On the diffuse interface method using a dual-resolution Cartesian grid

On the diffuse interface method using a dual-resolution Cartesian grid
复制标题

基于双分辨率笛卡尔网格的漫反射界面方法

DOI:
10.1016/j.jcp.2014.05.005
复制
发表时间:
2014-09-15
影响因子:
4.1
通讯作者:
Yuan, Cheng-jun
Yuan, Cheng-jun
中科院分区:
物理与天体物理2区
文献类型:
--
作者:
Ding, Hang;Yuan, Cheng-jun

文献摘要

被引文献

相似文献

We investigate the applicability and performance of diffuse interface methods on a dual-resolution grid in solving two-phase flows. In the diffuse interface methods, the interface thickness represents a cut-off length scale in resolving the interfacial dynamics, and it was found that an apparent loss of mass occurs when the interface thickness is comparable to the length scale of flows [24]. From the accuracy and mass conservation point of view, it is desirable to have a thin interface in simulations. We propose to use a dual-resolution Cartesian grid, on which a finer resolution is applied to the volume fraction C than that for the velocity and pressure fields. Because the computation of C field is rather inexpensive compared to that required by velocity and pressure fields, dual-resolution grids can significantly increase the resolution of the interface with only a slight increase of computational cost, as compared to the single-resolution grid. The solution couplings between the fine grid for C and the coarse grid (for velocity and pressure) are delicately designed, to make sure that the interpolated velocity is divergence-free at a discrete level and that the mass and surface tension force are conserved. A variety of numerical tests have been performed to validate the method and check its performance. The dual-resolution grid appears to save nearly 70% of the computational time in two-dimensional simulations and 80% in three-dimensional simulations, and produces nearly the same results as the single-resolution grid. Quantitative comparisons are made with previous studies, including Rayleigh Taylor instability, steadily rising bubble, and partial coalescence of a drop into a pool, and good agreement has been achieved. Finally, results are presented for the deformation and breakup of three-dimensional drops in simple shear flows. (C) 2014 Elsevier Inc. All rights reserved.