Solitonlike Coherent Structure Is a Universal Motion Form of Viscous Fluid

Preprint

Article

Solitonlike Coherent Structure Is a Universal Motion Form of Viscous Fluid

Altmetrics

Downloads

270

Views

124

Comments

Bohua Sun^*

This version is not peer-reviewed

Submitted:

22 October 2023

Posted:

23 October 2023

Read the latest preprint version here

Alerts

Abstract

The discovery of solitonlike coherent structure (SCS) in boundary layer flows is crucial for understanding turbulence origins and in particular for laminar-turbulence transition. However, the task of finding solutions for solitonlike coherent structure from the Navier-Stokes equations poses a significant challenge. In this paper, based on the author's previous work [Sun (2023) in Ref.31], we are able to study convergent flow boundary layers, whose solution encompasses both shock wave and solitary wave solutions, and their superposition gives rise to solitary-like waves, namely solitonlike coherent structure. It is found that the solitonlike coherent structure can only be obtained by combining the Navier-Stokes equations and mass conservation, since the combined equations will have the third order derivatives respect to coordinate x, in other words, without the mass conservation condition, the Navier-Stokes equations does not contain solitonlike coherent structure by itself. Finally proved solitonlike coherent structure is a universal motion form of viscous fluid.

Keywords:

Subject: Physical Sciences - Fluids and Plasmas Physics

1. Introduction

