IJMMSInternational Journal of Mathematics and Mathematical Sciences1687-04250161-1712Hindawi Publishing Corporation31926910.1155/2009/319269319269Research ArticleConstruction Solutions of PDE in Parametric FormVolosovaAlexandra K.VolosovKonstantin AlexandrovichGilanyiAttilaDepartment of Applied MathematicsFaculty of Computer SciencesMoscow State University of Railway Engineering141400 Himki, Moscow regionRussiamiit.ru200912082009200922072008310720092009Copyright © 2009This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The new important property of wide class PDE was found solely by K. A. Volosov. We make an arbitrary replacement of variables. In the case of two independent variables x,t, then it always gives the possibility of expressing all PDE second and more order as AX=b. This is a linear algebraic equations system with regards derivatives to old variables x(ξ,δ), t(ξ,δ) on new variables ξ,δ:(xξ',xδ',tξ',tδ'). This system has the unique solution. In the case of three and more independent variables x,y,t,, then it gives the possibility of expressing PDE second order as AX=b, if we do same compliment proposes. In the present paper, we suggest a new method for constructing closed formulas for exact solutions of PDE, then support on this important new property.

1. Introduction

A new method of construction of exact solutions for partial differential equations (PDE) is proposed in this article. The classical authors in mathematics used change of variables for classification of linear PDE. However, they did not notice an important property of broad class of PDE, which was discovered in Volosov's articles . (Literature reference  is available for review on http://eqworld.ipmnet.ru/) This property gives the possibility of expressing PDE as AX=b. This is a linear algebraic equations system with regards derivatives to old variables on new variables. This equations system has the unique solution. New identity was obtained which follows from solvability conditions of the system. The methods of the above calculations and their consequence are described in this article. Let us consider the:

Zt-(K(Z)Zx)x+F(Z)=0. Let's describe the proposed method based on (1.1). The proposed algorithm works providing that all functions are continuously differentiable functions.

Let's make an arbitrary replacement of variables:

Z(x,t)|x=x(ξ,δ),t=t(ξ,δ)=U(ξ,δ). The inverse replace of variables defines the function Z(x,t)of (1.1) from the function U(ξ,δ):

Z(x,t)=U(ξ,δ)|ξ=ξ(x,t),δ=δ(x,t).

We note that detJ=xξtδ-tξxδ0 is nonezero and not interminable, where

J=(xξtξxδtδ).

An inverse transformation exists, at least locally:

ξ=ξ(x,t),δ=δ(x,t).

The derivatives of the old independent variables on the new variables are determind as follows:

xξ=detJδt,tξ=-detJδx,xδ=-detJξt,tδ=detJξx. Let us introduce the following relation:

K(Z)Zx|x=x(ξ,δ),t=t(ξ,δ)=Y(ξ,δ),K(Z)Zt|x=x(ξ,δ),t=t(ξ,δ)=T(ξ,δ). Using (1.6) and (1.7), we obtain the formulas:

K(U(ξ,δ))(Uξtδ-Uδtξ)=Y(ξ,δ)[xξtδ-tξxδ],K(U(ξ,δ))(-Uξxδ+Uδxξ)=T(ξ,δ)[xξtδ-tξxδ]. Equation (1.1) takes the form

T(ξ,δ)-K(U)((Y/ξ)(t/δ)-(Y/δ)(t/ξ))[xξtδ-tξxδ]+K(U)F(U)=0. Let us rewrite (1.7) in the form

Z(x,t)x=Y(ξ,δ)K(U)|ξ=ξ(x,t),δ=δ(x,t),Z(x,t)t=T(ξ,δ)K(U)|ξ=ξ(x,t),δ=δ(x,t).

As Z is continuously differentiable function, with the necessary of (/t)Zx=(/x)Zt in the variables ξ,δ.

Taking into consideration on (1.6), (1.7), we can write this equality in the form

-xδξ[YK(U)]+xξδ[YK(U)]-tδξ[TK(U)]+tξδ[TK(U)]=0.

System (1.8)–(1.11) will be analyzed in two equationstages. At the first stage, we consider system (1.8)–(1.11) as a nonlinear algebraic equationsystem regarding the derivatives xξ,  xδ,  tξ,  tδ.

Theorem 1.1.

