Next Article in Journal
Estimates for the Commutators of p-Adic Hausdorff Operator on Herz-Morrey Spaces
Previous Article in Journal
Turning Hild’s Sculptures into Single-Sided Surfaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra

by
Terry E. Moschandreou
1,2
1
London International Academy, 365 Richmond Street, London, ON N6A 3C2, Canada
2
Department of Applied Mathematics, Faculty of Science, Western University, London, ON N6A 5C1, Canada
Mathematics 2019, 7(2), 126; https://doi.org/10.3390/math7020126
Submission received: 5 December 2018 / Revised: 18 January 2019 / Accepted: 23 January 2019 / Published: 27 January 2019

Abstract

:
A method of solution to solve the compressible unsteady 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity equation in cylindrical coordinates is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the tube. A reduction to a single partial differential equation is possible and integral calculus methods are applied for the case of a body force in the direction of gravity to obtain an integral form of the Hunter-Saxton equation.

1. Introduction

Compressible flow has many applications some of which are of physics, mathematics and engineering interest. In general we have two types of flows, internal and external. Internal flows in ducts are important in industry and nozzles and diffusers used in engines are also an applied area for these types of flows. In general, density changes are related to temperature changes. External flows can be important for airplanes and projectiles where compressiblilty effects are important. The Navier Stokes equations have been dealt with extensively in the literature for both analytical [1] and numerical solutions [2,3]. Some work converting the compressible Navier-Stokes equations to the Schrödinger equation in quantum mechanics by means of transformations has been carried out by Vadasz [4]. General mathematical and computational methods for compressible flow have been outlined in [5]. Methods in more general fluid mechanics are also addressed in [6]. In the context of functional analysis it has been shown in [7] that generally the motion of a compressible fluid with fixed initial velocity field and constant initial density converges to that of an incompressible fluid as its sound speed goes to infinity. Some analytical methods such as in [8] have been successfully carried out for one and two dimensional isentropic unsteady compressible flow. In the present work I introduce a new procedure to write the compressible unsteady Navier Stokes equations with a general spatial and temporal varying density term in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. This procedure can be used in physics and engineering to simplify a complex system of PDE to a simpler one such as complex multiphase flows [9] and other areas in physics using Geometric Algebra, such as for example the study of the Maxwell equations. In mathematics, the method leads to the present work which shows the necessary conditions for the full solution of the system of equations for compressible flow to blowup at a certain time for certain types of initial conditions and also implications of using assumptions for necessary and physically meaningful functional forms of density. Following the general procedure, a dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the tube. It is concluded that the total divergence of the flow can be expressed as the integral with respect to time of the line integral of the dot product of inertial and azimuthal velocity. The line integral can be evaluated on a contour that is annular and traces the boundary layer as time increases in the flow. A reduction to a single partial differential equation is possible and integral calculus methods are applied for the case of a body force directed downwards with gravity to obtain an integral form of the Hunter-Saxton equation. A specific density function that is exponential in time and having characteristics of a pulse that propagates radially and axially in the tube and can be interpreted as having a large density gradient in time and in space is used.

