Wetting boundary conditions for multicomponent pseudopotential lattice Boltzmann

Coelho, Rodrigo C. V.; Moura, Catarina B.; Gama, M. M. Telo da; Araújo, Nuno A. M.

doi:10.1002/fld.4988

Cited by 14 publications

(4 citation statements)

References 47 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Moreover, in both the two schemes, the local geometry information (e.g. local curvature) is not explicitly incorporated, indicating inconsistency of measured contact angles in different pore structures, as also noted in Coelho et al (2021).…”

Section: Geometric-function-based Two-component Contact Angle Schemementioning

confidence: 97%

Lattice Boltzmann modelling of isothermal two-component evaporation in porous media

et al. 2023

View full text Add to dashboard Cite

A mesoscopic lattice Boltzmann model is implemented for modelling isothermal two-component evaporation in porous media. The model is based on the pseudopotential multiphase model with two components to mimic the phase-change component (e.g. water and its vapour) and the non-condensible component (e.g. dry air), and the cascaded collision operator is used to enhance the numerical performance. The model is first analysed based on Chapman–Enskog expansion and then validated by the theoretical solution of an isothermal diffusive evaporation problem. We then discuss in detail the implementation of wettability based on a geometric function scheme and further validate the model with microfluidic evaporation experiments. We apply the method to simulate the convective evaporation of a dual-porosity medium and investigate the effects of inflow vapour concentration ( ${Y_{vapour,in}}$ ) and contact angle ( $\theta$ ) on the evaporation. Simulation results reproduce the typical transition from the constant evaporation regime (CRP) at large liquid saturation (S) to the receding front period (RFP) at small S, with an intermediate falling rate period in between. The dependence of the average evaporation rates on ${Y_{vapour,in}}$ and $\theta$ during CRP and RFP is investigated. A universal scaling formulation for the evaporation rate during CRP is found with respect to the concentration-related mass transfer number $B_Y$ , contact angle $\theta$ and inflow Reynolds number Re, i.e. $E{R_{CRP}} = {k_3}\ln \left ( {1 + {B_Y}} \right ) {\cdot } \left [ {\ln \left ( {1 + {Re}} \right ) + {k_2}} \right ]\left [ {\cos (\theta ) + {k_1}} \right ]$ , where ${k_1}$ , ${k_2}$ and ${k_3}$ are fitting parameters.

show abstract

Section: Geometric-function-based Two-component Contact Angle Schemementioning

confidence: 97%

Lattice Boltzmann modelling of isothermal two-component evaporation in porous media

et al. 2023

View full text Add to dashboard Cite

show abstract

“…We implement the wetting boundary condition as in references [19,20], which is summarized as follows. A field φ(x) is set to one at the solid nodes and zero at the fluid nodes.…”

Section: Wettingmentioning

confidence: 99%

Dynamics of two-dimensional liquid bridges

Coelho

Gazola

et al. 2022

J. Phys.: Condens. Matter

Self Cite

View full text Add to dashboard Cite

We have simulated the motion of a single vertical, two-dimensional liquid bridge spanning the gap between two flat, horizontal solid substrates of given wettabilities, using a multicomponent pseudopotential lattice Boltzmann method. For this simple geometry, the Young–Laplace equation can be solved (quasi-)analytically to yield the equilibrium bridge shape under gravity, which provides a check on the validity of the numerical method. In steady-state conditions, we calculate the drag force exerted by the moving bridge on the confining substrates as a function of its velocity, for different contact angles and Bond numbers. We also study how the bridge deforms as it moves, as parametrized by the changes in the advancing and receding contact angles at the substrates relative to their equilibrium values. Finally, starting from a bridge within the range of contact angles and Bond numbers in which it can exist at equilibrium, we investigate how fast it must move in order to break up.

show abstract

“…We implement the wetting boundary condition proposed in Ref. [17], which is summarized as follows. A field φ(x) is set to one at the solid nodes and zero at the fluid nodes.…”

Section: Wettingmentioning

confidence: 99%

Dynamics of two-dimensional liquid bridges

Coelho,

Cordeiro,

Gazola

et al. 2021

Preprint

Self Cite

View full text Add to dashboard Cite

We have simulated the motion of a single vertical, two-dimensional liquid bridge spanning the gap between two flat, horizontal solid substrates of given wettabilities, using a multicomponent pseudopotential lattice Boltzmann method. For this simple geometry, the Young-Laplace equation can be solved (quasi-)analytically to yield the equilibrium bridge shape under gravity, which provides a check on the validity of the numerical method. In steady-state conditions, we calculate the drag force exerted by the moving bridge on the confining substrates as a function of its velocity, for different contact angles and Bond numbers. We also study how the bridge deforms as it moves, as parametrized by the changes in the advancing and receding contact angles at the substrates relative to their equilibrium values. Finally, starting from a bridge within the range of contact angles and Bond numbers in which it can exist at equilibrium, we investigate how fast it must move in order to break up.

show abstract

Wetting boundary conditions for multicomponent pseudopotential lattice Boltzmann

Cited by 14 publications

References 47 publications

Lattice Boltzmann modelling of isothermal two-component evaporation in porous media

Lattice Boltzmann modelling of isothermal two-component evaporation in porous media

Dynamics of two-dimensional liquid bridges

Dynamics of two-dimensional liquid bridges

Contact Info

Product

Resources

About