The implicit linear algebraic equationsystem (1.8)–(1.11), with regarding the derivatives xξ,  xδ,  tξ,    tδ, has the unique solution: xξ=Ψ1(ξ,δ),xδ=Ψ2(ξ,δ),tξ=Ψ3(ξ,δ),tδ=Ψ4(ξ,δ), where Ψ1(ξ,δ)=def(TUδYξ2K[TUδYξ2-FKUξ(UδTξ-TδUξ)as-(-TUδYξ2-TTδUξ2+TUδTξUξ-YYδTξUξ+TYδYξUξ+YTδUξYξTUδYξ2)TTδUξ2]TUδYξ2)/P1(ξ,δ),Ψ2(ξ,δ)=defK[-FKUδ[UδTξ-TδUξ]-TTξUδ2+YYδTξUδ+TTδUξUδ-YTδYξUδ+TYδYξUδ-TYδ2Uξ]/P1(ξ,δ),Ψ3(ξ,δ)=defK[-YYξ+FKUξ+TUξ][UδYξ-YδUξ]P1(ξ,δ),Ψ4(ξ,δ)=defK[-YYδ+FKUδ+TUδ][UδYξ-YδUξ]P1(ξ,δ), where P1(ξ,δ)=FK[(TYξ-TξY)Uδ+(YTδ-TYδ)Uξ]+TY[-UδTξ+UξTδ]+Y2[YδTξ-TδYξ]+T2[UδYξ-YδUξ], and moreover detJ=K(U)2(YδUξ-UδYξ)2P1(ξ,δ).

Proof.

Equation (1.11) is linear. Let's divide the first equation (1.8) by Y, divide second equation by T deduct second equation from first equation, and obtain linear equation. Analogously, the linear equation arises from first equation (1.8) and equation (1.9).

We can express any three derivatives xξ,  xδ,  tξ, by means of one derivative. For example, we can express them by means of tδ.

We can substitute it to (1.9) and obtain a linear algebraic equation for tδ(!). This is valid for any PDE of the second order. This is the substance of the newly discovered property of PDE.

We obtained AX=b.

Matrix A has the form

A=(YUδ-YUξTUδ-TUξa21a22a23a2400a33a34000a44), where is a21=-K(U)Yδ+YK(U)Uδ,  a22=K(U)Yξ-YK(U)Uξ,  a23=-K(U)Tδ+TK(U)Uδ,  a24=K(U)Tξ-TK(U)Uξ,a33=-YYδ+F(U)K(U)Uδ+TUδ,  a34=YYξ-F(U)K(U)Uξ-TUξ, and     a44=P1(ξ,δ). Vectors X,  b have the form X=(xξ,xδ,tξ,tδ)T, b=(0,0,0,b4)T,   b4=Ψ4(ξ,δ)P1(ξ,δ). Vector symbol T means conjugation. Eigen values of matrix A have the form λ1=-YYδ+FKUδ+TUδ,λ2=Y2[TδYξ-YδTξ]-T[FK+T][UδYξ-YδUξ]+Y[FK+T][UδTξ-TδUξ],λ3,4=12[M±D],M=KYξ+Y(Uδ-K(U)Uξ),D=4YK(YδUξ-UδYξ)+[KYξ+Y(Uδ-K(U)Uξ)2]. The author proposes the alternative classification for PDE solutions on Eigen values.

Eigen pair can be discovered easily. We are not going to discuss here their interesting properties.

At the second stage, consider the new first-order system (1.12) with the functions x=x(ξ,δ),  t=t(ξ,δ).

It is well known that the solvability of a system of this type is verified by calculating the second mixed derivatives of the functions x=x(ξ,δ) and t=t(ξ,δ) on the arguments ξ and δ:

Theorem 1.2.

(1) one has the new identity (Ψ1/δ-Ψ2/ξ)T(Ψ3/δ-Ψ4/ξ)Y, where the functions Ψi,  i=1,,4 have the form (1.13)–(1.17).

(2) Two solvability conditions (1.21) of system (1.12) have multiply coefficient (or record monomial factor) of arbitrary functions U,Y,TδΨ3-ξΨ4=0, where Ψ3,  Ψ4 are the right-parts in (1.15), (1.16).

Corollary 1.3.

If some free functions U,  Y,  T satisfy the condition (1.23), then systems (1.12) and (1.8)– (1.11) are solvable.

This property (Theorems 1.1 and 1.2) of second-order partial differential equations was not known before.

The specific ways of satisfying the conditions of (1.23) are discussed below.

