Next Article in Journal
A Prototype of Speech Interface Based on the Google Cloud Platform to Access a Semantic Website
Previous Article in Journal
Cryptanalysis of an Image Encryption Algorithm Based on Combined Chaos for a BAN System, and Improved Scheme Using SHA-512 and Hyperchaos
Previous Article in Special Issue
Some Approaches to the Calculation of Conservation Laws for a Telegraph System and Their Comparisons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Second-Order Conditional Lie–Bäcklund Symmetries and Differential Constraints of Nonlinear Reaction–Diffusion Equations with Gradient-Dependent Diffusivity

1
Department of Information and Computational Science, He’nan Agricultural University, Zhengzhou 450002, China
2
Department of Applied Mathematics, Zhejiang University of Technology, Hangzhou 310023, China
*
Author to whom correspondence should be addressed.
Symmetry 2018, 10(7), 267; https://doi.org/10.3390/sym10070267
Submission received: 23 April 2018 / Revised: 1 July 2018 / Accepted: 2 July 2018 / Published: 7 July 2018

Abstract

:
The radially symmetric nonlinear reaction–diffusion equation with gradient-dependent diffusivity is investigated. We obtain conditions under which the equations admit second-order conditional Lie–Bäcklund symmetries and first-order Hamilton–Jacobi sign-invariants which preserve both signs (≥0 and ≤0) on the solution manifold. The corresponding reductions of the resulting equations are established due to the compatibility of the invariant surface conditions and the governing equations.

1. Introduction

The classical symmetry [1] for a general system of partial differential equations (PDEs) is originally created by Sophus Lie more than 100 years ago. A Lie group of transformation depends on continuous parameters and consists of either point transformations (point symmetries) acting on the independent and dependent variables of the system. The Lie group is completely characterized by the infinitesimal generator and the determining equation for the infinitesimal generator of such group is a linear homogeneous system of PDEs. The setting up of the determining equation for a given system of PDEs is entirely a matter of routine.
Over the years, many generalizations of the classical symmetry have been proposed. The possibility of the existence of Lie–Bäcklund symmetry (also known as generalized symmetry or higher-order symmetry) appears to have been first considered by Noether [2] and is characterized by an infinitesimal generator that acts on independent variables, dependent variables and their derivatives to some finite order. The procedure for constructing Lie–Bäcklund symmetries is about the same as that for the classical symmetry and the Lie–Bäcklund symmetries admitted by the system in question are found by solving a linear system of PDEs called the determining system. The method of conditional symmetry can be referred to [3,4,5] and the references therein.
Zhdanov [6], Fokas and Liu [7] independently introduced the method of conditional Lie–Bäcklund symmetry (CLBS) for scalar evolution equations. The reductions of the initial value problem for an evolution equation due to CLBS are considered by Zhdanov et al. in [8,9,10]. The studies of CLBS for an evolution system are discussed by Sergyeyev [11] and independently by Ji et al. in [12,13]. The CLBS can be regarded as a natural generalization of the conditional symmetry in a similar way to how the Lie–Bäcklund symmetry is a generalization of the classical symmetry. Therefore, the procedure for computing CLBS is about the same as that for the conditional symmetry. The crucial problem is that we need to determine the form of CLBS presumably. It has been proved that CLBS with the characteristic
η = f ( u ) l x + a 1 ( x ) f ( u ) ( l 1 ) x + + a l 1 ( x ) f ( u ) x + a l ( x ) f ( u )
is powerful to study classifications and reductions of second-order nonlinear diffusion equations [14,15,16]. In fact, the linear CLBS with the characteristic
σ 1 = v l x + a 1 ( x ) v ( l 1 ) x + + a l 1 ( x ) v x + a l ( x ) v
is responsible for the invariant subspace [17,18]
W l = W f 1 ( x ) , f 2 ( x ) , , f l ( x ) j = 1 l C j f j ( x ) , C j R
admitted by the considered equation. It is noted that f ( u ) j x = j f ( u ) x j , v j x = j v x j and j = 1 , 2 , , l . The CLBSs with the characteristic
η = u x x + H ( u ) u x 2 + G ( x , u ) u x + F ( x , u )
and
η = u x x + H ( u ) u x 2 + G ( x , u ) u x 2 m + F ( x , u ) u x 1 m
of nonlinear diffusion equations are considered in [19,20,21,22,23] and are related to first-order Hamilton–Jacobi sign-invariants (H–J SIs) of second-order nonlinear equations in question. In addition, the CLBS with the characteristic
η = [ f ( u ) ] x t
and
η = [ f ( u , u x ) ] x t
can respectively be used to construct functionally separable solutions [24] and derivative-dependent functionally separable solutions [25] of nonlinear diffusion equations. A large number of exact solutions for different types of nonlinear diffusion equations have been constructed due to the compatibility of the invariant surface conditions η = 0 and the governing equations. The base of the reductions of CLBS is due to the compatibility of the invariant surface condition η = 0 and the governing equation, that is, the invariant surface condition η = 0 is nothing but a differential constraint (DC) of the original equation.
The key idea of the method of DC is proposed by Yanenko in [26]. The survey of this method is presented by Sidorvo, Shapeev and Yanenko in [27], where the method of DC is successfully introduced into practice on gas dynamics. The general formulation of this method requires that the original differential equation be enlarged by appending additional equations such that the over-determined system satisfies some conditions of compatibility. The procedure to determine whether or not a given differential equation is compatible with the original equation is straightforward; however, for a given equation, the determination of the most general DC that is compatible with the original equation is a very difficult, if not impossible, problem. Therefore, it is better to make oneself content with finding constraints in some classes, and these classes must be chosen using additional considerations. It is proved that DCs corresponding to the invariant surface conditions of the above CLBSs related to invariant subspace [14,15,16], sign-invariant [23,28] and separation of variables [24,25] are all very effective to study classifications and reductions of nonlinear diffusion equations.
It is also interesting to note that second-order CLBSs of nonlinear diffusion equations can be translated to first-order H–J SIs of the considered equations [23,28]. The property of SI is a junction of qualitative and quantitative properties of parabolic equations. In [18], the SIs for second-order parabolic equations are extensively discussed from many points of view, including constructing SIs of the considered equations, deriving partial differential inequalities with physical interpretations and discussing the geometrical and physical properties of the exact solutions obtained.
In this paper, we are mainly concerned with the second-order CLBS
η = u r r + H ( u ) u r 2 + G ( r , u ) u r + F ( r , u )
and the first-order H–J SI
J = u t A ( u ) u r m + 1 B ( r , u ) u r m C ( r , u ) u r m 1 E ( r , u )
as well as the exact solutions of the radially symmetric nonlinear reaction–diffusion equation with gradient-dependent diffusivity
u t = 1 r n 1 [ r n 1 D ( u ) u r m ] r + Q ( r , u ) ,
where D ( u ) = u k and Q ( r , u ) are respectively the diffusion and reaction term. It is noted that m, k in (3) are constants, which may be positive or negative. This equation has been applied to describe several situations such as nonlinear heat conduction and nonlinear shear flows of nonlinear Newtonian fluids [29,30,31] (see more details, e.g., in [32]). The seminal work about the application of the nonlinear Equation (3) with D = 1 and Q = 0 to the turbulent gas flows can be referred to [33].
Lie symmetry of the simplest case of Equation (3) with n = m = 1 , D = 1 and Q = 0 , i.e., u t = u r r are discussed in [1] and [34]. Recently, Cherniha et al. [32] described all possible Lie symmetries of Equation (3) in the case n = 1 and D = 1 . For the case of m = 1 and n = 1 , Equation (3) has been studied extensively from many different points of view and a huge number of results on symmetries, exact solutions and qualitative analysis including existence, uniqueness, blow up, extinction, large time behavior and geometrical properties have been obtained (see [35,36,37,38] and references therein). Lie symmetries and reductions of the heat equation with a source (sink)
u t = d i v D ( u ) u
in the two-dimensional and three-dimensional cases are studied in [39]. The nonclassical symmetries and reductions of Equation (4) in two spatial dimensions can be referred to [40]. The nonclassical symmetries and reductions of Arrhenius reaction–diffusion in n dimensions are considered in [41].
It has been known that the approaches of SI [42,43,44] and CLBS [14,15,16,17,18,19,20,21,22,23,24,25] have been successfully applied to obtain solutions and to explore the properties of Equation (3) with n = 1 or m = 1 . Galaktionov [42,43,44] studied the existence of H–J SIs to the 1 + 1 -dimensional nonlinear diffusion equations and the higher-dimensional case. He showed that Equation (3) with n = 1 and m = 1 admits the H–J SI
J = u t + A ( u ) u r 2 + B ( u ) u r + C ( u )
for some coefficient functions. He also verified that the n + 1 -dimensional nonlinear reaction–diffusion equation with gradient-dependent diffusivity
u t = d i v ( | u | σ u ) + f ( u ) , t > 0
admits the H–J SIs with the form
J = u t h ( u ) .
Thus, it is of great interest to develop both methods to study Equation (3) with n 1 and m 1 . Notably, all possible Lie symmetries of Equation (6) for the case of u possessing two space variables were found in [32].
The remainder of this paper is organized as follows. In the following section, we review some necessary notations, definitions and fundamental theorems on CLBS, DC and H–J SI. The corresponding conditional invariance condition for (3) is also presented here. In Section 3, second-order CLBS (1) and first-order H–J SI (2) of Equation (3) are identified. The corresponding reductions of the resulting Equation (3) are presented in Section 4. The last section is devoted to conclusions and remarks.

