Research Article  Open Access
A Method for Generating Approximate Similarity Solutions of Nonlinear Partial Differential Equations
Abstract
Standard application of similarity method to find solutions of PDEs mostly results in reduction to ODEs which are not easily integrable in terms of elementary or tabulated functions. Such situations usually demand solving reduced ODEs numerically. However, there are no systematic procedures available to utilize these numerical solutions of reduced ODE to obtain the solution of original PDE. A practical and tractable approach is proposed to deal with such situations and is applied to obtain approximate similarity solutions to different cases of an initialboundary value problem of unsteady gas flow through a semiinfinite porous medium.
1. Introduction
The mathematical modeling of most of the physical processes in fields like diffusion, chemical kinetics, fluid mechanics, wave mechanics, and general transport problems is governed by such nonlinear PDEs whose analytic solutions are hard to find. Therefore, the approach of investigating nonlinear PDEs via reduction to ODEs becomes important and has been quite fruitful in analysis of many physical problems. The reader is referred to [1] for an introduction to different types of such reduction approaches and to have an idea about the advances made in the fields of nonlinear diffusion, fluid mechanics, and wave propagation from the utilization of reductiontoODE approach.
Lie symmetry method provides a powerful general technique for analyzing nonlinear PDEs and reducing them to ODEs. PDEs modeling physical processes naturally inherit symmetries from the underlying physical system. Lie symmetry method takes advantage of these natural symmetries in a PDE and provides similarity variables that lead to the reduction to ODEs; compare, for example, [2–4]. A large amount of the literature about the Lie symmetry theory, its applications, and its extensions is available; see, for example, [1, 5–12].
Although there have been some notable contributions in the applications of similarity method to boundary value problems of PDEs (cf. [1, 12–15]), in general the method has not been utilized in a great deal mainly because of a restriction explained below. The success of the method relies heavily on the success in solving the reduced ODE. Only in few cases, the reduced ODE is integrable in terms of elementary or tabulated functions, but in most cases it is not a simple matter and it is suggested to resort to numerical methods to solve the reduced ODE. This retreat from analytic calculations to numerical computations is not uncommon, but it has not proved practical due to the lack of practical and systematic procedures to utilize the numerical solution of reduced ODE to obtain the solution of original PDE. The aim of this work is to provide one practical way of dealing with such situations.
As a test problem to demonstrate our method, we consider a physical problem arising in the transient flow of gas through a semiinfinite porous medium. In the investigation of the unsteady flow of gas through a semiinfinite porous medium [16–19] initially filled with gas at a uniform pressure at time , the pressure at the outflow face is suddenly reduced from to ( is the case of diffusion into a vacuum) and is thereafter maintained at this lower pressure. The studies in [18, 20] show that the unsteady flow of gas in a porous medium is modeled by a nonlinear partial differential equation.
The nonlinear partial differential equation that describes the unsteady flow of gas through a semiinfinite porous medium has been derived by Muskat [20] in the form where is the pressure within porous medium, is the porosity, is the viscosity, and is the permeability. In the onedimensional medium extending from to , this reduces to with the boundary conditions where the constant is given by the properties of the medium.
In the next section we provide details about the method and its application to obtain approximate similarity solutions of BVPs of the above form.
2. The Method and Its Implementation
The main idea of the method rests on finding an approximation of solution of reduced ODE in the form of a function. This has clear advantages over the option of working only with numerical solution of reduced ODE. In the first place, it can be used to generate approximate solution of original PDE via the similarity variables. Secondly, having the approximate solution of PDE in function form can be more useful than numerical solution as it displays the parameters and variables of the problem, so it requires less processing time and can be used for applications in real time. A detailed description of the method is provided below.
To explain and illustrate our method we consider the IBVP with Without loss of generality we can assume because the change of variable leads to a similar problem with . So the IBVP under study here is with initial and boundary conditions
Step 1 (Reduction of IVBP to a BVP of ODE). The similarity transformations [16, 17] with , reduce the above IBVP (6)(7) to BVP of ODE as follows: The aim of the remaining steps is to find an approximate solution of BVP (9) in function form and then use the similarity transformations (8) to obtain approximate of PDE problem (6)(7).
Step 2. Find numerical solution of BVP (9) and use this as a benchmark for obtaining function form of the solution.
Step 3 (Obtain an initial guess for ). This is a crucial step and in our case we use the lower solution of BVP (9), obtained in [17], as our initial guess. As shown in numerical simulations below, lower solution provides a good initial guess that leads to an accurate enough approximate solution in few iterations. So the initial guess for approximate solution of BVP (9) for all the cases below is taken as follows, see [17, Example 2.1]:
where denotes the error function.
As an alternate, a solution obtained by homotopy analysis method can also be used as initial guess.
Step 4. Improve the initial approximation to get up to the desired level of accuracy.
Here we adopt the following procedure for improving the level of accuracy, starting from initial approximation. The lower solution is of the form
with giving the initial approximation , that is, the lower solution. Numerical simulations suggest that as the values of decrease from by a small decrement, the lower solution moves uniformly towards the numerical solution.
Given a function and a number , we say that lies within band of on an interval if
For a suitable value and numbers , , using the sequences of values
in (11) generates a sequence of curves that uniformly approach numerical solution, finally resulting in the curve
which lies in an band around the graph of numerical solution . The number is chosen according to the desired level of accuracy and the value of is approximated via numerical simulations.
Step 5. Use the similarity variables (8) to get the approximate solution of original IBVP (6)(7).
In the subsequent subsections we illustrate implementation of the above procedure, and we provide simulation results and approximate similarity solutions for different cases of values of parameter .
2.1. Approximate Similarity Solution of the IBVP (6)(7) for
In this case, and the lower solution of the ODE problem (9) becomes which serves as our initial approximation for approximating the solution of ODE problem. Solving the BVP (9) numerically to get and uniformly improving the approximations by simulating the procedure explained above we obtain an approximate solution of the ODE problem as with
The plots of the initial approximation , the numerical solution , the approximate solution , and the are given in Figure 1 where The dotted curves demonstrate some intermediary curves involved in simulations of the uniform approximation process from to , in a manner that as decreases from by a small decrement the curves move uniformly towards numerical solution .
Finally the similarity variables (8) provide the approximate solution of IBVP (6)(7) for as
2.2. Approximate Similarity Solution of the IBVP (6)(7) for
In this case, and hence the initial approximation determined by the lower solution of the ODE problem (9) becomes Proceeding as in Section 2.1 we obtain an approximate solution of the ODE problem, for the case , as with The plots of the initial approximation , the numerical solution , the approximate solution , and the , for the case , are given in Figure 2.
As in Section 2.1 the dotted curves demonstrate some intermediary curves involved in simulations of the uniform approximation process from to , in a manner that as decreases from by a small decrement the curves move uniformly towards numerical solution .
The approximate solution of IBVP (6)(7) for can be found using the similarity variables (8) and is given by
2.3. Approximate Similarity Solution of the IBVP (6)(7) for
Using gives the lower solution of the ODE problem (9) as which serves as our initial approximation for approximating the solution of ODE problem for this case. As above, solving the BVP (9) numerically to get and uniformly improving the approximations by simulating the procedure explained above, we obtain an approximate solution of the ODE problem as with
The plots of the initial approximation , the numerical solution , the approximate solution , and the , for the case , are given in Figure 3. As in the above sections, the dotted curves demonstrate some intermediary curves involved in simulations of the uniform approximation process from to , in a manner that as decreases from by a small decrement the curves move uniformly towards numerical solution .
Using the similarity variables, the approximate solution of IBVP (6)(7) for is found as
2.4. Discussion
The plots of the numerical solution , the lower solution , and the approximate solution for the ODE problem (9) are given in Figures 1, 2, and 3, respectively, for , , and . It can be seen that the lower solution provides a better initial approximation for the values of closer to . Yet a comparison of the plots of initial approximation and the accuracy benchmark curve , in all cases, clearly emphasizes the need of improving the initial approximation to get to get accurate enough analytic approximation of the benchmark curve . The plots in Figures 1, 2, and 3 clearly show that and curves overlap in each case and the corresponding plots of the error demonstrate that, in each case, provides an analytic approximation of the benchmark numerical solution curve . In all cases, the solution decreases with the increase in but descends faster for the case as compared to case .
A comparison of the approximate solution of BVP (9) obtained by our method with solutions obtained by other methods [21–23], for , is presented in Table 1.