Example 1.4.

Let's consider the Zel'dovich equation, which is well known in combustion theory ,(Kolmogorov, Petrovsky, Piskunov ): Zt-Z′′xx-Z2(1-Z)=0. Let's consider (1.1) with K(Z)=1, and F(Z)=Z2(1-Z).

Suppose that G(ξ,U) is a function of two variables and w(U),   v(U) are functions of one variable. We seek the functions Y(ξ,δ),  T(ξ,δ) in (1.7) in the form Y(ξ,δ)=G(ξ,U)+h(U),  T(ξ,δ)=w(U)+v(U)G(ξ,U), where U=U(ξ,δ).

In the articles [3, 4] the case h(U)=0 was considered.

Theorem 1.5.

Let oneG(ξ,U)C2[RR], w(U),  h(U),  v(U)C2[R]. Condition (1.23) takes the form [v′′(U)]G3+[2(h-v)v(U)+vh′′-2hv′′-w′′(U)]G2+[(3U2(-1+U)+2hv+2w-4hh)v(U)-2hvh′′+h2v′′+2hw′′(U)]G+U(-2+3U)hv+U(2-3U)w+(1-U)U2vh+2h(-U2+U3+w-hh)v+(U-1)U2w-h2vh′′+h2w′′=0.

We can try to satisfy this (1.25) by making equal the coefficients of the powers of G to zero. We obtain a system of four equations for the two functions w,and  v:

v′′(U)=0,2(h-v)v(U)+vh′′-2hv′′-w′′(U)=0,(3U2(-1+U)+2hv+2w-4hh)v(U)-2hvh′′+h2v′′+2hw′′(U)=0,U(-2+3U)hv+U(2-3U)w+(1-U)U2vh+2h(-U2+U3+w-hh)v+(U-1)U2w-h2vh′′+h2w′′=0.

It turns out that in a number of interesting cases, all the four equations can be solved. Moreover, the function G(ξ,U), as well as U(!), remains arbitrary.

Equation (1.25) can be solved by setting v(U)=3U/2-1/2,   w(U)=3U2(1-U)/2+(-1+3U)h(U)/2.

System (1.12) has the form

xξ=((w+vG)Gξ-(Fv+vw+v2G+G2v+Gw-(w+vG)Gξ)Uξ)P1,xδ=((Gv+w)Gδ-[Fv+v2G+vw+G2v+Gw-(w+vG)GU]Uδ)P1,tξ=(-GGξ+(F+Gv+w-GGU)Uξ)P1,tδ=(F+w+Gv-GGU)UδP1, where P1=Fw+w2+vwG-G2(Gv+w). The Jacobian has the form J=GξUδ/P1.

The exact solutions of the PDE are already constructed, since relations (1.21), (1.23) were satisfied. However, usually it require more detailed formulas. Let's extend reviewed situation.

Let's choose the function G(ξ,U) so that the last two equations (1.32) of the system take the elementary form

t(ξ,δ)=ξ,tξ=1,tδ=0.

Then, we have the following system of two equations for G(ξ,U),:

GU=[(1-U)U2+2(3U-1)G-2h+32Uh-2Gh-2hh]2(G+h),Gξ=3[2U(U-1)-2G-2h]2(2G+2h-2U2)82(G+h). Integrating this system, we obtain:

C1+3ξ2=2(1-U)2U(U-1)-2G-2h-ln(2U(U-1)-2G-2h)+ln(2U2-2G-2h). We still need to analyze the first pair of equations (1.31) for the function x(ξ,δ).

It turns out that, if we can do the variables replacement x(ξ,δ)=χ(U,R),  U=U(ξ,δ),  R=G(ξ,U), then system (1.30)–(1.32) was integrated. After that, we can return to the variables ξ,and  δ, which have the form

x=23(C2+2ln[U-1]-2ln[U])+[(U)+2(h)24(-ln[U-1]+ln[U]+ln(-2G+2(U-1)U-2h)-  ln(-2G+2U2-2h))×((1-3U)2+22(1-3U)h(U)+2(h)2)]/(3(2(1-3U)2+4(1-3U)h(U)+22(h)2))-[22(U-1)((-1+3U)3-32(1-3U)2h(U)+6(-1+3U)(h)2-22(h)3)22]/[3(2G-2(U-1)U+2h)(1-3U+2h)×(2(1-3U)2+4(1-3U)h(U)+22(h)2)]. We have (1.35), (1.36) exact solution equation (1.24) in parametric form.

