Transport Theorem for Spaces and Subspaces of Arbitrary Dimensions

: Using the apparatus of traditional differential geometry, the transport theorem is derived for the general case of a M -dimensional domain moving in a N -dimensional space, £ M N . The interesting concepts of curvatures and normals are illustrated with well-known examples of lines, surfaces and volumes. The special cases where either the space or the moving subdomain are material are discussed. Then, the transport at hypersurfaces of discontinuity is considered. Finally, the general local balance equations for continuum of arbitrary dimensions with discontinuities are derived.


Introduction
The transport theorem is a fundamental theorem used in formulating the basic conservation and balance laws in continuum mechanics (mass, momentum, and energy), which are adopted from classical mechanics and thermodynamics where the system approach is normally followed. Analogous to the classical Reynolds transport in continuum mechanics, the surface transport theorem is essential in the study of thin films undergoing large deformations, in epitaxial growth and in the study of phase boundary evolution. It is also important in the modeling of a singular surface which carries a certain structure of its own as it migrates. There is a vast literature on transport theorem and many references can be found in [1][2][3].
Betounes formulated and proved the general transport theorem associated with the motion of an arbitrary p -dimensional submanifold in a n -dimensional semi-Riemannian manifold [4]. He used the language and notation of modern differential geometry on manifolds e.g., [5,6], which is inconvenient for engineering and physics applications. Here, we formulate and prove the theorem using the language and concepts of traditional differential geometry and tensor calculus (e.g., [7,8]). Moreover, we apply the transport theorem to hypersurfaces of discontinuities and discuss the applications in continuum mechanics.
Petryk and Mroz derived the expressions for the first-and second-time derivatives of integrals and functionals defined on volume and surface domains which vary in time [9]. Their result is more general then the classical transport theorem as it pertains to piecewise regular surfaces and contains the edge terms. Cermelli et al. proved a transport theorem for smooth surfaces which evolve with time in Euclidean space, expressed in terms of the parameter-independent derivatives [10]. Recently, Sequin et al. extended the 3D transport theorem to rough domains of integration [11].
The need for the transport theorem arises in different contexts and consequently requires different derivation methods. The space-time approach [12] was used in [13] to derive the transport theorem for moving surface in a moving 3D region. A general transport theorem for moving surfaces based on the theory of generalized derivatives in n-dimensional space is presented in [14].
Two interesting attempts to present a unified approach to the topic of continuum mechanics on arbitrarily moving domains are given in [15,16]. They point out that it is desirable to formulate the transport theorem in a single unified way by using the classical approach expressed in terms of standards quantities from differential geometry and explicitly displaying the features that are common to all submanifolds, regardless of their finite dimensions.
This paper is organized as follows: In Sections 2 and 3, we consider geometry and kinematics in higher dimensions with special emphasis on the definitions of curvature and normals. Section 4 contains the derivation of the generalized transport theorem. In Section 5, we illustrate the concepts with well-known examples of lines, surfaces and volumes. In Section 6, we consider dependence of parametrization, i.e., on the choice of coordinates. In Sections 7-9, we consider the cases where the space and/or the moving subdomain are material in the sense of material in continuum mechanics. In Section 10, we consider a moving domain with hypersurfaces of discontinuity. Finally, in Section 11, we use the transport theorem to formulate the general local balance equations for continuum of arbitrary dimensions with discontinuities.

Geometry of
where or a u u we denote a typical point of ( ) Equivalently, ( , ) 0, 1, 2,..., where To illustrate dual representations (2) and (3), consider the representations of 1 ( ) V t in 3 E . The family of curves in 3 E , given in the parametric form may be represented as the intersection of circular cylinders: and as hyperbolic paraboloids: Each of the two representations has advantages and further we shall make use of both. While the representation (2) provides a convenient description of kinematics of ( ) M V t , it is dependent on parametrization. On the other hand, the representation (3) is independent of parameterization, i.e., independent of the choice of intrinsic coordinate system a u . It means that any transformation of intrinsic coordinate systems: does not change the representation (3). Consequently, the vectors independent of parameterization a u . Moreover, they are linearly independent because of (4), and may be taken as the base in Consider a different parametrization U such that From this, we have: are the basis vectors of ( ) M V t with respect to coordinate systems a u and U G , respectively.
Moreover, since ( , ) a k x u t in (2) satisfies (3), it follows that: whence we conclude that are also independent of parameterization. Moreover, from (8), it follows that: Further, we make use of vectors ( )   Then, vectors k g at the points x of ( ) M V t may be decomposed as: This relation is of crucial importance for decomposition of any tensor quantity defined on We will often make use of the relation between metric tensors of ( ) M V t defined in two coordinate systems a u and U L , i.e., the relation Then, The metric tensor