2. Basic Notations and Theorems

Let us review some basic facts on the CLBS method, the DC method and the H–J SI method of a nonlinear evolution equation.
Let
V = k = 0 D r k η ( r , t , u , u r , , u l r ) u k r
be an evolutionary vector field with the characteristic η and
u t = E ( r , t , u , u r , , u n r )
be a nonlinear evolution equation, where we use the following notations
D r = r + k = 0 u ( k + 1 ) r u k r , D r j + 1 = D r ( D r j ) , D r 0 = 1 , u k r = k u r k .
Definition 1 ([45,46,47]).
The evolutionary vector field (8) is said to be a Lie–Bäcklund symmetry of (9) if and only if
V ( u t E ) | L = 0 ,
where L is the set of all differential consequences of the equation, that is, D r j D t k ( u t E ) = 0 , j , k = 0 , 1 , 2 , .
Definition 2.
[6,7] The evolutionary vector field (8) is said to be a CLBS of (9) if and only if
V ( u t E ) | L M r = 0 ,
where M r denotes the set of all differential consequences of equation η = 0 with respect to r, that is, D r j η = 0 , j = 0 , 1 , 2 , .
The above conditional invariance criterion can be reduced to
D t η | L M r = 0 .
The fact that Lie–Bäcklund vector field (8) is conditional invariant with respect to evolution Equation (9) yields the compatibility of the invariant surface condition
η ( r , t , u , u r , , u l r ) = 0
and the governing Equation (9).
Definition 3.
[48] The differential constraint (11) and the evolution system (9) satisfy the compatibility condition if
D t η | L r M r = 0 ,
where L r denotes the set of all differential consequences of the Equation (9) with respect to r, that is, D r j ( u t E ) = 0 , j = 0 , 1 , 2 , .
Olver [49] put forward the view that conditional symmetries and their further generalizations are connected with DC. For evolution Equation (9), the compatibility condition (12) is nothing but the conditional invariance criterion (10).
Let
J [ u ] = J ( r , t , u , u t , u r )
be sufficiently smooth; then, J [ · ] is a general first-order Hamilton–Jacobi operator.
Definition 4.
[42]. (13) is said to be a H–J SI of (9) if it preserves both signs ≥0 and ≤0 on the solution manifold of (9). This means
J [ u ] 0 ( r e s p . 0 ) i n R f o r t = 0 ,
J [ u ] 0 ( r e s p . 0 ) i n R f o r t > 0 .
This property is understood in the sense of a formal application of the Maximum Principle of parabolic equations. It has been applied to obtain a large number of new exact solutions to nonlinear diffusion equations with source terms [42,43,44]. The equivalence classes of Lie–Bäcklund transformation [45] shows that the second-order CLBS (1) of Equation (3) degenerates into conditional symmetry if the corresponding first-order H–J SI (2) is of the form
J = u t R ( r , u ) u r E ( r , u ) .
Here, we do not consider these degenerated CLBS.
Since (3) admits the CLBS with the characteristic (1), a direct computation of the conditional invariance condition (10) will yield that
D t η | L M r [ D ( 3 m + 1 ) H D 3 m H D + m ( 3 m + 2 ) H 2 D m D H m 2 ( m + 1 ) D H 3 + m ( 3 m + 1 ) D H H ] u r m + 3 + { n 1 r D ( 2 m 1 ) H D + m ( m 1 ) D H 2 ( m 1 ) D H ( 3 m + 2 ) G D + m [ 3 ( 2 m + 1 ) H D G + ( 3 m + 1 ) D G H m ( 3 m + 1 ) D G H 2 + ( 3 m 1 ) D H G u 3 D G u D G u u ] } u r m + 2 { 2 ( n 1 ) r 2 D ( m 1 ) D H + n 1 r [ ( m 1 ) D G u + 2 m G D 2 m ( m 1 ) D H G ] m [ ( 3 m + 1 ) G 2 D m ( 3 m 1 ) D H G 2 + ( 3 m 1 ) D H G r + ( 3 m 1 ) D G G u + 2 ( 3 m + 1 ) H F D + ( 3 m + 1 ) D F H + 3 ( m 1 ) D H F u ( 3 m 2 m + 2 ) D F H 2 3 F u D D F u u 2 D G r u ] + ( 3 m + 1 ) D G r + 3 ( m + 1 ) F D } u r m + 1 { n 1 r 2 D ( 2 m 1 ) G + 2 r + n 1 r [ m D G r m ( m 1 ) D G 2 2 m ( m 1 ) D H F + ( m 1 ) D F u + ( 2 m + 1 ) F D ] + ( 3 m + 1 ) D F r + m [ 2 D F r u + D G r r ( 3 m 1 ) D G G r ( 3 m 1 ) D F G r + m ( m 1 ) D G 3 ( 6 m + 1 ) D G F + 2 ( 3 m 2 3 m + 2 ) D F G H 3 ( m 1 ) D G F u 3 ( m 1 ) D H F r ] } u r m m { n 1 r D [ F r 2 ( m 1 ) F G 2 r F ] 3 ( m 1 ) D G F r ( 3 m 1 ) D F G u 3 ( m 1 ) D F F u + D F r r + ( 3 m 2 5 m + 4 ) D F 2 H 3 m F 2 D + ( m 1 ) ( 3 m 2 ) D G 2 F } u r m 1 m ( m 1 ) D F [ ( 3 m 4 ) F G 3 F r n 1 r F ] u r m 2 m ( m 1 ) ( m 2 ) D F 3 u r m 3 + ( H Q + H Q u + Q u u ) u r 2 + ( 2 Q r u + 2 H Q r + Q G u ) u r + Q r r + Q F u Q u F + G Q r ,
where the prime denotes the derivative with respect to u and the subscripts denotes the partial derivative with respect to the indicated variables. In view of the above expression, we need to distinguish several cases corresponding to different m for further discussion.

