Compressible Navier-Stokes Equations in Cylindrical Passages and General Dynamics of Surfaces—(I)-Flow Structures and (II)-Analyzing Biomembranes under Static and Dynamic Conditions

: A new approach to solve the compressible Navier-Stokes equations in cylindrical co-ordinates using Geometric Algebra is proposed. This work was recently initiated by corresponding author of this current work, and in contrast due to a now complete geometrical analysis, particularly, two dimensionless parameters are now introduced whose correct deﬁnition depends on the scaling invariance of the N-S equations and the one parameter δ deﬁnes an equation in density which can be solved for in the tube, and a geometric Variational Calculus approach showing that the total energy of an existing wave vortex in the tube is made up of kinetic energy by vortex movement and internal energy produced by the friction against the wall of the tube. Density of a ﬂowing gas or vapour varies along the length of the tube due to frictional losses along the tube implying that there is a pressure loss and a corresponding density decrease. After reducing the N-S equations to a single PDE, it is here proven that a Hunter-Saxton wave vortex exists along the wall of the tube due to a vorticity argument. The reduced problem shows ﬁnite-time blowup as the two parameters δ and α approach zero. Finally we propose a CMS (Calculus of Moving Surfaces)–invariant variational calculus to analyze general dynamic surfaces of Riemannian 2-Manifolds in R 3 . Establishing ﬂuid structures in general compressible ﬂows and analyzing membranes in such ﬂows for example ﬂows with dynamic membranes immersed in ﬂuid (vapour or gas) with vorticity as, for example, in the lungs there can prove to be a strong connection between ﬂuid and solid mechanics.


Introduction
The fundamental lemma of Geometric-Algebra is that the geometric product can be decomposed into an inner product and outer product of two complex vectors. It can be proven that the bivector product of two complex vectors is a complex vector. With this in mind we introduce complex vectors for the representation of the inertial term f = a · ∇ a and the velocity b = a ρ . In the transformation of the governing Navier Stokes equations it can be proven that there is a density time dependent derivative present in the transformed equation which presents a pseudo-scalar in the transformed Navier-Stokes-Continuity equations. The stage is set to discuss solutions to two types of problems one which is in the setting of Geometric Algebra/Calculus using multivectors and the second type of solution being one in which the elements of the partial differential equations derived are not complex vectors(exponential function in time variable is not multivector dependent). We attempt to address one type of solution, the simpler one which leads to the existence of a Hunter-Saxton wave vortex in the large regimes of fluid flowing in a horizontal tube configuration. In the present work the full 3D compressible cylindrical Navier-Stokes equations are reduced using a new procedure to a Hunter-Saxton equation expressed in terms of the azimuthal angle in the flow. Also it is shown that there exists a complex vector F which can be expressed as F = f + Ic b with multivectors. Here we are interested in the angle in both the plane of rotation and outside the plane of rotation between the vectors b and f in the plane I.
Having introduced the goal of the paper for fluid dynamics of compressible flow in a tube at this stage of the paper, we are interested in further establishing a theoretical means to discuss the interaction of vortices with solid and flexible walls of tubes. Impingement of vortices on small biomembranes immersed inside tubes is also of general interest. These biological structures undergo large deformations under forces induced by vortices. Given this we now discuss a framework to analyze such membranes first in general and then with connection to fluid dynamics. It is important to do this in order to formulate a language to discuss fluid-membrane interaction problems in general. Variational Calculus is extremely useful in deriving laws regarding physical phenomena by minimizing the total energy required to form a system [17]. How such a system evolves may be determined by variating the integral of the total energy density known as the Action S of the system. This method of analyzing membranes has been often used in several physical fields such as physics and chemistry but is not seen as widespread in the Biomedical Sciences, particularly where it would have potential for analyzing Biological Phospholipid Bilayer Membranes [18]. The paper is organized into two parts, the first is (i)Flow Structures and the second is (ii) Analyzing Biomembranes under Static and Dynamic Conditions. The discussion consists of a relationship between these two parts with a conclusion connecting and summarizing these.

A New Composite Velocity Formulation
The 3D compressible cylindrical unsteady Navier-Stokes equations are written in expanded form, for each component, u r ,u θ and u z : where u r is the radial component of velocity, u θ is the azimuthal component and u z is the velocity component in the direction along tube, ρ is density, µ is dynamic viscosity ,Fg r , Fg θ , Fg z are body forces on fluid. The total gravity force vector is expressed as F T = (Fg r , Fg θ , Fg z ). The following relationships between starred and non-starred dimensional quantities together with a non-dimensional quantity δ are used: The density is defined in general as, αθ αθ ∈ (−π, 0), where i = √ −1, α > 0, ρ 1 is general function, and further below in this paper a geometric proof is shown that the sign of the radial velocity u r will change on a surface given by z * = r * 2 − θ * 2 .
It can also be observed that the density as defined above is real valued on (−π, π).

A Solution Procedure for δ Arbitrarily Small in Quantity
Multiplication of Equation (2.10) by ρ and Equation (2.12) below by L 1 , addition of the resulting equations [29], and using the ordinary product rule of differential multivariable calculus a form as in Equation (2.13) is obtained whereby a is given by a = ρ L 2 .
The continuity equation in cylindrical co-ordinates is Taking the geometric product in the previous equation with the inertial vector term, where b = a ρ is defined, where in the context of Geometric Algebra, the following scalar and vector grade equations arise, The geometric product of two vectors [30] is defined by A B = A · B + A × B. Taking the divergence of Equation (2.16) results in with the divergence of the following non-linear terms, , and which becomes where vorticity in the z * direction appears. It will be assumed for now that the vorticity is linear in r * and separable in r * from remaining independent variables. It can be generalized to a form r * F(r * , θ * , z * , t * ) with the restricton that F is increasing in r * towards the wall of the tube. For the moment, Div=0 will be set without justification and δ ≈ 0(but not zero) which imply, Substitution into the formula for u * r * above and solving for u * r * results in, The radial velocity, u * r * , depends on the terms u * θ * and hyperbolic part, r * 2 − θ * 2 − z * 2 . The sign of the radial velocity will change on the surface given by z * = r * 2 − θ * 2

