## Mesoscopic hydrodynamic model for spreading, sliding and coarsening compound drops

Jan Diekmann 1, ∗ and Uwe Thiele 1,2,3, †

1

Institute of Theoretical Physics, University of M¨nster, Wilhelm-Klemm-Str. 9, 48149 M¨ unster, Germany u 2 Center for Nonlinear Science (CeNoS), University of M¨ unster, Corrensstr. 2, 48149 M¨nster, Germany u 3 Center for Multiscale Theory and Computation (CMTC),

University of M¨ unster, Corrensstr. 40, 48149 M¨ unster, Germany

We revisit the mesoscopic hydrodynamic description of the dynamics of sessile partially wetting compound drops, i.e., of drops that consist of two immiscible nonvolatile partially wetting liquids and are situated on a smooth rigid solid substrate. We briefly discuss and complete existing dynamic models employing a gradient dynamics approach. Thereby, the underlying energy features capillarity and wettability contributions for all relevant interfaces in full-curvature formulation. Establishing transparent consistency relations between macroscopic and mesoscopic parameters, we obtain mesoscopic Neumann and Young laws that are also fully consistent with the macroscopic ones. In particular, we discuss the minimal requirements for the wetting energy that ensure the full spectrum of macroscopic parameters for partially wetting cases is addressed by the mesoscopic model. Subsequently, we distinguish long-wave and full-curvature variants of the dynamical model based on properties of the energy, and employ the latter to illustrate the usage of the mesoscopic model. As examples, we chose the spreading of individual compound drops on one-dimensional horizontal substrates, sliding compound drops on one-dimensional inclined substrates, and the coarsening of drop ensembles on one- and two-dimensional horizontal substrates. In each case, the discussion emphasizes occurring qualitative changes in the drop configurations.

## Publication note

This version of the article has been accepted for publication, after peer review but is not the Version of Record and does not reflect all post-acceptance improvements and corrections. The Version of Record is available online at

Diekmann and Thiele, Physical Review Fluids 10, 024002, 2025, doi: 10.1103/PhysRevFluids.10.024002.

## I. INTRODUCTION

Bilayer liquid films and compound drops of two immiscible liquids on solid substrates have over time attracted considerable interest. Studied phenomena include, for instance, the dewetting of two-layer films of immiscible liquids [1-10], the spreading of drops on another liquid [11, 12], the development of long-wave thermocapillary or electro-static instabilities in such bilayer films [13-18], the occurrence of different steady configurations of sessile compound drops of two (or more) liquids [19-24], and the capillary leveling of two-layer films [25]. More involved systems are also investigated as, for instance, heated bilayer systems with soluble or insoluble surfactants and evaporation [26-28], sliding drops on liquid layers [29], and bilayer films of viscoelastic and elastic materials [30]. The decomposition of a film of a liquid mixture into a two-layer film that may at the same time dewet is also considered [31, 32]. The influence of slip is discussed in [33, 34]. Finally note that the description of drops on liquid-infused substrates is also closely related to two-layer systems [35-37].

Although descriptions of the dynamics of two-layer films with the full Navier-Stokes equations are possible [38], most modeling efforts focus on thin-film (or lubrication, or long-wave) models [39, 40] that allow one to consider a reduced dynamics in the form of coupled time evolution equations for either the thickness profiles of the two liquid layers or the height profiles of the two free interfaces [1, 4, 6, 7, 9, 12, 24, 41, 42].

Here, we focus on the isothermal case of thin films and droplets of maximal thicknesses well below the capillary length, i.e., we neglect hydrostatic influences and fully concentrate on a mesoscopic hydrodynamic description that faithfully incorporates substrate wettability and fluid-fluid interface energies. Such a description is commonly employed to study, for instance, the rupture of bilayer films [4, 41], the initial dewetting together with the later drop coarsening and related morphology changes for one-dimensional [5] and two-dimensional [43] substrates, the dewetting of a liquid layer and the subsequent dynamics of liquid droplets on a stable layer of another liquid [9], and the behavior of liquid lenses [12]. As the underlying long-wave approximation assumes a small ratio of typical length scales orthogonal and parallel to the solid substrate, the resulting models

∗ jan.diekmann@uni-muenster.de; ORCID ID: 0009-0005-4741-3368

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

are normally limited to small slopes of all involved interfaces. However, below we will also discuss a full-curvature variant of such models that combine long-wave dynamics with full-curvature energetics. For one-layer systems this approach is discussed, e.g., in Refs. [44-46] as well as in appendix A of [47]. This is closely related to the possibility of formulating all mentioned mesoscopic models [1, 4, 6, 7, 9, 12, 41, 42] as gradient dynamics on an underlying energy functional [4, 5, 9, 10, 24, 43, 48]. Note, however, that some of these formulations need to be carefully scrutinized to uncover their gradient dynamics form.

For static situations, e.g., sessile compound drops, a macroscopic description is also commonly used [19, 20, 22, 23]. Then, the system can be completely characterized by the involved fluid-solid and fluid-fluid interface energies. All involved fluid-fluid interfaces are of constant curvature governed by respective Laplace laws. All three-phase contact lines where three fluid-fluid interfaces meet correspond to circular arcs, and the angles between the interfaces are governed by Neumann's law (chap. 6.1 and 6.2 of [49], also compare [50]). For three fluids including the ambient gas or liquid there are three fluid-fluid and three fluid-solid interface energies. That is, macroscopically there are six energetic parameters of which two can be eliminated by general arguments (see below). For a complete nondimensional macroscopic description, one has to add one volume ratio. 1

Note that other approaches exist that, e.g., employ a phase-field model [21] or use piecewise thin-film descriptions for the individual interfaces coupled at sharp contact lines [10, 29]. However, then the equilibrium Neumann law is imposed at the triple line [29] what might be problematic in dynamic situations. A corresponding comparison of a mesoscopic thin-film model (that allows for adsorption layers, i.e., sharp three-phase contact lines are replaced by smooth transition regions) and a piecewise model for relatively slow relaxational dynamics is provided by [48]. The deviation from Neumann's law in the dynamic situation of forced wetting is discussed in [52] for a polymer brush-covered substrate that is moved at given velocity into a liquid bath.

The aim of the present work is to present a mesoscopic hydrodynamic model for compound drops of two immiscible nonvolatile partially wetting liquids on smooth rigid solid substrates that is based on a gradient dynamics approach. It shall be sufficiently general to (i) describe the initial two-layer dewetting dynamics and the subsequent coarsening process, (ii) capture all types of steady compound drops described in the literature as well as their dynamics (e.g., spreading, sliding, coarsening), and (iii) allow for consistency with all relevant macroscopic Laplace, Neumann and Young laws in the case of partial wetting. Furthermore, a long-wave variant and a full-curvature variant will be discussed based on different treatments of the underlying energy functional.

We argue that such a re-visitation is necessary as the mesoscopic models in the literature are normally not of a form that can be easily related to the appropriate macroscopic description, and do not allow for a consistent analysis of sessile compound drops in the full range of macroscopic parameters. For instance, although the two-layer models of Ref. [5, 43] are able to capture a number of compound drops (see e.g., Fig. 14 of [5]), the employed form of the wetting and interface energy does not lend itself easily to a clear discussion of consistency conditions between macroscopic and mesoscopic modeling or, indeed, to the discussion of a mesoscopic equivalent of Neumann's law. The two-layer models that focus on the initial stages of the dewetting and rupture process are not able to describe sessile compound drops because no stabilizing short-range interactions are included in the wetting energy [4, 41]. Other works assume the liquid of the lower layer to completely wet the substrate, i.e., again, no compound drops can emerge [9, 10, 12]. These restrictions also apply to other more complicated models that involve more components beside the two liquids, e.g., surfactants, further liquids, and/or evaporation [26].

The mesoscopic dynamic model is provided in section II in the form of a gradient dynamics on a mesoscopic energy functional that incorporates wetting and interface energies in the full-curvature variant. A corresponding long-wave variant is also given. The subsequent section III discusses important aspects related to the form of the energy functional, and shows that the resulting mesoscopic Laplace, Neumann and Young laws all converge to the correct macroscopic expressions if certain consistency conditions are fulfilled. Details of the related energy minimizations and resulting comparisons of mesoscopic and macroscopic Laplace, Young and Neumann laws are given in appendix A 3. The full-curvature variant of the presented model is then employed for a few examples in section IV: First, we investigate in section IV A the spreading of compound drops, then in section IV B we consider how the configuration of sliding compound drops may change, and in section IV C we illustrate the dewetting and coarsening behavior of a two-layer film. Finally, section V concludes and gives an outlook.

## II. GRADIENT DYNAMICS MODEL

We employ the governing equations for the time evolution of the bottom and top layer thickness profiles h 1 ( r , t ) and h 2 ( r , t ) , respectively, or alternatively, for the interface height profiles h 12 ( r , t ) and h 23 ( r , t ) , see sketches in Fig. 1. They are written in standard gradient dynamics formulation for two scalar fields with conserved dynamics [4, 5, 9, 10, 24, 42, 43, 46, 48], e.g., for

1 Note that here we do not take into account that the macroscopic interface energies might be related among themselves by microscopic considerations, e.g., in the case of purely apolar fluids by sum rules as, e.g., discussed in chapter 10 of [51].

FIG. 1. Sketches illustrate (a,c) macroscopic compound drops and (b,d) their mesoscopic equivalents where thin adsorption layers are present. Shown are (a,b) a drop-on-drop and (c,d) a drop-covers-drop configuration where adsorption layers occur at (b) the liquid 1-gas and the solid-gas interface, and (d) at the solid-liquid 2 and solid-gas interface. Macroscopic three-phase contact-points in (a) and (c) translate into mesoscopic three-phase contact regions in (b) and (d), respectively (indicated by gray shading). Magnifications of the relevant regions are given and illustrate corresponding Neumann and Young constructions in their macroscopic [(a) and (c)] and mesoscopic [(b) and (d)] form. Here, mesoscopic contact angles are determined at inflection points of the respective profiles. Note that the mesoscopic profiles in (b) and (d) represent actual numerically obtained steady states and are connected to the macroscopic profiles in (a) and (c) by the consistency relations discussed in the main text. Furthermore, panel (b) indicates the two different choices for the considered profiles: either in terms of layer thickness profiles h 1 and h 2 (red arrows) or in terms of interface height profiles h 12 and h 23 (gray arrows).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-bffd0f01-de19-4acf-a138-c30a59c6dbf4.jpg)

the layer thicknesses

$$\partial _ { t } h _ { 1 } & = - \nabla \cdot \mathbf j _ { 1 } = \nabla \cdot \left ( Q _ { 1 1 } \nabla \frac { \delta \mathcal { F } } { \delta h _ { 1 } } + Q _ { 1 2 } \nabla \frac { \delta \mathcal { F } } { \delta h _ { 2 } } \right ) \\ \partial _ { t } h _ { 2 } & = - \nabla \cdot \mathbf j _ { 2 } = \nabla \cdot \left ( Q _ { 2 1 } \nabla \frac { \delta \mathcal { F } } { \delta h _ { 1 } } + Q _ { 2 2 } \nabla \frac { \delta \mathcal { F } } { \delta h _ { 2 } } \right ) \,, \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \end{cases}$$

where F is the underlying energy functional and Q is a symmetric, positive definite mobility matrix ensuring positive entropy production and Onsager's reciprocity relations. The positions on the two-dimensional planar solid rigid substrate are indicated by r = ( x, y ) T , time is t , ∂ t stands for the partial time derivative, and ∇ = ( ∂ , ∂ x y ) T . Furthermore, j 1 and j 2 are the lateral fluxes within the respective liquid layer. For the presently used variables ( h 1 and h 2 ) one has [9, 33]

$$\underline { \mathbf q } = \begin{pmatrix} Q _ { 1 1 } & Q _ { 1 2 } \\ Q _ { 2 1 } & Q _ { 2 2 } \end{pmatrix} = \frac { 1 } { 6 \eta _ { 1 } } \begin{pmatrix} 2 h _ { 1 } ^ { 3 } & 3 h _ { 1 } ^ { 2 } h _ { 2 } \\ 3 h _ { 1 } ^ { 2 } h _ { 2 } & 2 \frac { \eta _ { 1 } } { \eta _ { 2 } } h _ { 2 } ^ { 3 } + 6 h _ { 1 } h _ { 2 } ^ { 2 } \end{pmatrix},$$

where η 1 and η 2 are the dynamic viscosities of bottom and top liquid, respectively. These mobilities are determined from a long-wave (lubrication) approximation of the governing hydrodynamic equations and appropriate boundary conditions (no-slip and no-penetration at solid substrate, velocity and shear stress continuity at liquid-liquid interface, zero shear stress at liquid-gas interface) [4, 5, 42] following the general procedure as outlined in the reviews [39, 40, 53]. For a case with slip at the solid substrate see Ref. [5]. Slip at liquid-solid and liquid-liquid interface is incorporated in [33].

The employed energy functional is

$$\mathcal { F } = \int _ { A } \left [ \gamma _ { s 1 } + \gamma _ { 1 2 } \, \xi _ { 1 2 } + \gamma _ { 2 3 } \, \xi _ { 2 3 } + g ( h _ { 1 }, h _ { 2 }, \xi _ { 2 3 } ) \right ] \, \mathrm d ^ { 2 } r \,.$$

In the full-curvature variant we use the exact metric factors ξ 12 = √ 1 + |∇ h 1 | 2 and ξ 23 = √ 1 + |∇ ( h 1 + h 2 ) | 2 , while they are approximated in the long-wave variant. 2 The first three terms correspond to the interface energies of the solid-liquid, liquidliquid and liquid-gas interface, respectively (with energy densities γ s1 , γ 12 and γ 23 ), while the fourth term collects all parts of the wetting energy. The particular form of g h , h ( 1 2 , ξ 23 ) including its dependence on a metric factor is discussed below in section III. As in the following we consider an isothermal situation, the energy F may be seen as a thermodynamic free energy.

The variations of F [ h , h 1 2 ] in Eqs. (1) correspond to the (vertically uniform) pressures within the respective layers

$$\frac { \delta \mathcal { F } } { \delta h _ { 1 } } = p _ { 1 } ( \mathbf r, t ) \quad \text{and} \quad \frac { \delta \mathcal { F } } { \delta h _ { 2 } } = p _ { 2 } ( \mathbf r, t ) \,.$$

Their gradients act as thermodynamic forces that drive the dynamics (1). In Eqs. (1) with (2) and (4) we denote the full-curvature variant of the mesoscopic hydrodynamic model for a two-layer system written in gradient dynamics form. The long-wave variant uses the approximation of the ξ 's and is discussed in section III D.

An advantage of the presented gradient dynamics form is the resulting straightforward way of changing the describing variables [47, 54]: Instead of writing the model in layer thicknesses h 1 and h 2 [9, 10, 48], alternatively, one may use the height profiles of the two free interfaces h 12 = h 1 and h 23 = h 1 + h 2 [4, 5, 33, 42]. The two representations are related by a linear transformation

$$\tilde { h } _ { \alpha } = \sum _ { \beta } R _ { \alpha \beta } h _ { \beta }$$

where here ( ˜ h , h 1 ˜ 2 ) stands for ( h 12 , h 23 ) and the 2 × 2 transformation matrix is R = ((1 0) (1 1)) , , , . As the variations transform as

$$\frac { \delta } { \delta h _ { \alpha } } = \sum _ { \beta } R _ { \beta \alpha } \frac { \delta } { \tilde { \delta h _ { \beta } } } \,,$$

one obtains from Eqs. (1) the transformed set of governing equations

$$\partial _ { t } \tilde { h } _ { \alpha } = - \nabla \cdot \tilde { \mathbf j } _ { \alpha } = \nabla \cdot \sum _ { \beta = 1 } ^ { 2 } \left [ \widetilde { Q } _ { \alpha \beta } ( \tilde { h } _ { 1 }, \tilde { h } _ { 2 } ) \nabla \frac { \delta \mathcal { F } } { \delta \tilde { h } _ { \beta } } \right ] \quad \text{for} \quad \alpha = 1, 2$$

