## Gradient dynamics model for drops of volatile liquid on a porous substrate

Simon Hartmann 1, 2, ∗ and Uwe Thiele 1, 2, †

1 Institut f¨r u Theoretische Physik, Universit¨t a M¨ unster, Wilhelm-Klemm-Str. 9, D-48149 M¨ unster, Germany 2

Center for Nonlinear Science (CeNoS), Universit¨t M¨nster, Corrensstr. 2, 48149 M¨nster, Germany a u u

We present a mesoscopic hydrodynamic model for a spreading drop of volatile partially wetting liquid on a solid porous layer of small thickness. Thereby, evaporation takes place under strong confinement, i.e., we consider a drop that spreads on one of two parallel plates that form a narrow gap. Our gradient dynamics model describes the coupled dynamics of the vertically averaged vapor density profile, the drop height profile, and the vertically averaged saturation profile in the porous layer. The underlying free energy incorporates saturation-dependent wettability, capillarity of the drop and within the porous layer, air and vapor entropy, and an entropic contribution relevant close to complete filling of the porous layer. After developing the model we discuss the resulting sorption isotherm and illustrate typical drop spreading and imbibition behavior including the formation of a saturation halo within the porous substrate.

## I. INTRODUCTION

The behavior of liquids in and on porous media has been studied for a long time [1-6] due to the wide spectrum of natural and technical systems where this is relevant. Examples include the motion of particles with a porous outer layer in a liquid [7], mushy regions occurring during solidification of alloys [8], droplets of simple or complex liquids spreading on and within porous solids [5, 6, 9-17], reactive drops within a porous medium [18], the dewetting of a liquid film on a porous medium [19, 20], thermocapillary instabilities in systems involving porous layers [8, 21, 22], the infiltration of water into soil [23-25], the coating of porous substrates [26, 27], drops on polymer layers and brushes [17], and phase behavior in nanoporous materials without [28] and with light-switched properties [29].

Particular interest has focused on liquid flow past a porous medium [22, 30-38] as well as on the motion and stability of moving liquid-gas interfaces (imbibition fronts) within a porous medium [12, 23, 39-43] including cases where chemical reactions occur [44, 45]. Also the absorption of vapors into porous materials is of interest [46]. Note finally that the creation of liquid-infused slippery substrates is often based on porous substrate layers [47-51].

There are various approaches to modeling the dynamics of the imbibition process, i.e., the flow of liquid into and within porous material. In one approach the wet and dry regions within the porous medium are separated by a sharp imbibition front and the flow within the 'wet' part is described via Darcy's law v = -( κ/η ) ∇ p . It describes the velocity field v , which is driven by the pressure gradient ∇ p between the liquid reservoir and the Laplace pressure inside the pores of the porous medium [12, 13]. Here, κ denotes the permeability of the porous medium, and η the dynamic viscosity of the liquid.

A spreading droplet on a permeable membrane (porous medium with only vertical pores) is modeled by Davis and Hocking [6] using a thin-film equation with a sink term that is driven by the same pressure as the convective flux (see their section II). In particular, their Eqs. (1a-c) can be written as a gradient dynamics model [52] for the film height h :

$$\partial _ { t } h \, = \, \nabla \cdot \left [ Q \nabla \frac { \delta \mathcal { F } } { \delta h } \right ] - M \left ( \frac { \delta \mathcal { F } } { \delta h } - p _ { \text{vap} } \right ),$$

with the mobilities Q = h / η 3 3 , M = κ/η , the free energy F = ∫ [ σ ∂ h ( x ) 2 / 2 + ϱ liq gh / 2 2]d x , and p vap = 0. Here, σ is the interface tension, ϱ liq is the mass density of the liquid, and g is the gravitational acceleration. In this model, no lateral motion within the porous layer is accounted for.

Davis and Hocking [6] also describe in their section III a model for a droplet on a fully saturated porous substrate where lateral Darcy flows may occur. Furthermore, a model with an explicit wet-dry boundary within the porous layer is considered (their section IV). The process of imbibition into the vertical pores is included in Ref. [9], resulting in the coupled dynamics of height profiles on top and inside the porous medium, modeled through two coupled thin-film equations. Both height profiles of the form h i ( x, t ) are measured from the porous-liquid interface, i.e., every pore

∗ s.hartmann@uni-muenster.de; ORCID ID: 0000-0002-3127-136X

† u.thiele@uni-muenster.de; http://www.uwethiele.de; ORCID ID: 0000-0001-7989-9271

consists of a filled and an empty section. This approach does not allow for lateral transport within the porous layer and is also followed by Alleborn and Raszillier [11].

Partially filled pores can be accounted for in phase-field models of infiltration [23-25, 53, 54]. The governing equation describes the dynamics of a filling ratio within the porous medium and is mathematically closely related to a thin-film equation. Such models amend Richards equation [1], which itself extends Darcy's law to unsaturated porous media. In a related approach, Kap et al. [55] have recently presented a gradient dynamics model for a volatile droplet on a grafted polymer substrate, which describes the coupled dynamics of the droplet, the liquid within the substrate and the vapor, also cf. the review in Ref. [56].

Here, we adapt the latter approach by replacing the polymer brush substrate with a porous medium of small finite thickness, and by furthermore incorporating the driving by capillary pressure as employed by Cueto-Felgueroso and Juanes [25], Beltrame and Cajot [43]. The resulting model effectively extends Eq. (1) with additional equations for the vapor dynamics and the dynamics of the liquid within the pores. It can be applied to studies of spreading drops of volatile liquids on porous substrates that are initially unsaturated. The pursued gradient dynamics approach to mesoscopic hydrodynamics is based on the observation that in the Stokes limit, i.e., without inertia, the dynamics of a hydrodynamic system is simply driven by the minimization of a free energy functional. Such a thermodynamically consistent approach has for a long time been employed for the description of diffusive processes like the decomposition of binary mixtures via the Cahn-Hilliard equation [57-59], before it was noted by Mitlin [60] that the long-wave equation [61], that describes the dewetting of a thin layer of simple nonvolatile liquid, belongs to the same class of models [52]. Subsequently, the finding was extensively employed in the analysis of such systems and, more importantly, extended to more complex thin-film systems [56, 62-66]. Writing the hydrodynamic thin-film equations in gradient dynamics form automatically ensures that the derived equations respect the principles of thermodynamics, i.e., provides a convenient consistency check. This facilitates, in particular, the incorporation of additional effects via inclusion of further degrees of freedom and energy contributions [67].

The mesoscopic hydrodynamic model is presented in section II as a three-field gradient dynamics. Particular attention rests on the involved transport processes and the underlying energy functional. The subsequent section III discusses the resulting sorption isotherm of the porous medium, while section IV analyzes the coupled spreading, evaporation and imbibition dynamics of a sessile drop on a porous substrate. Finally, we conclude in section V with a brief discussion.

## II. MESOSCOPIC HYDRODYNAMIC MODEL

## A. Three-Field Gradient Dynamics

FIG. 1. Sketch of the considered geometry for a volatile liquid drop on a porous substrate within a shallow chamber of height d . The porous substrate has mostly vertical pores with some horizontal connections allowing for (slow) lateral diffusion of liquid. The drop profile is described by the height h ( x , t ), the substrate has a constant thickness z 0 , and the liquid contained in the pores is described by an effective height ζ ( x , t ). The particle densities of vapor ρ vap ( x , t ) and ambient air ρ air ( x , t ) together account for a constant total density ρ tot in the gas phase.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00179v1-figure-e6298726-f1cb-4505-a784-7146a2ede42e.jpg)