Characterization of the Sign of the Vorticity
The cylindrical function (z * ) 2 = (r * ) 2 − (θ * ) 2 defines interior and exterior regions where the radial velocity will be positive and conversely negative; the radial velocity will vanish on regions which belong to the cylindrical surface outlined by the cylindrical function. Graphing modalities encounter difficulties due to plotting softwares protocols in handling square root functions. This cylindrical function has a natural representation which can be found by warping the cylindrical coordinate system from (r * , θ * ) to (φ, ξ) which is outlined by the substitution {r * = φ cosh ξ, θ * = φ}. This transformation has the surface jacobian outlined by: In this natural representation, the Radial Velocity null-surface is outlined by: where ξ max is specified to be coincident with a cylinder of radius r * = r max oriented along the z-axis at the value z * = z max . Since it is known that the conditions specify the equalities z max = φ max sinh ξ max and r max = φ max cosh ξ max , the value of ξ max can be given by This image can be plotted continuously (See Figures 1-3). An important property of this surface is that taking the cylindrical representation of the function, in the large regimes of r, the function reduces to a cone. This represents the quality that far from the center of the tube, the induced vorticity dissipates near the wall. At infinity for z * , we approach a perfect circle as shown partially in Figure 1.

Non-Linear Further Reduction
The other nonlinear terms in Equation (2.16), that is, 2r The value of r * is chosen to be large as mentioned previously for δ small in value. This implies that the above expression is negligible and can be omitted from the subsequent analysis. This is true since ω 3 is assumed to be linear in r * as will be further clarified below. It is known that the three dimensional vorticity controls the breakdown of smooth solutions of the 3D Euler equations [31].
We denote the following in Equation (2.15) As it will be shown in this work, for r * not very large so that the flow is confined to a central core in the tube, the term ψ is significant but not unbounded on t * ∈ (0, ∞). It is only when r * → ∞ that solutions breakdown outside the central core in the vicinity of the wall of the tube. So ψ dependent on ω 3 does control the breakdown of smooth solutions of the 3D compressible equations as will be seen in this study.
The expression f · Ψ in Equation (2.15) will vanish further below when we take a dot product in the z * direction of flow or k direction downstream in tube. We therefore do not include it after Equation(2.24) below. Upon multiplication of Equation (2.17) by, the resulting equation is which results upon using Equation (2.15) in, The continuity equation is written in terms of b as, where the following compact expression is given, For the term ρ f · b ∂ρ ∂t in Equation (2.24) we obtain upon using Equations (2.15) again, (2.26) and multiplying by ( f · f ) −1 in Equation (2.24) and, using properties of third derivatives involving the gradient and in particular the fact that the Laplacian of the divergence of a vector field is equivalent to the divergence of the Laplacian of a vector field, leads to the following form, where We consider Ω term now. We use Equation (2.15) and the following identities, Further in this paper we will take curl of the desired equation and the curl of the previous expression above when dotted with b will be zero. As a result the term that is left over is, If the Laplacian of the vorticity vector is the triple (-2,-2,-2)( ie. vorticity is quadratic), then using the following identity, This involves the vector projection of b onto f which is written in the conventional form,

