## Flow problem in three-dimensional geometry

Maxim Zaytsev 1) , Vyacheslav Akkerman 2)

1) Nuclear Safety Institute of Russian Academy of Sciences

115191 Moscow, ul. Bolshaya Tulskaya, 52, Russia

2) West Virginia University Morgantown, WV 26506-6106, USA

## Abstract

It  is  shown  how  a  complete  set  of  hydrodynamic  equations  describing  an  unsteady  threedimensional viscous flow nearby a solid body, can be reduced to a closed system of surface equations using the method of dimension reduction of over-determined systems of differential equations.  These  systems  of  equations  allow  determining  the  surface  distribution  of  the resulting stresses on the surface of this body as well as all other quantities characterizing the hydrodynamic flow around it.

Keywords :  flow problem, hydrodynamic flow, viscosity, solids, stress on the body surface, differential equation on the surface

## Introduction

In  application,  nonlinear  equations  of  hydrodynamics  are  often  one  of  the  obstacles involved in the study of many phenomenas, such as turbulence, hydrodynamic discontinuities motion,  movement of the fronts of chemical reactions, etc. [1-4]. (During their modeling in practical problems often one has to produce a very fine numerical grid in space and time, which requires  large  computing  power  and  time-consuming)  When  modeling  such  equations  in practical problem, one has to often produce a very fine numerical grid in space and time which is time consumming and requires large computing power hence rendering the methode as not practical [5]. In this regard, various ways of reducing the complete system of hydrodynamic equations in terms of a system of surface equations are of a large scientific interest. [6-11]. Such a procedure can reduce the dimension of the problem by one (3 2 D D  , 2 1 D D  ),  which significantly  reduces  the  required  computational  power.  In  particular,  the  corresponding computer program could directly calculate the hydrodynamic discontinuities with viscosity, sound formation and other changes in the density of gases and liquids (using the information only on the surface). For example, the problem of describing the potential flow in a plane can be reduced to an integral equation on the boundary for a vanishingly small viscosity (Dirichlet, Neumann) [1, 4]. This equation relates the tangential and normal components of the velocity. If one of the components is known on the boundary of the body, one can define all the external flow.

In this paper we reduce the equations of a viscous medium "in the volume" to surface equations and demonstrate how the proposed method can reduce the dimension of the problems in viscous flow, which makes the equations easier to simplify the calculations in a variety of applications. General methods proposed for overdetemination of the Navier-Stokes equations in three-dimensional and two-dimensional case, which theoretically allows reducing them to a closed system of equations on any surface [12, 13]. No significant simplifying assumptions, which would limit the generality of the consideration of these phenomena, are considered.

## 1. Formulation of problem

Movement of solids in viscous media is characterized by the conditions on the boundary of a solid, where the particles stick to it. In other words, the rate of viscous medium on the surface of the body u is determined by the kinematic characteristics of the surface, for example,

it  can be expressed in terms of the speed of the body surface V .  Force F acting on a unit surface area is [1]

$$F _ { i } = P n _ { i } - \sigma _ { i k } ^ { / } n _ { k }$$

where n is normal to the surface, and

$$\sigma _ { i k } ^ { \prime } = \eta \left ( \frac { \partial u _ { i } } { \partial x _ { k } } + \frac { \partial u _ { k } } { \partial x _ { i } } \right )$$

the  viscous  stress  tensor  and  -  coefficient  of  dynamic  viscosity.  Therefore,  according  to equation (1), in this case, for any point on the body surface in a Cartesian coordinate system   2 1 , , τ τ n one has

$$F _ { n } = P + 2 \eta \left ( \frac { \partial u _ { 1 } } { \partial \tau _ { 1 } } + \frac { \partial u _ { 2 } } { \partial \tau _ { 2 } } \right ),$$

$$F _ { 1 } = \eta \left ( \frac { \partial u _ { 1 } } { \partial n } + \frac { \partial u _ { n } } { \partial \tau _ { 1 } } \right ) = \eta \left ( \omega _ { 2 } + 2 \frac { \partial u _ { n } } { \partial \tau _ { 1 } } \right ),$$

$$F _ { _ { 2 } } = \eta \left ( \frac { \partial u _ { _ { 2 } } } { \partial n } + \frac { \partial u _ { _ { n } } } { \partial \tau _ { _ { 2 } } } \right ) = \eta \left ( - \omega _ { _ { 1 } } + 2 \frac { \partial u _ { _ { n } } } { \partial \tau _ { _ { 2 } } } \right )$$

(see Fig. 1). In addition, the force F relates to the velocity of the body V via the laws of rigid body motion dynamics. Unsteady flow around fixed solids in viscous medium is different from the above by cancelling velocity conditions at the boundary [1]:

$$V = 0 \, \mathbf u \, \mathbf u = 0 \,.$$

## 2. Overdetermination of Navier-Stokes equations

We confine ourselves here to an incompressible viscous flow in three-dimensional space in  the  gravity  field g .  To  illustrate  this,  the  Navier-Stokes  equations  take  the  form  [1]

$$\frac { \partial \mathbf u } { \partial t } - \mathbf u \times \mathfrak { o } + \nabla \left ( P + \frac { 1 } { 2 } \mathbf u ^ { 2 } \right ) = - \nu \nabla \times \mathfrak { o } + g \,,$$

$$\omega = \nabla \times \mathfrak { u } \,,$$

div

0,



u

(9)

where  - is the kinematic viscosity coefficient. These equations allow us to find the unknown velocity vector   t , r u and the pressure profile   t P , r at any given point and time in space, i.e. for any   t , r . We introduce a "correction" α to the vector u and 'pseudo-density'   t , r  in order to generalize the system of equations (7) - (9) as follows