The gradient dynamics model that we develop in the following includes fields for the local vapor concentration above the substrate and the local amount of liquid within the porous medium. A sketch of the considered geometry is shown

in Fig. 1. Importantly, we assume that both, the gas layer and the porous substrate, are thin as compared to lateral extensions such that the local vapor concentration, as well as the local filling ratio of the pores, are approximately homogeneous in vertical direction. For the vapor layer this is ensured by considering a confined geometry, i.e., the gas phase is limited by a horizontal solid plate situated at a distance d from the substrate. Such a setup for the description of the diffuse vapor dynamics was introduced for a volatile drop on a smooth solid substrate by Hartmann et al. [68] and employed for such a drop on a brush-covered substrate in Ref. [55].

Given these assumptions, we use three variables ψ i ( x , t ) to describe the state of the system formed by the drop, porous substrate, and vapor layer. The general gradient dynamics framework reads [67]

$$\partial _ { t } \psi _ { i } = \nabla \cdot \left [ \sum _ { b = 1 } ^ { 3 } Q _ { i j } \, \nabla \frac { \delta \mathcal { F } } { \delta \psi _ { j } } \right ] - \sum _ { j = 1 } ^ { 3 } M _ { i j } \, \frac { \delta \mathcal { F } } { \delta \psi _ { j } }$$

where the subscripts i, j = 1 , 2 , 3 refer to the three fields and the 3 × 3 mobility matrices with components M ij and Q ij represent the conserved and nonconserved dynamics, respectively.

The variable ψ 1 describes the number of liquid particles within the drop per substrate area and relates to its profile height h as

$$\psi _ { 1 } ( x, t ) = \rho _ { \text{liq} } h ( x, t ).$$

where ρ liq is the constant and uniform particle number density of the liquid. Further, ψ 3 describes the number of vapor particles within the gas layer per substrate area and relates to the height-averaged vapor density ρ vap ( x , t ) via

$$\psi _ { 3 } ( x, t ) = \rho _ { \text{vap} } ( x, t ) [ d - h ( x, t ) ].$$

For simplicity, we consider the system as isothermal and assume that the total particle density ρ tot in the gas layer is constant and uniform. This assumption reflects that the velocity of sound is much larger than the velocity of all dynamics processes we are interested in.

Here, ρ tot is the sum of the vapor particles and of the remaining particles of the atmosphere constituents, simply labeled as 'air'. Therefore,

$$\rho _ { t o t } = \rho _ { \text{vap} } ( x, t ) + \rho _ { \text{air} } ( x, t ).$$

Finally, the variable ψ 2 ( x , t ) is the number of liquid particles within the porous layer per substrate area. It is proportional to the local pore filling ratio Θ( x , t ), i.e.,

$$\psi _ { 2 } ( x, t ) = \rho _ { l i q } z _ { 0 } b \, \Theta ( x, t ) = \rho _ { l i q } \zeta ( x, t )$$

where b is the uniform and constant porosity (void fraction) and z 0 is the thickness of the porous layer. We have also introduced ζ ( x , t ) = z b 0 Θ( x , t ) as the effective height of the substrate-contained liquid.

## B. Transport Processes

The conserved mobilities Q account for three processes: (i) viscous motion within the drop, (ii) Darcy-type diffusive transport of liquid particles within the substrate, and (iii) diffusive transport of vapor particles within the gas phase. This automatically neglects any dynamic coupling between the phases. In particular, there is no direct dynamic coupling between the liquid in the drop and the liquid in the porous layer, i.e., Q 12 = Q 21 = 0. This should be a good approximation at small porosity, but could in principle be improved along the lines of Thiele et al. [36], i.e., by considering Darcy-Brinkman dynamics in the porous layer. The resulting mobility matrix is diagonal, namely,

$$\mathbf Q = \begin{pmatrix} \frac { 1 } { \rho _ { 1 1 q } } \frac { \psi _ { 1 } ^ { 3 } } { 3 \eta } & 0 & 0 \\ 0 & \hat { D } _ { \text{port} } \psi _ { 2 } & 0 \\ 0 & 0 & \frac { 1 } { k _ { B T } } D _ { \text{avp} } \psi _ { 3 } \end{pmatrix}$$

with the vapor diffusion coefficient D vap . The mobility of the liquid within the porous substrate, ˆ D por , needs to be determined, e.g., from Darcy's equation [4] and will depend on the permeability κ of and the effective viscosity η eff = η/b in the porous layer. Following the phase-field approach of Cueto-Felgueroso and Juanes [25], we obtain

$$\hat { D } _ { \text{por} } = z _ { 0 } \kappa k _ { r } ( \Theta ) / \eta _ { \text{eff} },$$

where k r (Θ) is a relative permeability, i.e., an increasing convex function of the liquid saturation Θ. For an example of such a function see Fig. 2 of Ref. [25]. Darcy's law is recovered for k r (Θ) = const.

The nonconserved mobilities M account for all three transfer processes, namely, (i) evaporation (or condensation) between the film and the gas phase, (ii) evaporation (or condensation) directly between the substrate and the gas phase, and (iii) liquid transfer between the porous substrate and the drop. We assume that these processes are solely driven by the partial pressures of the phases and assign corresponding Onsager coefficients. The resulting matrix reads

$$\begin{pmatrix} M _ { \text{im} } + M _ { \text{ev} } & - M _ { \text{im} } & - M _ { \text{ev} } \\ - M _ { \text{im} } & M _ { \text{im} } + M _ { \text{ev} } ^ { \prime } & - M _ { \text{ev} } ^ { \prime } \\ - M _ { \text{ev} } & - M _ { \text{ev} } ^ { \prime } & M _ { \text{ev} } + M _ { \text{ev} } ^ { \prime } \end{pmatrix}. \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad \\.. \ \sim & \quad \sim.. \ \neq \. \ \neq \. \ \neq \. \ \neq \. \ \neq \. \ \neq \. \ \neq \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \$$

where M ev and M ′ ev are the Onsager coefficient of phase change for the drop and the substrate, respectively, and M im is an imbibition rate coefficient. The particular form of the matrix ensures that ψ 1 + ψ 2 + ψ 3 is conserved, i.e., ∂ t ( ψ 1 + ψ 2 + ψ 3 ) = -∇· j ψ 1 + ψ 2 + ψ 3 .

As the porous medium exposes only a fraction of the contained liquid to the surface via the pores, we ensure that the evaporation rate from the substrate surface scales with the porosity employing the simple assumption

$$M _ { e v } ^ { \prime } & = b M _ { e v }.$$

Keeping the overall structure of the matrix, in principle, more complicated dependencies may be incorporated into the M 's.

## C. Free Energy Functional

Next, we discuss the underlying free energy functional F [ ψ , ψ , ψ 1 2 3 ]. We employ the long-wave approximation [61, 69] and, for convenience, use the fields h , ρ vap , and ζ as abbreviations for the relations containing only ψ i . We have

$$\mathcal { F } & = \int _ { \Omega } \left [ \, \underbrace { \frac { \gamma _ { \text{g} } } { 2 } ( \nabla h ) ^ { 2 } } _ { \text{liquid-gas interface energy} } + \underbrace { \gamma _ { \text{sl} } ( \Theta ) } _ { \text{liquid-substrate interface energy} } + \underbrace { f _ { \text{wet} } ( h, \Theta ) } _ { \text{wetting energy} } + \underbrace { ( h + \zeta ) f _ { \text{liq} } ( \rho _ { \text{liq} } ) } _ { \text{liquid bulk energy} } \\ & \, \rightharpoonup \, ( d - h \, \hat { f } \, \left ( \alpha \, \right ) \, \rightharpoonup \, ( d - h \, \hat { f } \, \left ( \alpha \, \right ) \, \rightharpoonup \, \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde } \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde                                                                                                                                                                                                         }                                                                                                                                                                                                        ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }. } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }, } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ( } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } | } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } & } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } [ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \\ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }] } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }$$