3. Equation (3) Admitting CLBS (1) and H–J SI (2)

It is noted from the expression (15) that for m = 2 , 1 , 2 , the terms in (15) may be combined. For instance, the terms containing u r m + 3 and u r can be put together when m = 2 . The over-determined system of the coefficient functions H ( u ) , G ( r , u ) , F ( r , u ) and Q ( r , u ) is different for different m. Each case of m = 2 , 1 , 2 will be respectively considered.
For the case of m = 2 , it follows from (15) that the coefficient functions in (3) and (1) satisfy the system
D 7 H D 6 H D + 16 H 2 D 2 D H 12 D H 3 + 14 D H H = 0 , n 1 r ( D 3 H D + 2 D H 2 D H ) 8 G D 6 D G u + 30 H D G + 14 D G H 28 D G H 2 + 10 D H G u 2 D G u u = 0 , 2 ( n 1 ) r 2 ( D H D ) n 1 r ( D G u + 4 G D 4 D H G ) + 14 G 2 D 7 D G r 9 F D 20 D H G 2 + 10 D H G r + 10 D G G u + 28 H F D + 14 D F H + 6 D H F u 24 D F H 2 6 F u D 2 D F u u 4 D G r u = 0 , n 1 r ( 2 D G r 2 D G 2 4 D H F + D F u + 5 F D ) + 3 ( n 1 ) r 2 D G + 2 ( n 1 ) r 3 D 7 D F r 4 D F r u 2 D G r r + 10 D G G r 4 D G 3 + 6 D H F r + 10 D F G u + 26 D G F 32 D F G H + 6 D G F u + H Q + H Q u + Q u u = 0 , 4 ( n 1 ) r 2 D F + 2 ( n 1 ) r D ( 2 F G F r ) + 6 D G F r + 10 D F G r + 6 D F F u 2 D F r r 12 D F 2 H + 12 F 2 D 8 D G 2 F + 2 Q r u + 2 H Q r + Q G u = 0 , 2 ( n 1 ) r D F 2 + 6 D F F r 4 D G F 2 + Q r r + Q F u Q u F + G Q r = 0
with D = u k .
Setting
A ( u ) = D 2 D H , B ( r , u ) = n 1 r D 2 D G , C ( r , u ) = 2 D F , E ( r , u ) = Q ,
system (16) reads as
5 D D A 7 D A A 2 D 2 A + 2 A D D 5 A D 2 + 8 A 2 D 3 A 3 = 0 , 6 ( n 1 ) r D ( A D D A A 2 ) + 7 D B A + 5 D A B u 11 A B D + 2 B D 2 3 D D B u B D D + 2 D 2 B u u + 7 B A 2 = 0 , 3 n ( n 1 ) r 2 A D 2 2 ( n 1 ) r D ( 2 D B u B D + 4 A B ) 3 B 2 D 8 A C D + 5 D B B u + 3 A D C u + 7 C D A u + 5 A D B r + 2 D 2 C u u + 4 D 2 B r u D D C u 2 D D B r + 6 C A 2 + 5 A B 2 = 0 , n 1 r 2 D 2 [ ( n + 1 ) B + 5 C ] n 1 r D ( 3 D B r + 2 D C u + 6 C A + 2 B 2 ) + 5 C D B r + 3 B D C u A D Q u + 2 B C D + A Q D + 5 D B B r + 3 A D C r D Q A + D D Q u + 4 D 2 C r u Q D 2 + 2 D 2 Q u u + 2 D 2 B r r + D Q D + 8 A B C + B 3 = 0 , n 1 r 2 D 2 C n 1 r D ( D C r + 2 B C ) + 5 C D B r 2 A D Q r D Q B u + Q B D + 3 B D C r + 3 D C C u + 2 D D Q r + 2 D 2 C r r + 4 D 2 Q r u + 3 A C 2 + 2 C B 2 = 0 , n 1 r D 2 Q r + 3 C C r + D C Q u D B Q r D Q C u + C Q D + 2 D 2 Q r r + C 2 B = 0 .
Formula (2) is a H–J SI of Equation (3) for the case of m = 2 if the coefficient functions satisfy system (18). This result can be proved in the same manner as that for Proposition 2.1 in [43]. We omit the details here. It follows that the CLBS (1) is equivalent to the H–J SI (2) for Equation (3). In other words, the CLBS method provides a symmetry interpretation to the method of H–J SI.
Substituting D = u k into the first one of the system (16), we derive that H satisfies the second-order ordinary differential equation (ODE)
2 u 3 H + 2 ( 3 k 7 u H ) u 2 H + 12 u 3 H 3 16 k u 2 H 2 + 7 k ( k 1 ) u H k ( k 1 ) ( k 2 ) = 0 .
It seems impossible to obtain the general solution of this nonlinear differential equation, and so we restrict ourselves to considering the case of u H = h , where h is an arbitrary constant. Then, (19) reads as
( 2 h k ) ( 2 h k + 1 ) ( 3 h k + 2 ) = 0 ,
which has three solutions
h = k 2 , k 1 2 , k 2 3 .
Thus, Equation (19) has three particular solutions listed as
H ( u ) = k 2 u , k 1 2 u , k 2 3 u .
Substituting D ( u ) = u k and H ( u ) = k 2 u into the second one of (16), we deduce that G ( r , u ) satisfies
2 G u u + k u G u k u 2 G + k ( n 1 ) 2 r u 2 = 0 .
The corresponding solutions of this linear differential equation is easily derived, which is listed as follows:
( i )
For k 2 , 0 ,
G ( r , u ) = g 1 ( r ) u k 2 + g 2 ( r ) u + n 1 2 r .
( i i )
For k = 2 ,
G ( r , u ) = g 1 ( r ) u + g 2 ( r ) u ln u + n 1 2 r .
( i i i )
For k = 0 ,
G ( r , u ) = g 1 ( r ) u + g 2 ( r ) .
For case (i), the third one of (16) is simplified as
2 F u u + 3 k u F u + k ( k 2 ) u 2 F + k u k 1 g 1 2 ( r ) 2 ( 2 k + 5 ) u g 2 2 ( r ) ( 3 k + 10 ) u k 2 g 1 ( r ) g 2 ( r ) + 2 ( k + 2 ) g 2 ( r ) 2 ( k + 2 ) ( n 1 ) r g 2 ( r ) = 0 .
The form of F ( r , u ) can be determined by solving this linear differential equation, which includes the undetermined coefficients f 1 ( r ) and f 2 ( r ) . Similarly, the fourth one of (16) will determine the form of Q ( r , u ) which includes the unknown coefficients q 1 ( r ) and q 2 ( r ) . The last two ones of the system (16) yield that the undetermined parts f 1 ( r ) , f 2 ( r ) , g 1 ( r ) , g 2 ( r ) , q 1 ( r ) and q 2 ( r ) satisfy a system of nonlinear coupled ODEs. It is still very difficult to manipulate this nonlinear coupled system. However, we can solve the system explicitly for several special situations. The computations involved here are quite lengthy and it is impossible to solve the nonlinear system by hand. The software Maple is used to help us complete the tedious computations. The detailed procedure is omitted and we just present the corresponding results in Table 1. Similarly, we can identify the admissible functions D ( u ) , Q ( r , u ) and the corresponding CLBS (1) of Equation (3) for other cases. The obtained results are also listed in Table 1 and so does for the results for m = 1 and m = 2 . It is noteworthy that some solutions of (16) lead to the vanishing of all the coefficients of (15). Since the vanishing of all the coefficients of (15) yields the over-determined system for general m, which will be discussed in Section 4, we omit these results here. Thus, we only consider the results that ensure that there exists at least one non-zero coefficient in the polynomial (15). For case 1 of Table 1, it is noted that we can set a = 0 due to the transformation u u + ( 2 n + 3 ) a t .
For general m, it follows from Equation (15) that the coefficient functions H ( u ) , G ( r , u ) , F ( r , u ) and Q ( r , u ) satisfy the following over-determined system:
D ( 3 m + 1 ) H D 3 m H D + m ( 3 m + 2 ) H 2 D m D H m 2 ( m + 1 ) D H 3 + m ( 3 m + 1 ) D H H = 0 , n 1 r [ D ( 2 m 1 ) H D + m ( m 1 ) D H 2 ( m 1 ) D H ] ( 3 m + 2 ) G D + m [ 3 ( 2 m + 1 ) H D G m ( 3 m + 1 ) D G H 2 + ( 3 m + 1 ) D G H + ( 3 m 1 ) D H G u 3 D G u D G u u ] = 0 , n 1 r [ ( m 1 ) D G u + 2 m G D 2 m ( m 1 ) D H G ] + 2 ( n 1 ) r 2 [ D ( m 1 ) D H ] + ( 3 m + 1 ) D G r + 3 ( m + 1 ) F D m [ ( 3 m + 1 ) G 2 D + ( 3 m 1 ) D H G r + ( 3 m 1 ) D G G u m ( 3 m 1 ) D H G 2 + 2 ( 3 m + 1 ) H F D + ( 3 m + 1 ) D F H + 3 ( m 1 ) D H F u ( 3 m 2 m + 2 ) D F H 2 3 F u D D F u u 2 D G r u ] = 0 , n 1 r [ m D G r m ( m 1 ) D G 2 2 m ( m 1 ) D H F + ( m 1 ) D F u + ( 2 m + 1 ) F D ] n 1 r 2 D [ ( 2 m 1 ) G + 2 r ] + ( 3 m + 1 ) D F r + m [ 2 D F r u + D G r r ( 3 m 1 ) D G G r + m ( m 1 ) D G 3 ( 3 m 1 ) D F G u ( 6 m + 1 ) D G F + 2 ( 3 m 2 3 m + 2 ) D F G H 3 ( m 1 ) D G F u 3 ( m 1 ) D H F r ] = 0 , n 1 r D [ F r 2 ( m 1 ) F G 2 r F ] 3 ( m 1 ) D G F r ( 3 m 1 ) D F G r 3 ( m 1 ) D F F u + D F r r + ( 3 m 2 5 m + 4 ) D F 2 H 3 m F 2 D + ( m 1 ) ( 3 m 2 ) D G 2 F = 0 , D F [ ( 3 m 4 ) F G 3 F r n 1 r F ] = 0 , H Q + H Q u + Q u u = 0 , 2 Q r u + 2 H Q r + Q G u = 0 , Q r r + Q F u Q u F + G Q r = 0 , D F 3 = 0 .
It is noted that hereafter a , b , c and d in the above table are arbitrary constants.
The last one of (20) implies F = 0 . It is equivalent to the system
( 2 m + 1 ) D D A ( 3 m + 1 ) D A A m D 2 A + m A D D ( 2 m + 1 ) A D 2 + ( 3 m + 2 ) A 2 D ( m + 1 ) A 3 = 0 , 2 ( m + 1 ) ( n 1 ) r D ( D A A D + A 2 ) ( 3 m + 1 ) D B A + ( 6 m 1 ) A B D ( 3 m 1 ) A D B u 2 ( m 1 ) B D 2 + ( 2 m 1 ) D D B u m D 2 B u u + ( m 1 ) B D D ( 3 m + 1 ) A 2 B = 0 , ( m + 1 ) n ( n 1 ) r 2 A D 2 2 ( n 1 ) r D [ ( m 1 ) B D m D B u 2 m A B ] + 3 ( m 1 ) B 2 D ( 3 m 1 ) D B B u ( 3 m 1 ) A D B r + 2 ( m 1 ) D D B r 2 m D 2 B r u ( 3 m 1 ) A B 2 = 0 , ( n 1 ) ( m n n + 1 ) r 2 D 2 B n 1 r D [ 2 ( m 1 ) B 2 + ( 2 m 1 ) D B r ] + ( m 1 ) B 3 + ( 3 m 1 ) D B B r + m D 2 B r r = 0 , D Q A + A Q D A D Q u + D Q D Q D 2 + D D Q u + m D 2 Q u u = 0 , B Q D D Q B u + 2 D D Q r 2 A D Q r + 2 m D 2 Q r u = 0 , n 1 r D 2 Q r + D C Q u B D Q r D Q C u + C Q D + m D 2 Q r r = 0 ,
where A ( u ) = D m D H , B ( r , u ) = n 1 r m D G . It turns out that Equation (3) admits the H–J SI (2) with C ( r , u ) = 0 and E ( r , u ) = Q .
On substitution of D = u k into the first one of the system (20), we deduce that H satisfies the second-order ODE
m u 3 H + m [ 3 k u 2 ( 3 m + 1 ) u 3 H ] H + m 2 ( m + 1 ) ( u H ) 3 m ( 3 m + 2 ) k ( u H ) 2 + ( 3 m + 1 ) k ( k 1 ) u H k ( k 1 ) ( k 2 ) = 0 .
It seems impossible to obtain the general solution of Equation (22). However, it is ready to write down three special solutions, which are listed as
H ( u ) = k m u , H ( u ) = k 1 m u , H ( u ) = k 2 ( m + 1 ) u .
Substituting D ( u ) , H ( u ) and F ( r , u ) into (20), we can obtain the explicit forms of G ( r , u ) and Q ( r , u ) . The results are summarized in Table 2.
Due to the relation (17), we can derive the H–J SI (2) of Equation (3) from the results listed in Table 1 and Table 2.