$$\frac { \partial ( \mathfrak { u } + \mathfrak { a } ) } { \partial t } - \left [ ( \mathfrak { u } + \mathfrak { a } ) \times ( \mathfrak { w } + \nabla \times \mathfrak { a } ) \right ] + \nabla \left ( P + \frac { 1 } { 2 } \mathfrak { u } ^ { 2 } \right ) = \mathfrak { g } \,,$$

$$\frac { \partial \rho } { \partial t } + ( \mathfrak { u } + \mathfrak { a } ) \cdot \nabla \rho + \rho \cdot \text{div} ( \mathfrak { u } + \mathfrak { a } ) = 0 \,,$$

$$\frac { \partial \alpha } { \partial t } = [ \alpha \times \omega + \mathfrak { u } \times ( \nabla \times \mathfrak { a } ) + \mathfrak { a } \times ( \nabla \times \mathfrak { a } ) ] + \nu \nabla \times \omega$$

$$\omega = \nabla \times \mathfrak { u } \,,$$

$$\text{div} \, = \, 0 \,.$$

Actually, the equations (11) and (12) respectively determine   t , r  and α but the equation (10) is resulting from Equations (7) and (12). Vector field α u  can be formally considered as inviscid hydrodynamic flow with density   t , r  and pressure   2 2 1 1 2 2 P          u u α .

We go one to one to the pseudo Lagrangian variables of initial position of the gas particles and time  ( )  (i.e. marking) , 0 t r

d

r

dt





(

,

)

(

,

)

(

,

)

(

,

)

t

t

t

t











u

r

α

r

ω

r

α

r

,

$$\mathbf r = \mathbf r ( \mathbf r _ { 0 }, t ) \quad \mathbf r \quad \mathbf r _ { 0 } = \mathbf r _ { 0 } ( \mathbf r, t ), \ \mathbf r _ { 0 } \Big | _ { t = 0 } = \mathbf r \,.$$

We chose ( , )  so that t  r

$$\lambda ( \mathfrak { w } + \nabla \times \mathfrak { a } ) \cdot \nabla \rho + \rho ( \mathfrak { w } + \nabla \times \mathfrak { a } ) \cdot \nabla \lambda = 0 \,.$$

Hence, it can be shown that

$$\frac { \partial \mathbf r _ { 0 } } { \partial t } + \left ( \mathbf u + \mathbf a + \mathcal { A } \left ( \mathfrak w + \nabla \times \mathbf a \right ) \right ) \cdot \nabla \mathbf r _ { 0 } = 0 \,.$$

Then, according to appendix A, from the equations (10) and (11) one yields the following expressions:

$$\frac { ( \omega + \nabla \times \mathfrak { a } ) \cdot \nabla r _ { 0 } } { \rho } = \frac { ( \omega _ { 0 } ( \mathbf r _ { 0 } ) + \nabla \times \mathfrak { a } _ { 0 } ( \mathbf r _ { 0 } ) ) } { \rho _ { 0 } ( \mathbf r _ { 0 } ) },$$

$$\frac { 1 } { \rho } \frac { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } { \partial ( x, y, z ) } = \frac { 1 } { \rho _ { 0 } ( \mathbf r _ { 0 } ) } \,,$$

where 0 α , 0  , 0 0 u ω    -  initial  distribution  of  the  variables α ,  , ω .  Let 0 0 and  α 0 1. Then (18) and (19) can be transformed to the form [14]:  

$$\omega _ { x } + \nabla \times \alpha _ { x } - \left | \frac { \partial \mathbf r _ { 0 } } { \partial y } \left ( \frac { \partial \mathbf r _ { 0 } } { \partial z } \right ) \left ( \mathbf r _ { 0 } \right ) \right | = 0 \,,$$

$$\omega _ { y } + \nabla \times \alpha _ { y } - \left | \frac { \partial \mathbf r _ { 0 } } { \partial z } \left ( \frac { \partial \mathbf r _ { 0 } } { \partial x } \right ) \left ( \mathbf r _ { 0 } \right ) \right | = 0 \,,$$

$$\omega _ { z } + \nabla \times \alpha _ { z } - \left | \frac { \partial r _ { 0 } } { \partial x } \left ( \frac { \partial r _ { 0 } } { \partial y } \right ) = 0 \,.$$

Actually, from three equations (20) - (22) only two are independent [14]. Formula (19) takes the form

$$\frac { 1 } { \rho } \frac { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } { \partial ( x, y, z ) } = 1.$$

Thus, we have obtained an overdetermined system of 16 differential equations (7) - (9) (11) (12) (16) (20) (21) (23) plus any of the equations from (17) and of 15 unknowns α , u , ω , P ,  ,  и 0 r . Causation analysis shows that the system is consistent and valid [12-14]. It is possible to add to the equation (12) a nontrivial external force.

## 3. Viscous flow in a three-dimensional case

We derive a system of equations describing the unsteady flow around the fixed rigid body in three-dimensional incompressible viscous stream. Let us denote

$$& S _ { 1 } = \alpha _ { x } \,, \ S _ { 2 } = \alpha _ { y } \,, \ S _ { 3 } = \alpha _ { z } \,, \ S _ { 4 } = u _ { x } \,, \ S _ { 5 } = u _ { y } \,, \ S _ { 6 } = u _ { z } \,, \ S _ { 7 } = \omega _ { x } \,, \ S _ { 8 } = \omega _ { y } \,, \ S _ { 9 } = \omega _ { z } \,, \\ & S _ { 1 0 } = P \,, \ S _ { 1 1 } = \rho \,, \ S _ { 1 2 } = \lambda \,, \ S _ { 1 3 } = x _ { 0 } \,, \ S _ { 1 4 } = y _ { 0 } \,, \ S _ { 1 5 } = z _ { 0 } \,.$$

