Next Article in Journal
Revisiting Angular Momentum Conservation in Transport Simulations of Intermediate-Energy Heavy-Ion Collisions
Previous Article in Journal
How to Strengthen Constraints on Non-Newtonian Gravity from Measuring the Lateral Casimir Force
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Forms of Solutions to System of Nonlinear Partial Differential Equations

by
Alexander D. Bruno
1,*,† and
Alexander B. Batkhin
1,2,†
1
Keldysh Institute of Applied Mathematics of RAS, sq. Miusskaya 4, 125047 Moscow, Russia
2
Moscow Institute of Physics and Technology, National Research University, 9 Institutskiy per., 141701 Dolgoprudny, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2023, 9(1), 35; https://doi.org/10.3390/universe9010035
Submission received: 29 September 2022 / Revised: 14 December 2022 / Accepted: 29 December 2022 / Published: 3 January 2023
(This article belongs to the Section Mathematical Physics)

Abstract

:
Here, we considerably develop the methods of power geometry for a system of partial differential equations and apply them to two different fluid dynamics problems: computing the boundary layer on a needle in the first approximation and computing the asymptotic forms of solutions to the problem of evolution of the turbulent flow. For each equation of the system, its Newton polyhedron and its hyperfaces with their normals and truncated equations are calculated. To simplify the truncated systems, power-logarithmic transformations are used and the truncated systems are further extracted. Here, we propose algorithms for computing unimodular matrices of power transformations for differential equations. Results: (1) the boundary layer on the needle is absent in liquid, while in gas it is described in the first approximation; (2) the solutions to the problem of evolution of turbulent flow have eight asymptotic forms, presented explicitly.

1. Introduction

A universal asymptotic nonlinear analysis is formed, whose unified methods allow finding asymptotic forms and expansions of solutions to nonlinear equations and systems of different types:
  • Algebraic;
  • Ordinary differential equations (ODEs);
  • Partial differential equations (PDEs).
This calculus contains two methods:
1
Transformation of coordinates, bringing equations to normal form;
2
Separating truncated equations.
Two kinds of coordinate changes can be used to analyze the resulting equations:
1
Power;
2
Logarithmic.
In this paper, we consider systems of nonlinear partial differential equations in two variants:
1
With boundary conditions;
2
Without boundary conditions.
We show how to find asymptotic forms of their solutions using algorithms of power geometry. In this case, by asymptotic form of solution, we mean a simple expression in which each of the independent or dependent variables tends to zero or infinity.
Here, we consider two fluids problems: (1) boundary layer and (2) turbulence flow by methods of power geometry.
For problem (1), it was firstly given in [1] (Chapter 6, Section 6). The usual approach was in papers [2] and [3]; see also [4,5]. For the new approach via power geometry, see [6], and here in Section 3. A boundary layer on a needle has a stronger singularity than on a plane, and it was first considered in [6]. We are not sure that it is possible with the usual analysis. Our approach is, in some sense, opposite to the approach in [7].
For problem (2), we firstly make it here and we are not sure that it is possible with the usual analysis. Our approach is, in a sense, opposite to the approach in [8].
The structure of the paper is as follows. Section 2 outlines the basics of power geometry for partial differential equations. These are applied in Section 3 to calculate the boundary layer on the needle. In Section 4, the theory and algorithms are further developed to apply to variant 2 problems. In Section 5, Section 6 and Section 7, they are used to compute asymptotic forms of evolution of turbulent flow. Section 8 contains a summary of the computed asymptotics in the above sections.

2. Basics of Power Geometry

For more detail, see [1] (Chapters VI–VIII).
Let X = ( x 1 , , x m ) C m be independent and Y = ( y 1 , , y n ) C n be dependent variables. Place Z = ( X , Y ) C n + m . Differential monomial a ( Z ) is a product of an ordinary monomial c Z R = c z 1 r 1 z m + n r m + n , where c = const C , and a finite number of derivatives of the form
l y j x 1 l 1 l m x m l y j X L , l j 0 , j = 1 m l j = l , L = l 1 , , l m .
The differential monomial a ( Z ) corresponds to its vector exponent of degree Q ( a ) R m + n , formed by the following rules:
Q Z R = R , Q l y j / X L = ( L , E j ) ,
where E j is the unit vector. The product of monomials corresponds to the sum of their vector exponents of degree:
Q ( a b ) = Q ( a ) + Q ( b ) .
Differential sum is the sum of differential monomials:
f ( Z ) = a k ( Z ) .
The set S ( f ) of vector exponents Q ( a k ) is called support of sum f ( Z ) . The closure of the convex hull
Γ ( f ) = Q = λ j Q j , Q j S , λ j 0 , λ j = 1
of the support S ( f ) is called the polyhedron of the sum f ( Z ) . The boundary Γ of the polyhedron Γ ( f ) consists of generalized faces Γ j ( d ) , where d = dim Γ j ( d ) , 0 d m + n 1 . Each face Γ j ( d ) corresponds to:
  • Normal cone
    U j ( d ) = P R * m + n : P , Q = P , Q > P , Q , where Q , Q Γ j ( d ) , Q Γ \ Γ j ( d ) ,
    where the space R * m + n is conjugate to the space R m + n , · , · is a scalar product;
  • Truncated sum
    f ^ j ( d ) ( Z ) = a k ( Z ) over Q ( a k ) Γ j ( d ) S .
Consider a system of equations:
f i ( X , Y ) = 0 , i = 1 , , n ,
where f i are differential sums. Each equation f i = 0 corresponds to:
  • Its support S ( f i ) ;
  • Its polyhedron Γ ( f i ) with a set of faces Γ i j d i in the main space R m + n ;
  • Set of their normal cones U i j d i in the dual space R * m + n ;
  • Set of truncated equations f ^ i j d i ( X , Y ) = 0 .
The set of truncated equations
f ^ i j i d i ( X , Y ) = 0 , i = 1 , , n ,
is a truncated system if the intersection is not empty:
U 1 j i d 1 U n j n d n .
A truncated system is always a quasi-homogeneous system.
In the solution of the system (4),
y i = φ i ( X ) , i = 1 , , n ,
where φ i are series in powers of x k and their logarithms, each φ i corresponds to its support, polyhedron, normal cones u i , and truncations. Here, the logarithm ln x i has a zero exponent of degree on x i .
The set of truncated solutions y i = φ ^ i , i = 1 , , n , corresponds to the intersection of their normal cones:
u = i = 1 n u i R * m + n .
If it is not empty, it corresponds to truncated solution:
y i = φ ^ i , i = 1 , , n .
Theorem 1.
If the normal cone u intersects the normal cone (6), then the truncation y i = φ ^ ( X ) , i = 1 , , n , of this solution satisfies the truncated system (5).
Multiplying the differential sum (5) with the support S ( f ) by the monomial Z R gives the differential sum, g ( Z ) = Z R f ( Z ) , with the support S ( g ) = R + S ( f ) . Thus, the multiplication leads to a shift of supports. Multiplications by monomials form a group of linear transformations of supports, and they can be used to simplify supports, differential sums, and systems of equations.
Let S ( f ) be the support of the differential sum f ( Z ) and Q S ( f ) . The set
S ˜ ( f ) S ( f ) Q
is called shifted support of the sum f ( Z ) .
Each equation f i = 0 in the system (4) corresponds to a support S ( f i ) and a shifted support S ˜ ( f i ) . Let Γ ˜ be the convex hull of their union
S ˜ ( f 1 ) S ˜ ( f n ) ,
and d is the dimension of Γ ˜ . If  d < m + n , then the system (4) is quasi-homogeneous.
A similar technique is valid for systems of equations containing small or large parameters. Here, the exponents of degrees of these parameters are taken into account in the same way as the exponents of degrees of variables tending to zero or to infinity.

3. Boundary Layer on a Needle