Origins of turbulence is a ubiquitous phenomenon and one of the greatest unsolved mysteries of classical physics [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. Although the phenomenon of turbulence is still not yet fully understood, some universal structures of turbulence have been discovered, such as, in an incompressible boundary layer, typical coherent structures include ∧-vortices, solitonlike coherent structure (SCS), hairpin vortices (with their asymmetric and composite variations), pocket structure, typical eddies (TEs), high/low-speed streaks(LSSs), streamwise vortices, ringlike vortices, turbulent spots, dark spots, very large scale motions (VLSMs), and so on.

In particular, Lee [12,16,21,22] has successfully discovered solitonlike coherent structure (SCS) in different types of transitional boundary layers and turbulent boundary layers through experiments. They have utilized Tomographic Particle Image Velocimetry (Tomo-PIV) and their own developed Near-Wall Particle Image Velocimetry (Near-Wall PIV) to capture three-dimensional measurement data of randomly occurring strip-like structures in the turbulent boundary layer. Furthermore, they have employed the Lagrangian tracking method to invert the solitonlike coherent structure, revealing that the fundamental structure of near-wall flow is a solitary-like wave, and concluded that boundary layer flows all possess a common structure, namely solitonlike coherent structure (SCS, which has some typical characteristics as a soliton, such as, traveling downstream at a finite speed with a strong localized pattern and decreases in the later stage of boundary-layer transition, in particular, SCS has the longest living time compared with other coherent structures in wall-bounded flow [12,16].

For the significant discovery of the existence of solitonlike coherent structure or solitary-like wave structures in boundary layer flows, we believe that it is necessary to provide theoretical understanding and explanation. To do so, it is important to investigate whether there exist solitary wave solutions in the Navier-Stokes equations or in the non-steady boundary layer equations derived from the simplification of the Navier-Stokes equations using Prandtl’s boundary layer concept. Ideally, it would be best to find an analytical solution for the non-steady boundary layer equations that clearly includes solitary wave solutions. Unfortunately, to date, an analytical or exact solution to this problem has not been found [1,2,24,25,26,27,28,29,30,31].

Regarding water surface waves, researchers have observed the phenomenon of solitary waves and established the equation

u_{, t} + 6 u u_{, x} = ν u_{, x x x}

for water waves with a free boundary. They obtained the exact solution

u (x, t) = 2 ν c_{2}^{2} {tanh}^{2} (c_{1} + c_{2} x + c_{3} t) - (8 ν c_{2}^{3} + c_{3}) / (6 c_{2})

. However, the Navier-Stokes equations does not include third-order derivatives with respect to spatial coordinates, such as

u_{, x x x}

. Even in the 1D nonlinear Burgers equation

u_{, t} + u u_{, x} = ν u_{, x x}

, the solution

u (x, t) = - 2 ν c_{2} tanh (c_{1} + c_{2} x + c_{3} t) - c_{3} / c_{2}

is not a solitary wave solution but a shock wave solution.

This naturally raises a question: whether the Navier-Stokes equations or their simplified forms have solitary-like wave solutions? Our belief is that since the Navier-Stokes equations can describe the laws of fluid motion, if there is a common structure of solitonlike coherent structure in boundary layer flows, then the non-steady boundary layer equations of the Navier-Stokes together with mass conservation should have solitary-like wave solutions.

To achieve this belief, the best approach is to obtain exact solutions for the non-steady boundary layer equations of the Navier-Stokes (N-S) equations in certain cases. Due to the complexity of the N-S equations, only a few 1D non-steady problems have obtained exact solutions, such as the exact solution to Stokes’ problem, the flow due to an oscillating infinite plane, which was the first exact solution in fluid dynamics, and the solution to Rayleigh’s problem of an infinite flat plate started impulsively into motion in its own plane with velocity U, given by

u (y, t) = U erfc (y / (2 \sqrt{ν t}))

. Unfortunately, to date, no exact solutions have been obtained for the unsteady boundary layer equations of any 2D problem, and neither the solutions for Stokes nor Rayleigh include solitary-like wave solutions. This shattered our expectation of finding solitary-like wave solutions from historical exact solutions, leaving us with no choice but to face the problem and strive to find new ideas to explore new frontiers.

Lee et al [12,16,21,22] has found that boundary layer flows all exhibit a common structure, namely solitonlike coherent structure. This means that we only need to study one or a few specific cases, as universal laws exist within these special problems. By obtaining an analytical solution for any 2D unsteady boundary layer equations, we can analyze whether it contains solitary-like wave solutions.

In this paper, the author addresses the challenge of finding solitonlike coherent structures (SCS) in boundary layer flows, which are important for understanding turbulence origins and laminar-turbulence transition. The author builds upon their previous work and utilizes a combination of the Navier-Stokes equations and mass conservation to study convergent flow boundary layers.

By incorporating mass conservation, the combined equations now contain third-order derivatives with respect to the coordinate x. This inclusion is crucial, as it allows for the existence of solitonlike coherent structures. Without the mass conservation condition, the Navier-Stokes equations alone do not possess solitonlike coherent structures.

Through their analysis, the author demonstrates that the solution to convergent flow boundary layers encompasses both shock wave and solitary wave solutions. Furthermore, the superposition of these solutions gives rise to solitary-like waves, which are identified as solitonlike coherent structures.

The paper concludes by proving that solitonlike coherent structures exist in all kinds of flows, highlighting the significance of combining the Navier-Stokes equations with mass conservation to obtain these coherent structures. This finding contributes to our understanding of turbulence origins and the transition from laminar to turbulent flow.

After introduction in Section 1, the rest of this paper is organized as follows. In Section 2, based on Sun’s work [31], we formulate the 2D convergent boundary layers. In Section 3, an analytical solution is obtained. In Section 4, proved solitonlike coherent structure do exit in all kind of flows. Finally, in Section 5, conclusions and perspectives are drawn.

2. A thin plate immersed in unsteady convergent flow

A thin flat plate is immersed at zero incidence in a uniform stream, which flows with speed

U (x)

and is assumed not to be affected by the presence of the plate, except in the boundary layer. The fluid is supposed unlimited in extent, and the origin of coordinates is taken at the leading edge, with x measured downstream along the plate and y perpendicular to it.

For 2D unsteady boundary layers, Sun [31] introduced a novel similar transformation

ψ = U (x) δ (x) f (η, τ), η = \frac{y}{δ (x)}

and

τ = \frac{ν}{δ^{2}} t

, and universally converted the 2D unsteady boundary layers governing equation

\frac{\partial ψ}{\partial t \partial y} + \frac{\partial ψ}{\partial y} \frac{\partial ψ}{\partial x \partial y} - \frac{\partial ψ}{\partial x} \frac{\partial ψ}{\partial y^{2}} = U \frac{d U}{d x} + ν \frac{\partial^{3} ψ}{\partial y^{3}}

into a single partial differential equation as follows

\begin{matrix} f_{, η η η} & + α f f_{, η η} + β [1 - {(f_{, η})}^{2}] = f_{, τ η} + γ τ (f_{, τ} f_{, η η} - f_{, η} f_{τ η}), \end{matrix}

(1)

where the coefficients are

α = \frac{δ}{ν} \frac{d U δ}{d x}

β = \frac{δ^{2}}{ν} \frac{d U}{d x}

, and

γ = \frac{U}{ν} \frac{d δ^{2}}{d x} = 2 (α - β)

. If

α

and

β

are constants, we have main-stream flow

U (x) = C x^{m}

and boundary layer thickness

δ (x) = {[\frac{| (2 α - β) x | ν}{| U (x |)}]}^{1 / 2}

, in which the exponent

m = \frac{β}{2 α - β}

In case of boundary layer flow in a converging channel flow between two non-parallel planes, which has been studied by Pohlhausen [33], who’s solution can also be found in Page 162 in Fluid Mechanics by Landau and Lifshitz [34], the parameters

α = 0, β = 1

γ = - 2

and

m = - 1

, we have the convergent channel flow velocity function

U (x) = \frac{| Q |}{ρ θ x}

, which simply expresses the coacervation of the discharge

Q < 0

in the flow,

θ

being the angle between the planes. Now we have the boundary layer thickness

δ (x) = x {(\frac{ν ρ θ}{| Q |})}^{1 / 2}

, similarity variable

η = \frac{y}{δ} = \frac{y}{x} {(\frac{| Q |}{ν ρ θ})}^{1 / 2}

, and dimensionless diffusion time

τ = \frac{| Q |}{ρ θ} \frac{t}{x^{2}}

Referring to the basic problem is to find a function

f (η, τ)

satisfying Sun’s equation [26] as follows

\begin{matrix} f_{, η η η} + 1 - {(f_{, η})}^{2} = f_{, τ η} - 2 τ (f_{, τ} f_{, η η} - f_{, η} f_{τ η}) . \end{matrix}

(2)

Eq.2 can be solved numerically, some results are depicted in Figure 1.

3. Analytic solitonlike solution

Although we can solve the equation in Eq.2 numerically and find its solution, we do not know what kind of elements are contained in the solution since the solution is the final composition of those elements. This means that we need to split the solution into some elements. However, numerical methods would not be helpful, we need to perform some analytical analysis. Unfortunately, we can not find the exact solution of Eq. 2, so we have to find its approximate analytical solution.

Based on the arguments of Stewartson [24,25], Takuda[26] proposed a Figure 2, which can be used vividly to study this problem. From Figure 2, in the range of

0 < τ ≪ 1 / 2

, we can ignore the term

τ (f_{, τ} f_{, η η} - f_{, η} f_{τ η})

in Eq.2, hence the approximation of Eq.2 will be

\begin{matrix} f_{η η η} + 1 - {(f_{η})}^{2} \approx f_{τ η} . \end{matrix}

(3)

The exact solution of Eq.3 is obtained as follows

\begin{matrix} f (η, τ) & = - \sqrt{3} tanh ζ - 2 \sqrt{3} ln (cosh ζ) - \frac{\sqrt{3}}{2} ln \frac{tanh ζ - 1}{tanh ζ + 1} - \frac{η}{2} + g (τ), \end{matrix}

(4)

where

ζ = \frac{\sqrt{3}}{6} η - \frac{5}{6} τ - c

, and

g (τ)

is an arbitrary function, which can be determined by

f (0, τ) = 0

and gives

\begin{matrix} g (τ) & = - \sqrt{3} tanh (\frac{5}{6} τ + c) + 2 \sqrt{3} ln (cosh (\frac{5}{6} τ + c)) - \frac{\sqrt{3}}{2} ln \frac{1 - tanh (\frac{5}{6} τ + c)}{1 + tanh (\frac{5}{6} τ + c)} . \end{matrix}

(5)

Then we can get the first and the second partial derivative of

f (η, τ)

respect to

η

as follows

f_{, η} = - tanh ζ - \frac{1}{2} {sech}^{2} ζ

, and

f_{, η η} = - \frac{\sqrt{3}}{6} (1 - tanh ζ) {sech}^{2} ζ

The integral constant can be determined by

\frac{\partial f (η, τ)}{\partial η} |_{(0, 0)} = 0

and leads to

c = \frac{ln (1 + \sqrt{2})}{2}

The flow velocity component in x coordinate is obtained as follows

\begin{matrix} u & = \frac{| Q |}{ρ θ x} f_{, η} = \frac{ν R e}{x} (\frac{1}{2} {sech}^{2} ζ + tanh ζ), \end{matrix}

(6)

and since

U δ =

constant, so we have

\begin{matrix} v & = - [\underset{= 0}{\underset{︸}{{(U δ)}_{, x}}} f - U η δ_{, x} f_{, η} - 2 U τ δ_{, x} f_{, τ}] \\ = U η δ_{, x} f_{, η} + 2 U τ δ_{, x} f_{, τ} = \frac{θ}{x} \sqrt{R e} (η f_{, η} + 2 τ f_{, τ}) \\ = \frac{θ}{x} \sqrt{R e} (3 η - 10 \sqrt{3} τ) (\frac{1}{6} {sech}^{2} ζ + \frac{1}{3} tanh ζ) + \\ + \frac{5 \sqrt{3}}{3} \sqrt{R e} \frac{θ τ}{x} [{sech}^{2} (\frac{5}{6} τ + \frac{1}{2} ln (1 + \sqrt{2}) + tanh (\frac{5}{6} τ + \frac{1}{2} ln (1 + \sqrt{2})] \end{matrix}

(7)

where

ζ = \frac{\sqrt{3}}{6} η - \frac{5}{6} τ - \frac{1}{2} ln (1 + \sqrt{2})

and the Reynolds number is

R e = \frac{| Q |}{ρ θ ν}

. The velocity u and v are depicted in Figure 3 and Figure 4, respectively.

We can see clearly that solution inu and v contain solitary-like wave solutions

{sech}^{2} ζ

and shock wave solutions

tanh ζ

. They compete and suppress each other during different stages of development, eventually forming a solitonlike coherent structure as Lee discovered. The comparing is depicted in Figure 5 below.

The shear stress is given by

\begin{matrix} τ_{x y} & = μ \frac{\partial u}{\partial y} = \frac{ρ ν}{δ} \frac{\partial u}{\partial η} = {(\frac{ν}{ρ})}^{1 / 2} {(\frac{Q}{θ})}^{3 / 2} \frac{1}{x^{2}} f_{, η η} \\ = - \frac{\sqrt{3}}{6} {(\frac{ν}{ρ})}^{1 / 2} {(\frac{| Q |}{θ})}^{3 / 2} \frac{1}{x^{2}} (1 - tanh ζ) {sech}^{2} ζ, \end{matrix}

(8)

The shear stress is depicted in Figure 6.

The wall shear stress

τ_{w} = τ_{x y} |_{y = 0}

, we have

τ_{w} = - \frac{\sqrt{3}}{6} {(\frac{ν}{ρ})}^{1 / 2} {(\frac{| Q |}{θ})}^{3 / 2} \frac{1}{x^{2}} [1 + tanh (c + \frac{5}{6} τ)] {sech}^{2} (c + \frac{5}{6} τ),

(9)

The wall shear stress is depicted in Figure 7.

The expressions of shear stress in Eq.8 indeed indicate the presence of solitonlike coherent structure solutions. Solitary waves are localized disturbances that maintain their shape and amplitude as they propagate, while solitonlike waves are a superposition of different wave components that resemble solitons but may not exhibit all their characteristics. The identification of solitonlike coherent structure (SCS) in boundary layer flow is a significant step towards understanding turbulence origins.

4. Solitonlike coherent structure is a universal motion form of viscous fluid

In the previous section, we demonstrated the existence of solitonlike waves through specific examples. However, the question remains as to whether solittonlike waves are universally present in various types of flows. Analyzing specific problems individually is clearly not a viable approach, and therefore, it is necessary to explore this question theoretically in order to gain a rational understanding and obtain general conclusions.

By introducing a stream function

ψ (x, y, t)

, the velocity components can be expressed as follows

u = \frac{\partial ψ}{\partial y}, v = - \frac{\partial ψ}{\partial x}

. With the relation, the mass conservation is satisfied, and the momentum conservation of 2D boundary layers becomes [31]

\frac{\partial ψ}{\partial t \partial y} + \frac{\partial ψ}{\partial y} \frac{\partial ψ}{\partial x \partial y} - \frac{\partial ψ}{\partial x} \frac{\partial ψ}{\partial y^{2}} = U \frac{d U}{d x} + ν \underset{︸}{\frac{\partial^{3} ψ}{\partial y^{3}}},

(10)

We can easily figure out why the abive equation has solitonlike solution is that the presence of a third-order derivative term in the equation, namely

\frac{\partial^{3} ψ}{\partial y^{3}}

. The presence of the third-order derivative respect to spatial coordinate is the key factors that allows for the existence of soliton-like solutions. This term introduces dissipation in to the nonlinearity of the equation, which is essential for the formation and propagation of solitary waves.

In general, solitons arise in nonlinear systems where the nonlinearity counteracts the dispersion effects. Dispersion tends to spread out waves, causing them to disperse and lose their shape over time. However, nonlinearity can balance this dispersion, leading to the formation of solitary waves that maintain their shape as they propagate.

The specific form of the equation and the coefficients of the terms also play a crucial role in determining the existence and properties of solitons. Different types of equations may have different conditions for the existence of soliton-like solutions. Therefore, it is important to analyze the specific equation and its properties to determine if solitons are present.

Overall, the presence of a third-order derivative term is a key factor that allows for the existence of soliton-like solutions, but other factors such as nonlinearity and specific equation properties also contribute to the formation and behavior of solitary waves.

In the same way, we can confirm that the general 2D, 3D flows and turbulent flows have solitonlike solutions.

A. 2D laminar flows

A general 2D flows governing equations are

\begin{matrix} \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} & = 0, \end{matrix}

(11)

\begin{matrix} \frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} & = - \frac{1}{ρ} p_{, x} + ν (\frac{\partial^{2} u}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}}), \end{matrix}

