Calculus of moving surfaces

Short description: Extension of the classical tensor calculus

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

The calculus of moving surfaces (CMS) is an extension of the classical tensor calculus to deforming manifolds. Central to the CMS is the Tensorial Time Derivative [math]\displaystyle{ \dot{\nabla} }[/math] whose original definition was put forth by Jacques Hadamard. It plays the role analogous to that of the covariant derivative [math]\displaystyle{ \nabla _{\alpha } }[/math] on differential manifolds in that it produces a tensor when applied to a tensor.

Jacques Salomon Hadamard, French Mathematician, 1865–1963 CE

Suppose that [math]\displaystyle{ \Sigma_t }[/math] is the evolution of the surface [math]\displaystyle{ \Sigma }[/math] indexed by a time-like parameter [math]\displaystyle{ t }[/math]. The definitions of the surface velocity [math]\displaystyle{ C }[/math] and the operator [math]\displaystyle{ \dot{\nabla} }[/math] are the geometric foundations of the CMS. The velocity C is the rate of deformation of the surface [math]\displaystyle{ \Sigma }[/math] in the instantaneous normal direction. The value of [math]\displaystyle{ C }[/math] at a point [math]\displaystyle{ P }[/math] is defined as the limit

[math]\displaystyle{ C=\lim_{h\to 0} \frac{\text{Distance}(P,P^*)}{h} }[/math]

where [math]\displaystyle{ P^{*} }[/math] is the point on [math]\displaystyle{ \Sigma_{t+h} }[/math] that lies on the straight line perpendicular to [math]\displaystyle{ \Sigma_{t} }[/math] at point P. This definition is illustrated in the first geometric figure below. The velocity [math]\displaystyle{ C }[/math] is a signed quantity: it is positive when [math]\displaystyle{ \overline{PP^{*}} }[/math] points in the direction of the chosen normal, and negative otherwise. The relationship between [math]\displaystyle{ \Sigma_{t} }[/math] and [math]\displaystyle{ C }[/math] 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 [math]\displaystyle{ \delta/\delta t }[/math]-derivative of an invariant field F

The Tensorial Time Derivative [math]\displaystyle{ \dot{\nabla} }[/math] for a scalar field F defined on [math]\displaystyle{ \Sigma_{t} }[/math] is the rate of change in [math]\displaystyle{ F }[/math] in the instantaneously normal direction:

[math]\displaystyle{ \frac{\delta F}{\delta t}=\lim_{h\to 0} \frac{F(P^*)-F(P)}{h} }[/math]

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 [math]\displaystyle{ \dot{\nabla} }[/math] in terms of elementary operations from calculus and differential geometry.

Analytical definitions

For analytical definitions of [math]\displaystyle{ C }[/math] and [math]\displaystyle{ \dot{\nabla} }[/math], consider the evolution of [math]\displaystyle{ S }[/math] given by

[math]\displaystyle{ Z^i = Z^i \left( t ,S \right) }[/math]

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

[math]\displaystyle{ V^i =\frac{\partial Z^i \left( t ,S \right)}{\partial t } }[/math]

The velocity [math]\displaystyle{ C }[/math] can be computed most directly by the formula

[math]\displaystyle{ C=V^i N_i }[/math]

where [math]\displaystyle{ N_i }[/math] are the covariant components of the normal vector [math]\displaystyle{ \vec{N} }[/math].

Also, defining the shift tensor representation of the Surface's Tangent Space [math]\displaystyle{ Z^{\alpha}_{i}=\textbf{S}^\alpha\cdot\textbf{Z}_i }[/math] and the Tangent Velocity as [math]\displaystyle{ V^\alpha=Z^\alpha_iV^i }[/math] , then the definition of the [math]\displaystyle{ \dot{\nabla} }[/math] derivative for an invariant F reads

[math]\displaystyle{ \dot{\nabla}F=\frac{\partial F\left( t ,S \right)}{\partial t }-V^\alpha\nabla _{\alpha }F }[/math]

where [math]\displaystyle{ \nabla_\alpha }[/math] is the covariant derivative on S.

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

[math]\displaystyle{ \dot{\nabla}T^{i\alpha }_{j\beta }=\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 }+\dot{\Gamma}^\alpha_\eta T^{i\eta }_{j\beta }-\dot{\Gamma}^\eta_\beta T^{i\alpha }_{j\eta } }[/math]