The theory of the boundary layer on a plate for the flow of a viscous incompressible fluid was developed by Prandtl ([2], 1904) and Blasius ([3], 1908) (see [1] (Chapter 6, Section 6)). More developed their theory; see [4,5]. However, a similar theory for the boundary layer on a needle was not available until recently, for the sticking conditions on a needle correspond to a stronger singularity than on a plane. This theory has been constructed using power geometry [6].
Let there be an axis x in three-dimensional space, r is the distance from it, and a semi-infinite needle located on the semiaxis x > 0 , r = 0 . Stationary axisymmetric viscous fluid flows were studied, which at x = have a constant velocity parallel to the axis x, and on the needle satisfy the sticking conditions (Figure 1). Two variants were considered.
The first variant: an incompressible fluid. For this, the Navier–Stokes equations in independent variables x , r are equivalent to a system of two PDEs for the flow function ψ and pressure p:
g 1 1 r ψ x r 1 r ψ r + 1 r ψ r x 1 r ψ r + 1 ρ p x ν 1 r r r r 1 r ψ r + 2 x 2 1 r ψ r = 0 , g 2 1 r ψ x r 1 r ψ x 1 r ψ r x 1 r ψ x + 1 ρ p r + ν r 1 r 2 ψ x r + 2 x 2 1 r ψ x = 0 ,
where density ρ and viscosity ν = const , with boundary conditions
ψ = u 2 r 2 , p = p 0 at x = , u , p 0 = const ;
ψ x = ψ r = 2 ψ x r = 2 ψ r 2 = 0 in x 0 , r = 0 .
The system (8) has the form (4) with m = n = 2 and m + n = 4 . Thus, supports of Equation (8) should be considered in R 4 . It turns out that the polyhedra Γ ( g 1 ) and Γ ( g 2 ) of the Equation (8) are three-dimensional tetrahedrons which can be placed into one linear three-dimensional subspace by parallel transfer, which simplifies the separation of truncated systems. Analyzing the solutions of the truncated systems and the results of their jointing, it was possible to show that the system (8) has no solution with p 0 satisfying both boundary conditions (9) and (10).
The second variant: a compressible thermally conductive fluid and a nonthermally conductive needle. For this variant, the Navier–Stokes equations in independent variables x , r are equivalent to a system of three PDEs for the flow function ψ , density ρ , and enthalpy h (analog of temperature):
f 1 1 r ψ x r 1 ρ r ψ x + 1 r ψ r x 1 ρ r ψ x A r ( ρ h ) + + 2 3 C ν r h ν r r 1 ρ ψ x 2 3 C ν r h ν r x 1 ρ ψ r 2 C ν r r h ν r r 1 ρ r ψ x + C ν x h ν r 1 ρ r ψ r C ν x h ν x 1 ρ r ψ x + 2 C ν h ν ρ r 3 ψ x = 0 , f 2 1 r ψ x r 1 ρ r ψ r 1 r ψ r x 1 ρ r ψ r A x ( ρ h ) + + 2 3 C ν x h ν r r 1 ρ ψ x 2 3 C ν x h ν r x 1 ρ ψ r + + C ν r r h ν r r 1 ρ r ψ r C ν r r h ν r x 1 ρ r ψ x + + 2 C ν x h ν x 1 ρ r ψ r = 0 , f 3 1 r ψ x h r 1 r ψ r h x A ρ r ψ x ( ρ h ) r + A ρ r ψ r ( ρ h ) x + + 2 C ν h ν r 1 ρ r ψ x 2 + 2 C ν h ν 1 r 2 ρ ψ x 2 + 2 C ν h ν x 1 ρ r ψ r 2 + + C ν h ν x 1 ρ r ψ x 2 C ν h ν x 1 ρ r ψ x r 1 ρ r ψ r + + C ν h ν r 1 ρ r ψ r 2 2 3 C ν h ν 1 r r 1 ρ ψ x 2 + + 4 C ν h ν 3 r r 1 ρ ψ x x 1 ρ r ψ r 2 3 C ν h ν x 1 ρ r ψ r 2 + + C ν σ r r r h ν h r + C ν σ x h ν h x = 0 ,
where the parameters A , C , σ > 0 and ν [ 0 , 1 ] , with boundary conditions
ψ = ψ 0 r 2 , ρ = ρ 0 , h = h 0 at x = , ψ 0 , ρ 0 , h 0 = const
and (10). Here, x 1 = x , x 2 = r , y 1 = ψ , y 2 = ρ , y 3 = h , so m = 2 , n = 3 , m + n = 5 . In the space R 5 , all polyhedra Γ ( f 1 ) , Γ ( f 2 ) , Γ ( f 3 ) of equations (11) are three-dimensional and can be shifted parallel in one linear three-dimensional subspace. In the coordinates q ˜ 1 , q ˜ 2 , q ˜ 3 of this three-dimensional space, they are shown in Figure 2, Figure 3 and Figure 4, respectively.
It follows from the boundary conditions (12) that the boundary layer corresponds to a normal vector P = ( 2 , 1 , 2 , 0 , 0 ) . Thus, the truncated system corresponding to the boundary layer on the needle has the form:
f ^ 12 ( 0 ) A ( ρ h ) r = 0 or ( ρ h ) r = 0 , f ^ 22 ( 2 ) 1 r ψ x r 1 ρ r ψ r 1 r ψ r x 1 ρ r ψ r A ( ρ h ) x + C ν r r h ν r r 1 ρ r ψ r = 0 , f ^ 32 ( 2 ) 1 r ψ x h r 1 r ψ r h x A ρ r ψ x ( ρ h ) r + A ρ r ψ r ( ρ h ) x + C ν h ν r 1 ρ r ψ r 2 + + C ν σ r r r h ν h r = 0 ,
with self-similar variables
ψ = x G ( ξ ) , ρ = P ( ξ ) , h = H ( ξ ) , ξ = r 2 / x ,
and the boundary conditions
ψ = ψ 0 r 2 , ρ = ρ 0 , h = h 0 ; ψ 0 , ρ 0 , h 0 = const , r
and (10). In Figure 2, Figure 3 and Figure 4, the vertex and faces corresponding to the truncated system (13) are boldfaced. According to Equations  (13)–(15), the product P ( ξ ) H ( ξ ) = const = C 0 ρ 0 h 0 . Thus, P ( ξ ) = C 0 / H ( ξ ) and the system (13), for variables (14) is equivalent to a system of two ODEs:
F 2 G G H + 2 C ν ξ H ν G H = 0 , F 3 2 G H + 16 C ν C 0 2 ξ H ν G H 2 + 4 C ν σ 1 ξ H ν H = 0
where d / d ξ , with boundary conditions
G = ψ 0 ξ , H = h 0 for ξ + ,
G = d G / d ξ = 0 for ξ = 0 .
The problem (16), (17) and (18) has an invariant manifold G H = 0 , on which it reduces to one equation
Δ 2 ξ H ν H H 2 ξ H ν H 2 + ( ξ + c 2 ) H = 0 ,
where c 2 is an arbitrary constant, with boundary conditions
H 1 at ξ + , H + at ξ + 0 .
An analysis of the solutions of the last problem by methods of planar power geometry [9] shows that for ν ( 0 , 1 ) , it has solutions of the form
H c 3 | ln ξ 1 / ν , ξ 0 ,
where c 3 is an arbitrary constant.
Thus, at  ν ( 0 , 1 ) in the boundary layer r 2 / x < const at x + and ξ = r 2 / x 0 , the asymptotic form of the flow is obtained:
ψ c 1 r 2 ln ξ 1 / ν , ρ c 2 ln ξ 1 / ν , h c 3 ln ξ 1 / ν
i.e., near the needle, the density decreases to zero and the temperature increases to infinity as the distance from the tip of the needle tends to plus infinity.

4. Algorithms of Power Geometry

4.1. Euler’s Algorithm and a Generalization of Continued Fraction

A matrix α is called unimodular if all its elements are integer and det α = ± 1 .
Problem 1.
Let n-dimensional integer vector A = ( a 1 , a 2 , , a n ) be given. Find an n-dimensional unimodular matrix α such that the vector A α = C = ( c 1 , , c n ) contains only one coordinate c n different from zero.
To solve it, Euler [10] proposed the following algorithm. Firstly, let all coordinates of the vector A be non-negative. Using the permutation A α 0 = a ˜ 1 , a ˜ 2 , , a ˜ n , we order the coordinates
a ˜ j a ˜ j + 1 , j = 1 , , n 1 .
Here, α 0 is a unimodular permutation matrix. Let a ˜ k be the smallest of the numbers a ˜ j different from zero.
Let
b j = a ˜ j a ˜ k , j = 1 , , n ,
where [ x ] is the integer part of the number x. In this case, b 1 = = b k 1 = 0 , b k = 1 . Let us perform the transformation
d j = a ˜ j b j a ˜ k , 1 j n , j k , d k = α ˜ k .
It corresponds to a unimodular matrix α 1 that has ones on the diagonal, and in the kth row are elements
0 , 0 , , 0 , 1 , b k + 1 , , b n ,
i.e.,
A ˜ α 1 = D = ( d 1 , , d n ) .
Now we order the components of the vector D using the unimodular permutation matrix β 0 such that D β 0 = D ˜ = 0 , , 0 , d ˜ k , , d ˜ n , where d ˜ j d ˜ j + 1 .
Let d ˜ l be the smallest of d ˜ j different from zero, and  e j = d ˜ j / d ˜ l , j = 1 , , l . We perform the transformation
f j = d ˜ j e j d ˜ l , 1 j n , j l , f l = d ˜ l ,
and so on. At each step, the maximum of the coordinates of the vector decreases and is the nth coordinate. Thus, after a finite number of steps we obtain a vector with one nonzero coordinate which is the last one. Its value is the GCD of all initial coordinates a 1 , , a n . Each step consists of a permutation matrix and a triangular matrix with a unit diagonal:
A α 0 α 1 β 0 β 1 γ 0 γ 1 ω 0 ω 1 = A α = C = ( 0 , , 0 , c n ) .
Matrix
α = α 0 α 1 β 0 β 1 γ 0 γ 1 ω 0 ω 1
is the solution to Problem 1.
If not all coordinates a j of the original vector A are of the same sign, then we first order them by modulo
a ˜ j a ˜ j + 1
and suppose
b j = [ a ˜ j / a ˜ k ] sign a ˜ j sign a k ˜ .
Remark 1.
By multiplying the matrix α on the right by a unimodular permutation matrix, we can obtain a vector from vector C that has all but one coordinate equal to zero, and a single nonzero coordinate located at any position.

4.2. Power Transformations

