Conditional Lie–Bäcklund Symmetries and Functionally Generalized Separation of Variables to Quasi-Linear Diffusion Equations with Source

Abstract

The conditional Lie–Bäcklund symmetry method is applied to investigate the functionally generalized separation of variables for quasi-linear diffusion equations with a source. The equations and the admitted conditional Lie–Bäcklund symmetries related to invariant subspaces are identified. The exact solutions possessing the form of the functionally generalized separation of variables are constructed for the resulting equations due to the corresponding symmetry reductions.

1. Introduction

The separation of variables is the most widely used method for solving linear heat- and mass-transfer equations and other linear equations of mathematical physics. The productive separable solutions were presented in [1,2]. The additive separable solutions were constructed for some classes of nonlinear first-order partial differential equations (PDEs) in [3,4], and for some nonlinear heat- and mass-transfer and wave equations of the second and higher orders in [5,6]. The functional separation of variables (FSV)

$$\begin{array}{c}\hfill f\left(u\right)=\varphi \left(x\right)+\psi \left(t\right)\end{array}$$

(1)

is sought for some nonlinear wave and heat equations. The solutions of this form were also obtained [7] for the nonlinear diffusion equation $u_t={\left[D\left(u\right)u_x\right]}_x$.

The nonlinear separation of variables (NSV)

$$\begin{array}{c}\hfill u=\varphi \left(x\right)\psi \left(t\right)+\eta \left(t\right)\end{array}$$

(2)

was introduced [8,9] for the study of the nonlinear parabolic equation. The generalized separation of variables (GSV)

$$\begin{array}{c}\hfill u={\tilde{C}}_1\left(t\right){\tilde{f}}_1\left(x\right)+{\tilde{C}}_2\left(t\right){\tilde{f}}_2\left(x\right)+\cdots +{\tilde{C}}_n\left(t\right){\tilde{f}}_n\left(x\right)(n\ge 2)\end{array}$$

(3)

can also provide exact solutions involving more summands, which was extensively discussed by the invariant subspace method [10].

The additive separable solutions and the productive separable solutions are both the particular case of the FSV (1) and the GSV (3). A further extension to the separable solutions is of the form

$$\begin{array}{c}\hfill f\left(u\right)=C_1\left(t\right)f_1\left(x\right)+C_2\left(t\right)f_2\left(x\right)+\cdots +C_n\left(t\right)f_n\left(x\right)(n\ge 2),\end{array}$$

(4)

which we call the functionally generalized separation of variables (FGSV). It can be regarded as a natural generalization of GSV in a similar way to how FSV is a generalization of the additive separation of variables. The more general separable solutions

$$\begin{array}{c}\hfill f\left(u\right)={\varphi }_0\left(t\right)+{\varphi }_1\left(t\right)g_1(x,t)+{\varphi }_2\left(t\right)g_2(x,t)+\cdots +{\varphi }_n\left(t\right)g_n(x,t)\end{array}$$

are proposed due to the additional generating condition method in [11]. All these different types of separable solutions are very useful to study different kinds of PDEs. This is particularly true for nonlinear diffusion equations [10,11,12,13].

Conditional Lie–Bäcklund symmetry (CLBS) [14]

$$\begin{array}{c}\hfill \rho ={\left[f\left(u\right)\right]}_{xt}\end{array}$$

(5)

is introduced to study FSV (1). The key point is due to the compatibility of $\rho =0$ and the governing equation. As a consequence, CLBS

$$\begin{array}{c}\hfill \sigma ={\left[f\left(u\right)\right]}_{nx}+a_1\left(x\right){\left[f\left(u\right)\right]}_{(n-1)x}+\cdots +a_n\left(x\right)f\left(u\right)\end{array}$$

(6)

can be used to study FGSV (4). It is noted that ${\left[f\left(u\right)\right]}_{kx}={\partial }^{k}f\left(u\right)/\partial x^k$ and $k=1,2,\cdots ,n$.

A great number of publications have been devoted to study the quasi-linear diffusion equation with a source

$$\begin{array}{c}\hfill u_t={\left[D\left(u\right)u_x\right]}_x+Q\left(u\right),x\in R,t>0,\end{array}$$

(7)

where $D\left(u\right)$ and $Q\left(u\right)$ are respectively the diffusion and source term. The involved methods include the Lie-point symmetry method [15], the nonclassical symmetry method [16,17], the CLBS method [18,19,20,21], the nonlocal symmetry method [22], the truncated Painlevé approach [23], the differential constraints method [24], the sign-invariant and invariant subspace method [25,26], the transformation method [27,28], the ansatz-based method [29], the spectral volume method [30], etc.

The studies about FSV (1), NSV (2), and FGSV (4) of Equation (7) can respectively refer to [14,31], [31,32], and [10,31]. CLBS (6) is the key point to consider FGSV (4) of Equation (7). Equation (7) admits CLBS (6) is equivalent to saying that the equation

$$\begin{array}{c}\hfill v_t=A\left(v\right)v_{xx}+B\left(v\right)v_x^2+C\left(v\right)\end{array}$$

(8)

admits the CLBS

$$\begin{array}{c}\hfill \eta =v_{nx}+a_1\left(x\right)v_{(n-1)x}+\cdots +a_n\left(x\right)v.\end{array}$$

(9)

In fact, Equation (8) can be derived from Equation (7) due to the transformation $v=f\left(u\right)$, and these two equations are related as