$$+ \underbrace { ( d - h ) f _ { \varphi } ( \rho _ { \varphi } ) } _ { \text{var} \text{ energy} } + \underbrace { ( d - h ) f _ { \text{air} } ( \rho _ { \text{air} } ) } _ { \text{air energy} } + \underbrace { z _ { 0 } \frac { \Gamma } { 2 } ( \nabla \Theta ) ^ { 2 } + z _ { 0 } \Psi ( \Theta ) } _ { \text{substrate energy} } \right ] d x d y,$$

where f liq , f vap , and f air are bulk liquid, vapor, and air energies per volume, γ lg is the constant liquid-gas interface energy, γ sl is the (adaptive) filling ratio-dependent liquid-substrate interface energy, and f wet ( h, Θ) is the (per area) wetting energy.

The remaining terms (last two) represent the energetic description of the porous substrate layer that we base on the work of Cueto-Felgueroso and Juanes [25]. Namely, the substrate energy contains the capillary potential Ψ(Θ) and a nonlocal term that encodes the apparent interface energy Γ of imbibition fronts. The factor z 0 transforms the per-volume contributions into a per-area contribution in our height-averaged setting. It also relates the filling ratio Θ to the effective height of liquid in the pores, ζ = z b 0 Θ, see Eq. (6). The capillary potential is related to the capillary pressure function J (Θ) ∼ -∂ Θ Ψ. Analytical expressions for J and the relative permeability k r (Θ) are given, e.g., by the van Genuchten-Mualem model (vGM) [25, 70], where

$$k _ { r } & = \sqrt { \Theta } \left [ 1 - \left ( 1 - \Theta ^ { 1 / m } \right ) ^ { m } \right ] ^ { 2 }, \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

$$J = p _ { c } \, \left ( \overset { \mu } { \Theta } ^ { - 1 / m } - 1 \right ) ^ { 1 / n }.$$

Here, p c is the scale of the capillary pressure, n &gt; 1 is a measure for pore size distribution, and m = 1 -1 /n . The capillary pressure constant p c results from the Laplace pressure inside the pores and relates to the capillary rise h cap via Jurin's law [25, 71], i.e.,

$$p _ { c } = \frac { 2 \gamma _ { l g } \cos \theta _ { \text{eq} } } { r _ { 0 } } = \varrho _ { \text{liq} } \, g \, h _ { \text{cap} }$$

FIG. 2. Capillary pressure ∂ Θ Ψ and corresponding potential Ψ as functions of the pore filling ratio Θ according to Eq. (15) in units of f 0 = ρ liq k B T . We consider a low porosity b = 0 01, all other parameters are as in Tab. I, i.e., the dimensionless . capillary pressure constant is p /f c 0 ≈ 0 1. . Shown are the cases with (solid line) and without (dashed line) the entropic correction f entropic , which ensures that the capillary pressure diverges when the pores are completely filled, Θ → 1. Note that pressure and potential have identical units as Θ is dimensionless.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00179v1-figure-74b769ee-ea41-4d01-9e42-2d5e0b41a090.jpg)

with the typical pore radius r 0 ≈ √ κ/b , the liquid mass density ϱ liq , and the equilibrium contact angle θ eq of the liquid on the pure (non-porous) substrate material.

Note that the capillary pressure ∂ Θ Ψ is always negative for 0 &lt; Θ &lt; 1, zero for Θ = 1, and undefined for Θ &gt; 1. In other words, the energetic minimum of the capillary potential is at full saturation Θ = 1. The case Θ &gt; 1 corresponds to the (unphysical) oversaturation of the porous medium with liquid, i.e., the volume of the contained liquid would be larger than the space available in the pores. To prevent such unphysical behavior in our dynamical model, we additionally incorporate and entropic contribution arising for the two-phase system formed by the liquid and the gas within the pores [72]. For the nearly saturated liquid of interest, this gives the entropic contribution to the local free energy

$$f _ { \text{entropic} } ( \Theta ) = b \, \rho _ { \text{liq} } \, k _ { B } T \, ( 1 - \Theta ) \log ( 1 - \Theta ).$$

The complementary entropic contribution of the gas in the pores is negligible against the capillary potential that is part of the employed vGM model (12). We combine the entropic contribution (14) with the previously employed vGM expression (12) into the capillary potential Ψ. Here, we use this to directly define the pressure

$$\partial _ { \Theta } \Psi = - J + \partial _ { \Theta } f _ { \text{entropic} }.$$

At low porosity b , the entropic part only corresponds to a small correction to the capillary pressure (and resulting potential), as can be assessed in Fig. 2. There both the capillary pressure ∂ Θ Ψ and potential Ψ are illustrated. Nevertheless, including the entropy is crucial as it ensures that the capillary pressure diverges when approaching complete filling, Θ → 1, i.e., the minimum of the capillary potential is always at Θ &lt; 1.

To account for the adaptive wettability of the porous substrate we follow an approach similar to Kap et al. [55]. Namely, we assume a power law dependence of the Hamaker constant on the volume fraction c of the substrate material in the porous layer.

$$f _ { \text{wet} } ( h, \, \Theta ) = A _ { 0 } \, c ^ { a } \left ( \frac { h _ { p } ^ { 3 } } { 5 h ^ { 5 } } - \frac { 1 } { 2 h ^ { 2 } } \right ),$$

where A 0 is the Hamaker constant of a poreless substrate, and a ≥ 0 is some constant exponent. The substrate volume fraction is c = 1 -b Θ, i.e., the complement to the liquid concentration.

Furthermore, we assume that the interface energy between the drop and the porous layer adapts with the same power law as used in the wetting potential. This ensures consistency with the adaption of a macroscopic substrate-gas interface energy, similar to the case of a polymer-brush substrate discussed in Ref. [73]. Hence, we write for the substrate-liquid interfacial energy

$$\gamma _ { s l } ( \Theta ) = \gamma _ { s l, \, \text{dry} } \, c ^ { a }$$

with some positive exponent a . This implies that the interface energy decreases when the substrate fills with liquid.

Vapor and air are considered to correspond to ideal gases, cf. Ref. [68]. In consequence, their respective free energy densities are purely entropic

$$f _ { \varphi } = k _ { B } T \rho _ { \varphi } \left [ \log ( \Lambda ^ { 3 } \rho _ { \varphi } ) - 1 \right ] \text{ and } f _ { \text{air} } = k _ { B } T \rho _ { \text{air} } \left [ \log ( \Lambda ^ { 3 } \rho _ { \text{air} } ) - 1 \right ]$$

with the mean free path length Λ. This completes the discussion of all modeling elements.

## D. The kinetic model

The dynamic equations (2) involve the variations of the free energy functional, which compute as

