Scientific journal
European Journal of Natural History
ISSN 2073-4972

STUDY OF NAVIER – STOKES EQUATION SOLUTION I. THE GENERAL SOLUTION OF NONLINEAR ORDINARY DIFFERENTIAL EQUATION

Yakubovskiy E.G.
Values of large dimensionless unknown functions (for example, a large Reynolds number) can be found out as solutions of non-linear partial differential equation. In this case these equations can be brought to some number of non-linear ordinary differential equations. Turbulent solutions corresponding to large values of unknown function are complex. Transition from real solution to complex turbulent solution is realized through infinity of the right parts of ordinary differential equation system to which Navier – Stokes equations are brought. Thus, real solution of Navier – Stokes equation for turbulent mode yields function going to infinity. At the same time, complex solution for the turbulent mode is finite. Fluid flow resistance coefficient is calculated for round pipeline with different pipeline walls roughness.
fluid flow resistance coefficient
turbulent function
Navier – Stokes equation
solutions of non-linear partial differential equation
large dimensionless unknown functions
round pipeline

Problem Formulation

Let us consider Navier – Stokes problem and continuity equation for incompressible fluid. They are as following

yakubovsk001.wmf

Boundary conditions on the body boundary adjoining to fluid are V(t, r) = 0, r ∈ S where S is a body boundary. We will seek a solution in the form of series using Galerkin method (hereinafter N → ∞)

yakubovsk002.wmf

yakubovsk003.wmf

where space C2 is twice continuously differentiable function, Ψ0(r) is a defined external action which, in case of pipeline, is equal to

yakubovsk004.wmf

where z – direction of the pipeline longitudinal; P, P0 – pressure at the beginning and the end of the pipeline; L – pipeline length.

Now we substitute these functions into the differential equation, multiply by Ψm(r) and integrate over the volume, then we obtain following differential equations system:

yakubovsk005.wmf m = 1,..., 3N;

yakubovsk006.wmf m = 1,..., N. (A.1)

After we resolved the second equation (A.1), substituted

yakubovsk007.wmf n = 1,..., N

from the second equation (A.1) to the first one, we have

yakubovsk008.wmf m = 1,..., 2N;

yakubovsk009.wmf m = 2N + 1,..., 3N. (A.2)

Defining xn+3N(t), n = 1,..., N corresponding to pressure change, from the second equation (A.2) and substituting found out value into the first equation (A.2), we have equations system

yakubovsk010.wmf m = 1,..., 2N. (A.3)

At that

yakubovsk011.wmf yakubovsk012.wmf n = 1,..., N.

Where values

yakubovsk013.wmf l = 1,..., 3; yakubovsk014.wmf n = 1,..., N.

are known and coefficients Fmpq, Gmp, Hm, yakubovsk015.wmf yakubovsk016.wmf yakubovsk017.wmf yakubovsk018.wmf yakubovsk019.wmf yakubovsk020.wmf cnm, bnm are constants. This system of non-linear ordinary differential autonomous equations (A.3) is to be solved. Solution convergence issues will be discussed below in the text.

Finding of Solution of Ordinary Differential Equations in Complex Plane

Let us consider system of non-linear differential autonomous equations

yakubovsk021.wmf l = 1, ..., N. (1)

Navier – Stokes equation system and continuity equation can be brought to system of non-linear differential equations:

yakubovsk022.wmf m = 1,..., N, (2)

where three-dimensional velocity is defined by formula yakubovsk023.wmf At that, function φn(x1, x2, x3) is given in the form of sine or cosine. Then coefficients cn(t) for continuous function decrease not more rapidly than yakubovsk024.wmf when index increases and series reduction is possible, i.e. instead of infinite number of terms, finite terms number is used. At the same time, the infinite number of terms forms convergent series.

It was found out that a set of N + 1 coordinates for the system balance position exists (2). Indeed, let us assume that we have found several balance positions with coordinates yakubovsk025.wmf l = 1, ..., N. Let us seek the solution in the form yakubovsk026.wmf. For that we will substitute the solution into the right part of the differential equation (2) and will equate it to zero, then following equations system is obtained