$$\begin{array}{c}\hfill \begin{array}{c}A\left(v\right)=D\left[g\left(v\right)\right],B\left(v\right)=A^{\prime }\left(v\right)+\frac{g^{″}\left(v\right)}{g^{\prime }\left(v\right)}A\left(v\right),C\left(v\right)=\frac{Q\left[g\right(v\left)\right]}{g^{\prime }\left(v\right)},\hfill \end{array}\end{array}$$

(10)

where $u=g\left(v\right)$ is the inverse function of $v=f\left(u\right)$.

CLBS (9) is the key point to give a symmetry interpretation to the invariant subspace method proposed by Galaktionov and Svirshchevskii [10]. Since Equation (8) admits CLBS (9), it is easy to know that the corresponding group invariant solutions of Equation (8) are exactly defined on the linear solution space

$$\begin{array}{c}\hfill S_n=S\left\{f_1\left(x\right),f_2\left(x\right),\cdots ,f_n\left(x\right)\right\}\equiv \left\{{\Sigma }_{i=1}^n{C}_i{f}_i\left(x\right),C_i\in R\right\}\end{array}$$

(11)

determined by the linear ordinary differential equation (ODE) $\eta =0$. It is proved that the maximal dimension of the linear solution space is five for the second-order diffusion Equation (8) [33]. This confines oneself to study CLBS (9) with $2\le n\le 5$ of Equation (8).

The classifications and reductions of the nonlinear diffusion equations with convection and the source

$$\begin{array}{c}\hfill u_t={\left[D\left(u\right)u_x\right]}_x+P\left(u\right)u_x+Q\left(u\right)\end{array}$$

due to CLBS (6) were studied in [34], where the studies were concerned with the case of $P\left(u\right)\ne 0$. Indeed, some results can degenerate to the case of $P\left(u\right)=0$. However, most of these degenerated cases correspond to x-independent $a_1\left(x\right)$ and x-independent $a_2\left(x\right)$. Here, we are mainly concerned with x-dependent $a_i\left(x\right)(i=1,2)$.

The remainder of this paper is arranged as follows. In Section 2, the classification of Equation (7) due to CLBS (6) is presented. In the subsequent section, the FGSV (4) of Equation (7) is constructed. The last section is devoted to the conclusions.

2. Conditional Lie-Bäcklund Symmetry (6) of Equation (7)

Consider a nonlinear evolution equation

$$\begin{array}{c}\hfill u_t=F(t,x,u,u_x,u_{2x},\cdots ,u_{kx})\end{array}$$

(12)

with Lie–Bäcklund vector field (LBVF)

$$\begin{array}{c}\hfill V=\eta \frac{\partial }{\partial u}+D_x\eta \frac{\partial }{\partial u_x}+D_t\eta \frac{\partial }{\partial u_t}+D_x^2\eta \frac{\partial }{\partial u_{2x}}+D_x{D}_t\eta \frac{\partial }{\partial u_{xt}}\cdots ,\end{array}$$

(13)

where $\eta =\eta (x,t,u,u_x,u_t,u_{2t},u_{xt},\cdots )$ is the characteristic of LBVF (13) and the total differentiation operators respectively denote

$$\begin{array}{c}{D}_t=\frac{\partial }{\partial t}+u_t\frac{\partial }{\partial u}+u_{2t}\frac{\partial }{\partial u_t}+u_{xt}\frac{\partial }{\partial u_x}+\cdots ,\\ D_x=\frac{\partial }{\partial x}+u_x\frac{\partial }{\partial u}+u_{xt}\frac{\partial }{\partial u_t}+u_{2x}\frac{\partial }{\partial u_x}+\cdots .\end{array}$$

Definition 1.

The evolutionary vector field (13) is said to be a Lie–Bäcklund symmetry of the evolution Equation (12) if the condition

$$\begin{array}{c}\hfill V\left(u_t-F\right){\mid }_M=0\end{array}$$

holds, where M denotes the set of all differential consequences of the Equation (12).

Definition 2.

The evolutionary vector field (13) is said to be a CLBS of (12) if the following condition

$$\begin{array}{c}\hfill V\left(u_t-F\right){|}_{M\cap L_x}=0\end{array}$$

holds, where $L_x$ denotes the set of all differential consequences of equation $\eta =0$ with respect to x.

The following proposition is useful for further study.

Proposition 1

([18,19]). Equation (8) admits the CLBS (9) if there exists a function $W(t,x,v,\eta )$ such that

$$\begin{array}{c}\hfill \frac{\partial \eta }{\partial t}=[F,\eta ]+W(t,x,v,\eta ),W(t,x,v,0)=0,\end{array}$$

where $[F,\eta ]=F^{\prime }\eta -{\eta }^{\prime }F$ and $F\left[v\right]=A\left(v\right)v_{xx}+B\left(v\right)v_x^2+C\left(v\right)$, the prime denotes the Gateaux derivative, and W is an analytic function of t, x, v, $v_x$, ⋯ and η, $D_x\eta$, $D_x^2\eta$, ⋯.

A direct conclusion of this proposition is that Equation (8) admits the CLBS with the characteristic (9) if

$$\begin{array}{c}\hfill D_t{\eta |}_{L\cap M}=0,\end{array}$$

(14)

where M is the set of all differential consequences of the equation, that is,

$$\begin{array}{c}\hfill v_t-F\left[v\right]=0,D_x^j{D}_t^k\left(v_t-F\left[v\right]\right)=0,j,k=0,1,2,\cdots \end{array}$$