(12)

\begin{matrix} \frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} & = - \frac{1}{ρ} p_{, y} + ν (\frac{\partial^{2} v}{\partial x^{2}} + \frac{\partial^{2} v}{\partial y^{2}}), \end{matrix}

(13)

Introducing a stream function

ψ (x, y, t)

, the velocity components can expressed as follows

u = ψ_{, y}, v = - ψ_{, x},

(14)

we have a system of partial differential equations

\begin{matrix} ψ_{, y t} + ψ_{, y} ψ_{, y x} - ψ_{x} ψ_{, y y} & = - \frac{1}{ρ} p_{, x} + ν (ψ_{, y x x} + \underset{︸}{ψ_{, y y y}}), \end{matrix}

(15)

\begin{matrix} - ψ_{, x t} - ψ_{, y} ψ_{, x x} + ψ_{x} ψ_{, x y} & = - \frac{1}{ρ} p_{, y} + ν (- \underset{︸}{ψ_{, x x x}} - ψ_{, x y y}) . \end{matrix}

(16)

It is clear to see that the 3rd order derivatives,

ψ_{, x x x}

and

ψ_{, y y y}

, are in the Eq.15 and Eq.16, therefore the general 2D flows have solitonlike solutions.

B. 3D laminar flows

A general 3D laminar flows, whose governing equations are

