[논문 리뷰] The Consistent Multi-fluid Advection method

T. Plewa, Ewald Müller|ArXiv.org|1998. 07. 23.

Ionosphere and magnetosphere dynamics인용 수 13

한 줄 요약

이 논문은 질량 분율 정규화를 강제하고 수정된 보간 프로파일을 통해 수치적 확산을 줄임으로써 다중 유체 유동을 정확하게 운반하는 보존적이고 고차수의 Godunov 유형 방법인 일관된 다중유체운반(CMA) 방법을 소개한다. 이 방법은 보존성과 단조성을 유지하며, 특히 수치적 확산이 감소함으로써 중성자원자로 구성된 원소 농도가 몇 배 정도 변화할 수 있는 타입 II 초신성 모델에서 조성 해상도를 크게 향상시킨다.

ABSTRACT

Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to one in every computational zone throughout the simulation is fulfilled by renormalizing the mass fractions during the advection step. The proposed modification is appropriate for any difference scheme written in conservation form. Unlike other commonly used methods it does not violate the conservative character of the advection method. A new steepening mechanism, which is based on modification of interpolation profiles, is used to reduce numerical diffusion across composition discontinuities. Additional procedures are described, which are necessary to enforce monotonicity. Several numerical experiments are presented which demonstrate the capability of our Consistent Multi-fluid Advection (CMA) method in case of smooth and discontinuous distributions of fluid phases and under different hydrodynamic conditions. It is shown that due to the reduced diffusivity of the proposed scheme the abundance of some heavy elements obtained from hydrodynamic simulations of type II supernova explosions can change by a factor of a few in the most extreme cases.

연구 동기 및 목표

• 다중유체 유체역학 시뮬레이션에서 수치적 확산과 질량 분율의 일관성 없는 문제를 해결하기 위해.
• 여러 유체 성분을 운반할 때 고차수 Godunov 스킴의 보존적 성격을 유지하기 위해.
• 보존성이나 단조성을 위반하지 않으면서 조성 불연속성의 해상도를 향상시키기 위해.
• 천체물리학적 환경(예: 타입 II 초신성)에서 유체 혼합을 정확하게 모델링할 수 있도록 하기 위해.
• 보존 형식으로 작성된 임의의 차분 스킴에 적용 가능한 일반화 가능한 방법을 개발하기 위해.

제안 방법

• 모든 계산 영역에서 질량 분율의 합이 1이 되도록 하기 위해 운반 단계 동안 재정규화를 통해 제약 조건을 강제한다.
• 조성 인터페이스를 가로질러 수치적 확산을 줄이기 위해 보간 프로파일을 수정하여 새로운 급격한 변화 메커니즘을 창출한다.
• 보존 형식으로 작성된 임의의 고차수 Godunov 유형의 유한체적 방법과 호환된다.
• 해당 해법의 해가 물리적으로 비현실적인 진동을 일으키지 않도록 단조성을 보장하기 위해 추가 절차를 구현한다.
• 원래 스킴의 보존적 성격을 유지하며 인위적인 소스 항목을 피한다.
• 매끄럽고 불연속적인 유동 조건 하에서 성능을 테스트하기 위해 수치 실험을 여러 차례 수행한다.

실험 결과

연구 질문

• RQ1 조성 불연속성에서 정확성을 유지하면서도 보존적인 방식으로 다중유체 운반을 수행할 수 있는 방법은 무엇인가?
• RQ2 고차수 Godunov 스킴에 어떤 수정을 가해야 보존성을 위반하지 않으면서 수치적 확산을 줄일 수 있는가?
• RQ3 질량 분율의 재정규화가 다중유체 시뮬레이션의 안정성과 정확성에 어떤 영향을 미치는가?
• RQ4 새로운 보간 기반의 급격한 변화 메커니즘이 얼마나 높은 수준의 유체 인터페이스 해상도를 향상시키는가?
• RQ5 감소된 수치적 확산이 초신성 시뮬레이션에서 중원소 농도 예측에 어떤 영향을 미치는가?

주요 결과

• CMA 방법은 시뮬레이션 도메인 전반에서 질량 분율 정규화를 기계 정밀도 수준까지 유지한다.
• 수정된 보간 프로파일이 날카로운 조성 불연속성에서 수치적 확산을 크게 감소시킨다.
• 이 방법은 원래 Godunov 스킴의 보존적 성격을 유지하며 인위적인 소스 항목을 피한다.
• 추가적인 알고리즘 절차를 통해 단조성이 보장되어 비물리적 진동이 발생하지 않는다.
• 타입 II 초신성 시뮬레이션에서, 감소된 확산으로 인해 극단적인 경우 중원소 농도가 몇 배 정도 변화한다.
• 수치 실험 결과, 매끄럽고 불연속적인 유동 조건을 포함한 다양한 유체역학적 조건에서 안정적인 성능을 확인할 수 있었다.