$$\text{max equations} \left ( \varepsilon \text{ in vivue } \text{various on one one energy} \text{ quanto}, \text{mean compute as} \\ \frac { \delta \mathcal { F } } { \delta h } = \rho _ { \text{liq} } \frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } = - \gamma _ { \text{lg} } \Delta h + \partial _ { h } f _ { \text{wet} } ( h, \Theta ) + f _ { \text{liq} } \\ \frac { \delta \mathcal { F } } { \delta \zeta } = \rho _ { \text{liq} } \frac { \delta \mathcal { F } } { \delta \psi _ { 2 } } = \partial _ { \zeta } \left [ \gamma _ { \text{sl} } ( \Theta ) + f _ { \text{wet} } ( h, \Theta ) + z _ { 0 } \Psi ( \Theta ) \right ] + f _ { \text{liq} } - \frac { \Gamma } { z _ { 0 } b } \Delta \Theta \\ \frac { \delta \mathcal { F } } { \delta \psi _ { 3 } } = k _ { B } T \log \left ( \frac { \rho _ { \text{vap} } } { \rho _ { \text{tot} } - \rho _ { \text{vap} } } \right ), \\ \text{the final equation we have used that the vapor particle density is much smaller than the total gas particle}$$

where in the final equation we have used that the vapor particle density is much smaller than the total gas particle density, which itself is much smaller than the liquid density ρ vap ≪ ρ tot ≪ ρ liq , as discussed in more detail by Hartmann et al. [68].

Inserting the energy variations as well as the mobility matrices (Eqs. (7) and (9)) into the three-field gradient dynamics Eq. (2) gives the resulting dynamic equations

$$\circ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \partial _ { t } \psi _ { 1 } / \rho _ { \text{liq} } & = \nabla \cdot \left [ \frac { \psi _ { 3 } ^ { 3 } } { 3 \eta \rho _ { \text{liq} } ^ { 2 } } \nabla \frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } \right ] - j _ { \text{ev} } - j _ { \text{im} }, \\ \partial _ { t } \psi _ { 2 } / \rho _ { \text{liq} } & = \nabla \cdot \left [ \frac { D _ { \text{por} } } { \rho _ { \text{liq} } } \, \psi _ { 2 } \, \nabla \frac { \delta \mathcal { F } } { \delta \psi _ { 2 } } \right ] - j _ { \text{ev} } ^ { \prime } + j _ { \text{im} }, \\ \partial _ { t } \psi _ { 3 } / \rho _ { \text{liq} } & = \nabla \cdot [ D _ { \text{vap} } ( d - h ) \nabla \rho _ { \text{vap} } ] / \rho _ { \text{liq} } + j _ { \text{ev} } + j _ { \text{ev} } ^ { \prime }. \\ \text{s are}$$

The transfer fluxes are

$$j _ { e v } = \frac { M _ { e v } } { \rho _ { \text{Iiq} } } \left ( \frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } - \frac { \delta \mathcal { F } } { \delta \psi _ { 3 } } \right ),$$

$$j _ { e v } ^ { \prime } = \frac { \stackrel { \cdot \cdot \cdot \cdot } { M _ { e v } ^ { \prime } } } { \rho _ { \text{Liq} } } \left ( \frac { \delta \mathcal { F } } { \delta \psi _ { 2 } } - \frac { \delta \mathcal { F } } { \delta \psi _ { 3 } } \right ),$$

$$j _ { \lim } = \frac { \dot { M _ { \lim } } } { \rho _ { \text{liq} } } \left ( \frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } - \frac { \delta \mathcal { F } } { \delta \psi _ { 2 } } \right ),$$

namely, the film evaporation/condensation flux, the evaporation/condensation flux of the liquid contained in the porous layer, and the imbibition flux, respectively.

Finally, for an efficient numerical treatment, we reformulate the governing equation in terms of the dimensionless relative vapor concentration (or relative humidity) ϕ ( x , t ) = k B Tρ vap ( x , t ) /p sat and the dimensionless filling ratio of the substrate Θ( x , t ) = ζ/ z b ( 0 ).

$$\langle \varkappa, \varkappa / - s / \wedge \varkappa \rangle. \\ \partial _ { t } h & = \nabla \left [ \frac { h ^ { 3 } } { 3 \eta } \nabla \frac { \delta \mathcal { F } } { \delta h } \right ] - j _ { e$$

This completes the model development. The model is now employed to discuss in the next section the resulting sorption isotherm of the porous medium. The subsequent section IV analyzes the coupled spreading, evaporation and imbibition dynamics of a sessile drop on a porous substrate.

## III. SORPTION ISOTHERM OF A POROUS MEDIUM

Before we employ numerical tools to study the behavior of the coupled three-phase model, we consider spatially homogeneous equilibrium states, where h ( x , t ) = const, Θ( x , t ) = const, and ϕ ( x , t ) = const. These states can be assessed analytically, as the governing equations strongly simplify when all spatial and temporal derivatives are set to zero. In addition, we assume that the liquid film is either very thick or that it corresponds to a thin precursor layer, i.e., h ≫ h p or h = h p , respectively. In both cases the wetting potential is negligible, f wet ( h, Θ) ≈ 0.

Evaluating the partial pressures (Eq. (19)) for a homogeneous equilibrium state yields

$$\frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } & = \frac { 1 } { \rho _ { \text{liq} } } f _ { \text{liq} } \\ \frac { \delta \mathcal { F } } { \delta \psi _ { 2 } } & = \frac { 1 }$$

In thermodynamic equilibrium, the system is at rest. This implies that the transfer fluxes j ev , j ′ ev , and j im are zero as the partial pressures balance. In contrast, a larger pressure in the film than in the substrate causes an absorption of liquid into the substrate until the pressures equilibrate. Their balance results in coexistence conditions for the three phases.

When a thick homogeneous film is in equilibrium with an ambient vapor phase, the air is considered saturated with vapor. Hence, comparing the film and vapor partial pressures defines the saturation vapor density, i.e.,

$$\frac { \delta \mathcal { F } } { \delta \psi _ { 1 } } = \frac { \delta \mathcal { F } } { \delta \psi _ { 3 } } \Big | _ { \rho _ { \text{vag} } = \rho _ { \text{sat} } } \quad \quad \quad$$