Then,  a  system  of  16  equations  (7)  -  (9),  (11),  (12),  (16),  (20),  (21),  (23)  plus  any  of  the equations (17) can be written in form

$$H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial r }, \frac { \partial S _ { v } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 1 5 \,, \, k = 1 \dots 1 6 \,.$$

We now turn to the coordinate system   1 2 , , τ τ n at point M (see Fig. 1). Then equations (25) can be written as

$$H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial \tau _ { 1 } }, \frac { \partial S _ { v } } { \partial \tau _ { 2 } }, \frac { \partial S _ { v } } { \partial n }, \frac { \partial S _ { v } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 1 5 \,, \, k = 1 \dots 1 6 \,.$$

We differentiate equations (26) in the direction n 14 times. Then we obtain the 240 equations of the form

$$\left [ H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial \tau _ { 1 } }, \frac { \partial S _ { v } } { \partial \tau _ { 2 } }, \frac { \partial S _ { v } } { \partial n } \frac { \partial S _ { v } } { \partial t } \dots \right ) \right ] ^ { ( i ) } _ { n } = 0 \,, \, v = 1 \dots 1 5 \,, \, k = 1 \dots 1 6 \,, \, i = 0 \dots 1 4 \,.$$

We denote

$$S _ { v } ^ { j } = \partial ^ { j } S _ { v } / \partial n ^ { j }$$

where 0...15 , j  0 v v S S  . Then 240 equations (27) can be written as:

$$O _ { _ { h } } \left ( S _ { _ { v } } ^ { j }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial \tau _ { 1 } }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial \tau _ { 2 } }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 1 5 \,, \, h = 1 \dots 2 4 0 \,, \, j = 0 \dots 1 5 \,.$$

We extend the variable 0 12 12 S S    along the boundary of the body as follow:

$$( \alpha \cdot \mathbf n ) + \lambda ( \omega + \nabla \times \alpha ) \cdot \mathbf n = 0 \,.$$

Let 0  be the equation of the solid surface in the coordinate system ( ) , ( 0 0  t g r )  (15). Then , 0 t r from (6) and (30) its speed is equal to

$$\cdot = \frac { \dot { g } _ { t } + \left ( \left ( \mathbf u + \mathbf \alpha + \lambda \left ( \mathbf \omega + \nabla \times \mathbf \alpha \right ) \right ) \cdot \nabla g \right ) } { } =$$

$$\dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ & V _ { 0 } = \frac { g _ { 0 t } } { | \nabla _ { 0 } g _ { 0 } | } = \frac { g _ { t } + \left ( \left ( u + \alpha + \lambda \left ( \omega + \nabla \times \alpha \right ) \right ) \cdot \nabla g \right ) } { | \nabla _ { 0 } g _ { 0 } | } = \\ & = \frac { | \nabla g | } { | \nabla _ { 0 } g _ { 0 } | } \left ( V + u _ { n } + \left ( \alpha \cdot \mathbf n \right ) + \lambda \left ( \omega + \nabla \times \alpha \right ) \cdot \mathbf n \right ) = 0. \\ - \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \collections$$

In the coordinate system  ( ) , 0 t r the solid body is not moving.

Thus, the dependence of the derivative to n in  the  transformed system of differential equations (29) is missing, and we obtain a closed system of 240 surface differential equations strictly along the boundary of a solid and a corresponding number of variables. In fact, some of the  variables  are  determined  from  (6)  and  (30).  At  the  initial  time,  it  is  required  to  know nontrivial vorticity   0 0 r ω and some of its spatial derivatives on the surface (formula (20), (21)). It is always possible to find where the vorticity is non-trivial, because there is a viscous boundary layer on the surface.

The advantage is that the system actually describes the process of flow past the body surface  in  terms  of  the  surface  itself.  Calculations  show  [15,  16]  that  the  expressions  (29) obtained in an analytical form are extremely cumbersome. Their writing out is beyond the scope of this paper. However, one provides a detailed algorithm of their obtaining.

## 4. Overdetermination of the complete system of hydrodynamic equations

Complete  system  of  equations  describing  the  hydrodynamic  environment  in  a  nonuniform unsteady 3D case is [1]

$$\frac { \partial \rho } { \partial t } + \text{div} ( \rho \mathbf u ) = 0 \,,$$

$$\frac { \partial \mathbf u } { \partial t } - \mathbf u \times \mathbf v = - \nabla \left ( \frac { \mathbf u ^ { 2 } } { 2 } \right ) - \frac { \nabla P } { \rho } + \frac { \nabla \sigma } { \rho } + \frac { \mathbf F } { \rho } \,,$$

$$\omega = \nabla \times \mathfrak { u } \,,$$

$$\psi = & \text{div} \,,$$

$$\rho T \left [ \frac { \partial s } { \partial t } + ( \mathbf u \nabla ) s \right ] = Q - \text{div} q + \Phi \,,$$

$$T = T ( \rho, s ), \ P = P ( \rho, s ),$$

$$\nabla w = T \cdot \nabla s + \frac { \nabla P } { \rho },$$

where  , u , P , T , w и s - are density, velocity, pressure, temperature, specific heat function, and the entropy (per unit mass) of the medium, respectively. Here F и Q -  are volumetric force and heat generation respectively, T    q -  is  the  heat  flow  equation  (  -  thermal conductivity),  - the viscous stress tensor, for which we have

$$\nabla \sigma = \mu \left ( \Delta \mathfrak { u } + \frac { 1 } { 3 } \nabla ( \nabla \mathfrak { u } ) \right ) = \mu \left ( - \nabla \times \mathfrak { w } + \frac { 4 } { 3 } \nabla \psi \right ).$$

Here is  - shear dynamic viscosity. The dissipation function  in (36) is equal to

