Next Article in Journal
SCC-GPT: Source Code Classification Based on Generative Pre-Trained Transformers
Next Article in Special Issue
Analytical Computation of Hyper-Ellipsoidal Harmonics
Previous Article in Journal
Characterizing Finite Groups through Equitable Graphs: A Graph-Theoretic Approach
Previous Article in Special Issue
Existence and Uniqueness of Weak Solutions to Frictionless-Antiplane Contact Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unsteady Magnetohydrodynamics PDE of Monge–Ampère Type: Symmetries, Closed-Form Solutions, and Reductions

by
Andrei D. Polyanin
1,* and
Alexander V. Aksenov
2,*
1
Ishlinsky Institute for Problems in Mechanics, Russian Academy of Sciences, 101 Vernadsky Avenue, Bldg 1, 119526 Moscow, Russia
2
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, 1 Leninskie Gory, Main Building, 119991 Moscow, Russia
*
Authors to whom correspondence should be addressed.
Mathematics 2024, 12(13), 2127; https://doi.org/10.3390/math12132127
Submission received: 19 May 2024 / Revised: 2 July 2024 / Accepted: 4 July 2024 / Published: 6 July 2024

Abstract

:
The paper studies an unsteady equation with quadratic nonlinearity in second derivatives, that occurs in electron magnetohydrodynamics. In mathematics, such PDEs are referred to as parabolic Monge–Ampère equations. An overview of the Monge–Ampère type equations is given, in which their unusual qualitative features are noted. For the first time, the Lie group analysis of the considered highly nonlinear PDE with three independent variables is carried out. An eleven-parameter transformation is found that preserves the form of the equation. Some one-dimensional reductions allowing to obtain self-similar and other invariant solutions that satisfy ordinary differential equations are described. A large number of new additive, multiplicative, generalized, and functional separable solutions are obtained. Special attention is paid to the construction of exact closed-form solutions, including solutions in elementary functions (in total, more than 30 solutions in elementary functions were obtained). Two-dimensional symmetry and non-symmetry reductions leading to simpler partial differential equations with two independent variables are considered (including stationary Monge–Ampère type equations, linear and nonlinear heat type equations, and nonlinear filtration equations). The obtained results and exact solutions can be used to evaluate the accuracy and analyze the adequacy of numerical methods for solving initial boundary value problems described by highly nonlinear partial differential equations.

Contents

1.   
2.   
3.   
4.   
5.   
6.   
7.   
8.   
9.   
10. 
11. 
12. 
13. 
14. 
15. 
      
References

1. Introduction

1 ° . In plasma electron magnetohydrodynamics, an approach has been developed in which the electron component of the plasma is represented as a set of random or regularly distributed point electron vortices. Local unsteady motions in the corresponding two-dimensional vortex system are described by the nonlinear equation [1,2,3]:
u t = u x x u y y u x y 2 ,
in the right side of which an unessential constant factor is omitted.
In the special case, for u t = 0 , Equation (1) reduces to the stationary homogeneous Monge–Ampère equation
u x x u y y u x y 2 = 0 .
First integrals, exact solutions, symmetries, and invariant transformations of this equation and its multidimensional generalization can be found, for example, in [4,5,6,7,8,9,10]. It is important to note that the general solution of Equation (2) can be presented in parametric form [4,10].
The quadratic combination of second derivatives
u x x u y y u x y 2 det u x x u x y u x y u y y ,
which is included in Equations (1) and (2) and defined as the determinant of the matrix of second derivatives, will be called nonlinearity of the Monge–Ampère type.
For the first time, nonstationary equations with nonlinearity of the Monge–Ampère type for spatial variables appeared in [11], in which, in particular, the equation was given
u t = det [ u x i x j ] f ( x , t ) ,
where x = ( x 1 , , x n ) and u x i x j = 2 u x i x j . The matrix of second derivatives [ u x i x j ] included in Equation (3) describes the local curvature of a function of many variables and is called the Hessian matrix.
In the case of two spatial variables x 1 = x and x 2 = y in Equation (3), we should put x = ( x , y ) and det [ u x i x j ] = u x x u y y u x y 2 . For n = 2 and f 0 , Equation (3) becomes Equation (1).
Equations (1) and (3) and more complex related nonlinear PDEs containing the first time derivative u t and a combination of the second derivatives in spatial coordinates of the form det [ u x i x j ] are called parabolic Monge–Ampère equations. Such PDEs were considered in many works (see, for example, [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]), in which, mainly, questions of the existence and uniqueness of some classes of solutions for various internal and external initial-boundary value problems were studied, and geometric applications were also discussed.
It should be noted that in [11,12,14,16,18,22,23,26] the parabolic Monge–Ampère equation
u t det [ u x i x j ] = f ( x , t ) ,
which has a more complex nonlinearity than PDE (3), was studied.
In the case of two spatial variables, the parabolic Monge–Ampère equations have the form
u t = F ( x , y , t , u , u x , u y , u x x u y y u x y 2 ) .
2 ° . At u t = F t = 0 , the parabolic Equation (4) are reduced to the stationary Monge–Ampère equations, which are usually written in the form resolved with respect to the combination of the highest derivatives
u x x u y y u x y 2 = Φ ( x , y , u , u x , u y ) .
Very complete information about stationary Monge–Ampère equations, which are often encountered in differential geometry, can be found in the book [29]. In the special case, for Φ = a = const < 0 , the general solution of Equation (5) can be represented in parametric form [4,10]. Various transformations, exact solutions, and linearizable equations of the form (5) are given in [7,9,10,30,31,32,33,34].
Equation (5) is highly nonlinear (quadratic with respect to the higher derivatives) and has properties unusual for quasi-linear equations of mathematical physics that are linear with respect to the higher derivatives. In particular, even for the simplest right-hand side Φ = a = const , the qualitative features of this equation depend on the sign of the constant a, since for a > 0 , Equation (5) is an equation of elliptic type, and for a < 0 it is an equation of hyperbolic type [10,29,35].
The above-mentioned qualitative features of the stationary Monge–Ampère equations significantly complicate the procedure of their numerical integration. Therefore, sometimes for the numerical solution of the stationary Monge–Ampère equation (5), various modifications of the time-relaxation method are used using the parabolic Monge–Ampère equation, which is obtained by adding the nonstationary first derivative u t to the right side of (5). The indicated approach involving the parabolic Monge–Ampère equations (1) and (3) was used in [36,37,38,39] for the numerical integration of boundary value problems described by the corresponding stationary Monge–Ampère equations.
3 ° . The symmetries of the parabolic Monge–Ampère equation (1) have not been studied previously, and reductions and exact solutions of this highly nonlinear PDE have been little considered. In the works [40,41] by the classical method of separation of variables, exact solutions of Equation (1) were obtained in the form of a product and a sum of functions of different arguments, by analogy with how this is performed for linear second-order PDEs. A generalized separable solution of polynomial form in the spatial variables x and y with coefficients depending on time t is described in [42] (for more details on this solution, see Item 9 ° in Section 8).
Reductions of the equation under consideration are equations with fewer independent variables or equations of lower order, all solutions of which are solutions of the given equation. Reductions play a key role in constructing exact solutions to differential equations and usually lead to lower-order or lower-dimensional equations. The most important for nonlinear partial differential equations are one-dimensional reductions, using which it is possible to present their solutions in terms of solutions of much simpler ordinary differential equations.
In this paper, exact closed-form solutions of nonlinear partial differential equations are understood as solutions that are expressed in terms of:
(i)
elementary functions,
(ii)
elementary functions and indefinite integrals (solutions in quadrature),
(iii)
solutions of ODE or ODE systems.
In addition to exact closed-form solutions, we will also consider very meaningful solutions that can be expressed in terms of solutions of linear PDEs.
It is important to note that exact solutions of nonlinear partial differential equations of mathematical physics play an important role in mathematical standards, which are widely used to evaluate accuracy, verify and develop various numerical, asymptotic, and approximate analytical methods.
Symmetries, reductions, and exact solutions of nonlinear partial differential equations of mathematical physics are most often constructed using the classical method of symmetry reductions [5,9,43,44,45,46,47,48], the direct method of symmetry reductions [10,49,50,51,52,53,54,55], the nonclassical symmetries methods [50,52,56,57,58,59,60,61,62,63,64,65], methods of generalized separation of variables [10,42,55,66,67,68,69,70,71,72,73], methods of functional separation of variables [10,55,60,70,74,75,76,77,78,79,80,81,82,83], the method of differential constraints [10,51,55,84,85,86,87], and some other exact analytical methods [10,55,88,89,90,91,92,93,94,95]. On methods for constructing exact solutions of nonlinear delay PDEs and functional PDEs, see, for example, [96,97,98,99,100,101,102,103,104].
Let us note a couple of characteristic qualitative features of Equation (1) that distinguish it from the vast majority of other equations of mathematical physics that do not explicitly depend on independent variables. This equation has neither simple solutions of the traveling wave type u = U ( a x + b y + c t ) , nor more complex solutions of the form u = V ( a x + b y , t ) , where a, b, c are any constants, and U ( z ) and V ( ξ , t ) are some functions other than constants. This circumstance complicates to some extent the search for exact solutions of the PDE under consideration.
4 ° . In this paper, to find exact solutions to the nonlinear magnetohydrodynamics equation (1), we used the classical method of symmetry reductions [5,9,43,44], methods of generalized and functional separation of variables [10,42,55], as well as given in [10,34], exact solutions of simpler than the original PDE, intermediate reduced equations with fewer independent variables. In some cases, exact solutions were obtained by combining the classical method and methods of generalized and/or functional separation of variables (the order in which these methods are applied may vary, and the solutions found in this way will be non-invariant). Special attention is paid to the construction of simple exact solutions that are expressed in terms of elementary functions or in closed form (in quadratures). It is most convenient to use such solutions as test problems to evaluate the accuracy and verify the adequacy of numerical methods for solving highly nonlinear partial differential equations.

2. Symmetries of the Magnetohydrodynamics Equation, Admissible Operators, and Invariant Transformations