The authors believe that the solution (1.36) cannot be constructed based on the classic technique of the group analysis. However, if we consider h(U)=0 we have the possibility to return to the original variables x,  t. This equationsimplified solution can be constructed by other methods, for example, by the method of Satsuma-Hirota.

We can came back to the variables x,  t. Let's introduce the notation t=ξ,x(ξ,δ)=X(t,s(x,t)),δ=s(x,t),G=r(t,U(t,s(x,t))), and h(U)=0.

Let's consider (1.24) and the change of variables:

Z(x,t)=U(t,s(x,t)). Let a t>0, and let U(t,s)C2[R+R], have that the relation for the function r(t,U)=r(t,U(t,s(x,t))) takes the form

C1+t=22(1-U)32U(U-1)-6r-23ln(2U(U-1)-2r)+23ln(2U2-2r), (see (1.35), h=0). Also relation (1.36) takes the form

X(t,s)=C1+23[1-U-2r+2(U-1)U+2  ln(2r-2(U-1)U)-2ln(2r-2U2)].

We have a system of two equations for X(t,s),

Xt=12-3U2+3U2(U-1)2r+Utr,Xs=Usr. This is system (1.30)-(1.31) for h=0.

Theorem 1.6.

Suppose that (1.39), (1.40) are hold, and h(U)=0. In this case, (1.24) can be solved. Exact solution equation (1.24) has the form Z(x,t)=U(t,s(x,t)),Z(x,t)=2exp(t/2)-2exp(x/2)2exp(t/2)+exp(x/2)(-2+2t+x2).

Proof.

Let's express from (1.40) ln(2U(U-1)-2r) and substitute in (1.39): r(t,U)=(1-U)(-C1U+2(1+tU)+UX)/(2C1-2t-2X). We will substitute r in (1.40) and raise it to power 3/22 and proceed to exponential functions. Here have put a constant of shift on X (1.40),  C1=0. It always can be restored.

Consider the case h(U)0.

Theorem 1.7.

Suppose h(U)0. Exact solution equation (1.24) has the form (1.35), (1.36).

This is a new real family of solution to (1.24).

Example 1.8.

Consider the Fitz-Hugh-Nagumo-Semenov semilinear parabolic partial differential equation, which is well known in biophysics theory [6, 9]: Zt-Zxx-Z(1-Z2)=0. Here we show one more way of integration of system (1.30)–(1.32). Suppose that G(ξ,U)=ξ (!),   w(U),   v(U) are functions of one variable. We seek the functions   Y(ξ,δ),    T(ξ,δ) in (1.7) in the form Y(ξ,δ)=ξ,T(ξ,δ)=w(U)+v(U)ξ, where U=U(ξ,δ).

In this case, equation analogous (1.25) can be solved by setting v(U)=3U/2,  w(U)=3U(1-U2)/2. To avoid any doubts, let's express (1.30)–(1.32) in an obvious form:

xξ=[2U(1-U2+2ξ)+(-2U2+2U4-2ξ-22ξ2)Uξ]P1,xδ=(-2U2+2U4-2ξ-22ξ2)UδP1,tξ=[-2(2ξ+U(-1+U2-32ξ)Uξ)](3P1),tδ=[2(U-U3+32Uξ)Uδ](3P1).

Suppose G=ξ, and calculate integrals (!).

The denominator in the system of ODE (1.31)–(1.32) can be presented in the form of six multipliers:

P1=Fw+w2+vwG-G2(Gv+w)=U2-2U4+U6+32U2ξ-32U4ξ-2ξ2+6U2ξ2-22ξ3=(U+2+2ξ21/4)(U-2+2ξ21/4)(U-(-1-1+42ξ)2)×(U-(1-1+42ξ)2)(U-(-1+1+42ξ)2)×(U-(1+1+42ξ)2).

Then the exact solution of (1.43) can be written in parametric form as

x(ξ,δ)=12+82ξ(2+8ξ+2+82ξ-1-22ξ-1+42ξ×arctan[2U-1-22ξ-1+42ξ]+-2-8ξ+2+82ξ-1-22ξ+1+42ξ×  arctan[2U-1-22ξ+1+42ξ]),t(ξ,δ)=-23ln[U-2+2ξ21/4]-23ln[U+2+2ξ21/4]+13ln[U+(-1+1+42ξ)2]+13ln[U+(1+1+42ξ)2]+13ln[U-(-1+1+42ξ)2]+13ln[U-(1+1+42ξ)2].