$$\Phi = \frac { \mu } { 2 } \left ( \frac { \partial u _ { i } } { \partial x _ { k } } + \frac { \partial u _ { k } } { \partial x _ { i } } \right ) ^ { 2 } - \frac { 2 } { 3 } \, \mu \left ( \frac { \partial u _ { i } } { \partial x _ { i } } \right ) ^ { 2 }.$$

Equations (37) and (38) -are thermodynamic relations that depend on the properties of the medium.

We introduce a "correction"   t , r α to the vector )  and 'pseudo-density' , ( t r u   t , r   , which we determine from equations:

$$\frac { \hat { \partial } \rho ^ { \bullet } } { \hat { \partial } t } + ( \mathfrak { u } + \mathfrak { a } ) \cdot \nabla \rho ^ { \bullet } + \rho ^ { \bullet } \text{div} ( \mathfrak { u } + \mathfrak { a } ) = 0 \,,$$

$$\frac { \partial \alpha } { \partial t } = [ \alpha \times \omega + \mathfrak { u } \times ( \nabla \times \mathfrak { a } ) + \mathfrak { a } \times ( \nabla \times \mathfrak { a } ) ] - T \nabla s - \frac { \nabla \sigma } { \rho } - \frac { \mathrm F } { \rho } \,.$$

Termwise folded formula (42) and (33) and taking the operation «   »   on both sides, we find, using (38), that equation (33) is equivalent to the following expressions

$$\lambda & \lambda _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots \Delta _ { 0 } ^ { 2 } \cdots.$$

$$\frac { \partial ( \omega + \nabla \times \mathfrak { a } ) } { \partial t } = \nabla \times [ ( \mathfrak { u } + \mathfrak { a } ) \times ( \omega + \nabla \times \mathfrak { a } ) ].$$

We go one to one to the pseudo Lagrangian variables of initial position of the gas particles and time  ( )  (i.e. marking) , t r

0

$$\frac { d { \mathbf r } } { d t } = { \mathbf u } \left ( { \mathbf r }, t \right ) + { \mathbf a } \left ( { \mathbf r }, t \right ) + \lambda \left ( { \mathbf v } \left ( { \mathbf r }, t \right ) + \nabla \times { \mathbf a } \left ( { \mathbf r }, t \right ) \right ),$$

$$= \mathbf r ( \mathbf r _ { 0 }, \mathbf l ) \ \mathbf r _ { 0 } = \mathbf r _ { 0 } ( \mathbf r, \mathbf l ), \ \mathbf r _ { 0 } | _ { \mathbf l = 0 } = \mathbf r _ { 0 }.$$

$$\frac { \dots } { d t } = \mathfrak { u } \left ( \mathbf r, t \right ) + \mathbf a \left ( \mathbf r, t \right ) + \lambda \left ( \mathfrak { o } \left ( \mathbf r, t \right ) + \nabla \times \mathbf a \left ( \mathbf r, t \right ) \right ), \\ \mathbf r = \mathbf r \left ( \mathbf r _ { 0 }, t \right ) \quad \mathfrak { u } \quad \mathbf r _ { 0 } = \mathbf r _ { 0 } \left ( \mathbf r, t \right ), \ \mathbf r _ { 0 } \Big | _ { t = 0 } = \mathbf r \,.$$

We chose ( , )  so that t  r

$$\lambda ( \tilde { \omega } + \nabla \times \tilde { \alpha } \ ) \cdot \nabla \rho ^ { \bullet } + \rho ^ { \bullet } ( \tilde { \omega } + \nabla \times \tilde { \alpha } \ ) \cdot \nabla \lambda = 0.$$

Hence, it can be shown that

$$\frac { \partial r _ { 0 } } { \partial t } + \left ( \mathfrak { u } + \mathfrak { a } + \lambda \left ( \mathfrak { v } + \nabla \times \mathfrak { a } \ \right ) \right ) \cdot \nabla r _ { 0 } = 0 \,.$$

Equations (41) and (44) can be transformed to the form (see Appendix A)

$$\frac { ( \omega + \nabla \times \alpha ) \cdot \nabla r _ { 0 } } { \rho ^ { \bullet } } = \frac { ( \omega _ { 0 } ( r _ { 0 } ) + \nabla \times \alpha _ { 0 } ( r _ { 0 } ) ) } { \rho ^ { \bullet } _ { 0 } ( r _ { 0 } ) },$$

$$\frac { 1 } { \rho ^ { \bullet } } \frac { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } { \partial ( x, y, z ) } = \frac { 1 } { \rho ^ { \bullet } _ { 0 } ( \mathbf r _ { 0 } ) } \,,$$

where 0 α , 0   , 0 0 u ω    - initial distribution of the variables α ,   , ω . Let 0 0 and  α 0 0     . Then it can be shown [14], equations (48) - (49) can be transformed to

$$\omega _ { x } + \nabla \times \alpha _ { x } - \left | \frac { \partial \mathbf r _ { 0 } } { \partial y } \left ( \frac { \partial \mathbf r _ { 0 } } { \partial z } \right ) \left ( \mathbf r _ { 0 } \right ) \right | = 0 \,,$$

$$\omega _ { y } + \nabla \times \alpha _ { y } - \left | \frac { \left | \partial \mathbf r _ { 0 } } { \partial z } \left ( \frac { \partial \mathbf r _ { 0 } } { \partial x } \right ) \left | = 0 \,,$$

$$\omega _ { z } + \nabla \times \alpha _ { z } - \left | \frac { \partial r _ { 0 } } { \partial x } \left ( \frac { \partial r _ { 0 } } { \partial y } \right ) = 0 \,.$$

and the formula (49) takes the form

$$\frac { 1 } { \rho ^ { \bullet } } \frac { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } { \partial ( x, y, z ) } = \frac { 1 } { \rho _ { 0 } \left ( { \mathbf r } _ { 0 } \right ) }.$$

As a result, we get the following overdetermined system of 21 equations (32) - (38) (41) - (43), (46), (50), (51), (53), plus any of the equations from (47) and of 20 unknowns  , *  ,

$$\mathfrak { u } \,, \, \omega \,, \, \mathfrak { a } \,, \, P \,, \, T \,, \, s \,, \, w \,, \, \psi \,, \, \lambda \,, \, \mathfrak { r } _ { 0 }.$$

## 5. Three-dimensional general stream

We derive a system of equations describing the unsteady flow around the fixed rigid body in three-dimensional stream in the general case. Let

1 x S   , 2 y S   , 3 z S   , 4 x S u  , 5 y S u  , 6 z S u  , 7 x S   , 8 y S   , 9 z S   , 10 S P  , 11 S   , 12 S   , 13 0 S x  , 14 0 S y  , 15 0 S z  , 16 S    , 17 S T  , 18 S s  , 19 S   , 20 S w  , 21 x S q  , 22 y S q  , 23 z S q  . (54)

Then the system of 24 partial differential equations of the first order (32) - (38) (41) (43), (46), (50), (51), (53), plus any of the equations from (47) can be written as

$$H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial r }, \frac { \partial S _ { v } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 2 3 \,, \, k = 1 \dots 2 4 \,.$$

We now turn to the coordinate system   1 2 , , τ τ n at point M (see Fig. 1). Then equations (55) can be written as

$$H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial \tau _ { 1 } }, \frac { \partial S _ { v } } { \partial \tau _ { 2 } }, \frac { \partial S _ { v } } { \partial n }, \frac { \partial S _ { v } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 2 3 \,, \, k = 1 \dots 2 4 \,.$$

Differentiating equations (56) in the direction n 22 times, we obtain 552 equations of the form

$$\left [ H _ { k } \left ( S _ { v }, \frac { \partial S _ { v } } { \partial \tau _ { 1 } }, \frac { \partial S _ { v } } { \partial \tau _ { 2 } }, \frac { \partial S _ { v } } { \partial n } \frac { \partial S _ { v } } { \partial t } \dots \right ) \right ] _ { n } ^ { ( i ) } = 0 \,, \, v = 1 \dots 2 3 \,, \, k = 1 \dots 2 4 \,, \, i = 0 \dots 2 2 \,.$$

Denote

$$S _ { v } ^ { j } = \partial ^ { j } S _ { v } / \partial n ^ { j }$$

where 0...23, j  0 v v S S  . Then 552 equations (57) can be written as:

$$O _ { _ { h } } \left ( S _ { _ { v } } ^ { j }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial \tau _ { 1 } }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial \tau _ { 2 } }, \frac { \partial S _ { _ { v } } ^ { j } } { \partial t } \dots \right ) = 0 \,, \, v = 1 \dots 2 3 \,, \, h = 1 \dots 5 5 2 \,, \, j = 0 \dots 2 3 \,.$$