We are looking for the symmetry operators of Equation (1) in the form
X = ξ 1 ( x , y , t , u ) x + ξ 2 ( x , y , t , u ) y + ξ 3 ( x , y , t , u ) t + η ( x , y , t , u ) u .
Applying the invariance criterion [5], we obtain the following overdetermined linear homogeneous system of determining equations
ξ t 1 = 0 , ξ u 1 = 0 , ξ t 2 = 0 , ξ u 2 = 0 , ξ x 3 = 0 , ξ y 3 = 0 , ξ u 3 = 0 , η t = 0 , η u 2 ξ x 1 2 ξ y 2 + ξ t 3 = 0 , ξ x x 1 = 0 , ξ x y 1 = 0 , ξ y y 1 = 0 , ξ x x 2 = 0 , ξ x y 2 = 0 , ξ y y 2 = 0 , ξ t t 3 = 0 , η x x = 0 , η x y = 0 , η y y = 0 .
It is not difficult to show that the joint solution of the overdetermined system of Equation (6) is given by the formulas
ξ 1 = c 1 x + c 2 y + c 3 , ξ 2 = c 4 x + c 5 y + c 6 , ξ 3 = c 7 t + c 8 , η = ( 2 c 1 + 2 c 5 c 7 ) u + c 9 x + c 10 y + c 11 ,
where c j ( j = 1 , , 11 ) are arbitrary constants. From which two propositions follow.
Proposition 1.
The basis of the Lie algebra of symmetry operators of Equation (1) has the form
X 1 = x , X 2 = y , X 3 = t , X 4 = u , X 5 = y x , X 6 = x y , X 7 = x u , X 8 = y u , X 9 = x x + 2 u u , X 10 = y y + 2 u u , X 11 = t t u u .
Proposition 2.
The transformation
x ¯ = a 1 x + b 1 y + c 1 , y ¯ = a 2 x + b 2 y + c 2 , t ¯ = d 1 t + d 2 , u ¯ = ( a 1 b 2 a 2 b 1 ) 2 d 1 u + a 3 x + b 3 y + c 3 ( a 1 b 2 a 2 b 1 0 , d 1 0 )
transforms any solution u = φ ( x , y , t ) of Equation (1) into an eleven-parameter family of solutions
u = d 1 ( a 1 b 2 a 2 b 1 ) 2 φ ( a 1 x + b 1 y + c 1 , a 2 x + b 2 y + c 2 , d 1 t + d 2 ) a 3 x b 3 y c 3 ,
where a 1 , a 2 , a 3 , b 1 , b 2 , b 3 , c 1 , c 2 , c 3 , d 1 , d 2 are arbitrary constants.
It is important to note that in the Formula (8) the free parameters can be complex numbers if the solution obtained in this way is valid (for details, see [94]). Section 4 provides a specific example of using this approach.
Remark 1.
The transformation (7) leaves the form of Equation (1) invariant. The Formula (8) makes it possible to obtain multi-parameter families of solutions from simpler solutions.

3. Two-Dimensional Similarity Reductions

The regular procedure for constructing two-dimensional similarity reductions of partial differential equations is described in [5,44]. In our work, we will confine ourselves to the most informative examples of the construction of two-dimensional reductions of the parabolic Monge–Ampère equation (1) based on the symmetries described above.
1 ° . Passing in (1) to variables of the traveling wave type
u = U ( ξ , η ) , ξ = x a 1 t , η = y a 2 t ,
where a 1 and a 2 are arbitrary constants, we arrive at a two-dimensional Monge–Ampère type equation with constant coefficients
U ξ ξ U η η U ξ η 2 + a 1 U ξ + a 2 U η = 0 .
The solution of the form (9) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = a 1 X 1 + a 2 X 2 + X 3 = a 1 x + a 2 y + t .
Equation (10) admits exact solutions, quadratic in any independent variable, of the form
U 1 = f 1 ( ξ ) η 2 + g 1 ( ξ ) η + h 1 ( ξ ) , U 2 = f 2 ( η ) ξ 2 + g 2 ( η ) ξ + h 2 ( η ) ,
where the functions f i , g i , h i ( i = 1 , 2 ) are described by the corresponding one-dimensional systems of three ODEs, which are omitted here.
Remark 2.
The consistent use of Formulas (9) and (11) leads to exact non-invariant solutions obtained by combining the classical method and the method of generalized separation of variables.
Remark 3.
Generalized separable solutions that are quadratic in one or more variables, similar to (11), are often used to construct exact solutions to reaction-diffusion equations and some other nonlinear PDEs (see, for example, [10,42,55,71]).
2 ° . Passing in (1) to variables of self-similar type
u = t 2 α 2 β 1 U ( ξ , η ) , ξ = x t α , η = y t β ,
where α and β are arbitrary constants, we obtain a two-dimensional equation of the Monge–Ampère type with variable coefficients for the lower derivatives
U ξ ξ U η η U ξ η 2 α ξ U ξ β η U η + ( 2 α + 2 β + 1 ) U = 0 .
The solution of the form (12) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = α X 9 + β X 10 X 11 = α x x + β y y t t + ( 2 α + 2 β + 1 ) u u .
Equation (13) admits exact non-invariant solutions, quadratic with respect to any independent variable, of the form (11). The values of α = β = 0 in (12) correspond to a multiplicative separable solution.
3 ° . Passing in (1) to variables of limit self-similar type
u = e 2 ( α + β ) t U ( ξ , η ) , ξ = x e α t , η = y e β t ,
where α and β are arbitrary constants, we obtain another two-dimensional equation of the Monge–Ampère type with variable coefficients for the lower derivatives
U ξ ξ U η η U ξ η 2 α ξ U ξ β η U η + 2 ( α + β ) U = 0 .
The solution of the form (14) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = α X 9 + β X 10 X 3 = α x x + β y y t + 2 ( α + β ) u u .
Equation (15) admits exact non-invariant solutions, quadratic with respect to any independent variable, of the form (11).
4 ° . Passing in (1) to invariant variables
u = 1 t U ( ξ , η ) , ξ = x + λ 1 ln t , η = y + λ 2 ln t ,
where λ 1 and λ 2 are arbitrary constants, we obtain a two-dimensional equation of the Monge–Ampère type with constant coefficients for the lower derivatives
U ξ ξ U η η U ξ η 2 λ 1 U ξ λ 2 U η + U = 0 .
The solution of the form (16) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = λ 1 X 1 + λ 2 X 2 X 11 = λ 1 x + λ 2 y t t + u u .
Equation (16) admits traveling wave solutions, as well as non-invariant solutions quadratic in any independent variable of the form (11).
The values of λ 1 = λ 2 = 0 in (16) correspond to a multiplicative separable solution.
5 ° . Passing in (1) to other invariant variables
u = x 2 U ( ξ , η ) , ξ = t + α ln x , η = y + β ln x ,
where λ 1 and λ 2 are arbitrary constants, we obtain a two-dimensional equation of the Monge–Ampère type, which is omitted here. The values of α = β = 0 in (17) correspond to a multiplicative separable solution.
The solution of the form (17) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = β X 2 α X 3 + X 9 = x x β y α t + 2 u u .
6 ° . Equation (1) using invariant variables
u = e ( α 2 β ) x U ( ξ , η ) , ξ = t e α x , η = y e β x ,
is reduced to a two-dimensional Monge–Ampère type equation, which is omitted here.
The solution of the form (18) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = X 1 β X 10 α X 11 = x β y y α t t + ( α 2 β ) u u .
7 ° . Equation (1) using invariant variables
u = e 2 α x U ( ξ , η ) , ξ = x + β t , η = y e α x ,
is reduced to a two-dimensional Monge–Ampère type equation, which is omitted here.
The solution of the form (19) is invariant with respect to the one-parameter group of transformations given by the symmetry operator
Y = β X 1 X 3 α β X 10 + X 9 = β x α β y y t 2 α β u u .
Remark 4.
More complex two-dimensional reductions of Equation (1) can be obtained by applying to the solutions (9), (12), (14), (16), (17)–(19) the reproduction Formula (8).

4. One-Dimensional Similarity Reductions and Exact Solutions

The standard procedure for obtaining one-dimensional similarity reductions of PDEs is outlined in [5]. In our work, we will give only a few illustrative examples of constructing invariant solutions (in terms of ODEs) of the parabolic Monge–Ampère equation (1) based on the symmetries described above. We will also present several exact solutions of this PDE in elementary functions.
1 ° . The simplest invariant solution of Equation (1), which allows a scaling transformation, is a solution in the form of a product of the corresponding degrees of independent variables
u = x 2 y 2 12 t .
Below, we consider several more complex invariant solutions that can be obtained from the solution (20) using simple methods described in [93,94].
The solution (20) is a special case of a wider family of invariant solutions of the form
u = x 2 t f ( z ) , z = y + β ln t ,
where β is an arbitrary constant, and the function f = f ( z ) is described by the ODE
2 f f z z 4 ( f z ) 2 β f z + f = 0 .
The solution (20) is a special case of another broader family of invariant solutions of the form
u = x 2 t g ( ξ ) , ξ = y + λ ln x ,
where λ is an arbitrary constant, and the function f = f ( z ) satisfies the ODE
λ g ξ g ξ ξ 2 g g ξ ξ + 4 ( g ξ ) 2 = g .
The solution (20) is also a special case of another broader family of invariant solutions of the form
u = x 2 y 2 φ ( η ) , η = t + γ ln y ,
where γ is an arbitrary constant, and the function φ = φ ( η ) is described by the ODE
2 γ 2 φ φ η η 4 γ 2 ( φ η ) 2 10 γ φ φ η φ η 12 φ 2 = 0 .
In solutions (21)–(23), the spatial variables x and y can be swapped or use the Formula (8).
Let us apply the Formula (8) with a 1 = a 2 = b 2 = d 1 = 1 , b 1 = 1 , c 1 = c 2 = d 2 = a 3 = b 3 = c 3 = 0 to solution (20). As a result, we obtain a solution of a more complex form,
u = ( x 2 y 2 ) 2 48 t .
Following [93,94], we show how another exact solution can be obtained from solution (24) by using a complex parameter. Equation (1) does not change if we make the transformation u ¯ = u , x ¯ = i x , where i 2 = 1 (this is equivalent to choosing the complex parameter a 1 = i in Formula (8)). Having made the same transformation in (24), we arrive at the following simple solution with axial symmetry:
u = ( x 2 + y 2 ) 2 48 t .
2 ° . Equation (1) using invariant variables
u = x 2 2 β t 2 α 1 V ( ζ ) , ζ = x β t α y
is reduced to a second-order ODE
β ( β + 1 ) ζ V ζ 2 ( β 1 ) ( 2 β 1 ) V V ζ ζ + ( β 2 ) 2 ( V ζ ) 2 + α ζ V ζ ( 2 α + 1 ) V = 0 .
The solution of the form (26) is an invariant solution for the two-parameter transformation group given by symmetry operators
Y 1 = α X 9 β X 11 = α x x β t t + ( 2 α + β ) u u , Y 2 = α X 10 X 11 = α y y t t + ( 2 α + 1 ) u u .
3 ° . Equation (1) using invariant variables
u = e 2 α t x 2 2 β V ( ζ ) , ζ = e α t x β y
is reduced to a second-order ODE
β ( β + 1 ) ζ V ζ 2 ( β 1 ) ( 2 β 1 ) V V ζ ζ + ( β 2 ) 2 ( V ζ ) 2 + α ζ V ζ 2 α V = 0 .
The solution of the form (27) is an invariant solution for the two-parameter transformation group given by symmetry operators
Y 1 = X 3 + α X 10 = α y y t + 2 α u u , Y 2 = X 9 β X 10 = x x β y y 2 ( β 1 ) u u .
4 ° . Equation (1) using invariant variables
u = t 2 α 1 e 2 β x V ( ζ ) , ζ = t α e β x y
is reduced to a second-order ODE
β 2 ( ζ V ζ 4 V ) V ζ ζ + β 2 ( V ζ ) 2 + α ζ V ζ ( 2 α + 1 ) V = 0 .
The solution of the form (28) is an invariant solution for the two-parameter transformation group given by symmetry operators
Y 1 = X 1 β X 10 = x β y y 2 β u u , Y 2 = α X 10 X 11 = α y y t t + ( 2 α + 1 ) u u .
5 ° . Equation (1) using invariant variables
u = e 2 α t 2 β x V ( ζ ) , ζ = e α t + β x y
is reduced to a second-order ODE
β 2 ( ζ V ζ 4 V ) V ζ ζ + β 2 ( V ζ ) 2 + α ζ V ζ 2 α V = 0 .
The solution of the form (29) is an invariant solution for the two-parameter transformation group given by symmetry operators
Y 1 = X 1 β X 10 = x β y y 2 β u u , Y 2 = X 3 α X 10 = t α y y 2 α u u .
Remark 5.
More complex two-dimensional reductions of Equation (1) can be obtained by applying to the solutions (22), (23), (26)–(29) the reproduction Formula (8).

