Calculus of moving surfaces

The surface of a flag in the wind is an example of a deforming manifold.

The calculus of moving surfaces (CMS) [1] is an extension of the classical tensor calculus to deforming manifolds. Central to the CMS is the \delta /\delta t-derivative whose original definition [2] was put forth by Jacques Hadamard. It plays the role analogous to that of the covariant derivative \nabla _{\alpha } on differential manifolds. In particular, it has the property that it produces a tensor when applied to a tensor.

Jacques Salomon Hadamard, French Mathematician, 1865–1963 CE

Suppose that S_t is the evolution of the surface S indexed by a time-like parameter t. The definitions of the surface velocity C and the operator \delta /\delta t are the geometric foundations of the CMS. The velocity C is the rate of deformation of the surface S in the instantaneous normal direction. The value of C at a point P is defined as the limit

C=\lim_{h\to 0} \frac{\text{Distance}(P,P^*)}{h}

where P^{*} is the point on S_{t+h} that lies on the straight line perpendicular to S_{t} at point P. This definition is illustrated in the first geometric figure below. The velocity C is a signed quantity: it is positive when \overline{PP^{*}} points in the direction of the chosen normal, and negative otherwise. The relationship between S_{t} and C is analogous to the relationship between location and velocity in elementary calculus: knowing either quantity allows one to construct the other by differentiation or integration.

Geometric construction of the surface velocity C
Geometric construction of the \delta/\delta t-derivative of an invariant field F

The \delta /\delta t-derivative for a scalar field F defined on S_{t} is the rate of change in F in the instantaneously normal direction:

\frac{\delta F}{\delta t}=\lim_{h\to 0} \frac{F(P^*)-F(P)}{h}

This definition is also illustrated in second geometric figure.

The above definitions are geometric. In analytical settings, direct application of these definitions may not be possible. The CMS gives analytical definitions of C and \delta /\delta t in terms of elementary operations from calculus and differential geometry.

Analytical definitions

For analytical definitions of C and \delta /\delta t, consider the evolution of S given by

Z^i = Z^i \left( t ,S \right) \,

where Z^{i} are general curvilinear space coordinates and S^{\alpha } are the surface coordinates. By convention, tensor indices of function arguments are dropped. Thus the above equations contains S rather than S^\alpha.The velocity object v^{i} is defined as the partial derivative

v^i =\frac{\partial Z^i \left( t ,S \right)}{\partial t }

The velocity C can be computed most directly by the formula

C=v^i N_i \,

where N_i are the covariant components of the normal vector \vec{N}.

The definition of the \delta /\delta t-derivative for an invariant F reads

\frac{\delta F}{\delta t}=\frac{\partial F\left( t ,S \right)}{\partial t }-v^{i}Z^{\alpha }_{i}\nabla _{\alpha }F

where Z^\alpha_i is the shift tensor and \nabla_\alpha is the covariant derivative on S.

For tensors, an appropriate generalization is needed. The proper definition for a representative tensor T^{i\alpha }_{j\beta } reads

\frac{\delta T^{i\alpha }_{j\beta }}{\delta t}=\frac{\partial T^{i\alpha }_{j\beta }}{\partial t}-v^{\eta }\nabla _{\eta }T^{i\alpha }_{j\beta }+v^{m}\Gamma ^{i}_{mk}T^{k\alpha }_{j\beta }-v^{m}\Gamma ^{k}_{mj}T^{i\alpha }_{k\beta }+\nabla _{\eta }v^{\alpha }T^{i\eta }_{j\beta }-\nabla _{\beta }v^{\eta }T^{i\alpha }_{j\eta }

where \Gamma^k_{mj} are Christoffel symbols.

Properties of the δ/δt-derivative

The \delta /\delta t-derivative commutes with contraction, satisfies the product rule for any collection of indices