where the transformed (again symmetric positive definite) mobility matrix is ˜ Q = RQR T . The lateral fluxes ˜ j α transform as the profiles, i.e., as in Eq. (5). Specifically, one has [4, 5, 33, 42]

$$\underline { \widetilde { Q } } = \frac { 1 } { \eta _ { 1 } } \begin{pmatrix} \frac { h _ { 1 2 } ^ { 3 } } { 3 } & \frac { h _ { 1 2 } ^ { 2 } } { 2 } \begin{pmatrix} h _ { 2 3 } - \frac { h _ { 1 2 } } { 3 } \end{pmatrix} \\ \frac { h _ { 1 2 } ^ { 2 } } { 2 } \begin{pmatrix} h _ { 2 3 } - \frac { h _ { 1 2 } } { 3 } \end{pmatrix} & \frac { ( h _ { 2 3 } - h _ { 1 2 } ) ^ { 3 } } { 3 } \begin{pmatrix} \eta _ { 1 1 } - 1 \end{pmatrix} + \frac { h _ { 2 3 } ^ { 3 } } { 3 } \end{pmatrix}.$$

Note that the transformation is also applied within the energy functional. Then, the two alternative formulations are fully equivalent. After reformulating by using Eq. (6), i.e., δ/δh ˜ α = ∑ β R -1 βα δ/δh β , the variations are

$$\frac { \delta \mathcal { F } } { \delta h _ { 1 2 } } & = \frac { \delta \mathcal { F } } { \delta h _ { 1 } } - \frac { \delta \mathcal { F } } { \delta h _ { 2 } } = p _ { 1 } - p _ { 2 } \\ \frac { \delta \mathcal { F } } { \delta h _ { 2 3 } } & = \frac { \delta \mathcal { F } } { \delta h _ { 2 } } = p _ { 2 } - p _ { 3 } \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

and correspond to pressure jumps across interfaces. The pressure in the gas phase p 3 is normally neglected in the case of nonvolatile liquids.

Note that beside the two alternative gradient dynamics formulations, in terms of layer thicknesses h 1 and h 2 (Eqs. (1)) as used, e.g., in Refs. [9, 10, 24, 33, 48], and in terms of interface heights h 12 and h 23 (Eq. (7)) as used, e.g., in Refs. [4, 5, 33, 42, 55], there exists also a mixed formulation, see Refs. [6, 7, 12, 15, 16, 18, 56-58]. There, the time derivatives are written in terms of h 12 and h 23 , but the fluxes are given in terms of the gradients of the pressures p 1 and p 2 within the layers, and not in terms of gradients of pressure jumps p 1 -p 2 and p 2 -p 3 across the interfaces at h 12 and h 23 , respectively. Although the

2 The sometimes conceived conflict between mathematically proper asymptotic usage of smallness parameters as in the long-wave variant (that then in practice are often not small) and an approximation scheme where the thermodynamic structure of the model is employed (and always preserved as the 'higher good') in approximation schemes like the full-curvature variant is discussed, e.g., in Refs. [44-46] and appendix A of [47]. We advocate the latter and independently approximate mobilities and energies (like it is common for, e.g., the Cahn-Hilliard equation where mobilities are often set constant and energies are expanded up to fourth order).

resulting models are still mathematically correct, the mixed formulation should be avoided as it obfuscates the underlying gradient dynamics structure. For instance, the resulting mobility matrix is not symmetric anymore. 3 This is inconvenient as it can not be employed to easily express the linear stability analysis as a generalized eigenvalue problem (cf. [5]) rendering the analysis quite cumbersome. Furthermore, in this formulation it may appear that Derjaguin pressures do not result from a single underlying wetting potential.

## III. STEADY STATES

First, we consider static situations, i.e., steady states described by Eqs. (1) with ∂ h t 1 = ∂ h t 2 = 0 . Integrating the resulting equations once, the integration constants are set to zero for the considered situation without fluxes across the lateral domain boundaries. Taking then into account that the mobility matrix Q is positive definite for h , h 1 2 &gt; 0 , i.e., det Q &gt; 0 , one can integrate a second time. One obtains the condition for steady states, namely, that the variations of F given by Eq. (4) are uniform in space, i.e., the pressures are equal to constants P 1 and P 2 that appear in the second integration. Equations (4) then correspond to coupled nonlinear ordinary differential equations for h 1 and h 2 that may be solved, e.g., by numerical continuation as done in Ref. [5], see their Figs. 8, 9, 12-14.

Resulting steady compound drops may show three-phase contact regions where three media meet, see Fig. 1(b),(d). On the one hand, this may be the liquid-liquid-gas contact of three fluids where macroscopically the Neumann law holds. On the other hand, this may be a solid-liquid-gas or a solid-liquid-liquid contact where two fluids and the solid substrate meet. There, macroscopically a respective Young law holds. The derivation of these macroscopic laws via minimization of a macroscopic equivalent of the energy (3) is briefly reviewed in Appendix A. Here, we next derive the mesoscopic version of these laws and subsequently provide consistency relations connecting the mesoscopic and macroscopic descriptions.

## A. Mesoscopic Laplace, Neumann and Young laws

In this section, we restrict the consideration to one-dimensional substrates and use the formulation in interface height profiles h 12 and h 23 . We start with the mesoscopic energy functional

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta ( 1 0 ) \\ \mathcal { F } _ { \text{meso} } = \int _ { - \infty } ^ { \infty } & ( \gamma _ { \text{s1} } + \gamma _ { \text{1} } 2 \xi _ { \text{1} } 2 + \gamma _ { \text{2} } 3 \xi _ { \text{2} } 3 + f _ { \text{1} } ( h _ { \text{1} } _ { \text{2} } ) + \xi _ { \text{2} } 3 f _ { \text{2} } ( h _ { \text{2} } - h _ { \text{1} } _ { \text{2} } ) + f _ { \text{3} } ( h _ { \text{2} } _ { \text{3} } ) \\ & \quad - P _ { \text{1} } ( h _ { \text{1} } _ { \text{2} } - \bar { h } _ { \text{1} } ) - P _ { \text{2} } ( h _ { \text{2} } _ { \text{3} } - h _ { \text{1} } _ { \text{2} } - \bar { h } _ { \text{2} } ) \Big ) \, \text{d} x \,, \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ delta ( 1 0 )$$

where we have directly included the constraints for volume conservation

$$\int _ { - \infty } ^ { \infty } ( h _ { 1 2 }$$

$$\begin{array} {$$

with respective Lagrange multipliers P 1 and P 2 , and respective mean layer thicknesses ¯ h 1 and ¯ h 2 . The energy F meso may still be seen as a thermodynamic Helmholtz free energy as the constraints for volume conservation do not contribute when introducing (10) into the mass-conserving dynamics (1). However, when considering steady states the constraints matter as they control the volume and one may see the functional as a Gibbs free energy.

Note that in Eq. (10) we have specified the general wetting energy using

$$g ( h _ { 1 }, h _ { 2 }, \xi _ { 2 3 } )$$

i.e., the sum of the interactions of liquid 2 with the substrate across liquid 1 ( f 1 ), of the ambient gas with liquid 1 across liquid 2 ( f 2 ), and of the ambient gas with the substrate across liquids 1 and 2 ( f 3 ). To capture partially wetting liquids we assume that each of the three f i features a single minimum at some microscopically small thickness and approaches zero for h →∞ . Below, for specific calculations we further restrict this to the combination of two specific antagonistic power laws. Situations where one

3 The transformation h = R -1 ˜ h is applied only on the left hand side of Eq. (1), therefore, ∂ t ˜ h = [ ̂ Q δ F ] with the asymmetric mobility matrix

$$\left \{ \tilde { h } _ { 1 } ^ { 2 } ( 3 \t$$

∇· ∇ δ h ̂ Q = RQ = 1 6 η 1 ( 2 ˜ h 3 1 3 ˜ h 2 1 ( ˜ h 2 -˜ h 1 ) ˜ h 2 1 (3 ˜ h 2 ˜ h 1 ) ( ˜ h 2 ˜ h 1 )[3 ˜ h 2 1 +6 ˜ h 1 ( ˜ h 2 ˜ h 1 ) + 2 η 1 η ( ˜ h 2 ˜ h 1 ) 2 ] ) where again ( ˜ h , h 1 ˜ 2 ) stands for ( h 12 , h 23 ) .

of the two liquids is completely wetting are not considered, but see, e.g., Refs. [9, 10] for drops of partially wetting liquid on a liquid substrate, and Ref. [59] for partially wetting sessile drops underneath a second (wetting) liquid.

It is instructive to consider an analogy between steady state compound drops and Newtonian point mechanics: Interface heights of steady states that vary with position along the substrate are seen as coordinates of pseudo point particles that vary with time. Then the integrand of the energy functional F meso can be seen as a Lagrangian L . This then allows one to define and determine the corresponding Hamiltonian H within the mechanical analogon. It represents a quantity (not a usual thermodynamic potential) that is uniform in space when discussing the thermodynamic steady states. In the mechanical analogon the constraint terms in Eq. (10) represent external force fields.

The corresponding Hamiltonian is given by

$$\mathcal { H } = \sum _ { j } \frac { \partial \mathcal {$$

where ∂ L ∂ ∂ ( x ˜ h j ) and ˜ h j are generalized momenta and generalized coordinates, respectively (with ˜ h 1 = h 12 and ˜ h 2 = h 23 as before). The Hamiltonian reads explicitly

$$\mathcal { H } = \gamma _ { 1 2 } \frac { ( \partial _ { x } h _ { 1 2 } ) ^ { 2 } } { \xi _ { 1 2 } } + ( \gamma _ { 2 3 } + f _ { 2 } ) \frac { ( \partial _ { x } h _ { 2 3 } ) ^ { 2 } } { \xi _ { 2 3 } } - \mathcal { L }$$

$$= - \frac { \gamma _ { 1 2 } } { \xi _ { 1 2 } } - \frac { \overset { \overset { \overset { \overset { \overset { \overset { \overset { \overs { \overs { \overs { \overs { \overs { \overs { \ \ } } } } } } } } { \gamma _ { 2 3 } } + f _ { 2 } } } } { \xi _ { 2 3 } } - \gamma _ { 3 1 } - f _ { 1 } - f _ { 3 } + P _ { 1 } ( h _ { 1 2 } - \bar { h } _ { 1 } ) + P _ { 2 } ( h _ { 2 3 } - h _ { 1 2 } - \bar { h } _ { 2 } ) \,.$$

Note that from hereon we use the energy E = -H . Minimization of F meso w.r.t. h 12, and h 23 gives the Euler-Lagrange equations

$$- \gamma _ { 1 2 } \kappa _ { 1 2 } + f _ { 1 } ^ { \prime } - \xi _ { 2 3 } f _ { 2 } ^ { \prime } = & P _ { 1 } - P _ { 2 } \,,$$

and

$$- \left ( \gamma _ { 2 3 } + f _ { 2 } \right ) \kappa _ { 2 3 } + \left ( \xi _ { 2 3 } - \frac { ( \partial _ { x } h _ { 2 3 } - \partial _ { x } h _ { 1 2 } ) \partial _ { x } h _ { 2 3 } } { \xi _ { 2 3 } } \right ) f _ { 2 } ^ { \prime } + f _ { 3 } ^ { \prime } = & P _ { 2 },$$

respectively, with curvatures κ 12 = ( ∂ xx h 12 ) /ξ 3 12 and κ 23 = ( ∂ xx h 23 ) /ξ 3 23 . In the mechanical analogon they correspond to coupled Newton's equations of motion for particles with 'velocity- and interaction-dependent masses' γ 12 /ξ 3 12 and ( γ 23 + f 2 ) /ξ 3 23 .

To obtain mesoscopic Laplace, Neumann and Young laws we consider three regions (i), (ii) and (iii) at large (macroscopic) distance from the two distinguished contact regions (macroscopically located at L 1 and L 2 , cf. Fig. 1(a) and 1(b)) and consider the corresponding limiting cases of Eqs. (16), (17) and (18). First, we consider region (i) where x ≫ L 1 , i.e., x →∞ . There, the substrate is only covered by ultrathin adsorption layers of uniform thicknesses h a1 and h a2 , implying that Eqs. (17) and (18) become after rearrangement

$$f _ { 1 } ^ { \prime } ( h _ { \mathbf a 1 } ) + f _ { 3 } ^ { \prime } ( h _ { \mathbf a 1 } + h _ { \mathbf a 2 } ) = P _ { 1 },$$

$$f _ { 2 } ^ { \prime } ( h _ { \mathrm a 2 } ) + f _ { 3 } ^ { \prime } ( h _ { \mathrm a 1 } + h _ { \mathrm a 2 } ) = P _ { 2 } \,.$$

Similarly, Eq. (16) becomes

$$E = \gamma _ { s 1 } + \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 1 } ( h _ { \mathbf a 1 } ) + f _ { 2 } ( h _ { \mathbf a 2 } ) + f _ { 3 } ( h _ { \mathbf a 1 } + h _ { \mathbf a 2 } ) - P _ { 1 } ( h _ { \mathbf a 1 } - \bar { h } _ { 1 } ) - P _ { 2 } ( h _ { \mathbf a 2 } - \bar { h } _ { 2 } ) \,.$$

Second, we consider the intermediate region (ii) where L 2 &lt; x &lt; L 1 with | x -L i | ≫ h a where h a represents the order of magnitude of the adsorption layers. There, the thickness of liquid 1 is large and the wetting energies f 1 and f 3 and their derivatives all approach zero. Furthermore, we assume the adsorption layer of liquid 2 on liquid 1 to be of uniform thickness ˆ h a2 , i.e., slopes and curvatures of the two profiles are identical ( ∂ h x 12 = ∂ h x 23, κ 12 = κ 23 ). Therefore, ξ 12 = ξ 23 = 1 / cos(Φ 13 ) , where Φ 13 is the angle of the local tangent of the curved interface with the horizontal. From (17) and (18) we then obtain

$$\gamma _ { 1 2 } \kappa _ { 1 2 } + \xi _ { 2 3 } f _ { 2 } ^ { \prime } ( \hat { h } _ { \mathbf a 2 } ) = P _ { 2 } - P _ { 1 } \\ / \quad \cdot \hat { x } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{cases}$$

$$\left ( \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { \partial 2 } ) \right ) \kappa _ { 2 3 } - \xi _ { 2 3 } f _ { 2 } ^ { \prime } ( \hat { h } _ { \partial 2 } ) = - P _ { 2 } \,.$$

Their summation gives the mesoscopic Laplace law valid in region (ii):

$$\left ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) \right ) \kappa _ { 2 3 } = - P _ { 1 } \,.$$

Similarly, Eq. (16) becomes

$$E = & \gamma _ { s 1 } + \cos ( \Phi _ { 1 3 } ) ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) ) - P _ { 1 } ( h _ { 1 2 } - \bar { h } _ { 1 } ) - P _ { 2 } ( \hat { h } _ { a 2 } - \bar { h } _ { 2 } ) \,. \\.$$

Note that as discussed below, in general, ˆ h a2 = h a2 .

Third, we consider the central region (iii) where x &lt; L 2 with | x -L 2 | ≫ h a . There, the thicknesses of both layers are large and all wetting energies can be neglected. Furthermore, we write the metric factors as ξ 12 = 1 / cos(Φ 12 ) and ξ 23 = 1 / cos(Φ 23 ) where Φ 12 and Φ 23 are the respective angles of the local tangents of the curved interfaces h 12 and h 23 with the horizontal. From (17) and (18) we then obtain the Laplace laws

̸

$$\gamma _ { 1 2 } \kappa _ { 1 2 } & = P _ { 2 } - P _ { 1 } \,, & & ( 2 6 ) \\ \varepsilon _ { \infty } \kappa _ { \infty } & = - P _ { \infty } & & ( 2 7 )$$

$$\gamma _ { 2 3 } \kappa _ { 2 3 } = - P _ { 2 } \,.$$

Finally, Eq. (16) becomes

$$E = \gamma _ { s 1 } + \gamma _ { 1 2 } \cos ( \Phi _ { 1 2 } ) + \gamma _ { 2 3 } \cos ( \Phi _ { 2 3 } ) - P _ { 1 } ( h _ { 1 2 } - \bar { h } _ { 1 } ) - P _ { 2 } ( h _ { 2 3 } - h _ { 1 2 } - \bar { h } _ { 2 } ) \,. \quad \, \quad \, \quad \, \quad ( 2 8 ) \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

Having discussed the individual asymptotic regions (i), (ii) and (iii), we now use that E, P 1 and P 2 are constant across the system, thereby connecting the previously obtained relations. In principle, this is possible for arbitrary macroscopic liquid volumes, i.e., for finite pressures P 1 and P 2 . However, here we consider drop volumes that are sufficiently large to assume P , P 1 2 → 0 . In other words, we consider ranges that mesoscopically are far outside the contact regions but macroscopically are close to contact. In consequence, Eqs. (25) and (28) relate the angles that the interfaces at the macroscopic meeting point of regions (ii) and (iii) form with the horizontal by

$$\gamma _ { 1 2 } \cos ( \Psi _ { 1 2 } ) + \gamma _ { 2 3 } \cos ( \Psi _ { 2 3 } ) = \cos ( \Psi _ { 1 3 } ) ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) ).$$

This is the horizontal component of the mesoscopic Neumann law.

The thicknesses h a1 and h a2 of the adsorption layers on the solid substrate far away from the sessile compound drop [in region (i)] are obtained from the two relations (19) and (20) with P , P 1 2 → 0 , i.e., from f ′ 1 ( h a1 ) + f ′ 3 ( h a1 + h a2 ) = 0 and f ′ 2 ( h a2 ) + f ′ 3 ( h a1 + h a2 ) = 0 . Note that in general, h a2 differs from the thickness ˆ h a2 of the adsorption layer of liquid 2 on liquid 1 in region (ii). The latter is obtained from Eq. (22) or (23), i.e., in the considered limit f ′ 2 ( ˆ h a2 ) = 0 .

To obtain the vertical component of the mesoscopic Neumann law we employ the mechanical analogon of the motion of particles and consider the generalized momenta ∂ L ∂ ∂ ( x ˜ h j ) = π j in region (ii) at some x &gt; L 2 and in region (iii) at some x &lt; L 2 . In particular, we use the conservation of the total generalized momentum π tot = π 1 + π 2 across the liquid-liquid-gas contact region at x = L 2 . It actually corresponds to a completely inelastic collision of two pseudo particles: Arriving from region (iii) the particles approach each other in the contact region, interact inelastically due to their interaction potential (the wetting energy f 2 ( h 23 -h 12 ) ) and continue as a bound state in region (ii). In this interpretation, Eq. (29) represents energy conservation. Specifically, it describes how the generalized kinetic energy after the collision is related to the inner energy of the bound state. This is further discussed in section III D.

Under the above assumptions ( P , P 1 2 → 0 ) the total momentum in region (iii) is

$$\pi _ { \text{tot} } = \gamma _ { 1 2 } \frac { \partial _ { x } h _ { 1 2 } } { \xi _ { 1 2 } } + \gamma _ { 2 3 } \frac { \partial _ { x } h _ { 2 3 } } { \xi _ { 2 3 } } = \gamma _ { 1 2 } \sin ( \Psi _ { 1 2 } ) - \gamma _ { 2 3 } \sin ( \Psi _ { 2 3 } ).$$

In region (ii), one has ∂ h x 12 = ∂ h x 23 and obtains

$$\pi _ { \text{tot} } & = \gamma _ { 1 2 } \frac { \partial _ { x } h _ { 1 2 } } { \xi _ { 1 2 } } + \frac { \partial _ { x } h _ { 2 3 } } { \xi _ { 2 3 } } \left ( \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) \right ) \\ & = - \gamma _ { 1 2 } \sin ( \Psi _ { 1 3 } ) - \left ( \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) \right ) \sin ( \Psi _ { 1 3 } ) \,. \\ \dots \dots \dots \dots$$