5. Exact Solutions with Multiplicative Separation of Variables

1 ° . Equation (1) has the multiplicative separable solution [40],
u = t 1 U ( x , y ) ,
which is a special case of solution (12) for α = β = 0 . Here, the function U = U ( x , y ) is described by the stationary Monge–Ampère equation
U x x U y y U x y 2 + U = 0 .
If U = F ( x , y ) is a solution of Equation (30), then the function
U = 1 ( a 1 b 2 a 2 b 1 ) 2 F ( a 1 x + b 1 y + c 1 , a 2 x + b 2 y + c 2 ) ,
where a j , b j , c j ( j = 1 , 2 ) are arbitrary constants, is also a solution of this equation.
Invariant solutions of the reduced Equation (30) are described below in Items 2 ° and 3 ° .
2 ° . Equation (30), in turn, admits the multiplicative separable solution
U = x 2 φ ( y ) ,
where the function φ = φ ( y ) satisfies the autonomous ODE
2 φ φ y y 4 ( φ y ) 2 = φ .
The substitution Z ( φ ) = ( φ y ) 2 reduces this equation to the first-order linear ODE
φ Z φ 4 Z + φ = 0 ,
the general solution of which has the form
Z = C 1 φ 4 + 1 3 φ ,
where C 1 is an arbitrary constant. Replacing Z in (33) with ( φ y ) 2 , we obtain an ODE with separable variables, the general solution of which can be represented in implicit form.
The ODE (32) admits a particular solution φ = 1 12 ( y + C 2 ) 2 , which corresponds to the value C 1 = 0 in (33).
3 ° . Equation (30) has an invariant solution of the form
U = x 2 φ ( ξ ) , ξ = y + γ ln x ,
where γ is an arbitrary constant, and the function φ = φ ( ξ ) satisfies the autonomous ODE
γ φ ξ φ ξ ξ 2 φ φ ξ ξ + 4 ( φ ξ ) 2 φ = 0 .
4 ° . Equation (30) has a generalized separable solution that is quadratic with respect to any independent variable, e.g.,
U = x 2 φ ( y ) + x ψ ( y ) + χ ( y ) ,
where the functions φ = φ ( y ) , ψ = ψ ( y ) , and χ = χ ( y ) are described by the ODE system
2 φ φ y y 4 ( φ y ) 2 = φ , 2 φ ψ y y 4 φ y ψ y = ψ , 2 φ χ y y ( ψ y ) 2 = χ .
The non-invariant solution (34) is a generalization of solution (31). It can be seen that the first equation of system (35) coincides with the nonlinear ODE (32), and the second and third equations are linear with respect to the sought functions.
The following statement is true.
Let φ = φ ( y ) be a solution to the first equation (35). Then, the corresponding general solution to the second equation is given by the formula
ψ = C 1 φ + C 2 y φ ,
where C 1 and C 2 are arbitrary constants.
The proof of this statement is carried out by directly substituting expression (36) into the second equation of system (35) taking into account the fact that the function φ satisfies the first equation (35).
Note that the system (35) admits a particular solution,
φ = 1 12 ( y + A 1 ) 2 , ψ = 0 , χ = | y | [ A 2 cos ( 23 2 ln | y + A 1 | ) + A 3 sin ( 23 2 ln | y + A 1 | ) ] ,
where A 1 , A 2 , and A 2 are arbitrary constants.
Remark 6.
Using the Formula (36), it is possible to obtain a more complex solution of the system (34) than (37) with the same function φ = 1 12 ( y + A 1 ) 2 . In this case, the linear inhomogeneous third ODE of the system (34) is the Euler equation, and its particular solution is sought by the method of indefinite coefficients in the form of a sum of three terms of the power form B j y k j , where B j and k j are the sought ones coefficients.
5 ° . More complex solutions to Equation (30) are discussed further in Section 11 and Section 12.

6. Reductions with Additive and Generalized Separation of Variables Leading to Stationary Monge–Ampère Equations, Exact Solutions

1 ° . Equation (1) has additive separable solutions of the form [32]:
u = A t + w ( x , y ) ,
where A is an arbitrary constant, and the function w is described by the stationary nonhomogeneous Monge–Ampère equation with a constant right-hand side,
w x x w y y w x y 2 = A .
Remark 7.
Solution (38) can be obtained using a linear combination of the Lie symmetries X 3 and X 4 (see Section 2).
2 ° . It is easy to verify that Equation (1) admits an additive separable solution of the form (38), which is expressed in elementary functions
u = C 1 x 2 + C 2 x y + C 3 y 2 + C 4 x + C 5 y + ( 4 C 1 C 3 C 2 2 ) t + C 6 ,
where C 1 , ⋯, C 6 are arbitrary constants.
3 ° . Using the results of [10] one can obtain, for example, the following exact non-invariant solutions of the form (38) of Equation (1):
u = A t ± A C 2 x ( C 1 x + C 2 y ) + φ ( C 1 x + C 2 y ) + C 3 x + C 4 y + C 5 , u = A t + 1 x + C 1 C 2 y 2 + C 3 y + C 3 2 4 C 2 A 12 C 2 ( x 3 + 3 C 1 x 2 ) + C 4 y + C 5 x + C 6 , u = A t ± 2 A 3 C 1 C 2 ( C 1 x C 2 2 y 2 + C 3 ) 3 / 2 + C 4 x + C 5 y + C 6 ,
where C 1 , ⋯, C 6 are arbitrary constants, and φ = φ ( z ) is an arbitrary function.
Remark 8.
For A > 0 , the general solution of the nonhomogeneous Monge–Ampère equation (39) can be represented in parametric form [4,10].
4 ° . Equation (1) admits more complex solutions, than (38), with a generalized separation of variables of the form
u = ( a x + b y + c ) t + w ( x , y ) ,
where a, b, and c are arbitrary constants, and the function w is described by the stationary nonhomogeneous Monge–Ampère equation with a variable right-hand side,
w x x w y y w x y 2 = a x b y c .
For b = c = 0 , Equation (40) has, for example, the following generalized separable solutions:
w = ± 2 3 a y ( a x ) 3 / 2 + C 1 y + φ ( x ) , w = C 1 y 2 + C 2 x y + C 3 y a 12 C 1 x 3 + C 2 2 4 C 1 x 2 + C 4 x + C 5 , w = C 1 y 2 x + C 2 y a 24 C 1 x 4 + C 3 x + C 4 ,
where φ ( x ) is an arbitrary function, and C 1 , …, C 4 are arbitrary constants.

7. Reduction with Generalized Separation of Variables Leading to the Linear Heat Equation

1 ° . Equation (1) admits generalized separable solutions of the form
u = 1 2 a y 2 + b x y + 1 2 c x 2 + d y + ( a c b 2 ) t + U ( x , t ) ,
where a, b, c, and d are arbitrary constants, and the function U = U ( x , t ) is described by the linear heat equation
U t = a U x x .
2 ° . For solutions to Equation (42) under arbitrary initial and boundary conditions, see, for example, in [105,106]. In [106], there is a large list of exact solutions to the linear heat equation, which are expressed in elementary functions. To illustrate what has been said, we present several simple solutions to Equation (42), which contain exponential and trigonometric functions:
U = C 1 exp ( a μ 2 t ± μ x ) + C 2 , U = C 1 exp ( a μ 2 t ) cos ( μ x ) + C 2 , U = C 1 exp ( a μ 2 t ) sin ( μ x ) + C 2 , U = C 1 exp ( μ x ) cos ( μ x 2 a μ 2 t ) + C 2 , U = C 1 exp ( μ x ) sin ( μ x 2 a μ 2 t ) + C 2 ,
where C 1 , C 2 , and μ are arbitrary constants.
Formula (41) for b = c = 0 together with the second and third functions (43) gives two solutions periodic in the spatial variable x, and together with the fourth and fifth functions (43) describes two periodic solutions in time t.

8. Generalized Separable Solutions in the Form of Polynomials in One Spatial Variable