4. Symmetry Reductions of Equation (3)

Symmetry reductions of Equation (3) can be derived by using the compatibility of the invariant surface condition η = 0 for the admitted CLBS (1) and the governing Equation (3). One first solves the ODE η = 0 to obtain u as a function of r with r-independent integration constants, and then substitutes this solution into Equation (3) to determine the time evolution of these constants. We just present several examples to illustrate the reduction procedure.
Example 1.
Equation
u t = 1 r 2 r 2 1 u 2 u r 2 r + a u
admits the second-order CLBS
η = u r r 3 2 u u r 2 2 r u r + ( b 18 ) u r 2
and the first-order H–J SI
J = u t 1 u 3 u r 3 6 r u 2 u r 2 + 2 ( b + 8 ) r 2 u u r a u .
The compatibility of η = 0 and the governing equation yields the following results:
  • For b > 27 2 ,
    u ( r , t ) = 4 ( 2 b 27 ) α ( t ) sinh d 1 ln r + β ( t ) cosh d 1 ln r 2 r 3 ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = 3 8 α 3 + 1 4 d 1 α 2 β 3 8 α β 2 1 4 d 1 β 3 1 2 a α , β = 3 8 β 3 1 4 d 1 α β 2 + 3 8 α 2 β + 1 4 d 1 α 3 1 2 a β .
  • For b = 27 2 ,
    u ( r , t ) = 4 α ( t ) + β ( t ) ln r 2 r 3 ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = 3 2 α β 2 1 2 a α + β 3 , β = 3 2 β 3 1 2 a β .
    The solutions of this system of ODEs with a = 0 are listed below:
    α ( t ) = c 2 ln ( c 1 3 t ) 3 c 1 3 t , β ( t ) = ± 1 c 1 3 t .
  • For b < 27 2 ,
    u ( r , t ) = 4 ( 2 b 27 ) α ( t ) sin d 2 ln r + β ( t ) cos d 2 ln r 2 r 3 ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = 3 8 α 3 + 1 4 d 2 α 2 β 3 8 α β 2 + 1 4 d 2 β 3 1 2 a α , β = 3 8 β 3 1 4 d 2 α β 2 3 8 α 2 β 1 4 d 2 α 3 1 2 a β .
    It is noted that d 1 = 2 b 27 2 and d 2 = 27 2 b 2 .