Theorem 1.9.

Suppose that (1.46) is hold. Exact solution equation (1.43) has the form Z(x,t)=U(ξ,δ)|δ=δ(x,t),    ξ=ξ(x,t).

After that, we came back to the variable ξ from (1.46) which takes the form

ξ=-U2+exp(3t)(U2-1)±U2+exp(3t)(1-U2)2(-1+exp(3t)). Also δ is arbitrary parameter.

Let's substitute this expression in the first parity (1.46) x=x(ξ,U) and we will receive the solution in the parametrical form. This is a new real family of solution to (1.43).

Example 1.10.

In the three-dimensional case, consider the semilinear parabolic and the change of variables: Zt-Z′′xx-Z′′yy+f(t,Z)=0,Z(x,y,t)|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=U(ξ,δ,τ). Suppose that the Jacobian J0. Suppose also that there exists (at least locally) the inverse transformation ξ=ξ(x,y,t),δ=δ(x,y,t),τ=τ(x,y,t), and the derivatives an related as τt=yδxξ-xδyξJ,τy=xδtξ-tδxξJ, See . We set Zx|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=Y(ξ,δ,τ),Zy|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=M(ξ,δ,τ),Zt|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=T(ξ,δ,τ). The functions Y,M,and  T are unknown.

Let's supplement the relations written above by the equalities of mixed derivatives

Zx,y=Zy,x,Zx,t=Zt,x,Zy,t=Zt,y in the variables ξ,  δ,  and  τ.

The nonlinear algebraic system of seven equations in the variables ξ,  δ,  τ is similar to system (1.7)–(1.9), which follows from system (1.48)–(1.51) with respect to the nine variable-derivatives, xξ,  xδ,  xτ,  yξ,  yδ,  yτ,  tξ,  tδ,  tτ,which are underdetermined. Hence, there is much arbitrariness in the choice of functions. In this case we used "method C" from . Solutions of (1.51) have the form

yξ=-[T(ξ,δ,τ)tξ-Uξ+Y(ξ,δ,τ)xξ]M,yδ=-[T(ξ,δ,τ)tδ-Uδ+Y(ξ,δ,τ)xδ]M,yτ=-[T(ξ,δ,τ)tτ-Uτ+Y(ξ,δ,τ)xτ]M. By integrating the all resulting relations (1.51) and analogous equation (1.9) we have the theorem.

Theorem 1.11.

Suppose that G(ξ,δ,τ,U) is a twice continuously differentiable function of four variables     and h(t,U),  w(t,U),  r(t,U) are twice continuously differentiable functions of two variables, where U=U(ξ,δ,τ).

Then, the exact solution of (1.48) can by written in parametric form as

T(ξ,δ,τ)=c1w(t,U)+w(t,U)UoUwtdsw2(t,s),Y(ξ,δ,τ)=c2w(t,U),M(ξ,δ,τ)=w(t,U),x(ξ,δ,τ)=r(t,U)+h(t,U)G(ξ,δ,τ,U),y(ξ,δ,τ)=-c1t(ξ,δ,τ)-c2r(t,U)-c2G(ξ,δ,τ,U)h(t,U)+UoUdsw(t,s)+c3,ci0,i=1-3.

The function w(t,U) is determined from the differential

f(t,U)w(t,U)+c1+UoUwt(t,s)dsw2(t,s)-2wU(t,U)=0.

The twice differentiable functions G(ξ,δ,τ,U),r(t,U),h(t,U),t(ξ,δ,τ), as well as the function U(!), remain arbitrary.

If f=f(Z), in (1.48), then h(U),  w(U),  r(U) are functions of one variable. The function w(U) is determined from the first-order ordinary differential equation (known as the Abel equation)

f(U)=w(U)[(c22+1)w(U)-c1].

Example 1.12.

Let's consider the modification equation Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) : Zt-Zxx+F(Z)=0,Z(-,t)=a0,Z(,t)->a1, where F(Z)=2Z3λ2/9±(a0-Z)(a1-Z). In  there were parameters λ=0,  a0=0, a1=1,  a1>a0.

We can choose F(U) for (1.1), then v(U)=±3U/2λ+Uλ,   w(U)=-3F(U)/2 and the relations (1.23), (1.25) hold.