1 ° . Equation (1) has generalized separable solutions that are quadratic with respect to any spatial variable, e.g.,
u = y 2 f ( x , t ) + y g ( x , t ) + h ( x , t ) ,
where the functions f = f ( x , t ) , g = g ( x , t ) , and h = h ( x , t ) are described by the following system of PDEs with two independent variables:
f t = 2 f f x x 4 f x 2 , g t = 2 f g x x 4 f x g x , h t = 2 f h x x g x 2 .
Here the first equation for f is nonlinear and isolated (i.e., does not depend on other equations), and the other two equations are linear with respect to the desired functions g and h (moreover, the equation for g is homogeneous, and the equation for h is nonhomogeneous).
Solution (41), leading to the linear heat equation (42), is a special case of the generalized separable solution (44) with
f ( x , t ) = 1 2 a , g ( x , t ) = b x + d , h ( x , t ) = 1 2 c x 2 + ( a c b 2 ) t + U ( x , t ) .
2 ° . The last two equations of the system (45) have the obvious solution g = C 2 , h = C 3 , where C 2 and C 3 are arbitrary constants. A more general statement is also true.
Statement. Let f = f ( x , t ) be any solution to the first equation (45). Then, the last two equations of the system (45) admit partial solutions
g = C 1 f + C 2 , h = 1 4 C 1 2 f + C 3 ,
where C 1 , C 2 , and C 3 are arbitrary constants.
This statement is proven by substituting functions (46) into the last two equations of system (45) and comparing the resulting expressions with the first equation (45).
3 ° . The simplest solution to the first equation (45) is the constant
f = 1 2 a ,
where a is an arbitrary constant. In this case, the last two PDEs (45) are linear heat equations
g t = a g x x , h t = a h x x g x 2 .
the first of which is homogeneous, and the second is nonhomogeneous. Solutions to these equations can be obtained for arbitrary initial and boundary conditions (see, for example, [105,106]).
Note that the first equation (48) has a degenerate stationary solution g = b x , where b is an arbitrary constant. Some other simple solutions to this equation, which are expressed in elementary functions, are given by the right-hand sides of the Formula (43).
Let us present simple solutions of the traveling wave type to Equation (48) containing exponential functions:
g = C 1 exp ( k x + a k 2 t ) , h = C 2 exp ( k x + a k 2 t ) + C 1 2 2 a exp [ 2 ( k x + a k 2 t ) ] ,
where C 1 , C 2 , and k are arbitrary constants. Substituting (47) and (49) into Formula (44), we obtain an exact solution to the original equation (1).
4 ° . The first equation (45) has the stationary solution
f = 1 C 1 x + C 2 ,
as well as a more complex non-stationary three-parameter exact solution with generalized separation of variables
f = ( x + C 1 ) 2 12 ( t + C 2 ) + C 3 ( t + C 2 ) 1 / 3 ,
where C 1 , C 2 , and C 3 are arbitrary constants.
5 ° . The first equation of system (45) is simplified using the substitution f = 1 / θ , which leads to the following PDE
θ θ t = 2 θ x x ,
that admits a fairly obvious generalized separable solution, linear in t, of the form θ = φ ( x ) t + ψ ( x ) . It follows that the first equation (45) has a solution
f = 1 φ ( x ) t + ψ ( x ) ,
where the functions φ = φ ( x ) and ψ = ψ ( x ) are described by the ODE system
2 φ x x = φ 2 , 2 ψ x x = φ ψ .
System (51) has the following exact solution, which is expressed in elementary functions:
φ = 12 ( x + C 1 ) 2 , ψ = C 2 ( x + C 1 ) 3 + C 3 ( x + C 1 ) 2 ,
where C 1 , C 2 , and C 3 are arbitrary constants.
Let φ = φ ( x ) be a solution to the first ODE of the system (51). Then, the second ODE of this system has a particular solution ψ = φ ( x ) . Since the second ODE is linear with respect to ψ , its general solution is given by the formula [107]:
ψ = C 3 φ + C 3 φ d x φ 2 ,
where C 3 and C 4 are arbitrary constants. The general solution of the first (autonomous) ODE of system (51) is easy to represent in implicit or equivalent parametric form [107]. Taking into account the above and using Formula (51), we obtain the general solution of system (51) in parametric form
x = ± ξ d ζ 1 3 ζ 3 + C 1 1 / 2 + C 2 , φ = ξ , ψ = C 3 ξ + C 4 ξ ξ d ζ 1 3 ζ 3 + C 1 1 / 2 .
where 0 ξ < and C 1 0 .
6 ° . By replacing f = ξ 1 / 2 , the first equation (45) is reduced to a nonlinear heat equation of the special form
ξ t = 2 ( ξ 1 / 2 ξ x ) x ,
which admits a number of invariant and more complex solutions (a review of these solutions is given, for example, in [10,108]).
Remark 9.
In [109,110], a solution to Equation (52) of the form ξ = [ φ ( x ) t + ψ ( x ) ] 2 was obtained, which leads to solution (50).
7 ° . The system of PDEs (45) admits an exact solution of the traveling wave type
f = f ( z ) , g = g ( z ) , h = h ( z ) , z = k x λ t ,
where k and λ are arbitrary constants. The solution of the corresponding autonomous ODE for the function f (which reduces to a linear ODE of the first order) can be expressed in elementary functions.
8 ° . PDE (1) has generalized separable solutions in the form of a fourth-degree polynomial with respect to the spatial variable,
u = y 4 F ( x , t ) + y 2 G ( x , t ) + H ( x , t ) ,
where the functions F = F ( x , t ) , G = G ( x , t ) , and H = H ( x , t ) are described by the following overdetermined PDE system with two independent variables:
3 F F x x 4 F x 2 = 0 , F t = 12 F G x x + 2 G F x x 16 F x G x , G t = 12 F H x x + 2 G G x x 4 G x 2 , H t = 2 G H x x .
Here, the first equation for F is nonlinear and isolated (i.e., does not depend on other equations); it can be considered as an ODE that implicitly depends on t.
Note that the solutions in elementary functions (20), (24), and (25) are particular solutions of the form (53).
The first equation (54) can be satisfied, for example, by setting F = f 0 ( t ) . Then, from the second and other PDEs of this system, we will consistently obtain
G = g 2 x 2 + g 1 x + g 0 , H = h 4 x 4 + h 3 x 3 + h 2 x 2 + h 1 x + h 0 ,
where the functions f 0 = f 0 ( t ) , g i = g i ( t ) , and h j = h j ( t ) are described by a corresponding ODE system.
9 ° . If in Formulas (53) and (55), we remove the odd components in the variable x, setting g 1 = h 3 = h 1 = 0 , then, we arrive at the solution discussed in [42]:
u = f 0 y 4 + h 4 x 4 + g 2 x 2 y 2 + g 0 y 2 + h 2 x 2 + h 0 ,
whose time-dependent functional coefficients are described by the ODE system
f 0 = 24 f 0 g 2 , h 4 = 24 g 2 h 4 , g 2 = 144 f 0 h 4 12 g 2 2 , g 0 = 4 g 0 g 2 + 24 f 0 h 2 , h 2 = 4 g 2 h 2 + 24 g 0 h 4 , h 0 = 4 g 0 h 2 .
The first three ODEs of the system (57) form an independent subsystem, the general solution of which can be represented in closed form (expressed in quadratures). The last three ODEs of this system admit a trivial particular solution g 0 = h 2 = h 0 = 0 .
System (57) has a simple solution of the form
f 0 = a , h 4 = g 2 = 0 , g 0 = 24 a b t + c , h 2 = b , h 0 = 48 a b 2 t 2 + 4 b c t + d ,
where a, b, c, and d are arbitrary constants, which determines the following exact solution of the original equation (1):
u = a y 4 + ( 24 a b t + c ) y 2 + b x 2 + 48 a b 2 t 2 + 4 b c t + d .
Note that this solution was not considered in [42].
It can be shown that the ODE system (57) also has solutions when all the desired functional coefficients in the right part of (56) are proportional to ( t + C ) 1 , where C is an arbitrary constant. In particular, system (57) admits the following solution [42]:
f 0 = A 48 ( t + C ) , h 4 = 1 48 A ( t + C ) , g 2 = 1 24 ( t + C ) , g 0 = h 2 = h 0 = 0 ,
where A and C are arbitrary constants ( A 0 ).
Remark 10.
One can also look for more complex polynomial solutions of the form u = k = 0 4 F k ( x , t ) y k which generalize (53).

9. Reductions to a Monge–Ampère Type Equation in Traveling Wave Variables, Linearizable PDEs, and Two-Phase Solutions