Taken together, the two relations give

$$\gamma _ { 1 2 } \sin ( \Psi _ { 1 2 } ) + \left ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) \right ) \sin ( \Psi _ { 1 3 } ) = \gamma _ { 2 3 } \sin ( \Psi _ { 2 3 } ) \,,$$

corresponding to the vertical component of the mesoscopic Neumann law.

Next, we consider energy conservation between regions (ii) and (i), i.e., across the solid-liquid-gas contact region at x = L 1 (again for P , P 1 2 → 0 ). Then, Eqs. (21) and (25) give for the angle Θ 13 that the straight interface in region (ii) forms with the horizontal

$$\cos ( \Theta _ { 1 3 } ) ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { \alpha 2 } ) ) = \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 1 } ( h _ { \alpha 1 } ) + f _ { 2 } ( h _ { \alpha 2 } ) + f _ { 3 } ( h _ { \alpha 1 } + h _ { \alpha 2 } ) \,,$$

which is the mesoscopic Young law for the three-phase contact of substrate, liquid 1 and gas. Note that within the mechanical analogon, it corresponds to a completely inelastic collision of a particle with a resting particle of infinite mass. The total inner energy after collision is given by f 1 + f 2 + f 3 . Due to the infinite mass, the momentum of the impacting particle is simply absorbed, and no equation of momentum conservation in the mechanical analogon can be given.

In other configurations of the compound drop, the interface between liquid 1 and liquid 2, and the interface between liquid 2 and fluid 3 may both independently approach the substrate (cf. in Fig. 1 the bottom row the inner and outer contact line, respectively), then similar considerations result in the mesoscopic Young laws

$$\gamma _ { 1 2 } \cos \Theta _ { 1 2 } = \gamma _ { 1 2 } + f _ { 1 } ( \hat { h } _ { \text{a1} } ) \,,$$

and

$$\gamma _ { 2 3 } \cos \Theta _ { 2 3 } = \gamma _ { 2 3 } + f _ { 1 } ( h _ { \mathbf a 1 } ) + f _ { 2 } ( h _ { \mathbf a 2 } ) + f _ { 3 } ( h _ { \mathbf a 1 } + h _ { \mathbf a 2 } ) - f _ { 1 } ( \hat { h } _ { \mathbf a 1 } ) \,,$$

respectively. Here, ˆ h a1 corresponds to the adsorption layer thickness of liquid 1 under a thick layer of liquid 2. Similar to the case of the macroscopic description (Appendix A) one may combine the three mesoscopic Young laws from Eqs. (33), (34) and (35) into

$$\gamma _ { 1 2 } \cos ( \Theta _ { 1 2 } ) + \gamma _ { 2 3 } \cos ( \Theta _ { 2 3 } ) = \cos ( \Theta _ { 1 3 } ) ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } ) )$$

that is of exactly the same form as the horizontal component of Neumann's law but with the three Young angles.

## B. Meso-macro consistency condition

Before we embark on a detailed discussion of the wetting energy we first establish conditions for the consistency of the presented mesoscopic description with the macroscopic picture. For ease of comparison, Appendix A briefly reviews the macroscopic energy, its minimization, and the resulting Laplace, Neumann and Young laws.

TABLE I. Conditions that ensure consistency between mesoscopic and macroscopic descriptions. Within the mesoscopic description, the wetting energies are shown in their general form (third column) as well as in the specific form resulting from Eq. (37) (fourth column). Further details are given in section III C.

| interface energy   | macroscopic   | mesoscopic [general]                                                   | mesoscopic [specific with (37) and (38)-(41)]                                      |
|--------------------|---------------|------------------------------------------------------------------------|------------------------------------------------------------------------------------|
| solid-liquid 1     | γ s1          | γ s1                                                                   | γ s1                                                                               |
| solid-liquid 2     | γ s2          | γ s1 + γ 12 + f 1 ( ˆ h a1 )                                           | γ s1 + γ 12 - 3 A 1 10 h 2 a1                                                      |
| solid-gas          | γ s3          | γ s1 + γ 12 + γ 23 + f 1 ( h a1 ) + f 2 ( h a2 ) + f 3 ( h a1 + h a2 ) | γ s1 + γ 12 + γ 23 - 3 A 1 10 h 2 a1 - 3 A 2 10 h 2 a2 - 3 A 3 10( h a1 + h a2 ) 2 |
| liquid 1-liquid 2  | γ 12          | γ 12                                                                   | γ 12                                                                               |
| liquid 1-gas       | γ 13          | γ 12 + γ 23 + f 2 ( ˆ h a2 )                                           | γ 12 + γ 23 - 3 A 2 10 h 2 a2                                                      |
| liquid 2-gas       | γ 23          | γ 23                                                                   | γ 23                                                                               |

The macroscopic description is based on the six interface energies: liquid 1 to gas γ 13 , liquid 2 to gas γ 23 , liquid 1 to liquid 2 γ 12 , liquid 1 to solid substrate γ s1 , liquid 2 to substrate γ s2 , and gas to substrate γ s3 . In contrast, the mesoscopic description presented above is based on the same γ 23 , γ 12 , γ s1 , and, additionally, the wetting energy. In other words, three of the macroscopic parameters are replaced in the mesoscopic description by the wetting energy that is a function of the two layer thicknesses.

Comparing the mesoscopic Laplace, Neumann, and Young laws [(24), (26), (27); (29), (32); (33), (34), (35)] with their macroscopic equivalents [(A14), (A12), (A13); (A23), (A21); (A24), (A26), (A25), respectively], one finds that full consistency of macroscopic and mesoscopic description is ensured by the equivalences listed in Table I. Next, we discuss more in detail the conditions this imposes on the wetting energy.

## C. Wetting energies

When passing from the general gradient dynamics model in section II to the static considerations in section III we have in Eq. (13) specified the general form of the wetting energy to be the sum of three contributions. Thereby, f 1 ( h 12 ) encodes the

interaction of the bulk substrate with a bulk liquid 2 across a mesoscopic layer of liquid 1, f 2 ( h 23 -h 12 ) encodes the interaction of a bulk liquid 1 with the bulk gas phase across a mesoscopic layer of liquid 2, and f 3 ( h 23 ) represents a correction that ensures that the sum of all three terms encodes the interaction of the bulk substrate with the bulk gas phase across mesoscopic layers of liquid 1 and liquid 2. The derivation of the mesoscopic Neumann law has made it clear that f 2 has to be multiplied by the metric factor ξ 23 to correctly contribute to the macroscopic interface energy γ 13 .

̸

An obvious question is whether the term f 3 ( h 23 ) is indeed necessary. Pursuing this by setting f 3 = 0 , i.e., by assuming that wetting interactions can be completely encoded by adding one contribution only related to the lower layer and another contribution only related to the upper layer, would directly result in the identifications ˆ h a1 = h a1 and ˆ h a2 = h a2 . This in turn would imply the relation γ s2 -γ 12 = γ s3 -γ 13 between macroscopic quantities. As such a relation does, in general, not hold, the simplification f 3 = 0 should not be employed. In other words, the contribution f 3 ( h 1 + h 2 ) has to be incorporated for a faithful mapping of all macroscopic parameters for partially wetting liquids into the mesoscopic model. Note that one is still able to construct a model such that ˆ h a1 = h a1 and ˆ h a2 = h a2 as well as f 3 ( h a1 + h a2 ) = 0 . In terms of the general wetting energy g h , h ( 1 2 ) the condition is g h ( a1 , h a2 ) = g h ( ˆ a1 , ∞ ) + g ( ∞ , h ˆ a2 ) .

̸

Overall, the formulation of the mesoscopic model either with layer thicknesses h 1 and h 2 or interface height profiles h 12 and h 23 together with the above discussion of the three contributions to the wetting energy allows for a consistency check that should be applied to existing literature models: In the formulation with h 1 and h 2 one should have three terms f 1 ( h 1 ) + f 2 ( h 2 ) + f 3 ( h 1 + h 2 ) (here, for simplicity we dropped the metric factor). This implies that the variations of the energy w.r.t. h 1 , and h 2 should contain contributions f ′ 1 ( h 1 ) + f ′ 3 ( h 1 + h 2 ) , and f ′ 2 ( h 2 ) + f ′ 3 ( h 1 + h 2 ) , respectively. In contrast, in the formulation with h 12 and h 23 one should have three terms f 1 ( h 12 ) + f 2 ( h 23 -h 12 ) + f 3 ( h 23 ) . Then, the variations of the energy w.r.t. h 12 , and h 23 should contain contributions f ′ 1 ( h 12 ) -f ′ 2 ( h 23 -h 12 ) , and f ′ 2 ( h 23 -h 12 ) + f ′ 3 ( h 23 ) , respectively. Literature models that do not follow such a structure either unduly limit the macroscopic parameter space or break basic physical principles.

For the example calculations shown below, we employ a specific thickness-dependence in f 1 , f 2 and f 3 . We combine a long-range van der Waals interaction with a short-range power law, i.e.,

$$f _ { i } ( h ) = \frac { A _ { i } } { 2 h ^ { 2 } } \left [ \frac { 2 } { 5 } \left ( \frac { h _ { i } ^ { * } } { h } \right ) ^ { 3 } - 1 \right ]$$

where h stands for the relevant argument, the A i are Hamaker-type constants, and the h ⋆ i are parameters that represent typical length scales related to the equilibrium adsorption layer thicknesses h a1 , h a2 , ˆ h a1 and ˆ h a2 . For A &gt; i 0 the short-range interaction is repulsive (here, stabilizing) while the long-range one is attractive (here,, destabilizing) as appropriate for partially wetting liquids. The function f i ( h ) has a minimum f i ( h ⋆ i ) = -3 A / i 10( h ⋆ i ) 2 at h = h ⋆ i , diverges toward ∞ for h → 0 , and approaches zero from below at large h → ∞ . Note that this relatively simple choice is sufficient for our present purposes. However, to study more involved phenomena, e.g., involving first order wetting transitions, more complicated wetting energies should be employed, e.g., featuring a maximum at an intermediate thickness or several minima [60-63].

Assuming flat interfaces and large drops ( P , P 1 2 , κ 12 → 0 ), one can determine the thicknesses of all coexisting adsorption layers by solving Eqs. (19), (20), (22) and the equilibrium condition for an adsorption layer of liquid 1 underneath a bulk liquid 2, i.e., by solving the system