We extend the variable 0 12 12 S S    along the boundary of the body as follows:

$$( \alpha \cdot \mathbf n ) + \lambda ( \omega + \nabla \times \alpha ) \cdot \mathbf n = 0 \,.$$

Let 0 ) , ( 0 0  t g r be the equation of the solid surface in the coordinate system  ( )   (45). , 0 t r Similarly as in section 3, in the coordinate system  ( )  , the solid is static. , 0 t r

Similarly,  dependence  on  derivative  to n in  the  transformed  system  of  differential equations (59) is missing, and we obtain a closed system of 552 surface differential equations strictly along the boundary of a solid and the same number of variables (58). In fact, some

variables were already defined by (6) and (60), as well as by conditions of the heat transfer on surface of the body.

## Conclusion

In this paper, using a special transformation of variables we have reduced the full set of hydrodynamic equations describing the unsteady flow around a rigid body in terms of threedimensional flow to a system of equations on the surface. These equations can greatly simplify the numerical simulation and explore the deeper features of the flow process. First, they reduce the dimension of the task by a unit; there is no need to solve the hydrodynamic equations in the boundary layer. Second, besides the gas velocity at the surface they allow defining how all other parameters characterizing the flow (such as, u , P , ω etc.) change at the boundary.

It should be noted that not all of the obtained equations are with respect to time distinct systems. Therefore, in order to determine the stress distribution in addition to the initial data one should require examining boundary conditions. Through the previous, the information on the external flow is undoubtedly affecting the evolution of the whole effect of the flow on the body. They may need to be defined numerically using the additional code in space, the accuracy of which is important, but not in the entire volume, and only near the curves on the surface of the body along which the boundary conditions are determined.

The systems of equations obtained in this paper are essentially three-dimensional ones, ie, degenerate into the two-dimensional case. In the two-dimensional case (see Fig. 2) in order to override the hydrodynamic equations, it is necessary instead of transformations (15), (45) to use a different transformation similar to (A.12) (see Appendix A):

$$\frac { d \mathbf r } { d t } = \mathbf u + \mathbf \alpha + \frac { \left ( \omega + \nabla \times \mathbf \alpha \right ) } { \left | \omega + \nabla \times \mathbf \alpha \right | ^ { 2 } } \times \nabla \beta \,, \quad \mathbf r = \mathbf r ( \mathbf r _ { 0 }, t ) \,, \, \mathbf r _ { 0 } = \mathbf r _ { 0 } ( \mathbf r, t ) \,, \, \mathbf r _ { 0 } \right | _ { \mathbf r = 0 } = \mathbf r \,.$$

An additional condition should be implemented instead of (16)

$$\left ( \frac { \left ( \omega + \nabla \times \alpha \right ) } { \left | \omega + \nabla \times \alpha \right | ^ { 2 } } \times \nabla \beta \right ) \nabla \rho + \rho \nabla \beta \cdot \left ( \nabla \times \frac { \left ( \omega + \nabla \times \alpha \right ) } { \left | \omega + \nabla \times \alpha \right | ^ { 2 } } \right ) = 0 \,.$$

Parameter  is determined so that the boundary of the body did not move in the coordinates 0 r (see (30) and (60)).

Despite the complexity , our description can successfully be used in applications, such as less  expensive  calculation  of  lift  or  drag  coefficient  in  the  external  wing  of  impingement (including supersonic) flow, as well as the walls of the heat exchangers of nuclear reactors. The

Approach developed here is quite general and can be applied to other problems [14, 17]. If a more simple way is subsequently proposed to override the Navier-Stokes equations, the number of computations significantly reduced. This article may be of interest to researchers engaged in finding the analytical solutions of the equations of hydrodynamics. For example, according to the proposed method the hydrodynamic equations in the volume can be reduced to a system of equations  on  the  surface  and  even  get  an  overdetermined  system  of  surface  equations. Consequently, it can already be reduced to an overdetermined system of equations along a curve on the surface, etc. up to the analytical solution.

It should be noted that the override method itself is quite simple and already described by the authors [14] where a method of overriding the hydrodynamic equations described in Sec. 2 and Sec. 4 is provided. However, in this paper we apply it to a specific practical problem.

## Appendix A

Integrals of motion

Consider the Euler equations of an ideal incompressible fluid in the 3D case, in the form [1]

$$\frac { \partial \omega } { \partial t } = \nabla \times [ \mathfrak { u } \times \omega ],$$

$$\overline { \partial t }$$

$$\text{div} \, \iota = 0 \,.$$

We also consider the Lagrangian variables (i.e. markup) ) . Multiply both sides of , ( t r r r 

0

0 0 (A.1) to x  . Then

$$\nabla x _ { 0 } \cdot \frac { \partial \omega } { \partial t } = \nabla x _ { 0 } \cdot \nabla \times [ \mathfrak { u } \times \omega ] = - \text{div} ( \nabla x _ { 0 } \times [ \mathfrak { u } \times \omega ] ) = - \text{div} ( \mathfrak { u } ( \omega \cdot \nabla x _ { 0 } ) ) + \text{div} ( \omega ( \mathfrak { u } \cdot \nabla x _ { 0 } ) ) =$$

$$= - \mathfrak { u } \cdot \nabla ( \mathfrak { v } \cdot \nabla x _ { 0 } ) + \mathfrak { v } \cdot \nabla ( \mathfrak { u } \cdot \nabla x _ { 0 } ) = - \mathfrak { u } \cdot \nabla ( \mathfrak { v } \cdot \nabla x _ { 0 } ) - \mathfrak { v } \cdot \nabla \left ( \frac { \partial x _ { 0 } } { \partial t } \right ).$$

Here we have used the property of Lagrange variables:

$$\frac { \partial x _ { 0 } } { \partial t } + u _ { x } \, \frac { \partial x _ { 0 } } { \partial x } + u _ { y } \, \frac { \partial x _ { 0 } } { \partial y } + u _ { z } \, \frac { \partial x _ { 0 } } { \partial z } = 0 \,.$$

Consequently,

$$\frac { \partial } { \partial t } \left ( \mathfrak { w } \cdot \nabla x _ { 0 } \right ) + \mathfrak { u } \cdot \nabla ( \mathfrak { w } \cdot \nabla x _ { 0 } ) = \frac { d } { d t } \left ( \mathfrak { w } \cdot \nabla x _ { 0 } \right ) = 0 \Rightarrow \mathfrak { w } \cdot \nabla x _ { 0 } = \omega _ { 0 x } ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) \,,$$