Example 2.
Equation
u t = 1 r n 1 ( r n 1 u 2 u r 1 ) r + ( a + b u a n ) u ln r + ( c + d u a n ) u
admits the second-order CLBS
η = u r r n + a n u u r 2 + 1 r u r
and the first-order H–J SI
J = u t n u 2 r u r 1 a + b u a n u ln r + c + n a n + d u a n u .
The compatibility of η = 0 and the governing equation yields the following results:
  • For a 0 ,
    u ( r , t ) = a n α ( t ) ln r + β ( t ) n a ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = a ( c + 1 ) n a n 2 α + 1 n a 2 β + b , β = a ( c + 1 ) n a n 2 β + 1 n a 2 α 1 β 2 + d .
    The solutions blow up along the curves r = exp β ( t ) / α ( t ) when a > 0 and extinguish along the curves when a < 0 .
  • For a = 0 ,
    u ( r , t ) = r α ( t ) β ( t ) ,
    where α ( t ) and β ( t ) are given as
    (i) 
    For b 0 and b n ,
    α ( t ) = b t + c 1 , β ( t ) = c 2 ( b t + c 1 ) n b c + d + 1 b n ( b t + c 1 ) .
    (ii) 
    For b = 0 ,
    α ( t ) = c 1 , β ( t ) = c 2 exp n c 1 t ( c + d + 1 ) c 1 n .
    (iii) 
    For b = n ,
    α ( t ) = n t + c 1 , β ( t ) = c 2 + c + d + 1 n ln ( n t + c 1 ) ( n t + c 1 ) .
Example 3.
Equation
u t = 1 r n 1 ( r n 1 u k u r 2 ) r + a u + [ b r 2 ( k 3 ) ( 2 k n + k 6 n 4 ) r + c ] u k 2
admits the second-order CLBS
η = u r r + 2 k u u r 2 + 1 2 r u r
and the first-order H–J SI
J = u t ( 4 k ) u k 1 u r 1 a u n r u k u r 2 [ b r 2 ( k 3 ) ( 2 k n + k 6 n 4 ) r + c ] u k 2 .
The compatibility of η = 0 and the governing equation yields the following results.
  • For k 3 ,
    u ( r , t ) = ( 3 k ) 2 α ( t ) r + β ( t ) 1 3 k ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = 1 2 ( k 3 ) ( 4 n k + k 12 n 4 ) α 1 β + ( 3 k ) a α + b 2 , β = n ( k 3 ) 2 α 2 β 2 + ( 3 k ) a β + c .
    The solutions blow up along the curves r = β ( t ) / 2 α ( t ) 2 when k > 3 and extinguish along the curves when k < 3 .
  • For k = 3 ,
    u ( r , t ) = exp [ α ( t ) r + β ( t ) ] ,
    where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
    α = 2 α 1 + b , β = 4 n α 2 + a + c .
