# Construction Solutions of PDE in Parametric Form

Authors: Alexandra K. Volosova; Konstantin Alexandrovich Volosov
Journal: International Journal of Mathematics and Mathematical Sciences
 (2009)

---

## Abstract

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 variablesx,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.


---

## Body

## 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 [1–4]. (Literature reference [3] 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:(1.1)Z′t-(K(Z)Z′x)′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:(1.2)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(ξ,δ):(1.3)Z(x,t)=U(ξ,δ)|ξ=ξ(x,t),δ=δ(x,t).We note thatdetJ=x′ξt′δ-t′ξx′δ≠0 is nonezero and not interminable, where(1.4)J=(x′ξt′ξx′δt′δ).An inverse transformation exists, at least locally:(1.5)ξ=ξ(x,t),δ=δ(x,t).The derivatives of the old independent variables on the new variables are determind as follows:(1.6)∂x∂ξ=detJ∂δ∂t,∂t∂ξ=-detJ∂δ∂x,∂x∂δ=-detJ∂ξ∂t,∂t∂δ=detJ∂ξ∂x.
Let us introduce the following relation:(1.7)K(Z)∂Z∂x|x=x(ξ,δ),t=t(ξ,δ)=Y(ξ,δ),K(Z)∂Z∂t|x=x(ξ,δ),t=t(ξ,δ)=T(ξ,δ).
Using (1.6) and (1.7), we obtain the formulas:(1.8)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(1.9)T(ξ,δ)-K(U)((∂Y/∂ξ)(∂t/∂δ)-(∂Y/∂δ)(∂t/∂ξ))[x′ξt′δ-t′ξx′δ]+K(U)F(U)=0.
Let us rewrite (1.7) in the form(1.10)∂Z(x,t)∂x=Y(ξ,δ)K(U)|ξ=ξ(x,t),δ=δ(x,t),∂Z(x,t)∂t=T(ξ,δ)K(U)|ξ=ξ(x,t),δ=δ(x,t).AsZ is continuously differentiable function, with the necessary of (∂/∂t)Z′x=(∂/∂x)Z′t in the variables ξ,δ.Taking into consideration on (1.6), (1.7), we can write this equality in the form(1.11)-∂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:
(1.12)∂x∂ξ=Ψ1(ξ,δ),∂x∂δ=Ψ2(ξ,δ),∂t∂ξ=Ψ3(ξ,δ),∂t∂δ=Ψ4(ξ,δ),
where
(1.13)Ψ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(ξ,δ),(1.14)Ψ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(ξ,δ),(1.15)Ψ3(ξ,δ)=defK[-YY′ξ+FKU′ξ+TU′ξ][U′δY′ξ-Y′δU′ξ]P1(ξ,δ),(1.16)Ψ4(ξ,δ)=defK[-YY′δ+FKU′δ+TU′δ][U′δY′ξ-Y′δU′ξ]P1(ξ,δ),
where
(1.17)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
(1.18)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 derivativesx′ξ,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 obtainedAX=b.
MatrixA has the form
(1.19)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.20)λ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 functionsx=x(ξ,δ) and t=t(ξ,δ) on the arguments ξ and δ:
(1.21)x′′ξδ=x′′δξ,t′′ξδ=t′′δξ.
The central result of this article is as follows [1–4].Theorem 1.2.
(1) one has the new identity(1.22)(∂Ψ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(1.23)∂∂δΨ3-∂∂ξΨ4=0,
where Ψ3,Ψ4 are the right-parts in (1.15), (1.16).Corollary 1.3.
If some free functionsU,Y,T satisfy the condition (1.23), then systems (1.12) and (1.8)– (1.11) are solvable.This property (Theorems1.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 [5],(Kolmogorov, Petrovsky, Piskunov [6–8]):
(1.24)Z′t-Z′′xx-Z2(1-Z)=0.
Let's consider (1.1) with K(Z)=1,  and F(Z)=Z2(1-Z).Suppose thatG(ξ,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[R⊗R],  w(U),h(U),v(U)∈C2[R]. Condition (1.23) takes the form
(1.25)[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,andv:(1.26)v′′(U)=0,(1.27)2(h′-v)v′(U)+vh′′-2hv′′-w′′(U)=0,(1.28)(3U2(-1+U)+2hv+2w-4hh′)v′(U)-2hvh′′+h2v′′+2hw′′(U)=0,(1.29)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 functionG(ξ,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(1.30)∂x∂ξ=((w+vG)G′ξ-(Fv+vw+v2G+G2v′+Gw′-(w+vG)G′ξ)U′ξ)P1,(1.31)∂x∂δ=((Gv+w)G′δ-[Fv+v2G+vw+G2v′+Gw′-(w+vG)G′U]U′δ)P1,(1.32)∂t∂ξ=(-GG′ξ+(F+Gv+w-GG′U)U′ξ)P1,∂t∂δ=(F+w+Gv-GG′U)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 functionG(ξ,U) so that the last two equations (1.32) of the system take the elementary form(1.33)t(ξ,δ)=ξ,t′ξ=1,t′δ=0.Then, we have the following system of two equations forG(ξ,U),:(1.34)G′U=[(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:(1.35)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 replacementx(ξ,δ)=χ(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(1.36)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 variablesx,t. Let's introduce the notation
(1.37)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:(1.38)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(1.39)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(1.40)X(t,s)=C1+23[1-U-2r+2(U-1)U+2ln(2r-2(U-1)U)-2ln(2r-2U2)].We have a system of two equations forX(t,s),(1.41)∂X∂t=12-3U2+3U2(U-1)2r+U′tr,∂X∂s=U′sr.
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)),(1.42)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 caseh(U)≠0.Theorem 1.7.
Supposeh(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]:
(1.43)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:(1.44)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).SupposeG=ξ, and calculate integrals (!).The denominator in the system of ODE (1.31)–(1.32) can be presented in the form of six multipliers:(1.45)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(1.46)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(1.47)ξ=-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:(1.48)Z′t-Z′′xx-Z′′yy+f(t,Z)=0,Z(x,y,t)|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=U(ξ,δ,τ).
Suppose that the Jacobian J≠0. 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
(1.49)∂τ∂t=y′δx′ξ-x′δy′ξJ,∂τ∂y=x′δt′ξ-t′δx′ξJ,…
See [3]. We set
(1.50)∂Z∂x|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=Y(ξ,δ,τ),∂Z∂y|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=M(ξ,δ,τ),∂Z∂t|x=x(ξ,δ,τ),y=y(ξ,δ,τ),t=t(ξ,δ,τ)=T(ξ,δ,τ).
The functions Y,M,andT are unknown.Let's supplement the relations written above by the equalities of mixed derivatives(1.51)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 [3]. Solutions of (1.51) have the form(1.52)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 thatG(ξ,δ,τ,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
(1.53)T(ξ,δ,τ)=c1w(t,U)+w(t,U)∫UoUw′tdsw2(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,ci≠0,i=1-3.
The functionw(t,U) is determined from the differential
(1.54)f(t,U)w(t,U)+c1+∫UoUwt′(t,s)dsw2(t,s)-2w′U(t,U)=0.The twice differentiable functionsG(ξ,δ,τ,U),r(t,U),h(t,U),t(ξ,δ,τ),  as well as the function U(!), remain arbitrary.Iff=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)(1.55)f(U)=w(U)[(c22+1)w′(U)-c1].Example 1.12.
Let's consider the modification  equation  Fisher-Kolmogorov-Petrovsky-Piskunov(FKPP) [8]:
(1.56)Zt′-Zxx′′+F(Z)=0,Z(-∞,t)=a0,Z(∞,t)->a1,
where F(Z)=2Z3λ2/9±(a0-Z)(a1-Z). In [5] there were parameters λ=0,a0=0, a1=1,a1>a0.We can chooseF(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:
(1.57)∂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)G′U]U′δP1,(1.58)∂t∂ξ=-GG′ξ+(F+Gv+w-GG′U)U′ξP1,∂t∂δ=(F+w+Gv-GG′U)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 assumea0=0,a1=1,G(ξ,U)=ξ. The systems (1.57), (1.58) have the form:
(1.59)∂x∂ξ=[3[-6Uλ2(-9+9U+2U2λ2)+18λ(3+2Uλ2)ξ+Q1U′ξ]]λP2,∂x∂δ=[3[Q1U′δ]]λP2,(1.60)∂t∂ξ=-6[18λξ+Q2U′ξ]P2,∂t∂δ=-3Q2U′δP2,
where(1.61)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.DenominatorP2 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,andG(ξ,U)=ξ that one has system (1.59).
Exact solution equation of  (1.56) in parametric form has the form:
(1.62)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,andE1=exp[(3kt+4xλ)/(2η)],E2=exp[(3t/(2η)].Proof.
Exact solution of system (1.59) has the form:
(1.63)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
(1.64)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.(1)
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.(2)
There is a new identity which allows to record monomial factor (1.23) under the solvability conditions. It is sometimes called closure condition.(3)
Monomial factor (1.23) is the closure condition; the solvability condition is one between ness relation for three arbitrary functions U,Y,andT. This property was not noticed earlier. It is useful to consider the respective closure condition together with each differential equation.(4)
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.(5)
It is useful to study the Eigen pair's properties for matrixA. 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 Jacobi2×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.



---

*Source: 319269-2009-09-28.xml*