where         u t dt d / / -  time  derivative  in  the  Lagrangian  variables, , , ( z y x ω

vorticity distribution ω . Similarly 0 0 0 0 0 , , z y x y y ω   0 0 0 0 0 , , z y x z z     ω .

) -0 0 0 0 initial       and

From the continuity equation (A.2) one follows

$$& \text{from one community equation (A.z.) one follows} \\ & \text{div} \, = \, 0 \, \Rightarrow \, \left ( \frac { \partial ( u _ { x }, y, z ) } { \partial ( x, y, z ) } + \frac { \partial ( x, u _ { y }, z ) } { \partial ( x, y, z ) } + \frac { \partial ( x, y, u _ { z } ) } { \partial ( x, y, z ) } \right ) = 0 \,,$$

$$& \partial _ { \Delta } \cdot \frac { d \Delta } { d t } = 0 \,, \\ & \tau. e. \, \frac { \partial ( x, y, z ) } { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } = 1.$$

$$\ t. e. \, \frac { \partial ( x, y, z ) } { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } = 1.$$

Expressing ω using (A.3), we also have   x t z y x x 0 0 0 0 0 , , ,    ω  ,   y t z y x y 0 0 0 0 0 , , ,   ω  ,   z t z y x z 0 0 0 0 0 , , ,    ω  , где    0 0 0 0 , , x y z        (see [1] page 31).

We also consider the general equation of continuity

$$\frac { \hat { \partial } \rho } { \hat { \partial } t } + \mathfrak { u } \cdot \nabla \rho + \rho \cdot \text{div} \mathfrak { u } = 0 \,.$$

Passing to the Lagrangian variables, we find [6]

$$& \text{assing to one Lagrangian variances, we in} \, \{ 0 \} \\ & \frac { d \rho } { d t } + \rho \cdot \left ( \frac { \partial ( u _ { x }, y, z ) } { \partial ( x, y, z ) } + \frac { \partial \left ( x, u _ { y }, z \right ) } { \partial ( x, y, z ) } + \frac { \partial ( x, y, u _ { z } ) } { \partial ( x, y, z ) } \right ) = 0 \,,$$