Example 4.
Equation
u t = 1 r n 1 ( r n 1 u k u r m ) r + a u + b + c r m + 1 m u 1 k m , k m + 1
admits the second-order CLBS
η = u r r + k 1 m u u r 2 1 m r u r
and the first-order H–J SI
J = u t u k 1 u r m + 1 n r u k u r m a u b + c r m + 1 m u 1 k m .
The compatibility of η = 0 and the governing equation gives
u ( r , t ) = ( k + m 1 ) α ( t ) m + 1 r m + 1 m + β ( t ) m m m + k 1 ,
where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
α = ( m + k 1 ) n + m + 1 m α m + 1 + ( m + k 1 ) a m α + ( m + 1 ) c m , β = ( m + k 1 ) n m α m β + ( m + k 1 ) a m β + b .
Notably, exact solutions of this form were found by Cherniha et al. [32] in the special case n = 2 and k = a = b = c = 0 , which corresponds to the turbulent gas flows. The solutions generalize the instantaneous source solution of the porous medium equation. It is noted that the solutions blow up along the curves r = ( m + 1 ) β ( t ) / m α ( t ) + m m + 1 when m ( m + k 1 ) < 0 and extinguish along the curves when m ( m + k 1 ) > 0 , namely, the interfaces of the solutions are the curves r = ( m + 1 ) β ( t ) / m α ( t ) + m m + 1 .
Example 5.
Equation
u t = 1 r n 1 ( r n 1 u m + 1 u r m ) r + a ln u + b + c r m + 1 m u
admits the second-order CLBS
η = u r r 1 u u r 2 1 m r u r
and the first-order H–J SI
J = u t u m u r m + 1 n r u m + 1 u r m a ln u + b + c r m + 1 m u .
The compatibility of η = 0 and the governing equation gives
u ( r , t ) = exp m m + 1 α ( t ) r m + 1 m + β ( t ) ,
where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
α = m + 1 m α m + 1 + a α + ( m + 1 ) c m , β = n α m β + a β + b .
The solutions are also generalizations of the instantaneous source solution of the porous medium equation.
Example 6.
Equation
u t = 1 r n 1 r n 1 u ( m 1 ) 2 ( n + 2 ) + 5 m 3 m n + m n + 1 u r m r + b + c r m + 1 m 1 u m n + 2 m n + 2 m n + m n + 1 + a u
admits the second-order CLBS
η = u r r m n + 2 m n + 2 ( m n + m n + 1 ) u u r 2 2 ( m 1 ) r u r
and the first-order H–J SI
J = u t u m ( m n + 2 m n + 2 ) m n + m n + 1 u r m + 1 m n + m n + 1 ( m 1 ) r u ( m 1 ) 2 ( n + 2 ) + 5 m 3 m n + m n + 1 u r m b + c r m + 1 m 1 u m n + 2 m n + 2 m n + m n + 1 a u .
The compatibility of η = 0 and the governing equation gives
u ( r , t ) = m 1 m n + m n + 1 α ( t ) r m + 1 m 1 m + 1 m n + m n + 1 β ( t ) m n + m n + 1 m + 1 ,
where α ( t ) and β ( t ) satisfy the two-dimensional dynamical system
α = ( m + 1 ) 2 ( m 1 ) 2 α m β ( m + 1 ) a m n + m n + 1 α + ( m + 1 ) c m 1 , β = ( m + 1 ) a m n + m n + 1 β + b .
The solutions blow up along the curves r = ( m + 1 ) β ( t ) / ( m 1 ) α ( t ) + m 1 m + 1 when ( m + 1 ) ( m n + m n + 1 ) > 0 and extinguish along the curves when ( m + 1 ) ( m n + m n + 1 ) < 0 .

5. Conclusions

We have described the structure of the second-order CLBSs (1) and the first-order H–J SIs (2) for second-order nonlinear reaction–diffusion Equations (3) with gradient-dependent diffusivity. Equation (3) admitting CLBSs (1) and H–J SIs (2) are identified. The corresponding reductions of the resulting equations are established, which in general reduce Equation (3) to two-dimensional dynamical systems due to the compatibility of the invariant surface condition (11) and the governing Equation (3). Those solutions extend the known ones such as instantaneous source solutions of the porous medium equation with absorption term. The phenomena of blow up and extinguish of the obtained solutions are also described. In fact, the invariant surface condition (11) is just the second-order DC of Equation (3). For Equation (3), the first-order H–J SI (2) is a natural extension of the conditional symmetry with the characteristic (14).
For the second-order nonlinear diffusion equation, we found that the second-order CLBS (1) is very effective and may yield some interesting results. Are there other choices of CLBS? Higher-order CLBS may lead to some new results. The discussion about the study of second-order CLBS and first-order H–J SI can be also extended to consider multi-dimensional nonlinear diffusion equations or other types of evolution equations. The CLBSs related to invariant subspace [14,15,16], sign-invariant [23,28] and separation of variables [24,25] for scalar diffusion equation can be generalized to study classifications and reductions of diffusion system. All of these problems will be involved in our future research.

Author Contributions

All the authors contributed equally to the work. All authors read and approved the final manuscript.

Funding

This work was supported by the Chinese National Natural Science Foundation (No. 11501175 and No. 11401529), the Subsidy scheme for young backbone teachers of higher education institutions in He’nan Province (No. 2015GGJS-080) and the Natural Science Foundation of Zhejiang Province (No. LY18A010033).

Acknowledgments