$$f _ { 1 } ^ { \prime } ( h _ { \mathbf a 1 } ) & = - f _ { 3 } ^ { \prime } ( h _ { \mathbf a 1 } + h _ { \mathbf a 2 } ) \\ f _ { 2 } ^ { \prime } ( h _ { \mathbf a } \wedge & = - f _ { 3 } ^ { \prime } ( h _ { \mathbf a } + h _ { \mathbf a } \wedge$$

$$0 = f _ { 2 } ^ { \prime } ( \hat { h } _ { \mathrm a 2 } )$$

$$\ddot { f } ^ { \prime } _ { 2 } ( h _ { \mathbf a 2 } ) & = - f ^ { \prime } _ { 3 } ( h _ { \mathbf a 1 } + h _ { \mathbf a 2 } ) \\ \cap & = e ^ { \prime } / \hat { \imath } \quad \real$$

$$0 = f _ { 1 } ^ { \prime } ( \hat { h } _ { \mathbf a 1 } )$$

for h a1 , h a2 , ˆ h a1 and ˆ h a2 . The specific form (37) with independent parameters h ⋆ i may be numerically solved. However, we can further simplify and still access the complete relevant macroscopic parameter space discussed in Table I. Namely, we assume that f 1 , f 2 , and f 3 have their minima at h a1 , h a2 , and h a1 + h a2 , respectively. This holds if in Eq. (37) one sets h ⋆ 1 = h a1 , h ⋆ 2 = h a2 and h ⋆ 3 = h a1 + h a2 , and further implies ˆ h a1 = h a1 and ˆ h a2 = h a2 .

A typical resulting wetting potential is presented in nondimensional form in Fig. 2. Note that expressions (38) to (41) only hold for P , P 1 2 → 0 . In other words, the adsorption layer thicknesses slightly change for finite pressures. Normally, an increase in pressure results in a slight increase in the corresponding thickness. Geometrically, finite pressures within the two liquids (e.g., Laplace pressures within compound drops of finite volumes) result in a tilt of the energy in Fig. 2 that slightly shifts the minima.

Introducing the resulting f i ( h ⋆ i ) for P , P 1 2 → 0 into the expressions of column 3 of Table I one obtains the consistency relations given in column 4 of Table I. Vice versa, the three Hamaker constants A i can be directly related to the macroscopic

FIG. 2. Shown is a typical wetting potential [(13) and (37)] in the scaled form g/ h ( 2 a1 A 1 ) = [ f 1 ( h 1 ) + f 2 ( h 2 ) + f 3 ( h , h 1 2 )] / h ( 2 a1 A 1 ) with equal Hamaker constants and for homogeneous films, i.e., ξ 23 = 1 . Panel (a) gives the potential as a contour plot in the ( h , h 1 2 ) -plane where positive values are suppressed for the sake of clarity. Specific values of chosen equipotential lines are given in the color bar. Panels (b) and (c) present cuts at several fixed h 1 and h 2 , respectively, as indicated in (a) by respective vertical and horizontal lines of corresponding line style.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-2ee2a359-363c-48b8-bb59-a54315c02d1e.jpg)

interface energies by

$$A _ { 1 } = - \frac { 1 0 h _ { a 1 } ^ { 2 } } { 3 } ( \gamma _ { s 2 } - \gamma _ { s 1 } - \gamma _ { 1 2 } )$$

$$A _ { 2 } = - \frac { \lambda } { \stackrel { 1 0 h _ { a 2 } ^ { 2 } } { 3 } } ( \gamma _ { 1 3 } - \gamma _ { 1 2 } - \gamma _ { 2 3 } ) \\ \stackrel {.. } {.. } \, \stackrel {.. } {.. } \, \stackrel {.. } {. }$$

$$A _ { 3 } = - \frac { \lambda } { 3 } ( \gamma _ { 1 2 } - \gamma _ { 1 3 } + \gamma _ { s 3 } - \gamma _ { s 2 } ) \,.$$

As in the present setup, the Hamaker constants have to be positive to ensure short-range repelling and long-range attracting interactions between interfaces, Eqs. (42)-(44) imply three inequalities that the macroscopic interface energies have to fulfill to allow for a mesoscopic description by a wetting energy that combines expressions of the form (37). For instance, Eq. (43) shows that A &gt; 2 0 directly corresponds to the macroscopic Neumann triangle inequality γ 12 + γ 23 &gt; γ 13 .

Furthermore, the sequence of the two liquids has to be chosen as energetically favorable, i.e., γ s1 + γ 12 + γ 23 &lt; γ s2 + γ 12 + γ 13 , therefore, a further inequality to hold is γ s1 -γ s2 &lt; γ 13 -γ 23 . Note, that (44) well illustrates that neglecting f 3 (i.e., A 3 = 0 ) is only valid for the particular values that fulfill γ s2 -γ 12 = γ s3 -γ 13 , i.e., represents a nongeneric case.

In summary, the restrictions of the considered parameter space, on the one hand, result from our focus on the case of partial wettability. On the other hand, they result from our mesoscopic modeling approach where, e.g., all contact angles need to be smaller than π/ 2 , and situations like encapsulation of one of the liquids inside a drop of the other liquid are not considered.

Up to here, all relations are derived within the full-curvature variant of the mesoscopic model. In the next section, we consider how the various obtained laws and relations simplify when considering the long-wave variant.

## D. Long-wave approximation

To apply a long-wave approximation [39, 40] one assumes that all relevant length scales parallel to the substrate are large compared to relevant vertical scales. This implies that all interface slopes, in particular, Young angles and Neumann angles w.r.t. the horizontal, have to be small. In the energy functional (10) the approximated metric factors are ξ lw 12 = 1 + 1 2 |∇ h 12 | 2 and ξ lw 23 = 1 + 1 2 |∇ h 23 | 2 resulting as before in the variations (17) and (18) but with approximated curvatures κ lw 12 = ∂ xx h 12 and κ lw 23 = ∂ xx h 23 .

Considering as above the three regions (i), (ii) and (iii) (cf. Fig. 1) and subsequently the conservation of energy and generalized

momentum in the mechanical analogon, we obtain the two components of the long-wave Neumann law

$$\gamma _ { 1 2 } \Psi _ { 1 2 } + \left ( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { a 2 } )$$

$$\gamma _ { 1 2 } \Psi _ { 1 2 } ^ { 2 } + \gamma _ { 2 3 } \Psi _ { 2 3 } ^ { 2 } - ( \gamma _ { 1 2 } + \gamma$$

and the three Young laws

$$( \gamma _ { 1 2 } + \gamma _ { 2 3 } + f _ { 2 } ( \hat { h } _ { \mathbf a _ { 2 } } ) ) \Theta _ { 1 3 } ^ { 2 }$$

$$\gamma _ { 1 2 } \Theta _ { 1 2 } ^ { 2 } = - 2 f _ { 1 } ( \hat { h } _ { \mathfrak { a } 1 } ). \quad \quad \quad$$

$$\gamma _ { 2 3 } \Theta _ { 2 3 } ^ { 2 } & = - 2 ( f _ { 1 } ( h _ { \mathbf a 1 } ) + f _ { 2 } ( h _ { \mathbf a$$

Using the consistency conditions in Table I, Eqs. (45) to (49) are fully consistent with the long-wave variant (A28) to (A32), respectively, of the macroscopic laws in Appendix A.

Note that the mechanical analogon discussed above in the context of the vertical component of Neumann's law, in the longwave case corresponds to standard Newtonian mechanics of two point particles: For P , P 1 2 → 0 , Eqs. (46) and (45) correspond to energy and momentum conservation during a collision of point particles at positions h 12 and h 23 with interaction potential f 2 ( h 23 -h 12 ) , see above. The interface energies γ 12 and γ 23 then represent the corresponding inertial masses. Incorporating finite pressures into the relations would amount to adding an external gravitational potential with the respective gravitational masses being proportional to P 1 and P 2 .

Our discussion of the wetting energy g has made it clear that the Derjaguin (or disjoining) pressures that enter the pressure expressions for the two layers are obtained as derivatives of the total wetting energy w.r.t. the relevant layer heights or thicknesses (depending on the particular formulation, see section II). 4 Also note that none of the literature models cited in sections II and III takes the metric factor in front of f 2 ( h 23 -h 12 ) into account, cf. Eq. (10). In consequence, to our knowledge, none of the existing mesoscopic two-layer models provides exact mesoscopic Neumann and Young laws either within the full-curvature variant or the long-wave variant that are consistent with the corresponding macroscopic laws.

## IV. APPLICATIONS

In the previous section, we have introduced a mesoscopic hydrodynamic model for a two-layer system of two nonvolatile immiscible partially wetting liquids. The employed gradient dynamics framework ensures thermodynamic consistency. Furthermore, consistency with the macroscopic description of two-liquid compound drops is guaranteed by three conditions formulated for the wetting energy that relate macroscopic parameters to mesoscopic ones. Additionally, the vertical order of layering implies a fourth condition in the form of an inequality that connects four of the six macroscopic interface energies. Next, we apply the established model to a number of typical situations and discuss selected results obtained by direct numerical simulations of the dynamic model (1). We remind the reader that one may use either of the two given mesoscopic formulations, i.e., in terms of layer thicknesses (Eq. (1)-(2)) or in terms of interface heights (Eq. (7) and (8)). Both exist as full-curvature and long-wave variants. Here, we employ the full-curvature variant and effectively nondimensionalize the system by setting γ 12 = η 1 = h a1 = 1 . Additionally, we set the interface energy between liquid 1 and the substrate to γ s1 = 0 as only differences in substrate energies enter (macroscopic) Young laws. All numerical results are obtained by using the finite-element library oomph-lib providing space- and time-adaptivity [64].

## A. Spreading compound drops

First, we numerically investigate how sessile compound drops spread. We employ a symmetric initial configuration, namely, a drop of liquid 1 that is entirely covered by a drop of liquid 2, i.e., a drop-covers-drop configuration. We consider a onedimensional substrate, i.e., the liquid-liquid and liquid-gas interfaces are circular arcs. The volume of liquid 1 is fixed at V 1 = 10 5 and two different values for the volume V 2 of liquid 2 are used. Note that these volumes are defined as the excess volume above the adsorption layers, e.g., V 1 = ∫ ( h 1 -h a1 )d x . The dynamic viscosity of liquid 2 is η /η 2 1 = 0 8 . and its adsorption layer height is h a2 /h a1 = 1 . The interface energies are set to ( γ 13 , γ 23 , γ s2 , γ s3 ) = (2 , 1 3 . , 0 9 . , 1 6) . such that all four inequalities discussed at the end of Section III C hold. The resulting macroscopic equilibrium contact angles are ( θ 12 , θ 13 , θ 23 ) ≈ (25 8 . ◦ , 36 9 . ◦ , 57 4 . ◦ ) . As the initial droplets have contact angles θ init 12 = θ init 23 = 72 ◦ larger than the corresponding equilibrium values both liquids

4 If this is not taken into account Newton's third law is broken. In our opinion this occurs in Ref. [8] (their Eqs. (3)-(5) are inconsistent, possibly due to a missing term in Eq. (3)).

FIG. 3. Shown are snapshots of the relaxational spreading dynamics of a compound drop for (a) V /V 2 1 = 3 and (b) V /V 2 1 = 1 . In both cases, an initial drop-covers-drop configuration is used (left-most panels), i.e., a drop of liquid 2 entirely covers a drop of liquid 1. As reflection symmetry holds during the entire evolution, only the right hand side of the profiles is shown. Streamlines indicate velocity fields. Additionally, their absolute values are shown by colored shading. Further details and parameter values are given in the main text, further analysis is provided in Fig. 4. Also, consider the corresponding movies (file numbers 00 and 01) in the Supplementary Material [65].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-0f656bb3-083a-4324-bdfa-bd7db0a51c68.jpg)

spread. Snapshots from the simulations of the spreading processes are shown in Fig. 3 (a) and (b) for the volume ratios V /V 2 1 = 3 and V /V 2 1 = 1 , respectively.

For the compound drop of large volume ratio, the configuration does not change qualitatively during the relaxational spreading process, see Fig. 3(a). Both liquids spread till the contact angles have reached their equilibrium values, and always remain in drop-covers-drop configuration. A deeper insight into the dynamics of the deformable interfaces is attained by inspecting the vertically resolved velocity fields within the two liquids that are reconstructed from the film thickness profiles as detailed in appendix B 1. These fields are incorporated in the respective two central panels of Figs. 3 (a) and (b) as selected streamlines and coloring that encodes the absolute value of the local velocity. The larger velocities in the vicinity of the contact lines result from strong changes in pressure gradients (cf. B 1). 5

In contrast, at the smaller volume ratio employed in Fig. 3 (b) the configuration changes during relaxation. As in the previous case, both liquids initially spread, however, at t ≈ 3 2 . · 10 4 the contact line where the liquid 1-liquid 2 interface meets the substrate catches up with the one of the liquid 2-gas interface. In consequence, the substrate-liquid 2 interface shrinks to zero and a transient four-phase contact is formed. Afterward, the liquid 1-liquid 2-gas contact separates from the substrate and moves inwards, i.e., a drop of liquid 2 retracts onto the drop of liquid 1 (that still spreads) until equilibrium Neumann angles are reached. In the final state, liquid 2 forms a drop positioned symmetrically on top of a drop of liquid 1, i.e., a drop-on-drop configuration. During the spreading process the liquid flow can either form a convection roll-like structure in one of the liquids or across both liquids as in the third panels of Figs. 3 (a) and 3 (b) or be predominantly parallel as in the second panel of Fig. 3 (b) when both drops spread at similar speeds. Note that the nearly vertical structure visible in some panels at the location of the inner contact line reflects that, e.g., a faster spreading inner drop can drag the upper liquid along and 'push it' into the more quiescent outer liquid in front of the advancing inner contact line. A magnification of such a region is shown in the inset of the third panel of Fig. 3 (b). We point out that simulations with the long-wave variant of the model at identical parameters and initial conditions as Fig. 3 (not shown) result in identical qualitative behavior as obtained with the full-curvature variant. In particular, the equilibrated compound drops show the identical configurations in both variants. Instead of circular arcs in the full-curvature variant the final height profiles in the long-wave variant are composed by parabolas. Only pressure gradients and, therefore, driving forces and velocities quantitatively differ for the two variants.

Further quantitative information about the dynamic processes showcased in Fig. 3 can be found in Fig. 4 (a)-(c) for V /V 2 1 = 3 , and Fig. 4 (d)-(f) for V /V 2 1 = 1 . Thereby, the first column displays the dependence of the free energy (3) on time. When both drops continuously spread keeping the initial drop-covers-drop configuration [Fig. 4 (a)] the energy decreases rather smoothly till it approaches its equilibrium values. In the process, interface energy is viscously dissipated within liquid 1 (rate D 1 ) and liquid 2 (rate D 2 ), see Fig. 4 (b). The dissipation rates are determined using Eq. (B5) that is derived in Appendix B 2. As the viscous dissipation in the vicinity of a contact line dominates (cf. Ref. [66]) and is larger for smaller contact angles, in the considered case of similar viscosities the always smaller contact angle ( θ 12 &lt; θ 23 ) at and the larger velocity of the inner contact

5 Due to the limited exactness of the plotting routines the coloring should be interpreted as providing the relative strength of the local flow.

FIG. 4. The top and bottom row provide further analysis of the dynamics shown in Figs. 3 (a) and 3 (b), respectively. Given are the time dependencies of [(a),(d)] the free energy [Eq. (3)], [(b),(e)] the dissipation rates in liquid 1 ( D 1 ) and liquid 2 ( D 2 ), (cf. section B), and [(c),(f)] the contact line positions. The insets in panels (e) and (f) magnify the time interval when the two contact line regions cross.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-1d63c7b5-3dd7-4b63-86d3-76f596383701.jpg)

FIG. 5. Shown are equilibrated compound drops together with corresponding space-time diagrams of the contact line positions where the liquid 1-liquid 2 or the liquid 1-gas interface meet the substrate (blue line) and where the liquid 2-gas interface meets the substrate or the three fluids meet (red line). The volume ratios V /V 2 1 are (a) 0.5, (b) 0.8, and (c) 3.3. The interface energies are ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 3 . , 0 31 0 95 1 1) . , . , . while all other parameters are as in Fig. 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-a2342bd5-394f-4cc2-9286-2b1a62088cfa.jpg)

line implies D &gt; D 1 2 for the entire process as observed in Fig. 4 (b). Finally, the third column of Fig. 4 characterizes the overall dynamics by showing the two contact line positions over time. It clearly indicates that the configuration does not change (no change in sequence and/or number of contacts), that the inner contact line is always faster than the outer one (steeper slope of the red line), and that the outer contact angle is the first to reach equilibrium. The slight overshoot of the outer contact line is a dynamic effect of the still spreading inner drop.

In contrast to Fig. 4 (a), in the case of Fig. 4 (d) the energy seems to converge at t ≈ 10 4 but then abruptly decreases in a

nearly step-wise manner. This is related to the approach of the two contact lines (still as drop-covers-drop configuration) into the short-lived transient of four-phase coexistence and the subsequent configuration change to drop-on-drop configuration, i.e., the outer contact line has transformed into an inner (dynamic) Neumann region. This dramatic change is also reflected in features of the time dependencies of dissipation rates depicted in Fig. 4 (e). Both show a local minimum at the transient four-phase coexistence where contact regions 'cross' and their speeds are minimal, see Fig. 4 (f). Afterwards, D 2 and D 1 first sharply rise before decreasing again till finally equilibrium is reached. Close inspection of Fig. 4 (f) shows that before contact regions cross at t ≈ 10 4 the outer (liquid 2-gas) contact line slightly recedes. Then, directly after the crossing the same happens for the new outer (liquid 1-gas) contact line. This indicates a direct attractive interaction between contact regions. As the four-phase coexistence only exists for a very brief transient time span it most likely does not correspond to a metastable state (a local minimum of the energy). However, the seemingly convergent energy shortly before the step-change indicates that it might correspond to an unstable steady state representing a saddle point of the energy. Such a saddle is approached by the trajectory of the system along a stable direction in phase space before it is left behind when the trajectory is captured by an unstable direction, and ultimately converges to a stable state.

We further scrutinize the occurrence of four-phase contact by considering a different set of interface energies where, starting with initial conditions as in Fig. 3, for a range of volume ratios the compound drop may actually relax toward a four-phasecontact configuration. In particular, for ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 3 0 31 0 95 1 1) . , . , . , . (resulting in equilibrium contact angles ( θ 12 , θ 13 , θ 23 ) ≈ (18 2 . ◦ , 32 2 . ◦ , 61 1 . ◦ ) ), we find the equilibrated compound drops shown in the top row of Fig. 5, also see the corresponding space-time plots of the contact lines positions in the bottom row. For V /V 2 1 = 0 5 . [Fig. 5(a)] and V /V 2 1 = 3 3 . [Fig. 5(c)] the resulting equilibrium configurations are similar to those of Fig. 3(b) and (a), respectively, i.e., drop-on-drop and drop-covers-drop. However, in contrast to Fig. 3(b), in the case of Fig. 5(a) there is a finite time-span with four-phase contact where both contact regions move together. Subsequently, they separate again, the newly formed Neumann contact line region moves inwards while the liquid 1-gas-substrate contact line continues to move outwards before both reach their equilibrium positions in the drop-on-drop configuration. One may say that the motion, interaction and crossing of the two contact line regions parallels certain aspects of soliton-soliton interactions [67].