Equation (2.28) can be written compactly as
where U f ξ is the scalar projection for b, G = 1 ρ 2 ∂ρ ∂t 2 , Q (a differential operator defined by Equations (2.16) and (2.28) and hence for a constant positive function α, with solution in terms of a function B, At this stage of the analysis we introduce the vorticity equation for compressible flow, We consider the third component of the vorticity equation in z * . It is assumed that the vorticity is an exponential function of z * and t * , . This can be non-zero except at the center of the tube where there the flow is irrotational, also G(θ * ) = 0 further away from the center of tube. Substitution of this form of ω into the vorticity equation and letting a * 3 be an increasing linear function in the form βz * , results in ∇ · a being solved for and expanded in a series for small α and large β, resulting in, It is worthy to note that ρ is decreasing down the tube and a * 3 is increasing and this is physically true for gas pipe flows.
Here we note that due to the appearance of the z * -direction of vorticity in the cylindrical Navier Stokes equations that this expression's form can control the breakdown of smooth solutions for 3-D cylindrical compressible equations as we will see below in this paper. The compressible Navier-Stokes equations have regular solutions that blow up in finite time. The remaining part of paper is to prove this using integral calculus methods. It will also be seen that if ρ is independent of t * then there are classical global solutions in t * .
Recall from the very onset definition of ρ that it is changing exponentially in time. It may be proven that B is written in terms of a non-homogeneous Green's function in spatial and time variables. Next we have from Equation (2.31), Using Equation (2.34) that is, ∇ · b = O(1/r * ) and substituting in the left side of Equation (2.35) and taking the limit as r * → ∞ for ρ as defined in this work, the left side of Equation (2.35) vanishes with the exception of term, − bG(ρ, ∂ρ ∂t * ) (cancels with exact term that is part of Q) and we obtain the following for Q, . Ω in Equation (2.28) vanishes due to assumption on rate of change of density with respect to t and c = 1/δ. Also, We obtain the following upon taking the curl of Equation (2.36), where the conservative gravity force drops out. Multiply Equation (2.37) by the normal vector cos(θ) a which is the normal component of a at wall of moving control volume (CV) in the z * direction,(direction of flow downstream in tube)

Stokes Theorem Applied to Dynamic Surfaces
Recalling Divergence theorem and Stoke's theorem, for general F, where C is the contour of a circle in control volume of tube and S consists of all surfaces of control volume. The geometry of the present problem is shown as a cross section of a tube in Figure 4 below.
Defining the following vector field, wheren = cos(θ) a and Stoke's theorem has been used. Applying Stoke's theorem to W, hence b. It is evident that the derivative of the circulation The last integral in previous series of equalities is an integral of a perfect differential around a closed path and is therefore equal to zero. To use Stokes theorem for dynamic surfaces that change with time, it can be seen in Appendix B that using the Calculus of Moving Surfaces (CMS) [18,32], the only requirement is that the paraboloid surface considered here has the same boundary as a disk which is coincident to the boundary of the paraboloid. Under this observation, the boundary is stationary and so the surface integral over a dynamic surface can be reduced down to the line integral around the stationary path, and is thus equivalent to the path integral around a closed disk coincident with the mouth of the paraboloid. Therefore, the dynamic nature of the paraboloids studied in this paper is not a point of concern when applying Stokes Equation to obtain a line integral. Refer to Appendix A for general time derivatives of dynamic line integrals. Proceeding we obtain, where T is unit tangent vector to closed curve C and ds is arc length, The third term in the parenthesis in Equation (2.44) is integrated by parts for line integral and we obtain the following,˛C Parametrizing the circle as r = g(θ) in polar coordinates it can be proven that the line integral in Equation (2.45) is, The normal form of Green's theorem can be used for the line integral in Equation (2.46), setting first, The line integral in Equation ( where M and N are given by Equations (2.47) and (2.48) respectively and R is the open disk with boundary C. The gravity force F T 1 which includes multiplication by δ 3 (see Equations (2.10),(2.11)), is expressed as follows, where S i is the sine integral, φ is a potential function, h is the negative height in the direction of the vector b g in Figure 4, T = 1 α is a time constant and g is gravity constant. Upon substitution of M and N into Equation (2.49) and further substitution of, gives us a very complicated pde to solve. However it is noteworthy to see that it is in fact separable in time t (Maple verification). The form is

The Hunter-Saxton Equation
Vorticity is assumed to be in the form ω 3 = r * F (θ * , z * , t * ) and in Equation (??) ρ avg is taken to be negligible, Differentiate with respect to θ * , use Equation (2.51), substitute r * 2 − θ * 2 − z * 2 = 0, so that we are on the circle at infinity and letting A denotes the partial derivative ∂b 2 ∂θ * and B the azimuthal velocity b 2 . Letting the entire expression above for the derivative of F T 2 term with respect to θ be denoted by G and scaling G by dividing by 4 −1 r * 3 z * 6 and taking the limit as z * approaches infinity gives, Substituting A = r * 2 A 2 and B = r * 2 B 2 where A 2 and B 2 are independent of r * as previously assumed due to vorticity being in this form as well(note we could have chosen A 2 ,B 2 to be dependent on r * with similar results to follow since the form of ω 3 can be generalized), we obtain for Equation (2.54) in the limit as α → 0 the following, Finally taking the limit as α → 0, since γ = 1 − 1 δ and matching r * large for δ small (recall r = r * δ), we obtain, It follows that the right hand side of ∂N ∂θ * is precisely 1 2 A 2 2 = 1 2 ( ∂b 2 ∂θ * ) 2 Using Equations (2.47-2.49) and integrating out r * parts and setting all constants to unity gives us a Hunter-Saxton-like equation at z * = ∞ on a perfect circle, that is, This equation is a Hunter-Saxton Equation that describes the evolution of the vorticity at the boundary of the tube. For very large regimes of the tube the fluid propagates as a Hunter-Saxton-like wave vortex along the boundary of the tube. It can be seen that b 1 = 0 and Equation (2.47) vanishes since F T 1 = 0 on the surface r * 2 − θ * 2 − z * 2 = 0. To transform the Hunter-Saxton Equation into cartesian co-ordinates, we use the chain rule, with the following transformation, To obtain the following evaluation of the derivatives: By multiplying and subtracting(adding) x and y it can be shown that the following commutator relations hold, Upon further differentiation it can be shown that, Using Equations (2.48) and (2.63) we obtain a pde which is separable in time and exhibits finite time blowup. The spatial part of the solution has been solved numerically. It is a simple exercise to show that waves at infinity occur showing that there is a vortex wave stucture there.

Non-Blowup Result
In Equation (2.55) it can be seen that if δ approaches infinity then this term approaches zero and we do not have a Hunter-Saxton equation. We have from Equations (??) and (2.13) that ρ will tend to zero and solutions will exist since only the viscosity term, pressure gradient and gravitational force terms will survive in Equation (2.13). To give the definition of density replacing Equations (??) and (??) for large δ, we write, this would prove positive for the existence of regular solutions to the Incompressible Cylindrical Navier-Stokes equations. It is immediate that as δ → ∞, ρ 1 → ρ avg = 0 in general.

Geometrical and Variational Analysis of the Hunter-Saxton Equation
The Hunter-Saxton equation mentioned above for b 2 (θ * , t * ) contains geometrical significance. Primarily, to define the region upon which the Hunter-Saxton equation is prescribed in this case, let Σ ⊂ R 2 represent an infinite cylindrical tube of radius R oriented along the z axis given by the following embedding in R 3 : For such a surface, a Lagrangian upon it would usually utilize L(θ * , z) but in this case it is sufficient to describe the Lagrangian as only dependent on θ * (and possibly time, t * ). The Lagrangian for the above formulation of the Hunter-Saxton Equation can be seen to be described by the following Surface Integral over the coordinate space of Σ denoted dR 2 = dθ * dzdt * of an Energy Density prescribed by L: Using Variational Calculus reproduces the following: which is the Hunter-Saxton Equation. The Lagrangian has difficulty being interpreted in the significance of the terms and how they would correspond to the geometry of the Surface it is defined on. If one considers the prototypical example of a cylindrical surface, a cylinder as prescribed above with Radius 1, it can be shown that the Lagrangian can be expressed in terms of CMS (Calculus of Moving Surfaces) objects on the surface Σ (the cylinder with radius 1) and expressing them through differential forms on b 2 = b 2 (θ * , t) by the following using the conventions of Ivancevic & Ivancevic [33]: where ·, · Σ is the inner product induced by the metric tensor on Σ, S 1 is the tangent vector to Σ in the θ direction (thus curves around the cylinder),∇ is the Invariant Time Derivative from CMS [18] defined on Σ, d (Σ) is the differential of a function on Σ, (·) is the "sharp" notation taking 1-forms to their corresponding 1-vector through the metric tensor on Σ and N is the Normal on Σ. It can be shown that while this holds for Σ as defined, it also holds for small dynamic deformations of Σ. Defining the original Lagrangian for static surfaces as L 0 , if Σ is dynamic, the Lagrangian must be modified to the following: where the Vector V || denotes the tangential velocity of the Surface. Through integrating the Lagrangian over the surface, the additional term becomes converted into a line integral of b 2 V || , n ∂Σ Σ around the boundary of the cylindrical tube ∂Σ where n ∂Σ is the normal to the boundary in the z direction (since n ∂Σ is perpendicular to both S 1 and N). It can be noted that the normal to the boundary ∂Σ is in the longitudinal direction whereas the tangential velocity strictly curves around the boundary of the tube and thus, the additional value vanishes. So: And thus produces the same variation. Therefore the Dynamic Nature of the Surface represents a gauge invariance in the Lagrangian which is fixed by the choice of the definition of the Invariant Time Derivative from CMS∇ = ∂ t − V || · ∇ Σ which maps tensors to tensors on Σ whilst the surface is deforming. Defining the Inner Product Norm on Σ, (·, ·) Σ =´Σ ·, · Σ dΣ, the Energy Density may be written in an alternate form: This may be simplified if we defined the 1-form A = d (Σ) b 2 and the 1-vector B = √ b 2 R · NA :

From this expression for the Energy Density, it is clear what the Variational Result which is the Hunter-Saxton Equation is
communicating. The First integrand is a method to calculate the work that the Kinetic Energy of the moving vortex is contributing to the total energy, while the second integrand describes the internal energy produced by the friction against Σ induced by the Hunter Saxton vortex wave solution to the Variational Equation. If the Hodge Star operator is defined on Σ, integrating the function over the surface, produces the following equivalent expression of the energy density using the adjoint nature of the Hodge Star to the inner product on Σ: In this way, the Hunter Saxton Equation for the current setting (regarding fluid dynamics) with vorticity is intuitively obtained.

The Static Local Energy Density
A summary of the techniques used to analyze Biomembranes in steady regimes is presented here. This usually involves a notion of some sort of total energy density H. It can also be shown that Hamiltonian's Variational Equations are tensorial if the Hamiltonian Density is tensorial itself [17]. Therefore, it can be easily seen why a Hamiltonian must be tensorial in its terms and over all scalar-valued

Geometrical Preliminaries, Defining the Biomembrane Continuum
The analysis of the Local Energy Density is assumed to be on a Riemannian 2-Manifold Σ that is locally isomorphic to R 2 . Since it is a smooth Riemannian Manifold the manifold possesses a Vector Space at every pointS ∈ Σ referred to as the Tangent Space atS, TSΣ with a dual co-tangent space T * S Σ. Tensors and Differential Forms atS ∈ Σ may be constructed as elements of the Tensor Product Spaces for a Tensor of rank (p, q) and as an element Exterior Product Space for a κ-Differential Form. The space of Tensors of rank (p, q) atS ∈ Σ, is denoted by T (p,q) S (Σ) and the space of κ-Forms constructed from the Exterior Product Spaces of either the Tangent Space or its dual are denoted as Λ κS (Σ) or Λ κ * S (Σ), respectively. They are constructed using the conventions from Ivancevic and Ivancevic [33] as such: It can be seen that since the exterior product ∧ is a restriction of the tensor product ⊗, the space of Differential Forms from the tangent/co-tangent exterior product spaces is a subset of the space of Tensors such that Λ κS (Σ) ⊂ T (Σ). Also since the distinction is made between cotangent/tangent exterior product spaces and tensor product spaces, the musical isomorphisms are defined as (·) and (·) where (·) : Λ κ * S (Σ) → Λ κS (Σ) and (·) : Λ κS (Σ) → Λ κ * S (Σ). Accordingly, several geometric tensors and forms from Hodge-DeRham Theory and Differential Geometry may be defined on Σ such as Σ's Surface 2-Form dΣ ∈ Λ 2 * S (Σ) such that dΣ = 1 2 ε αβ dS α ∧ dS β = det gSdS 1 ∧ dS 2 = Σ 1. dΣ is the Unique Co-tangent 2-form defined on the Riemannian 2-Manifold Σ where ε αβ is the components of the Levi-Civita cotangent 2-Form on Σ, det gS is the determinant of the Bilinear Symmetric (0,2)-Tensor gS(·, ·) ∈ T (0,2) S (Σ) defined such that gS : TSΣ × TSΣ → R atS ∈ Σ. This tensor is the Metric Tensor atS ∈ Σ and Σ is the Hodge Star Operator on Σ defined such that Σ : Λ κ * S (Σ) → Λ 2−κ S (Σ) when operating on cotangent differential κ-forms or equivalently, Σ : Λ κS (Σ) → Λ (2−κ) * S (Σ) when operating on tangent differential κ-forms.
Differential Operators may also be defined on Σ. The exterior derivatived : Λ κ * S (Σ) → Λ (κ+1) * S (Σ) may be defined and the co-differential defined explicitly byδ = (−1) Σd Σ may also be defined on Σ. Since the exterior derivative operates strictly on cotangent κ-forms, the codifferential must necessarily act on tangent κ-forms given by A ∈ Λ κS (Σ) and is generally defined asδ : Λ κS (Σ) → Λ κ−1 S (Σ) . Using these, the nominal Laplace-DeRham operator∆ : Λ κ * S (Σ) → Λ κ * S (Σ) may be defined on a cotangent κ-form A ∈ Λ κ * S (Σ) explicitly by∆A =d(δA ) + (δ(dA) ) . Since the manifold is restricted to R 2 , the operator can only be defined on cotangent 0-forms and 1-forms. It can be shown that for a 0-form (a function) f , sinceδ is the components of the covariant derivative on Σ using the notiation of ι A as the insertion operator on Σ. With these preliminary definitions, the Energy Density Action on Σ may be analyzed.

Introducing the Energy Density Action
On such a manifold, the Local Energy to be variated is typically formulated in the following general form: S =ˆΣ H(S αβ , B αβ )dΣ (3.2) where S αβ are the components of the Metric Tensor gS(·, ·) and B αβ are the components of the curvature tensor on Σ. It is assumed that ∂Σ = 0 and thus an analogue of Stokes Theorem which resembles Gauss' Divergence Theorem on Σ holds for a 1-form A ∈ Λ 1S (Σ) [18] : where n ∂Σ ∈ T (1,0) S (Σ) is the surface vector which points outwards at the boundary ∂Σ. Consequentially, as per this constraint, since ∂Σ = 0, then for all vectors V,´ΣδVdΣ = 0 for all Σ .

Properties and Requirements of the Hamiltonian and Its variation
The General Hamiltonian outlined in Equation (3.2) contains several interesting properties. First, it is important to note that so long as the Hamiltonian Energy Density H(S αβ , B αβ ) is constructed from tensorial geometric objects on Σ, it will be invariant to changes of coordinates S α = S α (S α ). In addition, the Energy Density must be scalar-valued. Finally the Static Hamiltonian (as the name suggests) does not variate in time, in the sense that As mentioned before, this indicates that this Hamiltonian Formulation is not suited for analyzing moving surfaces and if it is to be reformulated the new formulation must preserve all the scalar and invariant properties.
In component form, the Variation of the Energy Density with respect to the surface configuration R explicitly given by R = R + δR results in the following form [24,34]: where f α is given by: Its worth to note that this can be abbreviated using the Normal Calculus of Moving Surface outlined in Reference [35] by considering the vectors ξ µ = {S 1 , S 2 , N} as an orthonormal vector basis:

Stress Tensor
It can be noted that if the deformation δR is allowed to correspond to constant translation, δR = a, where ∇ α a = 0, then the variation vanishes δS = 0. On a closed surface using the translation δR = a, the variation reduces to: This reflects that the quantity f α is conserved and is accordingly defined as the Stress for the system [25]. Accordingly, the tensor f αµ = f α · ξ µ is defined as the Stress Tensor of the System as per the Normal CMS Conventions outlined in Reference [35]. This makes physical intuitive sense; In equilibrium, the stress of the Surface should be conserved and the Surface's according Equilibrium Configuration is determined by the Value of ∇ α f α .
In practice, only the Normal Projection of ∇ α f α is taken; this is because the tangential projection can be eliminated by a sufficient reparametrization of the Surface [25]. Therefore the equilibrium equation is as follows: Defining the symmetric (2,0) tensorf S = f αβ S α ⊗ S β , the curvature tensor byB = B αβ S α ⊗ S β and the 1-tensor f N = f α3 S α , the conservation equation can be rewritten in terms of Hodge-Theory and also in a compressed conservation form through the Normal Calculus of Moving Surfaces [35] by the following: or even simple Mean-Curvature squared Energy Densities (i.e., H = (B α α ) 2 ) which have Biological Applications [18,37,38]. However these Energy Densities only permit the determination of equilibrium, static configurations of surfaces and do not allow for their Dynamic Modelling.
An apparent issue to ensuring the invariance of the Energy Density is shed light on by the recent field of CMS. By trying to incorporate terms which to capture the speed of the surface such the derivative with respect to time of the surface's configuration V = ∂ t R or higher order acceleration terms such as ∂ 2 t R [32] and projections of the surface's speed V · S α , as well as historical measurements of the speed such as R · V and 1 2 V · V [34,39], these terms violate of the Energy Density's need to be invariant under coordinate transformations of dynamic membranes who's transformations are time dependent S α (t, S α ) [18,32,40]. This presents issues in deriving physically realistic quantities to describe the motion of membranes which are obtained by methods in classical mechanics. For this reason deriving scalar quantities which capture the motion of surfaces and are CMS-invariant is essential.

Invariant Scalars of Motion, Linear and Quadratic Invariants and Normal Speed Gradients
Before attempting to construct an Energy Density which is dynamic, scalars of motion which are CMS-invariant must be developed. Now the surface previously denoted by Σ must not be treated as Σ t indicating that it is in motion. All the Tensor Product Space and Exterior Product space constructions remain the same but are denoted as originating from Σ t . The jewel of CMS is the Invariant Time Derivative Operator which is defined as∇ : T (Σ t ) and captures the time variance of the tensor operated on. It is defined on a scalar field Φ ∈ T (0,0) S (Σ t ) by the following [18,32,40]: where V ∈ Λ 1 * S (Σ t ) is the surface speed defined as the speed of the 2-Manifold in R 3 , pulled back to the surface: It is important to notice how if the surface is not moving, V = 0 guarantees that∇ = ∂ t [18]. Thus,∇ only is a simple time derivative when acting on fields defined on Σ t which is not in motion or has an entirely normal motion of the form V = |V|N (as is the case with a uniformly expanding sphere/cylinder surface) [32]. Otherwise, the partial time derivative operator ∂ t does not produce tensors on Σ t , so objects which are constructed from the partial time derivative operator must either be eliminated from a hypothetical Dynamic Energy Density or must be treated in a CMS-manner to ensure the invariant-scalar nature of the Energy Density and therefore the ensuing Variational Conditions. The Invariant Time Derivative operator can also be extended to preserve the property that: Ultimately we conclude that V = ∂ t R is not a tensor and thus cannot be used in an Energy Density. Therefore∇ and its combination with other operators on Σ t must be used to identify invariant scalars of motion. It is known from CMS through the definition given in Equation (84), if R is the embedded configuration of the surface in R 3 , then∇R = CN, where the value C = V · N is referred to as the Normal Speed of Σ t . It is a Linear Invariant of Motion with respect to R, and so any function of C is also invariant and provides a manner of extending Invariants of Motion to higher-orders. As defined earlier, the Laplace-Beltrami Operator which is obtained from the Laplace-DeRham Operator∆ is Invariant itself. So this another linear invariant of motion is∆C =δ(dC) . Therefore the two Linear Invariants of Motion I 1 and I 2 on Σ t are given by: To continue searching for higher order invariants, the Commutator [∇, ∇ α ] =∇∇ α − ∇ α∇ can find use [32,35,40]. It is known that the commutator acts on Φ ∈ T (0,0) S (Σ t ) as follows: Applying the commutator to R, we see that S α∇ S α −d(CN) = CB and sincedN = −B, then we see∇S α = N ι S α (dC) confirming the CMS identity that∇S α = N∇ α C. It also proves that the 1-formdC = (∇ α C)dS α is as expected invariant. Thus the norm of this norm may be obtained as ||(dC) || 2 = gS((dC) , (dC) ) = S αβ (∇ α C)(∇ β C). Noting that∇(N · N) = 0 and∇(S α · N), it can be confirmed with CMS that∇N = −(dC) reaffirming the tensorial nature ofdC.

Temporal Curvature Trace Tensors
Thus far, the commutator [∇, ∇ α ] was restricted to 0-forms and it was shown that S α ([∇, ∇ α ])Φ = ι CB (dΦ) . Defining the cotangent 1-form operatorΞ : T (n,m) S it is seen thatΞΦ = 0 This operator can easily be extended to the Co-tangent Exterior Product Space to operate on 1-forms asΞ : . This will be referred to as the Temporal Curvature Operator, as it can be seen that when applied to a section of the tangent space or cotangent space, the Temporal Curvature TensorṘ ∈ T (1,2) S (Σ t ) arises [32].
As expected, Ξ α S α may be decomposed into a tangential and normal component. The normal projection confirms the linear nature of∆C while the surface projection confirms the identity of the Temporal Curvature Tensor's componentsṘ γ αβ already established in other sources [32,40]: 14) The Tensor admits three tracesṘ γ γβ ,Ṙ γ αγ ,Ṙ γ αβ S αβ and by the antisymmetric structure within the Tensor, and the symmetric nature of the Metric Tensor, it can be easily seen thatṘ γ αγ = 0. Thus there are only two non-trivial traces in the Temporal Curvature Tensor. These will be denoted by a tangent 1-form and a cotangent 1-form T = T α S α =Ṙ γ αβ S αβ S γ &W = W α dS α =Ṙ γ γα S α . It can be easily seen that the equivalence exists as T = −W . Therefore, the two traces are related by an isomorphism. This originates in the identity of the Temporal Curvature Tensor thatṘ [α βγ S β]γ = 0. Thus we define the fundamental Temporal Curvature Trace 1-Form as A natural Invariant of Motion which may be obtained from this is the 1-form's norm given by gS(W ,W ). Also sinceW ∈ Λ 1 * S (Σ t ), then another Linear Invariant of Motion isδW.

Developing a CMS-Invariant Hamiltonian, Variating the Hamiltonian and Divergence form of the Variation
Reflecting on all the CMS methods of measuring speeds on a membrane, a Dynamic Energy Density may be constructed effectively. The resulting form of the Hamiltonian Action chosen is as follows: S =ˆtˆΣ H(S αβ , B αβ , C)dΣ dt (3.16) For the purposes of this introduction so far, the specific form of H will remain non-given however this term will incorporate one of the Linear Invariant Scalars of Motion, C derived above. It is worthy to note that the Energy Density must have units of [J/m 2 ] = [N/m], so the Action has units of [J · s]. As stated earlier in the Static Case regarding´Σ dΣ, the integral´t dt ranges over the interval {t ∈ R|t 0 ≤ t ≤ t f } meaning ∂t = {t 0 , t f }. Performing a Variation of the Hamiltonian much like the previous case [24], the total variation results in: where f α is the standard stress tensor from the previous Local Hamiltonian Energy Density [25]. This variation is fundamentally different from the Local Hamiltonian. It can be simplifed by defining the following linear operatorδ t =∇ − CB α α defined by the relation: In this case, the variation can be simplified to the following: This equation may be put into a divergence form as outlined in Reference [24] by defining the Spatiotemporal Operator∇ b where the latin lowercase index b takes on values {1, 2, 3}. In this case for b = 1..2,∇ b = ι S αd and for b = 3,∇ b = −δ t . In this case, by defining in a similar manner the vector F b where for b = 1..2, F b = f α and for b = 3, F b = ∂H ∂C N, then the variation can be rephrased in the following manner: In the definition of the Spatiotemporal Operator∇ b , there is an negative sign in the third component∇ 3 = −δ t . To account for this negative sign, a Positive Spatiotemporal Operator∇ b may be defined as∇ c = σ bc∇ c by defining the 3 × 3 Surface-Time Matrix using components from gS(·, ·): It is noteworthy to note that with this Matrix, defining its determinant as σ = det(σ bc ), then the following identity arises σ = −det(S αβ ). In this case, a 3 dimensional Surfacial-Temporal Space can be defined where the volume element in such a space is defined as dσ = √ −σdS 1 ∧ dS 2 ∧ dt. This is an interesting analogue to the Minkowski Space given in Elementary Physics Textbooks [34,41] And so an integral in this space is formulated as such:σ dσ =ˆtˆΣ dΣ dt (3.18) Therefore, the variation can be unformly expressed in this 3-dimensional Surface-Time: This 3-Dimensional Surface-Time shares several parallels with a 4-Dimensional Space Time. This has several applications:

•
The formulation of Surface-Time can be used for formulating and analyzing the geometry of a 2-Dimensional analogue of Space-Time [41].

•
The formulation of Surface-Time can also be used for introducing concepts of Relativity into moving Membranes, biological or otherwise [26,27].

•
Lastly, this formulation of Surface-Time can be used to describe a similar method for describing 4D Space-Time in a manner consistent with CMS [28].

Conservation Laws and the Stress Tensor
As with the Static Energy Density, if the deformation δR = a, then the variation vanishes δS = 0 [24,25]. On a closed surface using the translation δR = a and assuming that the Surface Σ evaluated at ∂t are equal ( that the Surface Geometry at the end of the deformation in time does not change from the beginning), then the Variation reduces to: (3.20) This reflect that the quantity F b is conserved and is defined as the Stress for the system. Accordingly, the tensor F bµ = ξ µ · F b , is defined as the Stress Tensor of the System in accordance with [35]. Just like the case with the Static Energy Density, the Stress of the Surface is Conserved and the Surface's Dynamic Evolution is determined by the value of σ bc∇ b F c .
Like the Local Static case, this paper will only examine the Normal Projection of the Stress' Divergence. So therefore, the Dynamic Equilibrium Eqution is as follows:

Specific Lagrangian: Quadratic Speed Lagrangian
It is well known that the Standard Energy Density for measuring the kinetic speed of a continuum (Surface or Volume) is given by [18,26,27,32,39]: It is well shown by several sources that the Energy Density H = 1/2ρV · V is not tensorial on a time-deforming surface Σ t [18,32], due to the non-tensorial identity of V belonging to Σ t [40]. Thus, a modified CMS-invariant version is given by the Energy Density: where ρ ∈ R represents the (assumed) homogeneous mass-density of the surface being analyzed; this definition can be made precise by using continuum methods outlined in Reference [39]. Using this Energy Density, it can be obtained that: This is the equation of motion for a surface which does not have any Energy Density on it except for its own motion. This is interesting in that the dynamics of the Surface are identical regardless of the homogeneous mass sensity, analogous to Newtonian Physics [42]. It is noteworthy that if C=0, the equation is satisfied. This can be satisfied by two ways:

•
If the surface is not moving, V = 0. • If the surface is only moving tangentially V = V α S α .
These two solutions are somewhat trivial but can be easily imagined. For example, a sphere moving purely tangentially would be rolling along its surface. A sheet of paper moving tangentially would be it moving in its plane. It is intuitive why these would minimize the Energy Density. In addition, a non-trivial solution can be derived for the Sphere. A sphere which has a variable Radius, R = R(t) simplifies the equation to the following: The solution to this equation is a straightforward one For the sphere, this result implies a linear increase in Area knowing that the area for a sphere is A = 4πR 2 . Therefore a sphere which minimizes its kinetic energy density Hamlitonian will increase and decrease its Area linearly.
It is worthy to note that if the Hamiltonian was constrained by having a constant Volume and a Constant Surface Tension, the Action would have been modified to [24,25]: S =ˆtˆΣ HdΣ dt + µˆtˆΣ dΣ dt + PˆtˆΩ dΩ dt (3.26) where Ω Represents the Volume the Surface Encloses. Obtaining the variational equations for this Action results in the following:∇ For a Sphere, the equation simplifies to the following equation: The equation dictates the way that a Sphere will expand or contract in response to the Surface Tension µ, and Pressure on the Membrane P. In equilibrium with the surrounding, the Time Derivatives of the Radius Vanish and the equation results in the Young Laplace Equation [43]. In cases of non equilibrium, the equation results in showing how the surface will evolve to reach an equilibium by dictating the form of R(t). The equilibrium solution for the Variational Equation may be found by setting C = 0. This results in the following: This is a compact way of notating the Young-Laplace Law.
In the case where λ = 0, the Equilibrium Condition may be found by writing the Variational Equation out in full and observing whereṘ = 0; this will yield Boyle's Law (P · Volume = constant) [44] and if the converse is done, the equilibrium solution when P = 0 yields the Minimal Surface Equation where stationary solutions have B α α = 0 which is the famous solution to Plateau's Minimal Surface Problem [18].

Work-Energy Theorem
The Form of the Lagrangian's Variational Equations yields insights into the definition of the various Tensors derived; this framework allows for an analogue of the Work-Energy Theorem [45] to be formed on Σ. Under special circumstances, using an Energy Density of H = U (C), the Variational Equation becomes∇ dU This equation can be restated in the following form: Therefore, the integrand is a constant of motion conserved on Σ t . Since the only Energy on a Surface is derived from U , the Equation can be interpreted as a form of a conservation of Kinetic Energy.
This can be generalized for an Energy Density H = G(S αβ , B αβ ) + U (C). Assuming the above energy density, the following results: where g α is the analogue of the f α stress tensor derived for the Energy Density G(S αβ , B αβ ). This equation outlines a relationship between the Kinetic Energy derived from U (C), and the work exerted on the surface by CN · ∇ α g α . This is a form of the Work Energy Theorem for Surfaces [45]. Therefore the equation is interpreted as the following: It is worthy to note that normally, Mechanical Power of a moving body is calculated using its local continuum velocity v(x) [39]; the analogue of this work for Surfaces V(S), has been shown by Grinfeld to be Non-tensorial for Surfaces [18,40]. However allowing the surface to move, and assuming a general Energy Density H(S αβ , B αβ , C), results in the identification of the expression CN. This is the only component of the Ambient Surface Velocity V = CN + V α S α that is Tensorial; therefore it makes intuitive sense that this is the only component taken into account for Mechanical Power which is commonly notated in Continuum Mechanics literature as´V v · b dV where b are all the body forces acting on a Body. The Dynamic Framework presented up to now supports some possible extensions.

Discussion
It has been earlier established in Section 2.8 that when modelling vorticity on a cylindrical tube, the evolution of the vorticity at the cylindrical surface Σ t can be obtained with the following Energy Density: where (A, B) =´Σ gS (A, B)dΣ. This Energy Density is seen to have the same stucture as the energy densities explored in Part II of this paper. When the Energy Density is variated with respect to b 2 (of Part I of paper), it results in the Hunter Saxton Equation.
Defining the 1-form A =db 2 , and the 1-vector B = √ b 2 R · NA : Through the Energy Density, it can be seen that the vorticity b 2 interacts with Σ t through the Metric Tensor, and through the Normal. In this particular interaction, Σ t = Σ is Stationary and referred to the static cylindrical tube. In the case where Σ t is allowed to variate, the motion of a membrane with a Vorticity that is pre-prescribed at Σ t such as b 2 = ω(R(S α )) can be obtained through Variating the Energy Density with respect to R. In this case, the Vorticity would affect the surface which minimizes the Energy Caused by the interaction. This Minimization can be explored in the future when desiring to observe the motion of a dynamic membrane immersed in a fluid with Vorticity.

Conclusions
A reduction of the compressible Navier Stokes equations coupled to the continuity equation in cylindrical co-ordinates to a simpler problem has been shown. Dimensionless parameters were introduced whereby in the small limit case for both of these a method of solution is sought for in the tube. A wave pulse that travels downwards towards the wall exists in Figure 1 [46]. A Dirac Delta density pulse exists exactly at the entrance of the tube and the density decreases downstream in the tube due to frictional and pressure losses. The density is a Gaussian-like wave function with an exponential in time. The result is finite time blowup for the velocity in the azimuthal direction. The Hunter-Saxton equation is a special case appearing at z * = ∞ in θ * ,t * . Future studies are required for parallel methods derived from this paper that can be applied to the problem of the Clay Institute for the Navier-Stokes Incompressible equations. The second objective of this paper was to develop a Variational Framework which is Dynamic in order to potentially model Biological Membranes. Several objectives were acheived:

•
The standard Variational Framework proposed in References [24,25] N · ∇ α f α = 0 was extended to Dynamic Surfaces by using the Calculus of Moving Surfaces (CMS) developed in Grinfeld's Textbook [18]. • Drawing, on previous extensions of CMS [32,35], several Invariants of Motion have been developed in this paper (C, Q α β ∇ α ∇ β C, P α β ∇ α C∇ β C) which may be used in the Dynamic extension of the Hamiltonian; however, in this paper, the only invariant of motion considered was C for simplicity. Line Integrals are difficult to evaluate when they are restricted to the boundary ∂Σ t of an open surface Σ t that is moving dynamically in time t; in several fields, but particularly in the Calculus of Moving Surfaces (CMS), Σ is usually parametrized by the vector function R Σ = R(S, t), where S = S α = {S 1 , S 2 } represent the surface variables chosen to parametrize the Surface. The parametrization of ∂Σ t is usually given by R λ = R(U, t), where U = U ψ = {U 1 } represents the variable used to parametrize ∂Σ t . Since the line is restricted to the boundary of the surface, R Σ is related to R λ through the composition R λ = R(S(U), t). In this case, S(U) = {S 1 (U), S 2 (U)} refers to the relation of the surface variables to the line variable. Consequentially, derivatives with respect to time of these line integrals are difficult to evaluate as well. The line integral of a vector function F taken over the line variable interval U 1 ∈ [a, b] is simplified to: where T is the unit tangent vector and U = U φψ is the metric tensor on ∂Σ t ; all naming conventions on the line and surface are outlined in Reference [18]. Using CMS, the time derivative can pass inside the integral since it is an single integral taken over a variable independent of t. Thus, the following develops: It is important to note that the partial derivative in the second term measures how much the tangent vector of the curve variates in time. This is evaluated further using tensors on the surface given in Reference [32] d if the tensorΓ =Γ β α (S α ⊗ S β ) is defined, then the following is given for the time derivative of a line integral:

Appendix B. Using Stokes Theorem on Dynamic Surfaces & Special Cases
In general, Stokes Theorem is not typically able to be extended to open surfaces who are dynamic in time. It is stated earlier that From CMS, it is known that an analogue of Gauss' Divergence Theorem exists on open surfaces: Framing the above as a classic line integral results in the familiar form where ∇ α is the component representation of the covariant derivative on the surface Σ. Knowing that ∇ α F = (S α · ∇)F, the above equation can be reduced to Stokes Equation, and so the equality can be viewed as an extension of Stokes Theorem to surfaces which are in motion. Utilizing the formula for the derivative of a line integral and substituting the right hand side in for the vector function, the following is obtained: Again recalling that ∇ α F = (S α · ∇)F, the surface integral can be broken up into two surface integrals: This is in general the relationship between the derivative of a line integral and the equivalent surface integral assuming the surface is in motion. It is worth noting as expected that the Tensor which captures the speed of the surface (and therefore the speed of the boundary),Γ is included in the surface integral as expected. A special case worth noting is that if in fact the speed on the surface evaluated at the boundary is stationary (i.e., a surface deforming with a fixed boundary ∂Σ t = ∂Σ), then it can be seen thatΓ = 0 and thus the equation reduces to the desired result: where the time derivative is taken over a stationary path ∂Σ, but the surface integral is taken over a dynamic surface Σ t . Therefore Stokes Theorem is valid for a surface that is moving in time, provided that the boundary is stationary. This is useful to evaluate integrals which are in motion as well.