Capillary surface

Last updated

In fluid mechanics and mathematics, a capillary surface is a surface that represents the interface between two different fluids. As a consequence of being a surface, a capillary surface has no thickness in slight contrast with most real fluid interfaces.

Contents

Capillary surfaces are of interest in mathematics because the problems involved are very nonlinear and have interesting properties, such as discontinuous dependence on boundary data at isolated points. [1] In particular, static capillary surfaces with gravity absent have constant mean curvature, so that a minimal surface is a special case of static capillary surface.

They are also of practical interest for fluid management in space (or other environments free of body forces), where both flow and static configuration are often dominated by capillary effects.

The stress balance equation

The defining equation for a capillary surface is called the stress balance equation, [2] which can be derived by considering the forces and stresses acting on a small volume that is partly bounded by a capillary surface. For a fluid meeting another fluid (the "other" fluid notated with bars) at a surface , the equation reads

where is the unit normal pointing toward the "other" fluid (the one whose quantities are notated with bars), is the stress tensor (note that on the left is a tensor-vector product), is the surface tension associated with the interface, and is the surface gradient. Note that the quantity is twice the mean curvature of the surface.

In fluid mechanics, this equation serves as a boundary condition for interfacial flows, typically complementing the Navier–Stokes equations. It describes the discontinuity in stress that is balanced by forces at the surface. As a boundary condition, it is somewhat unusual in that it introduces a new variable: the surface that defines the interface. It's not too surprising then that the stress balance equation normally mandates its own boundary conditions.

For best use, this vector equation is normally turned into 3 scalar equations via dot product with the unit normal and two selected unit tangents:

Note that the products lacking dots are tensor products of tensors with vectors (resulting in vectors similar to a matrix-vector product), those with dots are dot products. The first equation is called the normal stress equation, or the normal stress boundary condition. The second two equations are called tangential stress equations.

The stress tensor

The stress tensor is related to velocity and pressure. Its actual form will depend on the specific fluid being dealt with, for the common case of incompressible Newtonian flow the stress tensor is given by

where is the pressure in the fluid, is the velocity, and is the viscosity.

Static interfaces

In the absence of motion, the stress tensors yield only hydrostatic pressure so that , regardless of fluid type or compressibility. Considering the normal and tangential equations,

The first equation establishes that curvature forces are balanced by pressure forces. The second equation implies that a static interface cannot exist in the presence of nonzero surface tension gradient.

If gravity is the only body force present, the Navier–Stokes equations simplify significantly:

If coordinates are chosen so that gravity is nonzero only in the direction, this equation degrades to a particularly simple form:

where is an integration constant that represents some reference pressure at . Substituting this into the normal stress equation yields what is known as the Young-Laplace equation:

where is the (constant) pressure difference across the interface, and is the difference in density. Note that, since this equation defines a surface, is the coordinate of the capillary surface. This nonlinear partial differential equation when supplied with the right boundary conditions will define the static interface.

The pressure difference above is a constant, but its value will change if the coordinate is shifted. The linear solution to pressure implies that, unless the gravity term is absent, it is always possible to define the coordinate so that . Nondimensionalized, the Young-Laplace equation is usually studied in the form [1]

where (if gravity is in the negative direction) is positive if the denser fluid is "inside" the interface, negative if it is "outside", and zero if there is no gravity or if there is no difference in density between the fluids.

This nonlinear equation has some rich properties, especially in terms of existence of unique solutions. For example, the nonexistence of solution to some boundary value problem implies that, physically, the problem can't be static. If a solution does exist, normally it'll exist for very specific values of , which is representative of the pressure jump across the interface. This is interesting because there isn't another physical equation to determine the pressure difference. In a capillary tube, for example, implementing the contact angle boundary condition will yield a unique solution for exactly one value of . Solutions often aren't unique, this implies that there are multiple static interfaces possible; while they may all solve the same boundary value problem, the minimization of energy will normally favor one. Different solutions are called configurations of the interface.

Energy consideration

A deep property of capillary surfaces is the surface energy that is imparted by surface tension:

where is the area of the surface being considered, and the total energy is the summation of all energies. Note that every interface imparts energy. For example, if there are two different fluids (say liquid and gas) inside a solid container with gravity and other energy potentials absent, the energy of the system is

