# Second-Order Conditional Lie-Bäcklund Symmetry and Differential Constraint of Radially Symmetric Diffusion System

Authors: Jianping Wang; Huijing Ba; Yaru Liu; Longqi He; Lina Ji
Journal: Advances in Mathematical Physics
 (2021)

---

## Abstract

The classifications and reductions of radially symmetric diffusion system are studied due to the conditional Lie-Bäcklund symmetry method. We obtain the invariant condition, which is the so-called determining system and under which the radially symmetric diffusion system admits second-order conditional Lie-Bäcklund symmetries. The governing systems and the admitted second-order conditional Lie-Bäcklund symmetries are identified by solving the nonlinear determining system. Exact solutions of the resulting systems are constructed due to the compatibility of the original system and the admitted differential constraint corresponding to the invariant surface condition. For most of the cases, they are reduced to solving four-dimensional dynamical systems.


---

## Body

## 1. Introduction

In the latter part of the 19th century, Lie [1] introduced the notion of continuous groups, now known as Lie groups and classical symmetry, which at once unified a variety of special techniques designed to solve the ordinary differential equation (ODE) including the integrating factor, reduction of order, undetermined coefficients, Laplace transform, etc. These special methods were, in fact, all special cases of a general integration procedure based on the invariance of the differential equation (DE) under a continuous group of symmetry.To date, various generalizations to the concept of symmetry groups for nonlinear PDEs have been proposed. Noether [2] introduced the notion of the Lie-Bäcklund symmetry. Ovsiannikov [3] developed the method of partially invariant solutions. Bluman and Cole [4] suggested the so-called nonclassical method. Fushchych et al. [5] presented the conditional symmetry, which is an extension of the nonclassical method. Olver and Rosenau [6, 7] generalized the nonclassical method to the weak symmetry method and the side condition method. Nucci [8] gave the iteration of the nonclassical method. Fokas and Liu [9] and Zhdanov [10] independently carried out the conditional Lie-Bäcklund symmetry (CLBS). There are a lot of symmetry-related methods such as the direct and modified direct methods [11, 12], the ansatz-based method [13–15], the sign-invariant method [16–18], and the invariant subspace method [19, 20].The merging of the extension ideas of Lie-Bäcklund symmetry and conditional symmetry was independently carried out by Fokas and Liu [9] and Zhdanov [10]. In [9], Fokas and Liu introduced the notion of generalized conditional symmetry and constructed multishock and multisoliton solutions of certain nonintegrable equations. The reduction theorem of conditional Lie-Bäcklund symmetry (CLBS) which ensures that a PDE will be reduced to a system of ODEs was presented by Zhdanov in [10], where the corresponding symmetry reductions for the nonlinear heat conductivity equation were also displayed. Moreover, one-to-one correspondence relation between reducibility of a given evolution equation to a system of ODEs and its CLBS was established by Zhdanov in [21]. In fact, generalized conditional symmetry [9], CLBS [10], and higher conditional symmetry [21] are always the same.It is proven that CLBS are very effective to construct exact solutions and find symmetry reductions for different types of nonlinear diffusion equations [22–32]. CLBSs and symmetry reductions of nonintegrable equations can be referred to in [9, 33–35]. The reduction of initial-value problem for evolution equations to Cauchy problem for system of ODEs which can be fully characterized in terms of CLBS is shown by Zhdanov et al. in [21, 36–38]. The studies of CLBS for the evolution system are discussed by Andreytsev in [39] and by Sergyeyev in [40, 41], where the complete description of the evolution system which admits a given CLBS is provided. The reduction theorem for CLBS of the evolution system is constructed in [42], which generalizes the one for scalar evolution equation [10]. The studies of classifications and reductions for the two-component diffusion system can be referred to in [42–44]Olver [45] showed that the invariant surface condition corresponding to CLBS can be regarded as a differential constraint (DC) compatible with the initial equation within the framework of empiric compatibility theory. The term “conditional” is explained by the fact that attaching an additional differential equation called DC to the original PDE and the DC is right with the invariant surface equation corresponding to the admitted symmetry. In fact, conditional invariance criterion is nothing but a compatibility condition of the combined system including original PDE and additional invariant surface equation. The relations between conditional symmetry, reduction, and compatibility of the combined system were discussed in [45, 46]. The related papers addressing this problem include [33, 47, 48]. CLBS can be reformulated within the framework of the DC method. Olver [6, 7, 45], Kaptsov [49], Levi and Winternitz [50], and Pucci and Saccomandi [51] conclude that many reduction methods such as conditional symmetry, partial invariance, variable of separation, and direct method can be understood by using the technicalities of the method of DC. The method of CLBS provides an appropriate symmetry background for the method of DC. The base of symmetry reduction for CLBS is the fact that the corresponding invariant surface condition is formally compatible with the governing system, which is extensively discussed in [48], where it is shown that the problem of discussing the DC of the evolution system is equivalent to studying the CLBS of this system. The equivalence relation between CLBS, DC, and direct reduction is also discussed in our recent paper [43]The procedure for determining whether or not a given Lie-Bäcklund vector field is conditionally invariant of the considered equation is straightforward; however, the determination of the most general CLBSs admitted by a given PDE is a very difficult, if not impossible, problem since the associated determining system is an overdetermined nonlinear system of PDEs. Nevertheless, as is known, even finding particular CLBSs can lead to new explicit solutions of the considered equation. In practice, the principle direction of such research is to content oneself with finding CLBS in particular cases, and these cases must be chosen using additional considerations. It has been proven that CLBSs related with invariant subspaces [22–24], sign-invariants [25–30], and separation of variables [31, 32] are very effective to study the classifications and reductions of second-order nonlinear diffusion equations.In this paper, we will study the second-order CLBS with the characteristic(1)η1=urr+Huur2+Gr,uur+Fr,u,η2=vrr+Lvvr2+Mr,vv1+Nr,v,of a nonlinear radially symmetric diffusion system
(2)ut=1rn−1rn−1ukurr+Pu,v,vt=1rn−1rn−1vlvrr+Qu,v,which is equivalent to studying the second-order DC
(3)η1=urr+Huur2+Gr,uur+Fr,u=0,η2=vrr+Lvvr2+Mr,vv1+Nr,v=0,of Equation (2). It is noted that k and l in (2) are both arbitrary real constants. The form of CLBS (1) generalizes the one for a nonlinear diffusion system in [43], and the corresponding second-order CLBS for the scalar diffusion equation can provide symmetry interpretation for first-order Hamilton-Jacobi sign-invariant for the considered equation [25–30]. The discussion about second-order CLBS (1) of system (2) for the case of n=1 is referred to [43].The layout of this paper is listed here. Section2 is devoted to necessary definitions and notations about CLBS and DC of evolution system. Second-order CLBS (1) and DC (3) of system (2) are displayed in Section 3. Exact solutions of system (2) are constructed in Section 4. Conclusions and remarks are given in the last section.