yakubovsk027.wmf

For existence of non-zero solution of this differential equation, it is necessary that determinant is equal to zero:

yakubovsk028.wmf

Given zero determinant, coefficients yakubovsk029.wmf from linear equation will be defined up to a multiplier. This multiplier will be defined from equality to zero of determinant of non-linear equation system. Thus, we have N unknown multipliers, which will be defined from determinant equality to zero. I.e. set of N + 1 coordinates of the system balance position exists.

Differential equation system (2) for non-multiple balance positions can be expressed by means of yakubovsk030.wmf substitution. At that, the system (2) balance positions yakubovsk031.wmf l = 1, ..., N, s = 1, ..., S will be transformed into balance positions yakubovsk032.wmf l = 1, ..., N, s = 1, ..., S and eigen values and eigen vectors of the linearized system (2) will be defined as.

yakubovsk033.wmf

yakubovsk034.wmf

Equation system (2) will be

yakubovsk035.wmf (3а)

Values yakubovsk036.wmf satisfy condition yakubovsk037.wmf. k = 1, ..., N, s = 1, ..., S.

Equation system (3a) can be written as

yakubovsk038.wmf (3b)

where multiplier which can never be equal to 0 is used – yakubovsk039.wmf, and this multiplier is equal to

yakubovsk040.wmf

After this multiplier is substituted to (3b) we obtain (3a). Now we will demonstrate that this multiplier can never be equal to 0. When yakubovsk041.wmf l = 1, ..., N the limit of

yakubovsk042.wmf

is finite.

Here we canceled out a multiplier yakubovsk043.wmf, as we consider only not coincident roots being coordinates of balance position. So we showed that this multiplier does not equal to zero after infinite time.

Thus, the differential equation can be written as

yakubovsk044.wmf yakubovsk045.wmf l = 1, ..., N, (4)

where Hl(t, t0) – function which tends to infinity when coordinates tend to balance position. For real solutions, this function is monotonic. That is, we have obtained dependence of the solution on value Hl(t, t0), which is monotonic time-dependent function.

Lemma 1. Necessary and sufficient criterion for unknown function to tend to steady balance position coordinates is Hl(t, t0) → ∞ when t → ∞. At the same time, balance position coordinates have to have a real part.

So, we have

yakubovsk046.wmf (5)

at t → ∞ and hence Hl(t, t0) → ∞ , l = 1, ..., N as integral of constant. Inverse theorem is also valid, on condition Hl(t, t0) → ∞, l = 1, ..., N, one of steady balance positions is realized. This is a consequence of solution type; on condition Hl(t, t0) → ∞, l = 1, ..., N, according to Lemma 3, negative real part of value yakubovsk047.wmf exists in formula (6) and solution tends to balance position coordinate yakubovsk048.wmf in formula (4). If balance position coordinates have real parts, values yakubovsk049.wmf have real part. At that t → ∞.

Lemma 2. Solution of differential equation (1) is function xl(t) which satisfies to equation (6).

To obtain (6), let us divide equation (4) by product of multipliers yakubovsk050.wmf and multiply (4) by dHl(t, t0). Then we will decompose obtained fraction into sum of simple fractions and perform integration. The following equation is obtained

yakubovsk051.wmf

l = 1, ..., 2N.

Here for the case of sound energy emission in interval [t0, t] different branches of logarithm are obtained.

After the expression exponentiated, we have (6)

yakubovsk052.wmf

yakubovsk053.wmf (6)

where all values of balance position coordinates are not multiple and are not dependent on radiation process occurring in an interval [t0, t]. In case of laminar real solution, radiation will not appear, and in case of turbulent solution, followed by radiation, there will be energy transition. Really, presence of radiation yields the complex solution which describes turbulent pulsing mode. At that, at solution transformation, turbulent mode is followed by sound noise. Exponential multiplier does not affect balance position coordinates which define stationary solution. Existence of multiplier yakubovsk054.wmf changes calculated main branch of solution for coordinate xl, but will not affect balance position coordinate.