\begin{matrix} v_{, t} + v \cdot (\nabla v) & = \frac{1}{ρ} \nabla \cdot (p I) + ν \nabla^{2} v, \end{matrix}

(17)

\begin{matrix} \nabla \cdot v & = 0 . \end{matrix}

(18)

where the flow velocity

= v_{i} e_{i}

, the unit tensor

I = e_{i} \otimes e_{i}

, base vector

e_{k}

, the gradient operator

\nabla = e_{k} \frac{\partial}{\partial x^{k}}

, the Laplacian

\nabla^{2} = \nabla \cdot \nabla

Introducing a vectorial function

φ (x, y, z, t) = φ_{1} e_{1} + φ_{2} e_{2} + φ_{3} e_{3} = φ_{k} e_{k}

, and let

v = \nabla \times φ = (e_{i} \partial_{i}) \times (φ_{k} e_{k}) = φ_{k, i} (e_{i} \times e_{k})

, leads to the mass conservation

\nabla \cdot v = \nabla \cdot (\nabla \times φ) = φ_{k, i m} e_{m} \cdot (e_{i} \times e_{k}) = φ_{k, i m} ε_{i k ℓ} e_{m} \cdot e_{ℓ} = φ_{k, i m} ε_{i k ℓ} δ_{m ℓ} = φ_{k, i m} ε_{i k m}

