A contoured continuum surface force model for particle methods

A surface tension model is essential to simulate multiphase flows with deformed interfaces. This study develops a contoured continuum surface force (CCSF) model for particle methods. A color function that varies sharply across the interface to mark different fluid phases is smoothed in the transition region, where the local contour curvature can be regarded as the interface curvature. The local contour passing through each reference particle in the transition region is extracted from the local profile of the smoothed color function. The local contour curvature is calculated based on the Taylor series expansion of the smoothed color function, whose derivatives are calculated accurately according to the definition of the smoothed color function. Two schemes are proposed to specify the smooth radius: fixed scheme, where 2×re (re = particle interaction radius) is assigned to all particles in the transition region; and varied scheme, where re and 2×re are assigned to the central and edged particles in the transition region respectively. Numerical examples, including curvature calculation for static circle and ellipse interfaces, deformation of square droplet to a circle (2D and 3D), droplet deformation in shear flow, and droplet coalescence, are simulated to verify the CCSF model and compare its performance with those of other methods. The CCSF model with the fixed scheme is proven to produce the most accurate curvature and lowest parasitic currents among the tested methods.