Lemma 3. Sum of coefficients yakubovsk055.wmf by index s is equal to zero, i.e. yakubovsk056.wmf

In case if following fraction decomposed.

yakubovsk057.wmf

where QS–1(y) is S – 1-ordered polynomial. Equation yakubovsk058.wmf will remain satisfied,

yakubovsk059.wmf

Let us prove this. For this let us consider a sum

yakubovsk060.wmf

This sum is equal to P(y) = QS–1(y). We write formula for polynomial equal to QS–1(y), dividing the equation by product yakubovsk061.wmf we obtain

yakubovsk062.wmf

If suppose that yakubovsk063.wmf, equality yakubovsk064.wmf is satisfied when s + 1 balance position exists.

But to realize the solution, it is necessary to know balance positions of this non-linear equations system. Besides, balance positions can be multiple that changes the solution finding process, it becomes random or chaotic, but we are not going to consider this case. Nevertheless, it is possible to prove the following important theorem.

Theorem 1. Cauchy task is considered under arbitrary real initial conditions for system of orthogonal non-linear ordinary differential equations (1). If system (1) has complex conjugate balance positions with real parts then, for finite real argument t, Cauchy problem solution for the system (1), for real initial conditions, tends to infinity. Then this solution becomes a complex one, tending to balance position in case when complex balance position coordinates have real part. Here the right part of (1) is considered as being a regular function, real for real arguments. This function has finite number of non-multiple balance positions.

Proving

If the system (2) is resolved at non-multiple balance positions then, according to Lemma 2, we have

yakubovsk065.wmf (7)

where yakubovsk066.wmf selected complex balance position, yakubovsk067.wmf other balance positions. Coefficients yakubovsk068.wmf satisfy condition yakubovsk069.wmf, according to Lemma 3. At that, in sum yakubovsk070.wmf real part value yakubovsk071.wmf in case of complex solution yakubovsk072.wmf presents twice as all values yakubovsk073.wmf satisfy condition yakubovsk074.wmf, so we have formula

yakubovsk075.wmf

Let us substantiate solution (7). For that we will modify two complex conjugate terms of the solution (for expression simplicity, index l is omitted)

yakubovsk076.wmf (8)

where yakubovsk077.wmf. After integration (8) over argument x, we obtain formula (7)

yakubovsk078.wmf

The solution is yakubovsk079.wmf where

yakubovsk080.wmf

yakubovsk081.wmf

At that, value of yakubovsk082.wmf is real due to existence of complex conjugate balance positions. Thus, for yakubovsk083.wmf and finite t, we have equation

yakubovsk084.wmf (9)

Solution of this equation tends to infinity.

At that, solution of differential equation for rising Hl(t, t0), according to Lemma 1, can have complex roots

yakubovsk085.wmf

At that, as equation yakubovsk086.wmf is satisfied according Lemma 3 and balance positions have real parts, values with negative real part yakubovsk087.wmf exist, so convergence to one of the balance positions takes place. Real solution will tend to infinity at that existence, and uniqueness condition for Cauchy problem will be breached. According to Lemma 1, at Hl(t, t0) infinity, unknown function will tend to one of balance positions. This balance position cannot be real as the real solution is infinite. This means that the solution will have a branching point and will tend to complex balance position. That is, for balance complex positions, finite complex solution is obtained at Hl(t, t0) change. Thus, in some point a complex solution will begin.

End of the proof.

Now we will give an example describing this property of the differential equation, transition to the complex solution. So, for the differential equation, there can be a complex solution instead of infinite real one

yakubovsk088.wmf

And these balance positions are purely imaginary, that is, the solution cannot tend to balance position. And the real solution of this differential equation tends rapidly to infinity