, since

φ_{k, i m} = φ_{k, m i}

, while

ε_{i k m} = - ε_{m k i}

, therefore

\nabla \cdot v = φ_{k, i m} ε_{i k m} \equiv 0

, and momentum equations become

\nabla \times φ_{, t} + (\nabla \times φ) \cdot [\nabla (\nabla \times φ)] = - (1 / ρ) \nabla \cdot (p I) + ν \underset{︸}{\nabla^{2} (\nabla \times φ)} .

(19)

In Cartesian coordinates, if we express the flow velocity components as

u = φ_{3, y} - φ_{2, z}, v = φ_{1, z} - φ_{3, x}, w = φ_{2, x} - φ_{1, y} .

(20)

Eq.19 can be rewritten as follows

\begin{matrix} φ_{3, y t} - φ_{2, z t} + Q_{x} & = - \frac{1}{ρ} p_{, x} + ν \nabla^{2} (φ_{3, y} - φ_{2, z}), \\ φ_{1, z t} - φ_{3, x t} + Q_{y} & = - \frac{1}{ρ} p_{, y} + ν \nabla^{2} (φ_{1, z} - φ_{3, x}), \\ φ_{2, x t} - φ_{1, y t} + Q_{z} & = - \frac{1}{ρ} p_{, y} + ν \nabla^{2} (φ_{2, x} - φ_{1, y}) . \end{matrix}