$$\Leftrightarrow \quad \frac { \stackrel {.. } { f _ { \text{liq} } } } { \rho _ { \text{liq} } } = k _ { B } T \log \left ( \frac { \rho _ { \text{sat} } } { \rho _ { \text{tot} } - \rho _ {$$

Below, this result is used to relate the liquid free energy f liq to the vapor saturation pressure. Note that it is independent of the substrate state.

Similarly, equating the film and substrate pressures, δ F /δψ 1 = δ F /δψ 2 , allows one to define the equilibrium filling ratio for a substrate that is covered with a thick film

$$\partial _ { \zeta } \left [ \gamma _ { \text{sl} } ( \Theta ) + f _ { \text{wet} } ( h, \Theta ) + z _ { 0 } \Psi ( \Theta ) \right ] = 0.$$

Next, equating the substrate and vapor pressures, δ F /δψ 2 = δ F /δψ 3 , gives the equilibrium substrate filling ratio as a function of the ambient vapor density.

$$\partial _ { \zeta } \left [ \gamma _ { \text{sl} } ( \Theta ) + f _ { \text{wet} } ( h _ { p }, \Theta ) + z _ { 0 } \Psi ( \Theta ) \right ] + f _ { \text{liq} } = \rho _ { \text{liq} } k _ { B } T \log \left ( \frac { \rho _ { \text{vap} } } { \rho _ { \text{tot} } - \rho _ { \text{vap} } } \right ).$$

Using Eq. (27) to express the liquid bulk free energy through the saturation pressure, the equation can be written as

$$\phi ( \Theta ) = \exp \left \{ \frac { \partial _ { \zeta } \left [ \gamma _ { s l } ( \Theta ) + f _ { \text{w} } ( h _ { p }, \Theta ) + z _ { 0 } \Psi ( \Theta ) \right ] } { \rho _ { \text{liq} } k _ { B } T } \right \}.$$

The result implicitly defines the sorption isotherm Θ( ϕ ), i.e., the substrate state as a function of the vapor concentration. If the porosity is low (or the substrate is sufficiently thick), the term related to the (bulk) capillary potential dominates the interface terms. This is implied by their scaling behavior with respect to the substrate thickness z 0 or the porosity b . Thus, we simplify the definition of the isotherm by neglecting the interface terms and write

$$\phi ( \Theta ) = \exp \left \{ \frac { \partial _ { \Theta } \Psi ( \Theta ) } { b \rho _ { \text{liq} } k _ { B } T } \right \}.$$

Fig. 3 shows sorption isotherms Θ( ϕ ) of a porous substrate at fixed dimensionless capillary pressure constant p / ρ c ( liq k B T ) = 0 1 for various values of the porosity . b and the exponent n .

FIG. 3. Sorption isotherms Θ( ϕ ) of a porous medium in vapor according to Eq. (31) for various parameter configurations and fixed dimensionless capillary pressure constant p / ρ c ( liq k B T ) = 0 1. . In panel (a) we linearly vary the porosity b at fixed n = 5 and panel (b) shows the linear variation of the exponent n at fixed porosity b = 0 2. . Note that due to entropic influences some air remains in the pores even for a saturated gas phase, Θ( ϕ = 100%) &lt; 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00179v1-figure-966604e7-299e-473e-bf26-4559bc2955f1.jpg)

TABLE I. Model parameters used for the numerical simulation of the coupled spreading, evaporation and imbibition dynamics in Fig. 4.

| Parameter name                                    | Symbol           |                                                                       |
|---------------------------------------------------|------------------|-----------------------------------------------------------------------|
| viscosity contact angle parameter (dry            | η θ eq h p ρ liq | Value 0 . 89 mPas tan 30 ◦ 10 µ m 997 kg / m 3 18 g / mol N 2643 Pa ◦ |
| substrate)                                        |                  |                                                                       |
| precursor layer height                            |                  |                                                                       |
| liquid particle density                           |                  | A                                                                     |
| vapor saturation pressure                         | p sat            |                                                                       |
| temperature                                       | T                | 22 C                                                                  |
| initial drop volume                               | V 0              | 23 µ l                                                                |
| initial vapor concentration                       | ϕ lab            | 10 %                                                                  |
| liquid-gas interface energy                       | γ lg             | 72 . 8 mN / m                                                         |
| substrate-liquid interface energy (dry substrate) | γ bl , 0         | 3 mN / m                                                              |
| imbibition front interface energy                 | Γ                | 0 . 3 mN / m                                                          |
| substrate porosity                                | b                | 0 . 2                                                                 |
| substrate thickness                               | z 0              | 10 µ m                                                                |
| pore size distribution coefficient                | n                | 5                                                                     |
| adaption exponent (power law)                     | a                | 1                                                                     |
| permeability                                      | κ                | 2 × 10 - 17 m 2                                                       |
| vapor diffusion coefficient                       | D vap            | 2 . 82 × 10 - 5 m 2 / s                                               |
| imbibition rate coefficient                       | M im             | 10 - 13 m / Pas                                                       |
| bulk liquid evaporation rate coefficient          | M ev             | 10 - 13 m / Pas                                                       |
| substrate evaporation rate coefficient            | M ′ ev           | 10 - 13 m / Pas                                                       |
| simulation domain height                          | d                | 3 mm                                                                  |
| simulation domain width                           | L                | 20 mm                                                                 |

## IV. SPREADING DYNAMICS: HALO FORMATION

Next we consider an example for the coupled spreading, evaporation and imbibition dynamics of a sessile drop on a porous substrate. Fig. 4 shows four corresponding snapshots of a simulation of the developed three-field porous substrate model, Eq. (24), in a radial geometry. The boundary conditions are chosen such that the vapor particles are allowed to freely escape from the radial domain at the outer boundary (at r = 20mm), where the lab humidity ϕ lab = 10% is enforced. The full set of parameters is given in Tab. I. The parameters represent realistic values, however, our choice is not geared towards a particular experiment. The simulation is performed employing the finite element method implemented in the package oomph-lib [74]. Details regarding the model transformed into radial geometry can be found in Ref. 75.

At each of the four chosen times, Fig. 4 contains three panels: The top panel shows the three transfer fluxes between the phases, namely, liquid-to-vapor evaporation/condensation flux j ev , substrate-to-vapor evaporation/condensation flux j ′ ev , and liquid-to-substrate imbibition flux j im . The central panel visualizes the state of the system with the

FIG. 4. Snapshots at four selected times from a numerical simulation of the coupled spreading, evaporation and imbibition dynamics of a sessile drop on a porous substrate at parameters given in Tab. I. At each point in time (a)-(d) three subpanels are given: (top) local transfer rates j ev , j ev ', and j im on a symmetric log-scale with a linear scale between ± 10 -5 µ m s (volume / per area per time); (center) simulation domain with the drop in dark blue, the vapor concentration ϕ as blue shading, and the porous substrate filling state Θ as red-to-yellow shading; (bottom) relative vapor and substrate humidities, ϕ and Θ. For chosen times and various line styles see legends.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00179v1-figure-aff16e93-6971-4a1f-bb01-d58190cc9bfc.jpg)

drop given in dark blue, the ambient vapor concentration presented as a light blue shading, and the porous substrate state as a red-to-yellow shade below the droplet. The shading represents the local filling ratio, where darker colors correspond to a larger liquid content of the substrate. The bottom panel gives the humidity and filling ratio profiles of gas layer and porous substrate layer, respectively.

Inspection of Fig. 4 indicates that in the initial second, fast imbibition of liquid occurs underneath the drop (green line in bottom panels of Fig. 4(a) and (b)). In parallel, evaporation results in saturation of the vapor phase directly above the droplet that here has the properties of water (red line in bottom panels of Fig. 4(a) and (b)). Interestingly, there is some absorption of liquid from the vapor into the porous layer outside the drop (green line in top panels of Fig. 4(a) and (b)). Over the next minutes (Fig. 4(c) and (d)), lateral diffusion becomes the dominating process. On the one hand, the vapor diffuses what represent the limiting process for the evaporation of the drop as we are in the diffusion-limited regime, see Ref. [68] for an extensive discussion of the different regimes. On the other hand, also the substrate filling is extended far beyond the contact line of the drop. This occurs not only due to diffusion within the substrate but also because there is a significant condensation of liquid from the vapor phase into the substrate. We emphasize that the strength of this effect depends on the ratio between the various rates for the competing diffusive processes. For instance, if the substrate had a larger permeability, transport within the substrate could be faster than transport in the gas phase implying that one could in the first minutes observe evaporation of liquid from the

substrate outside the drop into the gas phase.

Here, the condensation continues for about a minute (between Fig. 4(c) and (d)). Then, the process reverses and liquid evaporates from the porous layer into the gas phase. This results in the emergence of a quasi-stationary macroscopic 'halo' in the substrate filling ratio. The system establishes a continuous flow of liquid from the drop radially outwards through the porous substrate. As it then evaporates thereby accelerating the transport through the gas phase the halo develops. The resulting halo persists for the full time until the droplet has completely evaporated. Similar results have been obtained for the case of a liquid drop evaporating on a polymer-brush covered substrate [55]. However, there the brush does not only adapt its 'filling ratio' but also its surface topography.

## V. CONCLUSION

We have presented a mesoscopic hydrodynamic model in gradient dynamics formulation for a sessile drop of volatile partially wetting liquid on a solid porous layer that is able to describe situations where several dynamic processes compete that are all driven by the same underlying energy functional. The latter incorporates saturation-dependent wettability, capillarity of the drop and within the porous layer, air and vapor entropy and an entropic contribution relevant close to complete filling of the porous layer that avoids unphysical behavior. The modeled transport and transfer processes are the spreading of the droplet, the evaporation of liquid from the drop into the ambient gas or condensation of vapor from the ambient gas, vapor diffusion in the gas phase, the vertical absorption of liquid from the drop or the vapor into the porous substrate, and the lateral imbibition dynamics of the liquid within the porous substrate layer. These coupled processes are captured by spatially two-dimensional partial differential equations for the particle densities per substrate area in the three regions: liquid drop, porous layer and gas phase. These particle densities are closely related to the drop height profile, the vertically averaged saturation profile in the porous layer, and the vertically averaged vapor density profile, respectively. The approach assumes strong confinement, i.e., we have (i) considered a sessile drop on one of two parallel plates that form a narrow gap, and (ii) a thin porous substrate layer. This has generalized the approach of Ref. [68] for drops of volatile liquids within a gap between smooth solid plates, and partially parallels such models for drops on polymer brushes [56]. In contrast, here, the imbibition process is incorporated adapting and extending models for the dynamics of the saturation profile in a porous medium that represent diffuse-interface generalizations of Darcy's law and Richards' equation used, e.g., to study gravity-driven saturation fronts within porous media [24, 25, 43].

The incorporation of drop and vapor dynamics now allows for the adaptation of the extended model to a number of practical problems. An example important in the context of rhizospheric soil in the presence of plant root-produced exopolysaccharides is the water transfer between a drop and an amphiphilic porous medium, studied for nonvolatile liquids in [76, 77]. There, one could introduce the interplay between ambient humidity and saturation dynamics. Other examples that one could now investigate with consideration of both, imbibition and vapor dynamics, include inkjet printing on layers of paper [78] and on microporous membranes [79], and sliding drops on porous substrates [80]. Furthermore, the approach could be expanded to cover various suspensions and solutions on porous substrates [81, 82] by combining the gradient dynamics approach presented here with elements of Refs. [64, 65].

In the present work we have given the mesoscopic hydrodynamic model in long-wave approximation [61]. However, as discussed in Ref. [67] in the spirit of the Cahn-Hilliard equation for phase separation one may also use a hybrid model where the mobilities are kept in long-wave approximation but the energy functional is made more precise by using a full-curvature formulation. It has been shown that this often results in more exact models, see e.g., Refs. [56, 83-85]. Here, the full-curvature formulation is obtained by replacing the long-wave metric factor 1 2 ( ∇ h ) 2 in Eq. (11) by the full expression √ 1 + ( ∇ h ) 2 . Note that the basic gradient dynamics can, in principle, also be adapted to situations like substrates with chemical or topographic heterogeneities [86, 87]. The employed gradient dynamics approach also provides a relatively simple, but by definition thermodynamically consistent approach to more complex situations [56, 65-67, 88]. In this way, our present approach could in the future be further expanded to cover situations where the underlying governing equations are not yet completely known. Examples include systems where imbibition results in the swelling of the porous medium as for thin deformable porous media [89], liquid-infused slippery substrates that are often based on porous substrate layers [48-51], imbibition into loose material as a powder bed [90], spreading and imbibition of drops of liquid mixtures, and photocatalytically active lubricant-impregnated surfaces [91].

## DATA AVAILABILITY STATEMENT

The underlying data needed to reproduce shown results are publicly available at the data repository zenodo , see: Hartmann, S. and Thiele, U. (2024). Data Supplement for: 'Gradient dynamics model for drops of volatile liquid on

## ACKNOWLEDGMENTS

The authors acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) via grant No. TH781/12-1 and TH781/12-2 within Priority Program (SPP) 2171 'Dynamic Wetting of Flexible, Adaptive, and Switchable Surfaces'; UT would further like to thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and hospitality during the program Anti-Diffusion in Multiphase and Active Flows (ADIW04) where some work on this paper was undertaken. The program was supported by EPSRC grant No. EP/R014604/1. We acknowledge fruitful discussions with Florian Cajot, Philippe Beltrame, Laura Gallardo Dominguez, Patrick Huber as well as many discussions on adaptive substrates with Jan Diekmann, Daniel Greve, and Christopher Henkel at the University of M¨ unster, and with participants of the events organized by SPP 2171.

- [1] L. A. Richards, Capillary conduction of liquids through porous mediums, Physics 1 , 318 (1931).
- [2] P. G. Saffman and G. Taylor, The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid, Proc. R. Soc. London Ser. A 245 , 312 (1958).
- [3] G. S. Beavers and D. D. Joseph, Boundary conditions at a naturally permeable wall, J. Fluid Mech. 30 , 197 (1967).
- [4] S. Whitaker, Flow in porous media I: A theoretical derivation of Darcy's law, Transp. Porous Media 1 , 3 (1986).
- [5] A. Marmur, Drop penetration into a thin porous-medium, J. Colloid Interface Sci. 123 , 161 (1988).
- [6] S. H. Davis and L. M. Hocking, Spreading and imbibition of viscous liquid on a porous base, Phys. Fluids 11 , 48 (1999).
- [7] S. I. Vasin, V. M. Starov, and A. N. Filippov, The motion of a solid spherical particle covered with a porous layer in a liquid, Colloid J. 58 , 282 (1996).
- [8] M. G. Worster, Instabilities of the liquid and mushy regions during solidification of alloys, J. Fluid Mech. 237 , 649 (1992).
- [9] S. H. Davis and L. M. Hocking, Spreading and imbibition of viscous liquid on a porous base. II, Phys. Fluids 12 , 1646 (2000).
- [10] V. M. Starov, S. A. Zhdanov, S. R. Kosvintsev, V. D. Sobolev, and M. G. Velarde, Spreading of liquid drops over porous substrates, Adv. Colloid Interface Sci. 104 , 123 (2003).
- [11] N. Alleborn and H. Raszillier, Spreading and sorption of a droplet on a porous substrate, Chem. Eng. Sci. 59 , 2071 (2004).
- [12] J. F. Xiao, H. A. Stone, and D. Attinger, Source-like solution for radial imbibition into a homogeneous semi-infinite porous medium, Langmuir 28 , 4208 (2012).
- [13] T. Gambaryan-Roisman, Liquids on porous layers: wetting, imbibition and transport processes, Curr. Opin. Colloid Interface Sci. 19 , 320 (2014).
- [14] O. Arjmandi-Tash, N. M. Koyalchuk, A. Trybala, I. V. Kuchin, and V. Starov, Kinetics of wetting and spreading of droplets over various substrates, Langmuir 33 , 4367 (2017).
- [15] Y. Yao, S. Alexandris, F. Henrich, G. Auernhammer, M. Steinhart, H. J. Butt, and G. Floudas, Complex dynamics of capillary imbibition of poly(ethylene oxide) melts in nanoporous alumina, J. Chem. Phys. 146 , 203320 (2017).
- [16] P. Johnson, A. Trybala, and V. Starov, Kinetics of spreading over porous substrates, Colloid Interfac. 3 , 38 (2019).
- [17] S. A. Etha, P. R. Desai, H. S. Sachar, and S. Das, Wetting dynamics on solvophilic, soft, porous, and responsive surfaces, Macromolecules 54 , 584 (2021).
- [18] P.-G. de Gennes, Running droplets in a random medium, C. R. Acad. Sci. II 327 , 147 (1999).
- [19] A. Aradian, E. Rapha¨ el, and P. G. de Gennes, Dewetting on porous media with aspiration, Eur. Phys. J. E 2 , 367 (2000).
- [20] L. Bacri and F. Brochard-Wyart, Dewetting on porous media, Europhys. Lett. 56 , 414 (2001).
- [21] T. Desaive, G. Lebon, and M. Hennenberg, Coupled capillary and gravity-driven instability in a liquid film overlying a porous layer, Phys. Rev. E 64 , 066304 (2001).
- [22] I. M. R. Sadiq, R. Usha, and S. W. Joo, Instabilities in a liquid film flow over an inclined heated porous substrate, Chem. Eng. Sci. 65 , 4443 (2010).
- [23] L. Cueto-Felgueroso and R. Juanes, Nonlocal interface dynamics and pattern formation in gravity-driven unsaturated flow through porous media, Phys. Rev. Lett. 101 , 244504 (2008).
- [24] L. Cueto-Felgueroso and R. Juanes, A phase field model of unsaturated flow, Water Resour. Res. 45 , W10409 (2009).
- [25] L. Cueto-Felgueroso and R. Juanes, Stability analysis of a phase-field model of gravity-driven unsaturated flow through porous media, Phys. Rev. E 79 , 036301 (2009).
- [26] O. Devauchelle, C. Josserand, and S. Zaleski, Forced dewetting on porous media, J. Fluid Mech. 574 , 343 (2007).
- [27] T. Goda, Y. Sasaki, M. Mizuno, K. Morizawa, H. Katakura, and S. Tomiya, Numerical analysis for predicting the operability window of slot-die coating onto porous media, J. Coat. Technol. Res. 14 , 1053 (2017).
- [28] P. Huber, Soft matter in hard confinement: Phase transition thermodynamics, structure, texture, diffusion and flow in nanoporous media, J. Phys.: Condens. Matter 27 , 103102 (2015).
- [29] Z. B. Wang, L. Heinke, J. Jelic, M. Cakici, M. Dommaschk, R. J. Maurer, H. Oberhofer, S. Grosjean, R. Herges, S. Brase, K. Reuter, and C. Woll, Photoswitching in nanoporous, crystalline solids: an experimental and theoretical study for

- azobenzene linkers incorporated in mofs, Phys. Chem. Chem. Phys. 17 , 14582 (2015).
- [30] G. Bayada and M. Chambat, On interface conditions for a thin film flow past a porous medium, SIAM J. Math. Anal. 26 , 1113 (1995).
- [31] J. A. Ochoa-Tapia and S. Whitaker, Momentum transfer at the boundary between a porous medium and a homogeneous fluid: I. Theoretical development, Int. J. Heat Mass Transfer 38 , 2635 (1995).
- [32] J. A. Ochoa-Tapia and S. Whitaker, Momentum transfer at the boundary between a porous medium and a homogeneous fluid: II. Comparison with experiment, Int. J. Heat Mass Transfer 38 , 2647 (1995).
- [33] J. T. Jeong, Slip boundary condition on an idealized porous wall, Phys. Fluids 13 , 1884 (2001).
- [34] B. Goyeau, D. Lhuillier, D. Gobin, and M. G. Velarde, Momentum transport at a fluid-porous interface, Int. J. Heat Mass Transf. 46 , 4071 (2003).
- [35] F. J. Valdes-Parada, B. Goyeau, and J. A. Ochoa-Tapia, Jump momentum boundary condition at a fluid-porous dividing surface: Derivation of the closure problem, Chem. Eng. Sci. 62 , 4025 (2007).
- [36] U. Thiele, B. Goyeau, and M. G. Velarde, Film flow on a porous substrate, Phys. Fluids 21 , 014103 (2009).
- [37] D. J. Knox, S. K. Wilson, B. R. Duffy, and S. Mckee, Porous squeeze-film flow, IMA J. Appl. Math. 80 , 376 (2015).
- [38] P. Angot, B. Goyeau, and J. A. Ochoa-Tapia, Asymptotic modeling of transport phenomena at the interface between a fluid and a porous layer: jump conditions, Phys. Rev. E 95 , 063302 (2017).
- [39] S. D. Fitzgerald and A. W. Woods, The instability of a vaporization front in hot porous rock, Nature 367 , 450 (1994).
- [40] A. D. Angelopoulos, V. N. Paunov, V. N. Burganos, and A. C. Payatakes, Lattice Boltzmann simulation of nonideal vapor-liquid flow in porous media, Phys. Rev. E 57 , 3237 (1998).
- [41] E. N. M. Cirillo, N. Ianiro, and G. Sciarra, Phase coexistence in consolidating porous media, Phys. Rev. E 81 , 061121 (2010).
- [42] Y. J. Chen, S. Watanabe, and K. Yoshikawa, Roughening dynamics of radial imbibition in a porous medium, J. Phys. Chem. C 119 , 12508 (2015).
- [43] P. Beltrame and F. Cajot, Model of hydrophobic porous media applied to stratified media: Water trapping, intermittent flow and fingering instability, Europhys Lett 138 , 53004 (2022).
- [44] A. de Wit and G. M. Homsy, Nonlinear interactions of chemical reactions and viscous fingering in porous media, Phys. Fluids 11 , 949 (1999).
- [45] A. De Wit, Fingering of chemical fronts in porous media, Phys. Rev. Lett. 8705 , 054502 (2001).
- [46] D. D. Do, H. D. Do, and K. N. Tran, Analysis of adsorption of gases and vapors on nonporous graphitized thermal carbon black, Langmuir 19 , 5656 (2003).
- [47] M. A. Samaha and M. Gad-el Hak, Polymeric slippery coatings: Nature and applications, Polymers 6 , 1266 (2014).

[48]

J. H. Guan, G. G. Wells, B. Xu, G. McHale, D. Wood, J. Martin, and S. Stuart-Cole, Evaporation of sessile droplets on slippery liquid-infused porous surfaces (slips), Langmuir

31

,

11781 (2015).

- [49] K. Manabe, T. Matsubayashi, M. Tenjimbayashi, T. Moriya, Y. Tsuge, K. H. Kyung, and S. Shiratori, Controllable broadband optical transparency and wettability switching of temperature activated solid/liquid-infused nanofibrous membranes, ACS Nano 10 , 9387 (2016).
- [50] J. T. Luo, N. R. Geraldi, J. H. Guan, G. McHale, G. G. Wells, and Y. Q. Fu, Slippery liquid-infused porous surfaces and droplet transportation by surface acoustic waves, Phys. Rev. Appl. 7 , 014017 (2017).
- [51] Z. Brabcova, G. McHale, G. G. Wells, C. V. Brown, and M. I. Newton, Electric field induced reversible spreading of droplets into films on lubricant impregnated surfaces, Appl. Phys. Lett. 110 , 121603 (2017).
- [52] U. Thiele, Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth, J. Phys.: Condens. Matter 22 , 084019 (2010).
- [53] L. Cueto-Felgueroso and R. Juanes, Adaptive rational spectral methods for the linear stability analysis of nonlinear fourthorder problems, J. Comput. Phys. 228 , 6536 (2009).
- [54] L. Cueto-Felgueroso and R. Juanes, Reply to comment by david a. dicarlo on 'a phase field model of unsaturated flow', Water Resour. Res. 46 , W12802 (2010).
- [55] ¨ . Kap, S. Hartmann, H. Hoek, S. de Beer, I. Siretanu, U. Thiele, and F. Mugele, Nonequilibrium configurations of swelling O polymer brush layers induced by spreading drops of weakly volatile oil, J. Chem. Phys. 158 , 174903 (2023).
- [56] S. Hartmann, J. Diekmann, D. Greve, and U. Thiele, Drops on polymer brushes - Advances in thin-film modelling of adaptive substrates, Langmuir 40 , 4001 (2024).
- [57] J. W. Cahn, On spinodal decomposition, Acta Met. 9 , 795 (1961).

[58]

J. W. Cahn, Phase separation by spinodal decomposition in isotropic systems, J. Chem. Phys.

42

,

93 (1965).

- [59] A. J. Bray, Theory of phase-ordering kinetics, Adv. Phys. 43 , 357 (1994).
- [60] V. S. Mitlin, Dewetting of solid surface: Analogy with spinodal decomposition, J. Colloid Interface Sci. 156 , 491 (1993).
- [61] A. Oron, S. H. Davis, and S. G. Bankoff, Long-scale evolution of thin liquid films, Rev. Mod. Phys. 69 , 931 (1997).
- [62] A. Pototsky, M. Bestehorn, D. Merkt, and U. Thiele, Morphology changes in the evolution of liquid two-layer films, J. Chem. Phys. 122 , 224711 (2005).
- [63] S. Bommer, F. Cartellier, S. Jachalski, D. Peschka, R. Seemann, and B. Wagner, Droplets on liquids and their journey into equilibrium, Eur. Phys. J. E 36 , 87 (2013).
- [64] U. Thiele, D. V. Todorova, and H. Lopez, Gradient dynamics description for films of mixtures and suspensions: dewetting triggered by coupled film height and concentration fluctuations, Phys. Rev. Lett. 111 , 117801 (2013).
- [65] U. Thiele, A. J. Archer, and L. M. Pismen, Gradient dynamics models for liquid films with soluble surfactant, Phys. Rev. Fluids 1 , 083903 (2016).

- [66] C. Henkel, M. Essink, T. Hoang, G. van Zwieten, E. van Brummelen, U. Thiele, and J. Snoeijer, Soft wetting with (a)symmetric Shuttleworth effect, Proc. R. Soc. A 478 , 20220132 (2022).
- [67] U. Thiele, Recent advances in and future challenges for mesoscopic hydrodynamic modelling of complex wetting, Colloid Surf. A 553 , 487 (2018).
- [68] S. Hartmann, C. Diddens, M. Jalaal, and U. Thiele, Sessile drop evaporation in a gap - crossover between diffusion-limited and phase transition-limited regime, J. Fluid Mech. 960 , A32 (2023).
- [69] U. Thiele, Structure formation in thin liquid films, in Thin Films of Soft Matter , edited by S. Kalliadasis and U. Thiele (Springer Vienna, Vienna, 2007) pp. 25-93.
- [70] M. T. van Genuchten, A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J. 44 , 892 (1980).
- [71] J. Jurin, An account of some experiments shown before the Royal Society; With an enquiry into the cause of the ascent and suspension of water in capillary tubes, Phil. Trans. 30 , 739 (1718).
- [72] M. Doi, Soft matter physics (Oxford University Press, Oxford, 2013).
- [73] D. Greve, S. Hartmann, and U. Thiele, Stick-slip dynamics in the forced wetting of polymer brushes, Soft Matter 19 , 4041 (2023).
- [74] M. Heil and A. L. Hazel, Oomph-lib - an object-oriented multi-physics finite-element library, in Fluid-Structure Interaction: Modelling, Simulation, Optimisation , edited by H.-J. Bungartz and M. Sch¨fer (Springer, Berlin, Heidelberg, 2006) pp. a 19-49.
- [75] S. Hartmann, Contact Line Dynamics on Rigid and Adaptive Substrates , Ph.D. thesis, University of M¨nster, M¨nster u u (2023).
- [76] F. Cajot, Modeling water transfer in soil in the presence of amphiphilic matter: application to the rhizosphere , Ph.D. thesis, University of Avignon, Avignon (2024).
- [77] F. Cajot, C. Doussan, and P. Beltrame, Model of water transfer through an amphiphilic porous medium, Preprint (2024).
- [78] A. Clarke, T. D. Blake, K. Carruthers, and A. Woodward, Spreading and imbibition of liquid droplets on porous surfaces, Langmuir 18 , 2980 (2002).
- [79] D. Bhattacharjee, H. Nazaripoor, B. Soltannia, M. F. Ismail, and M. Sadrzadeh, An experimental and numerical study of droplet spreading and imbibition on microporous membranes, Colloid Surf. A-Physicochem. Eng. Asp. 615 , 126191 (2021).
- [80] R. Gimenez, M. Mercuri, C. L. A. Berli, and M. G. Bellino, Sliding of drops on mesoporous thin films, Phys. Chem. Chem. Phys. 22 , 5915 (2020).
- [81] V. M. Starov, Surfactant solutions and porous substrates: spreading and imbibition, Adv. Colloid Interface Sci. 111 , 3 (2004).
- [82] R. K. Holman, M. J. Cima, S. A. Uhland, and E. Sachs, Spreading and infiltration of inkjet-printed polymer solution droplets on a porous substrate, J. Colloid Interface Sci. 249 , 432 (2002).
- [83] A. von Borries Lopes, U. Thiele, and A. L. Hazel, On the multiple solutions of coating and rimming flows on rotating cylinders, J. Fluid Mech. 835 , 540 (2018).
- [84] J. Diekmann, Macroscopic and Mesoscopic Modelling of Compound Drops , Master's thesis, University of M¨nster, M¨nster u u (2022).
- [85] J. Diekmann and U. Thiele, Mesoscopic hydrodynamic model for spreading, sliding and coarsening compound drops, Preprint (2024).
- [86] U. Thiele, L. Brusch, M. Bestehorn, and M. B¨r, Modelling thin-film dewetting on structured substrates and templates: a Bifurcation analysis and numerical simulations, Eur. Phys. J. E 11 , 255 (2003).
- [87] N. Savva and S. Kalliadasis, Influence of gravity on the spreading of two-dimensional droplets over topographical substrates, J. Eng. Math. 73 , 3 (2012).
- [88] J. Diekmann and U. Thiele, Drops of volatile binary mixtures on brush-covered substrates, Eur. Phys. J.-Spec. Top. 10.1140/epjs/s11734-024-01169-4 (2024).
- [89] M. Kvick, D. M. Martinez, D. R. Hewitt, and N. J. Balmforth, Imbibition with swelling: capillary rise in thin deformable porous media, Phys. Rev. Fluids 2 , 074001 (2017).
- [90] K. P. Hapgood, J. D. Litster, S. R. Biggs, and T. Howes, Drop penetration into porous powder beds, J. Colloid Interface Sci. 253 , 353 (2002).
- [91] S. Wooh and H. J. Butt, A photocatalytically active lubricant-impregnated surface, Angew. Chem. Int. Ed. 56 , 4965 (2017).
Simon Hartmann Uwe Thiele
Gradient dynamics model for drops of volatile liquid on a porous substrate
We present a mesoscopic hydrodynamic model for a spreading drop of volatile
partially wetting liquid on a solid porous layer of small thickness. Thereby,
evaporation takes place under strong confinement, i.e., we consider a drop that
spreads on one of two parallel plates that form a narrow gap. Our gradient
dynamics model describes the coupled dynamics of the vertically averaged vapor
density profile, the drop height profile, and the vertically averaged
saturation profile in the porous layer. The underlying free energy incorporates
saturation-dependent wettability, capillarity of the drop and within the porous
layer, air and vapor entropy, and an entropic contribution relevant close to
complete filling of the porous layer. After developing the model we discuss the
resulting sorption isotherm and illustrate typical drop spreading and
imbibition behavior including the formation of a saturation halo within the
porous substrate.