and L denotes the set of all differential consequences of equation $\eta =0$ with respect to x, that is,

$$\begin{array}{c}\hfill D_x^j\eta =0,j=0,1,2,\cdots .\end{array}$$

(14) is equivalent to

$$\begin{array}{c}\hfill {\left(F\left[v\right]\right)}_{nx}+a_1\left(x\right){\left(F\left[v\right]\right)}_{(n-1)x}+\cdots +a_n\left(x\right)\left(F\left[v\right]\right){|}_M=0,\end{array}$$

(15)

which is the right sufficient condition that the linear space (11) is invariant with respect to the operator $F[·]$ [10].

Different cases will be considered respectively according to different cases of $n=2,3,4,5$. For $n=2$, a direct computation from (15) gives a polynomial about $v_x$. The vanishing of all the coefficients yields the following over-determined system about the unknown functions in Equation (8) and the CLBS (9), which is listed as

$$\begin{array}{c}{B}^{″}=0,\hfill \\ a_1(4B^{\prime }+A^{″})=0,\hfill \\ 2(a_1^2-a_1^{\prime }-a_2)A^{\prime }-a_2[v(A^{″}+5B^{\prime })+B]+C^{″}+2(a_1^2-a_1^{\prime })B=0,\hfill \\ (3a_1{a}_2-2a_2^{\prime })v{A}^{\prime }-(a_1^{″}-2a_1{a}_1^{\prime }+2a_2^{\prime })A+2(2a_1{a}_2-a_2^{\prime })vB=0,\hfill \\ a_2^2{v}^2{A}^{\prime }+(2a_1^{\prime }{a}_2-a_2^{″})vA+2a_2^2{v}^{2}B+a_2(C-v{C}^{\prime })=0.\hfill \end{array}$$

(16)

We can not give the general solutions of the determining system (16) since it is a coupled nonlinear system of PDEs. The workable way is to find as many particular solutions as possible.

The first equation of the determining system (16) gives $B\left(v\right)=b_{1}v+b_2$. As a consequence, the second one of (16) can be simplified as

$$\begin{array}{c}\hfill a_1(4b_1+A^{″})=0.\end{array}$$

Thus, we consider two cases for further study.

• Case 1. $a_1\left(x\right)=0$.

Substituting $B\left(v\right)=b_{1}v+b_2$ and $a_1\left(x\right)=0$ into the system (16) and solving the resulting system, we can finally determine $A\left(v\right),B\left(v\right),C\left(v\right)$, and $a_2\left(x\right)$. The corresponding results are listed as

$$\begin{array}{l}\left(i\right)A\left(v\right)arbitrary,B\left(v\right)=b_{1}v+b_2,C\left(v\right)=c_{1}v+c_2,a_2\left(x\right)=0;\\ \left(ii\right)A\left(v\right)arbitrary,B\left(v\right)=b_{1}v,C\left(v\right)=[s(b_1{v}^2+A\left(v\right))+c_1]v,a_2\left(x\right)=s(s\ne 0);\\ \left(iii\right)A\left(v\right)=-\frac{1}{2}{b}_{2}v,B\left(v\right)=b_2,C\left(v\right)=c_{1}v,a_2^{″}\left(x\right)+3a_2^2\left(x\right)=0(a_2\left(x\right)\ne s);\\ \left(iv\right)A\left(v\right)=\frac{k_1}{v},B\left(v\right)=0,C\left(v\right)=c_{1}v+c_2,k_1[a_2^{″}\left(x\right)+a_2^2\left(x\right)]-c_2{a}_2\left(x\right)=0(a_2\left(x\right)\ne s).\end{array}$$

In [10], it is shown that the IS(11) is also invariant with respect to the operator $F\left[v\right]+pv$ if (11) is invariant with the operator $F\left[v\right]$. Thus, $v_t=F\left[v\right]+pv$ admits linear CLBS (9) if Equation (8) admits the linear CLBS (9). Hence, we can omit the linear term about v in $C\left(v\right)$. For the case $\left(i\right)$, the second one of (10) becomes

$$\begin{array}{c}\hfill \frac{g^{″}\left(v\right)}{g^{{}^{\prime }}\left(v\right)}=\frac{b_{1}v+b_2-A^{\prime }\left(v\right)}{A\left(v\right)}.\end{array}$$

It is well known that the diffusion coefficient $D\left(u\right)$ is a power function or exponential function. In general, we consider the following five cases to identify the transformation $u=g\left(v\right)$, which are presented as

$$\begin{array}{c}\left(1\right)A\left(v\right)=v^2,B\left(v\right)=kv(k\ne 1);\hfill \\ \left(2\right)A\left(v\right)=v^2,B\left(v\right)=v;\hfill \\ \left(3\right)A\left(v\right)=v^k(k\ne 1),B\left(v\right)=0;\hfill \\ \left(4\right)A\left(v\right)=v,B\left(v\right)=0;\hfill \\ \left(5\right)A\left(v\right)=v,B\left(v\right)=k(k\ne 0).\hfill \end{array}$$

For the case of $A\left(v\right)=v^2,B\left(v\right)=kv(k\ne 1)$, the second one of (10) is simplified as

$$\begin{array}{c}\hfill \frac{g^{″}\left(v\right)}{g^{{}^{\prime }}\left(v\right)}=\frac{k-2}{v}.\end{array}$$