To simplify a truncated system (5) and any quasi-homogeneous system, it is convenient to use a power transformation. Let α be a square real nondegenerate block matrix of dimension m + n of the form
α = α 11 α 12 0 α 22 ,
where α 11 and α 22 are square matrices of sizes m and n, respectively. We denote ln Z = ln z 1 , , ln z m + n , and by the asterisk * we denote transposition.
Transformation of the variables
ln W = ln Z α
is called power transformation.
Theorem 2
([1]). The power transformation (22) changes a differential monomial a ( Z ) with exponent of degree Q ( a ) into a differential sum b ( W ) with exponent of degree Q ( b ) :
R = Q ( b ) = Q ( a ) α 1 * .
Corollary 1.
The power transformation (22) changes the differential sum (3) with support S ( f ) into the differential sum g ( W ) with support
S ( g ) = S ( f ) α 1 * ,
i.e.,
S ( f ) = S ( g ) α *
Theorem 3
([1]). If the system (4) is a quasi-homogeneous system and d = dim Γ ˜ , then there exist a power transformation (22) and monomials Z T i , i = 1 , , n which change the system (4) into the system
g i ( W ) Z T i f i ( Z ) = 0 , i = 1 , , n ,
where all g i ( W ) are differential sums, and all their supports S ( g i ) have m + n d identical coordinates q j equal to zero.
Usually, the supports of differential equations are integer. For them, it is desirable to have power transformations that preserve the integrability of the supports. This property is possessed by power transformations (22) with a unimodular matrix α in which all elements are integers and det α = ± 1 . For a unimodular matrix α , its inverse α 1 and transpose α * matrices are also unimodular.
Let us compute the unimodular matrix (21) of the power transformation (22) in one important case. Suppose that in the system (4) all supports S ( f 1 ) , , S ( f n ) are integers and the normal to them is an integer vector N = ( ν 1 , , ν m + n ) 0 , i.e., for all Q S ( f i ) we have Q , N = λ i , i = 1 , , n . Split the vector N into two parts: N 1 = ( ν 1 , , ν m ) and N 2 = ( ν m + 1 , , ν m + n ) , and perform the same for the vector Q = ( Q 1 , Q 2 ) .
Consider three cases:
(1)
N 1 = 0 , N 2 0 , then Q 2 , N 2 = λ i ;
(2)
N 1 0 , N 2 = 0 , then Q 1 , N 1 = λ i ;
(3)
N 1 0 , N 2 0 .
Below, I is the unit matrix.
Lemma 1.
In the case 1, there exists a unimodular matrix α 22 of size n such that, after transformation (22) with α 11 = I , α 12 = 0 in each transformed differential sum g i ( W ) , the coordinate w m + n is contained only in a fixed degree λ ˜ i .
Proof of Lemma 1.
Using the Euler algorithm from Section 4.1 for the vector N 2 , we find such a unimodular matrix α 22 of size n that
N 2 α 22 = ( 0 , , 0 , μ 2 ) ,
and μ 2 is GCD of numbers ν m + 1 , , ν m + n . According to (23) and (24) for Q S ( f i ) we have
λ i = N 2 , Q 2 = N 2 , R 2 α 22 * = N 2 α 22 , R 2 = μ 2 r m + n .
Then, λ ˜ i = λ i / μ 2 . The proof is over.    □
Lemma 2.
In the case 2, there exists a unimodular matrix α 11 of size m such that after transformation (22) with α 22 = I , α 21 = 0 in each transformed differential sum g i ( W ) , the coordinate w m contains only a fixed degree λ ˜ i .
The proof is the same as the proof of Lemma 1.
Lemma 3.
In the case 3, if  gcd N 2 / gcd N 1 = ω is an integer, then there exists such a unimodular matrix α of (21) that every differential sum g i ( W ) contains the coordinate w m only in a fixed degree λ ˜ i .
Proof of Lemma 3.
Let
μ i = gcd N i , i = 1 , 2 .
By Euler’s algorithm, we obtain the representations
N i β i = ( 0 , , 0 , μ i ) , i = 1 , 2 ,
where β i are unimodular matrices of sizes m and n, respectively. In other words,
N β = ( 0 , , μ 1 , 0 , , μ 2 ) M = μ 1 ( 0 , , 1 , 0 , , ω ) ,
where β is a block unimodular matrix
β = β 1 0 0 β 2
Then, we have
M ( I γ ) = ( 0 , , μ 1 , 0 , , 0 ) ,
where I is a unit matrix of size m + n , and the matrix γ = ( γ i j ) has a single nonzero element γ m , m + n = ω . Then the matrix α = β ( I γ ) is unimodular, has a block structure (21), and each differential sum g i ( W ) contains the coordinate w m in degree λ ˜ = λ i / μ 1 . Reducing each of them by the value of w m λ ˜ i , we obtain a system in which the variable w m is contained with zero-degree exponent. The proof is over.    □
Remark 2.
If the relation ω is not integer, we can still perform a degree transformation of Lemma 3, but the support of the transformed system will not be integer.

4.3. Logarithmic Transformation

Let z j be one of the coordinates x k or y l according to the beginning of Section 2. Transformation:
ζ j = ln z j
Let us call this logarithmic transformation.
Theorem 4
([11]). Let f ( Z ) be such a differential sum that for all its monomials, jth component of q j vector degree exponent Q = ( q 1 , , q m + n ) is zero, then as a result of the logarithmic transformation (25), a differential sum f ( Z ) transforms into a differential sum from z 1 , , ζ j , , z n .
In the system
f i ( X , Y ) = 0 , i = 1 , , n ,
let all f i be differential sums. Let some of its truncated system be
f ^ i ( X , Y ) = 0 , i = 1 , , n .
It is quasi-homogeneous in dimension d < m + n . According to Theorem 3 there exists a power transformation (22) which reduces the system (27) to the system
g i ( W ) = 0 , i = 1 , , n ,
in which all supports of sums g i ( W ) have m + n d zero coordinates. A logarithmic transformation can be applied to these coordinates, which by theorem 4 will reduce the system (28) to the form
h i W ˜ = 0 , i = 1 , , n ,
where h i are differential sums, and  w ˜ j = w j or ln w j , j = 1 , , m + n . In the system (29) we can again select truncated systems and so on.
For z j 0 or , the coordinate ζ j = ln w j always tends to ± . If we are interested only in those solutions (7) which have a normal cone u intersecting a given cone K, then the cone K is called the cone of problem. Thus, after the logarithmic transformation (25) for the coordinate ζ j in the cone of the problem, we have p j 0 .
In the following, we will not consider all possible truncated systems (5), but only those in which one of the equations has dimension d i = m + n 1 . The calculations show that in this case the above procedure will cover all the truncated systems. Finally, it is convenient to combine the power and logarithmic transformations. Namely, the logarithmic transformation is performed for the coordinate w m + n in the case 1 and for the coordinate w m in the cases 2 and 3 of Section 4.2.

4.4. System of Notations

The original system is denoted by S, and its equations by E 1 S and E 2 S , respectively. For the equations of the system S, the polyhedrons, normal cones, are calculated and the corresponding shortened systems are found by them, which are denoted as S ( 1 ) , S ( 2 ) , etc. For the truncated system S ( k ) , a power and/or logarithmic transformation is applied, the result of which is the system P ( k ) . The corresponding truncations of the system P ( k ) are denoted by S ( k , 1 ) , S ( k , 2 ) , etc., and the results of their power-logarithmic transformations are denoted by P ( k , 1 ) , P ( k , 2 ) , etc. If new truncations are required, the corresponding systems are denoted as S ( k , l , m ) , and the results of the power-logarithmic transformations are denoted as P ( k , l , m ) . This branching procedure stops when one obtains a system that is solvable explicitly. Each system S ( m ) has its cone of problem K [ S ( m ) ] . In the following Section 5, Section 6 and Section 7, the vectors are denoted in square brackets [ x 1 , , x m ] , as is usual in Maple.

4.5. About the Computation of the Objects of Power Geometry

The computer algebra system Maple 2021 [12] was used for calculations in this work. A library of procedures based on the PolyhedralSets CAS Maple package was developed to implement the algorithms of power geometry. The library includes calculation procedures:
  • Vector degree exponent Q of the differential monomial a ( Z ) for a given order of independent and dependent variables.
  • Of the support S of a partial differential equation written as a sum of differential monomials.
  • Newton’s polyhedron Γ in the form of a graph of generalized faces Γ j ( d ) of all dimensions d for the given support of the equation (see below Figure 5 and Figure 6); the number j is given by the program; each generalized face has its own number j; each line of the graph contains all generalized faces Γ j ( d ) of the same dimension d, the first line contains the Newton’s polyhedron Γ , the next line contains all faces Γ j ( m + n 1 ) of dimension m + n 1 and so on; the last line contains the empty set; if Γ j ( d ) Γ k ( d + 1 ) , then they are connected by an arrow. In ([1], Ch. 1, Section 1), “the structural diagram” was used that is similar to the graph and differs from it in two properties: numeration of faces Γ j ( d ) is independent for each dimension d and arrows are replaced by segments (see also [13]).
  • Of the normal vector N j for the each generalized face Γ j ( m + n 1 ) for the second line of the graph;
  • Of the truncated equation f ^ j ( d ) = 0 by the given number j of the generalized face.
  • Of the truncated equation f ^ j ( d ) = 0 by a given normal vector N j , if d = m + n 1 .
  • Of the normal cone of the corresponding generalized face: if the face
    Γ j ( d ) = Γ i ( m + n 1 ) Γ k ( m + n 1 ) Γ l ( m + n 1 ) ,
    then the normal cone U j ( d ) is the conic hull of the normals N i , N k , , N l .
  • To calculate the power or logarithmic transformation of the original variables by a given normal N of the hyperface. For this purpose, the algorithms for constructing the unimodular matrix described in Section 4.1 are used.

5. The k ε Model of Evolution of Turbulent Bursts