## 2. Preliminaries

The one-parameter Lie-Bäcklund group of infinitesimal transformations(4)u~i=ui+εηir,t,ui,u1i,⋯,ulii+Oε2,u~ti=ui+εDtηir,t,ui,u1i,⋯,ulii+Oε2,u~ri=ui+εDrηir,t,ui,u1i,⋯,ulii+Oε2,⋯is generated by the Lie-Bäcklund vector field (LBVF)
(5)V=∑i=1mηir,t,ui,u1i,⋯,ulii∂∂ui+Drηir,t,ui,u1i,⋯,ulii∂∂u1i+Dtηir,t,ui,u1i,⋯,ulii∂∂uti+Dr2ηir,t,ui,u1i,⋯,ulii∂∂u2i+⋯.Definition 1 (see [3]).
The evolutionary vector field (5) is said to be a Lie-Bäcklund symmetry of the evolution system
(6)uti=Fit,r,u1,u11,⋯,uk11,⋯,um,u1m,⋯,ukmm,i=1,2,⋯,m,if
(7)Vuti−FiM=0,i=1,2,⋯,,m,where M denotes the set of all differential consequences of the system (6).Definition 2 (see [9, 10]).
The evolutionary vector field (5) is said to be a CLBS of (6) if
(8)Vuti−FiM∩Lr=0,i=1,2,⋯,m,where Lr denotes the set of all differential consequences of invariant surface condition
(9)ηir,t,ui,u1i,⋯,ulii=0,i=1,2,⋯,m,with respect to r.
In fact, the invariant criterion (8) is reduced to
(10)Vuti−FiM∩Lr=Dtηi−∑l=1m∑j=1klDxjηlFujliM∩Lx=DtηiM∩Lr=0.
The fact that LBVF (5) is a CLBS of system (1) leads to the compatibility of the invariant surface condition (9) and the governing system (1). For evolution system (1), the invariant criterion (10) is exactly the sufficient condition which is used to construct the DC of (1).Definition 3 (see [52]).
The differential constraints (9) and the evolution system (1) satisfy the compatibility condition if
(11)DtηiMr∩Lr=0,i=1,2,⋯,m,where Mr denotes the set of all differential consequences of the system (6) with respect to r.
The calculation of CLBS admitted can be divided into four steps. Firstly, we need to compute total derivatives ofηi=0 with respect to t for i=1,2,⋯,m, which yields
(12)ηti+ηuiiuti+ηu1iiu1ti+⋯+ηuliiiuliti=0.
It is noted that the subscripts denote partial derivatives with respect to the indicated variables. The next step is to eliminate all derivativesujti by using differential consequences Drjuti−Fi=0, which implies
(13)ηti+ηuiiFi+ηu1iiDxFi+⋯+ηuliiiDxliFi=0.
In fact, the left side of (13) can be simplified as a polynomial about u1i,u2i,⋯,uli−1i by substituting all higher-order derivatives of ui with respect to r by lower-order ones due to differential consequences Drjηi=0 for j=0,1,2,⋯. Consequently, equating the coefficients of resulting polynomials to zero will yield the so-called determining system for the undetermined parts in the governing system (6) and the corresponding characteristic system (9). Finally, solving this system leads to the form of CLBS (9) admitted by the system (6). In general, solving the resulting nonlinear determining system is as difficult as solving the original evolution system. However, even finding one particular CLBS admitted by the considered evolution system (6) will lead to the corresponding reductions of the original system. Thus, we will find particular solutions of the overdetermined system to determine the CLBS admitted by the system (6), which are proved to be very powerful to study classifications and reductions of scalar nonlinear diffusion equations [22–32].
In full analog with the classical symmetry reduction method, conditional symmetries can also be used to perform dimensional reductions of the governing system. The invariant solutions are defined by the characteristic system (9). The general solution to system (9) can be (locally) written as
(14)ui=fit,r,ϕ1it,ϕ2it,⋯,ϕliit,where ϕjit are arbitrary smooth functions. ηi in (9) depends only on ui and its higher-order derivatives, so the ansatz ui in (14) is independent with each other. Consequently, the studies of direct reductions, DCs and CLBSs for evolution system are about the same as that for scalar evolution equation. Thus, the equivalence relation [21, 45, 48] between direct reductions, DCs and CLBSs for scalar evolution equation holds true for the case of evolution system. It is concluded that reducibility of the evolution system (1) to the ODE system due to the ansatz (14), CLBS (5) and DC (9) of system (1) is, in some sense, equivalent [42, 43].