where the subscripts , , and respectively indicate the liquid–gas, solid–gas, and solid–liquid interfaces. Note that inclusion of gravity would require consideration of the volume enclosed by the capillary surface and the solid walls.

Illustration of distributed forces at the contact line, with the contact line perpendicular to the image. The vertical part of
g
L
G
{\displaystyle \gamma _{LG}}
is balanced by a slight deformation of the solid (not shown and inconsequential to this context). ForcesAtContactLine.png
Illustration of distributed forces at the contact line, with the contact line perpendicular to the image. The vertical part of is balanced by a slight deformation of the solid (not shown and inconsequential to this context).

Typically the surface tension values between the solid–gas and solid–liquid interfaces are not known. This does not pose a problem; since only changes in energy are of primary interest. If the net solid area is a constant, and the contact angle is known, it may be shown that (again, for two different fluids in a solid container)

so that

where is the contact angle and the capital delta indicates the change from one configuration to another. To obtain this result, it's necessary to sum (distributed) forces at the contact line (where solid, gas, and liquid meet) in a direction tangent to the solid interface and perpendicular to the contact line:

where the sum is zero because of the static state. When solutions to the Young-Laplace equation aren't unique, the most physically favorable solution is the one of minimum energy, though experiments (especially low gravity) show that metastable surfaces can be surprisingly persistent, and that the most stable configuration can become metastable through mechanical jarring without too much difficulty. On the other hand, a metastable surface can sometimes spontaneously achieve lower energy without any input (seemingly at least) given enough time.

Boundary conditions

Boundary conditions for stress balance describe the capillary surface at the contact line: the line where a solid meets the capillary interface; also, volume constraints can serve as boundary conditions (a suspended drop, for example, has no contact line but clearly must admit a unique solution).

For static surfaces, the most common contact line boundary condition is the implementation of the contact angle, which specifies the angle that one of the fluids meets the solid wall. The contact angle condition on the surface is normally written as:

where is the contact angle. This condition is imposed on the boundary (or boundaries) of the surface. is the unit outward normal to the solid surface, and is a unit normal to . Choice of depends on which fluid the contact angle is specified for.

For dynamic interfaces, the boundary condition showed above works well if the contact line velocity is low. If the velocity is high, the contact angle will change ("dynamic contact angle"), and as of 2007 the mechanics of the moving contact line (or even the validity of the contact angle as a parameter) is not known and an area of active research. [3]

See also

Related Research Articles

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

In vector calculus, the divergence theorem, also known as Gauss's theorem or Ostrogradsky's theorem, is a theorem relating the flux of a vector field through a closed surface to the divergence of the field in the volume enclosed.

In continuum mechanics, the infinitesimal strain theory is a mathematical approach to the description of the deformation of a solid body in which the displacements of the material particles are assumed to be much smaller than any relevant dimension of the body; so that its geometry and the constitutive properties of the material at each point of space can be assumed to be unchanged by the deformation.

<span class="mw-page-title-main">Four-vector</span> 4-dimensional vector in relativity

In special relativity, a four-vector is an object with four components, which transform in a specific way under Lorentz transformations. Specifically, a four-vector is an element of a four-dimensional vector space considered as a representation space of the standard representation of the Lorentz group, the representation. It differs from a Euclidean vector in how its magnitude is determined. The transformations that preserve this magnitude are the Lorentz transformations, which include spatial rotations and boosts.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

<span class="mw-page-title-main">Green's function</span> Impulse response of an inhomogeneous linear differential operator

In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.

A Newtonian fluid is a fluid in which the viscous stresses arising from its flow are at every point linearly correlated to the local strain rate — the rate of change of its deformation over time. Stresses are proportional to the rate of change of the fluid's velocity vector.

<span class="mw-page-title-main">Euler equations (fluid dynamics)</span> Set of quasilinear hyperbolic equations governing adiabatic and inviscid flow

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

<span class="mw-page-title-main">Scalar potential</span> When potential energy difference depends only on displacement

In mathematical physics, scalar potential, simply stated, describes the situation where the difference in the potential energies of an object in two different positions depends only on the positions, not upon the path taken by the object in traveling from one position to the other. It is a scalar field in three-space: a directionless value (scalar) that depends only on its location. A familiar example is potential energy due to gravity.

In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