According to [14,15,16,17], the model is described by the system
k t = k 2 ε k x x ε , ε t = k 2 ε ε x x γ ε 2 k .
Here, time t and coordinate x are independent variables, the turbulent density k and the dissipation rate ε are dependent variables, and γ is a real parameter. Here, m = n = 2 , m + n = 4 and x 1 = t , x 2 = x , y 1 = k , y 2 = ε .
The support of the first equation S 1 of the system (30) consists of points
Q 1 = [ 1 , 0 , 1 , 0 ] , Q 2 = [ 0 , 2 , 3 , 1 ] , Q 3 = [ 0 , 0 , 0 , 1 ] .
The support of the second equation S 2 of the system (30) consists of points
Q 4 = [ 1 , 0 , 0 , 1 ] , Q 5 = [ 0 , 2 , 2 , 0 ] , Q 6 = [ 0 , 0 , 1 , 2 ] .
The shifted supports S ˜ 1 = S 1 Q 3 and S ˜ 2 = S 2 Q 6 consist of three points:
R 1 Q 1 Q 3 = Q 4 Q 6 = [ 1 , 0 , 1 , 1 ] , R 2 Q 2 Q 3 = Q 5 Q 6 = [ 0 , 2 , 3 , 2 ] , 0 = Q 3 Q 3 = Q 6 Q 6 .
Therefore, d = 2 .
According to Theorem 3 let us introduce new dependent variables:
u = Z R 1 = t 1 k ε 1 , v = Z R 2 = x 2 k 3 ε 2 .
Then
k = x 2 v t 2 u 2 , ε = x 2 v t 3 u 3 .
This is a power transformation (22) with matrix (21), where
α 11 = 1 0 0 1 , α 12 = 1 0 0 2 , α 22 = 1 3 1 2 .
This power transformation is constructed directly on the support of the system such that it lies in the coordinate plane. The theory of Section 4.2 is not used here.
Change of the variables (31) leads the system (30) to the form
u t ( ln v ) t 2 u 2 t u t = v 6 12 U + 7 V + 6 U 2 7 U V + 2 V 2 2 x U x + x V x 1 , u t ( ln v ) t 3 u 3 t u t = v 6 17 U + 7 V + 12 U 2 10 U V + 2 V 2 3 x U x + x V x γ ,
where U = x ln u x , V = x ln v x .
Let us find the self-similar solutions of this system. Consider two cases.
The first case: u , v are constants. Then, the system (32) has the form
2 u = 6 v 1 , 3 u = 6 v γ .
Its solution
u = γ 1 , v = ( 3 2 γ ) / 6
has two critical values: γ = 1 and γ = 3 / 2 .
The second case: Let ζ = t σ x , where σ R . Now u and v are functions of ζ . In this case, in the matrix (21), the submatrix
α 11 = 1 σ 0 1
and the submatrices α 21 and α 22 are the same as before. For u ( ζ ) and v ( ζ ) , the system (32) after substitutions
t = d d ζ ζ t = σ ζ t d d ζ , x = d d ζ ζ x = ζ x d d ζ
generates a one-parameter by σ family of systems of two ODEs:
σ u ζ ( ln v ) ζ 2 u 2 u ζ u ζ = v 6 12 U + 7 V + 6 U 2 7 U V + 2 V 2 2 ζ U ζ + ζ V ζ 1 , σ u ζ ( ln v ) ζ 3 u 3 ζ u ζ = v 6 17 U + 7 V + 12 U 2 10 U V + 2 V 2 3 ζ U ζ + ζ V ζ γ ,
where U = ζ ln u ζ , V = ζ ln v ζ .
If u and v are functions only of t, then from (32) we obtain the system of ODES:
u t ln v t 2 u 2 t u t = 6 v 1 , u t ln v t 3 u 3 t u t = 6 v γ .
For its solutions
u = w / t
and, if γ 1 , γ 3 / 2 , then
v = 2 γ 3 6 + β w ( 3 2 γ ) / ( γ 1 ) ,
where w = ( γ 1 ) t + α and α , β are constants.
For γ = 1
v = 1 6 + β exp ( t / α ) .
For γ = 3 / 2
v = 1 β 12 ln ( α + t / 2 ) .
Let γ 1 , γ 3 / 2 , and u = γ 1 + α / t . Find solutions of the system (32) of the form v ( t , x ) = v p ( t ) x p with p 0 . For them, U = 0 , V = p , and V x = 0 and equations (32) reduce to one equation:
u t ln v p t = v p x p 6 + 7 p + 2 p 2 + 3 2 γ .
Here, the first and last terms are of order zero on x, and the middle term is of order p 0 . Consequently, this equation has a solution only if the middle term is zero, i.e.,
6 + 7 p + 2 p 2 = 0 .
This equation has two roots: p = 2 , p = 3 / 2 . Moreover, for these values of p it is possible to find a solution of the system (32) of the form
v ( t , x ) = v 0 ( t ) + v p ( t ) x p ,
where v 0 ( t ) is the solution (38) of the system (36). Here, for v p , we obtain the equation
u t ln v p t = 7 p + 3 p 2 v 0 + 2 γ 3 .
Here, the coefficient
7 p + 3 p 2 = 2 , if p = 2 , 15 / 4 , if p = 3 / 2 .
Thus, it is proven.
Theorem 5.
The system (32) reduces to a finite dimensional ODE system in three cases:
1. 
To a one-parameter family of two Equations (35);
2. 
To a system of three Equations (36) and (41) with p = 2 ;
3. 
To the system of three Equations (36) and (41) with p = 3 / 2 .
Up to now, only solutions to the case 2 at u = γ 1 have been known, i.e., solutions to the two-dimensional system of ODEs (see [15,16,17,18]).
Theorem 6.
The system (32) has a one-parameter by α family of solutions:
u = α x p t , p = 12 γ 17 ± 24 γ + 1 12 ( γ 2 ) , v = 25 + 12 γ 7 24 γ + 1 12 ,
where α is an arbitrary constant.
Proof of Theorem 6.
Here, v t = 0 , u + t u t = 0 , U = p , V = 0 , U x = 0 , V x = 0 . Thus, the Equation (32) take the form
v 6 12 p + 6 p 2 = 6 v ( p 1 ) 2 = 1 , v 6 17 p + 12 p 2 = v ( 4 p 3 ) ( 3 p 2 ) = γ .
Substituting the specified values of p and v here, we obtain two identities. □
If in (32) u and v are functions only of x, then they satisfy the system of ODEs:
2 u = v 6 12 U + 7 V + 6 U 2 7 U V + 2 V 2 2 U x x + V x x 1 , 3 u = v 6 17 U + 7 V + 12 U 2 10 U V + 2 V 2 3 U x x + V x x γ .
This is a particular case of the family (35) at σ = 0 .
Below, we assume that each intermediate variable is different from identical zero. Thus, we can consider its logarithm.
After the logarithmic transformation,
τ = ln t , ξ = ln x
the system (32) takes the form
u ( ln v ) τ 2 u 2 u τ = v 6 12 U + 7 V + 6 U 2 7 U V + 2 V 2 2 U ξ + V ξ 1 ,
u ( ln v ) τ 3 u 3 u τ = v 6 17 U + 7 V + 12 U 2 10 U V + 2 V 2 3 U ξ + V ξ γ ,
where U = ln u ξ , V = ln v ξ .
Below, all computations are performed for the system S consisting of a linear combination of the original equations:
1
Equation E 1 S is the difference of the Equations (44) and (45);
2
Equation E 2 S is the difference of the tripled Equation (44) and the doubled Equation (45).
As a result, the S system takes the form
u + u τ = 5 v U 7 v U 2 + 3 U v ξ + v U ξ + γ 1 ,
u ln v τ = 6 v 2 v U + 7 v ξ 6 v U 2 + U v ξ + v ξ V + v ξ ξ + 2 γ 3 .
To apply the Section 4.5 procedures, Equations (46) and (47) of the S system are rewritten as a sum of differential monomials:
E 1 S u 3 + u τ u 2 γ u 2 5 v u ξ u v u ξ , ξ u 3 u ξ v ξ u + 7 v u ξ 2 + u 2 = 0 ,
E 2 S u 3 v τ v ξ , ξ v u 2 6 v 2 u 2 7 v ξ v u 2 2 γ v u 2 v ξ 2 u 2 + + 2 v 2 u ξ u + u ξ v ξ v u + 6 v 2 u ξ 2 + 3 v u 2 = 0 .
The supports of Equations (48) and (49) are
S ( E 1 S ) = { [ 1 , 0 , 3 , 0 ] , [ 0 , 2 , 2 , 1 ] , [ 0 , 1 , 2 , 1 ] , [ 0 , 0 , 2 , 0 ] , [ 0 , 0 , 3 , 0 ] } , 111
S ( E 2 S ) = { [ 1 , 0 , 3 , 1 ] , [ 0 , 2 , 2 , 2 ] , [ 0 , 1 , 2 , 2 ] , [ 0 , 0 , 2 , 1 ] , [ 0 , 0 , 2 , 2 ] } .
To perform computations with a convex polyhedron of large dimension n, it is convenient to represent the latter as an oriented graph, all vertices of which have a unique number j (identifier) and correspond to a generalized face Γ j ( d ) of appropriate dimension d. The top vertex of the graph contains the polyhedron Γ itself, the next level contains generalized faces Γ k ( n 1 ) of dimension n 1 , below are generalized faces Γ k ( n 2 ) of dimension n 2 , and so on. The lowest vertex of the graph is an empty set. The segments connecting vertices of the graph mean that the lower element (the generalized edge) lies in the upper one (the generalized edge of higher dimension). The alternative sum of the number of vertices of the graph in the lines is equal to zero.
The graph of the polyhedron Γ ( E 1 S ) computed by support (50) is shown in Figure 5. The alternative sum of the numbers of elements in the rows is 1 5 + 10 10 + 5 1 = 0 . The polyhedron Γ ( E 1 S ) is a four-dimensional simplex and has five three-dimensional faces with identifiers 161 , 215 , 233 , 239 , 241 , computed by the program. They correspond to the external normals
N 161 ( 3 ) = [ 1 , 0 , 0 , 0 ] , N 215 ( 3 ) = [ 1 , 0 , 1 , 0 ] , N 233 ( 3 ) = [ 0 , 0 , 1 , 1 ] , N 239 ( 3 ) = [ 0 , 1 , 0 , 1 ] , N 241 ( 3 ) = [ 0 , 1 , 0 , 2 ] .
The graph of the polyhedron Γ ( E 2 S ) computed by support (51) is shown in Figure 6. The polyhedron Γ ( E 2 S ) lies in a three-dimensional plane with the normal
N 80 ( 3 ) ( E 2 S ) = [ 1 , 0 , 1 , 0 ]
and is a three-dimensional simplex, i.e., the Equation (49) is quasi-homogeneous.
Let us construct all truncations corresponding to the cone of problem K [ S ] = { p 1 , p 2 0 } according to change (43). The normals N 161 ( 3 ) , N 233 ( 3 ) , N 239 ( 3 ) , and N 80 ( 3 ) fall into the cone of problem K [ S ] . For each of the mentioned normals, we compute the truncations of the system (48), (49) and reject trivial, i.e., those consisting of a single algebraic monomial.
The truncation of Equation (49) corresponding to the normal N 239 ( 3 ) and the truncation of Equation (48) corresponding to the normal N 80 ( 3 ) consist of one algebraic monomial 6 u 2 v 2 and u 3 , respectively. There remain two nontrivial truncations, which we denote according to the notation system of Section 4.4 by S ( 1 ) and S ( 2 ) .
The truncated system S ( 1 ) depends on the variables ξ , u, v and is the system of ODEs, and cone of problem K [ S ( 1 ) ] = { p 1 0 } . The equations of the system have the form:
E 1 S ( 1 ) u 3 γ u 2 5 v u ξ u v u ξ , ξ u 3 u ξ v ξ u + 7 v u ξ 2 + u 2 = 0 ,
E 2 S ( 1 ) ( 3 2 γ ) v u 2 + 6 v 2 u ξ 2 + u ξ v ξ v u + 2 v 2 u ξ u v ξ 2 u 2 7 v ξ v u 2 6 v 2 u 2 v ξ , ξ v u 2 = 0 .
The truncated system of PDEs S ( 2 ) depends on the variables τ , ξ , u, v, and the cone of problem K [ S ( 2 ) ] = { p 1 , p 2 0 } . The equations of the system have the form:
E 1 S ( 2 ) u τ u 2 + u 3 5 v u ξ u v u ξ , ξ u 3 u ξ v ξ u + 7 v u ξ 2 = 0 , 1111
E 2 S ( 2 ) 6 v 2 u ξ 2 + u ξ v ξ v u + 2 v 2 u ξ u v ξ 2 u 2 7 v ξ v u 2 6 v 2 u 2 v ξ , ξ v u 2 + u 3 v τ = 0 .