## 3. CLBS (1) of System (2)

Since the diffusion system (2) admits of CLBS with the characteristic (1), the invariant criterion (10) can be simplified as
(15)Dtη1M∩Lr=−u3H′′+4uH−3ku2H′+4k1−kuH+5ku2H2−2u3H3+kk−1k−2uk−3ur4+−u2Guu+2uH−3kuGu+4u2H′+9kuH−4u2H2+5k−5k2G+kn−1rk−1−uHuk−2ur3+Puu+H′P+HPu−uk−2u2Fuu+3kuFu−4u2H′F−24kuH−2u2H2+3k−3k2F+22k−uHuGr−G2−2u2GGu+2u2Gru+2kn−1urG+1rur2+2Puv+HPvurvr+Pvv−PvLvr2+G−MPvvr+PGu+k7FG−4Fr−3n−1rFuk−1+2FGu+2GGr−Grr−2Fru−4HFG+n−1r2r2+1rG−Grukur+FuP−PuF−PvN+3kuk−1F2+2FGr−HF2−Frr+n−1r2rF−Fruk=0,Dtη2M∩Lr=−v3L′′+4vL−3lv2L′+4l1−lvL+5lv2L2−2v3L3−v3L′′+4vL−3lv2L′+4l1−lvL+5lv2L2−2v3L3+ll−1l−2vl−3vr4+−v2Mvv+2vL−3lvMv+4v2L′+9lvL−4v2L2+5l−5l2M+ln−1rl−1−vLvl−2vr3+Qvv+QL′+LQv−vl−2v2Nvv+3lvNv−4v2L′N−24lvL−2v2L2+3l−3l2N+22l−vLvMr−M2−2v2MMv+2v2Mrv+2ln−1vrM+1rvr2+2Quv+LQuurvr+Quu−QuHur2+M−GQuur+QMv+l7MN−4Nr−3n−1rNvl−1+2NMv+2MMr−Mrr−2Nrv−4LMN+n−1r2r2+1rM−Mrvlvr+NvQ−QuF−QvN+3lvl−1N2+2NMr−LN2−Nrr+n−1r2rN−Nrvl=0.The vanishing of all the coefficients of the two polynomials aboutur and vr yields the nonlinear system of determining equations
(16)−u3H′′+4uH−3ku2H′+4k1−kuH+5ku2H2−2u3H3+kk−1k−2=0,−u2Guu+2uH−3kuGu+4u2H′+9kuH−4u2H2+5k−5k2G+kn−1rk−1−uH=0,Puu+H′P+HPu−uk−2u2Fuu+3kuFu−4u2H′F−24kuH−2u2H2+3k−3k2F+22k−uHuGr−G2−2u2GGu+2u2Gru+2kn−1urG+1r=0,Puv+HPv=0,Pvv−PvL=0,G−MPv=0,PGu+k7FG−4Fr−3n−1rFuk−1+2FGu+2GGr−Grr−2Fru−4HFG+n−1r2r2+1rG−Gruk=0,FuP−PuF−PvN+3kuk−1F2+2FGr−HF2−Frr+n−1r2rF−Fruk=0−v3L′′+4vL−3lv2L′+4l1−lvL+5lv2L2−2v3L3+ll−1l−2=0,−v2Mvv+2vL−3lvMv+4v2L′+9lvL−4v2L2+5l−5l2M+ln−1rl−1−vL=0,Qvv+QL′+LQv−vl−2v2Nvv+3lvNv−4v2L′N−24lvL−2v2L2+3l−3l2N+22l−vLvMr−M2−2v2MMv+2v2Mrv+2ln−1vrM+1r=0,Quv+LQu=0,Quu−QuH=0,M−GQu=0,QMv+l7MN−4Nr−3n−1rNvl−1+2NMv+2MMr−Mrr−2Nrv−4LMN+n−1r2r2+1rM−Mrvl=0,NvQ−QuF−QvN+3lvl−1N2+2NMr−LN2−Nrr+n−1r2rN−Nrvl=0.It is impossible to present the general solutions of the nonlinear determining system since system (16) is a coupled nonlinear system of PDEs. The workable way is to construct particular solutions of the nonlinear determining system (16). The determined CLBS (1) will lead to symmetry reductions of the governing system (2).SettingHu=h/u into the first one of system (16), we can obtain the algebraic equation
(17)h−kh−k+1h−k2+1=0.Thus, we can list three special solutions(18)Hu=ku,k−1u,k/2−1u.The similar procedure gives(19)Lv=lv,l−1v,l/2−1v.We just consider six different cases(20)Hu=ku,Lv=lv,Hu=ku,Lv=l−1v,Hu=ku,Lv=l/2−1v,Hu=k−1u,Lv=l−1v,Hu=k−1u,Lv=l/2−1v,Hu=k/2−1u,Lv=l/2−1v,for further study due to the symmetrical form of system (2) and the admitted CLBS (1).Since system (2) is a coupled one, the sixth one and the fourteenth one of the determining system (16) implies
(21)Mr,u=Gr,v=Ar.As a consequence,Ar=n−1/r is finally determined due to the second one and the tenth one of system (16) for case i. Substituting Hu=k/u and Lv=l/v into the fifth one and the thirteenth one of the determining system (16), we get Pu,v and Qu,v, respectively, which satisfy the linear differential equation
(22)Pvv−lvPv=0,Quu−kuQu=0.Solving the differential equation aboutPu,v and Qu,v, we know that three subcases
(23)Pu,v=P1u+P2uv1+l,Qu,v=Q1v+Q2vu1+kk≠−1,l≠−1,Pu,v=P1u+P2uv1+l,Qu,v=Q1v+Q2vlnuk=−1,l≠−1,Pu,v=P1u+P2ulnv,Qu,v=Q1v+Q2vlnuk=−1,l=−1,will be considered. For the case of k≠−1,l≠−1, P2u=p1u−k and Q2v=p2v−l can be obtained by solving the fourth one and the twelfth one of system (16). Then, we conclude that Fr,u and Nr,v respectively satisfy the linear PDE
(24)Fru+2kuFr=0,Nrv−2lvNr=0.It is easy to know that(25)Fr,u=F1u+F2ru−2k,Nr,v=N1v+N2rv−2l.The eighth one of system (16) becomes
(26)ku−3k−1F22r+−u−2kP′1u−2ku−2k−1P1u+2ku−k−1F1uF2r−u−kF2′′r−n−1ru−kF2′r+F′1uP1u−P′1uF1u+kuk−1F12u+p1u−kF1′u−ku−3k−1F2r+ku−k−1F1uvl+1−p1l+1u−kN2rv−l−p1l+1u−kvlN1v=0.Consideringp1≠0 and l≠−1, we can set N1v=c1v+c2v−2l+c3v−l. As a consequence, we got
(27)ku−3k−1F22r+−u−2kP′1u−2ku−2k−1P1u+2ku−k−1F1uF2r−u−kF2′′r−n−1ru−kF2′r+F′1uP1u−P′1uF1u+kuk−1F12u−p1l+1c3u−k+p1u−kF1′u−ku−3k−1F2r+ku−k−1F1u−l+1c1u−kv1+l−p1l+1u−kN2r+c2v−l=0.It is noted that the above polynomial aboutv include two terms for the case of l=−1/2. Firstly, we consider the case of k≠−1/2. The polynomial is zero will yield that N2r=−c2. In addition, F2r=c4 and
(28)F1u=−c4u−2k+l+1c1k+1u+s1u−k,is derived due to
(29)u−kF1′u−ku−3k−1F2r+ku−k−1F1u−l+1c1u−k=0.Solving the third one and the eleventh one of system (16), we induce
(30)P1u=a1u+b1u−k+s1+l+1c1k+1u1+k,Q1v=a2v+b2v−l+c3+c1v1+l.The eighth one and the sixteenth one of the system (16) become
(31)ks1a1+a1s1−c1b1+p1c3−c1lb1+p1lc3=0,−c1b2+c3la2+p2s1k+p2s1+a2c3−b2lc1=0,which finally identify the undetermined functions in diffusion system (2) and the admitted CLBS (1). The corresponding results are listed as items 1-3 of Table 1. Similar discussion as above for the case of k=−1/2 will present the corresponding CLBS (1) and the governing system (2), which are listed as items 4-6 of Table 1. The results for k=−1,l≠−1 and k=−1,l=−1 are also listed in Table 1. We omit the tedious computational procedure for other cases of Hu and Lv and just list the corresponding results in Table 1. It is noted that the software Maple is used for calculations.Table 1
CLBS (1) of RD system (2).
No.RD system (2)CLBS (1)1ut=1rn−1rn−1ukurr+a1u+b1u−k+p1u−kv1+l,vt=1rn−1rn−1vlvrr+a2v+b2v−l+p2v−lu1+k.η1=urr+kuur2+n−1rur,η2=vrr+lvvr2+n−1rvr.2ut=1rn−1rn−1ukurr+a1u−k−1+lp1s21+ks1u+s1+p1u−kv1+l,vt=1rn−1rn−1vlvrr+a2v−l−1+kp2s11+ls2v+s2+p2v−lu1+k.η1=urr+kuur2+n−1rur+s1u−k,η2=vrr+lvvr2+n−1rvr+s2v−l.3ut=1rn−1rn−1ukurr+a1u+1+ka1s1+1+lp1s31+ls2u−k+1+ls21+ku1+k+s1+p1u−kv1+l,vt=1rn−1rn−1vlvrr+a2v+1+kp2s1+1+la2s31+ls2v−l+s2v1+l+s3+p2v−lu1+k.η1=urr+kuur2+n−1rur+1+ls21+ku+s1u−k,η2=vrr+lvvr2+n−1rvr+s2v+s3v−l.4ut=1rn−1rn−1u−1/2urr+a1u+b1uv,vt=1rn−1rn−1v−1/2vrr+a2v+b2uv.η1=urr−12uur2+n−1rur+4n−4r2u,η2=vrr−12vvr2+n−1rvr+4n−4r2v.5ut=1rn−1rn−1u−1/2urr+a1u+s1+2s3u+s32−a1s1s2uv,vt=1rn−1rn−1v−1/2vrr+a2v+s2+2s3v+s32−a2s2s1uv.η1=urr−12uur2+n−1rur+s3u+s1u+4n−4r2u,η2=vrr−12vvr2+n−1rvr+s3v+s2v+4n−4r2v.6ut=1rn−1rn−1u−1urr+a1u+b1ulnu+p1uv1+l,vt=1rn−1rn−1vlvrr+a2v+b2v−l+p2v−llnu.η1=urr−1uur2+n−1rur,η2=vrr+lvvr2+n−1rvr.7ut=1rn−1rn−1u−1urr+a1ulnu+1+lp1s1s2u+s2lnu+p1uv1+l,vt=1rn−1rn−1vlvrr+a2v+1+la2s1s2v−l+s21+lv1+l+s1+p2v−lu1+k.η1=urr−1uur2+n−1rur+s2ulnu,η2=vrr+lvvr2+n−1rvr+s21+lv+s1v−l.8ut=1rn−1rn−1u−1urr+a1u−1+lp1s1−a1s2s2ulnu+s2lnu+s3+p1uv1+l,vt=1rn−1rn−1vlvrr+a2v+p2v−l+s21+lv1+l+s1−1+la2s1−p2s2s3v−l.η1=urr−1uur2+n−1rur+s2ulnu+s3u,η2=vrr+lvvr2+n−1rvr+s21+lv+s1v−l.ut=1rn−1rn−1u−1urr+a1u+b1ulnu+p1ulnv,vt=1rn−1rn−1v−1vrr+a2v+b2vlnv+p2vlnu.η1=urr−1uur2+n−1rur,η2=vrr−1vvr2+n−1rvr.10ut=1rn−1rn−1u−1urr+a1u+b1ulnu−b1s3s1ulnv+s3,vt=1rn−1rn−1v−1vrr+a2v+b2vlnv−b2s1s3vlnu+s1.η1=urr−1uur2+n−1rur+s3u,η2=vrr−1vvr2+n−1rvr+s1v.11ut=1rn−1rn−1u−1urr+b1s3+p1s1s2u+b1ulnu+s2lnu+s3+p1ulnv,vt=1rn−1rn−1v−1vrr+b2s1+p2s3s2v+b2vlnv+s2lnv+s1+p2vlnu.η1=urr−1uur2+n−1rur+s2ulnu+s3u,η2=vrr−1vvr2+n−1rvr+s2vlnv+s1v.12ut=1rru−1/2urr+a1u+p1uv,vt=1rrv−1vrr+a2v+p2v2u.η1=urr−12uur2+1rur−8r2u,η2=vrr−2vvr2+1rvr+4r2v.13ut=1rn−1rn−1ukurr+a1u+b1u1−k+p1u1−kvl,vt=1rn−1rn−1vlvrr+a2v+b2v1−l+p2v1−lukη1=urr+k−1uur2−1rur,η2=vrr+l−1vvr2−1rvr.14ut=1rn−1rn−1u23−nn−4urr+p1s23−n+s12n−3s1u+p1u3n−10n−4v23−nn−4,vt=1rn−1rn−1v23−nn−4vrr+p2s13−n+s22n−3s2v+p2v3n−10n−4u23−nn−4.η1=urr−3n−10n−4uur2+3n−3rur+s1u3n−10/n−4−n−42r2u,η2=vrr−3n−10n−4vvr2+3n−3rvr+s2v3n−10/n−4−n−42r2v.15ut=1rn−1rn−1ulurr+a1u+p1u1−lvl,vt=1rn−1rn−1vlvrr+a2v+p2v1−lul.η1=urr+l−1uur2+n−3l−1l+1rur−2n−2l+1r2u,η2=vrr+l−1vvr2+n−3l−1l+1rvr−2n−2l+1r2v.16ut=1rru−1urr+a1u+p1u2v,vt=1rrv−1vrr+a2v+p2v2u.η1=urr−2uur2+srur+2s+1r2u,η2=vrr−2vvr2+srvr+2s+1r2v.17ut=1rn−1rn−1u23−nn−4urr+a1u+p1u3n−10n−4v23−nn−4,vt=1rn−1rn−1v23−nn−4vrr+a2v+p2v3n−10n−4u23−nn−4.η1=urr−3n−10n−4uur2+3n−3rur−n−42r2u,η2=vrr−3n−10n−4vvr2+3n−3rvr−n−42r2v.18ut=1rn−1rn−1ukurr+a1u+b1u1−k+p1u1−kv−2/n+2,vt=1rn−1rn−1v−4/n+2vrr+a2v+b2vn+4/n+2+p2vn+4/n+2uk.η1=urr+k−1uur2−1rur,η2=vrr−n+4n+2vvr2−1rvr.19ut=1rn−1rn−1u−4/n+2urr+a1u+b1un+4/n+2+p1un+4/n+2v−2/n+2,vt=1rn−1rn−1v−4/n+2vrr+a2v+b2vn+4/n+2+p2vn+4/n+2u−2/n+2.η1=urr−n+4n+2uur2−1rur,η2=vrr−n+4n+2vvr2−1rvr.20ut=1rn−1rn−1u23−n/n−4urr+a1u+p1u2n−7/n−4v3−n/n−4,vt=1rn−1rn−1v23−n/n−4vrr+a2v+p2v2n−7/n−4u3−n/n−4.η1=urr−2n−7n−4uur2+n−3rur+n−4r2u,η2=vrr−2n−7n−4vvr2+n−3rvr+n−4r2v.21ut=1rru−1urr+a1u+p1u3/2v,vt=1rrv−1vrr+a2v+p2v3/2u.η1=urr−32uur2−1rur+s1r2u,η2=vrr−32vvr2−1rvr+s1r2v.