where [math]\displaystyle{ \Gamma^k_{mj} }[/math] are Christoffel symbols and [math]\displaystyle{ \dot{\Gamma}^\alpha_\beta=\nabla_\beta V^\alpha-C B^\alpha_\beta }[/math] is the surface's appropriate temporal symbols ([math]\displaystyle{ B^\alpha_\beta }[/math] is a matrix representation of the surface's curvature shape operator)

Properties of the [math]\displaystyle{ \dot{\nabla} }[/math]-derivative

The [math]\displaystyle{ \dot{\nabla} }[/math]-derivative commutes with contraction, satisfies the product rule for any collection of indices

[math]\displaystyle{ \dot{\nabla}(S^i_\alpha T^\beta_j )=T^\beta_j\dot{\nabla}S^i_\alpha + S^i_\alpha\dot{\nabla}T^\beta_j }[/math]

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

[math]\displaystyle{ \dot{\nabla}F^j_k(Z,t) =\frac{\partial F^j_k}{\partial t}+CN^i \nabla _i F^j_k }[/math]

Chain rule shows that the [math]\displaystyle{ \dot{\nabla} }[/math]-derivatives of spatial "metrics" vanishes

[math]\displaystyle{ \dot{\nabla}\delta^i_j=0,\dot{\nabla}Z_{ij}=0,\dot{\nabla}Z^{ij}=0,\dot{\nabla}\varepsilon _{ijk}=0,\dot{\nabla}\varepsilon^{ijk}=0 }[/math]

where [math]\displaystyle{ Z_{ij} }[/math] and [math]\displaystyle{ Z^{ij} }[/math] are covariant and contravariant metric tensors, [math]\displaystyle{ \delta ^{i}_{j} }[/math] is the Kronecker delta symbol, and [math]\displaystyle{ \varepsilon _{ijk} }[/math] and [math]\displaystyle{ \varepsilon ^{ijk} }[/math] 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 [math]\displaystyle{ Z_{ij} }[/math].

Differentiation table for the [math]\displaystyle{ \dot{\nabla} }[/math]-derivative

The [math]\displaystyle{ \dot{\nabla} }[/math] derivative of the key surface objects leads to highly concise and attractive formulas. When applied to the covariant surface metric tensor [math]\displaystyle{ S_{\alpha \beta } }[/math] and the contravariant metric tensor [math]\displaystyle{ S^{\alpha \beta } }[/math], the following identities result

[math]\displaystyle{ \begin{align} \dot{\nabla}S_{\alpha \beta } &=0 \\[8pt] \dot{\nabla}S^{\alpha \beta } &=0 \end{align} }[/math]

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

[math]\displaystyle{ \begin{align} \dot{\nabla}B_{\alpha \beta }& = \nabla_\alpha \nabla_\beta C + CB_{\alpha \gamma }B^\gamma_\beta \\[8pt] \dot{\nabla}B^\alpha_\beta& = \nabla_\beta\nabla^\alpha C + CB^\alpha_\gamma B^\gamma_\beta \\[8pt] \dot{\nabla}B^{\alpha \beta }& = \nabla ^\alpha \nabla^\beta C + CB^{\gamma \alpha}B^\beta_\gamma \end{align} }[/math]

The shift tensor [math]\displaystyle{ Z^i_\alpha }[/math] and the normal[math]\displaystyle{ N^i }[/math] satisfy

[math]\displaystyle{ \begin{align} \dot{\nabla}Z^i_\alpha & = N^i\nabla _\alpha C \\[8pt] \dot{\nabla}N^i & = -Z^i_\alpha \nabla^\alpha C \end{align} }[/math]

Finally, the surface Levi-Civita symbols [math]\displaystyle{ \varepsilon _{\alpha \beta } }[/math] and [math]\displaystyle{ \varepsilon ^{\alpha \beta } }[/math] satisfy

[math]\displaystyle{ \begin{align} \dot{\nabla} \varepsilon _{\alpha \beta } & = 0 \\[8pt] \dot{\nabla} \varepsilon ^{\alpha \beta } & = 0 \end{align} }[/math]

Time differentiation of integrals

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