Methods for Constructing Complex Solutions of Nonlinear PDEs Using Simpler Solutions

: This paper describes a number of simple but quite effective methods for constructing exact solutions of nonlinear partial differential equations that involve a relatively small amount of intermediate calculations. The methods employ two main ideas: (i) simple exact solutions can serve to construct more complex solutions of the equations under consideration and (ii) exact solutions of some equations can serve to construct solutions of other, more complex equations. In particular, we propose a method for constructing complex solutions from simple solutions using translation and scaling. We show that in some cases, rather complex solutions can be obtained by adding one or more terms to simpler solutions. There are situations where nonlinear superposition allows us to construct a complex composite solution using similar simple solutions. We also propose a few methods for constructing complex exact solutions to linear and nonlinear PDEs by introducing complex-valued parameters into simpler solutions. The effectiveness of the methods is illustrated by a large number of speciﬁc examples (over 30 in total). These include nonlinear heat equations, reaction–diffusion equations, wave type equations, Klein–Gordon type equations, equations of motion through porous media, hydrodynamic boundary layer equations, equations of motion of a liquid ﬁlm, equations of gas dynamics, Navier–Stokes equations, and some other PDEs. Apart from exact solutions to ‘ordinary’ partial differential equations, we also describe some exact solutions to more complex nonlinear delay PDEs. Along with the unknown function at the current time, u = u ( x , t ) , these equations contain the same function at a past time, w = u ( x , t − τ ) , where τ > 0 is the delay time. Furthermore, we look at nonlinear partial functional-differential equations of the pantograph type, which, in addition to the unknown u = u ( x , t ) , also contain the same functions with dilated or contracted arguments, w = u ( px , qt ) , where p and q are scaling parameters. We propose an efﬁcient approach to construct exact solutions to such functional-differential equations. Some new exact solutions of nonlinear pantograph-type PDEs are presented. The methods and examples in this paper are presented according to the principle “from simple to complex”.