yakubovsk089.wmf

Using an implicit solution finding scheme, we obtain the following equation

yakubovsk090.wmf

Seeking solution in respect to unknown function x, we obtain the following implicit scheme

yakubovsk091.wmf

This implicit scheme with constant step correctly describes solution tendency to infinity. At reduced calculation step, this scheme yields larger value of variable t, that is, it yields larger value of unknown function. That is, it correctly describes behavior of the differential equation solution up to infinity. When infinity is reached, under condition yakubovsk092.wmf, the finite complex solution will be found. Numerical computation of this equation has validated this analysis of the solution obtained.

At that, the complex solution possesses new properties; it performs complex rotation around balance position. At the same time the real solution tends to infinity, i.e. right part of the differential equation tends to infinity and existence and uniqueness condition for Cauchy problem are breached, so additional complex solution is arisen.

The solution for complex initial data is given by formula

yakubovsk093.wmf

for any t. Thus, approximately we have

yakubovsk094.wmf

If we choose branch with positive β, we obtain converging series. At that, this fraction denominator never becomes zero.

That is, for real plane, finite solution does not exist. In complex plane, finite continuous solution exists in the case if balance positions are not multiple.

But there is a question – what is the physical meaning of imaginary part of the solution?

Physical Meaning of Exact Complex Solution

So, for turbulent solution corresponding to complex balance position coordinates, we have solution

yakubovsk095.wmf

The solution consists of step term in the form of delta-function and smooth part

yakubovsk096.wmf

As, at averaging over period, tangents sum without taking into account step term is equal to zero, we will study the step term of the solution. At that, this solution has singularity when condition yakubovsk097.wmf is satisfied. Step term of the solution is

yakubovsk098.wmf

That is, imaginary medium pulse is originated. Imaginary velocity means flow rotation or oscillation; flow step is originated which will be destroyed in time ΔD(hl) = π to originate repeatedly. Number of such steps is finite. But how to average this steps? You should pass to probabilistic interpretation of the description. That is, to average imaginary part over the period yakubovsk099.wmf Then we have local complex average solution

yakubovsk100.wmf

yakubovsk101.wmf

Continuous part of the solution has positive and negative parts which are compensated when averaged. To obtain a global average value it is necessary to average with respect to value k, so we have

yakubovsk102.wmf

We obtained complex velocity; imaginary part is defined up to multiplier. Real part of complex velocity corresponds to average value of velocity, and imaginary part is a mean square deviation. Simultaneously, there is a vortex motion consisting of positive and negative value of root from yakubovsk103.wmf.

Contribution of imaginary part to average value is equal to

yakubovsk104.wmf yakubovsk105.wmf

At that, module of average value, that defines real solution, is equal to

yakubovsk106.wmf yakubovsk107.wmf

where balance position coordinates and time are non-dimensional, then, as we calculate square root of imaginary part, we define branch yakubovsk108.wmf. Thus, average single-valued solution is found.

This multiplier Υl depends on the surface roughness and it is found from numerical experiment. As numerical experiment has shown, for round smooth pipeline the multiplier is equal to Υl = 1. At that, the smooth pipeline has a constant, minimum, average module of roughness inclination tangent equal to

yakubovsk109.wmf

that is associated with molecular roughness, see article II section 1. For this, one term of series which determines flow velocity is used. We calculate this value for one term of the series for smooth surface. The solution is

yakubovsk110.wmf

where yakubovsk111.wmf

yakubovsk112.wmf

yakubovsk113.wmf

This value exactly corresponds to experimental formula for round cross section pipeline if the solution roughness is taken into account. At that, to take roughness into account for internal problem, yakubovsk114.wmf is multiplied by yakubovsk115.wmf, here k/l is a constant average tangent of flowing surface inclination. From this we obtain yakubovsk116.wmf, see article II section 1. At constant average roughness height, the coefficient Υl is not a constant as k/l value is determined by other formula depending on dimensionless pressure, see article II section 1.