## 4. Exact Solutions of System (2)

Since the admitted CLBS (1) of system (2) rightly corresponds to the DC (3), the exact solutions of system (2) listed in Table 1 can be derived due to the compatibility of the governing system (2) and the admitted DC (3) corresponding to the invariant surface of CLBS (1). One first solves the two ODEs in (3) to determine the form of ur,t and vr,t, which are both functions about r with t-dependent integration constants. Substituting the resulting ur,t and vr,t into the original system (2), we can finally determine the time evolution integration constants. Here, we just present several examples to illustrate the reduction procedure.Example 4.
System(32)ut=1rn−1rn−1ukurr+a1u−k−1+lp1s21+ks1u+s1+p1u−kv1+l,vt=1rn−1rn−1vlvrr+a2v−l−1+kp2s11+ls2v+s2+p2v−lu1+k,admits of CLBS
(33)η1=urr+kuur2+n−1rur+s1u−k,η2=vrr+lvvr2+n−1rvr+s2v−l.
The exact solutions of (32) are listed as below.
(i)
Forn≠2,(34)ur,t=αtr2−n+βt−1+ks12nr21/1+k,vr,t=ϕtr2−n+ψt−1+ls22nr21/1+l,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(35)α′=−p1s21+ls1α+1+kp1ϕ,β′=−p1s21+ls1β+1+kp1ψ+1+ka1,ϕ′=−p2s11+ks2ϕ+1+lp2α,ψ′=−p2s11+ks2ψ+1+lp2β+1+la2.
The solutions of this linear system of PDEs are presented as follows.(i)
Forp2≠−p1s221+l/s121+k,
(36)αt=c1+c2exp−p1s221+l+p2s121+ks1s2t,βt=−s1s2c3p1s221+l+p2s121+kexp−p1s221+l+p2s121+ks1s2t+s11+kp1s2a21+l+p2s1a11+kp1s221+l+p2s121+kt+c4,ϕt=−p2s1c2s2p1exp−p1s221+l+p2s121+ks1s2t+s21+lc1s11+k,ψt=p2s12c3p1p1s221+l+p2s121+kexp−p1s221+l+p2s121+ks1s2t+1+lp1s2a21+l+p2s1a11+ks2p1s221+l+p2s121+kt+s21+la2s1−a1s2p1s221+l+p2s121+k+s21+lc4s11+k.(ii)
Forp2=−p1s221+l/s121+k,
(37)αt=c1t+c2,βt=p11+l1+ks1a2−s2a12s1t2+c3t+c4,ϕt=s2c11+ls11+kt+s21+lc2s11+k+c1p11+k,ψt=p1s21+l2a2s1−a1s22s12t2+1+l1+ka2s1−a1s2+s2c3s11+kt+s2p1c41+l−s1a11+k+s1c3s1p11+k.(iii)
Forn=2,
(38)ur,t=αtlnr+βt−1+ks14r21/1+k,vr,t=ϕtlnr+ψt−1+ls24r21/1+l,where αt,βt,ϕt, and ψt satisfy the linear system of ODEs (35).Example 5.
System(39)ut=1rn−1rn−1u−1/2urr+a1u+s1−a1s1s2uv,vt=1rn−1rn−1v−1/2vrr+a2v+s2−a2s2s1uv,admits of CLBS
(40)η1=urr−12uur2+n−1rur+s1u+4n−4r2u,η2=vrr−12vvr2+n−1rvr+s2v+4n−4r2v.
The exact solutions of this system are listed as follows.(i)
Forn≠2,6,
(41)ur,t=αtr4−n+βtr2−s18n−2r22,vr,t=ϕtr4−n+ψtr2−s28n−2r22,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(42)α′=12a1α−a1s12s2ϕ,β′=12a1β−a1s12s2ψ−2n+8,ϕ′=12a2ϕ−a2s22s1α,ψ′=12a2ψ−a2s22s1β−2n+8.
The solutions of this linear system of PDE are given as follows.(i)
Fora2≠−a1,
(43)αt=c1+c2exp12a1+a2t,βt=2c3a1+a2exp12a1+a2t+24−na1s1+a2s2a1+a2s2t+c4,ϕt=−a2s2c2a1s1exp12a1+a2t+s2c1s1,ψt=−2a2s2c3a1s1a1+a2exp12a1+a2t+24−na1a1s1+a2s2a1s1a1+a2t+a1+a2s2c4+44−ns2−s1a1+a2s1.(ii)
Fora2=−a1,
(44)αt=c1t+c2,βt=a14−ns2−s12s2t2+c3t+c4,ϕt=s2c1s1t+a1c2−2c1a1s1,ψt=a14−ns2−s12s1t2+2n−4s2−s1+a1s2c3a1s1t+2a1c4−2c3s2+44−ns2a1s1.
Forn=2,(45)ur,t=αtr2+βtr2+s1r24lnr−1322,vr,t=ϕtr2+ψtr2+s2r24lnr−1322,where αt,βt,ϕt, and ψt satisfy the linear system of (42).
Forn=6,(46)ur,t=αtlnrr2+βtr2+s132r22,vr,t=ϕtlnrr2+ψtr2+s232r22,where αt,βt,ϕt, and ψt satisfy the linear system of (42).Example 6.
System(47)ut=1rn−1rn−1u−1urr+a1u+b1ulnu+p1uv1+l,vt=1rn−1rn−1vlvrr+a2v+b2v−l+p2v−llnu,admits of CLBS
(48)η1=urr−1uur2+n−1rur,η2=vrr+lvvr2+n−1rvr.
The exact solutions of this system are listed as follows.
Forn≠2,(49)ur,t=expαtr2−n+βt,vr,t=ϕtr2−n+ψt1/1+l,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(50)α′=b1α+p1ϕ,β′=b1β+p1ψ+a1,ϕ′=a21+lϕ+p21+lα,ψ′=a21+lψ+p21+lβ+1+lb2.
The linear system is solvable. However, we do not list the corresponding results here because of the complexity form of the solutions.
Forn=2,(51)ur,t=rαtexpβt,vr,t=ϕtlnr+ψt1/1+l,where αt,βt,ϕt, and ψt satisfy the linear system of ODEs (50).Example 7.
System(52)ut=1rn−1rn−1u−1urr+a1u−s1p1s3ulnu+p1ulnv+s3,vt=1rn−1rn−1v−1vrr+a2v−s3p2s1vlnv+p2vlnu+s1,admits of CLBS
(53)η1=urr−1uur2+n−1rur+s3u,η2=vrr−1vvr2+n−1rvr+s1v.
Forn≠2,(54)ur,t=expαtr2−n+βt−s32nr2,vr,t=expϕtr2−n+ψt−s12nr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(55)α′=−s1p1s3α+p1ϕ,β′=−s1p1s3β+p1ψ+a1,ϕ′=−s3p2s1ϕ+p2α,ψ′=−s3p2s1ψ+p2β+a2.
The solutions of this linear system of ODEs are given as below.(i)
Forp2≠−s12p1/s32,
(56)αt=s3c1s1−s1p1c2s3p2exp−s32p2+s12p1s1s3t,βt=s3a2s1p1+a1s3p2ts12p1+s32p2−s1s3c3s12p1+s32p2exp−s32p2+s12p1s1s3t+c4,ϕt=c1+c2exp−s32p2+s12p1s1s3t,ψt=s1a2s1p1+a1s3p2ts12p1+s32p2+s32c3p2s12p1+s32p2p1exp−s32p2+s12p1s1s3t+s1a2s3−a1s1s12p1+s32p2s3+s1c4s3.(ii)
Forp2=−s12p1/s32,
(57)αt=s3c1s1t+s1c1+s3p2c2p2s1,βt=−p1a1s1−a2s32s3t2+c3t+c4,ϕt=c1t+c2,ψt=−s1p1a1s1−a2s32s32t2+a2s3−a1s1+s1c3s3t+c3−a1s3+s1p1c4s3p1.
Forn=2,(58)ur,t=expβt−s34r2rαt,vr,t=expψt−s14r2rϕt,where αt,βt,ϕt, and ψt satisfy the linear system of ODEs (68).Example 8.
System(59)ut=1rn−1rn−1ukurr+a1u+b1u1−k+p1u1−kvl,vt=1rn−1rn−1vlvrr+a2v+b2v1−l+p2v1−luk,admits of CLBS
(60)η1=urr+k−1uur2−1rur,η2=vrr+l−1vvr2−1rvr.
The solutions are listed as(61)ur,t=αtr2+βt1/k,vr,t=ϕtr2+ψt1/l,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(62)α′=2nk+2kα2+ka1α+kp1ϕ,β′=2nαβ+ka1β+kp1ψ+kb1,ϕ′=2nl+2lϕ2+la2ϕ+lp2α,ψ′=2nϕψ+la2ψ+lp2β+lb2.Example 9.
System(63)ut=1rn−1rn−1u23−n/n−4urr+p1s23−n+s12n−3s1u+p1u3n−10/n−4v23−n/n−4,vt=1rn−1rn−1v23−n/n−4vrr+p2s13−n+s22n−3s2v+p2v3n−10/n−4u23−n/n−4,admits of CLBS
(64)η1=urr−3n−10n−4uur2−3n−3rur+s1u3n−10/n−4−n−42r2u,η2=vrr−3n−10n−4vvr2−3n−3rvr+s1v3n−10/n−4−n−42r2v.
The solutions are listed as below.
Forn≠2,(65)ur,t=n−3s1n−22n−4r2+αtr4−n+βtr6−2n4−n/2n−3,vr,t=n−3s2n−22n−4r2+ϕtr4−n+ψtr6−2n4−n/2n−3,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(66)α′=2n−3p1n−4s1s2α−s1ϕ,β′=2n−4s12+n−3s2p1n−4s1β−n−22n−42n−3α2−2n−3p1n−4ψ,ϕ′=2n−3p2n−4s2s1ϕ−s2α,ψ′=2n−4s22+n−3s1p2n−4s2ψ−n−22n−42n−3ϕ2−2n−3p2n−4β.
Forn=2,(67)ur,t=2r2s1ln2r+αtlnr+βt,vr,t=2r2s2ln2r+ϕtlnr+ψt,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(68)α′=s2p1s1α−p1ϕ,β′=s2p1+2s12s1β−12α2−p1ψ,ϕ′=s1p2s2ϕ−p2α,ψ′=s1p2+2s22s2ψ−12ϕ2−p2β.Example 10.
System(69)ut=1rn−1rn−1ukurr+a1u+b1u1−k+p1u1−kv−2/n+2,vt=1rn−1rn−1v−4/n+2vrr+a2v+b2vn+4/n+2+p2vn+4/n+2uk,admits of CLBS
(70)η1=urr+k−1uur2−1rur,η2=vrr−n+4n+2vvr2−1rvr.
The solutions are listed as(71)ur,t=αtr2+βt1/k,vr,t=ϕtr2+ψt−n+2/2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(72)α′=2nk+2kα2+ka1α+kp1ϕ,β′=2nαβ+ka1β+kp1ψ+kb1,ϕ′=2nϕ2ψ−2a2n+2ϕ−2p2n+2α,ψ′=2nϕψ2−2a2n+2ψ−2p2n+2β−2b2n+2.Example 11.
System(73)ut=1rru−1urr+a1u+p1u2v,vt=1rrv−1vrr+a2v+p2v2u,admits of CLBS
(74)η1=urr−2uur2+srur+2s+1r2u,η2=vrr−2vvr2+srvr+2s+1r2v.
The solutions are listed as below.
Fors≠−3,(75)ur,t=3+sαtr−s+1+βtr2,vr,t=3+sϕtr−s+1+ψtr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(76)α′=s+3αβ−p1ϕ−a1α,β′=−a1β−p1ψ,ϕ′=s+3ϕψ−a2ϕ−p2α,ψ′=−a2ψ−p2β.
Fors=−3,(77)ur,t=1αtr2lnr+βtr2,vr,t=1ϕtr2lnr+ψtr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(78)α′=−p1ϕ−a1α,β′=−a1β−p1ψ−α2,ϕ′=−a2ϕ−p2α,ψ′=−a2ψ−p2β−ϕ2.Example 12.
System(79)ut=1rn−1rn−1u23−n/n−4urr+a1u+p1u2n−7/n−4v3−n/n−4,vt=1rn−1rn−1v23−n/n−4vrr+a2v+p2v2n−7/n−4u3−n/n−4,admits of CLBS
(80)η1=urr−2n−7n−4uur2+n−3rur+n−4r2u,η2=vrr−2n−7n−4vvr2+n−3rvr+n−4r2v.
The solutions are listed as below.
Forn≠2,(81)ur,t=αtr3−n+βtr4−n/n−3,vr,t=ϕtr3−n+ψtr4−n/n−3,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(82)α′=n−22n−3αβ2−n−3a1n−4α−n−3p1n−4ϕ,β′=n−22n−3β3−n−3a1n−4β−n−3p1n−4ψ,ϕ′=n−22n−3ϕψ2−n−3a2n−4ϕ−n−3p2n−4α,ψ′=n−22n−3ψ3−n−3a2n−4ψ−n−3p2n−4β.
Forn=2,(83)ur,t=1αtrlnr+βtr2,vr,t=1ϕtrlnr+ψtr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(84)α′=−α3−12a1+p1α,β′=−α2β−12a1+p1β,ϕ′=−ϕ3−12a2+p2ϕ,ψ′=−ϕ2ψ−12a2+p2ψ.Example 13.
System(85)ut=1rru−1urr+a1u+p1u3/2v,vt=1rrv−1vrr+a2v+p2v3/2u,admits of CLBS
(86)η1=urr−32uur2−1rur+s1r2u,η2=vrr−32vvr2−1rvr+s1r2v.
The solutions are listed as below.
Fors1+2<0,(87)ur,t=s1+2r2αtsin1/2−2s1+2lnr+βtcos1/2−2s1+2lnr2,vr,t=s1+2r2ϕtsin1/2−2s1+2lnr+ψtcos1/2−2s1+2lnr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(88)α′=12α3+12αβ2−12a1α−12p1ϕ,β′=12β3+12α2β−12a1β−12p1ψ,ϕ′=12ϕ3+12ϕψ2−12a2ϕ−12p2α,ψ′=12ψ3+12ϕ2ψ−12a2ψ−12p2β.
Fors1+2>0,(89)ur,t=s1+2r2αtsinh1/22s1+2lnr+βtcosh1/22s1+2lnr2,vr,t=s1+2r2ϕtsinh1/22s1+2lnr+ψtcosh1/22s1+2lnr2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(90)α′=−12α3+12αβ2−12a1α−12p1ϕ,β′=12β3−12α2β−12a1β−12p1ψ,ϕ′=−12ϕ3+12ϕψ2−12a2ϕ−12p2α,ψ′=12ψ3−12ϕ2ψ−12a2ψ−12p2β.
Fors1+2=0,(91)ur,t=1r2αtlnr+βt2,vr,t=1r2ϕtlnr+ψt2,where αt,βt,ϕt, and ψt satisfy the four-dimensional dynamical system
(92)α′=−α3−12a1α−12p1ϕ,β′=−α2β−12a1β−12p1ψ,ϕ′=−ϕ3−12a2ϕ−12p2α,ψ′=−ϕ2ψ−12a2ψ−12p2β.


## 5. Conclusions

The second-order CLBS (1) is used to study nonlinear radially symmetric diffusion system (2). Exact solutions and symmetry reductions are provided due to the compatibility of the governing system (2) and the admitted DC (3) corresponding to the invariant surface condition for CLBS (1). Those solutions extend the known ones such as instantaneous source solutions of the porous medium equation with absorption term. In addition, these results cannot be obtained in the framework within the classical symmetry and the nonclassical symmetry.The topics of exact solutions and symmetry reductions are interesting. Recently, Kumar and collaborators [53–60] have done many superior works about these topics for different kinds of multidimensional evolution equations. The studies of CLBS for multidimensional evolution equations will be involved in our future studies.



---

*Source: 8891750-2021-01-18.xml*