However, in a small range of intermediate volume ratios, e.g., at V /V 2 1 = 0 8 . shown in Fig. 5 (b), the contact line regions do not separate again and the four-phase-contact configuration corresponds to the equilibrium compound drop. The corresponding space-time plot reveals that this result is not just a coincidence for one specific volume ratio because the four-phase-contact is formed at a smaller x -value than the one of the final state, i.e., the combined contact region continues to move outwards towards the equilibrium state. The possibility of such an equilibrium configuration with four-phase contact that exists for a finite range of volume ratios is predicted in Ref. [19] based on macroscopic considerations as a state between drop-covers-drop and dropon-drop configuration. This prediction is here confirmed employing a dynamic mesoscopic hydrodynamic model. However, it remains a task for the future to identify the corresponding parameter ranges where the various states are stable, metastable or unstable steady states of the model.

## B. Sliding compound drop

In the previous section, we have considered the spreading of compound drops on horizontal substrates and discussed their relaxational dynamics toward various equilibrium configurations. Next, we investigate sliding compound drops on a homogeneous inclined substrate. As initial conditions we use the equilibrium drops obtained before. Thereby, we incorporate gravity while still considering small droplets. To focus on the sliding motion we neglect the hydrostatic pressure and only account for the symmetry-breaking downhill force parallel to the substrate. We introduce it into the free energy F meso [Eq. (10)] in analogy to Refs. [68, 69]. In the interface height formulation, the energy reads

$$\cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad\cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd� \quad x h _ { 2 3 } d x,$$

with the inclination parameter β = ρg sin α , where ρ is the common mass density of the liquids, g is the gravitational acceleration, and α is the inclination angle. Note that a generalization for different mass densities of the two liquids is straightforward but not pursued here. Inserting ˜ F meso into Eq. (7) the right-hand sides of the two dynamic equation acquire additional terms -β∂ x ( ˜ Q 12 ) and -β∂ x ( ˜ Q 22 ) , respectively.

Figs. 6 (a) and (c) show snapshots of resulting time simulations for the volume ratios V /V 2 1 = 3 and V /V 2 1 = 1 , respectively. The inclination parameter is β = 10 -6 while all other parameters are as in Fig. 3. Initially, at t = 0 , the center-of-mass (CoM) of the equilibrated compound drops are positioned at x = -1250 within a computational domain of total size L = 5000 . The respective next three panels display profiles at later times. Thereby the final two show converged stationary sliding states at times when the CoM of liquid 1 is at x = 0 and x = 1000 , respectively. Fig. 6 (b) and (d) give the CoM-velocities as a function of time in a log-log representation.

Inspecting the time evolution in Fig. 6 (a) one sees that the initial symmetric drop-covers-drop configuration at V /V 2 1 = 3 first becomes asymmetric under the influence of the lateral driving force (also see movie with file number 02 in the Supplementary

FIG. 6. Time evolution of sliding compound drops: Equilibrium compound drops with (a)-(b) V /V 2 1 = 3 and (c)-(d) V /V 2 1 = 1 obtained in section IV A are placed on an inclined plane with β = 1 · 10 -6 . The respective upper four panels (in (a) and (c)) show snapshots at selected times while (b) and (d) show a log-log representation of the time dependence of the center-of-mass (CoM) velocity individually for the two liquids. For further details see main text. Other parameters and shading are as in Fig. 3. Also consider the movies (file numbers 02 and 03) in the Supplementary Material [65].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-182c8201-2387-473a-8a10-63dccdf8c1c5.jpg)

Material [65]). In other words the smaller inner drop slides more slowly than the outer drop. As a result, their trailing contact regions approach till the liquid 2-gas-substrate contact has caught up with the liquid 1-liquid 2-substrate contact at t ≈ 1 2 . · 10 4 . Then, the configuration changes into a strongly asymmetric drop-beside-drop configuration, i.e., the trajectories of the contact regions cross. In the process, the liquid 1-liquid 2-substrate contact becomes the trailing liquid 1-gas-substrate contact while the liquid 2-gas-substrate contact becomes a liquid 2-liquid-1-gas contact, i.e., a dynamic Neumann region. In parallel, the trailing liquid 2-substrate interface vanishes and a liquid 1-gas interface emerges. This transition results in a small peak of the CoMvelocities. In the subsequent time span the drop of liquid 1 comes nearly to rest with a minimal CoM-velocity of v 1 ≈ 2 · 10 -4 at t ≈ 2 6 . · 10 4 while the drop of liquid 2 slides over it till the compound drop reaches the converged stationary sliding state that has a velocity of approximately 2 9 . · 10 -3 as obtained by a linear fit for the time interval bounded by liquid 1's CoM passing x = 0 and x = 1000 , and by averaging the CoM-velocities of the two liquids.

Inspecting next the time evolution in Fig. 6 (c) one notices that also the initial symmetric drop-on-drop configuration at V /V 2 1 = 1 becomes asymmetric, namely, the upper drop slides faster as it is 'lubricated' by the lower drop. Ultimately, also here, the configuration changes into a stationary sliding asymmetric compound drop in drop-beside-drop configuration. In contrast to the previous case, now the crossing of contact regions occurs at the front end, i.e., the liquid 1-liquid 2 interface becomes strongly tilted, a liquid 1-liquid 2-substrate contact emerges, the leading liquid 1-liquid 2-gas becomes a liquid 2-gassubstrate contact, and a liquid 2-substrate interface emerges. The process results in a local maximum of both CoM velocities at t ≈ 2 0 . · 10 5 . The converged stationary sliding drop has a velocity of 1 6 . · 10 -3 , i.e., it is slower than the one in Fig. 6 (a).

## C. Coarsening ensemble of drops

Having discussed individual compound drops on horizontal and inclined homogeneous substrates, next, we consider the coarsening behavior of drop ensembles on horizontal substrates. We set the interface energies to ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 8 . , 0 9 . , 0 7 . , 1 4) . while the remaining parameters are as in section IV A. The resulting equilibrium contact angles are ( θ 12 , θ 13 , θ 23 ) ≈ (45 1 . ◦ , 38 9 . ◦ , 38 9 . ◦ ) .

First, we consider a minimal example of coarsening, namely, a single initial compound drop consisting of a large and a small

FIG. 7. Shown is a minimal example of coarsening dynamics starting with a compound drop consisting of three subdrops (two of liquid 2 beside one of liquid 1). Panels (a) and (b) show selected snapshots of the profiles at times t marked by dashed horizontal lines in the space-time plot of panel (c). The latter is given in log-normal form with liquid 1 and liquid 2 indicated by red and transparent blue color to visualize relative positions and layering. Only thicknesses above the threshold value of 1 2 . are shown. The interface energies are ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 8 . , 0 9 . , 0 7 . , 1 4) . , and the domain size is L = 1600 , note that this is only partly shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-a59dcc01-71b6-4e2c-891a-a9b179cd81f0.jpg)

volume of liquid 2 to the left and right of a single intermediate-size volume of liquid 1, respectively, i.e., the compound drop consists of three subdrops, see Fig. 7. This is similar to states visible in Fig. 8 (d) of Ref. [70] or the lower panel in Fig. 23 of Ref. [71] where the phase decomposition of sessile drops of a mixture is considered. Here, coarsening occurs as expected based on the greater Laplace pressure in the smaller subdrop, i.e., liquid 2 is transported from the right to the left via the adsorption layer of liquid 2 on liquid 1. During this process of Oswald ripening (coarsening by mass transfer), interestingly, the whole compound drop moves by nearly one drop base length to the left. This is well visible in the space-time plot (log-normal representation) in Fig. 7 (c). The resulting equilibrated compound drop corresponds to an asymmetric drop-beside-drop configuration similar to the states in the top and central panels in Fig. 23 of Ref. [71].

Next, we consider a larger ensemble of compound drops: Two initially homogeneous layers of thicknesses ¯ h 1 = ¯ h 2 = 4 that are linearly unstable with respect to dewetting are perturbed by noise of small amplitude. This results in spinodal dewetting as shown in Fig. 8 similar to, e.g., [4-7, 56]. Initially, for the chosen domain size fluctuations grow into seven compound drops. Each is in a nearly symmetric drop-on-drop configuration, see Fig. 8 (a). Interestingly, before drops coarsen, each individual drop undergoes a morphological transition from the symmetric drop-on-drop configuration to an asymmetric drop-beside-drop configuration, see Fig. 8 (b). This is unlike the processes shown in [5, 7]. Subsequently, coarsening sets in and is normally accompanied by the above discussed lateral motion, see the corresponding space-time plot in Fig. 8 (e). Indeed, it shows that coarsening by translation is dominant while coarsening by mass transfer only scarcely occurs at later times. Translation seems to be of greater importance than in the coarsening of drops of a single simple liquid on rigid solid substrates [72-74]. This reminds coarsening on soft solid substrates where in an intermediate softness range the translation mode normally dominates [75]. Fig. 8 (f) gives the corresponding dependence of the total free energy on time. It shows the typical sequence of step-like transitions between plateaus [5]. This can be explained by the observation that transient drop states also correspond to steady states that are however linearly unstable. The plateaus correspond to linear phases where the system slowly leaves an unstable steady state (exponentially slow but accelerating) till strong nonlinear behavior controls the final collapse of one droplet finalizing mass transfer or the coalescence of drops finalizing translation. Ultimately, a single large compound drop in drop-beside-drop configuration remains (Fig. 8 (d)). It represents either the state of lowest energy or a metastable state that corresponds to a deep local minimum of the energy.

The coarsening pathway and emerging final state strongly depend on the system parameters, in particular the energies and mean layer thicknesses. Considering the case of ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 4 1 6 0 8 1 1) . , . , . , . and ¯ h 1 = 9 and ¯ h 2 = 3 , we show in Fig. 9 a scenario alternative to Fig. 8. The thinner upper layer results in a fast dewetting of the upper layer on the liquid substrate formed by the lower layer that is however deformed by the process, see Fig. 9 (a). The resulting many liquid lenses coarsen from nearly 30 to about 10 till dewetting of the lower layer is in full process at log t ≈ 4 5 . , see Fig. 9 (e). At t ≈ 10 5 two compound drops have formed, both in nearly symmetric drop-on-drop configuration, see Fig. 9 (c). As strong asymmetry that promotes translation is missing, they coarsen via mass transfer. Interestingly this occurs in two steps: first, the left compound drop transfers its droplet of liquid 2 to the right compound drop. Second, the remaining droplet of liquid 1 slowly shrinks by mass transfer before it collapses at t ≈ 10 9 . In contrast to Fig. 8, here, the final state is a symmetric drop-on-drop configuration, see Fig. 9 (d).

Finally, we illustrate in Fig. 10 the usage of the model for studying the coarsening of compound drops on two-dimensional substrates, thereby using the same parameters as in Fig. 8. After a short dewetting phase with a labyrinthine stripe-like structure (Fig. 10 (a), also cf. Fig. 8 of Ref. [43]), compound droplets emerge and start to coarsen (Fig. 10 (b)). Due to the large number of

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-1cf312d0-a6fb-409a-8760-1366f4899751.jpg)

F

FIG. 8. Shown is the coarsening dynamics of a small ensemble of compound drops. Panels (a)-(d) show selected snapshots at times t marked by dashed horizontal lines in the log-normal space-time plot of panel (e), and panel (f) gives the corresponding dependence of the total free energy on time. The mean thicknesses of the liquid layers are ¯ h 1 = 4 and ¯ h 2 = 4 and the domain size is L = 1000 . All remaining parameters are as in Fig. 7. Also, consider the movie (file number 04) in the Supplementary Material [65].

droplets, individual coarsening events can not be discerned in the smooth decrease of the total energy with advancing time, see Fig. 10 (e). Careful inspection of the intermediate panel Fig. 10 (b) already shows the prevalence of asymmetric drop-beside-drop configurations. However, also other (transient) configurations can be spotted: compound drops with more than two subdrops (two or three drops of liquid 1 covered by or adjacent to a larger drop of liquid 2), drop-on-ring structures mostly with partial rings unlike more symmetric transient states shown in Fig. 10 of [70]. These complex configurations have mostly disappeared in Fig. 10 (c) where drop-beside-drop configurations dominate. When only a few drops are left they mostly correspond to the ultimately surviving drop-beside-drop configuration, see Fig. 10 (d). This agrees with the case of the one-dimensional substrate in Fig. 8. Note that due to the large numerical effort caused by exponentially increasing coarsening times we do not follow the process till the state of minimal energy is reached.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-d4280a54-6aa9-4abc-be86-ede5fd0e909b.jpg)

F

FIG. 9. Shown is the coarsening dynamics at interface energies ( γ 13 , γ 23 , γ s2 , γ s3 ) = (1 4 1 6 0 8 1 1) . , . , . , . , mean thicknesses ¯ h 1 = 9 and ¯ h 2 = 3 , and domain size L = 1000 . Panels, remaining parameters and other details are as in Fig. 8. Also, consider the movie (file number 05) in the Supplementary Material [65].

FIG. 10. Shown is the coarsening dynamics for a large ensemble of compound drops on a two-dimensional substrate. Panels (a) to (d) present selected snapshots at times t marked by dashed vertical lines in the log-normal dependence of total energy on time given in panel (e). Mean thicknesses are ¯ h 1 = 2 and ¯ h 2 = 2 , and the domain size is L × L with L = 1000 . The remaining parameters are as in Fig. 7. For the visualization of the snapshots we use PYVISTA [76]. Also, consider the Supplementary Material [65] for more snapshots (file number 06) and movies (file numbers 07-11).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00820v2-figure-e625840a-20b9-44e7-90cb-4562a63b3483.jpg)

## V. CONCLUSION

̸