6. Asymptotic Forms of Solutions to the System S ( 1 )

Consider the computation of asymptotic forms of solutions to the system of ODEs S ( 1 ) in which Equations (52) and (53) depend on variables ξ , u, v, i.e., all corresponding objects of the power geometry are three-dimensional, and the cone of problem K [ S ( 1 ) ] = { p 1 0 } .
The supports of the Equations (52) and (53) are
S ( E 1 S ( 1 ) ) = [ 2 , 2 , 1 ] , [ 1 , 2 , 1 ] , [ 0 , 2 , 0 ] , [ 0 , 3 , 0 ] ,
S ( E 2 S ( 1 ) ) = [ 2 , 2 , 2 ] , [ 1 , 2 , 2 ] , [ 0 , 2 , 1 ] , [ 0 , 2 , 2 ] .
The convex polyhedron Γ ( E 1 S ( 1 ) ) is a tetrahedron, i.e., a three-dimensional simplex with normals to two-dimensional faces, computed by the program,
N 53 ( 2 ) = [ 0 , 1 , 1 ] , N 71 ( 2 ) = [ 0 , 1 , 0 ] , N 77 ( 2 ) = [ 1 , 0 , 1 ] , N 79 ( 2 ) = [ 1 , 0 , 2 ] .
The convex polyhedron Γ ( E 2 S ( 1 ) ) is a two-dimensional simplex, i.e., the left-hand side of the Equation (53) is a quasi-homogeneous differential sum. The corresponding normals are
N 71 ( 2 ) = [ 0 , 1 , 0 ] , N 72 ( 2 ) = [ 0 , 1 , 0 ] .
Suitable normals are those with numbers 53, 71, 77, 72. The corresponding truncated systems are S ( 1 , 1 ) , S ( 1 , 2 ) , S ( 1 , 3 ) , and S ( 1 , 4 ) .
The shortened system S ( 1 , 3 ) contains the trivial shortened equation E 2 S ( 1 , 3 ) 6 u 2 v 2 = 0 , and the shortened system S ( 1 , 4 ) contains the trivial equation E 1 S ( 1 , 4 ) u 3 = 0 . Therefore, we do not consider these systems below.

6.1. Analysis of the Truncated System S ( 1 , 1 )

Making truncation for the normal vector N 53 ( 2 ) = [ 0 , 1 , 1 ] , we obtain a system S ( 1 , 1 ) with equations
E 1 S ( 1 , 1 ) 7 v u ξ 2 3 u ξ v ξ u v u ξ , ξ u 5 v u ξ u + u 3 = 0 ,
E 2 S ( 1 , 1 ) 6 v 2 u ξ 2 + u ξ v ξ v u + 2 v 2 u ξ u v ξ 2 u 2 7 v ξ v u 2 6 v 2 u 2 v ξ , ξ v u 2 = 0 .
The normal vector N 53 ( 2 ) = [ 0 , 1 , 1 ] refers to the case 1 of Section 4.2 and by Lemma 1 defines a power-logarithmic substitution
u = r v , s = ln v ,
converting after reducing Equation (58) by v 3 and Equation (59) by v 4 of system S ( 1 , 1 ) into system P ( 1 , 1 ) with respect to variables ξ , r, s with equations
E 1 P ( 1 , 1 ) 3 r 2 s ξ 2 + r 3 5 r 2 s ξ r 2 s ξ , ξ + 9 r r ξ s ξ 5 r r ξ r r ξ , ξ + 7 r ξ 2 ,
E 2 P ( 1 , 1 ) 5 r 2 s ξ 2 5 r 2 s ξ r 2 s ξ , ξ + 13 r r ξ s ξ 6 r 2 + 2 r r ξ + 6 r ξ 2
with new cone of problem K [ S ( 1 , 1 ) ] = { p 1 , p 3 0 } . The supports of Equations (61) and (62) are
S ( E 1 P ( 1 , 1 ) ) = [ 2 , 2 , 0 ] , [ 2 , 2 , 1 ] , [ 2 , 2 , 2 ] , [ 1 , 2 , 0 ] , [ 1 , 2 , 1 ] , [ 0 , 3 , 0 ] , 1111
S ( E 2 P ( 1 , 1 ) ) = [ 2 , 2 , 0 ] , [ 2 , 2 , 1 ] , [ 2 , 2 , 2 ] , [ 1 , 2 , 0 ] , [ 1 , 2 , 1 ] , [ 0 , 2 , 0 ] .
They differ only in the last point of the support.
The normals to the two-dimensional faces of the convex polyhedron Γ ( E 1 P ( 1 , 1 ) ) of the support (63) are:
N 53 ( 2 ) = [ 1 , 1 , 0 ] , N 107 ( 2 ) = [ 1 , 2 , 0 ] , N 233 ( 2 ) = [ 0 , 1 , 0 ] , N 235 ( 2 ) = [ 0 , 0 , 1 ] , N 237 ( 2 ) = [ 1 , 0 , 1 ] ,
and the convex hull of the support (64) is a two-dimensional simplex with the normals:
N 232 ( 2 ) = [ 0 , 1 , 0 ] , N 233 ( 2 ) = [ 0 , 1 , 0 ] .
Only the normals N 53 ( 2 ) , N 232 ( 2 ) , N 233 ( 2 ) , and N 237 ( 2 ) are suitable, i.e., only they fall within the cone of problem. We denote the corresponding truncated systems by S ( 1 , 1 , 1 ) , S ( 1 , 1 , 2 ) , S ( 1 , 1 , 3 ) , and S ( 1 , 1 , 4 ) , respectively.
The truncated systems S ( 1 , 1 , 1 ) and S ( 1 , 1 , 2 ) are not considered below since they contain trivial equations in the form of a single monomial.

6.1.1. Asymptotic Forms of Solutions to the System S ( 1 , 1 , 3 )

The truncated ODE system S ( 1 , 1 , 3 ) has the form:
E 1 S ( 1 , 1 , 3 ) 7 r ξ 2 r r ξ , ξ 5 r r ξ + 9 r r ξ s ξ r 2 s ξ , ξ 5 r 2 s ξ + 3 r 2 s ξ 2 = 0 ,
E 2 S ( 1 , 1 , 3 ) 5 r 2 s ξ 2 5 r 2 s ξ r 2 s ξ , ξ + 13 r r ξ s ξ 6 r 2 + 2 r r ξ + 6 r ξ 2 = 0 .
The normal vector N 233 ( 2 ) = [ 0 , 1 , 0 ] belongs to the case 1 of Section 4.2 and by Lemma 1 defines the logarithmic transformation
T = ln r ,
translating, after reducing the Equations (65) and (66) by r 2 of the system S ( 1 , 1 , 3 ) into the system P ( 1 , 1 , 3 ) with respect to the variables ξ , T, s with the equations
E 1 P ( 1 , 1 , 3 ) 6 T ξ 2 9 T ξ 3 s ξ 2 + 5 T ξ + T ξ , ξ + 5 s ξ + s ξ , ξ ,
E 2 P ( 1 , 1 , 3 ) = 6 T ξ 2 13 T ξ 5 s ξ 2 2 T ξ + 5 s ξ + s ξ , ξ + 6 , 1111
and with new cone of problem K [ P ( 1 , 1 , 3 ) ] = { p 1 , p 2 , p 3 0 } . The supports S of the Equations (68) and (69) are
S ( E 1 P ( 1 , 1 ) ) = { [ 2 , 0 , 1 ] , [ 2 , 0 , 2 ] , [ 2 , 1 , 0 ] , [ 2 , 1 , 1 ] , [ 2 , 2 , 0 ] , [ 1 , 0 , 1 ] , [ 1 , 1 , 0 ] } ,
S ( E 2 P ( 1 , 1 ) ) = { [ 2 , 0 , 1 ] , [ 2 , 0 , 2 ] , [ 2 , 1 , 1 ] , [ 2 , 2 , 0 ] , [ 1 , 0 , 1 ] , [ 1 , 1 , 0 ] , [ 0 , 0 , 0 ] } .
Consistently computing the convex polyhedra Γ ( E 1 P ( 1 , 1 , 3 ) ) and Γ ( E 2 P ( 1 , 1 , 3 ) ) by supports (70) and (71), respectively, we find the corresponding external normals to their two-dimensional faces Γ ( E 1 P ( 1 , 1 , 3 ) ) and Γ ( E 2 P ( 1 , 1 , 3 ) ) , correspondingly:
N 161 ( 2 ) = [ 1 , 1 , 1 ] , N 185 ( 2 ) = [ 0 , 0 1 ] , N 209 ( 2 ) = [ 0 , 1 , 0 ] , N 233 ( 2 ) = [ 0 , 1 , 1 ] , N 241 ( 2 ) = [ 1 , 0 , 0 ] , N 161 ( 2 ) = [ 1 , 1 , 1 ] , N 209 ( 2 ) = [ 0 , 1 , 0 ] , N 77 ( 2 ) = [ 1 , 1 , 2 ] , N 241 ( 2 ) = [ 1 , 0 , 0 ] .
Only normal N 161 ( 2 ) is suitable, and its corresponding truncated system of ODEs has the form
E 1 S ( 1 , 1 , 3 ) T ξ + s ξ 6 T ξ + 3 s ξ 5 = 0 , E 2 S ( 1 , 1 , 3 ) 6 + 5 s ξ 2 T ξ 5 s ξ 2 13 T ξ s ξ 6 T ξ 2 = 0 .
This system is algebraic with respect to the quantities s ξ and T ξ , and its solutions are the following subsystems:
s ξ = 3 2 , T ξ = 3 2 ,
s ξ = 2 , T ξ = 2 ,
s ξ = 1 , T ξ = 4 3 .
Using the substitutions (60) and (67), we obtain that System (72) defines the asymptotic form
Asymp 1 S ( 1 , 1 , 3 ) : u = C 1 , v = C 2 e 3 ξ / 2 ,
System (73) defines the asymptotic form
Asymp 2 S ( 1 , 1 , 3 ) : u = C 1 , v = C 2 e 2 ξ ,
and System (74) defines asymptotic form
Asymp 3 S ( 1 , 1 , 3 ) : u = C 1 e ξ / 3 , v = C 2 e ξ ,
where C 1 and C 2 are arbitrary constants.

