[논문 리뷰] Linearly decoupled energy-stable numerical methods for multi-component two-phase compressible flow

Jisheng Kou, Shuyu Sun외 1명|arXiv (Cornell University)|2017. 12. 06.

Computational Fluid Dynamics and Aerodynamics참고 문헌 2인용 수 3

한 줄 요약

이 논문은 펑-로빈슨 상태방정식(예: Peng-Robinson)을 사용한 다성분 이상태 압축성 유체의 선형, 분리형, 에너지 안정성 수치 해법을 두 가지 제안한다. 중간 속도와 성분별 스칼라 보조변수(SAV) 방법을 도입함으로써 속도-몰 농도 상호작용을 분리하고, 각 시간 단계에서 선형 시스템만을 풀어, 무조건적인 에너지 소산을 보장하며, 복잡한 유체 시스템의 효율적이고 안정적인 시뮬레이션을 가능하게 한다.

ABSTRACT

In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure the consistency with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods preserve the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.

연구 동기 및 목표

• 실제 상태방정식(예: Peng-Robinson)을 사용하는 다성분 이상태 압축성 유체의 효율적인 수치 해법을 개발하여 열역학적으로 일관되며 압력 방정식 수립을 피하는 것.
• 질량 및 운동량 균형 방정식에서 속도와 몰 농도 간의 강한 결합을 해결하여 수치 해법의 복잡성을 줄이는 것.
• 무조건적인 이산 에너지 소산 법칙을 유지하는 선형, 분리형 해법을 구축하여 수치적 안정성을 확보하는 것.
• 다성분 유체에 대해 계산을 단순화하면서도 에너지 안정성을 유지할 수 있도록 성분별 SAV 방법을 도입하는 것.
• 이산 수치 설정에서 질량 균형 방정식과 운동량 방정식 간의 일관성을 확보하기 위해 중간 속도를 이산 형식에 통합하는 것.

제안 방법

• 비선형 헬름홀츠 자유 에너지 밀도를 다루기 위해 스칼라 보조변수(SAV) 방법을 활용하여 에너지 안정성 수치 해법을 선형 시스템으로 변환한다.
• 질량 균형 방정식의 대류 항에서 속도와 몰 농도 간의 밀접한 결합을 분리하기 위해 두 개의 중간 속도를 도입한다.
• 중간 속도를 포함한 이산 운동량 방정식을 설계하여 질량 균형 방정식과의 일관성을 확보하고 운동 에너지 변화를 유지한다.
• 각 성분에 대해 별도의 선형 질량 균형 방정식을 풀 수 있도록 성분별 SAV 수식을 제안하여 계산 효율을 크게 향상시킨다.
• 압력 방정식이 필요 없는 열역학적으로 일관된 모델링을 가능하게 하는 NVT 기반 모델링 프레임워크를 활용하며, 몰 수, 부피, 온도를 기본 변수로 사용한다.
• 이산 에너지 함수를 구성하고 시간 단계 간의 감쇠를 보여줌으로써 무조건적인 에너지 소산을 수학적으로 증명한다.

실험 결과

연구 질문

• RQ1 실제 상태방정식(예: Peng-Robinson)을 사용하는 다성분 이상태 압축성 유체에 대해 에너지 안정성, 선형성, 분리형 수치 해법을 어떻게 구성할 수 있는가?
• RQ2 질량 균형 방정식의 대류 항에서 속도-몰 농도 상호작용을 어떻게 분리하면서도 에너지 안정성을 유지할 수 있는가?
• RQ3 이산 설정에서 질량 균형 방정식과 운동량 방정식 양쪽과 일관성을 유지할 수 있도록 중간 속도를 어떻게 설계할 수 있는가?
• RQ4 계산 효율성을 향상시키기 위해 SAV 방법을 다성분 체계에 성분별로 어떻게 확장할 수 있는가?
• RQ5 성분별 SAV 방법의 사용이 다성분 압축성 유체의 수치 해법에 안정성과 정확성에 어떤 영향을 미치는가?

주요 결과

• 제안된 해법은 이산 에너지 함수를 통해 엄밀히 증명되어 각 시간 단계에서 감소하는 무조건적인 에너지 소산성을 보인다.
• 각 시간 단계에서 선형 방정식의 단순한 시퀀스만을 풀면 되므로, 비선형 결합 해법 대비 계산 효율성이 크게 향상된다.
• 성분별 SAV 방법을 통해 헬름홀츠 자유 에너지를 개별 성분 기여도로 분리하여 다성분 유체의 효율적 처리가 가능해진다.
• 수치 결과는 장기간 시뮬레이션 동안 복잡한 유동 패턴(예: 유동 쿼버, 속도 크기 등)을 효과적으로 시뮬레이션할 수 있음을 보여준다.
• 중간 속도의 사용은 운동량 방정식과 질량 균형 방정식 간의 일관성을 확보하여 이산 시스템의 물리적 일관성을 유지한다.
• 실제 Peng-Robinson 매개변수를 사용한 3성분 체계(메탄, 펜탄, 디케인)에서 해법이 다양한 시간 간격에 걸쳐 뛰어난 강건성과 안정성을 보였다.