Modeling surface tension of a two‐dimensional droplet using smoothed particle hydrodynamics

Surface tension plays a significant role at the dynamic interface of free‐surface flows especially at the microscale in capillary‐dominated flows. A model for accurately predicting the formation of two‐dimensional viscous droplets in vacuum or gas of negligible density and viscosity resulting from axisymmetric oscillation due to surface tension is solved using smoothed particle hydrodynamics composed of the Navier‐Stokes system and appropriate interfacial conditions for the free‐surface boundaries. The evolution of the droplet and its free‐surface interface is tracked over time to investigate the effects of surface tension forces implemented using a modified continuous surface force method and is compared with those performed using interparticle interaction force. The dynamic viscous fluid and surface tension interactions are investigated via a controlled curvature model and test cases of nonsteady oscillating droplets; attention is focused here on droplet oscillation that is released from an initial static deformation. Accuracy of the results is attested by demonstrating that (i) the curvature of the droplet that is controlled; (ii) uniform distribution of fluid particles; (iii) clean asymmetric forces acting on the free surface; and (iv) nonsteady oscillating droplets compare well with analytical and published experiment findings. The advantage of the proposed continuous surface force method only requires the use of physical properties of the fluid, whereas the interparticle interaction force method is restricted by the requirement of tuning parameters.     

A semi‐implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows

We present an implementation of Hysing's (Int. J. Numer. Meth. Fluids 2006; 51 :659–672) semi‐implicit method for treating surface tension, for finite volume models of interfacial flows. Using this method, the surface tension timestep restriction, which is often very stringent, can be exceeded by at least a factor of 5 without destabilizing the solution. The surface tension force in this method consists of an explicit part, which is the regular continuum surface force, and an implicit part which represents the diffusion of velocities induced by surface tension on fluids interfaces. The surface tension force is applied to the velocity field by solving a system of equations iteratively. Since the equations are solved only near interfaces, the computational time spent on the iterative procedure is insignificant. Copyright © 2008 John Wiley & Sons, Ltd.     

Issue Information

An integrated finite element method (FEM) is proposed to simulate incompressible two‐phase flows with surface tension effects, and three different surface tension models are applied to the FEM to investigate spurious currents and temporal stability. A Q2Q1 element is adopted to solve the continuity and Navier–Stokes equations and a Q2‐iso‐Q1 to solve the level set equation. The integrated FEM solves pressure and velocity simultaneously in a strongly coupled manner; the level set function is reinitialized by adopting a direct approach using interfacial geometry information instead of solving a conventional hyperbolic‐type equation. In addition, a consistent continuum surface force (consistent CSF) model is utilized by employing the same basis function for both surface tension and pressure variables to damp out spurious currents and to estimate the accurate pressure distribution. The model is further represented as a semi‐implicit manner to improve temporal stability with an increased time step. In order to verify the accuracy and robustness of the code, the present method is applied to a few benchmark problems of the static bubble and rising bubble with large density and viscosity ratios. The Q2Q1‐integrated FEM coupled with the semi‐implicit consistent CSF demonstrates the significantly reduced spurious currents and improved temporal stability. The numerical results are in good qualitative and quantitative agreements with those of the existing studies. Copyright © 2015 John Wiley & Sons, Ltd.     

Incompressible multiphase flow and encapsulation simulations using the moment‐of‐fluid method

A moment‐of‐fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment‐of‐fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi‐Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block‐structured adaptive grid. The new method is applied to multiphase problems illustrating contact‐line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two‐dimensional, three‐dimensional axisymmetric (R–Z), and three‐dimensional (X–Y–Z) coordinate systems. Copyright © 2015 John Wiley & Sons, Ltd.     

An improved method for calculation of interface pressure force in PLIC-VOF methods

A method for the application of interface force in the computational modeling of free surfaces and interfaces which uses PLIC-VOF methods is developed. This method is based on evaluation of the surface tension force only in the interfacial cells with out using the neighboring cells. The normal and the interface surface area needed for the calculation of the surface tension force are calculated more accurately. This method is applied on a staggered grid and it is referred to as Staggered Grid Interfere Pressure calculation method or SGIP. The present method is applied to a two-dimensional motionless liquid drop and a gas bubble. In addition, oscillations of a non-circular two-dimensional drop and a bubble due only to the surface tension forces are considered. It is shown that the new method predicts the pressure jump at the interface more accurately and produces less spurious currents compared to CSF, CSS and Meier's methods when applied to the same cases.     