1 ° . Equation (1) admits solutions with generalized separation of variables of the combined type
u = C 1 x 2 + C 2 x y + C 3 y 2 + C 4 x + C 5 y + C 6 t + U ( ξ , η ) , ξ = a 1 x + b 1 y + λ 1 t , η = a 2 x + b 2 y + λ 2 t ,
where C i , a j , b j , λ j ( i = 1 , , 6 ; j = 1 , 2 ) are arbitrary constants, and ξ and η are new traveling wave type variables, and the function U = U ( ξ , η ) is described by a two-dimensional Monge-Ampère type equation
C 6 + λ 1 U ξ + λ 2 U η = 4 C 1 C 3 C 2 2 + ( a 1 b 2 b 1 a 2 ) 2 ( U ξ ξ U η η U ξ η 2 ) + 2 ( a 1 2 C 3 + b 1 2 C 1 a 1 b 1 C 2 ) U ξ ξ + 2 ( a 2 2 C 3 + b 2 2 C 1 a 2 b 2 C 2 ) U η η + 2 [ ( 2 a 1 a 2 C 3 + 2 b 1 b 2 C 1 ( a 1 b 2 + b 1 a 2 ) C 2 ] U ξ η .
2 ° . Equation (59) admits exact solutions, quadratic with respect to any independent variable of the traveling wave type, of the form
U 1 = f 1 ( ξ ) η 2 + g 1 ( ξ ) η + h 1 ( ξ ) , U 2 = f 2 ( η ) ξ 2 + g 2 ( η ) ξ + h 2 ( η ) ,
where the functions f i , g i , h i ( i = 1 , 2 ) are described by the corresponding ODE systems, which are omitted here.
3 ° . Let us consider the special case (58) and (59), putting
a 1 = a , b 1 = b , λ 1 = λ , a 2 = b 2 = 0 , λ 2 = 1 , η = t ,
which corresponds to a solution of the form
u = C 1 x 2 + C 2 x y + C 3 y 2 + C 4 x + C 5 y + C 6 t + U ( ξ , t ) , ξ = a x + b y + λ t ,
where C i ( i = 1 , , 6 ), a, b, and λ are arbitrary constants. In this case, the function U = U ( ξ , t ) is described by the linear convective heat equation with constant source
U t = 2 ( a 2 C 3 + b 2 C 1 a b C 2 ) U ξ ξ λ U ξ + 4 C 1 C 3 C 2 2 C 6 .
Remark 11.
The considered Equation (1), according to the terminology introduced in [111], refers to conditionally integrable PDEs, since it admits solutions that are described by linear PDEs (42) and (61). Perhaps in such cases, it is more accurate to use the term partially linearizable equations.
4 ° . In particular, taking in (60) and (61), the function U with one argument ξ , and setting C 6 = 4 C 1 C 3 C 2 2 , we arrive at a linear ODE with constant coefficients,
2 ( a 2 C 3 + b 2 C 1 a b C 2 ) U ξ ξ λ U ξ = 0 ,
the general solution, which is expressed through the exponent
U ( ξ ) = A 1 exp λ 2 ( a 2 C 3 + b 2 C 1 a b C 2 ) ξ + A 2 ,
where A 1 and A 2 are arbitrary constants.
5 ° . Let us consider a two-phase solution (58) with additive separation of variables of the form
U ( ξ , η ) = V ( ξ ) + W ( η ) , a 1 = a , a 2 = k a , b 1 = b , b 2 = k b ,
where k is an arbitrary constant. As a result, we arrive at a functional differential equation with separable variables,
C 6 4 C 1 C 3 + C 2 2 + λ 1 V ξ + λ 2 W η = 2 ( a 2 C 3 + b 2 C 1 a b C 2 ) V ξ ξ + 2 k 2 ( a 2 C 3 + b 2 C 1 a b C 2 ) W η η ,
which reduces to two simple independent linear ODEs for determining the functions V = V ( ξ ) and W = W ( η ) :
2 ( a 2 C 3 + b 2 C 1 a b C 2 ) V ξ ξ λ 1 V ξ = B , 2 k 2 ( a 2 C 3 + b 2 C 1 a b C 2 ) W η η λ 2 W η = B + C 6 4 C 1 C 3 + C 2 2 ,
where B is an arbitrary constant.
In particular, setting B = 0 , C 6 = 4 C 1 C 3 C 2 2 in (63), we first find V = V ( ξ ) and W = W ( η ) , and then using Formula (62), we determine the required function U ( ξ , η ) included in the solution (58):
U ( ξ , η ) = A 1 exp λ 1 ξ 2 s + A 2 exp λ 2 η 2 k 2 s + A 3 , s = a 2 C 3 + b 2 C 1 a b C 2 , ξ = a x + b y + λ 1 t , η = k ( a x + b y ) + λ 2 t ,
where A 1 , A 2 , and A 3 are arbitrary constants.

10. Reduction Using a New Variable, Parabolic in Spatial Coordinates, Exact Solutions

1 ° . In variables, one of which is time, and the other is specified by a parabolic function in spatial variables,
u = U ( z , t ) , z = y + a x 2 ,
where a is an arbitrary constant, Equation (1) is reduced to the two-dimensional PDE
U t = 2 a U z U z z ,
which is found in the theory of nonlinear filtration of fluids in porous media.
Note that solutions of the form (64) refer to functional separable solutions (see, for example, [10,55]).
Some exact solutions of the reduced Equation (65) are described below.
2 ° . The reduced Equation (65) admits the additive separable solution
U = ± C 1 ( z + C 2 ) 3 / 2 + 9 4 a C 1 2 t + C 3 = ± C 1 ( y + a x 2 + C 2 ) 3 / 2 + 9 4 a C 1 2 t + C 3 ,
where C 1 , C 2 , and C 3 are arbitrary constants.
3 ° . Equation (65) has a simple solution in the form of a product of power functions of different arguments
U = ( z + C 2 ) 3 36 a ( t + C 1 ) = ( y + a x 2 + C 2 ) 3 36 a ( t + C 1 ) ,
where C 1 and C 2 are arbitrary constants.
4 ° . Equation (65) has the traveling wave solution
U = λ 4 a ξ 2 + C 1 ξ + C 2 , ξ = z + λ t y + a x 2 + λ t ,
where C 1 , C 2 , and λ are arbitrary constants.
Remark 12.
More general than (66), the solution of Equation (65) can be obtained if we look for a solution in the form
U = C t + W ( ξ ) , ξ = z + λ t y + a x 2 + λ t .
5 ° . Equation (65) admits the self-similar solution
U = t 3 β 1 V ( ζ ) , ζ = z t β ,
where β is an arbitrary constant, and the function V = V ( ξ ) is described by the non-autonomous ODE
2 a V ζ V ζ ζ = ( 3 β + 1 ) V + β ζ V ζ .
For β = 1 / 3 , the general solution of the last equation has the form
V = 1 36 a ζ 3 + C 1 ζ + C 2 ,
where C 1 and C 2 are arbitrary constants.
6 ° . Equation (65) has an invariant solution of the form
U = t 1 f ( η ) , η = z + λ ln t ,
where λ is an arbitrary constant, and the function f = f ( η ) is described by the autonomous ODE
2 a f η f η η = f + λ f η .
7 ° . Equation (65) also admits another invariant solution
U = e 3 β t g ( τ ) , τ = e β t z ,
where β is an arbitrary constant, and the function g = g ( τ ) is described by the non-autonomous ODE
2 a g τ g τ τ = 3 β g + β τ g τ .
8 ° . The reduced equation (65) admits an exact solution in the form of a cubic polynomial in z:
U = ψ 1 ( t ) + ψ 2 ( t ) z + ψ 3 ( t ) z 2 + ψ 4 ( t ) z 3 ,
where the functions ψ n ( t ) ( n = 1 , , 4 ) are described by the ODE system
ψ 1 = 4 a ψ 2 ψ 3 , ψ 2 = 4 a ( 3 ψ 2 ψ 4 + 2 ψ 3 2 ) , ψ 3 = 36 a ψ 3 ψ 4 , ψ 4 = 36 a ψ 4 2 .
This system is integrated in reverse order, starting with the last equation. Its general solution is described by the formulas
ψ 1 = 12 a C 2 C 3 ( t + C 1 ) 1 / 3 + 48 a 2 C 2 3 ( t + C 1 ) 1 + C 4 , ψ 2 = C 3 ( t + C 1 ) 1 / 3 12 a C 2 2 ( t + C 1 ) 1 , ψ 3 = C 2 ( t + C 1 ) 1 , ψ 4 = 1 36 a ( t + C 1 ) 1 ,
where C 1 , C 2 , C 3 , and C 4 are arbitrary constants.
9 ° . The reduced equation (65) also has an exact solution with a generalized separation of variables of a more exotic form
U = ϑ 1 ( t ) + ϑ 2 ( t ) z 3 / 2 + ϑ 3 ( t ) z 3 ,
where the functions ϑ n ( t ) ( n = 1 , 2 , 3 ) satisfy the ODE system
ϑ 1 = 9 4 a ϑ 2 2 , ϑ 2 = 45 2 a ϑ 2 ϑ 3 , ϑ 3 = 36 a ϑ 3 2 .
This system is integrated in reverse order, starting with the last equation. Its general solution has the form
θ 1 = 9 a C 2 2 ( t + C 1 ) 1 / 4 + C 3 , θ 2 = C 2 ( t + C 1 ) 5 / 8 , θ 3 = 1 36 a ( t + C 1 ) 1 ,
where C 1 , C 2 , and C 3 are arbitrary constants.
Remark 13.
To construct generalized separable solutions (67) and (68), we used invariant subspaces of the nonlinear differential operator F ( U ) = U z U z z , which is included in the right side of Equation (65) (for details see [42,55]).

11. Reduction Using a New Variable, Quadratic in Two Spatial Coordinates, Exact Solutions

1 ° . In variables, one of which is time and the other is quadratic with respect to both spatial variables,
u = U ( z , t ) , z = a x 2 + b x y + c y 2 + k x + s y ,
where a, b, c, k, and s are arbitrary constants, Equation (1) is reduced to a two-dimensional nonstationary equation
U t = 2 ( A z + B ) U z U z z + A U z 2 ; A = 4 a c b 2 , B = a s 2 + c k 2 b k s .
Solutions of the form (69) refer to functional separable solutions (see, for example, [10,55]).
Note that depending on the coefficients a, b, c of the quadratic form in (69), the curve z = const can be an ellipse (for A = 4 a c b 2 > 0 ), a hyperbola (for A < 0 ), or a parabola (for A = 0 ).
The degenerate case A = 0 leads to an equation of the form (65), which was discussed in detail in the previous section. Next, we will consider the non-degenerate case when A = 4 a c b 2 0 .
2 ° . The transformaation
t = t , z = | A | 2 ρ 2 B A , U = sign ( A ) W ( ρ , t ) ,
leads to Equation (70) to the canonical form
W t = ρ 1 W ρ W ρ ρ .
Let us consider some classes of exact solutions that Equation (71) admits.
3 ° . Equation (71) has an exact solution with additive separation of variables
W = C 1 C 2 2 t + C 2 C 1 ρ 2 + C 3 d ρ + C 4 ,
where C 1 , ⋯, C 4 are arbitrary constants. The integral on the right side is expressed in terms of different elementary functions
C 1 ρ 2 + C 3 d ρ = = 1 2 ρ C 1 ρ 2 + C 3 + C 3 2 C 1 ln C 1 ρ + C 1 ρ 2 + C 3 , if C 1 > 0 ; 1 2 ρ C 1 ρ 2 + C 3 + C 3 2 C 1 arctan C 1 ρ C 1 ρ 2 + C 3 , if C 1 < 0 ; C 3 ρ , if C 1 = 0 , C 3 > 0 .
4 ° . Equation (71) admits solutions in the form of a product of functions of different arguments
W = t 1 f ( ρ ) ,
where the function f = f ( ρ ) is described by the non-autonomous ODE
ρ 1 f ρ f ρ = f ,
which has the simple particular solution
f = 1 48 ρ 4 .
Returning now to the original variables, we obtain a five-parameter solution in elementary functions of Equation (1):
u = 1 12 ( 4 a c b 2 ) t a x 2 + b x y + c y 2 + k x + s y + a s 2 + c k 2 b k s 4 a c b 2 2 .
5 ° . Equation (71) admits the self-similar solution
W = t 4 γ 1 F ( z ) , z = t γ ρ ,
where γ is an arbitrary constant, and the function F = F ( z ) is described by the ODE
z 1 F z F z z = ( 4 γ + 1 ) F + γ z F z .
For γ = 1 / 1 , the general solution of the last equation has the form
F = 1 48 z 4 + C 1 z + C 2 ,
where C 1 and C 2 are arbitrary constants.
6 ° . Equation (71) has the invariant solution
W = e 4 λ t Φ ( ζ ) , ζ = ρ e λ t ,
where λ is an arbitrary constant, and the function Φ = Φ ( ζ ) satisfies the generalized homogeneous ODE
ζ 1 Φ ζ Φ ζ ζ = 4 λ Φ + λ ζ Φ ζ ,
whose order can be lowered by one [107].
7 ° . Equation (71) also has the generalized separable solution of polynomial form
W = θ 1 ( t ) + θ 2 ( t ) ρ 2 + θ 3 ( t ) ρ 4 ,
where the functions θ n ( t ) ( n = 1 , 2 , 3 ) satisfy the ODE system
θ 1 = 4 θ 2 2 , θ 2 = 32 θ 2 θ 3 , θ 3 = 48 θ 3 2 .
This system is integrated in reverse order, starting with the last equation. Its general solution is described by the formulas
θ 1 = 12 C 2 2 ( t + C 1 ) 1 / 3 + C 3 , θ 2 = C 2 ( t + C 1 ) 2 / 3 , θ 3 = 1 48 ( t + C 1 ) 1 ,
where C 1 , C 2 , and C 2 are arbitrary constants.

12. Reductions and Exact Solutions in Polar Coordinates

1 ° . At the point ( x 0 , y 0 ) , where x 0 and y 0 are arbitrary constants, we introduce polar coordinates r, φ using the formulas
x = x 0 + r cos φ , y = y 0 + r sin φ .
As a result, the original equation (1) is transformed to the form
u t = r 2 u r r ( u φ φ + r u r ) [ ( r 1 u φ ) r ] 2 .
This equation will be used in the future to construct reductions and exact solutions of the equation under consideration.
2 ° . From the results of group analysis of the transformed magnetohydrodynamics equation (72), it follows that the transformation
r ¯ = a r , φ ¯ = φ + b , t ¯ = p t + q , u ¯ = a 4 p u + c 1 r cos φ + c 2 r sin φ + c 3 ,
where a, b, c 1 , c 2 , c 3 , p, and q are arbitrary constants, leads Equation (72) to an equation of exactly the same form.
The seven-parameter invariant transformation (73) allows using simpler solutions of Equation (72) to construct its more complex exact solutions. Namely, if u = F ( r , φ , t ) is a solution to Equation (72), then the function
u = p a 4 F ( a r , φ + b , p t + q ) + s 1 r cos φ + s 2 r sin φ + s 3 ,
where s 1 = c 1 p / a 4 , s 2 = c 2 p / a 4 , and s 3 = c 3 p / a 4 are arbitrary constants, is also a solution of this equation.
3 ° . Equation (72) admits radially symmetric solutions independent of the angular variable, which are described by the two-dimensional equation
u t = r 1 u r u r r ,
which, up to obvious renotations, coincides with Equation (71). Therefore, it allows five exact solutions, described earlier in Items 3 ° 7 ° from Section 11.
4 ° . Passing in (72) to new variables of mixed type
u = t 4 α 1 U ( ξ , η ) , ξ = r t α , η = φ + β ln t ,
where α and β are arbitrary constants, we arrive at the two-dimensional equation
ξ 2 U ξ ξ ( U η η + ξ U ξ ) [ ( ξ 1 U η ) ξ ] 2 + ( 4 α + 1 ) U α ξ U ξ β U η = 0 .
The values of α = β = 0 in (74) correspond to a solution with multiplicative separation of variables of the form u = t 1 U ( r , φ ) .
The value β = 0 in (74) corresponds to a self-similar solution in the variables r and t.
5 ° . Passing in (72) to other variables of mixed type
u = e 4 γ U ( ξ , η ) , ξ = e γ t r , η = φ λ t ,
where γ and λ are arbitrary constants, we arrive at the two-dimensional equation
ξ 2 U ξ ξ ( U η η + ξ U ξ ) [ ( ξ 1 U η ) ξ ] 2 + 4 γ U γ ξ U ξ + λ U η = 0 .
The value γ = 0 in (75) corresponds to a solution of the traveling wave type in the variables φ and t.
6 ° . Equation (72) also has exact solutions with separation of variables of the form
u = r 4 U ( φ , t ) ,
where the function U = U ( φ , t ) is described by the two-dimensional equation
U t = 12 U U φ φ 9 U φ 2 + 48 U 2 .
7 ° . Since Equation (76) does not depend explicitly on independent variables, it has a solution of the traveling wave type
U = U ( η ) , η = φ λ t ,
where λ is an arbitrary constant, and the function U = U ( η ) is described by the autonomous ODE
12 U U η η 9 ( U η ) 2 48 U 2 + λ U η = 0 .
8 ° . Equation (76) also admits a multiplicative separable solution of the form
U = ( t + C ) 1 V ( φ ) ,
where C is an arbitrary constant, and the function V = V ( φ ) is described by the autonomous ODE
12 V V φ φ 9 ( V φ ) 2 + 48 V 2 + V = 0 .
The substitution Z ( V ) = ( V φ ) 2 leads (77) to a first-order linear ODE 6 V Z V 9 Z + 48 V 2 + V = 0 , the general solution of which is written as follows:
Z = C 1 V 3 / 2 16 V 2 + 1 3 V ,
where C 2 is an arbitrary constant. Integrating further, we find the general solution to Equation (77) in implicit form
d V C 1 V 3 / 2 16 V 2 + 1 3 V 1 / 2 = ± φ + C 2 .
9 ° . Equation (76) has generalized separable solutions of the form
v = f ( t ) + g ( t ) [ A cos ( 4 φ ) + B sin ( 4 φ ) ] ,
where A and B are arbitrary constants, and the functions f ( t ) and g ( t ) are described by the ODE system
f t = 48 f 2 144 ( A 2 + B 2 ) g 2 , g t = 96 f g .
Eliminating the variable t from (78), we arrive at the homogeneous first-order ODE
f g = 1 2 ( f / g ) + 3 2 ( A 2 + B 2 ) ( g / f ) ,
which, using the substitution w = f / g , is reduced to an ODE with separable variables. Integrating this equation, we obtain the general solution to ODE (79) in the form
f 2 = ( A 2 + B 2 ) g 2 + C 1 g ,
where C 1 is an arbitrary constant. Eliminating f from the second ODE (78) using (80), we arrive at an ODE with separable variables, the general solution of which can be represented in implicit form using elementary functions
z + z 2 + C 1 = C 2 exp [ ± 144 ( A 2 + B 2 ) t ] , z = g 3 / 2 .
Using this relation, it is possible to express the function g = g ( t ) explicitly.
Remark 14.
To construct generalized separable solutions of the form (77), we used invariant subspaces of the nonlinear differential operator F [ v ] = 12 v v φ φ 9 v φ 2 + 48 v 2 , included in the right side of Equation (76) (for details see [42,55]).

13. Reductions and Exact Solutions in Elliptic and Hyperbolic Coordinates

1 ° . At the point ( x 0 , y 0 ) , we introduce elliptic coordinates r, φ according to the formulas
x = x 0 + a r cos φ , y = y 0 + b r sin φ ,
where x 0 and y 0 are arbitrary constants, a and b are any positive constants. As a result, the original equation (1) takes the form
u t = ( a b ) 2 r 2 u r r ( u φ φ + r u r ) [ ( r 1 u φ ) r ] 2 .
Using a simple substitution u = ( a b ) 2 u ¯ (or t = ( a b ) 2 t ¯ ), this equation is reduced to Equation (72), the exact solutions of which are described in Section 12.
2 ° . In hyperbolic coordinates ζ , ψ , which are introduced by the formulas
x = x 0 + a ζ cosh ψ , y = y 0 + b ζ sinh ψ ,
where x 0 and y 0 are arbitrary constants, a and b are non-zero constants, Equation (1) takes the form
u t = ( a b ) 2 ζ 2 u ζ ζ ( u ψ ψ + ζ u ζ ) [ ( ζ 1 u ψ ) ζ ] 2 .
Using a simple substitution u = ( a b ) 2 u ¯ , this equation is reduced to Equation (72), the exact solutions of which are described in Section 12.

14. Using a Special Point Transformation to Construct Reductions and Exact Solutions

The special point transformation
x = ξ 1 + α ξ + β η , y = η 1 + α ξ + β η , u = w 1 + α ξ + β η ,
where α and β are arbitrary constants, brings Equation (1) to the form
w ξ ξ w η η w ξ η 2 = ( 1 + α ξ + β η ) 5 w t .
Note that transformation (81) brings the stationary Monge–Ampère equation u x x u y y u x y 2 = f ( x , y ) to a similar form with a different right-hand side [7,9].
Setting β = 0 in (81) and (82), we arrive at the equation
w ξ ξ w η η w ξ η 2 = f ( ξ ) w t ,
where f ( ξ ) = ( 1 + α ξ ) 5 .
Let us now describe some exact solutions of Equation (83) for a more general case, considering the function f ( ξ ) arbitrary.
1 ° . Equation (83) admits generalized separable solutions of the form
w = ( a ξ + b ) t + Z ( ξ , η ) ,
where a and b are arbitrary constants, and the function Z = Z ( ξ , η ) is described by the stationary Monge–Ampère equation
Z ξ ξ Z η η Z ξ η 2 = f ( ξ ) ( a ξ + b ) .
PDEs of this type were considered in [10]. Equation (84) has the following generalized separable solutions in closed form:
Z = ± η f ( ξ ) ( a ξ + b ) d ξ + φ ( ξ ) , Z = C 1 η 2 + C 2 ξ η + C 2 2 4 C 1 ξ 2 + 1 2 C 1 0 ξ ( ξ ζ ) f ( ζ ) ( a ζ + b ) d ζ + C 3 ξ + C 4 η , Z = 1 ξ + C 1 C 2 η 2 + C 3 η + C 3 2 4 C 2 + 1 2 C 2 0 ξ ( ξ ζ ) ( ζ + C 1 ) f ( ζ ) ( a ζ + b ) d ζ ,
where φ ( ξ ) is an arbitrary function, C 1 , …, C 4 are arbitrary constants.
2 ° . Passing in (83) to variables of self-similar type
w = t 2 k 1 W ( ξ , θ ) , θ = η t k ,
where k is a free parameter, we obtain a two-dimensional equation of the Monge–Ampère type with variable coefficients for the lower derivatives
W ξ ξ W θ θ W ξ θ 2 = f ( ξ ) [ ( 2 k + 1 ) W + k θ W θ ] .
3 ° . Passing in (83) to variables of limit self-similar type
w = exp ( 2 γ t ) W ( ξ , θ ) , θ = exp ( γ t ) η ,
where γ is a free parameter, we obtain a two-dimensional stationary equation of the Monge–Ampère type
W ξ ξ W θ θ W ξ θ 2 = γ f ( ξ ) ( 2 W + θ W θ ) .
4 ° . Equation (83) admits a one-dimensional invariant solution
w = t 1 η 2 φ ( ξ ) ,
where the function φ = φ ( ξ ) is described by the ODE
2 φ φ ξ ξ 4 ( φ ξ ) 2 = f ( ξ ) φ .

15. Brief Conclusions

We study a highly nonlinear partial differential equation of the form
u t = u x x u y y u x y 2 ,
which describes local unsteady plasma motions in electron magnetohydrodynamics. Invariant multiparameter transformations that preserve the form of this equation are considered. Some one-dimensional symmetry reductions leading to ordinary differential equations are described and corresponding self-similar and other invariant solutions were obtained. Many new exact solutions with additive, multiplicative, generalized, and functional separation of variables have been found. Special attention is paid to the construction of exact closed-form solutions. More than 30 solutions in elementary functions were obtained (it is such exact solutions that are most useful for testing numerical methods). Two-dimensional symmetry and non-symmetry reductions that reduce the original equation with three independent variables to simpler single PDE or PDE system with two independent variables are described. In particular, a class of solutions is described, which are expressed in terms of solutions to a linear diffusion equation. Note that in some cases, exact solutions were obtained by combining the classical method of symmetry reductions and methods of generalized and/or functional separation of variables (the order in which these methods are applied may vary).
Note that this paper did not use the nonclassical symmetries method and the method of differential constraints. These methods are based on the analysis of the compatibility of overdetermined nonlinear PDE systems and can also be used to construct exact solutions to the highly nonlinear PDE (1) (usually these methods are more labor-intensive and lead to a larger volume of calculations than the methods used in our work).

Author Contributions

Conceptualization and methodology, A.D.P. and A.V.A.; formal analysis and investigation, A.D.P. and A.V.A.; construction of symmetries and group classification, A.D.P. and A.V.A.; writing—original draft preparation, A.D.P. and A.V.A.; writing—review and editing, A.D.P. and A.V.A.; supervision, A.D.P. and A.V.A. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Ministry of Education and Science of the Russian Federation within the framework of the State Assignment (Reg. No. 124012500440-9).

Data Availability Statement

Data are contained within the present article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smirnov, V.V.; Chukbar, K.V. “Phonons” in two-dimensional vortex lattices. J. Exp. Theor. Phys. 2001, 93, 126–135. [Google Scholar] [CrossRef]
  2. Zaburdaev, V.Y.; Smirnov, V.V.; Chukbar, K.V. Nonlinear dynamics of electron vortex lattices. Plasma Phys. Rep. 2014, 30, 214–217. [Google Scholar] [CrossRef]
  3. Ohkitani, K.; Sultu, F.A. Singularity formation in the Smirnov–Chukbar–Zaburdaev equation for the deformation of vortex lattices. J. Phys. A Math. Theor. 2013, 46, 205501. [Google Scholar] [CrossRef]
  4. Goursat, E. A Course of Mathematical Analysis; Part 1; Gostekhizdat: Moscow, Russia, 1933; Volume 3. (In Russian) [Google Scholar]
  5. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: Boston, MA, USA, 1982. [Google Scholar]
  6. Fushchich, W.I.; Serov, N.I. Symmetry and some exact solutions of the multidimensional Monge–Ampère equation. Dokl. Acad. Nauk USSR 1983, 273, 543–546. (In Russian) [Google Scholar]
  7. Khabirov, S.V. Nonisentropic one-dimensional gas motions constructed by means of the contact group of the nonhomogeneous Monge–Ampère equation. Math. Sb. 1990, 181, 1607–1622. Erratum in Math. USSR-Sbornik 1992, 71, 447–462(In Russian) [Google Scholar] [CrossRef]
  8. Fushchich, W.I.; Shtelen, W.M.; Serov, N.I. Symmetry Analysis and Exact Solutions of Equations of Nonlinear Mathematical Physics; Kluwer Academic: Dordrecht, The Netherlands, 1993. [Google Scholar]
  9. Ibragimov, N.H. (Ed.) CRC Handbook of Lie Group Analysis of Differential Equations, Vol. 1, Symmetries, Exact Solutions and Conservation Laws; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
  10. Polyanin, A.D.; Zaitsev, V.F. Handbook of Nonlinear Partial Differential Equations, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2012. [Google Scholar]
  11. Krylov, N.V. Sequences of convex functions and estimates of the maximum of the solution of a parabolic equation. Sib. Math. J. 1976, 17, 226–236. [Google Scholar] [CrossRef]
  12. Spiliotis, J. Certain results on a parabolic type Monge–Ampère equation. J. Math. Anal. Appl. 1992, 163, 484–511. [Google Scholar] [CrossRef]
  13. Ivochkina, N.M.; Ladyzhenskaya, O.A. Parabolic equations generated by symmetric functions of the eigenvalues of the Hessian or by the principal curvatures of a surface. I. Parabolic Monge–Ampère equations. Algebra I Anal. 1994, 6, 141–160. (In Russian) [Google Scholar]
  14. Chen, L.; Wang, G.; Lian, S. Convex-monotone functions and generalized solution of parabolic Monge–Ampère equation. J. Differ. Equ. 2002, 186, 558–571. [Google Scholar] [CrossRef]
  15. Xiong, J.; Bao, J. On Jorgens, Calabi, and Pogorelov type theorem and isolated singularities of parabolic Monge–Ampère equations. J. Differ. Equ. 2011, 250, 367–385. [Google Scholar] [CrossRef]
  16. Tang, L. Regularity results on the parabolic Monge–Ampère equation with VMO type data. J. Differ. Equ. 2013, 255, 1646–1656. [Google Scholar] [CrossRef]
  17. Wang, J.; Yang, J.; Liu, X. The initial and Neumann boundary value problem for a class parabolic Monge–Ampère equation. Abstr. Appl. Anal. 2013, 2013, 535629. [Google Scholar] [CrossRef]
  18. Dai, L. Exterior problems of parabolic Monge–Ampère equations for n = 2. Comput. Math. Appl. 2014, 67, 1497–1506. [Google Scholar] [CrossRef]
  19. Dai, L. Exterior problems for a parabolic Monge–Ampère equation. Nonlinear Anal. Theory Methods Appl. 2014, 100, 99–110. [Google Scholar] [CrossRef]
  20. Tang, L. Boundary regularity on the parabolic Monge–Ampère equation. J. Differ. Equ. 2015, 259, 6399–6431. [Google Scholar] [CrossRef]
  21. Wang, B.; Bao, J. Asymptotic behavior on a kind of parabolic Monge–Ampère equation. J. Differ. Equ. 2015, 259, 344–370. [Google Scholar] [CrossRef]
  22. Zhang, W.; Bao, J.; Wang, B. An extension of Jorgens–Calabi–Pogorelov theorem to parabolic Monge–Ampère equation. Calc. Var. Partial Differ. Equ. 2018, 57, 90. [Google Scholar] [CrossRef]
  23. Dai, L.; Bao, J. Entire solutions of Cauchy problem for parabolic Monge–Ampère equations. Adv. Nonlinear Stud. 2020, 20, 769–781. [Google Scholar]
  24. Karatzas, I. Adaptive control of a diffusion to a goal and a parabolic Monge–Ampère-type equation. Asian J. Math. 1997, 1, 295–313. [Google Scholar] [CrossRef]
  25. Budd, C.J.; Galaktionov, V.A. On self-similar blow-up in evolution equations of Monge–Ampère type. IMA J. Appl. Math. 2013, 78, 338–378. [Google Scholar] [CrossRef]
  26. An, N.; Bao, J.; Liu, Z. Entire solutions to the parabolic Monge–Ampère equation with unbounded nonlinear growth in time. Nonlinear Anal. 2024, 239, 113441. [Google Scholar] [CrossRef]
  27. Ren, C. The first initial-boundary value problem for fully nonlinear parabolic equations generated by functions of the eigenvalues of the Hessian. J. Math. Anal. Appl. 2008, 339, 1362–1373. [Google Scholar] [CrossRef]
  28. Dai, L.; Guo, X. Parabolic Hessian equations outside a cylinder. Mathematics 2022, 10, 2839. [Google Scholar] [CrossRef]
  29. Pogorelov, A.V. Extrinsic Geometry of Convex Surfaces; American Math. Soc.: Providence, RI, USA, 1973. [Google Scholar]
  30. Leibov, O. Reduction and exact solutions of the Monge–Ampère equation. J. Nonlinear Math. Phys. 1997, 4, 146–148. [Google Scholar] [CrossRef]
  31. Arrigo, D.J.; Hill, J.M. On a class of linearizable Monge–Ampère equations. J. Nonlinear Math. Phys. 1998, 5, 115–119. [Google Scholar] [CrossRef]
  32. Rakhmelevich, I.V. On solutions of the Monge–Ampère equation with power-law non-linearity with respect to first derivatives. Tomsk State Univ. J. Math. Mech. 2016, 4, 33–43. [Google Scholar] [CrossRef]
  33. Feroze, T.; Umair, M. Optimal system and exact solutions of Monge–Ampère equation. Commun. Math. Appl. 2021, 12, 825–833. [Google Scholar] [CrossRef]
  34. Polyanin, A.D. Handbook of Exact Solutions to Mathematical Equations; CRC Press: Boca Raton, FL, USA; London, UK, 2024. [Google Scholar]
  35. Courant, R.; Hilbert, D. Methods of Mathematical Physics, Vol. 2; Wiley–Intersci. Publ.: New York, NY, USA, 1989. [Google Scholar]
  36. Sulman, M.M.; Williams, J.F.; Russell, R.D. An efficient approach for the numerical solution of the Monge–Ampère equation. Appl. Numer. Math. 2011, 61, 298–307. [Google Scholar] [CrossRef]
  37. Feng, X.; Lewis, T. Nonstandard local discontinuous Galerkin methods for fully nonlinear second order elliptic and parabolic equations in high dimensions. J. Sci. Comput. 2018, 77, 1534–1565. [Google Scholar] [CrossRef]
  38. Liu, H.; Glowinski, R.; Leung, S.; Qian, J. A finite element/operator-splitting method for the numerical solution of the three dimensional Monge–Ampère equation. J. Sci. Comput. 2019, 81, 2271–2302. [Google Scholar] [CrossRef]
  39. Caboussat, A.; Gourzoulidis, D.; Picasso, D. An adaptive method for the numerical solution of a 2D Monge–Ampère equation. In Proceedings of the 10th International Conference on Adaptive Modeling and Simulation (ADMOS 2021), Gothenburg, Sweden, 21–23 June 2021. [Google Scholar]
  40. Dubinov, A.E.; Kitayev, I.N. New exact solutions of the equation of non-linear dynamics of a lattice of electronic vortices in plasma in the framework of electron magnetohydrodynamics. Magnetohydrodynamics 2020, 56, 369–375. [Google Scholar]
  41. Rakhmelevich, I.V. Nonautonomous evolution equation of Monge–Ampère type with two space variables. Russ. Math. 2023, 67, 52–64. [Google Scholar] [CrossRef]
  42. Galaktionov, V.A.; Svirshchevskii, S.R. Exact Solutions and Invariant Subspaces of Nonlinear Partial Differential Equations in Mechanics and Physics; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  43. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  44. Olver, P.J. Application of Lie Groups to Differential Equations, 2nd ed.; Springer: New York, NY, USA, 2000. [Google Scholar]
  45. Andreev, V.K.; Kaptsov, O.V.; Pukhnachov, V.V.; Rodionov, A.A. Applications of Group-Theoretical Methods in Hydrodynamics; Kluwer: Dordrecht, The Netherlands, 1998. [Google Scholar]
  46. Vaneeva, O.O.; Johnpillai, A.G.; Popovych, R.O.; Sophocleous, C. Extended group analysis of variable coefficient reaction-diffusion equations with power nonlinearities. J. Math. Anal. Appl. 2007, 330, 1363–1386. [Google Scholar] [CrossRef]
  47. Vaneeva, O.O.; Popovych, R.O.; Sophocleous, C. Extended group analysis of variable coefficient reaction-diffusion equations with exponential nonlinearities. J. Math. Anal. Appl. 2012, 396, 225–242. [Google Scholar] [CrossRef]
  48. Shagolshem, S.; Bira, B.; Sil, S. Application of symmetry analysis to viscoelastic fluid model. Commun. Nonlinear Sci. Numer. Simulat. 2023, 125, 107417. [Google Scholar] [CrossRef]
  49. Clarkson, P.A.; Kruskal, M.D. New similarity reductions of the Boussinesq equation. J. Math. Phys. 1989, 30, 2201–2213. [Google Scholar] [CrossRef]
  50. Nucci, M.C.; Clarkson, P.A. The nonclassical method is more general than the direct method for symmetry reductions. An example of the Fitzhugh–Nagumo equation. Phys. Lett. A 1992, 164, 49–56. [Google Scholar] [CrossRef]
  51. Olver, P.J. Direct reduction and differential constraints. Proc. R. Soc. Lond. Ser. A 1994, 444, 509–523. [Google Scholar]
  52. Clarkson, P.A.; Ludlow, D.K.; Priestley, T.J. The classical, direct and nonclassical methods for symmetry reductions of nonlinear partial differential equations. Meth. Appl. Anal. 1997, 4, 173–195. [Google Scholar] [CrossRef]
  53. Hood, S. On direct, implicit reductions of a nonlinear diffusion equation with an arbitrary function—generalizations of Clarkson’s and Kruskal’s method. IMA J. Appl. Math. 2000, 64, 223–244. [Google Scholar] [CrossRef]
  54. Polyanin, A.D.; Zhurov, A.I. Direct functional separation of variables and new exact solutions to axisymmetric unsteady boundary-layer equations. Commun. Nonlinear Sci. Numer. Simulat. 2016, 31, 11–20. [Google Scholar] [CrossRef]
  55. Polyanin, A.D.; Zhurov, A.I. Separation of Variables and Exact Solutions to Nonlinear PDEs; CRC Press: Boca Raton, FL, USA; London, UK, 2022. [Google Scholar]
  56. Bluman, G.W.; Cole, J.D. The general similarity solution of the heat equation. J. Math. Mech. 1969, 18, 1025–1042. [Google Scholar]
  57. Levi, D.; Winternitz, P. Nonclassical symmetry reduction: Example of the Boussinesq equation. J. Phys. A 1989, 22, 2915–2924. [Google Scholar] [CrossRef]
  58. Pucci, E. Similarity reductions of partial differential equations. J. Phys. A Math. Gen. 1992, 25, 2631–2640. [Google Scholar] [CrossRef]
  59. Arrigo, D.; Broadbridge, P.; Hill, J.M. Nonclassical symmetry solutions and the methods of Bluman–Cole and Clarkson–Kruskal. J. Math. Phys. 1993, 34, 4692–4703. [Google Scholar] [CrossRef]
  60. Pucci, E.; Saccomandi, G. Evolution equations, invariant surface conditions and functional separation of variables. Physica D 2000, 139, 28–47. [Google Scholar] [CrossRef]
  61. Saccomandi, G. A remarkable class of non-classical symmetries of the steady two-dimensional boundary-layer equations. J. Phys. A Math. Gen. 2004, 37, 7005–7017. [Google Scholar] [CrossRef]
  62. Bradshaw-Hajek, D.H. Nonclassical symmetry solutions for non-autonomous reaction-diffusion equations. Symmetry 2019, 11, 208. [Google Scholar] [CrossRef]
  63. Pucci, E.; Saccomandi, G. Partial differential equations admitting a given nonclassical point symmetry. Stud. Appl. Math. 2020, 145, 81–96. [Google Scholar] [CrossRef]
  64. Cherniha, R.; Davydovych, V. Nonlinear Reaction-Diffusion Systems: Conditional Symmetry, Exact Solutions and Their Applications in Biology; Springer: Cham, Switzerland, 2017. [Google Scholar]
  65. Cherniha, R.; Davydovych, V.; King, J.R. The Shigesada–Kawasaki–Teramoto model: Conditional symmetries, exact solutions and their properties. Commun. Nonlinear Sci. Numer. Simulat. 2023, 124, 107313. [Google Scholar] [CrossRef]
  66. Galaktionov, V.A. Invariant subspaces and new explicit solutions to evolution equations with quadratic nonlinearities. Proc. R. Soc. Edinburgh Sect. A 1995, 125, 225–246. [Google Scholar] [CrossRef]
  67. Galaktionov, V.A.; Posashkov, S.A. Exact-solutions and invariant subspace for nonlinear gradient-diffusion equations. Comp. Math. Math. Phys. 1994, 34, 313–321. [Google Scholar]
  68. Galaktionov, V.A.; Posashkov, S.A.; Svirshchevskii, S.R. Generalized separation of variables for differential equations with polynomial nonlinearities. Diff. Equ. 1995, 31, 233–240. [Google Scholar]
  69. Polyanin, A.D. Exact solutions to the Navier–Stokes equations with generalized separation of variables. Dokl. Phys. 2001, 46, 726–731. [Google Scholar] [CrossRef]
  70. Polyanin, A.D.; Zhurov, A.I. Unsteady axisymmetric boundary-layer equations: Transformations, properties, exact solutions, order reduction and solution method. Int. J. Non-Linear Mech. 2015, 74, 40–50. [Google Scholar] [CrossRef]
  71. Kosov, A.A.; Semenov, E.I. Exact solutions of the generalized Richards equation with power-law nonlinearities. Diff. Equ. 2020, 56, 1119–1129. [Google Scholar] [CrossRef]
  72. Kosov, A.A.; Semenov, E.I. Reduction method and new exact solutions of the multidimensional nonlinear heat equation. Diff. Equ. 2022, 58, 187–194. [Google Scholar] [CrossRef]
  73. Svirshchevskii, S.R. Exact solutions of a nonlinear diffusion equation on polynomial invariant subspace of maximal dimension. Commun. Nonlinear Sci. Numer. Simulat. 2022, 112, 106515. [Google Scholar] [CrossRef]
  74. Grundland, A.M.; Infeld, E. A family of non-linear Klein-Gordon equations and their solutions. J. Math. Phys. 1992, 33, 2498–2503. [Google Scholar] [CrossRef]
  75. Miller, W.; Rubel, L.A. Functional separation of variables for Laplace equations in two dimensions. J. Phys. A 1993, 26, 1901–1913. [Google Scholar] [CrossRef]
  76. Zhdanov, R.Z. Separation of variables in the non-linear wave equation. J. Phys. A 1994, 27, L291–L297. [Google Scholar] [CrossRef]
  77. Doyle, P.W.; Vassiliou, P.J. Separation of variables for the 1-dimensional non-linear diffusion equation. Int. J. Non-Linear Mech. 1998, 33, 315–326. [Google Scholar] [CrossRef]
  78. Estevez, P.G.; Qu, C.; Zhang, S. Separation of variables of a generalized porous medium equation with nonlinear source. J. Math. Anal. Appl. 2002, 275, 44–59. [Google Scholar] [CrossRef]
  79. Polyanin, A.D. Construction of exact solutions in implicit form for PDEs: New functional separable solutions of non-linear reaction-diffusion equations with variable coefficients. Int. J. Non-Linear Mech. 2019, 111, 95–105. [Google Scholar] [CrossRef]
  80. Polyanin, A.D. Functional separable solutions of nonlinear reaction-diffusion equations with variable coefficients. Appl. Math. Comput. 2019, 347, 282–292. [Google Scholar] [CrossRef]
  81. Polyanin, A.D. Functional separable solutions of nonlinear convection-diffusion equations with variable coefficients. Commun. Nonlinear Sci. Numer. Simulat. 2019, 73, 379–390. [Google Scholar] [CrossRef]
  82. Polyanin, A.D. Construction of functional separable solutions in implicit form for non-linear Klein–Gordon type equations with variable coefficients. Int. J. Non-Linear Mech. 2019, 114, 29–40. [Google Scholar] [CrossRef]
  83. Polyanin, A.D.; Zhurov, A.I. Separation of variables in PDEs using nonlinear transformations: Applications to reaction-diffusion type equations. Appl. Math. Lett. 2020, 100, 106055. [Google Scholar] [CrossRef]
  84. Sidorov, A.F.; Shapeev, V.P.; Yanenko, N.N. Method of Differential Constraints and its Applications in Gas Dynamics; Nauka: Novosibirsk, Russia, 1984. (In Russian) [Google Scholar]
  85. Kaptsov, O.V.; Verevkin, L.V. Differential constraints and exact solutions of nonlinear diffusion equations. J. Phys. A Math. Gen. 2003, 36, 1401–1414. [Google Scholar] [CrossRef]
  86. Meleshko, S.V. Methods for Constructing Exact Solutions of Partial Differential Equations; Springer: New York, NY, USA, 2005. [Google Scholar]
  87. Kruglikov, B. Symmetry approaches for reductions of PDEs, differential constraints and Lagrange–Charpit method. Acta Appl. Math. 2008, 101, 145–161. [Google Scholar] [CrossRef]
  88. Calogero, F.; Degasperis, A. Spectral Transform and Solitons: Tolls to Solve and Investigate Nonlinear Evolution Equations; North-Holland Publ.: Amsterdam, The Netherlands, 1982. [Google Scholar]
  89. Ablowitz, M.J.; Clarkson, P.A. Solitons, Nonlinear Evolution Equations and Inverse Scattering; Cambridge Univ. Press: Cambridge, UK, 1991. [Google Scholar]
  90. Kudryashov, N.A. Simplest equation method to look for exact solutions of nonlinear differential equations. Chaos Solitons Fractals 2005, 24, 1217–1231. [Google Scholar] [CrossRef]
  91. Kudryashov, N.A. Methods of Nonlinear Mathematical Physics; Izd. Dom Intellekt: Dolgoprudnyi, Russia, 2010. (In Russian) [Google Scholar]
  92. Conte, R.; Musette, M. The Painlevé Handbook, 2nd ed.; Springer: Cham, Switzerland, 2020. [Google Scholar]
  93. Aksenov, A.V.; Polyanin, A.D. Methods for constructing complex solutions of nonlinear PDEs using simpler solutions. Mathematics 2021, 9, 345. [Google Scholar] [CrossRef]
  94. Aksenov, A.V.; Polyanin, A.D. Review of methods for constructing exact solutions of equations of mathematical physics based on simpler solutions. Theor. Math. Phys. 2022, 211, 567–594. [Google Scholar] [CrossRef]
  95. Bedrikovetsky, P.; Borazjani, S. Exact solutions for gravity-segregated flows in porous media. Mathematics 2022, 10, 2455. [Google Scholar] [CrossRef]
  96. Meleshko, S.V.; Moyo, S. On the complete group classification of the reaction–diffusion equation with a delay. J. Math. Anal. Appl. 2008, 338, 448–466. [Google Scholar] [CrossRef]
  97. Long, F.-S.; Meleshko, S.V. On the complete group classification of the one-dimensional nonlinear Klein–Gordon equation with a delay. Math. Methods Appl. Sci. 2016, 39, 3255–3270. [Google Scholar] [CrossRef]
  98. Polyanin, A.D.; Sorokin, V.G. A method for constructing exact solutions of nonlinear delay PDEs. J. Math. Anal. Appl. 2021, 494, 124619. [Google Scholar] [CrossRef]
  99. Polyanin, A.D.; Sorokin, V.G. Nonlinear pantograph-type diffusion PDEs: Exact solutions and the principle of analogy. Mathematics 2021, 9, 511. [Google Scholar] [CrossRef]
  100. Aibinu, M.O.; Thakur, S.C.; Moyo, S. Exact solutions of nonlinear delay reaction–diffusion equations with variable coefficients. Partial Diff. Equ. Appl. Math. 2021, 4, 100170. [Google Scholar] [CrossRef]
  101. Sorokin, V.G.; Vyazmin, A.V. Nonlinear reaction–diffusion equations with delay: Partial survey, exact solutions, test problems, and numerical integration. Mathematics 2022, 10, 1886. [Google Scholar] [CrossRef]
  102. Polyanin, A.D.; Sorokin, V.G. Reductions and exact solutions of nonlinear wave-type PDEs with proportional and more complex delays. Mathematics 2023, 11, 516. [Google Scholar] [CrossRef]
  103. Polyanin, A.D.; Sorokin, V.G. Exact solutions of reaction-diffusion PDEs with anisotropic time delay. Mathematics 2023, 11, 3111. [Google Scholar] [CrossRef]
  104. Polyanin, A.D.; Sorokin, V.G.; Zhurov, A.I. Delay Ordinary and Partial Differential Equations; CRC Press: Boca Raton, FL, USA, 2024. [Google Scholar]
  105. Tikhonov, A.N.; Samarskii, A.A. Equations of Mathematical Physics; Dover Publ.: New York, NY, USA, 1990. [Google Scholar]
  106. Polyanin, A.D.; Nazaikinskii, V.E. Handbook of Linear Partial Differential Equations for Engineers and Scientists, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  107. Polyanin, A.D.; Zaitsev, V.F. Handbook of Ordinary Differential Equations: Exact Solutions, Methods, and Problems; CRC Press: Boca Raton, FL, USA; London, UK, 2018. [Google Scholar]
  108. Ivanova, N.M. Exact solutions of diffusion-convection equations. Dyn. PDEs 2008, 5, 139–171. [Google Scholar] [CrossRef]
  109. Amerov, T.K. On conditional invariance of nonlinear heat equation. In Theoretical and Algebraic Analysis of Equations of Mathematical Physics; Inst. of Mathematics Kiev: Kiev, Ukraine, 1990; pp. 12–14. [Google Scholar]
  110. King, J.R. Some non-self-similar solutions to a nonlinear diffusion equations. J. Phys. A Math. Gen. 1992, 25, 4861–4868. [Google Scholar] [CrossRef]
  111. Broadbridge, P.; Bradshaw-Hajek, B.H.; Hutchinson, A.J. Conditionally integrable PDEs, non-classical symmetries and applications. Proc. R. Soc. A 2023, 479, 20230209. [Google Scholar] [CrossRef]
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

Polyanin, A.D.; Aksenov, A.V. Unsteady Magnetohydrodynamics PDE of Monge–Ampère Type: Symmetries, Closed-Form Solutions, and Reductions. Mathematics 2024, 12, 2127. https://doi.org/10.3390/math12132127

AMA Style

Polyanin AD, Aksenov AV. Unsteady Magnetohydrodynamics PDE of Monge–Ampère Type: Symmetries, Closed-Form Solutions, and Reductions. Mathematics. 2024; 12(13):2127. https://doi.org/10.3390/math12132127

Chicago/Turabian Style

Polyanin, Andrei D., and Alexander V. Aksenov. 2024. "Unsteady Magnetohydrodynamics PDE of Monge–Ampère Type: Symmetries, Closed-Form Solutions, and Reductions" Mathematics 12, no. 13: 2127. https://doi.org/10.3390/math12132127

APA Style

Polyanin, A. D., & Aksenov, A. V. (2024). Unsteady Magnetohydrodynamics PDE of Monge–Ampère Type: Symmetries, Closed-Form Solutions, and Reductions. Mathematics, 12(13), 2127. https://doi.org/10.3390/math12132127

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