The authors are grateful to the unknown reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lie, S. Über die integration durch bestimmte integrale von einer klasse linearer partieller differentialgleichungen. Arch. Math. 1881, 6, 328–368. [Google Scholar]
  2. Noether, E. Invariante Variationsprobleme. Nachr. König. Gesell. Wissen. Göttingen, Math. Phys. Kl. 1918, 1918, 235–257. [Google Scholar]
  3. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  4. Fushchych, W.I.; Shtelen, W.M.; Serov, M.I. Symmetry Analysis and Exact Solutions of Equations of Nonlinear Mathemtical Physics; Kluwer: Dordrecht, The Netherlands, 1993. [Google Scholar]
  5. Cherniha, R.; Serov, M.; Pliukhin, O. Nonlinear Reaction–Diffusion-Convection Equations: Lie and Conditional Symmetry, Exact Solutions and Their Applications; CRC Press: New York, NY, USA, 2018. [Google Scholar]
  6. Zhdanov, R.Z. Conditional Lie–Bäcklund symmetry and reduction of evolution equation. J. Phys. A: Math. Gen. 1995, 28, 3841–3850. [Google Scholar] [CrossRef]
  7. Fokas, A.S.; Liu, Q.M. Nonlinear interaction of traveling waves of nonintegrable equations. Phys. Rev. Lett. 1994, 72, 3293–3296. [Google Scholar] [CrossRef] [PubMed]
  8. Basarab-Horwath, P.; Zhdanov, R.Z. Initial-value problems for evolutionary partial differential equations and higher-order conditional symmetries. J. Math. Phys. 2000, 42, 376–389. [Google Scholar] [CrossRef]
  9. Zhdanov, R.Z.; Andreitsev, A.Y. Non-classical reductions of initial-value problems for a class of nonlinear evolution equations. J. Phys. A: Math. Gen. 2000, 33, 5763–5781. [Google Scholar] [CrossRef] [Green Version]
  10. Zhdanov, R.Z. Higher conditional symmetry and reductions of initial-value problems. Nonlinear Dyn. 2004, 28, 17–27. [Google Scholar] [CrossRef]
  11. Sergyeyev, A. Constructing conditionally integrable evolution systems in (1+1)-dimensions: A generalization of invariant modules approach. J. Phys. A: Math. Gen. 2002, 35, 7563–7660. [Google Scholar] [CrossRef]
  12. Ji, L.N.; Qu, C.Z. Conditional Lie–Bäcklund symmetry of evolution system and application for reaction–diffusion system. Stud. Appl. Math. 2014, 133, 118–149. [Google Scholar] [CrossRef]
  13. Wang, J.P.; Ji, L.N. Conditional Lie–Bäcklund symmetry, second-order differential constraint and direct reduction of diffusion systems. J. Math. Anal. Appl. 2015, 427, 1101–1118. [Google Scholar] [CrossRef]
  14. Ji, L.N.; Qu, C.Z. Conditional Lie–Bäcklund symmetries and invariant subspaces to nonlinear diffusion equations. IMA J. Appl. Math. 2011, 76, 610–632. [Google Scholar] [CrossRef]
  15. Ji, L.N.; Qu, C.Z. Conditional Lie–Bäcklund symmetries and invariant subspaces to nonlinear diffusion equaitons with convection and source. Stud. Appl. Math. 2013, 131, 266–301. [Google Scholar]
  16. Qu, C.Z.; Ji, L.N. Invariant subspaces and Conditional Lie–Bäcklund symmetries of inhomogeneous nonlinear diffusion equations. Sci. China Math. 2013, 56, 2187–2203. [Google Scholar] [CrossRef]
  17. Svirshchevskii, S.R. Invariant linear spaces and exact solutions of nonlinear evolution equations. J. Nonlinear Math. Phys. 1996, 3, 164–169. [Google Scholar] [CrossRef]
  18. Galaktionov, V.A.; Svirshchevskii, S.R. Exact solutions and Invariant Subspaces of Nonlinear Partial Differential Equations in Mechanics and Physics; Chapman and Hall: London, UK, 2007. [Google Scholar]
  19. Qu, C.Z. Exact solutions to nonlinear diffusion equations obtained by generalized conditional symmetry. IMA J. Appl. Math. 1999, 62, 283–302. [Google Scholar]
  20. Qu, C.Z. Group classification and generalized conditional symmetry reduction of the nonlinear diffusion-convection equation with a nonlinear source. Stud. Appl. Math. 1997, 99, 107–136. [Google Scholar] [CrossRef]
  21. Ji, L.N.; Qu, C.Z. Conditional Lie–Bäcklund symmetries and solutions to (n+1)-dimensional nonlinear diffusion equations. J. Math. Phys. 2007, 484, 103509. [Google Scholar] [CrossRef]
  22. Ji, L.N.; Qu, C.Z.; Ye, Y.J. Solutions and symmetry reductions of the n-dimensional nonlinear convection-diffusion equations. IMA J. Appl. Math. 2010, 75, 17–55. [Google Scholar] [CrossRef]
  23. Qu, C.Z.; Ji, L.N.; Wang, L.Z. Conditional Lie–Bäcklund symmetries and sign-invarints to quasi-linear diffusion equations. Stud. Appl. Math. 2007, 119, 355–391. [Google Scholar] [CrossRef]
  24. Qu, C.Z.; Zhang, S.L.; Liu, R.C. Separation of variables and exaxt solutions to quasilinear diffusion equations with nonlinear source. Phys. D: Nonlinear Phenom. 2000, 144, 97–123. [Google Scholar] [CrossRef]
  25. Zhang, S.L.; Lou, S.Y.; Qu, C.Z. New variable separation approach: Application to nonlinear diffusion equations. J. Phys. A: Math. Gen. 2003, 36, 12223–12242. [Google Scholar] [CrossRef]
  26. Yanenko, N.N. Theory of consistency and methods of integrating systems of nonlinear partial differential equations. In Proceedings of the Fourth All-Union Mathematical Congress, Leningrad, Russia; 1964; pp. 247–259. (In Russian). [Google Scholar]
  27. Sidorov, A.F.; Shapeev, V.P.; Yanenko, N.N. Method of Differential Constraints and Its Application to Gas Dynamics; Nauka: Novosibirsk, Russia, 1984. [Google Scholar]
  28. Qu, C.Z.; Estévez, P.G. On nonlinear diffusion equations with x-dependent convection and absorption. Nonlinear Anal. Theory Methods Appl. 2004, 37, 549–557. [Google Scholar] [CrossRef]
  29. Pascal, J.P.; Pascal, H. On some diffusive waves in nonlinear heat conduction. Int. J. Nonlinear Mech. 1993, 28, 641–649. [Google Scholar] [CrossRef]
  30. Pascal, J.P.; Pascal, H. On some nonlinear shear flows of non-Newtonian fluids. Int. J. Nonlinear Mech. 1995, 30, 487–500. [Google Scholar] [CrossRef]
  31. Assis, P.C.; Silva, L.R. Nonlinear diffusion equations, Tsallis formalism and exact solutions. J. Math. Phys. 2005, 46, 1–7. [Google Scholar] [CrossRef]
  32. Cherniha, R.; King, J.R.; Kovalenko, S. Lie symmetry properties of nonlinear reaction–diffusion equations with gradient-dependent diffusivity. Commun. Nonlinear Sci. Numer. Simul. 2016, 36, 98–108. [Google Scholar] [CrossRef]
  33. Barenblatt, G.I. On self-similar motions of compressible fluid in a porous medium. Appl. Math. Mech. 1952, 16, 679–698. [Google Scholar]
  34. Akhatov, I.S.; Gazizov, R.K.; Ibragimov, N.H. Group classification of the equations of nonlinear filtration. Sov. Math. Dokl. 1987, 35, 384–386. [Google Scholar]
  35. Fridman, A. Partial Differential Equations of Parabolic Type; PrenticeHall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  36. Galaktionov, V.A. Geometric Sturmian Theory of Nonlinear Parabolic Equations and Applications; Chapman & Hall/CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  37. Samarskii, A.A.; Galaktionov, V.A.; Kurdyumov, S.P.; Mikhailov, A.P. Blow up in Quasilinear Parabolic Equations; Nauka: Moscow, Russia, 1987. [Google Scholar]
  38. Galaktionov, V.A.; Vázquez, J.L. A Stability Technique for Evolution Partial Differential Equations: A Dynamical Systems Approach; Birkhäuser: Boston, MA, USA, 2004. [Google Scholar]
  39. Galaktionov, V.A.; Dorodnitsyn, V.A.; Elenin, G.G.; Kurdyumov, S.P.; Samarskii, A.A. A quasilinear heat equation with a source: Peaking, localization, symmetry, exact solutions, asymptotics, structures. J. Sov. Math. 1988, 41, 1222–1292. [Google Scholar] [CrossRef]
  40. Goard, J.M.; Broadbridge, P. Nonclassical symmetry analysis of nonlinear reaction–diffusion equations in two spatial dimensions. Nonlinear Anal. Theory Methods Appl. 1996, 26, 735–754. [Google Scholar] [CrossRef]
  41. Broadbridge, P.; Bradshaw-Hajek, B.H.; Triadis, D. Exact nonclassical symmetry solutions of Arrhenius reaction–diffusion. Proc. R. Soc. A 2015, 471, 20150580. [Google Scholar] [CrossRef]
  42. Galaktionov, V.A. Quasilinear heat equations with first-order sign-invariants and new explicit solutions. Nonlinear Anal. Theory Methods Appl. 1994, 23, 1595–1621. [Google Scholar] [CrossRef]
  43. Galaktionov, V.A.; Posashkov, S.A. New explicit solutions of quasilinear heat equations with general first-order sign-invariants. Phys. D: Nonl. Phenom. 1996, 99, 217–236. [Google Scholar] [CrossRef]
  44. Galaktionov, V.A.; Posashkov, S.A. Maximal sign-invariants of quasilinear heat equations with gradient diffusivity. J. Math. Phys. 1998, 39, 4948–4964. [Google Scholar] [CrossRef]
  45. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar]
  46. Olver, P.J. Applications of Lie Groups to Differential Equations, 2nd ed.; Springer: New York, NY, USA, 1993. [Google Scholar]
  47. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  48. Andreev, V.K.; Kaptsov, O.V.; Pukhnachev, V.V.; Rodionov, A.A. Applications of Group-Theoretic Methods in Hydrodynamics; Kluwer: Dordrecht, The Netherland, 1998. [Google Scholar]
  49. Olver, P.J. Direct reduction and differential constraints. Proc. R. Soc. A 1994, 444, 509–523. [Google Scholar] [CrossRef]