A fast method for solving fluid–structure interaction problems numerically

The paper presents a semi‐implicit algorithm for solving an unsteady fluid–structure interaction problem. The algorithm for solving numerically the fluid–structure interaction problems was obtained by combining the backward Euler scheme with a semi‐implicit treatment of the convection term for the Navier–Stokes equations and an implicit centered scheme for the structure equations. The structure is governed either by the linear elasticity or by the non‐linear St Venant–Kirchhoff elasticity models. At each time step, the position of the interface is predicted in an explicit way. Then, an optimization problem must be solved, such that the continuity of the velocity as well as the continuity of the stress hold at the interface. During the Broyden, Fletcher, Goldforb, Shano (BFGS) iterations for solving the optimization problem, the fluid mesh does not move, which reduces the computational effort. The term ‘semi‐implicit’ used for the fully algorithm means that the interface position is computed explicitly, while the displacement of the structure, velocity and the pressure of the fluid are computed implicitly. Numerical results are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Semi‐Lagrangian semi‐implicit time‐splitting scheme for the shallow water equations

Time‐splitting technique applied in the context of the semi‐Lagrangian semi‐implicit method allows the use of extended time steps mainly based on physical considerations and reduces the number of numerical operations at each time step such that it is approximately proportional to the number of the points of spatial grid. To control time growth of the additional truncation errors, the standard stabilizing correction method is modified with no penalty for accuracy and efficiency of the algorithm. A linear analysis shows that constructed scheme is stable for time steps up to 2h. Numerical integrations with actual atmospheric fields of pressure and wind confirm computational efficiency, extended stability and accuracy of the proposed scheme. Copyright © 2006 John Wiley & Sons, Ltd.     

A variationally consistent finite element approach to the two‐fluid internal contact‐line problem

A new method for simulating incompressible viscous fluid flow involving moving internal contact lines is presented. The steady state interface shape is determined by a variationally consistent formulation of the surface tension contribution to the equations of motion adapted to the case of internal contact lines through the application of a global force balance compatibility condition that consistently removes the pressure indeterminacy. The Crouzeix–Raviart element is chosen to capture the pressure discontinuity at the two‐fluid interface. The resulting discrete equations are solved by an iterative procedure which displays fast convergence characteristics for small capillary numbers. Numerical results for the case of the steady movement of a fluid meniscus in a two‐dimensional channel are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Simulating gas–liquid multiphase flow using meshless Lagrangian method

A multiphase flow model has been established based on a moving particle semi‐implicit method. A surface tension model is introduced to the particle method to improve the numerical accuracy and stability. Several computational techniques are employed to simplify the numerical procedure and further improve the accuracy. A particle fraction multiphase flow model is developed and verified by a two‐phase Poiseuille flow. The multiphase surface tension model is discussed in detail, and an ethanol drop case is introduced to verify the surface tension model. A simple dam break is simulated to demonstrate the improvements with various modifications in particle method along with a new boundary condition. Finally, we simulate several bubble rising cases to show the capacity of this new model in simulating gas–liquid multiphase flow with large density ratio difference between phases. The comparisons among numerical results of mesh‐based model, experimental data, and the present model, indicate that the new multiphase particle method is acceptable in gas–liquid multiphase fluids simulation. Copyright © 2014 John Wiley & Sons, Ltd.     

Improvement of stability in moving particle semi‐implicit method

Moving particle semi‐implicit (MPS) method is one of the particle methods, which can be used to analyze incompressible free surface flow without surface tracking by a mesh or a scalar quantity. However, MPS causes unphysical numerical oscillation of pressure with high frequencies. We proposed a new formulation for the source term of Poisson equation of pressure. The proposed source term consists of three parts, one main part and two error‐compensating parts. With proper selection of the coefficients for the error‐compensating parts, we can suppress the unphysical pressure oscillation. Smoother pressure distributions are obtained in hydrostatic pressure and dam break problems. Copyright © 2010 John Wiley & Sons, Ltd.     

On surface tension modelling using the level set method

