[论文解读] Spray formation: a numerical closeup

Yue Ling, Daniel Fuster+2 more|arXiv (Cornell University)|Nov 13, 2015

Fluid Dynamics and Heat Transfer被引用 40

一句话总结

本研究采用高达40亿个计算单元的三维直接数值模拟（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）以实现高分辨率模拟。
• 设置四种网格分辨率等级（M0至M3），评估收敛性及分辨率对液滴尺寸统计的影响。
• 采用动态自适应时间步长方案，以在界面快速变形与破碎事件中保持数值稳定性。
• 对下游采样区域的液滴尺寸数据进行统计分析及概率分布拟合（Log-Normal与Gamma分布）。

实验结果

研究问题

• RQ1 当前计算资源能否实现气液剪切层中初级雾化过程的完全解析直接数值模拟？
• RQ2 网格分辨率如何影响丝状物、液滴的形成及最终液滴尺寸分布？
• RQ3 液滴形成的主要机理是什么——是端缘丝状物的撕裂还是薄膜中孔洞的扩展？这些机理在模拟中如何被解析？
• RQ4 经典液滴尺寸分布模型（Log-Normal与Gamma）在多大程度上能准确描述模拟数据？
• RQ5 分辨率不足是否导致人为破碎，从而产生虚假的小液滴及错误的尺寸统计？

主要发现

• 最细网格（M3）约40亿个计算单元，可完整解析所有物理尺度，结果表明波的演化与液膜形成是丝状物与液滴生成的关键因素。
• 粗网格（M0–M2）无法解析波峰曲率与剪切层结构，导致错误的、过早的破碎及虚假小液滴的生成。
• 丝状物通过两种机制形成：Taylor-Culick端缘的指状撕裂，以及薄液膜中孔洞的扩展，其中孔洞演化过程与Taylor-Culick理论相符。
• Log-Normal分布对液滴尺寸概率密度函数（PDF）的拟合效果优于Gamma分布，表明模拟中具备足够的动态范围以支持尺度律假设。
• M3中PDF的统计振荡源于采样时间较短，而非双峰过程，且在尺寸分布中未发现明显的“边缘液滴”或“薄膜液滴”特征。
• 液滴的聚并现象极为罕见，表明在此系统中聚集对最终液滴尺寸分布的影响可忽略。