(21)

where

\begin{matrix} Q_{x} & = (φ_{3, y} - φ_{2, z}) {(φ_{3, y} - φ_{2, z})}_{, x} + (φ_{1, z} - φ_{3, x}) {(φ_{3, y} - φ_{2, z})}_{, y} + (φ_{2, x} - φ_{1, y}) {(φ_{3, y} - φ_{2, z})}_{, z}, \\ Q_{y} & = (φ_{3, y} - φ_{2, z}) {(φ_{1, z} - φ_{3, x})}_{, x} + (φ_{1, z} - φ_{3, x}) {(φ_{1, z} - φ_{3, x})}_{, y} + (φ_{2, x} - φ_{1, y}) {(φ_{1, z} - φ_{3, x})}_{, z}, \\ Q_{z} & = (φ_{3, y} - φ_{2, z}) {(φ_{2, x} - φ_{1, y})}_{, x} + (φ_{1, z} - φ_{3, x}) {(φ_{2, x} - φ_{1, y})}_{, y} + (φ_{2, x} - φ_{1, y}) {(φ_{2, x} - φ_{1, y})}_{, z}, \end{matrix}

(22)

and

\begin{matrix} \nabla^{2} (φ_{3, y} - φ_{2, z}) & = φ_{3, y x x} - φ_{2, z x x} + \underset{︸}{φ_{3, y y y}} - φ_{2, z y y} + φ_{3, y z z} - \underset{︸}{φ_{2, z z z}}, \\ \nabla^{2} (φ_{1, z} - φ_{3, x}) & = φ_{1, z x x} - \underset{︸}{φ_{3, x x x}} + φ_{1, z y y} - φ_{3, x y y} + \underset{︸}{φ_{1, z z z}} - φ_{3, x z z}, \\ \nabla^{2} (φ_{2, x} - φ_{1, y}) & = \underset{︸}{φ_{2, x x x}} - φ_{1, y x x} + φ_{2, x y y} - \underset{︸}{φ_{1, y y y}} + φ_{2, x z z} - φ_{1, y z z} . \end{matrix}

(23)

It is clear to see that the above equations contain the 3rd order differential respect to x, y and z, therefore the general 3D flows also have solitonlike solutions.

C. Turbulent flows

The Reynolds averaged Navier-Stokes equation [33] is

{\bar{v}}_{, t} + \bar{v} \cdot (\nabla \bar{v}) = - (1 / ρ) \nabla \bar{p} + ν \nabla^{2} \bar{v} + \frac{1}{ρ} \nabla \cdot (- ρ \bar{v^{'} \otimes v^{'}}),

(24)

and mass conservation

\nabla \cdot \bar{v} = 0,

(25)

where

\bar{v}

is mean velocity,

\bar{p}

is mean pressure,