The paper describes and compares the performance of two options for numerically representing the surface tension force in combination with the level set interface‐tracking method. In both models, the surface tension is represented as a body force, concentrated near the interface, but the technical implementation is different: the first model is based on a traditional level set approach in which the force is distributed in a band around the interface using a regularized delta function, whereas in the second, the force is partly distributed in a band around the interface and partly localized to the actual computational cells containing the interface. A comparative study, involving analysis of several two‐phase flows with moving interfaces, shows that in general the two surface tension models produce results of similar accuracy. However, in the particular case of merging and pinching‐off of interfaces, the traditional level set model of surface tension produces an error that results in non‐converging solutions for film‐like interfaces (i.e. ones involving large contact areas). In contrast, the second model, based on the localized representation of the surface tension force, displays consistent first‐order convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

Development of an improved divergence‐free‐condition compensated coupled framework to solve flow problems with time‐varying geometries

In this article a coupled version of the improved divergence‐free‐condition compensated method will be proposed to simulate time‐varying geometries by direct forcing immersed boundary method. The proposed method can be seen as a quasi‐multi‐moment framework due to the fact that the momentum equations are discretized by both cell‐centered and cell‐face velocity. For simulating time‐varying geometries, a semi‐implicit iterative method is proposed for calculating the direct forcing terms. Treatments for suppressing spurious force oscillations, calculating drag/lift forces, and evaluating velocity and pressure for freshly cells will also be addressed. In order to show the applicability and accuracy, analytical as well as benchmark problems will be investigated by the present framework and compared with other numerical and experimental results.     

An improved three-dimensional model for interface pressure calculations in free-surface flows

A three-dimensional method for the calculation of interface pressure in the computational modeling of free surfaces and interfaces is developed. The methodology is based on the calculation of the pressure force at the interfacial cell faces and is mainly designed for volume of fluid (VOF) interface capturing approach. The pressure forces at the interfacial cell faces are calculated according to the pressure imposed by each fluid on the portion of the cell face that is occupied by that fluid. Special formulations for the pressure in the interfacial cells are derived for different orientations of an interface. The present method, referred to as pressure calculation based on the interface location (PCIL), is applied to both static and dynamic cases. First, a three-dimensional motionless drop of liquid in an initially stagnant fluid with no gravity force is simulated as the static case and then two different small air bubbles in water are simulated as dynamic cases. A two-fluid, piecewise linear interface calculation VOF method is used for numerical simulation of the interfacial flow. For the static case, both the continuum surface force (CSF) and the continuum surface stress (CSS) methods are used for surface tension calculations. A wide range of Ohnesorge numbers and density and viscosity ratios of the two fluids are tested. It is shown that the presence of spurious currents (artificial velocities present in case of considerable capillary forces) is mainly due to the inaccurate calculation of pressure forces in the interfacial computational cells. The PCIL model reduces the spurious currents up to more than two orders of magnitude for the cases tested.

Also for the dynamic bubble rise case, it is shown that using the numerical solver employed here, without PCIL, the magnitude of spurious currents is so high that it is not possible to simulate this type of surface tension dominated flows, while using PCIL, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.     

Two‐dimensional modeling for stability analysis of two‐phase stratified flow

The effect of wavelength and relative velocity on the disturbed interface of two‐phase stratified regime is modeled and discussed. To analyze the stability, a small perturbation is imposed on the interface. Growth or decline of the disturbed wave, relative velocity, and surface tension with respect to time will be discussed numerically. Newly developed scheme applied to a two‐dimensional flow field and the governing Navier–Stokes equations in laminar regime are solved. Finite volume method together with non‐staggered curvilinear grid is a very effective approach to capture interface shape with time. Because of the interface shape, for any time advancement, a new grid is performed separately on each stratified field, liquid, and gas regime. The results are compared with the analytical characteristics method and one‐dimensional modeling. This comparison shows that solving the momentum equation including viscosity term leads to physically more realistic results. In addition, the newly developed method is capable of predicting two‐phase stratified flow behavior more precisely than one‐dimensional modeling. It was perceived that the surface tension has an inevitable role in dissipation of interface instability and convergence of the two‐phase flow model. Copyright © 2009 John Wiley & Sons, Ltd.     

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