Systems (1.12) has the form: xξ=(w+vG)Gξ-(Fv+vw+v2G+G2v+Gw-(w+vG)Gξ)UξP1,xδ=(Gv+w)Gδ-[Fv+v2G+vw+G2v+Gw-(w+vG)GU]UδP1,tξ=-GGξ+(F+Gv+w-GGU)UξP1,tδ=(F+w+Gv-GGU)UδP1. The functions P1,J have the form P1=Fw+w2+vwG-G2(Gv+w),  J=GξUδ/P1.

The exact solutions of the PDE are already constructed, since relations (1.23), (1.25) were satisfied. However, usually it require more detailed formulas.

Let's assume a0=0,  a1=1,  G(ξ,U)=ξ. The systems (1.57), (1.58) have the form: xξ=[3[-6Uλ2(-9+9U+2U2λ2)+18λ(3+2Uλ2)ξ+Q1Uξ]]λP2,xδ=[3[Q1Uδ]]λP2,tξ=-6[18λξ+Q2Uξ]P2,tδ=-3Q2UδP2, where

Q1=-27Uλ+27U2λ-18U2λ3+24U3λ3+4U4λ5-81ξ-54λ2ξ-36λ3ξ2,Q2=-9Uλ+9U2λ+2U3λ3-27ξ-18Uλ2ξ,P2=(-9+9U+2U2λ2-6λξ)(-9U2λ+9U3λ+2U4λ3-27Uξ-12U2λ2ξ+18λξ2)=(-9+9U+2U2λ2-6λξ)*(3Uk2+k(9U+4U2-12λξ2))×-3Uk2+k(9U+4U2-12λξ2)8λk2,k=9+8λ2.

Denominator P2 has three square multipliers. We can integration second and fourth equations (1.59).

Theorem 1.13.

The exact solution of (1.56) Z(x,t)[a0,a1] can be written in parametric form of systems (1.57), and (1.58), where v(U)=±3U/2λ+Uλ, w(U)=-3F(U)/2, and if a0=0,    a1=1,  and  G(ξ,U)=ξ that one has system (1.59).

Exact solution equation of (1.56) in parametric form has the form:

Z(x,t)=U(ξ,δ)|ξ=ξ(x,t)=[-9E1-3ηkE2±[(9E1+3ηkE2)2+24(E1-2kE2)λ2(-4kλξE2+E1(3+2λξ))](1/2)]/(4λ2(E1-2kE2))|ξ=ξ(x,t), where k=9+8λ2,η=k-3,  and  E1=exp[(3kt+4xλ)/(2η)],  E2=exp[(3t/(2η)].

Proof.

Exact solution of system (1.59) has the form: t=-23ln|-9+9U+2U2λ2-6λξ|--k+33kln|-3Uk2+k(9U+4U2λ2-12λξ)|+3+k3kln|3Uk2+k(9U+4U2λ2-12λξ)|,x=1λln|-9+9U+2U2λ2-6λξ|+-k+3+4λ22λkln|-3Uk2+k(9U+4U2λ2-12λξ)|-2λ(1-k)9+8λ2-3kln|3Uk2+k(9U+4U2λ2-12λξ)|. Hence, from (1.63) we have exp[4xλ+3tk-6+2k](-9+9U+2U2λ2-6λξ)+exp[3t-6+2k][3Uk2-k(9U+4U2λ2-12λξ)]=0.

2. Conclusion

The following new facts complementing PDE theory were discovered in the study of Mr. K. A. Volosov.

The system of four first-order equations (1.8)–(1.11) which is replacing the second-order PDE (1.1) is a system of linear algebraic equations regarding derivatives xξ,  xδ,  tξ,  tδ and has unique solution.

There is a new identity which allows to record monomial factor (1.23) under the solvability conditions. It is sometimes called closure condition.

Monomial factor (1.23) is the closure condition; the solvability condition is one between ness relation for three arbitrary functions U,  Y,  and  T. This property was not noticed earlier. It is useful to consider the respective closure condition together with each differential equation.

Using this approach it is possible to construct new families of exact solutions in the parametric form, the examples of which were presented in this article. These solutions cannot be constructed using the classical method of batch properties analysis.

It is useful to study the Eigen pair's properties for matrix A. The alternative classification of PDE solutions is possible through matrix A eigenvalues.

Below you can find the extract from the comments to the article of Mr. M.V.Karasev, Doctor of Physics and Mathematics, Professor, Head of Chair of Applied Mathematics of Moscow State Institute of Electronics and Mathematics, Laureate of State award of Russian Federation.

The article contains several interesting solutions and whole classes of solutions for nonlinear differential equations with partial secondary derivatives (PDE) important in applications. Sometimes in the theory of PDE the approach is used when by making certain transformation (e.g. by replacement of variables) which is converted in order to bring it to another equation whose solution is already known. Thus, the known solution generates nontrivial solution of initial equation. In Mr. K.A. Volosov's approach it is suggested a priori not fixing the type of replacement of variables, leaving that arbitrary on the first stage. The system derived after replacement of variables contains both required vector-function (sought after solution and its derivatives) and unknown coordinate transformation (Jacobi's matrix of variables replacements). The number of equations in the system is always less than the number of variables which opens possibilities for construction of solution by introduction of different interconnections between the components of sought after vector-function. Limitation for the choice of connections requires consistency of the equation for coordinate replacement function co-ordinate replacement function must be consistent. A condition of consistency is mandatory for the remaining independent components of sought after vector-function.

