Comments on the compositional boundary condition for diffuse interface model of a contact line
Abstract
The diffuse interface model of Cahn-Hilliard-van der Waals is often used to study various aspects of multi-phase flows such as droplets coalescence and contact line dynamics. The original model of Cahn-Hilliard-van der Waals uses an approximation which neglects some surface contributions to the free energy of the system and is justified by the large distance between the region of observation and the outer surface of the system. This is not a priori accurate when considering the dynamics of the contact line on a solid surface. A modification of the so called natural compositional boundary conditions obtained by minimization of the surface contribution to the free energy of the system is proposed for the vapor-liquid interface of a single component and for a quasi-incompressible binary mixture of partially miscible liquids. The results are obtained for a stagnant fluid in thermodynamic equilibrium but should be valid in a flow if local thermodynamic equilibrium can be assumed.
pacs:
I Introduction
A sharp interface continuum model cannot describe the movement of a contact line on a solid surface when the so called no-slip boundary condition is imposed at the wall. Various modifications of the sharp interface model have been described in review papers ((Bonn et al., 2009), (Sui et al., 2014), (Snoeijer and Andreotti, 2013)). Seppecher Seppecher (1996) and Jacqmin Jacqmin (2000) demonstrated that a diffuse interface model based on the modification of the theory of van der Waals van der Waals (1979) by Cahn and Hilliard Cahn and Hilliard (1958) does predict a motion of a contact line when a no-slip boundary condition is applied at the wall. The diffuse interface model uses a local thermodynamic equilibrium assumption in which the specific free energy of the fluid depends not only on the local thermodynamic properties of the fluid but also on the spatial gradients in these quantities. In the Cahn-Hiliard-van der Waals approximation the gradient length is assumed to be large compared to the molecular attraction length scale . When considering isotropic fluids the first gradient term vanishes in a Taylor expansion of the specific free energy around its value for a uniform fluid. The first non-vanishing non-local terms are a term proportional to the Laplacian of the fluid properties and a term proportional to the square of the gradient. When considering the total free energy of the system, the Laplacian term is split by partial integration into a bulk term proportional to the square of the gradient and a surface term. This surface term is neglected ((van der Waals, 1979),(Cahn and Hilliard, 1958)) using the argument that the region of interest is far removed from walls. Obviously this is questionable when considering contact line dynamics. The application of the diffuse interface theory to wetting in the seminal paper of Cahn Cahn (1977) uses this approximation as do the more recent studies ((Seppecher, 1996), (Jacqmin, 2000), (Qian et al., 2006), (Khatavkar et al., 2007), (Yue et al., 2010), Yue and Feng (2011), (Carlson et al., 2011), (Lee and Kim, 2011), (Magaletti et al., 2013), (Sibley et al., 2013)). To the authors knowledge this problem has not yet been discussed in the literature.
An advantage of the gradient theory is that in principle one could use a single set of equations through a two phase system. However as explained by Yue et al. Yue et al. (2010) in practice the spatial discretisation cannot be sufficient to resolve typical length scales in interfaces. Hence in principle only results independent from the exact value of the gradient length are physically relevant. While Yue et al. Yue et al. (2010) propose a method to compensate for errors due to too large interface thickness, the multi-scale approach of Seppecher Seppecher (1996) is promising as well. He matches an inner diffuse interface model to a far-field sharp interface model. Such a multi-scale model could be further extended to include local molecular phenomena near the contact line. An example of this is including a disjoining potential term making the interaction energy between the fluid and the solid dependent on the local thickness of
the liquid film. Also one could introduce stochastic source terms in the equations as proposed by Snoeijer and Andreotti Snoeijer and Andreotti (2013). We however limit our discussion to the most commonly used model proposed by Cahn Cahn (1977) in which the surface free energy is a function of the local bulk fluid properties adjacent to the wall. Hence the model remains essentially a continuum model. We will show that in the case of Cahn’s Cahn (1977) model for a semi-infinite two phase fluid in contact with a wall the extra surface free energy term cannot be neglected.
As proposed by Anderson et al. Anderson et al. (1998) we limit our discussion to a stagnant fluid at thermodynamic equilibrium. The results obtained can then be used to derive equations of motion if one assumes local thermodynamic equilibrium Anderson et al. (1998); Mauri (2013).
Firstly, the Cahn-Hilliard-van der Waals theory applied to the interface between a vapor and a liquid of a single component system will be discussed with focus on modified wall compositional boundary conditions. Secondly, results for the commonly studied quasi-incompressible binary mixture of partially miscible fluids (Lowengrub and Truskinovsky, 1998) are summarized. Finally the incompressible binary regular solution model of Cahn Cahn (1977) is discussed.
Ii vapor liquid interface
As stated before we limit our discussion to a stagnant fluid in thermodynamic equilibrium. The condition for thermodynamic equilibrium for a single component isolated system of volume enclosed by a surface is found by van der Waals (van der Waals, 1979) by seeking for a maximum of entropy at constant mass and internal energy . Ignoring surface effects, this corresponds to a maximum of the entropy:
(1) |
where is the specific entropy and is the fluid density, with the constrains:
(2) |
and
(3) |
where the specific internal energy. The constrained maximum of the entropy can be written using the Lagrange multipliers method as:
(4) |
where is the specific Helmholtz free energy. By definition the Lagrangian multipliers and are constants. The Lagrangian multiplier is the inverse of the temperature, which in equilibrium is uniform within . The Lagrange multiplier is the chemical potential divided by the temperature, which in equilibrium is uniform over the system.
We consider conditions such that two uniform phases phases coexist: a liquid phase of density and its vapor of density . The densities in the range (between the spinodal points) are for uniform systems either unstable or meta-stable. The key idea of van der Waals (van der Waals, 1979) is that in the interface of finite thickness between the liquid and the vapor all the densities in the range are present and stable. Thus the non-uniformity of the fluid in the interface stabilizes the fluid, allowing static equilibrium. Hence, for the closed system considered, one should minimise the total free energy of the system taking into account the fact that the specific free energy of the fluid depends not only on the local density but also on the gradients of the density: . The next step is to assume that the contribution of the gradients is small, in other words that the molecular length scale is small compared to the gradient length (the interface thickness) i.e. . For an isotropic fluid the leading order terms in a Taylor expansion of around a uniform state of density are given by:
(5) |
The index indicates that the coefficients are functions of only. As proposed by Cahn Cahn (1977), the interaction of the fluid with the wall surface is assumed to be described by a free energy per unit surface , which is a function of the density of the fluid at the wall (in a material element, large compared to , adjacent to the wall). The molecular interaction between the fluid and the wall is assumed to be limited to a length scale (monolayer), short compared to other length scales in the problem (Cahn, 1977). This leads to:
(6) | |||||
where it the outer normal on . In this expression we have carried out a partial integration to express the term with the Laplacian in the bulk of the fluid in terms proportional to and a surface contribution proportional to . We neglected the mass of the fluid monolayer on the wall surface compared to the mass of the fluid in .
Using the fact that in a stagnant fluid , by carrying a second partial integration and considering that the perturbations are arbitrary, we obtain:
(7) |
with defined as is usual in the literature:
(8) |
In the literature one often assumes that is a constant but in general it is a function of . After the second partial integration and application of the variation the surface integral becomes:
(9) |
Assuming that is a function of and consequently using one obtains the compositional natural boundary condition in the form of a differental equation for :
(10) |
We should now provide an initial condition for the integration of this equation. Note that this integration is carried out in the thermodynamic space and can be carried out independently of any flow simulations. This integration determines as a function of . For a wetting wall () we expect that at very high densities the wall will be saturated. Let be the maximum thermodynamically allowable liquid phase density at the temperature considered (limit of very high pressures). In principle is outside the range of that will be found in the actual flow. We expect that the derivative of with respect to vanishes for (i.e. ) and that . This ad hoc assumption should in principle be confirimed on the basis of physical models for the interaction of the fluid with the wall. Such an assumption could, however, also be imposed for convenience (without justification).
Note that Mauri (Mauri, 2013) obtained a similar result but stated that the perturbation can be assumed to vanish. This results into an explicit equation for as a function of , similar to the compositional boundary condition used by Seppecher Seppecher (1996) and Jacqmin Jacqmin (2000). Consequently the assumption of Mauri Mauri (2013) that is independent from is contradicted by the his conclusion that is a function of .
Iii Quasi-incompressible mixture of partially miscible liquids
We now consider a quasi-incompressible binary mixture of liquids (Lowengrub and Truskinovsky, 1998). The mixture has two equilibrium concentrations and ( is the mass fraction of one of the components) corresponding to liquid-liquid equilibrium of two uniform phases. The definition of a quasi-incompressible mixture is that the density of the mixture is a function of the mass fraction , but is independent of the pressure. The wall free energy density is function of as well. While the bulk specific free energy is a function of and gradients of viz. . The free energy of the system should now be minimized with the constraint of total mass conservation ( constant) and mass conservation for one of the two components ( constant). We furthermore assume a Taylor expansion similar to equation (5) and follow the same procedure as for a single component. Results concerning the chemical potential and reversible stress tensor are identical to the results given by Lowengrub and Truskinovsky (1998). For the natural compositional boundary condition we now find the differential equation for the normal gradient of composition at the surface :
(11) | |||||
An initial condition for integration of this differential equation could be obtained by considering the limit of a fully saturated wall (with one of the components). This corresponds to the mass fraction for . For the limit it seems rather logical to assume that because the mass fraction cannot become larger than one (). As in the single component case the integration is carried out in the thermodynamic variable space. The actual range in which varies in a physical system will in general not include . This model could be applied for convenience. An example of the integration in thermodynamic variable space is provided in the next section.
In the case of a compressible binary mixture we have . Following a similar procedure one obtains a set of partial differential equations for the gradients at the wall and as functions of and .
Iv The influence of the extra surface free energy term
Here we derive a model for the excess free energy of a semi infinite two phase binary mixture of partially miscible fluids in contact with a planar surface that takes into account the surface term neglected by Cahn Cahn (1977). For convenience in comparing the results obtained here to the original paper Cahn (1977) we consider the free energy per unit volume to b e a function of the molar fraction of one of the components. We do this for the particular case of an incompressible regular solution model as described in the paper of Cahn and Hilliard Cahn and Hilliard (1958). This will allow us to compare the term neglected with other terms in the model giving us some insight into whether or not this term is justifiably neglected by Cahn Cahn (1977).
We define the excess free energy per unit area as the free energy per unit area of a semi infinite two phase binary mixture of partially miscible fluids in contact with a planar surface minus that of a semi infinite uniform reference system . Hence this is the energy needed to create the fluid in contact with the wall out of a uniform reference fluid of concentration , the molar fraction far away from the wall. The reference state is one for which the fluid has no interaction with the wall. Thus we have
(12) |
The wall free energy per unit surface is a function of the molar fraction of the fluid adjacent to the wall. The free energy per unit volume can be expanded in a Taylor series around the uniform state which for an isotropic fluid yields Cahn and Hilliard (1958)
(13) |
Where and , the subscript specifies functions of (viz. not of and ). Substituting equation (13) into equation (12) yields
(14) |
Performing partial integration on the second term in the integrant and using and yields
(15) |
Where . For the particular case considered i.e. the regular solution model of Cahn and Hilliard Cahn and Hilliard (1958) we have where is a positive constant and . The free energy of the reference state is
(16) |
where . Thus defining the excess free energy can be written as:
(17) |
with . Taking the variation of we have
(18) | |||||
In equilibrium we have where surface term of the variation and the bulk term should vanish independently Courant and Hilbert (1963). Setting and we find the equilibrium condition for the bulk
(19) |
We can integrate this ordinary differential equation using the boundary condition and in the limit to find:
(20) |
Which can be rewritten to find
(21) |
where we have chosen the negative root because we take in this model (the wall is partially wetting) n.b. .
The surface tension is the minimum in excess free energy of the surface per unit area Cahn (1977) viz.
(22) |
Which except for the extra surface term corresponds to equation (10a) of Cahn Cahn (1977). Using we can rewrite this equation as follows
(23) |
This equation corresponds to equation (10b) of Cahn Cahn (1977) with an added surface free energy term . Setting and in equation (18) yields
(24) |
the so called natural boundary condition at the wall.
We assume i.e. the surface is partially wetting. For simplicity we assume . The regular solution model of Cahn-Hilliard Cahn and Hilliard (1958) results in and . Substituting this into equation (24) and defining , we have
(25) |
Which is a first order ordinary differential equation. Solving this equation and assuming for (as proposed for the quasi-incompressible case), yields
(26) |
We will use this to investigate the importance of the surface term in equation (23). That is we have
(27) |
where for convenience we have used the assumption that is a constant. Substituting equation (26) into equation (27) yields
(28) |
The second term on the right hand side of equation (28) is a extra surface free energy term omitted by Cahn Cahn (1977).
We can now go compare the term neglected by Cahn Cahn (1977) to the third term on the right hand side of equation (28), which is a bulk term also present in Cahn’s expression Cahn (1977), taking the fraction of the extra surface term with respect to the bulk term yields
(29) |
Which in general will not be small for and . We conclude that one cannot make the a priori assumption that the results obtained by Cahn and Hilliard Cahn and Hilliard (1958) for cases where boundary effects are neglected can be extended to problems that involve surface effects e.g. critical point wetting Cahn (1977).
V Discussion and conclusion
The square gradient model always implies a partial integration of the Laplacian term in the van der Waals expansion of the free energy taking non-local effects into account. This generates a surface term that is neglected in the literature without justification. For binary regular solutions Cahn and Hilliard (1958) and Verschueren (1999) show that the Laplacian term term is the leading order term in the van der Waals expansion. Hence one cannot use the argument that the term in (or ) is negligible compared to the square gradient term in the bulk. Therefore one cannot avoid the partial integration of the Laplacian term as a step in the diffuse interface theory. When considering contact line motion or wetting phenomena one should therefore quantify the surface contribution to the wall free energy resulting from this partial integration of the Laplacian term.
We have obtained a compositional natural compositional boundary condition at the wall for the Cahn-Hilliard-van der Waals diffuse interface model in the case of a single component vapor-liquid system and for a two phase system in a quasi-incompressible binary mixture of partially miscible fluids. Specifically we did not neglect the surface term proportional to in the free energy as done by Cahn and Hilliard ((Cahn and Hilliard, 1958), (Cahn, 1977)). This compositional boundary condition yields a differential equation for the wall normal gradient (or ) that should be integrated in the thermodynamic space. Determining a physically relevant initial condition for the integration of this differential equation (10) or (11) will depend on the thermodynamic model used for the bulk of the fluid and for the fluid-wall interaction. We proposed a possible initial condition.
At the present time the use of a near critical equation of state (Ginsburg-Landau (Jacqmin, 2000)) has become common practice. In order to avoid critical point wetting (Cahn, 1977) a particular expression for proposed by Jacqmin (2000) is used, imposing for and . These assumptions considerably simplify the numerical solution of the problem but make it difficult to use the theory for quantitative prediction of the behavior of real physical systems. Improving the natural compositional boundary condition for such simplified models has probably little added value. However for more realistic physical models this could be significant and deserves to be investigated.
References
- D. Bonn, J. Eggers, J. Indekeu, J. Meunier, and E. Rolley, Rev. Modern Physics, 81, 739 (2009).
- Y. Sui, H. Ding, and P. Spelt, Ann. Rev. Fluid Mech., 46, 97 (2014).
- J. Snoeijer and B. Andreotti, Ann. Rev. Fluid Mech., 45, 269 (2013).
- P. Seppecher, Int. J. Engng Sci., 34, 977 (1996).
- D. Jacqmin, J. Fluid Mech., 402, 57 (2000).
- J. D. van der Waals, Journal of Statistical Physics, 20, 200 (1979).
- J. Cahn and J. Hilliard, J. Chemical Physics, 28, 258 (1958).
- J. Cahn, J. Chemical Physics, 66, 3667 (1977).
- T. Qian, X. Wang, and P. Sheng, J. Fluid Mech., 564, 333 (2006).
- V. Khatavkar, P. Anderson, and H. Meijer, J. Fluid Mech., 572, 367 (2007).
- N. Yue, C. Zhou, and J. Feng, J. Fluid Mech., 645, 279 (2010).
- P. Yue and J. J. Feng, Physics of Fluids, 23, 012106 (2011).
- A. Carlson, M. Doquand, and G. Amberg, J. Fluid Mech., 682, 213 (2011).
- H. Lee and J. Kim, Computers and Fluids, 44, 178 (2011).
- F. Magaletti, F. Picano, M. Chinappi, Marino, and C. M. Casciola, J. Fluid Mech., 714, 95 (2013).
- D. Sibley, A. Nold, N. Savva, and S. Kalliadasis, Phys. Fluids, 25, 082105 (2013).
- D. Anderson, G. McFadden, and A. Wheeler, Ann. Rev. Fluid Mech., 30, 139 (1998).
- R. Mauri, Non-Equilibrium Thermodynamics in Multiphase Flows, Soft and Biological Matter No. ISBN 978-94-007-5461-4 (Springer, Dordrecht, The Netherlands, 2013).
- J. Lowengrub and L. Truskinovsky, Proc. R. Soc. Lond. A, 454, 2617 (1998).
- R. Courant and D. Hilbert, Methods of Mathematical Physics, 4th ed. (Inter science publishers, inc., 1963).
- M. Verschueren, A diffuse-interface model for structure development in flow, Ph.D. thesis, Technische Universiteit Eindhoven, The Netherlands (1999).