We have revisited the mesoscopic hydrodynamic description of the dynamics of sessile compound drops consisting of two immiscible nonvolatile partially wetting liquids on horizontal and inclined homogeneous smooth rigid solid substrates. We have discussed that there exist two formulations (in terms of either layer thickness or interface height profiles) that transparently correspond to a (thermodynamic) gradient dynamics form. A third one that mixes the two previous formulations should be avoided as it obfuscates the gradient dynamics form and complicates a linear stability analysis because the resulting mobility matrix is asymmetric. The subsequent discussion of the wetting energy has clarified that in contrast to common praxis three contributions are needed (interaction across lower layer, interaction across upper layer, and interaction across both layers - resulting in three Hamaker constants) to be able to reflect the full spectrum of macroscopic parameters for partially wetting liquids in the mesoscopic model. A condition for general wetting energies, g h , h ( 1 2 ) = g 1 ( h 1 ) + g 2 ( h 2 ) , has also been given. Based on these considerations, we have determined corresponding consistency relations between macroscopic and mesoscopic parameters that ensure the full equivalence of all mesoscopic Laplace, Neumann and Young laws with their macroscopic counterparts. They are most conveniently established by employing the full-curvature variant of the model but fully carry over to the long-wave variant. In this way, we have corrected a number of oversights in some previous models that were intended to capture the full parameter spectrum for partially wetting liquids but did not. Here we remark that our focus on partially wetting liquids mesoscopically described by a wetting energy that consists of three terms of identical functional form allows for transparent consistency conditions between mesoscopic and macroscopic description. These considerations have to be amended if a wider class of wetting energies is considered. This includes, e.g., situations where one of the liquids is wetting or situations where first order wetting transitions can occur.

The resulting dynamic model now allows one to consider dynamic processes involving two-layer films and compound drops consisting of two immiscible partially wetting liquids on solid substrates. Our approach ensures that these processes ultimately converge to equilibrium states that can be exactly predicted based on macroscopic parameters. It can also be employed to investigate how the behavior is influenced by mesoscale effects when considering very small drops. Here, we have considered three examples: the spreading of individual compound drops on a one-dimensional horizontal substrate, sliding compound drops on one-dimensional inclined substrates (by incorporating the 'downhill-component' of gravity into the model), and the coarsening of drop ensembles on one- and two-dimensional horizontal substrates.

First, we have investigated the spreading of compound drops employing a symmetric initial configuration where a drop of liquid 2 completely covers a drop of liquid 1 (drop-covers-drop configuration). We have found that even at fixed interface energies, the dynamic behavior may dramatically depend on the volume ratio of the two liquids and the initial contact angles: The outer drop may spread faster than the inner one resulting in no change of configuration. However, if the inner drop spreads faster the compound drop can change configuration. When the inner contact line reaches the vicinity of the outer one, the latter can either slightly recede before a Neumann region develops and the configuration becomes one where a drop of liquid 2 only covers the central top part of the drop of liquid 1 (drop-on-drop configuration), or the contact lines meet forming a fourphase contact. This four-phase configuration can represent the final stable state or only form a short or long transient intermediate configuration before a drop-on-drop configuration emerges. The discussed transitions all occur between configurations discussed as possible equilibrium states in Refs. [5, 19-22, 24]. Note, in particular, Ref. [19] for the four-phase configuration that our simulations consistently confirm as the final equilibrium state in a finite range of control parameters. This seems to indicate that morphological phase diagrams as in Ref. [20] might need to be reconsidered in the relevant parameter ranges. Note that the gradient dynamics formulation has allowed us to study the importance of the different dissipation channels and corresponding changes at morphological transitions.

Second, we have considered drops sliding on a homogeneous incline. Also here, transitions to other configurations may occur: Obviously an initially symmetric drop-covers-drop configuration will become asymmetric due to the down-hill driving force that breaks the left-right symmetry of the system. However, we have not encountered an asymmetric drop-covers-drop configuration as stationary state, i.e., a compound drop that slides with constant speed and shape. Instead, the final configuration might be a drop of liquid 1 adjacent in front of a drop of liquid 2 or vice versa (drop-beside-drop configuration). Also an initial drop-on-drop configuration will transition into a drop-beside-drop configuration. Here, a more extensive parametric study will most likely uncover more possible transitions. Note that also here the gradient dynamics formulation allows one to study the influence of the different dissipation channels and their dependence on parameters. We have left this as a task for the future.

Finally, we have returned to horizontal substrates and studied the coarsening of droplet ensembles on oneand twodimensional substrates. We find that the two basic coarsening models of mass transfer between drops and translation of drops both contribute. Their relative importance depends not only on drop size but also on configuration, an effect that requires further future scrutiny. Note that also in the two-dimensional case, asymmetric drops do not only appear as long-lived transients but may also feature as final equilibrium configuration - in agreement with [20, 21], termed 'Janus drops' in [20]. Note that the existence of such states is negated in Ref. [19] based on a torque balance argument that, however, does not apply in situations without an external force field like gravity. In [5, 43] all drop-like long-lived transients and final states are symmetric while in [24] asymmetric intermediate states are found. Also here future studies should expand on [24, 77] and provide morphological phase diagrams that clarify the ranges of existence and stability of such states in terms of physical control parameters. We are

not aware of recent experimental work on large-scale dewetting and subsequent coarsening phenomena for bi- or multi-liquid systems of partially wetting liquids beyond the already older works on decomposing and dewetting films reviewed in [31, 40]. We hope our study motivates further experiments on dewetting phenomena with a focus on complex liquids and multi-liquid set-ups. An interesting avenue could be the pulsed-laser-induced dewetting of multilayer metal films similar to the single-layer set-ups reviewed in [78]. Also see [70] for decomposing dewetting nanoscopic films of alloys.

Furthermore, our gradient dynamics approach lends itself to a number of future amendments and extensions. Beside the relatively simple incorporation of topographically and chemically heterogeneous substrates like for simple liquids done in [79, 80], e.g., to investigate the stick-slip motion [81] of compound drops, one may in a straightforward manner incorporate hydrostatic effects including different mass densities for the liquids. Beside gravity or electrical fields [14, 82, 83] one might incorporate volatility [84], soft elastic substrates [85], or surfactants [47]. Although the systems then become rather complicated, employing the gradient dynamics approach to mesoscopic hydrodynamics [46] allows for the combination of building blocks known from simpler systems. However, one has to go beyond the gradient dynamics framework if thermocapillary effects shall be added, e.g., to study the influence of internal convection and long-wave thermal Marangoni instabilities [86, 87]. For two-layer films such influences were incorporated into the two-layer model (but not employed) in [5], for further studies see [15-18]. Here, building on the present model would allow one to assess how thermocapillarity influences Young and Neumann angles.

## ACKNOWLEDGMENTS

The authors acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) via Grant No. TH781/12-1 and TH781/122 within Priority Program (SPP) 2171 'Dynamic Wetting of Flexible, Adaptive, and Switchable Surfaces'; We acknowledge preliminary work of Kevin Mitas on earlier models for two-layer systems published in his PhD thesis. Furthermore, we acknowledge fruitful discussions with Daniel Greve, Simon Hartmann, Christopher Henkel and Dominik Thy at the University of M¨ unster, and with participants of the events organized by SPP 2171. Part of the calculations were performed on the HPC cluster PALMA II of the University of M¨ unster, subsidized by the DFG via Grant INST 211/667-1.

This version of the article has been accepted for publication, after peer review but is not the Version of Record and does not reflect all post-acceptance improvements and corrections. The Version of Record is available online at 10.1103/PhysRevFluids.10.024002.

## DATA AVAILABILITY STATEMENT

The underlying data and source code required to reproduce the shown results is publicly available at the data repository zenodo with the corresponding DOI:10.5281/zenodo.13143462.

## Appendix A: Macroscopic description

For convenience of the reader, here, we review the macroscopic description of static compound drops on rigid solid substrates. In particular, we introduce the energy functional F macro based on the six interface energies (liquid 1 to fluid 3 (here ambient gas, but may also be an ambient liquid) γ 13 , liquid 2 to fluid 3 γ 23 , liquid 1 to liquid 2 γ 12 , liquid 1 to substrate γ s1 , liquid 2 to substrate γ s2 , and fluid 3 to substrate γ s3 ), and summarize how minimization of the total energy results in Laplace, Neumann, and Young laws and point out relations between them. This forms the basis of establishing the consistency conditions discussed in section III B.

To discuss the energy, we consider a drop-on-drop configuration consisting of a drop of liquid 2 (of volume V 2 and lateral extension 2 L 2 ) on top of a drop of liquid 1 (of volume V 1 and lateral extension 2 L 1 ) on a one-dimensional rigid smooth solid substrate as sketched in Fig. 1 (a). The configuration is symmetric and is embedded in an ambient fluid 3. The configuration features five of the six possible interfaces, indicated by the respective interface energies γ ij , e.g., γ 12 at the interface 1-2. The three static thickness profiles h 12 , h 13, and h 23 correspond to circular arcs. An alternative drop-covers-drop configuration (with a much larger V 2 ) would also feature the sixth energy γ s2 (cf. Fig. 1(c)).

and ensuring geometric consistency and