For example, in classic case, for PDE second order with two independent variables, the sought after vector-function has three components; coordinate replacement contains two more functions; and all of them together subordinate to four differential equations of first-order. It was found that all the components of Jacobi 2×2 matrix of co-ordinate replacement cab be explicitly expressed through the components of the sought after vector-function. On this stage it is not required to introduce any additional connections. Then the received formula for the Jacobi's matrix is considered the system of first-order equations relative to two functions of coordinate replacement, and the condition of this system consistency is noted. It is reduced only to one equation, but contains three components of sought after vector-function. As a result, a big freedom appears to meet the consistency condition. On this stage it is useful to introduce connections, which allows in many cases to reduce the condition of consistency to a standard differential equation or even make it explicitly solvable. After that it is necessary to return to construction of coordinate replacement, that is, to integration of the system of first order, but this time with right part already known. To the extent this integration can be executed explicitly in quadrature, to the same extent it is possible to construct the explicit solution of the initial nonlinear with partial derivatives.

This algorithm is not tied to any group or symmetry attributes. The proposed transformation of nonlinear differential multi variable equations has common nature with the mechanism of integration. The coefficients and type of nonlinearity in the equation solved are not specified, but remain general functions.

This work opens very interesting line in many areas of mathematics and its applications dealing with nonlinear equations with higher partial derivatives.

Acknowledgments

The author is grateful to V. P. Maslov, M. V. Karasev, S. Yu. Dobrokhotov, V. G. Danilov, A. S. Bratus for attention and advice.

VolosovK. A.Gertsenovskie chteniyaProceedings of Gertsen Readings ConferenceApril 2006St. Petersburg, Fla, USAVolosovK. A.Construction of solutions of quasilinear parabolic equations in parametric formDifferential Equations2007434507512MR2358692ZBL1131.35313VolosovK. A. , Doctoral dissertation2007Moscow, RussiaMIEMhttp://eqworld.ipmnet.ruVolosovK. A.Formulas for the exact solutions of quasilinear partial differential equations in implicit formRossiĭskaya Akademiya Nauk. Doklady Akademii Nauk200841811114MR2459484ZeldovichYa. B.BarenblattG. I.LibrovichV. B.MakhviladzeG. M.Mathematical Theory of Combustion and Explosion1980Moscow, RussiaNauka479MR613984MaslovV. P.DanilovV. G.VolosovK. A.Mathematical Modeling of Heat and Mass Transfer Processes (Evolution of Dissipative Structures)1987Moscow, RussiaNauka352MR938746DanilovV. G.MaslovV. P.VolosovK. A.Mathematical Modelling of Heat and Mass Transfer Processes1995348Dordrecht, The NetherlandsKluwer Academic Publishersx+316Mathematics and Its ApplicationsMR1369574KolmogorovA. N.PetrovskiiI. G.PiskunovI. S.A study of the diffusion equation with increase in the quantity of matter, and its application to a biological problemMoscow State University Bulletin A193716125VolosovK. A.Construction of solutions of quasilinear parabolic equations in parametric formSibirskii Journal of Industrial'noi Matematiki2008112(34)2939English translation in Journal of Applied and Industrial Mathematics