[논문 리뷰] Parallel variable-density particle-laden turbulence simulation

Hadi Pouransari, Milad Mortazavi외 1명|arXiv (Cornell University)|2016. 01. 20.

Particle Dynamics in Fluid Flows참고 문헌 14인용 수 13

한 줄 요약

이 논문은 저밀도 수流, 입자 랜딩 터뷸런스를 위한 병렬 MPI 기반 C++ 코드를 제시하며, 복사加열을 고려한 두 방향 결합 라그랑주-오일러 프레임워크를 사용한다. 주요 기여는 저마하 수류, 변동 밀도 유동에서 발생하는 변수 계수 푸아송 방정식을 해결하기 위한 새로운 병렬 선형 해법으로, 복사에 의한 난류 및 비중력 효과의 정확한 시뮬레이션을 가능하게 한다.

ABSTRACT

We have developed a fully parallel C++/MPI based simulation code for variable-density particle-laden turbulent flows. The fluid is represented through a uniform Eulerian staggered grid, while particles are modeled using a Lagrangian point-particle framework. Spatial discretization is second-order accurate, and time integration has a fourth-order accuracy. Two-way coupling of the particles with the background flow is considered in both momentum and energy equations. The code is fully modular and abstracted, and easily can be extended or modified. We have considered two different boundary conditions. We have also developed a novel parallel linear solver for the variable density Poisson equation that arises in the calculation.

연구 동기 및 목표

• 변동 밀도, 입자 랜딩 난류에 대한 확장 가능하고 모듈화되며 확장 가능한 시뮬레이션 프레임워크를 개발한다.
• 저마하 수류, 변동 밀도 유동에서 입자와 유체 간의 상호작용을 두 방향으로 정확히 모델링한다.
• 선호 농축, 선호 스위핑, 복사에 의한 난류와 같은 복잡한 현상을 시뮬레이션할 수 있도록 한다.
• 변동 밀도 유동에서 발생하는 변수 계수 푸아송 방정식을 해결하기 위한 새로운 병렬 선형 해법을 구현한다.
• 부소니스크 근사의 완화를 통해 균일한 등방성 난류와 비중력 유동을 모두 지원한다.

제안 방법

• 스토크스 항력과 복사 및 열전도를 통한 열 에너지 교환을 고려한 점입자 라그랑주 프레임워크를 사용한다.
• 밀도, 운동량, 온도와 같은 유체 변수를 두 번째 차수 정확도를 갖는 스타거드 격자로 표현하며, 네 단계 룬게-쿠타 시간 적분을 적용한다.
• 오일러 프레임과 라그랑주 프레임 간의 두 방향 결합을 위해 선형 보간 및 투영 연산자(I 및 P)를 적용한다.
• 유동 압력에 대한 변수 계수 푸아송 방정식을 해결하기 위해, 변동 밀도 효과를 처리하기 위한 새로운 병렬 반복 해법을 사용한다.
• 유체의 압력 일관성과 수치 오차로 인한 비정상 문제를 방지하기 위해, 수정된 룬게-쿠타 방법을 사용해 시간에 따라 변화하는 배경 압력 $ P_0 $ 를 계산한다.
• 스태거드 격자에서 두 번째 차수 보존형 유한차분 이산화를 통해 운동 에너지 보존을 확보한다.

실험 결과

연구 질문

• RQ1 입자에 대한 복사加열이 입자 랜딩 유동에서 난류의 발달에 어떤 영향을 미치는가?
• RQ2 유체 밀도의 변동이 두 방향 결합 입자 랜딩 유동에서 운동량 및 에너지 이송에 어떤 영향을 미치는가?
• RQ3 부소니스크 근사를 적용하지 않을 경우, 비중력 유도 난류 시뮬레이션에 어떤 영향을 미치는가?
• RQ4 병렬 보존형 수치 해법이 저마하, 변동 밀도 유동에서 변수 계수 푸아송 방정식을 정확히 해결할 수 있는가?
• RQ5 복사에 의한 열 기울기가 입자 선호 농축과 유체 난류 구조에 어떤 영향을 미치는가?

주요 결과

• 코드는 입자 선호 농축과 선호 스위핑과 같은 고전적 현상을 정확히 재현하여, 두 방향 결합 프레임워크의 타당성을 검증한다.
• 기존의 일정 밀도, 부소니스크 근사 기반 코드에서 관측된 바 있었던 복사에 의한 난류 효과를, 이제는 변동 밀도 조건에서도 정확히 재현한다 (Zamansky 등, 2014).
• 부소니스크 근사를 적용하지 않은 상태에서 균일한 비중력 유도 난류를 정확히 시뮬레이션하며, Batchelor 등 (1992)의 결과와 일치함을 확인한다.
• 변수 계수 푸아송 방정식을 위한 새로운 병렬 선형 해법을 개발하였으며, 변동 밀도 유동에서 수치적 안정성과 일관성을 유지함을 입증하였다.
• 시간에 따라 변화하는 배경 압력 $ P_0 $ 는 수정된 RK4 방법을 통해 계산되어 푸아송 방정식과의 일관성을 확보하고, 수치 오차로 인한 비정상 문제를 방지한다.
• 코드는 공간적으로 두 번째 차수 정확도와 시간적으로 네 번째 차수 정확도를 확보하였으며, 스타거드 격자에서 보존형 이산화를 통해 운동 에너지 보존을 달성한다.