$$\frac { d \rho } { d t } + \frac { \rho } { \Delta } \cdot \frac { d \Delta } { d t } = 0 \,,$$

i.e.

$$\rho \cdot \Delta = \rho \cdot \frac { \partial ( x, y, z ) } { \partial ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) } = \rho _ { 0 } ( x _ { 0 }, y _ { 0 }, z _ { 0 } ) \text{ - integral of motion.}$$

In the case const   , but const s  , there is a clear connection between the pressure and density    S P P  the corresponding integrals are

$$\frac { \omega \cdot \nabla x _ { _ { 0 } } } { \rho } \,, \ \frac { \omega \cdot \nabla y _ { _ { 0 } } } { \rho } \, \bar { u } \ \frac { \omega \cdot \nabla z _ { _ { 0 } } } { \rho } \,.$$

Similarly, we obtain

$$\frac { \nu \left ( \mathbf r _ { 0 }, t \right ) } { \rho } = \frac { \omega _ { 0 } \cdot \nabla _ { 0 } \mathbf r } { \rho _ { 0 } },$$

The same results can be obtained if a more general transformation of the variables is used instead of Lagrangian variables

$$\frac { d \mathbf r } { d t } = \mathbf u ( \mathbf r, t ) + \lambda ( t ) \cdot \mathbf v ( \mathbf r, t ) \,,$$