v^{'}

is fluctuation velocity, and

- ρ \bar{v^{'} \otimes v^{'}}

is the Reynolds stress tensor.

Introducing a vectorial function

ϕ (x, y, z, t) = ϕ_{k} e_{k} = ϕ_{1} e_{1} + ϕ_{2} e_{2} + ϕ_{3} e_{3}

, in and let

\bar{v} = \nabla \times ϕ

. The mass conservation,

\nabla \cdot \bar{v} = \nabla \cdot (\nabla \times ϕ) \equiv 0

, is satisfied. The momentum equations in Eq.24 become

\begin{matrix} \nabla \times ϕ_{, t} + (\nabla \times ϕ) \cdot [\nabla (\nabla \times ϕ)] = - (1 / ρ) \nabla \bar{p} + ν \underset{︸}{\nabla^{2} (\nabla \times ϕ)} + \frac{1}{ρ} \nabla \cdot (- ρ \bar{v^{'} \otimes v^{'}}) \end{matrix}

(26)

Except the turbulent term

\frac{1}{ρ} \nabla \cdot (- ρ \bar{v^{'} \otimes v^{'}}

, Eq.26 has a same structure as its corresponding 3D laminar flow in Eq.19 and Eq.23, namely, containing the 3rd order derivatives respect to x, y and z, therefore the general 3D turbulent flows also have solitonlike solutions.

5. Conclusions and perspectives

The approach of transforming the 2D unsteady laminar boundary layer equations into a single partial differential equation, as described in Ref. [31], is commendable. This transformation allows for an approximate analytical solution to be obtained for convergent laminar boundary layers.

Through this analysis, it is revealed that the solution for convergent laminar boundary layers includes shock wave solutions and solitary wave solutions. These two types of waves can combine to form solitary-like waves. This analytical prediction of solitary-like waves in boundary layer flow provides evidence for the existence of solitonlike coherent structures in boundary layer flows. It is worth noting that a similar prediction can also be made for turbulent boundary layers.

Furthermore, the study concludes that solitonlike coherent structure exist in all kinds of viscous flows, and reveals the flow motion is particle-SCS motion. However, it is important to emphasize that these coherent structures can only be obtained by combining the Navier-Stokes equations with mass conservation. The inclusion of mass conservation results in the presence of third-order derivatives with respect to the coordinate x. Without the mass conservation condition, the Navier-Stokes equations alone do not possess solitonlike coherent structures.

This research significantly contributes to our understanding of solitonlike coherent structures in boundary layer flows and highlights the importance of considering the combined equations of the Navier-Stokes equations and mass conservation in their study.

Acknowledgments

This work was supported by Xi’an University of Architecture and Technology (Grant No. 002/2040221134). The author would like to express his gratitude to Academician Xiaogang DENG of the Chinese Academy of Sciences and Prof. Cunbiao Lee from Peking University for professional discussions.

Availability of data

The data supporting the findings of this study are available from the corresponding author upon reasonable request.

Declaration of Competing Interest

The author declares that he has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

H. Schlichting, Boundary Layer Theory, 9th (Springer-Verlag Berlin, Heidelberg, 2017).
L. Rosenhead, Laminar Boundary Layers (Dover Publications, New York, 1988).
H. Tennekes and J.L. Lumley, A First Course of Turbulence (The MIT Press, Cambridge, 1972).
U. Frisch, Turbulence (Cambridge University Press, Cambridge,1995).
M. V. Morkovin, On the many faces of transition, in Viscous Drag Reduction, edited by C. S. Wells (Springer, New York, 1969), pp. 1–31.
E. Reshotko, Boundary-layer stability and transition, Annu. Rev. Fluid Mech. 8, 311–349 (1976). [CrossRef]
B. J. Bayly, S. A. Orszag, and T. Herbert, Instability mechanisms in shear-flow transition, Annu. Rev. Fluid Mech. 20, 359–391 (1988). [CrossRef]
I. E. Beckwith and C. G. Miller, Aerothermodynamics and transition in high-speed wind tunnels at NASA Langley, Annu. Rev. Fluid Mech. 22, 419–439 (1990). [CrossRef]
Y. S. Kachanov, Physical mechanisms of laminar-boundary-layer transition, Annu. Rev. Fluid Mech. 26, 411–482 (1994). [CrossRef]
H. L. Reed, W. S. Saric, and D. Arnal, Linear stability theory applied to boundary layers, Annu. Rev. Fluid Mech. 28, 389–428 (1996). [CrossRef]
T. Herbert, “Parabolized stability equations,” Annu. Rev. Fluid Mech. 29, 245–283 (1997). [CrossRef]
C. B. Lee, Possible universal transitional scenario in a flat plate boundary layer: Measurement and visualization, Phys. Rev. E 62, 3659 (2000). [CrossRef]
W. S. Saric, H. L. Reed, and E. B. White, Stability and transition of threedimensional boundary layers, Annu. Rev. Fluid Mech. 35, 413–440 (2003). [CrossRef]
P. Durbin and X. Wu, Transition beneath vortical disturbances, Annu. Rev. Fluid Mech. 39, 107–128 (2007). [CrossRef]
B. Eckhardt, T. M. Schneider, B. Hof, and J. Westerweel, Turbulence transition in pipe flow, Annu. Rev. Fluid Mech. 39, 447–468 (2007).
C. B. Lee and J. Z. Wu, Transition in wall-bounded flows, Appl. Mech. Rev. 61, 030802 (2008).
A. Fedorov, Transition and stability of high-speed boundary layers, Annu. Rev. Fluid Mech. 43, 79–95 (2011). [CrossRef]
X. Zhong and X. Wang, Direct numerical simulation on the receptivity, instability, and transition of hypersonic boundary layers, Annu. Rev. Fluid Mech. 44, 527–561 (2012). [CrossRef]
J. Jiménez, Coherent structures in wall-bounded turbulence, J. Fluid Mech. 842, 513–531 (2018).
X. Wu, Nonlinear theories for shear flow instabilities: Physical insights and practical implications, Annu. Rev. Fluid Mech. 51, 451–485 (2019). [CrossRef]
C.B. Lee and X.Y. Jiang, Flow structures in transitional and turbulent boundary layer, Phys. Fluids 31, 111301 (2019). [CrossRef]
X. Y. Jiang, C. B. Lee, X. Chen, C. R. Smith and P. F. Linden, Structure evolution at early stage of boundary layer transition: simulation and experiment, J. Fluid Mech., 890, A11(2020).
B.H. Sun, Turbulent Poiseuille flow modeling by modified Prandtl-van Driest mixing length, Acta Mech. Sin., Vol. 39, 322066 (2023). [CrossRef]
K. Stewartson, On the impulsive motion of a flat plate in a viscous fluid, Quart. J. Mech. Appl. Math., 4, 182-198(1951).
K. Stewartson, On the impulsive motion of a flate plate in a viscous fluid, II, Quart. J. Mech. Appl. Math., 26 (1973), pp. 143-152.
K. Takuda, On the impulsive motion of a flat plate in a viscous fluid, J. Fluid Mech. 33 (4), 657-675(1968).
M. G. Hall, The boundary layer over an impulsively started flat plate, Proc. Roy. Soc. Ser. A, 310, 401-414(1969). [CrossRef]
S. C. R. Dennis, Motion of a viscous.fluid past an impulsively started semi-infinite flat plate, J. Inst. Math. Appl., 10, 105-117 (1972).
K. Stewartson, The theory of unsteady laminar boundary layers, Advances in Applied Mechanics, vol. VI, Academic Press, New York, 1-37(1960).
N. Riley, Unsteady Laminar Boundary Layers, SIAM Review, 17(4), 1975.
B.H. Sun, Exact similarity solutions of unsteady laminar boundary layer flows, Preprints 2023, 2023092117. [CrossRef]
S.B. Pope, Turbulent Flows (Cambridge University, Cambridge, 2000).
K. Pohlhausen, Zur näherungsweisen Integration der Differentialgleichung der laminaren Grenzschicht, Z.Angew.Math.Mech.,1,252-268. [CrossRef]
L.D. Landau and E.M. Lifshitz, Fluid Mechanics (2nd edtion) (Elservier, Butterworth-Heinemann,1976).

Figure 1.

f_{, η}

Figure 1.

f_{, η}

Figure 2. Region of unsteady laminar boundary layers past a flat plate: Reylaigh solution and Blasius solution.

Figure 3.

\frac{x u}{ν R e} = \frac{1}{2} {sech}^{2} ζ + tanh ζ

Figure 3.

\frac{x u}{ν R e} = \frac{1}{2} {sech}^{2} ζ + tanh ζ

Figure 4.

\frac{x v}{θ \sqrt{R e}}

Figure 4.

\frac{x v}{θ \sqrt{R e}}

Figure 5. Red line:

{sech}^{2} ζ

, blue line:

tanh ζ

Figure 5. Red line:

{sech}^{2} ζ

, blue line:

tanh ζ

Figure 6.

K_{x y} = x^{2} {(\frac{ρ}{ν})}^{1 / 2} {(\frac{θ}{| Q |})}^{3 / 2} τ_{x y}

Figure 6.

K_{x y} = x^{2} {(\frac{ρ}{ν})}^{1 / 2} {(\frac{θ}{| Q |})}^{3 / 2} τ_{x y}

Figure 7.

K_{w} = x^{2} {(\frac{ρ}{ν})}^{1 / 2} {(\frac{θ}{| Q |})}^{3 / 2} τ_{w}

Figure 7.

K_{w} = x^{2} {(\frac{ρ}{ν})}^{1 / 2} {(\frac{θ}{| Q |})}^{3 / 2} τ_{w}

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Solitonlike Coherent Structure Is a Universal Motion Form of Viscous Fluid

Abstract

1. Introduction

2. A thin plate immersed in unsteady convergent flow

3. Analytic solitonlike solution

4. Solitonlike coherent structure is a universal motion form of viscous fluid

5. Conclusions and perspectives

Acknowledgments

Availability of data

Declaration of Competing Interest

References

MDPI Initiatives

Important Links

Subscribe