Navier-Stokes equations are partial differential equations governing the incompressible fluids motion [1] . These equations constitute the basic equations of fluid mechanics [2] . They are named after Claude Navier and George Stokes, two 19th century physicists. These equations are described as below read [3] [4] [5] :

∂ v ∂ t + v ⋅ ∇ v = f + α ∇ 2 v − 1 ρ ∇ P (1)

It is difficult to model a fluid as a single-phase continuous medium at the resolution of these equations. The mathematical existence of Navier-Stokes equations solutions was not demonstrated because, it is such a difficult problem to solve and so important (because its equations govern the water flow in a pipe, the ocean currents, the air movements in the atmosphere) [6] [7] . The complication comes from v ∇ v term. This term varies as velocity field square, and it has done the equation mathematically inextricable. We approach this by trying to find an analytical way of solving v ∇ v . After numerical verification, the solution and v ∇ v converge. Then we finish by resolving by this solution the Navier-Stokes equations as a differential equation. Finally we find an analytical solution that can help to resolve many fluid mechanics problems.

After a series of numerical and experimental studies [8] [9] [10] , we have found that v ∇ v converge as:

F ( v ) = ( N − 10 ) ∑ n = 1 + ∞ v 500 n + 2 n ! − 10 − 2 ∑ p = 0 N − 10     p v (2)

After simulation we obtain this Area 1 for F ( v ) with v random Figure 1: To compare with the Area 2 generates numerically by v ∇ v with v random, we simulate it with python Figure 2: We can see that the two areas are the same, which confirms that:

v ∇ v = ( N − 10 ) ∑ n = 1 + ∞ v 500 n + 2 n ! − 10 − 2 ∑ p = 0 N − 10     p v (3)

Then Navier Stokes equations become with ν = α the viscosity constant:

∂ v ∂ t + ( N − 10 ) ∑ n = 1 + ∞ v 500 n + 2 n ! − 10 − 2 ∑ p = 0 N − 10     p v = f + α ∇ 2 v − 1 ρ ∇ P (4)

∂ v ∂ t + ( N − 10 ) ∑ n = 1 + ∞ v 500 n + 2 n ! − ( N − 10 ) ( N − 9 ) 2 × 10 2 v = f + α ∇ 2 v − 1 ρ ∇ P (5)

Let v = h ( t ) g ( r ) with r ∈ { x , y , z } a coordinate:

1 h ∂ h ∂ t + ( N − 10 ) ∑ n = 1 + ∞ v 500 n + 1 n ! − ( N − 10 ) ( N − 9 ) 2 × 10 2 = f v + α 1 g ∇ 2 g − 1 v ρ ∇ P (6)

with 1 h ∂ h ∂ t = − K 2 ; K 2 = 1 τ c the temporal decay, we obtain:

h = h 0 exp ( − t τ c ) (7)

Let’s take ω 2 = 1 τ c + ( N − 10 ) ( N − 9 ) 2 × 10 2 we get:

∬ f d r 2 + α ∬ d 2 v − 1 ρ ∬ ∇ P d r 2 − ( N − 10 ) ∑ n = 1 + ∞ ∬ v 500 n + 2 n ! d r 2 + ω 2 ∬ v d r 2 = 0 (8)

f r 2 2 + α v 2 2 − P r ρ − ( N − 10 ) ∑ n = 1 + ∞ 1 n ! ∬ v 500 n + 2 d r 2 + ω 2 ∬ v d r 2 = 0 (9)

with g = exp ( a r ) we have d v = a v d r that implies d r = 1 a d v v , then we get:

a 2 f r 2 2 + a 2 α v 2 2 − a 2 P r ρ − ( N − 10 ) ∑ n = 1 + ∞ 1 n ! ∬ v 500 n d v 2 + ω 2 ∬ 1 v d v 2 = 0 (10)

a 2 f r 2 2 + a 2 α v 2 2 − a 2 P r ρ − ( N − 10 ) ∑ n = 1 + ∞ 1 n ! ( 500 n + 1 ) ( 500 n + 2 ) v 500 n + 2 − ω 2 v t τ c + a 2 ω 2 ∫ r v d r = 0 (11)

a 2 f r 2 2 + a 2 α v 2 2 − a 2 P r ρ − ( N − 10 ) ∑ n = 1 + ∞ 1 n ! ( 500 n + 1 ) ( 500 n + 2 ) v 500 n + 2 + ω 2 [ − v t τ c + a 2 ∫ r v d r ] = 0 (12)

1 ω 2 [ a 2 f r 2 2 + a 2 α v 2 2 − a 2 P r ρ − ( N − 10 ) ∑ n = 1 + ∞ 1 n ! ( 500 n + 1 ) ( 500 n + 2 ) v 500 n + 2 ]   − v t τ c + a 2 ∫ r v d r = 0 (13)

D ρ D t = − a v (14)

( ∂ ∂ t + v . ∇ ) ρ = − a v (15)

( ∂ ∂ t + v . ∇ ) ρ v = − a (16)

− ( ∂ ∂ t + v . ∇ ) ρ v = a (17)

with λ = a , we get:

{ λ l i q u i d = − ( ∂ ∂ t + v . ∇ ) ρ         foraliquid   ρ l i q u i d ≫ v f l u i d λ g a z = − ∇ ρ           foragaz   ρ g a z ≤ v g a z

τ c must be determined experimentally according the type fluid.

v ( r ∈ { x , y , z } ) = [ λ 2 α ∑ n = 1 + ∞ [ 2 ( λ r − t τ c ) ] n 2 n ! − ( N − 10 ) ∑ n = 1 + ∞     ∑ j = 1 + ∞ ( 500 n + 2 ) j − 1 n ! j ! ( 500 n + 1 ) ( λ r − t τ c ) j + λ 2 f r 2 2 − λ 2 P r ρ ] τ c ω 2 t     + λ 2 τ c t ∫ ∑ n = 1 + ∞     ( λ r 1 − n − t r − n τ c ) n n ! d r (18)

The analytical method only compensates solutions in which non-linear and complex structures in the Navier-Stokes equations are ignored within several assumptions [11] [12] [13] . Our analytical solution is only valid for fundamental cases such as fluids dynamics under an evenly distributed force and constant Volumic mass.

The TOUGMA’s solution is crucial to understand the fluid properties. The range of applications is extraordinarily broad [14] [15] :

1) External aerodynamics: airplane, glider, missile, launcher, space probe, automobile, flying insects, buildings and bridges;

2) Enternal aerodynamics: aircraft engines, gas turbines, rocket engines, air conditioning;

3) Hydrodynamics: boat, submarine, marine propulsion;

4) Hydraulics: pipes, open channels, waves, rivers, blood circulation; meteodynamics: meteorology, climatology;

5) Astrophysics: formation of stars, stellar systems, galaxies, the universe, stellar jets;

6) Armament: explosions;

7) Safety: table fires, forest fires and confined fires;

The present investigation will be very helpful to the researchers who are engaged for those area research works in earth and in Universe.

The present study received no financial aid from any government/non-government agency.

There are no conflicts of interests for this study.

Tougma, J.L.W. (2023) Method of Analytical Resolution of the Navier-Stokes Equations. Open Journal of Fluid Dynamics, 13, 226-231. https://doi.org/10.4236/ojfd.2023.135017

P: Position

v: Velocity

α : is the fluid coefficient viscosity

τ c : Experimental time

K 2 : Temporal decay

λ : Flux wave length

l c : Experimental critical length

ρ : Volumic mass

D D t : is a material derivative, stated as ∂ ∂ t + v ⋅ ∇

∇ : Divergence