#### 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 initial-boundary value problem of unsteady gas flow through a semi-infinite 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 reduction-to-ODE 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 semi-infinite porous medium. In the investigation of the unsteady flow of gas through a semi-infinite 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 semi-infinite 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 one-dimensional 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 semi-infinite 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 initial-boundary value problem of unsteady flow of gas through a semi-infinite 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 reduced-via-similarity 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.