Kinematics of
Making use of (3), we obtain: . By d dt we denote the time derivative along a = u const . Note that: represents the scalar-normal velocity in ( ) p n direction which is independent of parameterization. Then, we can write: We will consider the coordinate system A U as a convected system in ( ) M V t . In continuum mechanics, the terms material or Lagrangian coordinates are used for convected coordinates, since the material particles can be labeled by these coordinates and, as such, they do not change their values during the motion of the material body. Note that, for a moving non-material domain, the choice of convected coordinates is arbitrary. However, once chosen, they remain fixed. The final result-the transport theorem-will be expressed in terms independent of the choice of coordinates.
At any time t , we consider the relation as the coordinate transformation between the convected coordinate system U L and intrinsic coordinate system a u . Further, we consider where C L are some constants. Generally, we denote by d dt the time derivative of any quantity In particular: From (5), we obtain the relation between velocities v and V of the point From (12), (10), and (7), we conclude that: thus proving that the scalar-normal velocity in ( ) p n direction is independent of parameterization. Instead of (13), a more compact representation is given by: Alternatively, We may also write: In particular, for the hypersurface

Generalized Transport Theorem in
, then according to [8] (p.

262):
Next, making use of (5), we write: Thus, in view of this and (15) we have: where M R is defined with respect to convected coordinates U L , and therefore the integration with respect to U L is independent of t . Now, we show in Appendix A, that: where The geometric meaning of ( ) which is independent of parameterization. Further, The expression

Remark 1.
To justify its name, we write in a unique compact invariant form as: where n is a unit vector orthogonal to ( ) M R t . Then, we have: . Accordingly, we have simple expression for V : The has been discussed previously in [10,15].
, it easy to see that: Therefore, After substituting (18) and (19) into (16), we obtain: Further, using divergence theorem we obtain (when the boundary consists of "material" points defined by convected coordinates, see Appendix B): is intrinsic to the motion of ). Finally, the transport theorem with respect to convective coordinate reads: We emphasize that the first integral is intrinsic to ( ) M R t , while the second integral is intrinsic with respect to ( ) ¶ M R t . Since the first integral on the right-hand side of (22) is invariant to any parameterization we may equally apply it for parameters a u . In this case, the scalar normal velocity of ( ) ¶ M R t with respect to parameters a u , i.e., needs to be used (Appendix B). Then, we write the generalized transport theorem with respect to non-convective coordinates a u in the form:

Examples
In this section, we consider familiar special cases and motivate the subsequent analysis of material domains and propagating discontinuity fronts. The section ends with the analysis of the capillary flow problem, which encompasses many of the special cases.
In all the cases below, the transport theorems holds also for the coordinates a u , i.e., (22), when we substitute , as defined in Appendix B.

3D Domain Moving in 3D Space
On the other hand, the familiar form of the transport theorem for the material body V t ( ) and the material field ( , ) t F x given per unit mass [ ( , ) t r x is mass density], reads: The volume V t ( ) in (24) is the material volume and the motion of the each material point

Closed Line Moving in a 3D Space
Define the unit tangent, normal, and binormal, , and , t m b in the standard way, the curvature As a special case, we compute the change of the total line energy arising from the line energy density In anticipation of the next example, we establish the relationship between the normal line velocity m V m and one of the surfaces. Let the surface with the normal 2 ( ) n be a material surface, moving with the velocity (2) w . The line then glides on this surface with the velocity 2 ν x ( ) relative to the surface (2). The unit normal-tangent vectors point outside the respective surfaces: The normal component of the line velocity the normal-tangent component relative to the surface (1) are then:

Interface Energies for Liquid Drop on a Solid Surface
The free liquid-gas interface ( ) a t has the unit normal 1  Then, from (25) and (28): where K is twice the mean curvature of the surface a t ( ) .