2. 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 :
u r t + u r u r r + u θ r u r θ u θ 2 r + u z u r z μ ρ u r r 2 + 2 u r r 2 + 1 r u r r + 1 r 2 2 u r θ 2 2 r 2 u θ θ + 2 u r z 2 μ 3 ρ r u r r + 1 r u r θ + u r z + 1 ρ p r F g r = 0
u θ t + u r u θ r + u θ r u θ θ + 1 r u r u θ + u z u θ z μ ρ u θ r 2 + 2 u θ r 2 + 1 r u θ r + 1 r 2 2 u θ θ 2 + 2 r 2 u r θ + 2 u θ z 2 1 r μ 3 ρ θ u θ r + 1 r u θ θ + u θ z + 1 ρ r p θ F g θ = 0
u z t + u r u z r + u θ r u z θ + u z u z z μ ρ 2 u z r 2 + 1 r u z r + 1 r 2 2 u z θ 2 + 2 u z z 2 μ 3 ρ z u z r + 1 r u z θ + u z z + 1 ρ p z F g z = 0
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, F g r , F g θ , F g z are body forces on fluid. The total gravity force vector is expressed as F T = ( F g r , F g θ , F g z ) .
The following relationships between starred and non-starred dimensional quantities together with a non-dimensional quantity δ are used:
u r = 1 δ 2 u r *
u θ = 1 δ u θ *
u z = 1 δ u z *
r = r * δ
θ = θ *
z = 1 δ z * , t = δ t *
c a constant is defined below, δ = 1 r * ρ θ * ρ r * + ρ r * 1 as part of the continuity equation (Equation (12)). The solution of this equation in terms of ρ is of the form f ρ θ δ + ln r δ , z , t r , for general function f ρ , I first set f ρ = B ( r ) 1 + e ( θ δ + ln ( r ) ) 2 δ r e c t mz , which is a Dirac type pulse as δ gets arbitrarily large. For the second type of density function of small amplitude, f ρ = B ( r ) θ 1 e 1 ln r 2 r e c t mz , where B ( r ) is a function of compact support on an interval [ 0 , r w ] . For the second choice of f ρ the density ρ is defined everywhere in tube when θ = ϵ is small for δ large and in addition f ρ solves the pde in δ and ρ above for both types.This will eliminate two of the nonlinear terms in the Navier Stokes equations below for δ sufficiently large and negative.
Using Equations (4)–(9),multiplying Equations (1)–(3) by cartesian unit vectors e r * = ( 1 , 0 , 0 ) , e θ * = ( 0 , 1 , 0 ) and k = ( 0 , 0 , 1 ) respectively and adding Equations (1)–(3) gives the following equation, for the resulting composite vector L 1 = 1 δ 2 u r * * e r * + 1 δ u θ * * e θ * + 1 δ u z * * k ,
δ L 1 t + u r * * δ 2 L 1 r * + u θ * * δ r * L 1 θ * + 1 δ u z * * L 1 z * 1 δ 2 r * u θ * * 2 e r * + 1 δ 3 r * u r * * u θ * * e θ * μ ρ δ 2 L 1 r * 2 + δ 2 2 L 1 r * 2 + δ 2 1 r * L 1 r * + δ 2 1 r * 2 2 L 1 θ * 2 + 2 δ 2 r * 2 1 δ 2 u r * * θ * e θ * 1 δ u θ * * θ * e r * + δ 2 2 L 1 z * 2 μ 3 ρ r * u r * * r * + 1 r * u r * * θ * + u r * * z * e r * δ r * μ 3 ρ θ * u θ * * r * + 1 r * u θ * * θ * + u θ * * z * e θ * δ μ 3 ρ z * u z * * r * + 1 r * u z * * θ * + u z * * z * k + δ 1 ρ p r * e r * + δ 1 r * ρ p θ * e θ * + δ 1 ρ p z * k F T = 0
Furthermore one can take the cross product of Equation (10) with e r * and e θ * to form two equations respectively where in the second case multiplication by δ to cancel squared nonlinear terms, 1 δ 2 r * u θ * * 2 e r * , 1 δ 3 r * u r * * u θ * * e θ * , leads to,
δ L A t + u r * * δ 2 L A r * + u θ * * δ r * L A θ * + 1 δ u z * * L A z * + δ 2 u θ * * r * ρ L ¯ ρ t u z * * ρ z * ρ r * k μ ρ δ 2 L A r * 2 + δ 2 2 L A r * 2 + δ 2 1 r * L A r * + δ 2 1 r * 2 2 L A θ * 2 + 2 δ 2 r * 2 1 δ 2 u r * * θ * k 1 δ u θ * * θ * k + δ 2 2 L A z * 2 + + δ r * μ 3 ρ θ * u θ * * r * + 1 r * u θ * * θ * + u θ * * z * k + δ μ 3 ρ z * u z * * r * + 1 r * u z * * θ * + u z * * z * e θ * δ μ 3 ρ r * u r * * r * + 1 r * u r * * θ * + u r * * z * k + δ 2 μ 3 ρ z * u z * * r * + 1 r * u z * * θ * + u z * * z * e r * δ 1 r * ρ p θ * k δ 1 ρ p z * e θ * + δ 2 1 ρ p r * k δ 2 1 ρ p z * e r * + F T = 0
where L A = u z * * e r * 1 δ u z * * e θ * + 1 δ ( u r * * u θ * * ) k , L A = L 1 × e r * + δ L 1 × e θ * and continuity equation, Equation (12) below has been used. Equation (10) is sufficient to carry out the subsequent analysis.