$g\left(v\right)=v^{k-1}$ can be obtained by solving this equation, which is of the canonical coordinate form without considering the translation and scaling transformations. Consequently, $D\left(u\right)=u^{\frac{2}{k-1}}$ and $Q\left(u\right)=(k-1)q{u}^{\frac{k-2}{k-1}}$ are given from the first one and third one of (10) due to the transformation $u=v^{k-1}$. A similar discussion for Cases 2–5 will identify the corresponding governing Equation (7) and the admitted CLBS (6). The procedure of solving the system (16) for the case of (ii) is almost the same as that for the case of (i). Since the obtained $a_1\left(x\right)$ and $a_2\left(x\right)$ for these two cases are both x-independent and the resulting CLBS (6) and the governing Equation (7) are all special cases of the ones in [34], we do not consider these special cases for further study. So is the case for all the similar ones appearing below.

For the case of (iii), Equation (8) is invariant with respect to scaling transformation about x. Therefore, we consider $A\left(v\right)=v,B\left(v\right)=-2$ without loss of generality. The transformation $u=1/v^2$ is obtained by solving the second one of (10). Consequently, $D\left(u\right)=u^{-\frac{1}{2}}$ and $Q\left(u\right)=0$ are given. $a_2\left(x\right)$ satisfy the nonlinear ODE

$$\begin{array}{c}\hfill a_2^{″}\left(x\right)+3a_2^2\left(x\right)=0.\end{array}$$

Now that the general solution of this ODE is in terms of the Weierstrass function, the solution of second-order linear ODE $v_{xx}+a_2\left(x\right)v=0$ is not analytic. Hence, the analytic solution of Equation (8) can not be constructed due to the corresponding symmetry reductions. However, the particular CLBS (9) will lead to symmetry reductions of the governing Equation (8). Thus, we can satisfy our study by finding CLBS (9) in particular cases. Considering the form of the equation about $a_2\left(x\right)$, we apply the power law function constraint. $a_2\left(x\right)=-\frac{2}{x^2}$ is determined, which will yield the analytic solution of Equation (8). The corresponding result is listed as Case 9 for $s=-2$ in

Table 1. Conditional Lie–Bäcklund symmetry (CLBS) (6) for Equation (7).

No.	Equation (7)	CLBS (6)
1	$u_t={\left(u^{-\frac{1}{2}}{u}_x\right)}_x$	$\sigma ={\left(u^{\frac{1}{2}}\right)}_{xx}-\frac{6}{x^2}{u}^{\frac{1}{2}}$
2	$u_t={\left(u^{\frac{1}{k}}{u}_x\right)}_x+kq{u}^{\frac{k-1}{k}}$	$\sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\frac{1}{x}{\left(u^{\frac{1}{k}}\right)}_x$
3	$u_t={\left[exp\left(u\right)u_x\right]}_x+qexp(-u)$	$\sigma ={\left[exp\left(u\right)\right]}_{xx}-\frac{1}{x}{\left[exp\left(u\right)\right]}_x$
4	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x-\frac{3}{2}q{u}^{\frac{5}{3}}$	$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}-\frac{1}{x}{\left(u^{-\frac{2}{3}}\right)}_x$
5	$u_t={\left(u^{\frac{1}{k}}{u}_x\right)}_x+k(k+1)s{u}^{\frac{k+1}{k}}+kq{u}^{\frac{k-1}{k}}$	$\sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}+\sqrt{s}tan\left(\sqrt{s}x\right){\left(u^{\frac{1}{k}}\right)}_x$
		$\sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\sqrt{s}cot\left(\sqrt{s}x\right){\left(u^{\frac{1}{k}}\right)}_x$
		$\sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\sqrt{-s}tanh\left(\sqrt{-s}x\right){\left(u^{\frac{1}{k}}\right)}_x$
		$\sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\sqrt{-s}coth\left(\sqrt{-s}x\right){\left(u^{\frac{1}{k}}\right)}_x$
6	$u_t={\left[exp\left(u\right)u_x\right]}_x+sexp\left(u\right)+qexp(-u)$	$\sigma ={\left[exp\left(u\right)\right]}_{xx}+\sqrt{s}tan\left(\sqrt{s}x\right){\left[exp\left(u\right)\right]}_x$
		$\sigma ={\left[exp\left(u\right)\right]}_{xx}-\sqrt{s}cot\left(\sqrt{s}x\right){\left[exp\left(u\right)\right]}_x$
		$\sigma ={\left[exp\left(u\right)\right]}_{xx}-\sqrt{-s}tanh\left(\sqrt{-s}x\right){\left[exp\left(u\right)\right]}_x$
		$\sigma ={\left[exp\left(u\right)\right]}_{xx}-\sqrt{-s}coth\left(\sqrt{-s}x\right){\left[exp\left(u\right)\right]}_x$
7	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x-\frac{3}{4}s{u}^{-\frac{1}{3}}-\frac{3}{2}q{u}^{\frac{5}{3}}$	$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}+\sqrt{s}tan\left(\sqrt{s}x\right){\left(u^{-\frac{2}{3}}\right)}_x$
		$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}-\sqrt{s}cot\left(\sqrt{s}x\right){\left(u^{-\frac{2}{3}}\right)}_x$
		$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}-\sqrt{-s}tanh\left(\sqrt{-s}x\right){\left(u^{-\frac{2}{3}}\right)}_x$
		$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}-\sqrt{-s}coth\left(\sqrt{-s}x\right){\left(u^{-\frac{2}{3}}\right)}_x$