where ) , ( 0 t r r r  и ) , , ( 0 0 t r r r  r r   0 0 t , а   t  - some function on t . Similar arguments can be established that in these variables, 0 x   ω , 0 y   ω , 0 z   ω and the Jacobian of this transformation  do not depend explicitly on time in an ideal fluid. Instead of   t  we can take   t z y x , , ,  , including      0 , , , , , ,  t z y x t z y x ω to satisfy   0. div   ω u  Expressing also ω , we obtain   0 0 0 , x t x     r ω ,   0 0 0 , y t y     r ω ,   0 0 0 , z t z     r ω

.

In the case const   , but const s  , there is a clear connection between the pressure and density    S P P  it is possible for the same reason in (A.5) to take   t z y x , , ,  , including 0 . Then the corresponding integrals are         ω ω

$$\frac { \omega \cdot \nabla x _ { 0 } } { \rho } = \frac { \omega _ { 0 _ { x } } ( { \mathbf r } _ { 0 } ) } { \rho _ { 0 } ( { \mathbf r } _ { 0 } ) } \,, \ \frac { \omega \cdot \nabla y _ { 0 } } { \rho } = \frac { \omega _ { 0 _ { y } } ( { \mathbf r } _ { 0 } ) } { \rho _ { 0 } ( { \mathbf r } _ { 0 } ) } \,, \ \frac { \omega \cdot \nabla z _ { 0 } } { \rho } = \frac { \omega _ { 0 _ { z } } ( { \mathbf r } _ { 0 } ) } { \rho _ { 0 } ( { \mathbf r } _ { 0 } ) } \, \mathfrak { u } \ \rho \Delta = \rho _ { 0 } ( { \mathbf r } _ { 0 } ) \,.$$

In fact, it requires only a negligible term 2 /    P  in equation (A.1), taking into account the compressibility of the fluid, which is performed at a characteristic velocity of its motion being much less than the speed of sound. [1].

$$\nabla \times [ \mathbf F \times \mathfrak { o } ] = 0 \,.$$

$$& \text{Going further. We solve the equation for the vector } F \\ & \nabla \times \left [ F \times \omega \right ] = 0 \,.$$

Consequently,

$$[ F \times \omega ] = \nabla \beta \,.$$

Then

$$\omega \nabla \beta = 0 \ u \ \omega \times [ F \times \omega ] = \omega \times \nabla \beta.$$

The solution of this equation has the form

$$F = \frac { \omega } { | \omega | ^ { 2 } } \times \nabla \beta + \lambda \omega \,.$$

We now change the variables

$$\frac { d \mathbf r } { d t } = \mathbf u + \frac { \omega } { \left | \omega \right | ^ { 2 } } \times \nabla \beta + \lambda \omega \,, \, \mathbf r = \mathbf r ( \mathbf r _ { 0 }, t ) \,, \, \mathbf r _ { 0 } = \mathbf r _ { 0 } \left ( \mathbf r, t \right ) \,, \, \mathbf r _ { 0 } \right | _ { t = 0 } = \mathbf r$$

and thus implement the conditions

$$\omega \nabla \beta = 0 \text{ and div} \left ( \frac { \omega } { | \omega | ^ { 2 } } \times \nabla \beta + \lambda \omega \right ) = \nabla \beta \cdot \left ( \nabla \times \frac { \omega } { | \omega | ^ { 2 } } \right ) + \left ( \omega \nabla \right ) \lambda = 0 \,.$$

Then these variables 0 x  ω , 0 y  ω , 0 z  ω and the Jacobian of this transformation  do not depend explicitly on time. Similarly, one can take into account compressibility.

In the two-dimensional case, the change of variables (A.10) has the form

$$\frac { d \mathbf r } { d t } = \mathbf u + \frac { \omega } { \left | \omega \right | ^ { 2 } } \times \nabla \beta \,, \, \mathbf r = \mathbf r ( \mathbf r _ { 0 }, t ) \,, \, \mathbf r _ { 0 } = \mathbf r _ { 0 } ( \mathbf r, t ) \,, \, \mathbf r _ { 0 } \right | _ { t = 0 } = \mathbf r \,,$$

where

$$\begin{array} { c } \psi \\ \frac { \mathfrak { v } } { \left | \mathfrak { v } \right | ^ { 2 } } \times \nabla \beta = \frac { 1 } { \omega } \begin{pmatrix} - \frac { \partial \beta } { \partial y } \\ \frac { \partial \beta } { \partial x } \end{pmatrix}. \\ \dots \quad \dots \dots. \end{array}$$

Condition (A.11) becomes

$$\frac { \partial \omega } { \partial x } \frac { \partial \beta } { \partial y } - \frac { \partial \omega } { \partial y } \frac { \partial \beta } { \partial x } = 0 \, \quad \text{ or } \quad \beta = G ( \omega, t ) \,.$$

As a result, the change of variables (A.12) can be written as

$$\frac { d \mathbf r } { d t } = \mathbf u + \mu \left ( \omega, t \right ) \nabla \times \omega \,, \, \mathbf r = \mathbf r \left ( \mathbf r _ { 0 }, t \right ) \,, \, \mathbf r _ { 0 } = \mathbf r _ { 0 } \left ( \mathbf r, t \right ) \,, \, \mathbf r _ { 0 } \Big | _ { t = 0 } = \mathbf r \,.$$

( А .14)

These variables allow

$$\omega = \omega _ { 0 } \left ( x _ { 0 }, y _ { 0 } \right ) \text{ and } \frac { \hat { \partial } ( x _ { 0 }, y _ { 0 } ) } { \hat { \partial } ( x, y ) } = 1 \,.$$

## REFERENCES

- 1. L. D. Landau and E. M. Lifshitz , Course of Theoretical Physics, Vol. 6: Fluid Mechanics (Nauka, Moscow,1986; Butterworth-Heinemann, Oxford, 1987).
- 2. Williams F.A. Combustion Theory. -Benjamin Cummings, Menlo Park, CA, 1985.
- 3. Ya.  B.  Zel'dovich,  G.  I.  Barenblatt,  V.  B.  Librovich,  and  G.  M.  Makhviladze , The Mathematical  Theory  of  Combustion  and  Explosion  (Nauka,  Moscow,  1980;  Consultants Bureau, New York, 1985).
- 4. A. N. Tikhonov and A. A. Samarskii, Equations of Mathematical Physics (Nauka, Moscow, 1966; Dover, New York, 1990).
- 5. A.  A.  Samarskii  and  Yu.  P.  Popov, Difference  Methods  for  Solving  Problems  of  Gas Dynamics (Nauka, Moscow, 1980) [in Russian].
- 6. M.L.  Zaytsev,  V.B.  Akkerman, A  Nonlinear  Theory  for  the  Motion  of  Hydrodynamic Discontinuity Surfaces,  JETP, 108 , No 4, pp. 699-717 (2009)
- 7. M.L. Zaytsev and V.B. Akkerman, Method for Describing the Steady-State Reaction Front in a Two-Dimensional Flow, JETP Letters, 92 , No 11, pp. 731-734 (2010)
- 8. Bychkov V., Zaytsev M. and Akkerman V. Coordinate-free description of corrugated flames with realistic gas expansion // Phys. Rev. E. - 2003. - V. 68. -P. 026312.
- 9. Joulin G., El-Rabii H. and K. Kazakov K. On-shell description of unsteady flames // J. Fluid Mech. - 2008. - V. 608. - P. 217.
- 10. Ott E. Nonlinear evolution of the Rayleigh-Taylor instability of thin layor // Phys. Rev. Lett. -1972. - Vol. 29.  - P. 1429-1432.
- 11. Book  D.L,  Ott E., Sulton A.L. Rayleigh-Taylor  instability in the "shallow-water" approximation // Phys. Fluids. -1974. - Vol.17. N4. - P.676-678.
- 12. Зайцев M. Л ., Аккерман В . Б . Метод снижения размерности в задачах арогидродинамики // Труды 52ой научной конференции МФТИ " Современные проблемы фундаментальных и прикладных наук ". -2009. Часть VI. С . 10-13.
- 13. Зайцев M. Л ., Аккерман В . Б . Метод снижения размерности в переопределенных системах дифференциальных уравнений в частных производных // Труды 53ой научной конференции МФТИ " Современные проблемы фундаментальных и прикладных наук ". -2010. Часть VII, Т . 3. С . 41-42.
- 14. V.B.  Akkerman  and  M.L.  Zaytsev, Dimension  Reduction  in  Fluid  Dynamics  Equations, Computational Mathematics and Mathematical Physics, Vol. 8 No. 51, pp. 1418-1430 (2011)

- 15. G. V. Korenev, Tensor Calculus (Moscow Institute of Physics and Technology, Moscow, 1996) [in Russian].
- 16. J. A. Schouten, Tensor Analysis for Physicists (Clarendon, Oxford, 1954; Nauka, Moscow, 1965).
- 17. M.L. Zaytsev, V.B. Akkerman, Free Surface and Flow Problem for a Viscous Liquid, JETP, 113 , No 4, pp. 709-713 (2011)

## FIGURES

Fig. 1. Rigid body in a viscous medium

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0153v1-figure-cc25f551-d92f-44ca-b810-74a4dcbb8b94.jpg)

Fig. 2. A solid in viscous fluid in the two-dimensional case

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0153v1-figure-07104ed1-4d11-4894-bea7-8af9cf8f60ab.jpg)
Maxim Zaytsev Vyacheslav Akkerman
Flow problem in three-dimensional geometry
It is shown how a complete set of hydrodynamic equations describing an
unsteady three-dimensional viscous flow nearby a solid body, can be reduced to
a closed system of surface equations using the method of dimension reduction of
over-determined systems of differential equations. These systems of equations
allow determining the surface distribution of the resulting stresses on the
surface of this body as well as all other quantities characterizing the
hydrodynamic flow around it.