A consistent adaptive level set framework for incompressible two-phase flows with high density ratios and high Reynolds numbers
A consistent adaptive level set framework for incompressible two-phase flows with high density ratios and high Reynolds numbers
复制标题
DOI：
10.1016/j.jcp.2023.111971
复制
发表时间：
2023-01
期刊：
影响因子：
--
通讯作者：
Yadong Zeng;Han Liu;Q. Gao;A. Almgren;A. Bhalla;Lian Shen
中科院分区：
文献类型：
--
作者：
Yadong Zeng;Han Liu;Q. Gao;A. Almgren;A. Bhalla;Lian Shen
我们在多级并置网格布局中开发了一个一致的自适应框架，用于通过自适应网格细化（AMR）来模拟两相流。保守动量方程和质量方程在当前一致的框架下求解。这种一致的质量和动量传递处理极大地提高了模拟高密度比和高雷诺数的两相流的准确性和鲁棒性。界面捕获水平集方法与保守形式的纳维-斯托克斯方程相结合，并应用多级重新初始化技术来实现质量守恒。这种自适应框架允许我们使用子循环或非子循环方法逐级推进所有变量，以解耦每个级别上的数据推进。该框架的准确性和鲁棒性通过各种典型的两相流问题得到了验证。我们证明，一致的方案会在高密度比（高达 106）和高雷诺数（高达 106）的流动中产生数值稳定的解，而不一致的方案在这些测试中表现出非物理流体行为。此外，结果表明，子循环和非子循环方法提供了一致的结果，并且两者都可以准确地解析具有表面张力效应的两相流的界面。最后，模拟了 3D 破碎波问题，以显示所提出的使用 AMR 的框架的效率和显着加速。
We develop a consistent adaptive framework in a multilevel collocated grid layout for simulating two-phase flows with adaptive mesh refinement (AMR). The conservative momentum equations and the mass equation are solved in the present consistent framework. This consistent mass and momentum transport treatment greatly improves the accuracy and robustness for simulating two-phase flows with a high density ratio and high Reynolds number. The interface capturing level set method is coupled with the conservative form of the Navier–Stokes equations, and the multilevel reinitialization technique is applied for mass conservation. This adaptive framework allows us to advance all variables level by level using either the subcycling or the non-subcycling method to decouple the data advancement on each level. The accuracy and robustness of the framework are validated by a variety of canonical two-phase flow problems. We demonstrate that the consistent scheme results in a numerically stable solution in flows with high density ratios (up to 106) and high Reynolds numbers (up to 106), while the inconsistent scheme exhibits nonphysical fluid behaviors in these tests. Furthermore, it is shown that the subcycling and non-subcycling methods provide consistent results and that both of them can accurately resolve the interfaces of the two-phase flows with surface tension effects. Finally, a 3D breaking wave problem is simulated to show the efficiency and significant speedup of the proposed framework using AMR.