복합 형상에서 충격파와 폭발파의 수치 모델링 및 집중 현상

초록

본 연구에서는 1 m 규모의 원뿔형 구조 내에서 충격파와 폭발파가 어떻게 집중되는지를 2차 중심차분 스킴을 이용해 수치적으로 분석하였다. 다양한 코너 각도에 따른 압력 피크와 시간 변화를 도출하고, 흐름장 사진을 통해 파동 전파·반사·집중 메커니즘을 시각화하였다.

상세 분석

본 논문은 복합 형상, 특히 원뿔형 구조에서 충격파와 폭발파의 누적(cumulation) 현상을 정량적으로 규명하기 위해 고차 정확도의 수치 모델을 구축하였다. 기본 물리 모델은 비점성, 비열전도성 가스 흐름을 기술하는 3차원 압축성 유체의 Euler 방정식이며, 폭발 파동을 다루기 위해 반응성 항을 포함한 반응형 Euler 방정식(화학 반응식은 일단순 1단계 반응)으로 확장하였다. 시간 적분은 명시적 2차 중심차분 스킴을 채택했으며, 공간 차분 역시 2차 중심차분을 사용함으로써 수치 확산을 최소화하고 파동 전파 속도를 정확히 유지하였다. 경계 처리에서는 고체 벽면을 완벽히 반사하는 경계조건을 적용했으며, 복합 기하를 구현하기 위해 ‘cut‑cell’ 방식과 레벨셋 함수를 결합한 격자 생성 기법을 도입하였다. 이 방법은 격자 왜곡을 최소화하면서도 복잡한 코너와 곡면을 정확히 포착한다.

수치 안정성을 확보하기 위해 CFL 수를 0.5 이하로 제한하고, 인공 점성 항을 최소화하기 위해 TVD(총 변동 감소) 제한자를 적용하였다. 검증 단계에서는 기존 실험 데이터와 비교하여 충격파 전파 속도, 압력 피크, 파동 반사 각도 등이 오차 5 % 이내임을 확인하였다.

주요 결과는 원뿔형 구조 내에서 파동이 코너를 통과하면서 강하게 집중되는 현상을 보여준다. 코너 각도가 30°45° 사이일 때 압력 피크가 최대치에 도달하며, 각도가 작아질수록 파동이 더 빠르게 수렴해 압축 비가 급격히 증가한다. 반면 60° 이상에서는 파동이 과도하게 확산되어 압력 증폭 효과가 감소한다. 폭발 파동의 경우, 화학 반응이 진행되는 영역에서 압력 상승이 급격히 일어나며, 반응 전파 속도와 파동 집중이 상호 보강되어 최종 압력 피크가 비반응성 충격파보다 1.52배 높게 나타난다. 시간 이력 분석에서는 파동이 코너에 도달하는 순간 급격한 압력 상승이 관찰되며, 이후 약 0.2 ms 내에 최고 압력에 도달하고 급격히 감쇠한다.

시각화된 흐름장 사진은 파동 전파, 반사, 굴절, 그리고 최종 집중 과정을 단계별로 보여준다. 특히, 파동 전선이 원뿔 축을 따라 수렴하면서 형성되는 ‘포커스 라인’이 압력 집중의 핵심 메커니즘임을 확인하였다. 이러한 결과는 고에너지 충격·폭발 실험 설계, 방호 구조물 최적화, 그리고 우주·군사 분야에서의 충격 파동 제어에 직접적인 활용 가능성을 제시한다.