6.1.2. Asymptotic Forms of Solutions to the System S ( 1 , 1 , 4 )

According to Equations (61)–(64), the truncated ODE system, corresponding to N 237 ( 2 ) = [ 1 , 0 , 1 ] , is the following:
r 2 ( ξ ) 3 s ξ 2 5 s ξ + r ξ = 0 , r 2 ( ξ ) 5 s ξ 2 5 s ξ 6 = 0 .
The truncated ODE system S ( 1 , 1 , 4 ) after reduction by r 2 ( ξ ) has the form:
E 1 S ( 1 , 1 , 4 ) 3 s ξ 2 5 s ξ + r ξ = 0 ,
E 2 S ( 1 , 1 , 4 ) 5 s ξ 2 5 s ξ 6 = 0 .
This system is algebraic with respect to the quantities r , s ξ , and its solutions are the following subsystems:
r = a 1 , s ξ = b 1 ,
r = a 2 , s ξ = b 2 ,
where
a 1 , 2 = 13 ± 145 5 , b 1 , 2 = 5 145 10 .
Using the substitution (60), we obtain that the system (81) defines the asymptotic form
Asymp 1 S ( 1 , 1 , 4 ) : u = C 1 a 1 e b 1 ξ , v = C 1 e b 1 ξ ,
and the system (82) defines the asymptotic form
Asymp 2 S ( 1 , 1 , 4 ) : u = C 1 a 2 e b 2 ξ , v = C 1 e b 2 ξ .

6.2. Analysis of the Truncated System S ( 1 , 2 )

Now consider the truncated system S ( 1 , 2 ) for the normal N 71 ( 2 ) = [ 0 , 1 , 0 ] from the system S ( 1 ) with equations:
E 1 S ( 1 , 2 ) γ u 2 5 v u ξ u v u ξ , ξ u 3 u ξ v ξ u + 7 v u ξ 2 + u 2 , 11177
E 2 S ( 1 , 2 ) 3 v u 2 + 6 v 2 u ξ 2 + u ξ v ξ v u + 2 v 2 u ξ u v ξ 2 u 2 v ξ , ξ v u 2 2 γ v u 2 7 v ξ v u 2 6 v 2 u 2 .
The normal vector N 71 ( 2 ) belongs to the case 1 of Section 4.2 and by Lemma 1 defines the logarithmic transformation
r = ln u ,
which, after reducing the Equations (86a) and (86b) of the system S ( 1 , 2 ) by the factor u 2 to the system P ( 1 , 2 ) with the equations
E 1 P ( 1 , 2 ) 6 r ξ 2 v 5 r ξ v r ξ , ξ v 3 r ξ v ξ γ + 1 = 0 , E 2 P ( 1 , 2 ) 6 v 2 r ξ 2 + 2 v 2 r ξ + v r ξ v ξ 6 v 2 7 v ξ v 2 v γ v v ξ , ξ v ξ 2 + 3 v = 0 .
We calculate the supports of the equations of the system P ( 1 , 2 )
S ( E 1 P ( 1 , 2 ) ) = [ 2 , 2 , 1 ] , [ 1 , 1 , 1 ] , [ 2 , 1 , 1 ] , [ 0 , 0 , 0 ] , S ( E 2 P ( 1 , 2 ) ) = [ 2 , 2 , 2 ] , [ 1 , 1 , 2 ] , [ 2 , 1 , 2 ] , [ 0 , 0 , 2 ] , [ 1 , 0 , 2 ] , [ 0 , 0 , 1 ] , [ 2 , 0 , 2 ]
their polyhedra Γ ( E 1 P ( 1 , 2 ) ) , Γ ( E 2 P ( 1 , 2 ) ) and the normals to the two-dimensional faces:
for Γ ( E 1 P ( 1 , 2 ) ) : N 53 ( 2 ) = [ 1 , 1 , 0 ] , N 71 ( 2 ) = [ 0 , 1 , 1 ] , N 77 ( 2 ) = [ 0 , 0 , 1 ] , N 79 = [ 1 , 0 , 2 ] , for Γ ( E 2 P ( 1 , 2 ) ) : N 53 ( 2 ) = [ 1 , 1 , 0 ] , N 72 ( 2 ) = [ 0 , 1 , 0 ] , N 77 ( 2 ) = [ 0 , 0 , 1 ] , N 79 = [ 1 , 0 , 2 ] .
In cone of problem K [ P ( 1 , 2 ) ] = p 1 , p 2 0 only two normals, N 53 ( 2 ) = [ 1 , 1 , 0 ] and N 77 ( 2 ) = [ 0 , 0 , 1 ] , fall in.

6.2.1. Asymptotic Forms of Solutions to the System S ( 1 , 2 , 1 )

The truncation corresponding to the normal N 53 ( 2 ) = [ 1 , 1 , 0 ] gives the system S ( 1 , 2 , 1 ) with the equations:
E 1 S ( 1 , 2 , 1 ) 1 γ 5 r ξ v + 6 r ξ 2 v = 0 , E 2 S ( 1 , 2 , 1 ) 3 v 2 v γ 6 v 2 + 2 v 2 r ξ + 6 v 2 r ξ 2 = 0 ,
which we solve as an algebraic system with respect to the functions r ξ and v:
Sol 1 S ( 1 , 2 , 1 ) : r ξ = a 1 , v = b 1 , Sol 2 S ( 1 , 2 , 1 ) : r ξ = a 2 , v = b 2 ,
where
a 1 , 2 = 12 γ 17 ± 24 γ + 1 12 γ 24 , b 1 , 2 = γ ± 7 24 γ + 1 12 + 25 12 .
Returning to the original variables by (87) and (60), we obtain the asymptotic forms of the solutions
Asymp 1 ( 1 , 2 , 1 ) : u = C 1 e a 1 ξ , v = b 1 , Asymp 2 ( 1 , 2 , 1 ) : u = C 1 e a 2 ξ , v = b 2 .

6.2.2. Asymptotic Forms of Solutions to the System S ( 1 , 2 , 2 )

The truncation corresponding to the normal N 77 ( 2 ) = [ 0 , 0 , 1 ] gives the system S ( 1 , 2 , 2 ) with equations:
E 1 S ( 1 , 2 , 2 ) 5 r ξ v 3 v ξ r ξ r ξ , ξ v + 6 r ξ 2 v = 0 ,
E 2 S ( 1 , 2 , 2 ) 6 v 2 7 v ξ v v ξ 2 v ξ , ξ v + 2 r ξ v 2 + + v ξ r ξ v + 6 r ξ 2 v 2 = 0 .
The normal vector N 77 ( 2 ) belongs to the case 1 of Section 4.2; hence, by Lemma 1 we have a logarithmic transformation
T = ln r ,
which, after reducing Equation (89) by v and Equation (90) by v 2 leads to the system P ( 1 , 2 , 2 ) with cone of problem K [ P ( 1 , 2 , 2 ) ] = { p 1 , p 2 , p 3 0 } :
E 1 P ( 1 , 2 , 2 ) 3 r ξ T ξ 5 r ξ r ξ , ξ + 6 r ξ 2 , E 2 P ( 1 , 2 , 2 ) 2 T ξ 2 7 T ξ 6 T ξ , ξ + r ξ T ξ + 2 r ξ + 6 r ξ 2
The supports of these equations of the system P ( 1 , 2 , 2 ) are
S ( E 1 P ( 1 , 2 , 2 ) ) = { [ 2 , 1 , 1 ] , [ 1 , 1 , 0 ] , [ 2 , 1 , 0 ] , [ 2 , 2 , 0 ] } , S ( E 2 P ( 1 , 2 , 2 ) ) = { [ 2 , 0 , 2 ] , [ 1 , 0 , 1 ] , [ 0 , 0 , 0 ] , [ 2 , 0 , 1 ] , [ 2 , 1 , 1 ] , [ 1 , 1 , 0 ] , [ 2 , 2 , 0 ] } .
Both supports have the following normals:
N 53 ( 2 ) = [ 1 , 1 , 1 ] , N 71 ( 2 ) = [ 0 , 0 , 1 ] , N 77 ( 2 ) = [ 0 , 1 , 0 ] , N 79 ( 2 ) = [ 1 , 0 , 0 ] ,
of which the only normal N 53 ( 2 ) is suitable. The corresponding truncated ODE system has the form
6 r ξ 2 5 r ξ 3 r ξ T ξ = 0 , 6 r ξ 2 6 7 T ξ 2 T ξ 2 + 2 r ξ + r ξ T ξ = 0 .
We obtain an algebraic system with respect to the functions r ξ , T ξ , which has the following solutions:
Sol 1 S ( 1 , 2 , 2 ) : r ξ = 0 , T ξ = 2 , Sol 2 S ( 1 , 2 , 2 ) : r ξ = 0 , T ξ = 3 / 2 , Sol 3 S ( 1 , 2 , 2 ) : r ξ = 1 / 3 , T ξ = 1 .
According to (91), (87), and (60), these solutions correspond to asymptotic forms:
Asymp 1 S ( 1 , 2 , 2 ) : { u = C 1 , v = C 2 e 2 ξ } , Asymp 2 S ( 1 , 2 , 2 ) : { u = C 1 , v = C 2 e 3 ξ / 2 } , Asymp 3 S ( 1 , 2 , 2 ) : { u = C 1 e ξ / 3 , v = C 2 e ξ } .
It is not difficult to see that they correspond to the previously found asymptotic forms in Section 6.1.1.

