权重积分方法计算界面张力消除表面“寄生”流

1 INTRODUCTION In simulation of multiphase flow in fixed Eule- rian grids by the level set approach and similar meth- ods, the accurate computation of surface tension force is critical. Two popular models were developed to model surface tension force: one is referred to as the continuum surface force (CSF) model[1―3] and the other is the continuum surface stress (CSS) model[4]. Both methods are commonly used for advecting inter- faces in the level set approach[5] and VOF (vol- ume-of-…

A Simple Weighted Integration Method for Calculating Surface Tension Force to Suppress Parasitic Flow in the Level Set Approach

Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.