Capillary Flow
We consider an incompressible liquid flowing over a rigid solid and the surrounded by a gas with negligible viscosity and mass density, so that, without loss of generality, we assume the uniform vanishing pressure in the gas. The motion of the solid surface (2) w is prescribed and the triple line glides on the solid surface with the relative velocity The total energy of the system can be written as the sum of the bulk energy (kinetic + gravitational potential), the interface energies and the line energy: where y( ) x is the gravitational potential. Although the liquid-gas interface is not a material surface, the normal component of its velocity is identical to corresponding component of the material velocity of the fluid: Using (24), (26), (27) and (29), the rate of change of the total energy is where we have taken into account that for the rigid body motion (2) w : The last term in (30) is the correction arising from the rigid body motion of the whole assembly. To illustrate that point, consider a uniform translation . const (2) = w The last term in (30) is then: The liquid-gas surface contribution in (30) (for pure translation of the solid substrate) can then be written as: The energy balance requires that the total power input P be equal to the sum of the rate of total energy and dissipation rate D . We include the incompressibility condition, div , with the Lagrange multiplier field p which is recognized as the pressure in the fluid: where T represents the traction vector exerted by the solid on the liquid. We assume, without loss of generality, that the pressure in the (inviscid) gas vanishes, so that the traction on ( ) a t vanishes.
The flow at the re-entrant corner where the two surfaces intersect is singular (Taylor 1960). Within the sharp interface model the typical solution is to allow slip in some vicinity of the triple line [17]. Thus, at the solid surface we allow for slip, but not separation/penetration: The dissipation includes viscous dissipation in the bulk and the dissipation at the triple line: where τ is the viscous stress (deviatoric for incompressible fluid: tr 0 τ ( ) = ) and Q is the power conjugate of the triple line glide. Although we allow the slip at the solid surface, the dissipation (33) does not include this slip, owing to the choice of the external power input: ⋅ T v . Had we represented the external power as 2 ( ) ⋅ T w , the slip dissipation would have to be included. We will shortly see that the triple line dissipation is necessary if the experimentally observed difference between advancing and receding triple lines is to be described. The 2nd law of thermodynamics for isothermal processes reduces to the requirement that dissipation be positive. Assuming that the constitutive law for the fluid satisfies this requirement on its own (as it does for Newtonian fluids), it follows that the triple line force Q must have the same sign as the triple line glide x . The simplest linear constitutive law Upon substitution of (30) and (33) into (31) and some manipulation, we obtain the power balance equation: The simplest method for deriving the strong form of governing equations is to formulate the weak form directly via the Principle of virtual power (PVP) [18,19]. Application of the PVP yields the following governing equations and boundary conditions:

•
The natural boundary condition on ( ) A t : The governing (Cauchy) equations of motion in v t ( ) : The capillary jump in normal stress across ( ) a t : τ p K g -= I I. The deviatoric nature of the viscous stress τ then implies that on ( ) a t : The sign in the pressure jump condition is the consequence of the choice of 1 Moreover, from (35) for the advancing contact line Conversely, for the receding contact line the contact angle is smaller than the equilibrium contact angle. This is consistent with experimental observations [20]. Moreover, failure to include the triple line dissipation implies that the contact angle is always equal to the equilibrium one [21], which is in direct contradiction to experimental results.
Finally, we note that to complete the formulation, a definition of the slip constitutive law on

( )
A t is needed. This is beyond the scope of the current paper. We only note that such definition will not change any of the derived equations.

Transport Theorem Depending on Parametrization
In some cases, it is useful to use transport theorem in terms depending on parameterization. This can be done simply if we substitute (19), now written as: (20). It is convenient to put it in more familiar form: However: ( )

M V t in N-Dimensional Fluid
The transport relation given by (23) permits us to consider several particular cases of importance in continuum mechanics in a unified way. We consider N-dimensional fluid in analogy to 3dimensional fluid. First, we denote by ( , ) w x t the velocity of the fluid. We may write it as: where w a a = w tan a . Then, represent is the relative velocities of motion of ( ) M V t with respect to the fluid. We may decompose it as: Now, looking at the transport theorem the only term which may be influenced by motion of the fluid is . For definiteness we write: for the relative velocities of motion and normal migrational velocity of and hence: Upon substituting this into (23), we obtain: The corresponding transport theorem reads: ( ) Then, we note that: Therefore, the transport theorem becomes:

General Balance Laws
The basic laws of mechanics in 3D can all be expressed, in general, in the following form, for any bounded regular subregion P of the body B , and the vector field n , the outward unit normal to the boundary of the region t P in the current configuration. The quantities y and y s are tensor fields of order m , and y F is a tensor field of order 1 + m .
The relation (38) − the general balance of y in integral form − asserts that the rate of increase of the quantity y in a part P of a body is affected by the inflow of y through the boundary of P and the production of y within P . y F is the flux of y , and y s is the source of y . In general, the source y s may include external and internal sources.
We state general balance laws for where, according to (12), we have Thus, (A5) Formally, we may write (A4) as here, are the components of symmetric 2nd order tensor for each value of p [22]. Alternatively: Next, we calculate (A3) and obtain where we have made use of (A4) and (A5). Therefore (A2) is given by Next, we use (A4) and obtain Hence, Then, Grad F and its unit normal,