7. Asymptotic Forms of Solutions to the System S ( 2 )

Now consider the computation of the asymptotic forms of the solutions to the PDE system S ( 2 ) , in which Equations (54) and (55) depend on variables τ , ξ , u, v, and cone of problem K [ S ( 2 ) ] = { p 1 , p 2 0 } .
The normal vector N 233 ( 3 ) ( E 1 S ) = [ 0 , 0 , 1 , 1 ] refers to the case 1 of Section 4.2 and by Lemma 1 defines the power-logarithmic transformation
u = r v , s = ln v ,
reducing the system S ( 2 ) to the system P ( 2 ) with respect to the variables τ , ξ , r, and s with equations:
E 1 P ( 2 ) r 3 s τ + 3 r 2 s ξ 2 + r 3 + r 2 r τ 5 r 2 s ξ r 2 s ξ , ξ + + 9 r r ξ 5 r r ξ r r ξ , ξ + 7 r ξ 2 = 0 ,
E 2 P ( 2 ) r 3 s τ + 5 r 2 s ξ 2 5 r 2 s ξ r 2 s ξ , ξ + 13 r r ξ s ξ 6 r 2 + + 2 r r ξ + 6 r ξ 2 = 0 .
The cone of problem of the system P ( 2 ) is K = { p 1 , p 2 , p 4 0 } .
The supports of Equations (93) and (94) of the system P ( 2 ) are:
S ( E 1 P ( 2 ) ) = { [ 1 , 0 , 3 , 0 ] , [ 1 , 0 , 3 , 1 ] , [ 0 , 2 , 2 , 0 ] , [ 0 , 2 , 2 , 1 ] , [ 0 , 2 , 2 , 2 ] , [ 0 , 1 , 2 , 1 ] , [ 0 , 0 , 3 , 0 ] } , S ( E 2 P ( 2 ) ) = { [ 1 , 0 , 3 , 1 ] , [ 0 , 2 , 2 , 0 ] , [ 0 , 2 , 2 , 1 ] , [ 0 , 2 , 2 , 2 ] , [ 0 , 1 , 2 , 0 ] , [ 0 , 1 , 2 , 1 ] , [ 0 , 0 , 2 , 0 ] }
The normals to the three-dimensional faces of the convex polyhedron Γ ( E 1 P ( 2 ) ) are
N 485 ( 3 ) = [ 0 , 1 , 2 , 0 ] , N 647 ( 3 ) = [ 0 , 1 , 1 , 0 ] , N 701 ( 3 ) = [ 1 , 0 , 1 , 0 ] , N 707 ( 3 ) = [ 0 , 0 , 0 , 1 ] , N 713 ( 3 ) = [ 1 , 1 , 0 , 1 ] , N 727 ( 3 ) = [ 1 , 0 , 0 , 0 ] .
The convex polyhedron Γ ( E 2 P ( 2 ) ) is a three-dimensional simplex, i.e., the support of the equation E 2 P ( 2 ) lies in the hyperplane with normals N 700 ( 3 ) = [ 1 , 0 , 1 , 0 ] and N 701 ( 3 ) .
The normals with numbers 647, 700, 713, and 727 are suitable, and we denote the corresponding systems by S ( 2 , 1 ) , S ( 2 , 2 ) , S ( 2 , 3 ) , and S ( 2 , 4 ) .
The shortened system S ( 2 , 1 ) contains the trivial shortened equation E 2 S ( 2 , 1 ) 6 r 2 = 0 , and the shortened system S ( 2 , 2 ) contains the trivial equation E 1 S ( 2 , 2 ) r 3 = 0 . Therefore, we do not consider these systems below.

7.1. Analysis of the Truncated System S ( 2 , 3 )

The PDE system S ( 2 , 3 ) corresponding to the normal N 713 ( 3 ) = [ 1 , 1 , 0 , 1 ] consists of equations:
E 1 S ( 2 , 3 ) s τ r + 3 s ξ 2 5 s ξ + r = 0 ,
E 2 S ( 2 , 3 ) s τ r + 5 s ξ 2 5 s ξ 6 = 0 ,
derived from the corresponding equations of the system P ( 2 ) after reduction by the multiplier r 2 . Excluding the function r from E 2 S ( 2 , 3 ) and substituting it into E 1 S ( 2 , 3 ) , we obtain the equation:
E 1 S ( 2 , 3 ) 2 s ξ 2 s τ 5 s ξ 2 + 5 s ξ + 6 s τ + 6 = 0 ,
which we consider as one PDE. It can be solved by the method of separation of variables, considering the required function s ( τ , ξ ) in the form of
s ( τ , ξ ) = s 1 ( τ ) + s 2 ( ξ ) .
Then, after substitution, it turns out that Equation (97) can be considered as the equation of an algebraic curve of genus 0 with respect to the derivatives ( s 1 ) τ and ( s 2 ) ξ . This curve allows a rational parametrization
( s 1 ) τ = 5 C 1 2 5 C 1 6 2 ( C 1 2 3 ) , ( s 2 ) ξ = C 1 ,
where C 1 is an arbitrary constant. Hence, the solution of the system S ( 2 , 3 ) is the following:
Sol S ( 2 , 3 ) : r τ , ξ = 2 ( C 1 2 3 ) , s τ , ξ = 5 C 1 2 5 C 1 6 τ 2 ( C 1 2 3 ) + C 1 + C 2 ξ
which, according to (92), in the u , v variables is written as
u = 2 C 2 C 1 2 3 e w , v = C 2 e w ,
where w = 5 C 1 2 5 C 1 6 2 ( C 1 2 3 ) τ + C 1 ξ , and C 2 is an arbitrary constant.

7.2. Analysis of the Truncated System S ( 2 , 4 )

The truncated ODE system S ( 2 , 4 ) is
E 1 S ( 2 , 4 ) 7 r ξ 2 r r ξ , ξ 5 r r ξ + 9 r r ξ s ξ r 2 s ξ , ξ 5 r 2 s ξ + r 3 + 3 r 2 s ξ 2 = 0 ,
E 2 S ( 2 , 4 ) 6 r ξ 2 + 2 r r ξ 6 r 2 + 13 r r ξ s ξ r 2 s ξ , ξ 5 r 2 s ξ + 5 r 2 s ξ 2 = 0 .
Note that Equation (100) differs from Equation (65) of system S ( 1 , 1 , 3 ) only by monomial r 3 , and Equation (101) is exactly the same as Equation (66). Moreover, the variable derivatives τ of the functions r ( τ , ξ ) and s ( τ , ξ ) are not included in the system S ( 2 , 4 ) , which allows us to consider the latter as a ODE system of functions r ( ξ ) and s ( ξ ) that depend on one variable, ξ . Consequently, the objects of power geometry related to the system S ( 2 , 4 ) become three-dimensional in this case. The cone of problem corresponding to the system S ( 2 , 4 ) is K [ S ( 2 , 4 ) ] = { p 1 , p 3 0 } .
The supports of Equations (100) and (101) are
S ( E 1 S ( 2 , 4 ) ) = { [ 2 , 2 , 0 ] , [ 2 , 2 , 1 ] , [ 2 , 2 , 2 ] , [ 1 , 2 , 0 ] , [ 1 , 2 , 1 ] , [ 0 , 3 , 0 ] } , S ( E 2 S ( 2 , 4 ) ) = { [ 2 , 2 , 0 ] , [ 2 , 2 , 1 ] , [ 2 , 2 , 2 ] , [ 1 , 2 , 0 ] , [ 1 , 2 , 1 ] , [ 0 , 2 , 0 ] } ,
and the corresponding vectors of external normals are
N 53 ( 2 ) = [ 1 , 1 , 0 ] , N 107 ( 2 ) = [ 1 , 2 , 0 ] , N 233 ( 2 ) = [ 0 , 1 , 0 ] , N 234 ( 2 ) = [ 0 , 1 , 0 ] , N 235 ( 2 ) = [ 0 , 0 , 1 ] , N 237 ( 2 ) = [ 1 , 0 , 1 ] .
Only the normals with numbers 53, 233, 234, and 237 are suitable.
The truncations corresponding to the first and the third normals are trivial systems.
The truncated system corresponding to the normal N 233 ( 2 ) differs only by the sign from the system P ( 1 , 1 , 3 ) with Equations (68) and (69) from Section 6.1.1. Hence, it defines the same asymptotic forms of the solutions given by the Formulas (75)–(77).
A similar match takes place for the truncated system corresponding to the normal N 237 ( 2 ) , only in this case, the Equations (79) and (80) of the system S ( 1 , 1 , 4 ) from Section 6.1.2 are obtained. Hence, it defines the same asymptotic forms of solutions given by the Formulas (84) and (85).

8. Summary of Results for the System (30)

In this section, we present the final results in the form of exact solutions and asymptotic forms of the solutions to the original system (30) in the initial functions k ( t , x ) and ε ( t , x ) .

8.1. Self-Similar Solutions