\frac{\delta }{\delta t}\left( S^i_\alpha T^\beta_j
 \right)=\frac{\delta S^i_\alpha}{\delta t}T^\beta_j + S^i_\alpha \frac{\delta T^\beta_j }{\delta t}

and obeys a chain rule for surface restrictions of spatial tensors:

\frac{\delta F^j_k}{\delta t} =\frac{\partial F^j_k}{\partial t}+CN^i \nabla _i F^j_k

Chain rule shows that the \delta /\delta t-derivative of spatial "metrics" vanishes

\frac{\delta \delta^i_j}{\delta t},\frac{\delta Z_{ij}}{\delta t},\frac{\delta Z^{ij}}{\delta t},\frac{\delta \varepsilon _{ijk}}{\delta t},\frac{\delta \varepsilon^{ijk}}{\delta t}=0

where Z_{ij} and Z^{ij} are covariant and contravariant metric tensors, \delta ^{i}_{j} is the Kronecker delta symbol, and \varepsilon _{ijk} and \varepsilon ^{ijk} are the Levi-Civita symbols. The main article on Levi-Civita symbols describes them for Cartesian coordinate systems. The preceding rule is valid in general coordinates, where the definition of the Levi-Civita symbols must include the square root of the determinant of the covariant metric tensor Z_{ij}.

Differentiation table for the δ/δt-derivative

The \delta /\delta t-derivative of the key surface objects leads to highly concise and attractive formulas. When applied to the covariant surface metric tensor S_{\alpha \beta } and the contravariant metric tensor S^{\alpha \beta }, the following identities result

\begin{align}
\frac{\delta S_{\alpha \beta }}{\delta t} & = -2CB_{\alpha \beta } \\[8pt]
\frac{\delta S^{\alpha \beta }}{\delta t} & = 2CB^{\alpha \beta }
\end{align}

where B_{\alpha \beta } and B^{\alpha \beta } are the doubly covariant and doubly contravariant curvature tensors. These curvature tensors, as well as for the mixed curvature tensor B^\alpha_\beta, satisfy

\begin{align}
\frac{\delta B_{\alpha \beta }}{\delta t}& = \nabla _\alpha \nabla_\beta C - CB_{\alpha \gamma }B^\gamma_\beta \\[8pt]
\frac{\delta B^\alpha_\beta}{\delta t}& = \nabla^\alpha \nabla_\beta C + CB^\alpha_\gamma B^\gamma_\beta \\[8pt]
\frac{\delta B^{\alpha \beta }}{\delta t}& = \nabla ^\alpha \nabla^\beta C + 3CB^\alpha_\gamma B^{\gamma \beta}
\end{align}

The shift tensor Z^i_\alpha and the normal N^i satisfy

\begin{align}
\frac{\delta Z^i_\alpha}{\delta t} & = \nabla _\alpha \left( CN^i \right) \\[8pt]
\frac{\delta N^i}{\delta t} & = -Z^i_\alpha \nabla^\alpha C
\end{align}

Finally, the surface Levi-Civita symbols \varepsilon _{\alpha \beta } and \varepsilon ^{\alpha \beta } satisfy

\begin{align}
\frac{\delta \varepsilon _{\alpha \beta }}{\delta t} & = -\varepsilon _{\alpha \beta }CB^{\gamma }_{\gamma } \\[8pt]
\frac{\delta \varepsilon ^{\alpha \beta }}{\delta t} & = \varepsilon ^{\alpha \beta }CB^\gamma_\gamma
\end{align}

Time differentiation of integrals

The CMS provides rules for time differentiation of volume and surface integrals.

References

  1. Grinfeld, P. (2010). "Hamiltonian Dynamic Equations for Fluid Films". Studies in Applied Mathematics. doi:10.1111/j.1467-9590.2010.00485.x. ISSN 00222526.
  2. J. Hadamard, Lecons Sur La Propagation Des Ondes et Les Equations De l’Hydrodynamique. Paris: Hermann, 1903.