Table 1. Admissible functions D ( u ) , Q ( r , u ) and CLBS (1) of Equation (3) for m = 2 , 1 , 2 .
Table 1. Admissible functions D ( u ) , Q ( r , u ) and CLBS (1) of Equation (3) for m = 2 , 1 , 2 .
No.mn D ( u ) Q ( r , u ) H ( u ) G ( r , u ) F ( r , u )
12 1 2 ( n + 3 ) a 1 2 u 1 2 r a r u
223 1 u 3 2 9 b c r 3 + 2 3 a r 3 2 + d u 2 + b u 6 c 2 u 1 2 r c r u 2
323 1 u 2 a u 3 2 u 2 r ( b 18 ) u r 2
42 1 u a r 3 2 ( n 6 ) ( 2 n 3 ) 2 27 r 3 + b u 1 u 1 2 r ( 3 2 n ) u 2 r 2
526 1 u 0 1 u 1 2 r 9 u 2 a 2 r 2
6 1 u k a u + ( b + c ln r ) u k 2 k 2 u k 1 k 2 u 1 r 0
7 1 u 2 ( a + b u a n ) u ln r + ( c + d u a n ) u n + a n u 1 r 0
8 1 2 u 2 2 ( h + 1 ) ln r + a u + ( b + c ln r ) u h h u 1 r 0
9 2 u k [ b r 2 ( k 3 ) ( 2 k n + k 6 n 4 ) r + c ] u k 2 + a u 2 k u 1 2 r 0
10 2 u 3 n r ( 2 n + 1 ) 2 + b u + ( b + c r ) u 4 n + 3 2 ( 2 n + 1 ) 4 n + 3 2 ( 2 n + 1 ) u 1 2 r 0
11 2 u 2 ( 9 n + 5 ) 3 ( 2 n + 1 ) b r 1 3 + ( 3 n + 1 ) r 3 ( 2 n + 1 ) 2 + c u 2 ( 3 n + 2 ) 3 ( 2 n + 1 ) + a u 2 ( 3 n + 2 ) 3 ( 2 n + 1 ) u 2 3 r 0
Table 2. Admissible Functions D ( u ) , Q ( r , u ) and CLBS (1) of Equation (3) for general m with F ( r , u ) = 0 .
Table 2. Admissible Functions D ( u ) , Q ( r , u ) and CLBS (1) of Equation (3) for general m with F ( r , u ) = 0 .
No. D ( u ) Q ( r , u ) H ( u ) G ( r , u )
1 u k a u + b + c r m + 1 m u 1 k m k 1 m u 1 m r
k m + 1
2 u m + 1 a ln u + b + c r m + 1 m u 1 u 1 m r
3 u ( m 1 ) 2 ( n + 2 ) + 5 m 3 m n + m n + 1 b + c r m + 1 m 1 u m n + 2 m n + 2 m n + m n + 1 + a u m n + 2 m n + 2 ( m n + m n + 1 ) u 2 ( m 1 ) r
4 u m 2 n m n + 2 m 2 + 2 m m n + m + 1 b + c r m + 1 m u m n + 2 m + 2 m n + m + 1 + a u m n + 2 m + 2 ( m n + m + 1 ) u 1 m r

Share and Cite

MDPI and ACS Style

Ji, L.; Feng, W. Second-Order Conditional Lie–Bäcklund Symmetries and Differential Constraints of Nonlinear Reaction–Diffusion Equations with Gradient-Dependent Diffusivity. Symmetry 2018, 10, 267. https://doi.org/10.3390/sym10070267

AMA Style

Ji L, Feng W. Second-Order Conditional Lie–Bäcklund Symmetries and Differential Constraints of Nonlinear Reaction–Diffusion Equations with Gradient-Dependent Diffusivity. Symmetry. 2018; 10(7):267. https://doi.org/10.3390/sym10070267

Chicago/Turabian Style

Ji, Lina, and Wei Feng. 2018. "Second-Order Conditional Lie–Bäcklund Symmetries and Differential Constraints of Nonlinear Reaction–Diffusion Equations with Gradient-Dependent Diffusivity" Symmetry 10, no. 7: 267. https://doi.org/10.3390/sym10070267

APA Style

Ji, L., & Feng, W. (2018). Second-Order Conditional Lie–Bäcklund Symmetries and Differential Constraints of Nonlinear Reaction–Diffusion Equations with Gradient-Dependent Diffusivity. Symmetry, 10(7), 267. https://doi.org/10.3390/sym10070267

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