In this paper, we present a finite element method for two‐phase incompressible flows with moving contact lines. We use a sharp interface Navier–Stokes model for the bulk phase fluid dynamics. Surface tension forces, including Marangoni forces and viscous interfacial effects, are modeled. For describing the moving contact lines, we consider a class of continuum models that contains several special cases known from the literature. For the whole model, describing bulk fluid dynamics, surface tension forces, and contact line forces, we derive a variational formulation and a corresponding energy estimate. For handling the evolving interface numerically, the level‐set technique is applied. The discontinuous pressure is accurately approximated by using a stabilized extended finite element space. We apply a Nitsche technique to weakly impose the Navier slip conditions on the solid wall. A unified approach for discretization of the (different types of) surface tension forces and contact line forces is introduced. Results of numerical experiments are presented, which illustrate the performance of the solver. Copyright © 2016 John Wiley & Sons, Ltd.     

The Floating Ball Paradox

In capillary theory there are two kinds of surface tension. There is the surface tension at the interface between two immiscible fluids. Thomas Young [9] also allowed for there to be a surface tension associated with a liquid-solid interface. He proceeded to use a balance of forces argument to derive the well-known contact angle condition along a liquid-liquid-solid intersection. The validity of this argument has recently been called into question by R. Finn [6]. A floating ball experiment discussed in that paper leads to an apparent paradox. We address this issue.      

An assessment of a parallel,finite element method for three‐dimensional,moving‐boundary flows driven by capillarity for simulation of viscous sintering

A parallel, finite element method is presented for the computation of three‐dimensional, free‐surface flows where surface tension effects are significant. The method employs an unstructured tetrahedral mesh, a front‐tracking arbitrary Lagrangian–Eulerian formulation, and fully implicit time integration. Interior mesh motion is accomplished via pseudo‐solid mesh deformation. Surface tension effects are incorporated directly into the momentum equation boundary conditions using surface identities that circumvent the need to compute second derivatives of the surface shape, resulting in a robust representation of capillary phenomena. Sample results are shown for the viscous sintering of glassy ceramic particles. The most serious performance issue is error arising from mesh distortion when boundary motion is significant. This effect can be severe enough to stop the calculations; some simple strategies for improving performance are tested. Copyright © 2001 John Wiley & Sons, Ltd.     

A filter‐based,mass‐conserving lattice Boltzmann method for immiscible multiphase flows

Some issues of He–Chen–Zhang lattice Boltzmann equation (LBE) method (referred as HCZ model) (J. Comput. Physics 1999; 152 :642–663) for immiscible multiphase flows with large density ratio are assessed in this paper. An extended HCZ model with a filter technique and mass correction procedure is proposed based on HCZ's LBE multiphase model. The original HCZ model is capable of maintaining a thin interface but is prone to generating unphysical oscillations in surface tension and index function at moderate values of density ratio. With a filtering technique, the monotonic variation of the index function across the interface is maintained with larger density ratio. Kim's surface tension formulation for diffuse–interface method (J. Comput. Physics 2005; 204 :784–804) is then used to remove unphysical oscillation in the surface tension. Furthermore, as the density ratio increases, the effect of velocity divergence term neglected in the original HCZ model causes significant unphysical mass sources near the interface. By keeping the velocity divergence term, the unphysical mass sources near the interface can be removed with large density ratio. The long‐time accumulation of the modeling and/or numerical errors in the HCZ model also results in the error of mass conservation of each dispersed phase. A mass correction procedure is devised to improve the performance of the method in this regard. For flows over a stationary and a rising bubble, and capillary waves with density ratio up to 100, the present approach yields solutions with interface thickness of about five to six lattices and no long‐time diffusion, significantly advancing the performance of the LBE method for multiphase flow simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulations of droplet formation in a cross‐junction microchannel by the lattice Boltzmann method

An immiscible liquid–liquid multiphase flow in a cross‐junction microchannel was numerically studied using the lattice Boltzmann method. An improved, immiscible lattice BGK model was proposed by introducing surface tension force based on the continuum surface force (CSF) method. Recoloring step was replaced by the anti‐diffusion scheme in the mixed region to reduce the side‐effect and control the thickness of the interface. The present method was tested by the simulation of a static bubble. Laplace's law and spurious velocities were examined. The results show that our model is more advantageous for simulations of immiscible fluids than the existing immiscible lattice BGK models. Computational results of multiphase flow in a cross‐junction microchannel were obtained and analyzed based on dimensionless numbers. It is found that the flow pattern is decided mostly by the capillary number at a small inlet flux. However, at the same capillary number, a large inlet flux will lead to much smaller droplet generation. For this case, the flow is determined by both the capillary number and the Weber number. Copyright © 2007 John Wiley & Sons, Ltd.