$$\delta \mathcal { F } _ { \text{macro} } & = \left ( \gamma _ { 1 3 } \frac { \partial _ { x } h _ { 1 3 } } { \xi _ { 1 3 } } + \lambda _ { 1 } \right ) \delta h _ { 1 3 } \left ( L _ { 1 } \right ) - \left ( \gamma _ { 1 3 } \frac { \partial _ { x } h _ { 1 3 } } { \xi _ { 1 3 } } + \lambda _ { 3 } \right ) \delta h _ { 1 3 } \left ( L _ { 2 } \right ) \\ & \quad - \int _ { L _ { 2 } } ^ { L _ { 1 } } \left ( \gamma _ { 1 3 } \frac { \partial _ { x x } h _ { 1 3 } } { \xi _ { 1 3 } ^ { 3 } } + P _ { 1 } \right ) \delta h _ { 1 3 } ( x ) \text{d} x, \\ \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{d} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad  \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{d} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad  \text{e} \ \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{t} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad\text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e}\quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text` \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \tilde { \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{u} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad \text{e} \quad$$

respectively, while minimization w.r.t. the contact line positions results in

$$\delta \mathcal { F } _ { \text{macro} } = ( \gamma _ { 1 3 } \xi _ { 1 3 } + \gamma _ { s 1 } - \gamma _ { s 3 } - P _ { 1 } h _ { 1 3 } + \lambda _ { 1 } \partial _ { x } h _ { 1 3 } ) \, \delta L _ { 1 } \text{ \ at \ } x = L _ { 1 } \text{ \quad } ( A 1 0 )$$

and

$$\begin{array} { c } \delta \mathcal { F } _ { \text{macro} } = & ( \gamma _ { 1 2 } \xi _ { 1 2 } + \gamma _ { 2 3 } \xi _ { 2 3 } - \gamma _ { 1 3 } \xi _ { 1 3 } - P _ { 2 } \left ( h _ { 2 3 } - h _ { 1 2 } \right ) - P _ { 1 } \left ( h _ { 1 2 } - h _ { 1 3 } \right ) \\ & + \lambda _ { 2 } \left ( \partial _ { x } h _ { 2 3 } - \partial _ { x } h _ { 1 2 } \right ) + \lambda _ { 3 } \left ( \partial _ { x } h _ { 2 3 } - \partial _ { x } h _ { 1 3 } \right ) \right ) \delta L _ { 2 } \text{ at } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \end{array}$$

We identify the expression ∂ xx h ij ξ 3 ij = κ ij as the local curvature of h ij ( x ) .

## 1. Energy functional

The macroscopic energy functional for a compound drop consisting of two immiscible liquids for the configuration in Fig. 1(a) is given by

$$\cdot \varepsilon \cdot \dots \cdot \cdot \cdot \\ \mathcal { F } _ { \text{macro} } & = \int _ { 0 } ^ { L _ { 1 } } \gamma _ { s 1 } \mathrm d x + \int _ { 0 } ^ { L _ { 2 } } ( \gamma _ { 1 2 } \xi _ { 1 2 } + \gamma _ { 2 3 } \xi _ { 2 3 } ) \mathrm d x + \int _ { L _ { 2 } } ^ { L _ { 1 } } \gamma _ { 1 3 } \xi _ { 1 3 } \mathrm d x + \int _ { L _ { 1 } } ^ { D } \gamma _ { s 3 } \mathrm d x \\ & \quad - P _ { 1 } \left ( \int _ { 0 } ^ { L _ { 2 } } h _ { 1 2 } \mathrm d x + \int _ { L _ { 2 } } ^ { L _ { 1 } } h _ { 1 3 } \mathrm d x - V _ { 1 } \right ) - P _ { 2 } \left ( \int _ { 0 } ^ { L _ { 2 } } ( h _ { 2 3 } - h _ { 1 2 } ) \mathrm d x - V _ { 2 } \right ) \\ & \quad \pm \vee h \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \, \vee \, \vec { h } \dots \, \vec { r } \, \vee \cdot \end{cases}$$

+ λ h 1 13 ( L 1 ) + λ 2 ( h 23 ( L 2 ) -h 12 ( L 2 )) + λ 3 ( h 23 ( L 2 ) -h 13 ( L 2 )) , where due to the reflection symmetry we only consider half the domain, i.e., x ∈ [0 , D ] , assume that all interface profiles can be written in Monge representation h ij ( x ) , and use the one-dimensional equivalents of the metric factors in the mesoscopic energy (10), i.e., ξ ij = √ 1 + ( d h ij d x ) 2 . The constants P , P 1 2 , λ 1 , λ 2 and λ 3 represent the Lagrange multipliers for various constraints, namely, controlling the volumes of liquids 1 and 2

$$\int _ { 0 } ^ { L _ { 2 } } h _ { 1 2 } \text{d} x + \int _ { L _ { 2 } } ^ { L _ { 1 } } h _ { 1 3 } \text{d} x = V _ { 1 },$$

$$\int _ { 0 } ^ { L _ { 2 } } ( h _ { 2 3 } - h _ { 1 2 } ) \mathrm d x = V _ { 2 },$$

$$h _ { 1 3 } \left ( L _ { 1 } \right ) & = 0, & \left ( A 4 \right )$$

$$h _ { 1 2 } \left ( L _ { 2 } \right ) = h _ { 2 3 } \left ( L _ { 2 } \right ),$$

$$h _ { 1 3 } \left ( L _ { 2 } \right ) = h _ { 2 3 } \left ( L _ { 2 } \right ),$$

respectively. Next, the energy (A1) is minimized w.r.t. variations in the profiles h ij ( x ) and the positions of the contact lines L 1 and L 2 .

## 2. Minimization of F macro

Minimization of the energy functional (A1) w.r.t. the profiles h 12 , h 23 and h 13 gives

$$\delta \mathcal { F } _ { \text{macro} } = \left ( \gamma _ { 1 2 } \frac { \partial _ { x } h _ { 1 2 } } { \xi _ { 1 2 } } - \lambda _ { 2 } \right ) \delta h _ { 1 2 } \left ( L _ { 2 } \right ) - \int _ { 0 } ^ { L _ { 2 } } \left ( \gamma _ { 1 2 } \frac { \partial _ { x x } h _ { 1 2 } } { \xi _ { 1 2 } ^ { 3 } } + \left ( P _ { 1 } - P _ { 2 } \right ) \right ) \delta h _ { 1 2 } ( x ) \mathrm d x,$$

$$\delta \mathcal { F } _ { \text{macro} } = \left ( \gamma _ { 2 3 } \frac { \partial _ { x } h _ { 2 3 } } { \xi _ { 2 3 } } + \lambda _ { 2 } + \lambda _ { 3 } \right ) \delta h _ { 2 3 } \left ( L _ { 2 } \right ) - \int _ { 0 } ^ { L _ { 2 } } \left ( \gamma _ { 2 3 } \frac { \partial _ { x x } h _ { 2 3 } } { \xi _ { 2 3 } ^ { 3 } } + P _ { 2 } \right ) \delta h _ { 2 3 } ( x ) \text{d} x,$$

## 3. Laplace, Neumann and Young laws

We then obtain from Eqs. (A7) to (A9) the Laplace laws

γ

12

κ

12

=

P

2

P

1

γ

-

23

κ

23

=

P

2

γ

-

13

κ

13

=

-

P

1

$$& \text{for } x \in [ 0, L _ { 2 } ] & & ( \text{A12} ) \\ & \text{for } r \in [ \mathfrak { m } \ I _ { \cap } ] & & \text{for } \mathbb { m } \ I \mathbb { m }$$

$$\text{for } x \in [ L _ { 2 }, L _ { 1 } ] \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text$$

and the relations for the Lagrange multipliers

$$\gamma _ { 1 2 } \frac { \partial _ { x } h _ { 1 2 } } { \xi _ { 1 2 } } & = \lambda _ { 2 } & \text{ at } x = L _ { 2 } \,, & \text{(A15)}$$

$$\gamma _ { 2 3 } \frac { \partial _ { x } ^ { \mu } h _ { 2 3 } } { \xi _ { 2 3 } } = - \lambda _ { 2 } - \lambda _ { 3 } & \text{ at } x = L _ { 2 } \,,$$

$$\gamma _ { 1 3 } \frac { \partial _ { x } ^ { \dots } \bar { h } _ { 1 3 } } { \xi _ { 1 3 } } = - \lambda _ { 3 } & \text{ at } x = L _ { 2 } & \text{(A17)}$$

$$\gamma _ { 1 3 } \frac { \partial _ { x } ^ { \dots } h _ { 1 3 } } { \xi _ { 1 3 } } = - \lambda _ { 1 } & \text{ at } x = L _ { 1 } \,.$$

The Laplace laws state that the interfaces have constant curvatures κ ij , i.e., circular arcs in the considered geometry. Thereby, the κ ij are determined by the pressure jumps across the corresponding interface (identifying P 1 and P 2 as the pressure in liquid 1 and 2, respectively). Care has to be taken as different sign conventions exist. Combining Eqs. (A12) to (A14) gives the relation

$$0 & = - \gamma _ { 1 3 } \kappa _ { 1 3 } + \gamma _ { 1 2 } \kappa _ { 1 2 } + \gamma _ { 2 3 } \kappa _ { 2 3 }. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &   & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & . & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

Next, we consider the three-phase contact of the three fluids 1, 2, and 3 as shown in Fig. 1 (a), i.e., Neumann angles are defined as angles between the tangent to the respective interface at the contact point and the horizontal. We relate the angles to the interface slopes by ∂ h x ij | x = L 2 = ± tan( ψ ij ) the '+' is for interface 1-2, and the '-' is for interfaces 2-3 and 1-3, i.e., all angles shown in Fig. 1 (a) are defined positive.

The minimizations give Young's law for the three-phase contact of substrate, liquid 1, and fluid 3 at x = L 1 and Neumann's law at the three-fluid contact point at x = L 2 . Combining (A15), (A16) and (A17) gives

$$\gamma _ { 2 3 } \frac { \partial _ { x } h _ { 2 3 } } { \xi _ { 2 3 } } = - \gamma _ { 1 2 } \frac { \partial _ { x } h _ { 1 2 } } { \xi _ { 1 2 } } + \gamma _ { 1 3 } \frac { \partial _ { x } h _ { 1 3 } } { \xi _ { 1 3 } } \,.$$

Employing ξ -1 ij | x = L 2 = cos( ψ ij ) , this becomes the vertical component of Neumann's law

$$\gamma _ { 2 3 } \sin ( \psi _ { 2 3 } ) = \gamma _ { 1 2 } \sin ( \psi _ { 1 2 } ) + \gamma _ { 1 3$$

Further, we use Eq. (A11) in the limit P , P 1 2 → 0 (i.e., for large drops) together with Eqs. (A15) and (A17) to get

$$\begin{smallmatrix} \cdot \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot$$

$$+ \gamma _ { 1 2 } \sin ( \psi _ { 1 2 } ) ( - \tan ( \psi _ { 2 3 } ) - \tan ( \psi _ { 1 2 } ) )$$

Using Eq. (A21) and rearranging the trigonometric functions we obtain the horizontal component of Neumann's law

$$\gamma _ { 1 3 } \cos ( \psi _ { 1 3 } ) = & \gamma _ { 1 2 } \cos ( \psi _ { 1 2 } ) + \gamma _ { 2$$

Finally, Young's law is obtained rearranging Eqs. (A10) and (A18) into

$$\gamma _ { 1 3 } \cos ( \theta _ { 1 3 } ) = \gamma _ { s 3 } - \gamma _ { s 1 } \,. \quad \quad \quad$$

Considering other compound drop configurations than the one in Fig. 1(a), e.g., the one in Fig. 1(c), one obtains similar Young laws for contact points of the substrate with on the one hand liquid 2 and fluid 3 (cf. Fig. 1(c), outer contact line),

$$\gamma _ { 2 3 } \cos ( \theta _ { 2 3 } ) = \gamma _ { s 3 } - \gamma _ { s 2 } \,, \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots$$

and on the other hand with liquid 1 and liquid 2 (cf. Fig. 1(c), inner contact line),

$$\gamma _ { 1 2 } \cos ( \theta _ { 1 2 } ) = \gamma _ { s 2 } - \gamma _ { s 1 } \,. \\ \cdot \cdot \cdot \cdot \,.$$

Note that Eqs. (A24)-(A26) can be combined into

$$0 = - \gamma _ { 1 3 } \cos ( \theta _ { 1 3 } ) + \gamma _ { 1 2 } \cos ( \theta _ { 1 2 } ) + \gamma _ { 2 3 } \cos ( \theta _ { 2 3 } ) \,, & & ( A 2 7 ) \\ \therefore \, \, \cdots \, \cdots \, \cdots$$

a relation somewhat similar to (A19) and (A23).

$$& \text{for $x\in [0, L_{2}]$} \\ & \text{for $r\in [T_{\wedge} T_{\wedge} ]$}$$

and

$$v _ { 2 } = & \frac { 1 } { n _ { 2 } } \left \lceil \frac { z ^ { 2 } } { 2 } \partial _ { x } p _ { 2 } - z ( h _ { 1 } + h _ { 2 } ) \partial _ { x } p _ { 2 } + I \right \rceil$$

$$v _ { 2 } = & \frac { 1 } { \eta _ { 2 } } \left [ \frac { z ^ { 2 } } { 2 } \partial _ { x } p _ { 2 } - z ( h _ { 1 } + h _ { 2 } ) \partial _ { x } p _ { 2 } + I \right ] \\ w _ { 2 } = & - \frac { 1 } { \eta _ { 2 } } \left [ \frac { z ^ { 3 } } { 6 } \partial _ { x x } p _ { 2 } - \frac { z ^ { 2 } } { 2 } \partial _ { x } \left [ ( h _ { 1 } + h _ { 2 } ) \partial _ { x } p _ { 2 } \right ] + z \partial _ { x } I \\ & - \frac { h _ { 1 } ^ { 3 } } { 6 } \partial _ { x x } p _ { 2 } + \frac { h _ { 1 } ^ { 2 } } { 2 } \partial _ { x } \left [ ( h _ { 1 } + h _ { 2 } ) \partial _ { x } p _ { 2 } \right ] - h _ { 1 } \partial _ { x } I - \eta _ { 2 } w _ { 1 } | _ { z = h _ { 1 } } \right ] \,. \\. ^ { \cdot }$$

Here, I = h 2 1 2 ∂ x ( p 2 -p η /η 1 2 1 ) + (1 -η /η 2 1 ) h h ∂ p 1 2 x 2 , and the pressures p 1 and p 2 correspond to the ones in Eq. (4).

## 2. Dissipation channels

In order to analyze how dissipation occurs via the different channels, we consider the total time derivative of the free energy in Eq. (10), i.e., the negative of the total dissipation D and split it into parts D 1 and D 2 corresponding to liquid 1 and liqud 2,

## 4. Long-wave approximation

Finally, we consider the long-wave approximation of the macroscopic relations obtained in section A 3. To do so, we assume that all Young and Neumann angles are small, i.e., we approximate cos ϕ ≈ 1 -ϕ / 2 2 and sin ϕ ≈ ϕ . In this way, the components (A21) and (A23) of Neumann's law become

$$\gamma _ { 1 2 } \psi _ { 1 2 } + \gamma _ { 1 3 } \psi _ { 1 3 } = \gamma _ { 2 3 } \psi _ { 2 3 },$$

$$\text{and} \quad \gamma _ { 1 2 } \psi _ { 1 2 } ^ { 2 } + \gamma _ { 2 3 } \psi _ { 2 3 } ^ { 2 } - \gamma _ { 1 3 } \psi _ { 1 3 } ^ { 2 } = 2 ( \gamma _ { 1 2 } + \gamma _ { 2 3 } - \gamma _ { 1 3 } ) \,,$$

respectively. The Young laws (A24), (A25) and (A26) become

$$\gamma _ { 1 3 } \theta _ { 1 3 } ^ { 2 } & = 2 ( \gamma _ { s 1 } + \gamma _ { 1 3 } - \gamma _ { s 3 } ) \,, \\ \omega _ { 1 } \, \ a ^ { 2 } \, = \, \alpha / \omega _ { 1 } \, + \, \omega _ { 1 } \, \quad \omega _ { 1 } \, \quad \omega _ { 1 } \, \quad \omega _ { 1 } \, \quad$$

$$\text{and} \quad \gamma _ { 1 2 } \theta _ { 1 2 } ^ { 2 } = 2 ( \gamma _ { s 1 } + \gamma _ { 1 2 } - \gamma _ { s 2 } ) \,,$$

$$\gamma _ { 2 3 } \theta _ { 2 3 } ^ { 2 } & = 2 ( \gamma _ { s 2 } + \gamma _ { 2 3 } - \gamma _ { s 3 } ) \,, \\ \tilde { \tau } \tilde { \Delta } \tilde { \Delta } ^ { 2 } & = \tilde { \Delta } ( \tilde { \Delta$$

respectively. By combining them one obtains the relation

$$\gamma _ { 1 2 } \theta _ { 1 2 } ^ { 2 } + \gamma _ { 2 3 } \theta _ { 2 3 } ^ { 2 } - \gamma _ { 1 3 } \theta _ { 1 3 } ^ { 2 } = 2 ( \gamma _ { 1 2 } + \gamma _ { 2 3 } - \gamma _ { 1 3 } )$$

again a structure very similar to the horizontal component of Neumann's law (A29).

## Appendix B: Related physical quantities for the mesoscopic model

## 1. Velocitiy fields

To reconstruct the vertically resolved velocity fields v 1 = ( v , w 1 1 ) T and v 2 = ( v , w 2 2 ) T of the lower and upper liquid layer, respectively, from the film height profiles one may perform the standard derivation of the thin-film equation for the two-layer system [5, 33, 39, 40]. Explicitly, one uses no-slip and no-penetration boundary conditions at the solid substrate, continuity of shear stress and velocity at the liquid 1-liquid 2 interface, and zero shear stress and the kinematic boundary condition at the liquid 2-gas interface. Using the formulation in terms of layer thicknesses, one ultimately obtains the form in Eq. (2). However, as an intermediate result, one directly obtains the height-resolved x -components of the velocities v i ( x, z, t ) in the two liquid layers. The z -components one obtains using the incompressibility condition ∂ v x i + ∂ w z i = 0 and integration in z (similar to Ref. [66]) to finally obtain

$$v _ { 1 } = & \frac { 1 } { \eta _ { 1 } } \left [ \frac { z ^ { 2 } } { 2 } \partial _ { x } p _ { 1 } - z ( h _ { 1 } \partial _ { x } p _ { 1 } + h _ { 2 } \partial _ { x } p _ { 2 } ) \right ] \quad & \quad \quad \quad$$

$$w _ { 1 } = - \frac { \ddots } { 2 \eta _ { 1 } } \left [ \frac { z ^ { 3 } } { 3 } \partial _ { x x } p _ { 1 } - z ^ { 2 } \partial _ { x } ( h _ { 1 } \partial _ { x } p _ { 1 } + h _ { 2 } \partial _ { x } p _ { 2 } )$$

respectively. The approach is similar to the one employed in [52, 69] for drops on adaptive substrates. Using the formulation in terms of layer thicknesses in Eqs. (1)-(2) allows us to directly determine D 1 and D 2 . For the one-dimensional domain we obtain

$$\frac { \mathrm d } { \mathrm d t } \mathcal { F } = \int _ { \Omega } \sum _ { \alpha = 1 } ^ { 2 } \frac { \delta \mathcal { F } } { \delta h _ { \alpha } } \partial _ { t } h _ { \alpha } \mathrm d x = \int _ { \Omega } \sum _ {$$

$$= - \int _ { \Omega } \left ( \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { 1 } } \right ) \left [ \sum _ { \beta = 1 } ^ { 2 } Q _ { 1 \beta } \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { \beta } } \right ] \mathrm d x - \int _ { \Omega } \left ( \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { 2 } } \right ) \left [ \sum _ { \beta = 1 } ^ { 2 } Q _ { 2 \beta } \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { \beta } } \right ] \mathrm d x \,.$$

$$\frac { \mathrm d } { \mathrm d t } \mathcal { F } & = \int _ { \Omega } \tilde { \varprojlim } \frac { \delta \mathcal { F } } { \delta h _ { \alpha } } \partial _ { t h _ { \alpha } } d x = \int _ { \Omega } \tilde { \varprojlim } \frac { \delta \mathcal { F } } { \delta h _ { \alpha } } ( - \partial _ { x } j _ { \alpha } ) d x = \tilde { \varprojlim } \int _ { \Omega } \left ( \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { \alpha } } \right ) j _ { \alpha } \mathrm d x \\ & = - \int _ { \Omega } \left ( \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { 1 } } \right ) \left [ \sum _ { \lambda _ { - 1 } } ^ { 2 } Q _ { 1 \beta } \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { \beta } } \right ] \mathrm d x - \int _ { \Omega } \left ( \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { 2 } } \right ) \left [ \sum _ { \lambda _ { - 1 } } ^ { 2 } Q _ { 2 \beta } \partial _ { x } \frac { \delta \mathcal { F } } { \delta h _ { \beta } } \right ] \mathrm d x.$$

In the final transformation of Eq. (B5) we have used partial integration and eliminated the surface terms by assuming periodic boundary conditions. The negative of the first and second term on the right hand side of (B6) represent the two dissipation channels, namely, the overall viscous dissipation D 1 and D 2 in liquid 1 and in liquid 2, respectively.

- [1] F. Brochard-Wyart, P. Martin, and C. Redon, Liquid/liquid dewetting, Langmuir 9 , 3682 (1993).
- [2] A. Faldi, R. J. Composto, and K. I. Winey, Unstable polymer bilayers. 1. Morphology of dewetting, Langmuir 11 , 4855 (1995).
- [3] M. Sferrazza, M. Heppenstall-Butler, R. Cubitt, D. Bucknall, J. Webster, and R. A. L. Jones, Interfacial instability driven by dispersive forces: The early stages of spinodal dewetting of a thin polymer film on a polymer substrate, Phys. Rev. Lett. 81 , 5173 (1998).
- [4] A. Pototsky, M. Bestehorn, D. Merkt, and U. Thiele, Alternative pathways of dewetting for a thin liquid two-layer film, Phys. Rev. E 70 , 025201(R) (2004).
- [5] 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).
- [6] D. Bandyopadhyay, R. Gulabani, and A. Sharma, Stability and dynamics of bilayers, Ind. Eng. Chem. Res. 44 , 1259 (2005).
- [7] D. Bandyopadhyay and A. Sharma, Nonlinear instabilities and pathways of rupture in thin liquid bilayers, J. Chem. Phys. 125 , 054711 (2006).
- [8] S. Yadavali, H. Krishna, and R. Kalyanaraman, Morphology transitions in bilayer spinodal dewetting systems, Phys. Rev. B 85 , 235446 (2012).
- [9] 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).
- [10] D. Peschka, S. Bommer, S. Jachalski, R. Seemann, and B. Wagner, Impact of energy dissipation on interface shapes and on rates for dewetting from liquid substrates, Sci. Rep. 8 , 13295 (2018).
- [11] R. L¨ udtge, Ueber die Ausbreitung der Fl¨ssigkeiten auf einander, Ann. Phys. (Poggendorf) u 137 , 362 (1869).
- [12] R. V. Craster and O. K. Matar, On the dynamics of liquid lenses, J. Colloid Interface Sci. 303 , 503 (2006).
- [13] L. G. Badratinova, Non-linear long-wave equations in the thermocapillary convection problem for a two-layered liquid, Inst. Gifrodinamiki, So Akad. Nauk SSSR, Novosibirsk 69 , 3 (1985).
- [14] D. Bandyopadhyay and A. Sharma, Electric field induced instabilities in thin confined bilayers, J. Colloid Interface Sci. 311 , 595 (2007).
- [15] A. A. Nepomnyashchy and I. B. Simanovskii, Marangoni instability in ultrathin two-layer films, Phys. Fluids 19 , 122103 (2007).
- [16] A. Nepomnyashchy and I. Simanovskii, Instabilities and ordered patterns in nonisothermal ultrathin bilayer fluid films, Phys. Rev. Lett. 102 , 164501 (2009).

[17]

A. A. Nepomnyashchy and I. B. Simanovskii, Novel criteria for the development of monotonic and oscillatory instabilities in a two-layer film, Phys. Fluids

29

, 092104 (2017).

- [18] A. Nepomnyashchy and I. Simanovskii, Droplets on the liquid substrate: Thermocapillary oscillatory instability, Phys. Rev. Fluids 6 , 034001 (2021).
- [19] L. Mahadevan, M. Adda-Bedia, and Y. Pomeau, Four-phase merging in sessile compound drops, J. Fluid Mech. 451 , 411 (2002).
- [20] M. J. Neeson, R. F. Tabor, F. Grieser, R. R. Dagastine, and D. Y. C. Chan, Compound sessile drops, Soft Matter 8 , 11042 (2012).
- [21] M. Ben Said, M. Selzer, B. Nestler, D. Braun, C. Greiner, and H. Garcke, A phase-field approach for wetting phenomena of multiphase droplets on solid surfaces, Langmuir 30 , 4033 (2014).
- [22] Y. Zhang, D. Chatain, S. L. Anna, and S. Garoff, Stability of a compound sessile drop at the axisymmetric configuration, J. Colloid Interface Sci. 462 , 88 (2016).
- [23] R. Iqbal, S. Dhiman, A. K. Sen, and A. Q. Shen, Dynamics of a water droplet over a sessile oil droplet: Compound droplets satisfying a Neumann condition, Langmuir 33 , 5713 (2017).
- [24] G. Kitavtsev, Composite solutions to a liquid bilayer model, Journal of Nonlinear Science 35 , 10.1007/s00332-024-10108-5 (2024).
- [25] V. Bertin, C. L. Lee, T. Salez, E. Raphael, and K. Dalnoki-Veress, Capillary levelling of immiscible bilayer films, J. Fluid Mech. 911 , A13 (2021).
- [26] K. D. Danov, V. N. Paunov, N. Alleborn, H. Raszillier, and F. Durst, Stability of evaporating two-layered liquid film in the presence of surfactant - I. The equations of lubrication approximation, Chem. Eng. Sci. 53 , 2809 (1998).
- [27] K. D. Danov, V. N. Paunov, S. D. Stoyanov, N. Alleborn, H. Raszillier, and F. Durst, Stability of evaporating two-layered liquid film in the presence of surfactant - II. Linear analysis, Chem. Eng. Sci. 53 , 2823 (1998).
- [28] V. N. Paunov, K. D. Danov, N. Alleborn, H. Raszillier, and F. Durst, Stability of evaporating two-layered liquid film in the presence of surfactant - III. Non-linear stability analysis, Chem. Eng. Sci. 53 , 2839 (1998).

- [29] J. J. Kriegsmann and M. J. Miksis, Steady motion of a drop along a liquid interface, SIAM J. Appl. Math. 64 , 18 (2003).
- [30] R. Mukherjee and A. Sharma, Instability, self-organization and pattern formation in thin soft films, Soft Matter 11 , 8717 (2015).
- [31] M. Geoghegan and G. Krausch, Wetting at polymer surfaces and interfaces, Prog. Polym. Sci. 28 , 261 (2003).
- [32] L. V. Govor, G. Reiter, G. H. Bauer, and J. Parisi, Bilayer formation in thin films of a binary solution, Phys. Lett. A 353 , 198 (2006).
- [33] S. Jachalski, D. Peschka, A. M¨nch, and B. Wagner, Impact of interfacial slip on the stability of liquid two-layer polymer films, J. Eng. u Math. 86 , 9 (2014).
- [34] L. Xu, D. Bandyopadhyay, P. D. S. Reddy, A. Sharma, and S. W. Joo, Giant slip induced anomalous dewetting of an ultrathin film on a viscous sublayer, Sci. Rep. 7 , 14776 (2017).
- [35] J. H. Guan, G. G. Wells, B. Xu, G. McHale, D. Wood, J. Martin, and S. Stuart-Cole, Evaporation of sessile droplets on slippery liquidinfused porous surfaces (slips), Langmuir 31 , 11781 (2015).
- [36] C. Semprebon, G. McHale, and H. Kusumaatmaja, Apparent contact angle and contact angle hysteresis on liquid infused surfaces, Soft Matter 13 , 101 (2017).
- [37] M. Tress, S. Karpitschka, P. Papadopoulos, J. H. Snoeijer, D. Vollmer, and H. J. Butt, Shape of a sessile drop on a flat surface covered with a liquid film, Soft Matter 13 , 3760 (2017).
- [38] P. Valluri, L. O. N´raigh, H. Ding, and P. D. M. Spelt, Linear and nonlinear spatio-temporal instability in laminar two-layer flows, J. Fluid ´ a Mech. 656 , 458 (2010).
- [39] A. Oron, S. H. Davis, and S. G. Bankoff, Long-scale evolution of thin liquid films, Rev. Mod. Phys. 69 , 931 (1997).
- [40] R. V. Craster and O. K. Matar, Dynamics and stability of thin liquid films, Rev. Mod. Phys. 81 , 1131 (2009).
- [41] M. H. Ward, Interfacial thin films rupture and self-similarity, Phys. Fluids 23 , 062105 (2011).
- [42] S. Jachalski, R. Huth, G. Kitavtsev, D. Peschka, and B. Wagner, Stationary solutions of liquid two-layer thin-film models, SIAM J. Appl. Math. 73 , 1183 (2013).
- [43] A. Pototsky, M. Bestehorn, D. Merkt, and U. Thiele, Evolution of three-dimensional interface patterns in dewetting two-layer liquid films, Europhys. Lett. 74 , 665 (2006).
- [44] J. H. Snoeijer, Free-surface flows with large slopes: Beyond lubrication theory, Phys. Fluids 18 , 021701 (2006).
- [45] 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).
- [46] U. Thiele, Recent advances in and future challenges for mesoscopic hydrodynamic modelling of complex wetting, Colloid Surf. A 553 , 487 (2018).
- [47] U. Thiele, A. J. Archer, and L. M. Pismen, Gradient dynamics models for liquid films with soluble surfactant, Phys. Rev. Fluids 1 , 083903 (2016).
- [48] R. Huth, S. Jachalski, G. Kitavtsev, and D. Peschka, Gradient flow perspective on thin-film bilayer flows, J. Eng. Math. 94 , 43 (2015).
- [49] F. E. Neumann and A. Wangerin, Vorlesungen ¨ uber die Theorie der Capillarit¨t, gehalten an der Universit¨t K¨nigsberg von Franz a a o Neumann , Vorlesungen ¨ber mathematische Physik. (7. Hft.) (B. G. Teubner, Leipzig, 1894). u
- [50] P. d. Bois-Reymond, De aequilibrio fluidorum , Ph.D. thesis, Humboldt Universit¨ at, Berlin (1859).
- [51] J. N. Israelachvili, Intermolecular and Surface Forces , 3rd ed. (Academic Press, San Diego, 2011).
- [52] D. Greve, S. Hartmann, and U. Thiele, Stick-slip dynamics in the forced wetting of polymer brushes, Soft Matter 19 , 4041 (2023).
- [53] 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.
- [54] X. Xu, U. Thiele, and T. Qian, A variational approach to thin film hydrodynamics of binary mixtures, J. Phys.: Condens. Matter 27 ,
- 085005 (2015).
- [55] K. Kostourou, D. Peschka, A. M¨ unch, B. Wagner, S. Herminghaus, and R. Seemann, Interface morphologies in liquid/liquid dewetting, Chem. Eng. Process. 50 , 531 (2011).
- [56] L. S. Fisher and A. A. Golovin, Nonlinear stability analysis of a two-layer thin liquid film: Dewetting and autophobic behavior., J. Colloid Interface Sci. 291 , 515 (2005).
- [57] L. S. Fisher and A. A. Golovin, Instability of a two-layer thin liquid film with surfactants: dewetting waves, J. Colloid Interface Sci. 307 , 203 (2007).
- [58] A. A. Nepomnyashchy and I. B. Simanovskii, Dynamics of ultra-thin two-layer films under the action of inclined temperature gradients, J. Fluid Mech. 631 , 165 (2009).
- [59] D. Merkt, A. Pototsky, M. Bestehorn, and U. Thiele, Long-wave theory of bounded two-layer films with a free liquid-liquid interface: Short- and long-time evolution, Phys. Fluids 17 , 064104 (2005).
- [60] G. F. Teletzke, H. T. Davis, and L. E. Scriven, Wetting hydrodynamics, Rev. Phys. Appl. (Paris) 23 , 989 (1988).
- [61] V. M. Starov and M. G. Velarde, Surface forces and wetting phenomena, J. Phys.-Condens. Matter 21 , 464121 (2009).
- [62] A. P. Hughes, U. Thiele, and A. J. Archer, Influence of the fluid structure on the binding potential: comparing liquid drop profiles from density functional theory with results from mesoscopic theory, J. Chem. Phys. 146 , 064705 (2017).
- [63] H. Yin, D. N. Sibley, U. Thiele, and A. J. Archer, Films, layers and droplets: The effect of near-wall fluid structure on spreading dynamics, Phys. Rev. E 95 , 023104 (2017).
- [64] 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¨ afer (Springer, Berlin, Heidelberg, 2006) pp. 19-49.
- [65] See Supplemental Material at link.aps.org/supplemental/10.1103/PhysRevFluids.10.024002 for additional movies and snapshots for the spreading, sliding, and coarsening simulations.
- [66] S. Engelnkemper, M. Wilczek, S. V. Gurevich, and U. Thiele, Morphological transitions of sliding drops - dynamics and bifurcations, Phys. Rev. Fluids 1 , 073901 (2016), http://arxiv.org/abs/1607.05482.
- [67] E. Infeld and G. Rowlands, Nonlinear Waves, Solitons and Chaos (University Press, Cambridge, 1990).

- [68] M. Galvagno, D. Tseluiko, H. Lopez, and U. Thiele, Continuous and discontinuous dynamic unbinding transitions in drawn film flow, Phys. Rev. Lett. 112 , 137803 (2014).
- [69] 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).
- [70] J. A. Diez, A. G. Gonz´ alez, D. A. Garfinkel, P. D. Rack, J. T. McKeown, and L. Kondic, Simultaneous decomposition and dewetting of nanoscale alloys: A comparison of experiment and theory, Langmuir 37 , 2575 (2021).
- [71] M. Areshi, D. Tseluiko, U. Thiele, B. D. Goddard, and A. J. Archer, Binding potential and wetting behavior of binary liquid mixtures on surfaces, Phys. Rev. E 109 , 024801 (2024).
- [72] K. B. Glasner and T. P. Witelski, Collision versus collapse of droplets in coarsening of dewetting thin films, Physica D 209 , 80 (2005).
- [73] F. Otto, T. Rump, and D. Slepcev, Coarsening rates for a droplet model: rigorous upper bounds, SIAM J. Math. Anal. 38 , 503 (2006).
- [74] A. A. Pahlavan, L. Cueto-Felgueroso, A. E. Hosoi, G. H. McKinley, and R. Juanes, Thin films in partial wetting: stability, dewetting and coarsening, J. Fluid Mech. 845 , 642 (2018).
- [75] C. Henkel, J. H. Snoeijer, and U. Thiele, Gradient-dynamics model for liquid drops on elastic substrates, Soft Matter 17 , 10359 (2021).
- [76] C. B. Sullivan and A. Kaszynski, PyVista: 3d plotting and mesh analysis through a streamlined interface for the visualization toolkit (VTK), Journal of Open Source Software 4 , 1450 (2019).
- [77] J. Diekmann, Macroscopic and Mesoscopic Modelling of Compound Drops , Master's thesis, University of M¨nster, M¨ unster (2022), u master thesis.
- [78] L. Kondic, A. G. Gonz´ alez, J. A. Diez, J. D. Fowlkes, and P. Rack, Liquid-state dewetting of pulsed-laser-heated nanoscale metal films and other geometries, Ann. Rev. Fluid Mech. 52 , 235 (2020).
- [79] U. Thiele, L. Brusch, M. Bestehorn, and M. B¨ ar, Modelling thin-film dewetting on structured substrates and templates: Bifurcation analysis and numerical simulations, Eur. Phys. J. E 11 , 255 (2003).
- [80] N. Savva and S. Kalliadasis, Influence of gravity on the spreading of two-dimensional droplets over topographical substrates, J. Eng. Math. 73 , 3 (2012).
- [81] U. Thiele and E. Knobloch, Driven drops on heterogeneous substrates: Onset of sliding motion, Phys. Rev. Lett. 97 , 204501 (2006).
- [82] Z. Lin, T. Kerle, T. P. Russell, E. Sch¨ affer, and U. Steiner, Structure formation at the interface of liquid-liquid bilayer in electric field, Macromolecules 35 , 3971 (2002).
- [83] G. Amarandei, P. Beltrame, I. Clancy, C. O'Dwyer, A. Arshak, U. Steiner, D. Corcoran, and U. Thiele, Pattern formation induced by an electric field in a polymer-air-polymer thin film system, Soft Matter 8 , 6333 (2012).
- [84] 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).
- [85] 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).
- [86] A. Oron and P. Rosenau, Formation of patterns induced by thermocapillarity and gravity, J. Physique II France 2 , 131 (1992).
- [87] U. Thiele and E. Knobloch, Thin liquid films on a slightly inclined heated plate, Physica D 190 , 213 (2004).
Jan Diekmann Uwe Thiele
Mesoscopic hydrodynamic model for spreading, sliding and coarsening compound drops
We revisit the mesoscopic hydrodynamic description of the dynamics of sessile
partially wetting compound drops, i.e., of drops that consist of two immiscible
nonvolatile partially wetting liquids and are situated on a smooth rigid solid
substrate. We briefly discuss and complete existing dynamic models employing a
gradient dynamics approach. Thereby, the underlying energy features capillarity
and wettability contributions for all relevant interfaces in full-curvature
formulation. Establishing transparent consistency relations between macroscopic
and mesoscopic parameters, we obtain mesoscopic Neumann and Young laws that are
also fully consistent with the macroscopic ones. In particular, we discuss the
minimal requirements for the wetting energy that ensure the full spectrum of
macroscopic parameters {for partially wetting cases} is addressed by the
mesoscopic model. Subsequently, we distinguish long-wave and full-curvature
variants of the dynamical model based on properties of the energy, and employ
the latter to illustrate the usage of the mesoscopic model. As examples, we
chose the spreading of individual compound drops on one-dimensional horizontal
substrates, sliding compound drops on one-dimensional inclined substrates, and
the coarsening of drop ensembles on one- and two-dimensional horizontal
substrates. In each case, the discussion emphasizes occurring qualitative
changes in the drop configurations.