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:
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:
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
be independent and
be dependent variables. Place
.
Differential monomial is a product of an ordinary monomial
, where
, and a finite number of derivatives of the form
The differential monomial
corresponds to its vector exponent of degree
, formed by the following rules:
where
is the unit vector. The product of monomials corresponds to the sum of their vector exponents of degree:
Differential sum is the sum of differential monomials:
The set
of vector exponents
is called
support of sum . The closure of the convex hull
of the support
is called the
polyhedron of the sum . The boundary
of the polyhedron
consists of generalized faces
, where
,
. Each face
corresponds to:
Normal cone
where the space
is conjugate to the space
,
is a scalar product;
Consider a system of equations:
where
are differential sums. Each equation
corresponds to:
Its support ;
Its polyhedron with a set of faces in the main space ;
Set of their normal cones in the dual space ;
Set of truncated equations .
The set of truncated equations
is a
truncated system if the intersection is not empty:
A truncated system is always a quasi-homogeneous system.
In the solution of the system (
4),
where
are series in powers of
and their logarithms, each
corresponds to its support, polyhedron, normal cones
, and truncations. Here, the logarithm
has a zero exponent of degree on
.
The set of truncated solutions
,
, corresponds to the intersection of their normal cones:
If it is not empty, it corresponds to
truncated solution:
Theorem 1. If the normal cone intersects the normal cone (6), then the truncation , , of this solution satisfies the truncated system (5). Multiplying the differential sum (
5) with the support
by the monomial
gives the differential sum,
, with the support
. 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
be the support of the differential sum
and
. The set
is called
shifted support of the sum
.
Each equation
in the system (
4) corresponds to a support
and a shifted support
. Let
be the convex hull of their union
and
d is the dimension of
. If
, 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
,
. Stationary axisymmetric viscous fluid flows were studied, which at
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
are equivalent to a system of two PDEs for the flow function
and pressure
p:
where density
and viscosity
, with boundary conditions
The system (
8) has the form (
4) with
and
. Thus, supports of Equation (
8) should be considered in
. It turns out that the polyhedra
and
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
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
are equivalent to a system of three PDEs for the flow function
, density
, and enthalpy
h (analog of temperature):
where the parameters
and
, with boundary conditions
and (
10). Here,
,
,
,
,
, so
,
,
. In the space
, all polyhedra
,
,
of equations (
11) are three-dimensional and can be shifted parallel in one linear three-dimensional subspace. In the coordinates
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
. Thus, the truncated system corresponding to the boundary layer on the needle has the form:
with self-similar variables
and the boundary conditions
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
. Thus,
and the system (
13), for variables (
14) is equivalent to a system of two ODEs:
where
, with boundary conditions
The problem (
16), (
17) and (
18) has an invariant manifold
, on which it reduces to one equation
where
is an arbitrary constant, with boundary conditions
An analysis of the solutions of the last problem by methods of planar power geometry [
9] shows that for
, it has solutions of the form
where
is an arbitrary constant.
Thus, at
in the boundary layer
at
and
, the asymptotic form of the flow is obtained:
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.
5. The – Model of Evolution of Turbulent Bursts
According to [
14,
15,
16,
17], the model is described by the system
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, , and , , , .
The support of the first equation
of the system (
30) consists of points
The support of the second equation
of the system (
30) consists of points
The shifted supports
and
consist of three points:
Therefore, .
According to Theorem 3 let us introduce new dependent variables:
This is a power transformation (
22) with matrix (
21), where
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
where
,
.
Let us find the self-similar solutions of this system. Consider two cases.
The first case:
are constants. Then, the system (
32) has the form
Its solution
has two critical values:
and
.
The second case: Let
, where
. Now
u and
v are functions of
. In this case, in the matrix (
21), the submatrix
and the submatrices
and
are the same as before. For
and
, the system (
32) after substitutions
generates a one-parameter by
family of systems of two ODEs:
where
,
.
If
u and
v are functions only of
t, then from (
32) we obtain the system of ODES:
For its solutions
and, if
,
, then
where
and
,
are constants.
Let
,
, and
. Find solutions of the system (
32) of the form
with
. For them,
,
, and
and equations (
32) reduce to one equation:
Here, the first and last terms are of order zero on
x, and the middle term is of order
. Consequently, this equation has a solution only if the middle term is zero, i.e.,
This equation has two roots:
,
. Moreover, for these values of
p it is possible to find a solution of the system (
32) of the form
where
is the solution (
38) of the system (
36). Here, for
, we obtain the equation
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 ; - 3.
To the system of three Equations (36) and (41) with .
Up to now, only solutions to the case 2 at
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:where α is an arbitrary constant. Proof of Theorem 6. Here,
,
,
,
,
,
. Thus, the Equation (
32) take the form
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:
This is a particular case of the family (
35) at
.
Below, we assume that each intermediate variable is different from identical zero. Thus, we can consider its logarithm.
After the logarithmic transformation,
the system (
32) takes the form
where
,
.
Below, all computations are performed for the system S consisting of a linear combination of the original equations:
- 1
Equation
is the difference of the Equations (
44) and (
45);
- 2
Equation
is the difference of the tripled Equation (
44) and the doubled Equation (
45).
As a result, the
S system takes the form
To apply the
Section 4.5 procedures, Equations (
46) and (
47) of the
S system are rewritten as a sum of differential monomials:
The supports of Equations (
48) and (
49) are
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 of appropriate dimension d. The top vertex of the graph contains the polyhedron itself, the next level contains generalized faces of dimension , below are generalized faces of dimension , 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
computed by support (
50) is shown in
Figure 5. The alternative sum of the numbers of elements in the rows is
. The polyhedron
is a four-dimensional simplex and has five three-dimensional faces with identifiers
, computed by the program. They correspond to the external normals
The graph of the polyhedron
computed by support (
51) is shown in
Figure 6. The polyhedron
lies in a three-dimensional plane with the normal
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
according to change (
43). The normals
,
,
, and
fall into the cone of problem
. 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
and the truncation of Equation (
48) corresponding to the normal
consist of one algebraic monomial
and
, respectively. There remain two nontrivial truncations, which we denote according to the notation system of
Section 4.4 by
and
.
The truncated system
depends on the variables
,
u,
v and is the system of ODEs, and cone of problem
. The equations of the system have the form:
The truncated system of PDEs
depends on the variables
,
,
u,
v, and the cone of problem
. The equations of the system have the form:
6. Asymptotic Forms of Solutions to the System
Consider the computation of asymptotic forms of solutions to the system of ODEs
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
.
The supports of the Equations (
52) and (
53) are
The convex polyhedron
is a tetrahedron, i.e., a three-dimensional simplex with normals to two-dimensional faces, computed by the program,
The convex polyhedron
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
Suitable normals are those with numbers 53, 71, 77, 72. The corresponding truncated systems are , , , and .
The shortened system contains the trivial shortened equation , and the shortened system contains the trivial equation . Therefore, we do not consider these systems below.
6.1. Analysis of the Truncated System
Making truncation for the normal vector
, we obtain a system
with equations
The normal vector
refers to the case 1 of
Section 4.2 and by Lemma 1 defines a power-logarithmic substitution
converting after reducing Equation (
58) by
and Equation (
59) by
of system
into system
with respect to variables
,
r,
s with equations
with new cone of problem
. The supports of Equations (
61) and (
62) are
They differ only in the last point of the support.
The normals to the two-dimensional faces of the convex polyhedron
of the support (
63) are:
and the convex hull of the support (
64) is a two-dimensional simplex with the normals:
Only the normals , , , and are suitable, i.e., only they fall within the cone of problem. We denote the corresponding truncated systems by , , , and , respectively.
The truncated systems and 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
The truncated ODE system
has the form:
The normal vector
belongs to the case 1 of
Section 4.2 and by Lemma 1 defines the logarithmic transformation
translating, after reducing the Equations (
65) and (
66) by
of the system
into the system
with respect to the variables
,
T,
s with the equations
and with new cone of problem
. The supports
of the Equations (
68) and (
69) are
Consistently computing the convex polyhedra
and
by supports (
70) and (
71), respectively, we find the corresponding external normals to their two-dimensional faces
and
, correspondingly:
Only normal
is suitable, and its corresponding truncated system of ODEs has the form
This system is algebraic with respect to the quantities
and
, and its solutions are the following subsystems:
Using the substitutions (
60) and (
67), we obtain that System (
72) defines the asymptotic form
System (
73) defines the asymptotic form
and System (
74) defines asymptotic form
where
and
are arbitrary constants.
6.1.2. Asymptotic Forms of Solutions to the System
According to Equations (
61)–(
64), the truncated ODE system, corresponding to
, is the following:
The truncated ODE system
after reduction by
has the form:
This system is algebraic with respect to the quantities
, and its solutions are the following subsystems:
where
Using the substitution (
60), we obtain that the system (
81) defines the asymptotic form
and the system (
82) defines the asymptotic form
6.2. Analysis of the Truncated System
Now consider the truncated system
for the normal
from the system
with equations:
The normal vector
belongs to the case 1 of
Section 4.2 and by Lemma 1 defines the logarithmic transformation
which, after reducing the Equations (
86a) and (
86b) of the system
by the factor
to the system
with the equations
We calculate the supports of the equations of the system
their polyhedra
,
and the normals to the two-dimensional faces:
In cone of problem only two normals, and , fall in.
6.2.1. Asymptotic Forms of Solutions to the System
The truncation corresponding to the normal
gives the system
with the equations:
which we solve as an algebraic system with respect to the functions
and
v:
where
Returning to the original variables by (
87) and (
60), we obtain the asymptotic forms of the solutions
6.2.2. Asymptotic Forms of Solutions to the System
The truncation corresponding to the normal
gives the system
with equations:
The normal vector
belongs to the case 1 of
Section 4.2; hence, by Lemma 1 we have a logarithmic transformation
which, after reducing Equation (
89) by
v and Equation (
90) by
leads to the system
with cone of problem
:
The supports of these equations of the system
are
Both supports have the following normals:
of which the only normal
is suitable. The corresponding truncated ODE system has the form
We obtain an algebraic system with respect to the functions
,
, which has the following solutions:
According to (
91), (
87), and (
60), these solutions correspond to asymptotic forms:
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
Now consider the computation of the asymptotic forms of the solutions to the PDE system
, in which Equations (
54) and (
55) depend on variables
,
,
u,
v, and cone of problem
.
The normal vector
refers to the case 1 of
Section 4.2 and by Lemma 1 defines the power-logarithmic transformation
reducing the system
to the system
with respect to the variables
,
,
r, and
s with equations:
The cone of problem of the system is .
The supports of Equations (
93) and (
94) of the system
are:
The normals to the three-dimensional faces of the convex polyhedron
are
The convex polyhedron is a three-dimensional simplex, i.e., the support of the equation lies in the hyperplane with normals and .
The normals with numbers 647, 700, 713, and 727 are suitable, and we denote the corresponding systems by , , , and .
The shortened system contains the trivial shortened equation , and the shortened system contains the trivial equation . Therefore, we do not consider these systems below.
7.1. Analysis of the Truncated System
The PDE system
corresponding to the normal
consists of equations:
derived from the corresponding equations of the system
after reduction by the multiplier
. Excluding the function
r from
and substituting it into
, we obtain the equation:
which we consider as one PDE. It can be solved by the method of separation of variables, considering the required function
in the form of
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
and
. This curve allows a rational parametrization
where
is an arbitrary constant. Hence, the solution of the system
is the following:
which, according to (
92), in the
variables is written as
where
, and
is an arbitrary constant.
7.2. Analysis of the Truncated System
The truncated ODE system
is
Note that Equation (
100) differs from Equation (
65) of system
only by monomial
, and Equation (
101) is exactly the same as Equation (
66). Moreover, the variable derivatives
of the functions
and
are not included in the system
, which allows us to consider the latter as a ODE system of functions
and
that depend on one variable,
. Consequently, the objects of power geometry related to the system
become three-dimensional in this case. The cone of problem corresponding to the system
is
.
The supports of Equations (
100) and (
101) are
and the corresponding vectors of external normals are
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
differs only by the sign from the system
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
, only in this case, the Equations (
79) and (
80) of the system
from
Section 6.1.2 are obtained. Hence, it defines the same asymptotic forms of solutions given by the Formulas (
84) and (
85).