3. A Solution Procedure for δ Arbitrarily Large in Quantity

Multiplication of Equation (10) by ρ δ and Equation (12) below by L 1 δ , addition of the resulting equations [9], and using the ordinary product rule of differential multivariable calculus a form as in Equation (13) is obtained whereby I set a = ρ L 1 .
The continuity equation in cylindrical co-ordinates is
ρ t + 1 δ u r * * ρ r * + u θ * * r * ρ θ * + u z * * ρ z * = ρ 1 δ u r * * r * + 1 δ u r * * r * + 1 r * u θ * * θ * + u z * * z * = ρ L ¯
ρ a t + a · a + ρ 2 b · b = μ ( δ ) 2 b + 1 3 ( · b ) + P + F T
where μ ( δ ) = μ δ .
Taking the geometric product in the previous equation with the inertial vector term,
f = a · a
where b = a ρ is defined, where in the context of Geometric Algebra, the following scalar and vector grade equations arise,
f · ρ 2 b t + ρ b ρ t + f 2 + b · f ρ 2 · b = μ ( δ ) f · 2 b + μ ( δ ) 3 f · ( · b ) + f · P
ρ 2 b t + ρ b ρ t + b ρ 2 · b = μ ( δ ) 2 b + 1 3 ( · b ) + P + F T
The geometric product of two vectors [10], is defined by A B = A · B + A × B .
Taking the divergence of Equation (16) results in
ρ 2 t · b + ρ ρ t · b + b · ( ρ ρ t ) + · ( b ρ 2 · b ) = · μ ( δ ) 2 ( b + μ ( δ ) 3 · ( · b ) + · P + · F T
Upon multiplication of Equation (17) by,
H = ρ b · f ρ t
the resulting equation is
ρ 2 H t ( · b ) + ρ 2 b · f · b + H b · ( ρ ρ t ) + H · ( b ρ 2 · b ) = μ ( δ ) H · 2 b + μ ( δ ) 3 H · ( ( · b ) ) + H · ( P ) + H · F T
which results upon using Equation (15) in,
ρ 2 H t ( · b ) f · ρ 2 b t + ρ b ρ t f 2 μ | δ | f · 2 b μ | δ | 3 f · ( · b ) f · P + H b · ( ρ ρ t ) + H · ( b ρ 2 · b ) = μ | δ | H ( · 2 b ) μ | δ | 3 H · ( ( · b ) ) + H · P + H · F T
The continuity equation is written in terms of b as,
1 δ ρ t + ρ · b + b · ρ = 0
and
· b = 1 δ ρ ρ t ρ · b
where the following compact expression is given,
Y * = · b
Multiplying by ( f · f ) 1 in Equation (20) 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,
W * Y * t G ( ρ , ρ t ) W * F ( ρ , ρ t ) b f ( 1 + f · P ) ρ 2 V ( μ ) | δ | W * 2 ( Y * ) b f ρ t ρ 1 f 2 f · b t + Ω + b f ρ t ρ 3 f 2 μ | δ | f · 2 b + b f ρ t ρ 3 f 2 μ | δ | 3 f · ( · b ) ρ 2 b f 1 f 2 b · f · b ρ 2 · b + P + F T = 0
where Ω = ρ 2 H b · ( ρ ρ t ) in Equation (24) and,
W * = f · b f 2 f b = f · b f 2 f · b + f · b f 2 f × b = ξ + f × f · b f 2 b
This involves the vector projection of b onto f which is written in the conventional form,
proj f b = f · b f 2 f
Equation (24) can be written compactly as
Y * t G ( ρ , ρ t ) ρ 2 V ( μ ) | δ | 2 Y * ρ 2 2 P = ρ 2 U f Q ( ρ , ρ t , b , b t , P , F T ) + f U f b
where U f ξ is the scalar projection for b , G = 1 ρ 2 ρ t 2 , Q (a differential operator defined by Equations (16) and (24)) and hence for a constant α ,
Y * t G ( ρ , ρ t ) ρ 2 V ( μ ) | δ | 2 Y * ρ 2 2 P = ρ 2 Q ( ρ , ρ t , b , b t , P , F T ) + f b = α ρ 2 0
with solution in terms of a function B ,
Y * = · b = B ( α , r , θ , z , t )
If the rate of change of ρ with respect to t is changing exponentially in time as f ρ = N ( r ) e c t mz ( ρ = f ρ and c = 1 δ < 0 , m < 0 ) , where N ( r ) is either of two density type functions defined in the beginning of this paper, it may be proven that B is written in terms of a separable function, if ρ 2 2 P is chosen to be negligible. and I obtain,
Y r , θ , z , t = F 1 r F 2 θ F 3 z F 4 t 1 / 4 α δ V ( μ ( δ ) ) δ 3 C 2 C 1 r + C 2 r 2 + C 3 r
where F i are the separable parts of the function, C 1 , C 2 , C 3 , are arbitrary constants, and as C 2 ,
M ( r ) = 1 / 4 α δ V ( μ ( δ ) ) δ 3 C 2 C 1 r + C 2 r 2 + C 3 r = 1 / 4 α δ r 2 δ 3 V ( μ ( δ ) )
For c 2 approaching zero and using the properties of the norm,
Y * r , θ , z , t F 1 r F 2 θ F 3 z F 4 t b Q ρ , ρ t , b , b t , P , F T + f = 0
Furthermore using the property for a norm on a vector space of continuous functions (since away from r = 0 ), A = 0 iff A = 0 , and Equation (16) and the second line of Equation (24) gives,
Y * r , θ , z , t F 1 r F 2 θ F 3 z F 4 t b + 2 b t 1 ρ b ρ t F ( ρ , ρ t ) ( f + F T ) ρ 2 b · b ρ 2 · b + F T = 0
where Y * r , θ , z , t F 1 r F 2 θ F 3 z F 4 t 0 , for δ , with μ chosen such that V ( μ ( δ ) ) = α ( δ ) (here it can be proved that the viscosity is a function of density which is a Dirac pulse function, also in general it can be proven that Equation (28) with viscosity as a function of small amplitude type density function can have a similar form as in Equation (30)) for a new form of M ( r ) . Next continuing with Equation (33),
2 b t 1 ρ b ρ t ( c f b + F T ) ρ 2 b · F T = 0
where F ( ρ , ρ t ) = ρ 3 ρ t . Ω in Equation (24) vanishes due to assumption on rate of change of density with respect to t and c 2 0 , and finally Equation (22) has been used with a calculation done to show that the expression reduces in going from Equation (33) to Equation (34).
Also, F ( ρ , ρ t ) f = ρ 3 ρ t ρ 2 b · b = ρ 1 ρ t b · b = c f b .
For large r it can be proven that F 1 ( r ) is decaying as a J 0 Bessel function and F 0 is an increasing function radially in tube ( r > 0 ) with the following resulting upon taking the curl of Equation (34),
2 t × b c × f b × ρ 2 b · F T = 0
where the conservative gravity force drops out.
Multiply Equation (35) by the normal vector cos ( θ ) a which is the normal component of a at wall of moving control volume (CV) in Figure 1,
cos ( θ ) a · 2 t × b c × f b × ρ 2 b · F T = 0
Recalling Divergence theorem and Stoke’s theorem, for general F ,
V ( · F ) d V = S ( V ) F · n ^ d S
S ( × F ) · d S = C F ( r ) · d r
where C is the contour of a circle in control volume of tube and S consists of all surfaces of control volume.
Defining the following vector field,
W = t × b 1 2 × ρ 2 b · F T ,
S ( V ) W · n ^ d S = V ( · W ) d V
= c 2 C f b ( r ) · d r
where n ^ = cos ( θ ) a , and Stoke’s theorem has been used. Applying Stoke’s theorem to W , hence b .
c 2 C f b ( r ) · T d s = C ( b t + 1 2 ρ 2 b · F T ) · d r
where T is unit tangent vector to closed curve C and d s is arc length,
C ( c 2 f b ( r ) + b t + 1 2 ρ 2 b · F T ) · d r = 0
The third term in the parenthesis in Equation (42) is integrated by parts for line integral and I obtain the following,
C ( c 2 f b ( r ) + b t 1 2 F T · ( ρ 2 b ) ) · d r = 0
Parametrizing the circle as r = g ( θ ) in polar coordinates it can be proven that the line integral in Equation (43) is,
C ( c f b 1 2 + b 1 t 1 2 F T 1 · ( ρ 2 b ) ) · d θ C ( c f b 2 2 + b 2 t 1 2 F T 2 · ( ρ 2 b ) ) · d r = 0
The normal form of Green’s theorem can be used for the line integral in Equation (44), setting first,
M = c f b 1 2 + b 1 t 1 2 F T 1 · ( ρ 2 b )
N = c f b 2 2 + b 2 t 1 2 F T 2 · ( ρ 2 b )
The line integral in Equation (44) is equal to the following,
R M r + N θ d r d θ
where M and N are given by Equations (45) and (46) respectively and R is the open disk with boundary C. The gravity force F T 1 is expressed as follows,
F T 1 = ϕ = g h = g ( b T cos ( π 4 θ ) )
where ϕ is a potential function, h is the negative height in the direction of the vector b g in Figure 2, T is a time constant and g is gravity constant. Also cos ( π 4 θ ) b = b g at the vertical vector b g in Figure 2. It follows that Equation (47) becomes,
R 1 2 g T r 1 r 2 b g θ 2 ρ 2 + ( c r b g b g θ + b g t ) θ d r d θ
R 1 2 g T r b g θ 2 1 r 2 1 + e 1 ln r 2 r 2 + ( c r b g b g θ + b g t ) θ d r d θ
R 1 2 g T b g θ 2 4 e ln r 1 2 ln r 4 e ln r 1 2 2 r r r + e ln r 1 2 3 + ( c r b g b g θ + b g t ) θ d r d θ
R 1 2 g T 4 e ln r 1 2 ln r 4 e ln r 1 2 2 r + θ r r + e ln r 1 2 3 c r b g b g θ + r r + e ln r 1 2 3 b g t d r d θ
For δ the function b 2 = b g ( θ , t ) and b 1 = b 1 ( r , t ) 0 , and thus for c < 0 small, integration in r on the interval [ 0 , r w ] , where r w is near the wall in control volume,
0 2 π W 2 1 2 g T c 1 b g θ 2 + θ W 1 b g b g θ + r c 1 W 3 b g t ) θ d r d θ
In the above analysis I have integrated the second term in second line of Equation (52) wiith respect to r and have obtained,
W 1 ( r ) = 3 / 2 A π erf 5 / 2 + ln r + 3 / 4 B π 2 erf 1 / 2 2 ln r 3 2 + 1 / 6 C 3 π erf 1 / 6 6 ln r 7 3 + 1 / 4 r 4 + C
where A, B and C are constants. Letting r = r w δ and δ , I obtain a general non-zero constant W 1 = 3 A ( π ) ( 1 / 2 ) 3 B ( π ) ( 2 ) ( 1 / 4 ) ( 1 / 6 ) C ( 3 ) ( π ) + C . Integrating the first term in Equation (52) with respect to r following the same procedure as in previous step, the result is that I obtain a general nonzero constant W 2 . The same is true for the third term in Equation (52) where I have a non-zero constant W 3 .
Choosing constants in W 1 , W 2 and W 3 such that W 2 1 2 g T c 1 = 1 / 2 , W 1 = 1 and c 1 W 3 = 1, the resulting equation obtained is the Hunter-Saxton PDE.

