[논문 리뷰] Spray formation: a numerical closeup

Yue Ling, Daniel Fuster외 2명|arXiv (Cornell University)|2015. 11. 13.

Fluid Dynamics and Heat Transfer인용 수 40

한 줄 요약

이 연구는 40억 개의 계산 격자까지 사용하여 기체-액체 혼합층에서의 스프레이 형성에 대한 대규모 3차원 직접 시뮬레이션(DNS)을 수행하여 주 원자화에서의 전체 물리 스케일을 해석한다. 결과적으로, 정밀한 방울 크기 분포를 정확히 기록하기 위해 미세한 격자 해상도가 필수적임을 밝혀냈으며, 이는 Log-Normal 및 Gamma 분포가 데이터에 적합하지만, Log-Normal 분포가 충분한 스케일 범위를 제공함에 따라 더 나은 적합도를 보임을 시사한다.

ABSTRACT

Spray formation and atomization in a gas-liquid mixing layer is an important fundamental problem of multiphase flows. It is highly desirable to visualize the detailed atomization process and to analyze the instabilities and mechanisms involved, and massive numerical simulations are required, in addition to experiment. Rapid development of numerical methods and computer technology in the past decades now allows large-scale three-dimensional direct numerical simulations of atomization to be performed. Nevertheless, the fundamental question, whether all the physical scales involved in the primary breakup process are faithfully resolved, remains unclear. In the present study, we conduct direct numerical simulations of spray formation in a gas-liquid mixing layer with state-of-the-art computational resources (using up to 4 billion cells and 16384 cores), in order to obtain a high-fidelity numerical closeup of the detailed mechanisms of spray formation. We also aim to examine whether present computational resources are sufficient for a fully resolved direct numerical simulation of atomization.

연구 동기 및 목표

• 현재의 계산 자원이 기체-액체 혼합층에서의 주 원자화 과정에서 모든 물리 스케일을 완전히 해상할 수 있는지 조사하기 위해.
• 격자 해상도가 방울 크기 분포 및 파손 역학의 정확도에 미치는 영향을 검토하기 위해.
• 고해상도 시뮬레이션을 활용하여 얇은 액체 필름에서 리간드, 방울 및 구멍 형성 메커니즘을 분석하기 위해.
• 기존의 방울 크기 분포 모델(Log-Normal 및 Gamma)이 해상된 원자화 과정의 맥락에서 타당한지 평가하기 위해.
• 粗격자 해상도로 인해 표면파 발달 및 시어층 곡률을 제대로 해상하지 못해 인위적인 작은 방울이 생성되는지 여부를 판단하기 위해.

제안 방법

• 최대 16,384개의 CPU 코어와 40억 개의 계산 격자를 사용하여 준평면적인 기체-액체 혼합층에 대한 직접 시뮬레이션(DNS)을 수행하였다.
• 액체-기체 계면을 추적하고 얇은 액체 필름 및 리간드를 해상하기 위해 가변 격자 해상도를 적용한 체적 분율(VOF) 방법을 사용하였다.
• 평면적 제트를 생성하기 위해 고체 분리판을 사용한 설정을 시뮬레이션하였으며, 높이 H = 0.8 mm로 물리 척도를 축소하여 고해상도 시뮬레이션을 가능하게 하였다.
• 4개의 격자 해상도 수준(M0부터 M3까지)을 변화시켜 수렴성과 해상도가 방울 크기 통계에 미치는 영향을 평가하였다.
• 빠른 계면 변형 및 파손 사건 동안 수치적 안정성을 유지하기 위해 동적 가변 시간 스텝 스킴을 적용하였다.
• 하류 샘플링 박스에서 확득한 방울 크기 데이터에 대해 통계 분석 및 확률 분포 적합(로그정규 및 감마 분포)을 수행하였다.

실험 결과

연구 질문

• RQ1 현재의 계산 자원으로 기체-액체 혼합층에서의 주 원자화 과정을 완전히 해상된 직접 시뮬레이션을 달성할 수 있는가?
• RQ2 격자 해상도가 리간드, 방울 형성 및 그로 인한 방울 크기 분포에 어떤 영향을 미치는가?
• RQ3 주요 방울 형성 메커니즘은 무엇인가—끝면 림에서의 리간드 파손인가, 아니면 얇은 필름에서의 구멍 팽창인가—그리고 이는 어떻게 해상되는가?
• RQ4 기존의 방울 크기 분포 모델(Log-Normal 및 Gamma)이 시뮬레이션된 데이터를 얼마나 정확히 기술하는가?
• RQ5 해상도가 열악할 경우 인위적인 분열이 발생하여 인위적인 작은 방울이 생성되고 잘못된 크기 통계가 도출되는가?

주요 결과

• 약 40억 개의 격자로 구성된 가장 미세한 격자(M3)는 전체 물리 스케일을 해상하였으며, 파동 발달 및 필름 형성이 리간드 및 방울 생성에 핵심적임을 보여주었다.
• 粗격자(M0–M2)는 파동 정점의 곡률과 시어층의 구조를 제대로 해상하지 못해 잘못된, 조기 발생하는 파손과 인위적인 작은 방울 생성을 유도하였다.
• 리간드 형성은 두 가지 메커니즘을 통해 발생한다: 타일러-쿠리클 끝면에서의 패턴 형성과 얇은 액체 필름 내의 구멍 팽창이며, 구멍의 진화는 타일러-쿠리클 이론과 일치한다.
• 확률 밀도 함수(PDF)에 대한 Log-Normal 분포는 Gamma 분포보다 더 나은 적합도를 보이며, 스케일링 가설에 필요한 충분한 동적 범위를 제공함을 시사한다.
• M3의 PDF에서 관찰된 통계적 진동은 이분산 과정이 아니라 짧은 샘플링 시간으로 인한 것으로 판단되며, 크기 분포에 명백한 '림 방울' 또는 '필름 방울'의 이중성 증거는 발견되지 않았다.
• 방울의 융합은 거의 관찰되지 않아 이 시스템에서 응집이 최종 방울 크기 분포에 미치는 영향은 미미하다는 것을 시사한다.