Preliminary Remarks
Exact solutions of nonlinear partial differential equations and methods for their construction are necessary for the development, analysis, and verification of various mathematical models used in natural and engineering sciences, as well as for testing approximate analytical and numerical methods. There are several basic methods for finding exact solutions and constructing reductions of nonlinear partial differential equations: the method of group analysis of differential equations (the method of searching for classical symmetries) [1][2][3][4][5][6], methods for finding for nonclassical symmetries [7][8][9][10], the direct Clarkson-Kruskal method [11][12][13][14], methods for generalized separation of variables [13][14][15], methods for functional separation of variables [14,[16][17][18], the method of differential constraints [13,14,19], the method of truncated Painlevé expansions [13,20,21], and use of conservation laws to obtain exact solutions [22][23][24]. The application of these methods requires considerable special training and, as a rule, is accompanied by time-consuming analysis and a large volume of analytical transformations and intermediate calculations.
This paper describes a number of simple, but quite effective, methods for constructing exact solutions of nonlinear partial differential equations, which do not require much special training and lead to a relatively small amount of intermediate calculations. These methods are based on the following two simple, but very important, ideas: • simple exact solutions can serve as a basis for constructing more complex solutions of the equations under consideration, • exact solutions to some equations can serve as the basis for constructing solutions to other more complex equations.
The effectiveness of the proposed methods is illustrated by a large number of specific examples. Nonlinear heat equations, reaction-diffusion equations, wave type equations, Klein-Gordon type equations, equations of motion in porous media, hydrodynamic boundary layer equations, equations of motion of a liquid film, equations of gas dynamics, Navier-Stokes equations, and some other PDEs are considered. In addition to exact solutions of 'ordinary' partial differential equations, some exact solutions of more complex nonlinear delay PDEs with constant and variable delay and pantograph-type functional-differential equations with partial derivatives are also described.
The methods and examples in the article are presented according to the principle "from simple to complex". For the convenience of a wide audience with different mathematical backgrounds, the authors tried to do their best, wherever possible, to avoid special terminology.

Concept of 'Exact Solution' for Nonlinear PDEs
In this article, the term 'exact solution' for nonlinear partial differential equations will be used in cases where the solution is expressed: (i) in terms of elementary functions, functions included in the equation (this is necessary when the equation contains arbitrary functions), and indefinite or/and definite integrals; (ii) through solutions of ordinary differential equations or systems of such equations.
Combinations of cases (i) and (ii) are also allowed. In case (i), the exact solution can be presented in explicit, implicit, or parametric form.

Some Definitions. Simplest Transformations
We say that a partial differential equation, F(x, t, u, u x , u t , u xx , u xt , u tt , . . .
It is important to note that the free parameter a, which can take values in a certain interval (a 1 , a 2 ), is not included in Equation (1).
Transformations that preserve the form of Equation (1)  A solution u = Φ(x, t) of Equation (1) is called invariant if, under the transformation (2), it transforms into exactly the same solutionū = Φ(x,t).
Further, we will consider only one-parameter transformations of the form x =x + b 1 , t =t + b 2 , u =ū + b 3 (translation); x = c 1x , t = c 2t , u = c 3ū (scaling), and the composition of these transformations. Here b n and c n (n = 1, 2, 3) are constants depending on the free parameter a. Such transformations will be called the simplest transformations.
The following example shows how to determine the invariants of the simplest transformations and the form of the corresponding invariant solutions. Example 1. Consider a transformation consisting of the translation in x and the scaling in t and u: where k and m are some constants. Excluding the parameter a, we find two functionally independent invariants: If the considered equation is invariant under transformation (3), then it admits an invariant solution, which can be represented as I 2 = ϕ(I 1 ) [1,4] or u = t k ϕ(z), where z = x + m ln t. Substituting the resulting expression into the original equation we arrive at an ordinary differential equation for the function ϕ = ϕ(z).

Construction of Complex Solutions from Simpler Solutions. Examples
Simple one-term solutions in the form of a product of functions of different variables are most easily found by the method of separation of variables (the simplest solutions of this type u = Ax α t β are easily determined from the equations under consideration by the method of undefined coefficients). The methods for constructing more complex solutions based on such solutions are described below.
First, we will consider a simple multiplicative separable solution of the special form where k is some constant and ϕ 1 (x) is some function. Such solutions do not change (are invariant) under the scaling transformation Below, in the form of a proposition, we describe a method that allows us to construct more complex solutions based on one-term solutions of the form (5).

Proposition 1. Let the equation
which does not explicitly depend on the x, has a simple solution of the form (5), and does not change under the scaling transformation (6) (i.e., Equation (7) has the same property as the original solution (5)). Then this equation also has a more complex solution of the form where m is an arbitrary constant.
Proof. Equation (7) does not explicitly depend on the spatial variable and is invariant under the translation in x. Consider transformation (3), which is a composition of the translation in x and scaling in t and u (see (6)). Transformation (3) has two functionally independent invariants (4). Therefore, the solution invariant under transformation (3) has the form (8) (see Example 1).
Consider now a simple multiplicative separable solution of the special form where n is some constant and ψ 1 (t) is some function. Such solutions do not change (are invariant) under the scaling transformation A more complex solution than (9) can be obtained by using Proposition 1, redefining the constants and variables in (5)-(8) accordingly. A different, but equivalent method for constructing a more complex solution is described below, which is sometimes more convenient to use in practice.

Proposition 2. Let the equation
which does not explicitly depend on t, has a simple solution of the form (9), and does not change under the scaling transformation (10) (i.e., Equation (11) has the same property as the original solution (9)). Then this equation also has a more complex solution of the form where p is an arbitrary constant. (11) is invariant under the translation in t. Consider a transformation that is a composition of the translation in t and the scaling in x and u (see (10)):

Proof. Equation
where p is an arbitrary constant (p = 0). Transformation (13) has two functionally independent invariants I 1 = y = xe pt and I 2 = e npt u. Therefore, a solution invariant with respect to transformation (13) has the form (12).

Remark 3.
In (12), the function argument y is linear in x. Therefore, solution (12) is easy to differentiate with respect to x. This solution representation should be used for equations which contain partial derivatives with respect to x of a higher order than with respect to t.

Example 2. Consider the Boussinesq equation
which describes the unsteady flow of groundwater in a porous medium with a free surface [45]. Equation (14) has a simple exact solution, which is simultaneously a solution of two types (5) and (9). Let us consider in order both possibilities of constructing more complex solutions based on solution (15). (15) and Equation (14) retain their form under the scaling transformation t = ct, u =ū/c. Therefore, by virtue of Proposition 1 Equation (14) admits a more complex exact solution, where the function ϕ = ϕ(z) satisfies the ordinary differential equation (hereinafter ODE): Note that Equation (16) for k = 0 admits a one-parameter family of solutions in the form of a quadratic polynomial, where C is an arbitrary constant. For C = 0 this solution coincides with the original solution (15). 2 • . Solution (15) and Equation (14) retain their form also under the scaling transformation x = cx, u = c 2ū . Therefore, by virtue of Proposition 2, Equation (14) admits another exact solution where p is an arbitrary constant and the function ψ = ψ(y) is described by the ODE:

Example 3. Consider now the Guderley equation
which is used to describe transonic gas flows [46]. Equation (17) admits a simple exact solution, which is a special case of two types of solutions (5) and (9). Let us consider in order both possibilities of constructing more complex solutions based on solution (18). 1 • . Solution (18) and Equation (17) retain their form under the scaling transformation x = cx, u = c −2ū . Therefore, by virtue of Proposition 1, Equation (17) has a more complex exact solution of the form where the function ϕ = ϕ(z) is described by the second-order ODE: For m = 0, this equation admits a one-parameter family of solutions in the form of a cubic polynomial, where C is an arbitrary constant. For C = 0 this solution coincides with the original solution (18). 2 • . Solution (18) and Equation (17) retain their form also under the scaling transformation y = cȳ, u = c 3ū . Therefore, by virtue of Proposition 2, one can also obtain another more complex exact solution, where p is an arbitrary constant and the function ψ = ψ(z) is described by the ODE: Example 4. In gas dynamics, there is a nonlinear wave equation, which admits a simple exact solution of the form This solution belongs to both classes of solutions (5) and (9). Therefore, based on solution (20), we can construct two more complex solutions described below. 1 • . Solution (20) and Equation (19) are invariant under the scaling transformation t = ct, u = c −2/bū . By virtue of Proposition 1, Equation (19) has a more complex solution of the form where the function ϕ = ϕ(z) is described by the ODE: (20) and Equation (19) are also invariant under the scaling transformation x = cx, u = c 2/bū . Therefore, by virtue of Proposition 2, Equation (19) admits another solution where p is an arbitrary constant and the function ψ = ψ(y) satisfies the second-order nonlinear ODE: Example 5. The system of boundary layer equations on a flat plate by introducing a stream function is reduced to one nonlinear third-order PDE: where ν is the kinematic viscosity of the fluid [47]. Equation (21) has a simple solution of the form which generates two more complex solutions. 1 • . Solution (22) and Equation (21) do not change under the scaling transformation x = ax, u = aū. Therefore, by virtue of Proposition 1, Equation (21) has a more complex solution of the form u = xϕ(z), z = y + k ln x, where k is an arbitrary constant and the function ϕ = ϕ(z) satisfies the ODE: −ϕϕ zz + (ϕ z ) 2 = νϕ zzz .
2 • . Solution (22) and Equation (21) do not change also under the scaling transformation y = aȳ, u =ū/a. Therefore, by virtue of Proposition 2, Equation (21) admits another solution where p is an arbitrary constant and the function ψ = ψ(z) is described by the ODE: Example 6. Consider a fourth-order nonlinear evolution equation describing the change in the film thickness of a heavy viscous liquid moving along a horizontal superhydrophobic surface with a variable surface tension coefficient, where a, b, and c are some constants [48,49]. Equation (23) has a simple solution of the form where the function f = f (t) is described by the first-order ODE with separable variables Solution (24) and Equation (23) are invariant under the scaling transformation x = kx, u = k 2/3ū . Therefore, by virtue of Proposition 2, Equation (23) also has a more complex solution of the form where the function ψ = ψ(y) satisfies the fourth-order nonlinear ODE [48,49]: where p is an arbitrary constant.

Example 7.
Consider a nth-order nonlinear PDE of the form Equation (25) has a simple solution, where the function ϕ = ϕ(x) is described by the ODE: Solution (26) and Equation (25) are invariant under the scaling transformation t = at, u = a 1/(1−s)ū . Therefore, by virtue of Proposition 1, Equation (25) also has a more complex solution of the form where m is an arbitrary constant and the function θ = θ(z) satisfies the ODE: Proposition 3. Let the equation which does not explicitly depend on x and t (and therefore admits the traveling-wave solution [13]) does not change under scaling of the unknown function where c > 0 is an arbitrary constant. Then this equation admits an exact solution (more complicated than the traveling-wave solution) of the form where k, p, and q are arbitrary constants (pq = 0).

Proof. Consider a transformation that is a composition of translations in
x and t and scaling of the unknown function (28): where a > 0 is an arbitrary constant (c = a −k/q ), p and q are some constants (pq = 0). The transformation (30) preserves the form of equation (27) and has two functionally independent invariants I 1 = z = px + qt and I 2 = e −kt u. Therefore, a solution that is invariant with respect to transformation (30), can be represented as (29). Solution of the form (29) is obtained from the invariant solution by applying the scaling transformation (28).

Example 8.
Consider the nonlinear heat-type equation where f = f (ξ) is an arbitrary function. Equation (31) is invariant under the scaling transformation (28). Therefore, by virtue of Proposition 3, this equation has a solution of the form (29), where the function ϕ = ϕ(z) satisfies the nonlinear ODE:

Example 9.
Consider a more complex nonlinear PDE of order n, Equation (32) is invariant under the scaling transformation (28). Therefore, by virtue of Proposition 3, this equation has a solution of the form (29), where the function ϕ = ϕ(z) satisfies the nonlinear ODE: where F(w 1 , w 2 , . . . , w n ) is an arbitrary function.

Generalization to Nonlinear Multidimensional Equations
The above Propositions 1-3 allow obvious generalizations to the case of an arbitrary number of spatial variables.

Example 10.
Consider the nonlinear heat equation with n spatial variables Equation (33) admits a simple multiplicative separable solution, where the function ϕ = ϕ(x 1 , . . . , x n ) satisfies the stationary equation Solution (34) and Equation (33) are invariant under the scaling transformation t = ct, u = c −kū . Therefore, by virtue of Proposition 1, Equation (33) also has a more complex solution of the form where m i are arbitrary constants, and the function θ = θ(z 1 , . . . z n ) satisfies the stationary equation

Generalization to Nonlinear Systems of Coupled Equations
The above Propositions 1-3 can also be used to find exact solutions of systems of coupled PDEs. Example 11. Consider the nonlinear system consisting of two coupled reaction-diffusion where a and b are some constants (b = 0), and f (z) and g(z) are arbitrary functions. System of Equation (35) has a simple solution of the form where the functions ϕ = ϕ(t) and ψ = ψ(t) are described by the system of first-order ODEs: Solution (36) and system of Equation (35) are invariant under the scaling transformation x = cx, u = c 2/bū , v = c 2/bv . Therefore, by virtue of Proposition 2, the system of Equation (35) also has a more complex solution of the form where the functions Φ = Φ(z) and Ψ = Ψ(z) are described by the ODE system: where m is an arbitrary constant.
Example 12. Consider another nonlinear system consisting of two coupled reaction-diffusion equations where a and b are some constants (b = 0), and f (z) and g(z) are arbitrary functions. System of Equation (37) has a simple solution of the form where the functions ϕ = ϕ(x) and ψ = ψ(x) are described by the second-order ODE system Solution (38) and system of Equation (37) are invariant under the scaling transformation t = ct, u = c −1/bū , v = c −1/bv . By virtue of Proposition 1, the system of Equation (37) also has a more complex solution of the form where m is an arbitrary constant, and the functions Φ = Φ(z) and Ψ = Ψ(z) satisfy the system ODE:

Construction of Complex Solutions by Adding Terms to Simpler Solutions
In some cases, simple solutions can be generalized by adding one or more additional terms to them, which leads to more complex solutions with generalized separation of variables [13][14][15]. We demonstrate the possible course of reasoning in such cases using the examples of the Boussinesq Equation (14) and the Guderley Equation (17).

Example 13.
As mentioned earlier, the Boussinesq Equation (14) has a solution with a simple separation of variables (quadratic in x, see (15)), which we write as Let us try to find a more complex solution in the form of the sum whose first term coincides with the solution (39). The second term of formula (40) includes the function ψ(t) and the coefficient k, which must be found. Substituting (40) in (14), after elementary transformations we get Since this equality must hold identically for any x, the functional coefficients for various powers of x in (41) must be zero. Thus, there are two possible cases k = 0 and k = 1/2 (both correspond to the vanishing of the coefficient at x 2k−2 ), which must be considered separately. 1 • . The first case. Substituting k = 0 into (41), to define the functions ϕ = ϕ(t) and ψ = ψ(t), we have the system of ODEs: the general solution of which is determined by the formulas where C 1 and C 2 are arbitrary constants. 2 • . The second case (the Barenblatt-Zeldovich dipole solution [50]). Substituting k = 1/2 into (41), we obtain a system of ODEs for determining the functions ϕ = ϕ(t) and Given the formulas (40), (42), (43), as a result, we obtain two three-parameter generalized separable solutions of Equation (14): where for the sake of greater generality, an arbitrary translation in x is additionally added.
We will look for more complex solutions (with generalized separation of variables) Equation (17) in the form where the functions ϕ(x) and ψ(x) and the constant k = 0 are determined in the subsequent analysis (solution (18) is a particular case of solution (44) for k = 3 and ψ = 0).
It is important to note that binomial solutions of the form (44) are quite often encountered in practice and are the simplest generalized separable solutions of nonlinear PDEs.
Substituting (44) in (17), after rearranging the terms, we come to the relation which contains the power functions y k and y 2k−3 and must be satisfied identically for any y.
Consider two cases: ψ xx = 0 and ψ xx = 0. 1 • . The first case. When ψ xx = 0 we get a binomial equation with separable variables, which can be satisfied if we set The general solution of the autonomous ODE (46) can be represented in the implicit form Moreover, this equation admits a particular solution of the power form ϕ = 1 3a (x + C 1 ) −2 , which leads to a three-parameter exact solution of Equation (17): where C 1 , C 2 , and C 3 are arbitrary constants. 2 • . Second case. To balance the function ψ xx = 0 with second term in equality (45), we must set k = 3/2. As a result, we obtain a binomial equation, which can be satisfied by setting ϕ xx = 0, ψ xx = 9 8 aϕ 2 . These equations are easily integrated and lead to a four-parameter exact solution of Equation (17): where C 1 , C 2 , C 3 , and C 4 are arbitrary constants.

Example 15.
Let us return to the hydrodynamic boundary layer Equation (21). It is easy to verify that this equation admits the self-similar solution [51]: where the function F = F(ξ) satisfies the third-order ODE: −(F z ) 2 = νF zzz . We look for a more general solution of Equation (21) by adding the function ϕ(x) to (49): Simple calculations show that ϕ(x) = a ln x, where a is an arbitrary constant. As a result, we obtain a non-self-similar solution of the boundary layer equation (21) of the form [13]: where the function F = F(ξ) is described by the third-order ODE: −(F z ) 2 − aF zz = νF zzz .

Construction of Compound Solutions (Nonlinear Superposition of Solutions)
In some cases, two similar but different solutions of the considered nonlinear PDE can be combined to obtain a more general composite solution. We demonstrate the possible course of reasoning in such cases by examples of the Goderley equation and the nonlinear diffusion equation with a second-order volume reaction. (47) and (48), it follows that the Guderley Equation (17) has two solutions of the same type u 1 = ϕy 3/2 + ψ and u 2 = ϕy 3 + ψ, which differ from each other by the exponent y. This circumstance suggests an attempt to construct a more general solution of Equation (17), that includes both terms with different exponents at once. In other words, we are looking for a composite solution of the form

Example 16. From expressions
Substitute it in Equation (17). After combining the functional factors for power-functions y 3n/2 (n = 0, 1, 2), we get For this equality to hold for any y, it is necessary to equate the functional factors of y 3n/2 to zero. As a result, we arrive at the system of ODEs: Thus, it is constructively proved that Equation (17) admits a solution of the form (50) (this solution was obtained in [52]).
It can be shown that the system (51) admits the exact solution

Example 17.
Let us now consider a nonlinear diffusion equation with the second-order volume reaction u t = a(uu x ) − bu 2 .
The procedure for constructing a composite solution of this equation will be carried out in two stages: first, we will find two fairly simple solutions, and then, using these solutions, we will construct a composite solution. 1 • . Solutions of exponential form in x. Exact generalized separable solutions of Equation (52) are sought in the form where functions ϕ = ϕ(t) and ψ = ψ(t) and the constant λ are to be determined in the subsequent analysis. Substituting (53) in (52) and collecting similar terms at exponents e nλx (n = 0, 1, 2), we get Since this equality must be satisfied identically for any x, the functional factors of e nλx must be equated to zero. As a result, we come to the differential-algebraic system which allows two solutions where C 1 and C 2 are arbitrary constants. 2 • . Composite solution of exponential form in x. From relations (53) and (54), it follows that Equation (52) has two solutions u 1,2 = ϕe ±λx + ψ. They differ in structure from each other only by the sign of the exponent λ.
This circumstance suggests trying to construct a more general solution of Equation (52), which includes both exponential terms at once. In other words, we are looking for a composite solution of the form Substituting (55) in (52), after elementary transformations we have Equating the functional factors of e nλx (n = 0, ±1) to zero, we arrive at the first-order ODE system (ϕ 1 ) t + 3 2 bϕ 1 ψ = 0, (ϕ 2 ) t + 3 2 bϕ 2 ψ = 0, Thus, it has been proved that Equation (52) admits the solution of the form (55).
By excluding ψ from the first two equations in (56), we obtain the equality (ϕ 1 ) t /ϕ 1 = (ϕ 2 ) t /ϕ 2 . This implies that ϕ 1 = Aϕ(t), ϕ 2 = Bϕ(t), where A and B are arbitrary constants. Therefore, the generalized separable solution (55) is reduced to the form where the functions ϕ = ϕ(t) and ψ = ψ(t) are described by the nonlinear system of two ODEs: By excluding t, this autonomous system is reduced to one ODE, which is homogeneous and therefore can be integrated [53]. Note that the system of equations (58) for AB > 0 admits two simple solutions which define the solution (57) in the form of a product of functions of different arguments.
3 • . Solution of trigonometric type in x. When writing formulas (55) and (57) implicitly, it was assumed that ab > 0. For ab < 0, we have In this case, in solution (57) instead of exponential functions, trigonometric functions appear, i.e., it can be represented in the form where A 1 and B 1 are arbitrary constants. Substituting (59) into Equation (52) and performing calculations similar to those in Item 2 • , we obtain the following nonlinear system of ODEs for the functions ϕ = ϕ(t) and ψ = ψ(t): This system allows for two simple solutions which determine the solution (59) in the form of a product of functions of different arguments.

Linear Partial Differential Equations
In the case of linear partial differential equations, the following proposition can be used to construct more complex solutions from simpler solutions. Proposition 4. Let a linear homogeneous PDE with two independent variables x and t have a one-parameter solution of the form u = ϕ(x, t, c), where c is a parameter that is not included in the original equation. Then the considered equation also has two two-parameter solutions where a and b are arbitrary real constants, Re z and Im z are real and imaginary parts of complex number z.
Proof. The validity of the proposition follows from the linearity of the equation and from the fact that the solution u = ϕ(x, t, c) is also a solution for c = a + ib.
Proposition 4 implies the validity of the following two consequences: Corollary 1. Let a linear homogeneous PDE not depend explicitly on the independent variable t and have a solution u = ϕ(x, t). Then this equation also has two one-parameter families of solutions where a is an arbitrary real constant.

Corollary 2. Let a linear homogeneous PDE not depend explicitly on the independent variable x and have a solution u = ϕ(x, t). Then this equation also has two one-parameter families of solutions
where a is an arbitrary real constant.

Example 18. Consider the linear heat equation
It is easy to verify that this equation admits an exact solution of the exponential form where c is an arbitrary parameter. Using Proposition 4, we obtain two more complicated two-parameter families exact solutions of Equation (62):

Example 19. Consider the linear wave equation
It is easy to verify that Equation (63) admits translation transformations for both independent variables and has the particular solution Making a translation in solution (64) with an imaginary parameter in t and using Corollary 1, we find two more complicated one-parameter families of solutions to Equation (63): Making a translation in solution (64) with an imaginary parameter in x and using Corollary 2, we find two other one-parameter families of solutions: Example 20. Consider the linear heat equation which describes two-dimensional processes with axial symmetry, where x is the radial coordinate. It is easy to verify that Equation (65) admits a translation transformation with respect to the variable t and has the particular solution Making a translation in solution (66) with an imaginary parameter in the variable t and using Corollary 1, we find two more complicated one-parameter families of solutions: t cos ax 2 4(t 2 + a 2 ) + a sin ax 2 4(t 2 + a 2 ) , .

Example 21. Consider the linear wave equation with variable coefficients
This equation admits a translation transformation with respect to t and has the exact solution where C is an arbitrary constant. By making a translation in solution (68) an imaginary parameter in t and using the Corollary 1, we can find two more complicated one-parameter families of solutions by formulas The final form of these solutions is not presented here, due to the cumbersomeness of their recording. The solution u 1 was obtained in [54] and was used to describe the propagation of localized disturbances in one-dimensional shallow water over an inclined bottom. Note, that in [55], another exact solution of the Equation (67) was obtained by integrating the parameter a. ϕ(x, t, c), where c is a real parameter that is not included in the equation. Then, by n-fold differentiation or integration of this solution, one can obtain other exact solutions of the considered equation [56,57].

Corollary 3.
Exact solutions of linear PDEs, which do not explicitly depend on the independent variable t, can be constructed by differentiating or/and integrating with respect to parameters a and b in solutions (61) that are obtained by introducing the complex parameter c = a + ib.

Remark 5.
In [55], by integrating formulas (69) with parameter a, it was obtained a new solution of Equation (67).

Nonlinear Partial Differential Equations
In some cases, it is possible to obtain another solution from one solution, passing from real parameters to complex ones in such a way that the transformed equation and the solution remain real. Let us explain this with a few examples.

Example 22.
Let us return again to Equation (52). It is easy to verify that its trigonometric solution (59) and the system of equations (60) can be obtained from the solution exponential form (57) and system of equations (58), if in the latter we formally set e λx = e iβx = cos(βx) + i sin(βx), e −λx = e −iβx = cos(βx) − i sin(βx), Example 23. Consider the equation where f (w) is an arbitrary function.
It is easy to verify that Equation (71) has the simple multiplicative separable solution (exponential in x): where the parameter λ and the function ψ(t) are to be determined in the subsequent analysis. Substituting (72) in (71), we obtain two solutions of the form (72), where The presence of two solutions of the same type corresponding to ±λ suggests trying them 'combined' and looking for a more general composite solution of the form where A and B are some constants. Substituting (74) in Equation (71), we obtain a solution of the form (74), where A and B are arbitrary constants, and the function ψ = ψ(t) satisfies the nonlinear ODE: Substituting (70) into (74) and (75), we arrive at a new solution containing already trigonometric functions in x, where A 1 and B 1 are arbitrary constants, and the function ψ = ψ(t) is described by a nonlinear ODE: Remark 6. Exact solutions of the nonlinear hyperbolic equation are constructed in the same way.
2 • . Exact solutions of some nonlinear PDEs can be obtained using the proposition below. Proposition 6. Let a nonlinear PDE have an exact solution involving trigonometric functions of the form u = F(x, t, A cos(βx) + B sin(βx), β 2 ), (76) where A, B, and β are free real parameters that are not included in the considered equation. Then this equation also has the exact solution involving hyperbolic functions: whereĀ,B, λ are free real parameters. The converse is also true: if an equation has the exact solution (77), then it also has the exact solution (76).
Solution (77) is obtained from (76) by renaming the parameters β = iλ, A =Ā, Example 24. Consider the fourth-order nonlinear equation to which the stationary Navier-Stokes equations are reduced in the planar case [58]. Equation (78) has the exact solution Therefore, this equation also has the exact solution These solutions and other examples of this kind can be found in [13].

Using Solutions of Simpler Equations for Construct Solutions to Complex Equations
Preliminary remarks. It is often possible to use solutions of simpler equations to construct exact solutions to complex differential equations. In this section, we will illustrate the reasoning in such cases for nonlinear PDEs (see Section 5.1), as well as for more complex nonlinear partial functional differential equations (see Sections 5.2-5.4).

Nonlinear Partial Differential Equations
The following example shows how precise solutions of nonlinear reaction-diffusion equations can be used to generate exact solutions to wave type equations.
Let us now consider a nonlinear wave type equation with a quadratic nonlinearity of the form Equations (79) and (84) differ only in the order of the derivative with respect to t in the left parts of the equations. Since the right-hand sides of these equations involving derivatives with respect to x are the same, it is natural to assume that the power structure of solutions with respect to x of both equations will also be the same, and only the functional factors that depend on t will change for different powers of x.
In other words, we look for exact solutions of wave type PDE (84) in the same form as solutions of reaction-diffusion PDE (79). As a result, we get the following four exact solutions of PDE (84): 1 • . The additive separable solution of the form (80), where the function ψ = ψ(t) is described by the ODE: The multiplicative separable solution of the form (81), where the function ψ = ψ(t) is described by the ODE: The generalized separable solution of the form (82), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the ODEs: The generalized separable solution of the form (83), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the ODEs: The considered example is a good illustration of a rather general fact, which is a consequence of the results of [15] (see also [13,14]) and can be formulated as the following proposition.

Proposition 7.
Let the evolution partial differential equation where F[u] ≡ F(u, u x , . . . , u (n) x ) is the nonlinear differential operator in x, has a generalized separable solution of the form Then, the more complex partial differential equation where L 1 and L 2 are any linear differential operators in t, also has the generalized separable solution of the form (86) with the same functions ϕ k (x) (but with other functions ψ(t)).

Remark 7.
In the equations (85) and (87), the nonlinear operator F can explicitly depend on the variables x and t.
Let us now give an example of constructing an exact solution that cannot be obtained by using Proposition 7.

Example 26.
Consider the nth-order nonlinear PDE: which differs from (32) only in the order of the derivative with respect to t on the left part of the equation. We look for the solution of Equation (88) in the same form as the solution of Equation (32). Substituting (29) in (88), for the function ϕ = ϕ(z) we obtain a nonlinear ODE: k 2 ϕ + 2kqϕ z + q 2 ϕ zz = ϕF(pϕ z /ϕ, p 2 ϕ zz /ϕ, . . . , p n ϕ (n) z /ϕ).

Partial Differential Equations with Delay
In biology, biophysics, biochemistry, chemistry, medicine, control theory, climate model theory, ecology, economics, and many other areas there are nonlinear systems, the rate of change of parameters of which depends not only on the current state of the system at a given time, but also on the state system at some previous time [59]. The differential equations that describe such processes, in addition to the unknown function u = u(x, t) also include the function w = u(x, t − τ), where t is time, τ > 0 is the constant delay. In some cases, we consider situations where the delay depends on the time, τ = τ(t).
The presence of a delay significantly complicates the analysis of such equations. Although nonlinear PDEs with constant delay allow solutions of the traveling wave type u = u(z), where z = x + λt (see, for example, [59][60][61][62]), they do not allow self-similar solutions of the form u = t β ϕ(xt λ ), which often have simpler PDEs without delay.

Example 27. Let us consider a nonlinear reaction-diffusion equation with a constant delay,
Equation (89) is more complicated than the ODE without delay (79) and goes into it at τ = 0. The presence of the delay in (89) does not affect the nonlinear term containing derivatives in x. Therefore, we can assume that the power structure of solutions in x of both equations will be the same, and only the functional factors that depend on t will change.
In other words, we look for exact solutions PDE with delay (89) in the same form, as solutions simpler PDE without delay (79). As a result, we get the following four exact solutions of the nonlinear delay PDE (89): 1 • . The additive separable solution of the form (80), where the function ψ = ψ(t) is described by the linear delay ODE: The multiplicative separable solution of the form (81), where the function ψ = ψ(t) is described by the nonlinear delay ODE: The generalized separable solution of the form (82), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the delay ODEs: The generalized separable solution of the form (83), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the delay ODEs:

Example 28.
More complex than (89), nonlinear PDE with variable delay where τ(t) is an arbitrary function, also admits four exact solutions of the form (80)- (83).

Example 29. The reaction-diffusion equation with logarithmic nonlinearity
admits the exact functional separable solution [15]: where the functions ψ n = ψ n (t) are described by the nonlinear system of ODEs: Let us now consider a more complex nonlinear reaction-diffusion equation with a constant delay, PDE with delay (92) in the special case τ = 0 passes into the simpler PDE without delay (90). For τ = 0 the solution of delay PDE (92), as for Equation (90), is sought in the form (91). As a result, for the functions ψ n = ψ n (t), we obtain the nonlinear system delay ODEs: Example 30. More complex than (92), nonlinear PDE with variable delay where τ(t) is an arbitrary function, also admits a solution with functional separation of variables of the form (91) [64].

Remark 8.
In [71,72], methods for constructing generalized traveling-wave solutions for nonlinear PDEs with delay by using exact solutions more simple PDEs without delay were proposed. This method is suitable for constructing exact solutions in both explicit and implicit form.
Below, with specific examples, it is shown that the exact solutions of nonlinear pantograph-type PDEs can be found by using simpler 'ordinary' PDEs, which do not contain the unknown function with dilated or contracted arguments.

Example 31.
Let us first consider a pantograph-type reaction-diffusion equation with quadratic nonlinearity, Equation (93) is more complicated than the ODE without argument scaling (79) and passes into it at p = q = 1. The presence in (93) of dilation in w does not affect the nonlinear term containing derivatives with respect to x. Therefore, we can assume that the power structure of solutions in x of both equations will be the same, and only the functional factors that depend on t will change (and for the additive separate solution, a factor in x 2 will change).
In other words, we look for exact solutions of pantograph-type PDE (93) in the same form as solutions of 'ordinary' simpler PDE (79). As a result, we get the following four exact pantograph-type solutions (93): 1 • . The additive separable solution is The multiplicative separable solution has the form (81), where the function ψ = ψ(t) is described by the pantograph-type ODE: The generalized separable solution has the form (82), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the pantograph-type ODEs: The generalized separable solution has the form (83), where the functions ψ 1 = ψ 1 (t) and ψ 2 = ψ 2 (t) are described by the pantograph-type ODEs: Example 32. The reaction-diffusion equation with power-law nonlinearity for k = 1 admits a self-similar solution [25]: where the function U = U(z) is described by the nonlinear ODE: Let us now consider a much more complex nonlinear partial functional-differential equation of the pantograph-type where p and q are free parameters (p > 0, q > 0). Parameter values 0 < p < 1 and 0 < q < 1 correspond to equations with proportional delay in two arguments. The functional-differential Equation (96) in the special case p = q = 1 passes into the 'ordinary' partial differential Equation (94). For k = 1 the solution of the pantograph-type PDE (96), as for Equation (94), is sought in the form (95). As a result, for the function U = U(z), we obtain a nonlinear ODE of the pantograph-type [91]: Remark 9. Equation (96) with proportional delays for 0 < p, q < 1 in the special case p = q 1/2 has an exact solution expressed in terms of the solution of the ODE without delay (97) with s = 1; for p < q 1/2 , Equation (96) reduces to the delay ODE with s < 1; and for p > q 1/2 , to the ODE with contracted argument for s > 1. Moreover, a solution of the ODE (96) for p, q > 1 for appropriate values of the parameters p and q can also be expressed in terms of the solution of the ODE with delay (s < 1), without delay (s = 1), and with contracted argument (s > 1).

Example 33.
Let us now consider the reaction-diffusion equation with exponential nonlinearity which for λ = 0 admits the invariant solution [25]: where the function U = U(z) is described by the nonlinear ODE: Let us now consider a much more complex nonlinear functional-differential equation of the pantograph-type where p and q are free parameters (p > 0, q > 0). Partial functional-differential Equation (100) for p = q = 1 passes into the 'ordinary' partial differential Equation (98). For λ = 0 the solution of the pantograph-type Equation (100), as for Equation (98), is sought in the form (99). As a result, for the function U = U(z), we obtain a nonlinear ODE of the pantograph-type [91]: In the special case for p = q 1/2 this equation is a standard ODE (without dilated or contracted arguments).

Example 34.
It is easy to show that the nonlinear Klein-Gordon type equation allows the multiplicative separable solution u(x, t) = ϕ(x)ψ(t).

Approach for Constructing Exact Solutions of Functional Partial Differential Equations
Below, a rather general approach for constructing exact solutions of functional partial differential equations of the pantograph-type is formulated in the form following principle.
The principle of analogy of solutions. Structure of exact solutions to partial functionaldifferential equations of the form F(u, w, u x , u t , u xx , u xt , u tt , . . . ) = 0, w = u(px, qt) (105) often (but not always) is determined by the structure of solutions to simpler partial differential equations: F(u, u, u x , u t , u xx , u xt , u tt , . . . ) = 0.
Equation (106) does not contain the unknown functions with dilated or contracted arguments; it is obtained from (105) by formally replacing w by u.
The solutions discussed in Examples 29-32 were constructed by using the principle of analogy of solutions. Below are two more complex examples.

Example 36.
Consider the pantograph-type reaction-diffusion equation with power-law nonlinearities u t = au xx + bu m w k , w = u(px, qt), that is more complex than (96). Following the principle of analogy of solutions, we set w = u in Equation (107). As a result, we arrive at the equation u t = au xx + bu m+k , which, after renaming m + k to k coincides with Equation (94). Taking into account that the solution of Equation (94) is determined by formula (95), the solution of Equation (107) (by renaming k to m + k in Equation (95)) is sought in the form [91]: u(x, t) = t 1 1−m−k U(z), z = xt −1/2 , k = 1 − m.
As a result, for the function U = U(z) we get the nonlinear pantograph-type ODE:

Remark 10.
In [91], a number of exact solutions of nonlinear pantograph PDEs of diffusion and wave types are obtained, which well confirm the principle of analogy of solutions for pantographtype PDEs.

Brief Conclusions
A number of simple, but quite effective, methods for constructing exact solutions of nonlinear partial differential equations that require a relatively small amount of intermediate calculations are described. These methods are based on the following two main ideas: (i) simple exact solutions can serve as the basis for constructing more complex solutions of the considered equations, (ii) exact solutions of some equations can serve as a basis for constructing solutions of more complex equations. The effectiveness of the proposed methods is illustrated by a large number of specific examples of constructing exact solutions of nonlinear heat equations, reaction-diffusion equations, wave type equations, hydrodynamics equations and some other PDEs. In addition to exact solutions to partial differential equations, some exact solutions to nonlinear delay PDEs and pantograph-type PDEs are also described. The principle of analogy of solutions is formulated, which allows us to constructively find exact solutions to such partial functional-differential equations.