8	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x$	$\sigma ={\left(u^{-\frac{2}{3}}\right)}_{xx}-\frac{2}{x}{\left(u^{-\frac{2}{3}}\right)}_x+\frac{2}{x^2}{u}^{-\frac{2}{3}}$
		$\sigma ={\left(u^{-\frac{4}{3}}\right)}_{xx}-\frac{6}{x}{\left(u^{-\frac{4}{3}}\right)}_x+\frac{12}{x^2}{u}^{-\frac{4}{3}}$
9	$u_t={\left(u^{\frac{s}{2-s}}{u}_x\right)}_x$	$\sigma ={\left(u^{\frac{s}{2-s}}\right)}_{xx}-\frac{s+2}{2x}{\left(u^{\frac{s}{2-s}}\right)}_x+\frac{s}{x^2}{u}^{\frac{s}{2-s}}$
10	$u_t={\left[exp\left(u\right)u_x\right]}_x$	$\sigma ={\left[exp\left(u\right)\right]}_{xx}-\frac{2}{x}{\left[exp\left(u\right)\right]}_x+\frac{2}{x^2}exp\left(u\right)$
11	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x-\frac{3}{4}q{u}^{\frac{7}{3}}$	$\sigma ={\left(u^{-\frac{4}{3}}\right)}_{xxx}-\frac{3}{x}{\left(u^{-\frac{4}{3}}\right)}_{xx}+\frac{3}{x^2}{\left(u^{-\frac{4}{3}}\right)}_x$
12	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x$	$\sigma ={\left(u^{-\frac{4}{3}}\right)}_{xxx}-\frac{6}{x}{\left(u^{-\frac{4}{3}}\right)}_{xx}+\frac{18}{x^2}{\left(u^{-\frac{4}{3}}\right)}_x-\frac{24}{x^3}{u}^{-\frac{4}{3}}$
13	$u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x-\frac{3}{4}s{u}^{-\frac{1}{3}}-\frac{3}{4}q{u}^{\frac{7}{3}}$	$\begin{array}{ll}\sigma =& {\left(u^{-\frac{4}{3}}\right)}_{xxx}+3\sqrt{s}tan\left(\sqrt{s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}\\ & +\left[1+3{tan}^2\left(\sqrt{s}x\right)\right]s{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$
		$\begin{array}{ll}\sigma =& {\left(u^{-\frac{4}{3}}\right)}_{xxx}-3\sqrt{s}cot\left(\sqrt{s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}\\ & +\left[1+3{cot}^2\left(\sqrt{s}x\right)\right]s{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$
		$\begin{array}{ll}\sigma =& {\left(u^{-\frac{4}{3}}\right)}_{xxx}-3\sqrt{-s}tanh\left(\sqrt{-s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}\\ & +\left[1-3{tanh}^2\left(\sqrt{-s}x\right)\right]s{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$
		$\begin{array}{ll}\sigma =& {\left(u^{-\frac{4}{3}}\right)}_{xxx}-3\sqrt{-s}coth\left(\sqrt{-s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}\\ & +\left[1-3{coth}^2\left(\sqrt{-s}x\right)\right]s{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$

. The corresponding results for Case (iv) is also presented as Case 1 of Table 1.

• Case 2. $a_1\left(x\right)\ne 0$.

$A\left(v\right)=-2b_1{v}^2+k_{1}v+k_2$ is given by solving $4b_1+A^{″}=0$. Since the power or exponential law diffusion is the usual phenomenon, $A\left(v\right)$ and $B\left(v\right)$ take the form of

$$\begin{array}{c}\left(1\right)A\left(v\right)=v,B\left(v\right)=k(k\ne 0);\hfill \\ \left(2\right)A\left(v\right)=v,B\left(v\right)=0;\hfill \\ \left(3\right)A\left(v\right)=v^2,B\left(v\right)=-\frac{1}{2}v.\hfill \end{array}$$

Accordingly,

$$\begin{array}{c}\left(1\right)D\left(u\right)=u^{\frac{1}{k}},g\left(v\right)=v^k;\hfill \\ \left(2\right)D\left(u\right)=exp\left(u\right),g\left(v\right)=ln\left(u\right);\hfill \\ \left(3\right)D\left(u\right)=u^{-\frac{4}{3}},g\left(v\right)=v^{-\frac{3}{2}}\hfill \end{array}$$

will be derived due to the the first two equation of (10). For the case of $A\left(v\right)=v,B\left(v\right)=k(k\ne 0)$, $C\left(v\right)=c_1{v}^2+c_2+c_3$ can be derived from the third one of the determining system (16). Then, the last one of the system becomes

$$\begin{array}{c}\hfill \left[-a_2^{″}\left(x\right)+(2k+1)a_2^2\left(x\right)+2a_2\left(x\right)a_2^{{}^{\prime }}\left(x\right)-c_1{a}_2\left(x\right)\right]v^2+c_3{a}_2\left(x\right)=0.\end{array}$$

$a_2\left(x\right)=0$ and $a_2\left(x\right)\ne 0$ will be respectively considered.

For the case of $a_2\left(x\right)=0$, the determining system (16) is finally simplified as

$$\begin{array}{c}(k+1)\left[a_1^2\left(x\right)-a_1^{{}^{\prime }}\left(x\right)\right]+c_1=0,\hfill \\ a_1^{″}\left(x\right)-2a_1\left(x\right)a_1^{{}^{\prime }}\left(x\right)=0.\hfill \end{array}$$

It is easy to know that $a_1\left(x\right)=s$ and $c_1=-(k+1)s^2$ satisfy the system. In addition, this system can be rewritten as $c_1=-(k+1)s$ and

$$\begin{array}{c}\hfill a_1{\left(x\right)}^2-a_1^{{}^{\prime }}\left(x\right)=s.\end{array}$$

The form of $a_1\left(x\right)$ can be represented as

$$\begin{array}{c}•a_1\left(x\right)=-\frac{1}{x}fors=0;\hfill \\ •a_1\left(x\right)=\sqrt{s}tan\left(\sqrt{s}x\right)or{a}_1\left(x\right)=-\sqrt{s}cot\left(\sqrt{s}x\right)fors>0;\hfill \\ •a_1\left(x\right)=-\sqrt{-s}tanh\left(\sqrt{-s}x\right)or{a}_1\left(x\right)=-\sqrt{-s}coth\left(\sqrt{-s}x\right)fors<0.\hfill \end{array}$$

The corresponding Equation (7) and the admitted CLBS (6) are listed in Table 1.

For the case of $a_2\left(x\right)\ne 0$, $c_3=0$ must hold. The third one of (16) is simplified as

$$\begin{array}{c}\hfill -(k+2)a_2\left(x\right)-2(k+1)\left[a_1^{\prime }\left(x\right)-a_1^2\left(x\right)\right]+2c_1=0.\end{array}$$

It becomes $c_1+a_1^{\prime }\left(x\right)-a_1^2\left(x\right)=0$ in view of $k=-2$. $a_1\left(x\right)a_2\left(x\right)=0$ can be derived from the fourth one of (16), which is contrary to $a_1\left(x\right)=\ne 0$ and $a_2\left(x\right)=\ne 0$. For the case of $k\ne -2$, the fourth equation of (16) becomes

$$\begin{array}{c}\hfill (4k+3)\left[(k+2)a_1^{″}\left(x\right)+2(k+1)a_1^3\left(x\right)-2(2k+3)a_1\left(x\right)a_1^{\prime }\left(x\right)+2c_1{a}_1\left(x\right)\right]=0\end{array}$$

due to

$$\begin{array}{c}\hfill a_2\left(x\right)=\frac{-2(k+1)\left[a_1^{\prime }\left(x\right)-a_1^2\left(x\right)\right]+2c_1}{k+2}.\end{array}$$

Thus, the fourth one of (16) holds for the case of $k=-\frac{3}{4}$. In addition, the last one becomes

$$\begin{array}{c}-5a_1^{‴}\left(x\right)+10a_1\left(x\right)a_1^{″}\left(x\right)+21{a_1^{\prime }}^2\left(x\right)-12a_1^{\prime }\left(x\right)a_1^2\left(x\right)\hfill \\ -53c_1{a}_1^{\prime }\left(x\right)+13c_1{a}_1^2\left(x\right)+a_1^4\left(x\right)=0.\hfill \end{array}$$

Although it is impossible to present the general solution of this nonlinear ODE, we can find that

$$\begin{array}{c}•a_1\left(x\right)=l,a_2\left(x\right)=\frac{2}{9}{l}^2;\hfill \\ •a_1\left(x\right)=-\frac{5}{x}{a}_2\left(x\right)=\frac{8}{x^2}and{c}_1=0;\hfill \\ •a_1\left(x\right)=-\frac{6}{x}{a}_2\left(x\right)=\frac{12}{x^2}and{c}_1=0\hfill \end{array}$$

satisfy the system (16). For the case of $k\ne -\frac{3}{4}$, $a_1\left(x\right)=-\frac{s+2}{2x},a_2\left(x\right)=\frac{s}{x^2}$ and $k=\frac{2-s}{s}$ can be obtained by solving the determining system (16). The corresponding Equation (7) and the admitted CLBS (6) for these cases are all listed in Table 1. We omit the tedious computational procedure for Cases (2) and (3) and just list the corresponding results in Table 1. It is noted that $c,k,s$, and q in Table 1 are arbitrary constants hereafter.

Although the determining system (16) is exactly the special case of the one in [34], we find many particular solutions of System (16) that are not special cases of the results listed in [34]. As a consequence, new forms of Equation (7) and the admitted CLBS (6) are obtained. A similar procedure as for $n=2$ will yield CLBS (6) with $n=3,4,5$ of Equation (7), and these results are also listed in Table 1.

It is interesting to note that some particular equations (including the ones in Cases 1, 8, and 10 of Table 1) admit several CLBSs. These different CLBSs of the same governing equation will naturally lead to different forms of exact solutions.

3. Exact Solutions of Equation (7)

In this section, we construct exact solutions of Equation (7) due to the compatibility of $\sigma =0$ and the governing Equation (7). $\sigma =0$ is exactly the invariant surface condition of the corresponding CLBS (6) admitted by Equation (7). Since $\sigma =0$ can be linearized due to the transformation $f\left(u\right)=v$, the corresponding solution is exactly defined on the space (11) and possesses the form

$$\begin{array}{c}\hfill u=g\left[C_1\left(t\right)f_1\left(x\right)+C_2\left(t\right)f_2\left(x\right)+\cdots +C_n\left(t\right)f_n\left(x\right)\right](2\le n\le 5),\end{array}$$

which is known as FGSV. Substituting the resulting solution into (7) will finally yield that the t-dependent coefficients $C_i\left(t\right)(i=1,2,\cdots ,n)$ satisfy the finite-dimensional dynamical system. Here, we just present several examples to illustrate the reduction procedure.

Example 1.

Equation

$$\begin{array}{c}\hfill u_t={\left(u^{\frac{1}{k}}{u}_x\right)}_x+kq{u}^{\frac{k-1}{k}}\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\frac{1}{x}{\left(u^{\frac{1}{k}}\right)}_x.\end{array}$$

The corresponding solutions are given by

$$\begin{array}{c}\hfill u(x,t)={\left[\alpha \left(t\right)+\beta \left(t\right)x^2\right]}^k,\end{array}$$

where $\alpha \left(t\right)$ and $\beta \left(t\right)$ are listed as below.

(i) For $k\ne -1,$

$$\begin{array}{c}\hfill \alpha \left(t\right)=\frac{[2(2k+1)t-c_1]q}{4(k+1)}+c_2{[2(2k+1)t-c_1]}^{-\frac{1}{2k+1}},\beta \left(t\right)=\frac{1}{c_1-2(2k+1)t}.\end{array}$$

(ii) For $k=-1,$

$$\begin{array}{c}\hfill \alpha \left(t\right)=\left[\frac{1}{2}qln(2t+c_1)\right](2t+c_1),\beta \left(t\right)=\frac{1}{2t+c_1}.\end{array}$$

For the case of $k<0$, the solutions have the asymptotical behavior $u⟶0$ as $x⟶\infty$ and exhibit singularity along the curves $x=\pm {\left[{(-\alpha /\beta )}_{+}\right]}^{1/2}$.

Example 2.

Equation

$$\begin{array}{c}\hfill u_t={\left(u^{\frac{1}{k}}{u}_x\right)}_x+k(k+1)s{u}^{\frac{k+1}{k}}+kq{u}^{\frac{k-1}{k}}\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\sqrt{s}cot\left(\sqrt{s}x\right){\left(u^{\frac{1}{k}}\right)}_x\end{array}$$

for the case of $s>0$ and

$$\begin{array}{c}\hfill \sigma ={\left(u^{\frac{1}{k}}\right)}_{xx}-\sqrt{-s}coth\left(\sqrt{-s}x\right){\left(u^{\frac{1}{k}}\right)}_x\end{array}$$

for the case of $s<0$.

The corresponding separable solutions are respectively

$$\begin{array}{c}\hfill u(x,t)={\left[\alpha \left(t\right)+\beta \left(t\right)cos\left(\sqrt{s}x\right)\right]}^k\end{array}$$

and

$$\begin{array}{c}\hfill u(x,t)={\left[\alpha \left(t\right)+\beta \left(t\right)cosh\left(\sqrt{-s}x\right)\right]}^k,\end{array}$$

where $\alpha \left(t\right)$ and $\beta \left(t\right)$ satisfy two-dimensional dynamical system

$$\begin{array}{c}\hfill {\alpha }^{\prime }=(k+1)s{\alpha }^2+ks{\beta }^2+q,{\beta }^{\prime }=(2k+1)s\alpha \beta .\end{array}$$

For $s>0$, the solutions are x-periodic with the period $2\pi /\sqrt{s}$.

Example 3.

Equation

$$\begin{array}{c}\hfill u_t={\left[exp\left(u\right)u_x\right]}_x+sexp\left(u\right)+qexp(-u)\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left[exp\left(u\right)\right]}_{xx}+\sqrt{s}tan\left(\sqrt{s}x\right){\left[exp\left(u\right)\right]}_x\end{array}$$

for the case of $s>0$ and

$$\begin{array}{c}\hfill \sigma ={\left[exp\left(u\right)\right]}_{xx}-\sqrt{-s}tanh\left(\sqrt{-s}x\right){\left[exp\left(u\right)\right]}_x\end{array}$$

for the case of $s<0$.

The corresponding separable solutions are respectively

$$\begin{array}{c}\hfill u(x,t)=ln\left[\alpha \left(t\right)+\beta \left(t\right)sin\left(\sqrt{s}x\right)\right]\end{array}$$

and

$$\begin{array}{c}\hfill u(x,t)=ln\left[\alpha \left(t\right)+\beta \left(t\right)sinh\left(\sqrt{-s}x\right)\right],\end{array}$$

where $\alpha \left(t\right)$ and $\beta \left(t\right)$ satisfy two-dimensional dynamical system

$$\begin{array}{c}\hfill {\alpha }^{\prime }=s{\alpha }^2+q,{\beta }^{\prime }=s\alpha \beta .\end{array}$$

$\alpha \left(t\right)$ and $\beta \left(t\right)$ are listed as below.

(i) For $qs>0,$

$$\begin{array}{c}\hfill \alpha \left(t\right)=\frac{\sqrt{qs}}{s}tan\left[\sqrt{qs}(t+c_1)\right],\beta \left(t\right)=c_{2}sec\left[\sqrt{qs}(t+c_1)\right].\end{array}$$

(ii) For $qs<0,$

$$\begin{array}{c}\hfill \alpha \left(t\right)=-\frac{\sqrt{-qs}}{s}tanh\left[\sqrt{-qs}(t+c_1)\right],\beta \left(t\right)=c_2\mathrm{sech}\left[\sqrt{-qs}(t+c_1)\right].\end{array}$$

(iii) For $q=0,$

$$\begin{array}{c}\hfill \alpha \left(t\right)=\frac{1}{c_1-st},\beta \left(t\right)=\frac{c_2}{c_1-st}.\end{array}$$

Example 4.

Equation

$$\begin{array}{c}\hfill u_t={\left(u^{\frac{s}{2-s}}{u}_x\right)}_x\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{\frac{s}{2-s}}\right)}_{xx}-\frac{s+2}{2x}{\left(u^{\frac{s}{2-s}}\right)}_x+\frac{s}{x^2}{u}^{\frac{s}{2-s}}.\end{array}$$

The corresponding separable solutions are given as follows.

(i) For $s=4,$

$$\begin{array}{c}\hfill u(x,t)={\left[\frac{c_2+\frac{1}{2}ln(c_1-t)}{c_1-t}{x}^2+\frac{1}{c_1-t}{x}^{2}lnx\right]}^{-\frac{1}{2}}.\end{array}$$

(ii) For $s\ne 4,$

$$\begin{array}{c}\hfill u(x,t)={\left\{\frac{1}{\frac{2(s-4)}{s}t+c_1}{x}^2+c_2{\left[\frac{2(s-4)}{s}t+c_1\right]}^{\frac{s(s-6)}{8}}{x}^{\frac{s}{2}}\right\}}^{\frac{2-s}{s}}.\end{array}$$

Example 5.

Equation

$$\begin{array}{c}\hfill u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{-\frac{4}{3}}\right)}_{xxx}-\frac{6}{x}{\left(u^{-\frac{4}{3}}\right)}_{xx}+\frac{18}{x^2}{\left(u^{-\frac{4}{3}}\right)}_x-\frac{24}{x^3}{u}^{-\frac{4}{3}}.\end{array}$$

The corresponding solutions are given by

$$\begin{array}{c}\hfill u(x,t)={\left[\frac{x^2+c_2{x}^3}{c_1+t}+\left(\frac{c_2^2}{4(c_1+t)}+c_3{(c_1+t)}^2\right)x^4\right]}^{-\frac{3}{4}}.\end{array}$$

It is easy to see that the solutions have the asymptotical behavior $u⟶0$ as $x⟶\infty$ and $u⟶0$ as $t⟶+\infty$.

Example 6.

Equation

$$\begin{array}{c}\hfill u_t={\left(u^{-\frac{4}{3}}{u}_x\right)}_x-\frac{3}{4}s{u}^{-\frac{1}{3}}-\frac{3}{4}q{u}^{\frac{7}{3}}\end{array}$$

admits the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{-\frac{4}{3}}\right)}_{xxx}-3\sqrt{s}cot\left(\sqrt{s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}+s\left[1+3{cot}^2\left(\sqrt{s}x\right)\right]{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$$

for $s>0$ and the CLBS

$$\begin{array}{c}\hfill \sigma ={\left(u^{-\frac{4}{3}}\right)}_{xxx}-3\sqrt{-s}coth\left(\sqrt{-s}x\right){\left(u^{-\frac{4}{3}}\right)}_{xx}+s\left[1-3{coth}^2\left(\sqrt{-s}x\right)\right]{\left(u^{-\frac{4}{3}}\right)}_x\end{array}$$

The corresponding solutions are given respectively as follows.

(i) For $s>0$,

$$\begin{array}{c}\hfill u(x,t)={\left[\alpha \left(t\right)+\beta \left(t\right)cos\left(\sqrt{s}x\right)+\gamma \left(t\right)cos\left(2\sqrt{s}x\right)\right]}^{-\frac{3}{4}}.\end{array}$$

(ii) For $s<0$,

$$\begin{array}{c}\hfill u(x,t)={\left[\alpha \left(t\right)+\beta \left(t\right)cosh\left(\sqrt{-s}x\right)+\gamma \left(t\right)cosh\left(2\sqrt{-s}x\right)\right]}^{-\frac{3}{4}}.\end{array}$$

The unknown functions in the solutions satisfy three-dimensional dynamical system

$$\begin{array}{c}{\alpha }^{\prime }=s{\alpha }^2-\frac{3}{8}s{\beta }^2-3s{\gamma }^2+q,\hfill \\ {\beta }^{\prime }=s\alpha \beta -3s\beta \gamma ,\hfill \\ {\gamma }^{\prime }=\frac{3}{8}s{\beta }^2-2s\alpha \gamma .\hfill \end{array}$$

It is noted that the constant $c_1$ can be removed from all the resulting solutions because Equation (7) is invariant with respect to time translation. The exact solutions of Example 1 were found in [35] for the first time and later rediscovered in [11]. The exact solutions presented in Example 3 are also constructed due to conditional symmetry [36]. In addition, the exact solutions of Examples 2 and 3 are both particular cases of those derived in [11,31]. The exact solutions listed in Example 5 are obtainable from the dynamical system presented in [13].

4. Conclusions

We obtained the classification of the quasi-linear diffusion Equation (7) due to the CLBS method, which provided a symmetry interpretation to the FGSV. The variant forms of (7) admitting the CLBS (6) were presented. As a consequence, we constructed exact solutions in separable form (4) to the resulting equations due to the corresponding symmetry reductions. Most of these solutions could not be obtained by the other symmetry reduction methods.

This method could also be used to consider other types of PDEs, including KdV-type equations, fourth-order equations, etc. Moreover, it was also effective to deal with PDE systems, such as all kinds of diffusion systems. In addition, the discussion of multi-dimensional evolution equations, including all kinds of multi-dimensional diffusion equations in 2D, will be involved in our further study.