In physics, the Young–Laplace equation is an algebraic equation that describes the capillary pressure difference sustained across the interface between two static fluids, such as water and air, due to the phenomenon of surface tension or wall tension, although use of the latter is only applicable if assuming that the wall is very thin. The Young–Laplace equation relates the pressure difference to the shape of the surface or wall and it is fundamentally important in the study of static capillary surfaces. It is a statement of normal stress balance for static fluids meeting at an interface, where the interface is treated as a surface :

The derivation of the Navier–Stokes equations as well as their application and formulation for different families of fluids, is an important exercise in fluid dynamics with applications in mechanical engineering, physics, chemistry, heat transfer, and electrical engineering. A proof explaining the properties and bounds of the equations, such as Navier–Stokes existence and smoothness, is one of the important unsolved problems in mathematics.

In fluid dynamics, aerodynamic potential flow codes or panel codes are used to determine the fluid velocity, and subsequently the pressure distribution, on an object. This may be a simple two-dimensional object, such as a circle or wing, or it may be a three-dimensional vehicle.

In fluid dynamics, Airy wave theory gives a linearised description of the propagation of gravity waves on the surface of a homogeneous fluid layer. The theory assumes that the fluid layer has a uniform mean depth, and that the fluid flow is inviscid, incompressible and irrotational. This theory was first published, in correct form, by George Biddell Airy in the 19th century.

<span class="mw-page-title-main">Stokes' theorem</span> Theorem in vector calculus

Stokes' theorem, also known as the Kelvin–Stokes theorem after Lord Kelvin and George Stokes, the fundamental theorem for curls or simply the curl theorem, is a theorem in vector calculus on . Given a vector field, the theorem relates the integral of the curl of the vector field over some surface, to the line integral of the vector field around the boundary of the surface. The classical theorem of Stokes can be stated in one sentence: The line integral of a vector field over a loop is equal to the surface integral of its curl over the enclosed surface. It is illustrated in the figure, where the direction of positive circulation of the bounding contour ∂Σ, and the direction n of positive flux through the surface Σ, are related by a right-hand-rule. For the right hand the fingers circulate along ∂Σ and the thumb is directed along n.

Curvilinear coordinates can be formulated in tensor calculus, with important applications in physics and engineering, particularly for describing transportation of physical quantities and deformation of matter in fluid mechanics and continuum mechanics.

In mathematical physics, the Gordon decomposition of the Dirac current is a splitting of the charge or particle-number current into a part that arises from the motion of the center of mass of the particles and a part that arises from gradients of the spin density. It makes explicit use of the Dirac equation and so it applies only to "on-shell" solutions of the Dirac equation.

Morris Muskat et al. developed the governing equations for multiphase flow in porous media as a generalisation of Darcy's equation for water flow in porous media. The porous media are usually sedimentary rocks such as clastic rocks or carbonate rocks.

In fluid mechanics, the thin-film equation is a partial differential equation that approximately predicts the time evolution of the thickness h of a liquid film that lies on a surface. The equation is derived via lubrication theory which is based on the assumption that the length-scales in the surface directions are significantly larger than in the direction normal to the surface. In the non-dimensional form of the Navier-Stokes equation the requirement is that terms of order ε2 and ε2Re are negligible, where ε ≪ 1 is the aspect ratio and Re is the Reynolds number. This significantly simplifies the governing equations. However, lubrication theory, as the name suggests, is typically derived for flow between two solid surfaces, hence the liquid forms a lubricating layer. The thin-film equation holds when there is a single free surface. With two free surfaces, the flow must be treated as a viscous sheet.

The streamline upwind Petrov–Galerkin pressure-stabilizing Petrov–Galerkin formulation for incompressible Navier–Stokes equations can be used for finite element computations of high Reynolds number incompressible flow using equal order of finite element space by introducing additional stabilization terms in the Navier–Stokes Galerkin formulation.

References

  1. 1 2 Robert Finn (1999). "Capillary Surface Interfaces" (PDF). American Mathematical Society.
  2. Landau, Lev D.; Lifshitz, Evgeny M. (1987). Fluid Mechanics. Vol. 6 (2nd ed.). Butterworth-Heinemann. ISBN   978-0-08-033933-7.
  3. E. B. Dussan V; Enrique Ramé; Stephen Garoff (September 1991). "On identifying the appropriate boundary conditions at a moving contact line: an experimental investigation". Journal of Fluid Mechanics. 230. CJO: 97–116. Bibcode:1991JFM...230...97D. doi:10.1017/S0022112091000721.