4. The Hunter-Saxton Equation

Substituting f 2 and F T 1 into Equations (45) and (46) since F T 2 = 0 and f 1 , b 1 terms are negligible in the boundary layer(Recall u r = u r * / δ 2 then one can obtain an integral form of the Hunter-Saxton equation using Equation (53),
b 2 t * + b 2 * b 2 θ * θ = 1 2 b 2 θ * 2
It is of interest that a more complicated partial differential equation is obtained upon taking F T 2 to be non zero and thus defining a rotational force related to the vorticity of the fluid elements in the boundary layer. In the interval [ 0 , t ] the boundary layer starts to form at t = 0 and reaches maximum height at time t. See Figure 1 and Figure 3 for control volume over the growth of the boundary layer and contours between boundary layer and tube wall. The time dependence is shown in the inertial term f . The right side of Equation (40) consists of nonlinear inertial term f and Equation (39) shows the dependence on these gradients and rate of change of curl of b with respect to t.

5. Discussion

It is known that there exist sufficient conditions for the blowup in time of smooth solutions to the compressible Navier-Stokes equations in arbitrary space dimensions with initial density of compact support [11]. The present work focuses on the compressible Navier Stokes equations in 3D in cylindrical coordinates. It has been shown that this system is reducible to the Hunter-Saxton equation which is used to model nematic liquid crystals and this nonlinear PDE also blows up in finite time if the initial boundary condition is set equal to a non-monotone increasing function [12].

6. Conclusions

An attempt has been made to reduce the compressible Navier Stokes equations coupled to the continuity equation in cylindrical co-ordinates to a simpler problem in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the tube. By seeking an additive solution using the continuity equation and the simplified vector equation obtained by a similar procedure, using the product rule of differentiable calculus, and using Geometric Algebra as a starting point, it follows through analysis that the integral of total divergence of a specific vector field over time can be expressed as the integral with respect to time of the line integral of the dot product of inertial and azimuthal velocity. The line integral can be evaluated on a contour that is annular and traces the boundary layer as time increases in the flow. It has been shown that a reduction of the 3D compressible unsteady Navier-Stokes equations to a single partial differential equation is possible and integral calculus methods are applied for the case of a body force in the direction of gravity to obtain an integral form of the Hunter-Saxton equation. The radial flow due to the gravitational pull at the centre of the disk in Figure 2 does not result in an immediate radial flow downwards towards the wall of the tube. Rather, the fluid compresses near the centre of the disk where the gravitational force exists (and is greatest there due to a nonuniform mass distribution where it is assumed that the mass is greatest at the center of disk) and hence the density increases locally in response to the force. The compressed fluid expands against the neighbouring fluid particles causing the neighboring fluid particles to compress and set in motion a wave pulse that travels downwards towards the wall in Figure 2 [13]. It is worthy to note that for this reason the density must be not-monotone increasing, i.e., it must either be a Dirac pulse or Gaussian-like wave function (i.e., of compact support), or a pulse of small amplitude.The reduction obtained in this paper allows one to analyze the Hunter-Saxton PDE which gives great insight into the nature of nematic liquid crystals which can be formed due to large temporal and spatial gradients in density. The viscosity can thus be used in describing the flow with boundary layer formation near the surface of the tube. Both a dynamic and rotational viscosity are related to dissipation of energy and on reorientation of the NLC director. The Hunter-Saxton system is known to blow up in finite time when the function used as an initial condition is a not-monotone increasing function. Since the density is necessarily to be either an infinite Dirac function or a pulse of finite and small amplitude, it is concluded that the governing system of compressible non-linear 3D Navier Stokes equations coupled to the space-time density dependent continuity equation will blow up in finite time. A restriction of F T 2 = 0 was made and one can consider a more general body force than has been shown where in addition there is a rotational force applied. Future studies are required to examine this problem.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Referring to Figure 4, let,
d = k a
d = m c
It follows that,
d + d 2 = d = k a + m c
If a and c are in the same direction then k m > 0 and k , m are either both positive or both negative. Also according to Figure 4 it follows that,
d 2 = c 2 2 + d 2 2
d 2 = a 1 2 + d 1 2
If the circle shrinks to an arbitrary small circle in diameter then
d 2 d 2 2
d 2 d 1 2
Summing up all the d i components and letting b 2 i e θ represent the tangential velocity for each point on the circumference of disk R, it follows for N distinct points that,
i = 1 N d i = i = 1 N b 2 i
Thus in the limit as the circle diameter gets small, introduce a function g defined as,
g = b 2 θ
where b = b 2 e θ and b θ is the vector pointing towards the center of the disk in R.

References

  1. Pereira, L.M.; Guerrero, J.S.P.; Cotta, R.M. Integral transformation of the Navier Stokes equations in cylindrical geometry. Comput. Mech. 1998, 21, 60–70. [Google Scholar] [CrossRef]
  2. Taylor, C.; Hughes, T.G. Finite Element Programming of the Navier Stokes Equations; Pineridge Press: Swansea, UK, 1981; Available online: https://trove.nla.gov.au/work/25314368?q&versionId=30518049 (accessed on 24 January 2018).
  3. MacCormack, R.W. A Numerical Method for Solving the Equations of Compressible Viscous Flow. AIAA J. 1982, 20, 1275–1281. [Google Scholar] [CrossRef]
  4. Vadasz, P. Rendering the Navier-Stokes equations for a compressible fluid into the Schrödinger Equation for quantum mechanics. Fluids 2016, 1, 18. [Google Scholar] [CrossRef]
  5. Feistauer, M.; Felcman, J.; Straskraba, I. Mathematical and Computational Methods for Compressible Flow; Oxford University Press: Oxford, UK, 2004. [Google Scholar]
  6. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge University Press: Cambridge, UK, 1967. [Google Scholar]
  7. Ebin, D.G. Motion of a slightly compressible fluid. Proc. Natl. Acad. Sci. USA 1975, 72, 539–542. [Google Scholar] [CrossRef] [PubMed]
  8. Cai, R. Some Explicit analytical solutions of unsteady compressible flow. J. Fluids Eng. 1998, 120, 760–764. [Google Scholar] [CrossRef]
  9. Moschandreou, T. A New Analytical Procedure to Solve Two Phase Flow in Tubes. Math. Comput. Appl. 2018, 23, 26. [Google Scholar] [CrossRef]
  10. Doran, C.; Lasenby, A. Geometric Algebra for Physicists; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  11. Xin, Z. Blowup of Smooth Solutions to the Compressible Navier-Stokes Equation with Compact Density. Commun. Pure Appl. Math. 1998, 51, 229–240. [Google Scholar] [CrossRef]
  12. Hunter, J.K.; Saxton, R. Dynamics of Director Fields. SIAM J. Appl. Math. 1991, 51, 1498–1521. [Google Scholar] [CrossRef]
  13. Morisson, F.A. Compressible Fluids; Michigan Technological University: Houghton, MI, USA, 2004; pp. 94–98. [Google Scholar]
Figure 1. Compressible Viscous Flow in a Tube, z b l is the distance to achieve maximum boundary layer height.
Figure 1. Compressible Viscous Flow in a Tube, z b l is the distance to achieve maximum boundary layer height.
Mathematics 07 00126 g001
Figure 2. Vector b g which is pointing in the direction of increasing gravitational force. The right angle triangles show vector addition as expressed by solution, L 1 = b = 1 δ 2 u r * * e r * + 1 δ u θ * * e θ * .
Figure 2. Vector b g which is pointing in the direction of increasing gravitational force. The right angle triangles show vector addition as expressed by solution, L 1 = b = 1 δ 2 u r * * e r * + 1 δ u θ * * e θ * .
Mathematics 07 00126 g002
Figure 3. A typical set of contour lines between edge of boundary layer and tube wall.
Figure 3. A typical set of contour lines between edge of boundary layer and tube wall.
Mathematics 07 00126 g003
Figure 4. Vectors a , c , d as well as d 1 and d 2 which are pointing in the direction of increasing gravitational force. (Appendix A).
Figure 4. Vectors a , c , d as well as d 1 and d 2 which are pointing in the direction of increasing gravitational force. (Appendix A).
Mathematics 07 00126 g004

Share and Cite

MDPI and ACS Style

Moschandreou, T.E. A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra. Mathematics 2019, 7, 126. https://doi.org/10.3390/math7020126

AMA Style

Moschandreou TE. A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra. Mathematics. 2019; 7(2):126. https://doi.org/10.3390/math7020126

Chicago/Turabian Style

Moschandreou, Terry E. 2019. "A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra" Mathematics 7, no. 2: 126. https://doi.org/10.3390/math7020126

APA Style

Moschandreou, T. E. (2019). A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra. Mathematics, 7(2), 126. https://doi.org/10.3390/math7020126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop