User:Hh73wiki/Dynamic Fluid Film Equations Update



Fluid films, such as soap films, are commonly encountered in everyday experience. A soap film can be formed by dipping a closed contour wire into a soapy solution as in the figure on the right. Alternatively, a catenoid can be formed by dipping two rings in the soapy solution and subsequently separating them while maintaining the coaxial configuration.

Stationary fluid films form surfaces of minimal surface area, leading to the Plateau problem.

On the other hand, fluid films display fascinating and rich dynamic properties. They can undergo enormous deformations away from the equilibrium configuration. Furthermore, they display several orders of magnitude variations in thickness from nanometers to millimeters. Thus, a fluid film can simultaneously display nanoscale and macroscale phenomena.

In the study of the dynamics of free fluid films, such as soap films, it is common to model the film as two dimensional manifolds. Then the variable thickness of the film is captured by the two dimensional density $$\rho $$.

The dynamics of fluid films can be described by the following system of exact nonlinear Hamiltonian equations which, in that respect, are a complete analogue of Euler's inviscid equations of fluid dynamics. In fact, these equations reduce to Euler's dynamic equations for flows in stationary Euclidean spaces.

The foregoing relies on the formalism of tensors, including the summation convention and the raising and lowering of tensor indices.

The full dynamic system
Consider a thin fluid film $$S$$ that spans a stationary closed contour boundary. Let $$C$$ be the normal component of the velocity field and $$V^{\alpha }$$ be the contravariant components of the tangential velocity projection. Let $$\nabla _{\alpha }$$ be the covariant surface derivative, $$B_{\alpha \beta }$$ be the covariant curvature tensor, $$B^{\alpha }_{\beta }$$ be the mixed curvature tensor and $$B^{\alpha }_{\alpha }$$ be its trace, that is mean curvature. Furthermore, let the internal energy density per unit mass function be $$e\left( \rho \right)$$ so that the total potential energy $$E$$ is given by


 * $$E=\int_{S}\rho e\left( \rho \right) \, dS.$$

This choice of $$e\left( \rho \right)$$ :


 * $$e\left( \rho \right) =\frac{\sigma }{\rho }$$

where $$\sigma$$ is the surface energy density results in Laplace's classical model for surface tension:


 * $$E = \sigma A \, $$

where A is the total area of the soap film.

The governing system reads


 * $$\begin{align}

\frac{\delta \rho }{\delta t}+\nabla _{\alpha }\left( \rho V^{\alpha } \right) &= \rho CB^{\alpha }_{\alpha } \\ \\  \rho \left( \frac{\delta C}{\delta t} + 2V^\alpha \nabla_\alpha C+B_{\alpha \beta }V^\alpha V^\beta \right) &= -\rho^2 e_\rho B^\alpha_\alpha \\ \\  \rho \left( \frac{\delta V^\alpha}{\delta t} + V^\beta \nabla_\beta V^\alpha - C\nabla^\alpha C - 2CV^\beta B^\alpha_\beta \right) &= -\nabla^\alpha \left( \rho^2 e_\rho \right) \end{align}$$

where the $${\delta}/{\delta} t$$-derivative is the central operator, originally due to Jacques Hadamard, in the The Calculus of Moving Surfaces. Note that, in compressible models, the combination $$\rho^2 e_\rho $$ is commonly identified with pressure $$p$$. The governing system above was originally formulated in reference 1.

For the Laplace choice of surface tension $$\left( e\left( \rho \right)=\sigma /\rho  \right)$$ the system becomes:


 * $$\begin{align}

\frac{\delta \rho }{\delta t} + \nabla_\alpha \left( \rho V^{\alpha } \right) &= \rho CB^\alpha_\alpha \\ \\  \rho \left( \frac{\delta C}{\delta t}+2V^\alpha \nabla_\alpha C+B_{\alpha \beta }V^\alpha V^\beta  \right) &= \sigma B^\alpha_\alpha \\ \\  \frac{\delta V^\alpha }{\delta t} + V^\beta \nabla_\beta V^\alpha - C\nabla ^\alpha C - 2V^\beta B^\alpha _\beta & = 0 \end{align}$$

Note that on flat ($$B_{\alpha \beta }=0$$) stationary ($$C=0$$) manifolds, the system becomes


 * $$\begin{align}

\frac{\partial \rho }{\partial t} + \nabla_\alpha \left( \rho V^{\alpha } \right) &= 0 \\ & \\   \rho \left( \frac{\partial V^\alpha}{\partial t} + V^\beta \nabla _\beta V^\alpha \right) &= -\nabla ^\alpha \left( \rho^2 e_\rho \right) \end{align}$$

which is precisely classical Euler's equations of fluid dynamics.

A simplified system
If one disregards the tangential components of the velocity field, as frequently done in the study of thin fluid film, one arrives at the following simplified system with only two unknowns: the two dimensional density $$\rho$$ and the normal velocity $$C$$:


 * $$\begin{align}

{\frac{\delta\rho }{\delta t}} &= \rho CB_\alpha^\alpha \\ &\\ \rho\frac{\delta C}{\delta t} &=\sigma B_\alpha^\alpha \\ \end{align}$$

Reformulated Fluid film equations with $$\textstyle \dot{\nabla}$$
The fluid film equations, first introduced in \cite{FluidFilmsJGSP}, are formulated in terms of the $$\textstyle \delta/\delta t$$-derivative. The last equation in the system, governs the evolution of the contravariant components $$\textstyle V^{\alpha}$$ of the tangential velocity field and features the term $$\textstyle \delta V^{\alpha }/\delta t$$ on the left hand side. Therefore, the corresponding equation for the evolution of the covariant components $$\textstyle V_{\alpha}$$ is different. The new operator $$\textstyle \dot{\nabla}$$ eliminates this incovenience and yields a system that is valid for covariant and contravariant components. In the new form, the system reads:


 * $$\begin{align}

\dot{\nabla}\rho+\nabla_{\alpha}\left( \rho V^{\alpha}\right) &= \rho CB^{\alpha }_{\alpha } \\ \\  \rho \left( \frac{\delta C}{\delta t} + 2V^\alpha \nabla_\alpha C+B_{\alpha \beta }V^\alpha V^\beta \right) &= -\rho^2 e_\rho B^\alpha_\alpha \\ \\  \rho\left(  \dot{\nabla}V_{\alpha}+V_{\beta}\nabla^{\beta}V_{\alpha}-C\nabla_{\alpha}C-CV_{\beta}B_{\alpha}^{\beta}\right)  &=- \nabla_{\alpha }\left(  \rho^{2}e_{\rho}\right) \end{align}$$

Note we in fact mixed the covariant and contraviant components $$\textstyle V_{\alpha}$$ and $$\textstyle V^{\alpha}$$ which we can now do more freely than before.