The exact solution (34) in variables u , v corresponds to the solution
k = x 2 3 + 2 γ 6 t 2 γ 1 2 , ε = x 2 3 + 2 γ 6 t 3 γ 1 3
The solutions to the system (36) take the following form:
  • For γ 1 , 3 / 2 :
    k = x 2 3 + 2 γ 6 + β t ( γ 1 ) + α δ t ( γ 1 ) + α 2 , ε = x 2 3 + 2 γ 6 + β t ( γ 1 ) + α δ t ( γ 1 ) + α 3 ,
    where δ = ( 3 + 2 γ ) / ( γ 1 ) .
  • For γ = 1 :
    k = x 2 6 + β e t / α α 2 , ε = x 2 6 + β e t / α α 3 .
  • For γ = 3 / 2 :
    k = 4 x 2 β + 12 ln t / 2 + α t + 2 α 2 , ε = 8 x 2 β + 12 ln t / 2 + α t + 2 α 3 .

8.2. Asymptotic Forms of Solutions to the System S ( 1 )

In Section 6, four groups of asymptotics were found, two of which coincided with each other.
The asymptotic forms of the system S ( 1 , 1 , 3 ) :
Asymp 1 S ( 1 , 1 , 3 ) : k = x C 2 t 2 C 1 2 , ε = x C 2 t 3 C 1 3 , Asymp 2 S ( 1 , 1 , 3 ) : k = C 2 t 2 C 1 2 , ε = C 2 t 3 C 1 3 , Asymp 3 S ( 1 , 1 , 3 ) : k = x 1 / 3 C 2 t 2 C 1 2 , ε = C 2 t 3 C 1 3 .
Asymptotic forms of the system S ( 1 , 1 , 4 ) :
Asymp 1 , 2 S ( 1 , 1 , 4 ) : k = x 2 C 1 a 1 , 2 2 x b 1 , 2 t 2 , ε = x 2 C 1 2 a 1 , 2 3 x 2 b 1 , 2 t 3 ,
where a 1 , 2 and b 1 , 2 are given by the Formula (83).
Asymptotic forms of the system S ( 1 , 2 , 1 )
Asymp 1 , 2 S ( 1 , 2 , 1 ) : k = x 2 b 1 , 2 t 2 C 1 2 x 2 a 1 , 2 , ε = x 2 b 1 , 2 t 3 C 1 3 a 1 , 2
where a 1 , 2 and b 1 , 2 are given by the Formula (88).
The asymptotic forms of the system S ( 1 , 2 , 2 ) coincide with the asymptotic forms of the system S ( 1 , 1 , 3 ) .

8.3. Asymptotic Forms of Solutions to the System S ( 2 )

The solution found for the truncated system S ( 2 , 3 ) gives the two-parameter asymptotic form
Asymp S ( 2 , 3 ) : { k = x ( 2 C 1 ) t C 1 2 C 1 3 / 2 C 1 2 6 4 C 2 C 1 2 3 2 , ε = x 2 ( 1 C 1 ) t ( 2 C 1 3 ) ( C 1 1 ) / C 1 2 3 8 C 2 2 C 1 2 3 3 } ,
defined for all parameter values C 1 ± 3 , C 2 0 .
The truncated system S ( 2 , 4 ) does not define new asymptotic forms.

Author Contributions

Conceptualization, A.D.B. and A.B.B.; methodology, A.D.B.; software, A.B.B.; validation, A.D.B. and A.B.B.; writing—original draft preparation, A.B.B.; writing—review and editing, A.D.B.; visualization, A.B.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PDEPartial differential equation
ODEOrdinary differential equation
CASComputer algebra system

References

  1. Bruno, A.D. Power Geometry in Algebraic and Differential Equations; Elsevier Science: Amsterdam, The Netherlands, 2000. [Google Scholar]
  2. Tollmien, W.; Schlichting, H.; Görtler, H.; Riegels, F.W. Über Flüssigkeitsbewegung bei sehr kleiner Reibung. In Ludwig Prandtl Gesammelte Abhandlungen: Zur angewandten Mechanik, Hydro- und Aerodynamik; Riegels, F.W., Ed.; Springer: Berlin/Heidelberg, Germany, 1961; pp. 575–584. [Google Scholar] [CrossRef]
  3. Blasius, H. Grenzschichten in Flüssigkeiten mit kleiner Rebung. Zeit. Math. Phys. 1908, 56, 1–37. [Google Scholar]
  4. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Course of Theoretical Physics; Pergamon Press: Oxford, UK, 1987; Volume 6. [Google Scholar]
  5. Boyd, J.P. The Blasius function: Computations before computers, the value of tricks, undergraduate projects, and open research problems. SIAM Rev. 2008, 50, 791–804. [Google Scholar] [CrossRef]
  6. Bruno, A.D.; Shadrina, T.V. Axisymmetric boundary layer on a needle. Trans. Mosc. Math. Soc. 2007, 68, 201–259. [Google Scholar] [CrossRef] [Green Version]
  7. Dalibard, A.L.; Saint-Raymond, L. Mathematical study of degenerate boundary layers: A Large Scale Ocean Circulation Problem. arXiv 2018, arXiv:1203.5663. [Google Scholar] [CrossRef] [Green Version]
  8. Kuksin, S. Kolmogorov’s theory of turbulence and its rigorous 1d model. arXiv 2021, arXiv:2101.03761. [Google Scholar] [CrossRef]
  9. Bruno, A.D. Asymptotics and expansions of solutions to an ordinary differential equation. Russ. Mathem. Surv. 2004, 59, 429–480. [Google Scholar] [CrossRef]
  10. Euler, L. De relatione inter ternas pluresve quantitates instituenda. Opera Omnia 1785, 4, 136–145. [Google Scholar]
  11. Bruno, A.D. Algorithms of the nonlinear analysis. Russ. Math. Surv. 1996, 51, 956. [Google Scholar]
  12. Thompson, I. Understanding Maple; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  13. Conradi, C.; Feliu, E.; Mincheva, M.; Wiuf, C. Identifying parameter regions for multistationarity. PLoS Comput. Biol. 2017, 13, e1005751. [Google Scholar] [CrossRef] [Green Version]
  14. Kolmogorov, A.N. Local structure of turbulence in an incompressible viscous fluid at very large Reynolds numbers. In Selected Works of A.N. Kolmogorov; Tikhomirov, V.M., Ed.; Kluwer Academic Publisher: Dordrecht, The Netherlands, 1991; Vol. I. Mathematics and Mechanics, pp. 312–318. [Google Scholar]
  15. Kolmogorov, A.N. On the degeneration of isotropic turbulence in an incompressible viscous fluid. In Selected Works of A.N. Kolmogorov; Tikhomirov, V.M., Ed.; Kluwer Academic Publisher: Dordrecht, The Netherlands, 1991; Vol. I. Mathematics and Mechanics, pp. 319–323. [Google Scholar]
  16. Kolmogorov, A.N. The equation of turbulent motion of an incompressible fluid. In Selected Works of A.N. Kolmogorov; Tikhomirov, V.M., Ed.; Kluwer Academic Publisher: Dordrecht, The Netherlands, 1991; Vol. I. Mathematics and Mechanics, pp. 328–330. [Google Scholar]
  17. Bertsch, M.; Dal Pazo, R.; Kersner, R. The evolution of turbulent bursts: The b-epsilon model. Europ. J. Appl. Math. 1994, 5, 537–557. [Google Scholar] [CrossRef]
  18. Galaktionov, V.A. Invariant solutions of two models of evolution of turbulent bursts. Europ. J. Appl. Math. 1999, 10, 237–249. [Google Scholar] [CrossRef]
Figure 1. Streamline of a needle by the filling flow.
Figure 1. Streamline of a needle by the filling flow.
Universe 09 00035 g001
Figure 2. The polyhedron Γ ( f 1 ) of the first equation of the system (11). The boldface point corresponds to the first equation of the truncated system (13).
Figure 2. The polyhedron Γ ( f 1 ) of the first equation of the system (11). The boldface point corresponds to the first equation of the truncated system (13).
Universe 09 00035 g002
Figure 3. The polyhedron Γ ( f 2 ) of the second equation of the system (11). The selected face corresponds to the second equation of the truncated system (13).
Figure 3. The polyhedron Γ ( f 2 ) of the second equation of the system (11). The selected face corresponds to the second equation of the truncated system (13).
Universe 09 00035 g003
Figure 4. The polyhedron Γ ( f 3 ) of the third equation of the system (11). The selected face corresponds to the third equation of the truncated system (13).
Figure 4. The polyhedron Γ ( f 3 ) of the third equation of the system (11). The selected face corresponds to the third equation of the truncated system (13).
Universe 09 00035 g004
Figure 5. Graph of the polyhedron Γ ( E 1 S ) of Equation (48).
Figure 5. Graph of the polyhedron Γ ( E 1 S ) of Equation (48).
Universe 09 00035 g005
Figure 6. Graph of the polyhedron Γ ( E 2 S ) of Equation (49).
Figure 6. Graph of the polyhedron Γ ( E 2 S ) of Equation (49).
Universe 09 00035 g006
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.

Share and Cite

MDPI and ACS Style

Bruno, A.D.; Batkhin, A.B. Asymptotic Forms of Solutions to System of Nonlinear Partial Differential Equations. Universe 2023, 9, 35. https://doi.org/10.3390/universe9010035

AMA Style

Bruno AD, Batkhin AB. Asymptotic Forms of Solutions to System of Nonlinear Partial Differential Equations. Universe. 2023; 9(1):35. https://doi.org/10.3390/universe9010035

Chicago/Turabian Style

Bruno, Alexander D., and Alexander B. Batkhin. 2023. "Asymptotic Forms of Solutions to System of Nonlinear Partial Differential Equations" Universe 9, no. 1: 35. https://doi.org/10.3390/universe9010035

APA Style

Bruno, A. D., & Batkhin, A. B. (2023). Asymptotic Forms of Solutions to System of Nonlinear Partial Differential Equations. Universe, 9(1), 35. https://doi.org/10.3390/universe9010035

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

Article Metrics

Back to TopTop