2.5. Comments on the Conservation Laws
The conservation laws of can be found using the direct multiplier method [24, 25]. Since the procedure is standard, only the results are stated. Implementing the direct multiplier procedure for multipliers of the form generates the determining equations which yield the solution . Hence there are two linearly independent conservation laws arising from the multipliers and . Next, elementary direct calculations [26] give the corresponding conservation laws as listed below.
For , the conservation law is and, for , the conservation law is where For a given conservation law of (28), if the spatial flux vanishes on the boundary and of the semiinfinite medium, then integration from to provides conserved quantity of the boundary value problem.
As an example, the conservation law (30) is applied to derive conserved quantity of the BVP The conservation law (30) can be written as Integrating with respect to from to and using the boundary conditions imply that which gives the time independent conserved quantity of the BVP.
In general, the boundary conditions will determine which conservation law is to be applied to obtain conserved quantities of the BVP of (28).
3. Conclusion
We present a practical way of obtaining approximate solution, in function form, for the class of PDEs where the PDE can be reduced to an ODE through similarity variables but the reduced ODE is not easily integrable in terms of elementary or tabulated functions. The idea presented here to get approximate solution of PDE, that is, approximating the surface , practically involves approximating a curve which is a tractable problem in comparison to increasingly complex and intractable problem of approximating the surface itself. A combination of simulations, initial approximation, and numerical solution of reduced ODE is utilized to obtain approximation of solution curve which readily generates, via similarity variables, the approximate solution surface of the PDE. This makes it a promising approach especially when reasonably accurate initial approximation of the solution of ODE can be obtained, as was the case here in the form of lower solution. The approach is applied to obtain approximate solutions for some cases of an initialboundary value problem of unsteady flow of gas through a semiinfinite porous medium. The approach can be adapted for obtaining approximate analytic solutions for the class of PDEs where the PDE can be reduced to an ODE through similarity variables. For instance, the approach can be directly applied to all the reducedviasimilarity BVPs of ODEs in [27]. For further application, the approach can be extended to obtain approximate solutions where the reduction is a system of ODEs like the reduced flow problems in [28, 29].
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 W. F. Ames, Nonlinear Partial Differential Equations in Engineering, vol. I, Academic Press, New York, NY, USA, 1965.
 H. Azad, M. T. Mustafa, and M. Ziad, “Group classification, optimal system and optimal reductions of a class of Klein Gordon equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 5, pp. 1132–1147, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bokhari, M. T. Mustafa, and F. D. Zaman, “An exact solution of a quasilinear Fisher equation in cylindrical coordinates,” Nonlinear Analysis, vol. 69, no. 12, pp. 4803–4805, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 M. T. Mustafa and K. Masood, “Symmetry solutions of a nonlinear elastic wave equation with thirdorder anharmonic corrections,” Applied Mathematics and Mechanics, vol. 30, no. 8, pp. 1017–1026, 2009. View at: Publisher Site  Google Scholar
 G. W. Bluman and J. D. Cole, Similarity Methods for Differential Equations, Springer, New York, NY, USA, 1974. View at: MathSciNet
 G. W. Bluman and S. Kumei, Symmetries and Differential Equations, Springer, New York, NY, USA, 1989. View at: Publisher Site  MathSciNet
 A. G. Hansen, Similarity Analyses of Boundary Value Problems in Engineering, Prentice Hall, Englewood Cliffs, NJ, USA, 1964. View at: MathSciNet
 P. E. Hydon, Symmetry Methods for Differential Equations, Cambridge University Press, Cambridge, UK, 2000. View at: Publisher Site  MathSciNet
 N. H. Ibragimov, Elementary Lie Group Analysis and Ordinary Differential Equations, John Wiley & Sons, Chichester, UK, 1999. View at: MathSciNet
 W. Miller, Symmetry and Separation of Variables, AddisonWesley, Reading, Mass, USA, 1977. View at: MathSciNet
 P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, NY, USA, 1986. View at: Publisher Site  MathSciNet
 L. V. Ovsiannikov, Group Analysis of Differential Equations, Academic Press, New York, NY, USA, 1982. View at: MathSciNet
 H. Azad, M. T. Mustafa, and A. F. M. Arif, “Analytic solutions of initialboundaryvalue problems of transient conduction using symmetries,” Applied Mathematics and Computation, vol. 215, no. 12, pp. 4132–4140, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. W. Bluman, “Applications of the general similarity solution of the heat equation to boundaryvalue problems,” Quarterly of Applied Mathematics, vol. 31, pp. 403–415, 1974. View at: Google Scholar  MathSciNet
 B. J. Cantwell, Introduction to Symmetry Analysis, Cambridge University Press, Cambridge, UK, 2002. View at: MathSciNet
 R. P. Agarwal and D. O'Regan, “Infinite interval problems modeling the flow of a gas through a semiinfinite porous medium,” Studies in Applied Mathematics, vol. 108, no. 3, pp. 245–257, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. P. Agarwal and D. O'Regan, “Nonlinear boundary value problems on the semiinfinite interval: an upper and lower solution approach,” Mathematika, vol. 49, no. 12, pp. 129–140, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 R. E. Kidder, “Unsteady flow of gas through a semiinfinite porous medium,” Journal of Applied Mechanics, vol. 27, pp. 329–332, 1957. View at: Google Scholar  MathSciNet
 T. Y. Na, Computational Methods in Engineering Boundary Value Problems, vol. 145 of Mathematics in Science and Engineering, Academic Press, New York, NY. USA, 1979. View at: MathSciNet
 M. Muskat, Flow of Homogeneous Fluids, J. Edwards, Ann Arbor, Mich, USA, 1946.
 Y. Khan, N. Faraz, and A. Yildirim, “Series solution for unsteady gas equation via Mldm—Pade Technique,” World Applied Sciences Journal, vol. 10, pp. 1452–1456, 2010. View at: Google Scholar
 M. A. Noor and S. T. MohyudDin, “Variational iteration method for unsteady flow of gas through a porous medium using He's polynomials and Pade approximants,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2182–2189, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. Wazwaz, “The modified decomposition method applied to unsteady flow of gas through a porous medium,” Applied Mathematics and Computation, vol. 118, no. 23, pp. 123–132, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. C. Anco and G. Bluman, “Direct constructi on method for conservation laws of partial di ff erential equations part I: examples of c onservation law classifications,” European Journal of Applied Mathematics, vol. 13, no. 5, pp. 545–566, 2002. View at: Google Scholar
 S. C. Anco and G. Bluman, “Direct construction method for conservation laws of partial differential equations part II: general treatment,” European Journal of Applied Mathematics, vol. 13, no. 5, pp. 567–585, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 T. Wolf, “A comparison of four approaches to the calculation of conservation laws,” European Journal of Applied Mathematics, vol. 13, no. 2, pp. 129–152, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Dresner, Similarity Solutions of Nonl inear Partial Differential Equations, vol. 88 of Research Notes in Mathematics Series, Longman, 1983. View at: MathSciNet
 M. M. Rashidi, T. Hayat, E. Erfani, S. A. M. Pour, and A. A. Hendi, “Simultaneous effects of partial slip and thermaldiffusion and diffusionthermo on steady MHD convective flow due to a rotating disk,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4303–4317, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M. M. Rashidi, S. A. Mohimanian Pour, T. Hayat, and S. Obaidat, “Analytic approximate solutions for steady flow over a rotating disk in porous medium with heat transfer by homotopy analysis method,” Computers & Fluids, vol. 54, no. 1, pp. 1–9, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Mazhar Iqbal et al. This 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.