[논문 리뷰] Numerical Solution of Cylindrically Converging Shock Waves

M. Shabouei, Reza Ebrahimi외 1명|arXiv (Cornell University)|2016. 02. 08.

Computational Fluid Dynamics and Aerodynamics참고 문헌 28인용 수 4

한 줄 요약

이 논문은 TVD-MUSCL 보간법과 Roe의 근사 리만 해법을 사용한 유한체적 방법으로 원통좌표계에서의 원통형 수축 충격파를 수치 시뮬레이션한다. 연구는 축 근처에서 충격파 역학을 정확하게 해석하고, 기존의 분석적 및 수치적 기준과의 비교를 통해 강력한 수축 충격파를 다룰 수 있음을 입증한다.

ABSTRACT

The cylindrically converging shock wave was numerically simulated by solving the Euler equations in cylindrical coordinates with TVD scheme and MUSCL approach, using Roe's approximate Riemann solver and super-bee nonlinear limiter. The present study used the in house code developed for this purpose. The behavior of the solution in the vicinity of axis is investigated and the results of the numerical solution are compared with the computed data given by Payne, Lapidus, Abarbanel, and Goldberg, Sod, and Leutioff et al.

연구 동기 및 목표

• 원통형 수축 충격파의 수치적 안정성과 정확성을 확보하여 시뮬레이션하는 것.
• 특이점과 높은 기울기를 포함하는 축 근처에서의 수치 해법의 거동을 조사하는 것.
• Payne, Sod 등이 발표한 분석적 및 수치적 자료와의 비교를 통해 자체 개발한 계산 코드를 검증하는 것.
• TVD-MUSCL 스킴에 초등비선형 제한기(super-bee limiter)를 적용하여 원통 기하학에서 강한 충격파를 얼마나 잘 해상도를 확보할 수 있는지 평가하는 것.
• 관성 압축 융합 및 고에너지 밀도 물리학 분야에서 충격파 수렴 현상을 연구하기 위한 신뢰할 수 있는 수치 프레임워크를 제공하는 것.

제안 방법

• 2차 정밀도를 확보하기 위해 MUSCL 재구성 기법을 사용한 원통좌표계에서의 오일러 방정식을 유한체적 방법으로 해석하는 것.
• 불연속성 근처에서의 임의의 진동을 방지하고 안정성을 확보하기 위해 총 변동 감소(Total Variation Diminishing, TVD) 기법을 적용하는 것.
• 세포 경계에서의 수치적 플럭스를 계산하기 위해 Roe의 근사 리만 해법을 사용하여 충격파를 정확하게 포착하는 것.
• 수치적 확산을 통제하고 충격파 전면을 뚜렷하게 유지하기 위해 초등비선형 제한기(super-bee)를 사용하는 것.
• 원통형 충격파 시뮬레이션에 특화된 자체 개발한 계산 코드를 구현하는 것.
• 특이점 처리와 보존 법칙 유지의 목적으로 축에서 적절한 경계 조건을 적용하는 것.

실험 결과

연구 질문

• RQ1 TVD-MUSCL 유한체적 스킴은 축 근처에서 원통형 수축 충격파의 역학을 얼마나 정확하게 포착할 수 있는가?
• RQ2 초등비선형 제한기와 함께 사용된 Roe의 리만 해법은 원통 기하학에서 강한 충격파 전면을 얼마나 잘 해상도를 확보하는가?
• RQ3 수치적 결과는 Payne, Lapidus, Abarbanel, Goldberg, Sod, Leutioff 등의 기존 기준 솔루션과 어떻게 비교되는가?
• RQ4 축 근처에서 충격파 수렴을 시뮬레이션할 때 발생하는 수치적 과제는 무엇이며, 본 연구에서 어떻게 완화되는가?
• RQ5 자체 개발한 코드는 알려진 충격파 수렴의 물리적 거동을 어느 정도 정확하게 재현하는가?

주요 결과

• 수치적 해법은 기하학적 특이점이 존재하는 축 근처에서도 충격파 수렴을 성공적으로 포착하고 안정성을 유지한다.
• Payne, Lapidus, Abarbanel, Goldberg, Sod, Leutioff 등의 기준 자료와 양호한 일치를 보이며 시뮬레이션의 정확성을 입증한다.
• 초등비선형 제한기의 적용으로 수치적 확산이 효과적으로 감소하고 충격파 전면의 선명함이 유지된다.
• TVD-MUSCL 접근법은 강한 기울기가 존재하는 상황에서도 단조성을 보장하고 물리적으로 비합리적인 진동을 방지한다.
• 자체 개발한 코드는 원통좌표계에서의 복잡한 충격파 역학을 시뮬레이션하는 데 있어 강력하고 신뢰할 수 있는 성능을 보여준다.
• 시뮬레이션 결과는 충격파 집중 현상에 따른 압력 및 밀도 증가와 같은 예상되는 물리적 거동을 확인한다.