In the present study, a robust and conservative numerical scheme is proposed to simulate the violent two-phase flows with high-density ratios. In this method, the mass conservation equation and the momentum equation are solved in a consistent manner. The tangent of hyperbola for interface capturing scheme is extended for the computation of the mass flux by which the sharpness and conservation property of density field is preserved. Compared with other recently proposed methods, no geometrical computation is involved in deriving the mass flux and the spurious velocity in the interfacial region can be completely avoided. To improve the computational efficiency, the present method is implemented on a parallel block-structured adaptive mesh refinement method with a staggered layout of variables. High-fidelity numerical simulation of plunging jet through the liquid surface is performed. A bubble detection algorithm is developed to track bubbles generated in air entrainment process. The evolution of the bubble cloud, air concentration, bubble-size, and bubble-velocity distributions are predicted and compared quantitatively with the experiment. Numerical results show the air entrainment and penetration depth are highly correlated with the upstream disturbance. The growing interfacial roughness of the jet yields more entrained air in the final stage of jet impingement. It is found that when the initial perturbation is introduced, the overall size of the equivalent bubble radius will expand, and the penetration depth of the bubble cloud will decrease, while a larger volume of air is entrained.
Surface tension in two-phase flow problems plays a dominant role in many micro-flow phenomena and has an important influence on the development of flow instability phenomena that contain free surfaces. In this study, the multi-moment finite volume method is extended for direct numerical simulation of two-phase flow problems. A constraint interpolation profile–CSL (semi-Lagrangian) scheme is used for discretization of the advection part in the momentum equation. A compact volume of fluid method–approximated piecewise linear calculation method without flux limiter is proposed for capturing the moving interface. For modeling the surface tension accurately, the logic in curvature estimation is redesigned based on the height function (HF) method. The isolated volumetric fractions that may reduce accuracy in HF integration are excluded, and the numerical solution shows that the accuracy in the curvature estimation is improved for a coarse mesh. The present method is implemented with a parallel block-structured adaptive mesh refinement (BAMR) strategy; thus, the computational cost can be reduced significantly. Numerical tests show that the present BAMR solver is capable of reproducing the theoretical predictions of capillary wave instability problems with high accuracy. The simulation of droplet collisions further demonstrates the accuracy of the surface tension model. Finally, we extend it to the liquid jet atomization. The wavy disturbance, film breakup, liquid filament pinch-off, and droplet generation are well reproduced. The droplet size distribution satisfies the experimental measurement and theoretical predictions power-law. BAMR shows a huge advantage in computational efficiency than the traditional Cartesian grid. The findings of this study can help for a better understanding of the micro-mechanism of surface tension driven flows.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.