Numerical Simulation of Moving Contact Line Problems Using a Volume-of-Fluid Method

Renardy, Michael; Renardy, Yuriko; Li, Jie

doi:10.1006/jcph.2001.6785

Cited by 257 publications

(159 citation statements)

References 19 publications

Supporting

Mentioning

159

Contrasting

Order By: Relevance

“…We first determine the value of the contact angle to apply at the wall. This value is then imposed as a boundary condition using relation (8) for the calculation of the capillary contribution (7) in the momentum balance (5). One objective of this work is to compare different possible modeling to the dynamic modeling introduced in our code JADIM (model Dyn2 in the following).…”

Section: Numerical Modeling Of the Contact Anglementioning

confidence: 99%

“…For the transport of the interface, we give here some examples showing that almost all the classical methods are concerned: Boundary Integral methods [1][2][3], adaptive grid methods [4,5], Level-Set methods [6,7], Volume of Fluid methods [8][9][10], Front Tracking Methods [11,12] and coupled Level set and Volume-of-Fluid (CLSVOF) methods [13]. Different methods have also been developed for the modeling of moving contact lines.…”

Section: Introductionmentioning

confidence: 99%

“…The simplest situation is then to impose a constant angle corresponding to the static angle, i.e. h W ¼ h S [8,11,7,6,9]. When a no-slip condition is imposed on the wall, the stress generated by a contact line moving at velocity U cl , can be estimated as…”

Section: Introductionmentioning

confidence: 99%

“…The stress at the contact line is clearly diverging when refining the grid size (see for example [14] where the evolution of the viscous stress at the wall is reported). Several authors [8,6,14] have dealt with the ''stress singularity'' paradox by introducing the Navier slip condition that gives a relation between the fluid velocity at the wall U W and a Navier slip length k N :…”

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

Comparison between numerical models for the simulation of moving contact lines

Legendre

Maglio

2015

Computers & Fluids

View full text Add to dashboard Cite

OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. a b s t r a c tThe aim of this study is to discus different numerically models for the simulation of moving contact lines in the context of a Volume of Fluid-Continuum Surface Force (VoF-CSF) method. We focus on the particular situation of spreading drops. We first present the numerical methods used for the simulation of moving contact line i.e. static contact angle versus dynamic contact angle, no slip condition versus slip condition. A grid and time convergence is performed for the different models. We show that the integration of the Continuum Surface Force using the finite volume method results in a grid dependence at the onset of the spreading. The static and dynamic models are compared to experiments. It is shown that the dynamic models based on the Cox's relation for the dynamic contact angle are able to reproduce experiments while static models overestimate the spreading time and are not able to reproduce the Tanner regime. The difference between static and dynamic models is shown to increase with the Ohnesorge number.

show abstract

Section: Numerical Modeling Of the Contact Anglementioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Comparison between numerical models for the simulation of moving contact lines

Legendre

Maglio

2015

Computers & Fluids

View full text Add to dashboard Cite

show abstract

“…The most important physics here is the motion of the contact line, which presents a well-known stress singularity that is conventionally removed by assuming ad hoc conditions such as Navier slip or numerical slip [49,50]. In recent years, the diffuse-interface model has emerged as a promising alternative that offers a more rational approach to this issue [33,[51][52][53].…”

Section: Drop Spreading On Partially Wetting Substratementioning

confidence: 99%

3D phase-field simulations of interfacial dynamics in Newtonian and viscoelastic fluids

Zhou

Yue

Feng

et al. 2010

Journal of Computational Physics

119

View full text Add to dashboard Cite

-This work presents a three-dimensional finite-element algorithm, based on the phase-field model, for computing interfacial flows of Newtonian and complex fluids. A 3D adaptive meshing scheme produces fine grid covering the interface and coarse mesh in the bulk. It is key to accurate resolution of the interface at manageable computational costs.The coupled Navier-Stokes and Cahn-Hilliard equations, plus the constitutive equation for non-Newtonian fluids, are solved using second-order implicit time stepping. Within each time step, Newton iteration is used to handle the nonlinearity, and the linear algebraic system is solved by preconditioned Krylov methods. The phase-field model, with a physically diffuse interface, affords the method several advantages in computing interfacial dynamics.One is the ease in simulating topological changes such as interfacial rupture and coalescence.Another is the capability of computing contact line motion without invoking ad hoc slip conditions. As validation of the 3D numerical scheme, we have computed drop deformation in an elongational flow, relaxation of a deformed drop to the spherical shape, and drop spreading on a partially wetting substrate. The results are compared with numerical and experimental results in the literature as well as our own axisymmetric computations where appropriate. Excellent agreement is achieved provided that the 3D interface is adequately resolved by using a sufficiently thin diffuse interface and refined grid. Since our model involves several coupled partial differential equations and we use a fully implicit scheme, the matrix inversion requires a large memory. This puts a limit on the scale of problems that can be simulated in 3D, especially for viscoelastic fluids.

show abstract

3‐D numerical simulations on flow and mixing behaviors in gas–liquid–solid microchannels

Tang

Liu

2013

AIChE Journal

View full text Add to dashboard Cite

Three-dimensional (3-D) gas-liquid-solid flow and mixing behaviors in microchannels were simulated by coupled volume of fluid and discrete phase method and simulations were validated against observations. The detachment time and length of gas slug are shortened in liquid-solid flow, compared with those in liquid flow due to higher superficial viscosity of liquid-solid mixture, which will move the bubble formation toward the dripping regime. Solid particles mainly distribute in liquid slug and particle flow shows obvious periodicity. With the increase of contact angle of the inner wall, gas slug (0-50 ), stratified (77-120 ), and liquid drop (160 ) flows are observed. The residence time distributions of solid and liquid phases are similar because particles behave as tracers. The backmixing of solid and liquid phases in liquid drop flow is the weakest among the three flow patterns, and the backmixing of gas phase in slug flow is weaker than that in both stratified and liquid drop flows. The results can provide a theoretical basis for the design of microreactors.

show abstract

Numerical Simulation of Moving Contact Line Problems Using a Volume-of-Fluid Method

Cited by 257 publications

References 19 publications

Comparison between numerical models for the simulation of moving contact lines

Comparison between numerical models for the simulation of moving contact lines

3D phase-field simulations of interfacial dynamics in Newtonian and viscoelastic fluids

3‐D numerical simulations on flow and mixing behaviors in gas–liquid–solid microchannels

Contact Info

Product

Resources

About