Free Boundary Formulation for BVPs on a SemiInfinite Interval and NonIterative Transformation Methods
Abstract
This paper is concerned with two examples on the application of the free boundary formulation to BVPs on a semiinfinite interval. In both cases we are able to provide the exact solution of both the BVP and its free boundary formulation. Therefore, these problems can be used as benchmarks for the numerical methods applied to BVPs on a semiinfinite interval and to free BVPs. Moreover, we emphasize how for two classes of free BVPs, we can define noniterative initial value methods, whereas BVPs are usually solved iteratively. These noniterative methods can be deduced within Lie’s group invariance theory. Then, we show how to apply the noniterative methods to the two introduced free boundary formulations in order to obtain meaningful numerical results. Finally, we indicate several problems from the literature where our noniterative transformation methods can be applied.
Key Words: BVPs on a semiinfinite interval, free boundary formulation, free boundary problems, noniterative numerical method.
AMS Subject Classifications: 65L10, 34B15, 65L08.
1 Introduction
Usually when dealing for the first time with a boundary value problem (BVP) defined on a semiinfinite interval the applied scientist does not know the exact or even an approximate solution. As a consequence, he often is tempted to try for a numerical solution of the problem. Therefore, along the years several approaches have been proposed in order to solve BVPs defined on a semiinfinite interval numerically.
The oldest and simplest approach is to replace the original problem by one defined on a finite interval, where a finite value, the socalled truncated boundary, is used instead of infinity. This approach was used, for instance, by Horwarth [1] and by Goldstein [2, p. 136] to get the tabulated numerical solution of the Blasius problem [3]. However, to get an accurate solution a comparison of numerical results obtained for several values of the truncated boundary is necessary as suggested by Fox [4, p. 92] or by Collatz [5, pp. 150151]. Moreover, in some cases accurate solutions can be found only by using very large values of the truncated boundary. This is, for instance, the case for the branches of the von Karman swirling flows where values of truncated boundaries up to several hundreds were used by Lentini and Keller [6].
The limitation of the above classical approach has lead de Hoog and Weiss [7], Lentini and Keller [8] and Markowich [9, 10] to develop a theory for defining the asymptotic boundary conditions to be imposed at a given value of the truncated boundary. Those asymptotic boundary conditions have to be derived by a preliminary asymptotic analysis involving the Jacobian matrix of the righthand side of the governing equations evaluated at infinity. The main idea of this asymptotic boundary conditions approach is to project the solution into the manifold of bounded solutions. By using the same value of the truncated boundary, a more accurate numerical solution can be found by this approach than that obtained by the classical approach, because in the first case the imposed boundary conditions are obtained from the asymptotic behaviour of the solution. However, we should note that this approach is not straightforward, see the remarks in [11], and that for nonlinear problems highly nonlinear asymptotic boundary conditions usually result. Asymptotic boundary conditions have been applied successfully to the numerical approximation of the socalled “connecting orbits” problems of dynamical systems, see Beyn [12, 13, 14]. Those problems are of interest, not only in connection with dynamical systems, but also in the study of travelling wave solutions of partial differential equations of parabolic and hyperbolic type as shown by Beyn [13], Friedman and Doedel [15], Bai et al. [16], and Liu et al. [17].
A different approach, for the numerical solution of BVPs defined on a semiinfinite, is to consider a free boundary formulation of the given problem, where the unknown free boundary can be identified with a truncated boundary. In this approach the free boundary is unknown and has to be found as part of the solution. This free boundary approach overcomes the need for a priori definition of the truncated boundary. Free BVPs represent a numerical challenge because they are always nonlinear as pointed out first by Landau [18]. However, a free boundary formulation has been successfully applied to several problems in the applied sciences: namely, the Blasius problem [19], a twodimensional stagnation point flow [20], the FalknerSkan model [21, 22, 23], and the model describing a fluid flowing around a slender parabola of revolution [24] in boundary layer theory, the computation of a twodimensional homoclinic connecting orbit [25], and a problem related to the deflection of a semiinfinite pile embedded in soft soil [26]. The last problem is of interest in foundation engineering, for instance, in the design of drilling rigs above the ocean floor, see Lentini and Keller [8] and the references quoted therein.
A different way to avoid the definition of a truncated boundary is to apply coordinate transforms. Coordinate transforms have been applied successfully to the numerical solution of ordinary and partial differential equations on unbounded domains, see Grosch and Orszag [27], Koleva [28] or Fazio and Jannelli [29].
This paper is concerned with two examples on the application of the free boundary formulation of BVPs on a semiinfinite interval. In both cases we are able to provide the exact solution of both the BVP and its free boundary formulation. Therefore, these problems can be used as benchmarks for the numerical methods applied to BVPs on a semiinfinite interval and to free BVPs. In this context, sometimes, it is possible to solve a given free BVP noniteratively, see the survey by Fazio [30], whereas BVPs are usually solved iteratively. Here, for two classes of free BVPs, we define noniterative initial value methods which are referred in literature as noniterative transformation methods (ITMs). Indeed, nonITMs can be defined within Lie’s group invariance theory. For the group invariance theory, the interested reader is referred to Bluman and Cole [31], Bluman and Kumei [32], Barenblatt [33], or Dresner [34].
Let us remark here that the first application of a noniterative initial value method for the numerical solution of Blasius problem of boundary layer theory was given by Töpfer [35]. The algorithm devised by Töpfer for Blasius problem was redefined, and extended to a class of problems, within group invariance theory by Klamkin [36]. The relationship between the invariance of the Blasius problem with respect to a linear group of transformations, the scaling group, and the applicability of a noniterative initial value method was point out by Na [37, 38]. Moreover, in the same papers Na considered BVPs on finite intervals and the invariance with respect to a nonlinear group of transformations: the spiral group. Fazio and Evans, in [39], showed how to apply the scaling group to solve noniteratively free BVPs. The translation group of transformation was first used for the noniterative solution of free BVPs by Fazio [40], see also the recent contribution by Fazio and Iacono [41]. Fazio [42], using the scaling group, defined a noniterative initial value method for the numerical solution of a free boundary value problem governed by a system of first order differential equations. For this author’s knowledge no other Lie’s group of point transformations has been applied to define a noniterative initial value method.
In the past, the main drawback of nonITMs was that they were considered not widely applicable: see the critical considerations by Fox, Erickson and Fan [43], Meyer [44, pp. 9798], Na [45, p. 137] or Sachdev [46, p. 218]. In fact, the simplest way in order to verify if a nonITM is applicable to a particular problem is to use an inspectional analysis as shown by Seshadri and Na [47, pp. 157168], cf. also the discussion on inspectional analysis by Birkhoff [48, pp. 99103].
If we consider any possible extension of nonITMs, then we have to recall the extension of scaling invariance involving physical parameters by Na [49], see also Na [45, Chapters 89]. Moreover, as shown by Fazio [50, 51, 24], within scaling invariance theory, it is possible to define an iterative extension of our approach that can be applied to the most general class of free BVPs.
2 Free boundary formulation
In this paper we provide two examples that support the idea that we can deal with BVPs defined on a semiinfinite interval via their free boundary formulation. In particular, in order to explain the main idea behind our free boundary formulation, we consider the class of BVPs
where is a given function of its arguments, and and are given constants. If we can assume that the first derivative of goes monotonically to zero at infinity, then we replace the problem (2) with its free boundary formulation
where is an unknown free boundary and is a parameter.
We have to remark here that monotonic properties of the solution, its first and second derivative have been demonstrated by Countyman and Kannan [52], for the class of problems in (2) where depends exclusively on .
The following theorem provides, under suitable smoothness conditions, the order of convergence (and the uniform convergence) of the solution of (2) to the solution of (2).
Theorem 1
Suppose and are continuous functions with respect to (and also with respect to in the related free boundary domain ) and that at least in a nonempty interval including , then
where is the maximum norm on and is a positive constant independent on .
The proof of this Theorem can be obtained along the lines of the proof for the convergence Theorem stated in [24] for a free boundary formulation for a class of problems governed by a third order differential equation.
The free boundary formulation allows us to embed a BVP in (2) into a class of problems involving the control parameter . When we solve the free boundary formulation (2) numerically, we can fix a very small value of and apply a grid refinement to verify whether the numerical results agree within a prefixed number of significant digits. Also, it is possible to fix a step size and let goes to zero and verify whether together with . Usually, it suffices to take , , , , , , and compare the obtained numerical results. Let us remark here that sometimes it is possible to solve the free boundary formulation noniteratively, see the survey by Fazio [30], whereas the BVP (2) is usually solved iteratively.
3 Two examples for the free boundary formulation
As a first example we consider the linear problem
where is a positive constant. The solution of (3) is easily found to be
(4) 
so that the missing initial condition is equal to , that is . Figure 1 shows the solution (4) of the BVP (3) for different values of . The bigger is the value of , the harder is to solve the BVP numerically. In fact, for large values of the solution has a boundary layer near .
Let us consider now the free boundary formulation for (3)
with . The solution of (3) is given by
(6) 
Therefore, we can easily verify that as goes to zero the solution of the free boundary formulation (3) converges to the solution of the original problem (3) and the free boundary goes to infinity. Moreover, we realize that the obtained approximation becomes the more accurate the more is near zero, see figure 2.
Let us remark here that the same exact solutions (4) and (6) are still valid if we replace the governing differential equation, in the BVP (3) and its free boundary formulation (3) with the nonautonomous one
(7) 
where we substitute in the free boundary case.
Replacing a linear problem with a nonlinear one can be justified, from a numerical viewpoint, only by considering that in this way we overcome the singularity related to the boundary condition prescribed at infinity. Of course, when the original problem is a nonlinear one a free boundary formulation for it can be really convenient to solve numerically.
As a second example we consider the nonlinear problem
where, again, is a positive constant. The solution of (3) is given by
(9) 
and, again, . Figure 3 shows the solution (9) of the BVP (3) for different values of . Again, for large values of the solution has a boundary layer near .
It can be easily verified that, for instance by comparing figure (1) with figure (3), for the same value of the parameter , the BVP (3) is more challenging than the BVP (3).
Let us consider now the free boundary formulation for (3)
with . The positive solution of (3) is given by
(11) 
where . Also in this case, as goes to zero the solution of the free boundary formulation (3) converges to the solution of the original problem (3) and the free boundary goes to infinity. Moreover, also in this case the obtained approximation becomes the more accurate the more is close to zero, figure 4.
4 NonITMs for the free boundary formulation
As we mentioned in the introduction numerical TMs can be defined within Lie’s group invariance theory.
4.1 Spiral and translation invariance
Let us define a nonITM for the class of free BVPs
(12)  
where is an arbitrary function of its arguments, , , and are arbitrary constants, and is the unknown free boundary. The governing differential equation and the two free boundary conditions are left invariant by the spiral group
(13) 
where is the group parameter.
Let us remark here that by setting we can recover from (4.1) the class of problems invariant with respect to the translation group of transformation defined by (13) for .
Using the transformation (13) we can define the following noniterative algorithm for the numerical solution of (4.1):

Input .

Solve backwards from to (to be defined below) the following auxiliary IVP
using an event locator in order to find such that
(15) 
Compute the group parameter
(16) the free boundary
(17) and the missing initial condition
(18) 
Compute the transformed solution
(19) where, of course, .
We define now a simple event locator which is suitable to be applied with the nonITM for the numerical solution of (4.1). Let us consider the case , the case can be treated in a similar way. We integrate the auxiliary IVP (2) until we get at a mesh point where , and we compute
(20) 
Then, we repeat the last step with the smaller step size given by
(21) 
In defining the last step size in equation (20)(21), we apply a first order (linear) Taylor formula at where we have replaced the first derivative with a backward finite difference approximation.
Let us notice here that this nonITM generalize the one proposed by Fazio and Iacono [41] for the numerical solution of free boundary problems with governing equations invariant with respect to a translation group.
4.2 Scaling invariance
Let us define a nonITM for the class of free BVPs
(22)  
where is an arbitrary function of its arguments, , and are arbitrary constants, and is the unknown free boundary. The governing differential equation and the boundary condition at are invariant with respect to the scaling group:
(23) 
where is the group parameter. Using the invariance properties, we can define the following noniterative algorithm for the numerical solution of (4.2):

Input , and .

Solve the auxiliary IVP

Compute by
(25) and the free boundary
(26) 
Rescale the numerical solution to get and according to (23). In particular, we find
If the computed value of is bigger than , then set a lager value for and repeat the computation. It is evident that there is no need to rewrite the given free BVP in standard form, as suggested by Ascher and Russell [53], because we can choose at our convenience and therefore we usually set .
The algorithm presented above is an original variant of the nonITM defined in [39], where the boundary conditions in (4.2) were replaced by
(28) 
where , and are arbitrary constants and we had to integrate backwards on because these free boundary values are invariant, but the boundary condition at zero is not invariant, under the scaling group (23).
5 Numerical tests
In this section we report on the application of the two nonITMs defined in the previous section. To this end we solve numerically the free boundary formulations defined in section 3.
5.1 Translation invariance
As an example we consider here the noniterative numerical solution of the free boundary formulation (3) with , . With a simple change of variables, the problem (3) clearly belongs to (4.1) with . In table 1 we list sample numerical results. We used a uniform grid, with and , and applied the classical fourth order RungeKutta method (RK4) [54, 55], the sixth order RungeKutta method (RK6), and the eighth order RungeKutta method (RK8) as reported by Butcher in [56] on page 178 and 180, respectively.
Method  

RK4  13.  8152456  2.  6660  0.  999734403 
RK6  13.  8152449  2.  6659  0.  999734408 
RK8  13.  8152449  2.  6659  0.  999734408 
Since the numerical results obtained by RK6 and RK8 are very close, we can infer that in order to improve the numerical accuracy we need to consider a grid refinement.
In table 2 we list the numerical results obtained by RK6 and a grid refinement with the reported step sizes.
0.1  13.  8149  6.  48  0.  999353 
0.05  13.  8152  2.  67  0.  999734 
0.025  13.  8154  7.  37  0.  999927 
0.0125  13.  8155  1.  42  0.  999987 
0.00625  13.  8155  4.  88  0.  999996 
0.003125  13.  8155  1.  71  1.  000001 
The reported numerical results clearly indicate that we are able to get an accurate numerical approximation of the BVP (3) noniteratively. This can be also realized by comparing the exact solution, for , plotted on figure 1 with the numerical solution shown on figure 5. Figure 5 is a frame, on , of the numerical solution, computed by the RK6 solver with .
From figure 5 we can notice how the last step, which is the one close to the origin, is smaller with respect to the previous ones. This is due to the application of our simple event locator defined by equations (20)(21).
In table 3 we report the numerical results obtained by fixing a value for the step size and repeating the computation for several decreasing values of .
1.0  2.  39790  5.  16  1.  099999948 
1.0  4.  61512  5.  35  1.  009999946 
1.0  6.  90875  9.  26  1.  000999907 
1.0  9.  21044  1.  23  1.  000099877 
1.0  11.  5129  3.  02  1.  000009970 
1.0  13.  8155  1.  25  1.  000000875 
1.0  16.  1181  4.  33  1.  000000057 
1.0  18.  4207  1.  09  0.  999999901 
1.0  20.  7233  9.  76  0.  999999903 
As it is easily seen, we can verify numerically that, the smaller the value of the larger the free boundary value .
5.2 Scaling invariance
As an example we consider here the noniterative numerical solution of the free boundary formulation (3). With a simple change of variables, the problem (3) clearly belongs to (4.2). The governing differential equation is invariant on condition that , that is .
In table 4 we list the numerical results obtained by setting , , and . We used a uniform grid with and the same RungeKutta solvers of the previous subsection.
Method  

RK4  8.  24683  9.  999999539  1.  000000092 
RK6  8.  24462  10.  000000058  0.  999999988 
RK8  8.  24461  9.  999999959  1.  000000008 
The reported numerical results clearly indicate that we are able to get an accurate numerical approximation of the BVP (3) noniteratively. This can be also realized by comparing the exact solution, for , plotted on figure 3 with the numerical solution shown on the bottom frame of figure 6. In figure 6 we plot the numerical solution of the initial value problem (2), obtained by RK6, as well as the rescaled solution for the free boundary formulation of the BVP (3).
We would like to remark that, in order to solve the free boundary formulation noniteratively we can use a whatever large value of . For instance, by setting and again we have computed with RK6 the values , , and . Indeed, we got a less accurate value for the first derivative of at , and this depends on the use of the same grid on although we have to cope with a faster transitory of near the origin, see figure 7. However, we verify numerically that the smaller the value of the larger the free boundary value .
6 Concluding remarks
In this paper we have reported on the free boundary formulation for BVPs defined on a semiinfinite interval. In this context, we discussed two simple test problems where we are able to get the exact analytical solution of both the original BVPs and of their free boundary formulation. Therefore, these problems can be used as benchmarks for numerical methods applied to BVPs on a semiinfinite interval and to free BVPs. Moreover, we emphasized how for two classes of free BVPs we can define noniterative initial value methods, whereas BVPs are usually solved iteratively. These noniterative methods can be deduced within Lie’s group invariance theory. We also applied the noniterative methods to the two introduced free boundary formulations and we obtained meaningful numerical results.
Let us indicate a few problems in the applied sciences where our nonITMs can be used. For instance, the nonITM defined using the translation group can be applied to the free boundary formulation of the BVP
where is a positive constant, arising in the theory of colloids, see Countryman and Kannan [52] and the references quoted therein. Moreover, the nonITM defined using the scaling group can be applied to the free boundary formulation of the BVP, studied by Seshadri and Na [57],
where and are physical parameters, arising in the study of longitudinal impact to a nonlinear viscoplastic thin rod. We would like to mention that the nonITMs defined above can be easily generalized to the free boundary formulation of BVPs governed by systems of differential equations like those derived by Dresner [58, pp. 6997] within the study of wave propagation problems.
Finally, let us report on a possible way to extend the noniterative methods defined in this paper to free boundary problems that are invariant with respect to a generic Lie group. To this end we assume that a given problem is invariant under the oneparametr group of point transformations
(31) 
where is the group parameter and and are the socalled group generators. The infinitesimal generator of (31) is given by
(32) 
If we introduce the variable transformation and , then we can rewrite the infinitesimal generator (32) as
(33) 
Now we can choose between two possible alternative: the socalled canonical variables, see Bluman and Kumei [32], that transforms as
(34) 
or the so called normal variables, see Fazio [40], that obey the relations
(35) 
We note that (35) is the infinitesimal form of a scaling group. In the case of the normal variable, taking into account the equations (33) and (35), we get
The general solution of (6) can be found by integrating the characteristic equations
(37) 
By integrating the first equation in (37) we get the socalled first group invariant , whereupon we have
where and are arbitrary functions of the first invariant. Fazio [59] gives an application of the idea developed above to the numerical length estimation of tubular flow reactors.
Acknowledgements: This work was partially supported by GNCS of INDAM and by University of Messina.
References
 L. Horwarth. On the solution of the laminar boundary layer equations. Proc. Roy. Soc. London A, 164:547–579, 1938.
 S. Goldstein. Modern Developments in Fluid Dynamics. Clarendon Press, Oxford, 1938.
 H. Blasius. Grenzschichten in Flüssigkeiten mit kleiner Reibung. Z. Math. Phys., 56:1–37, 1908.
 L. Fox. Numerical Solution of Twopoint Boundary Value Problems in Ordinary Differential Equations. Clarendon Press, Oxford, 1957.
 L. Collatz. The Numerical Treatment of Differential Equations. Springer, Berlin, 3rd edition, 1960.
 M. Lentini and H. B. Keller. The von Karman swirling flows. SIAM J. Appl. Math., 38:52–64, 1980.
 F. R. de Hoog and R. Weiss. An approximation theory for boundary value problems on infinite intervals. Computing, 24:227–239, 1980.
 M. Lentini and H. B. Keller. Boundary value problems on semiinfinite intervals and their numerical solutions. SIAM J. Numer. Anal., 17:577–604, 1980.
 P. A. Markowich. A theory for the approximation of solution of boundary value problems on infinite intervals. SIAM J. Math. Anal., 13:484–513, 1982.
 P. A. Markowich. Analysis of boundary value problems on infinite intervals. SIAM J. Math. Anal., 14:11–37, 1983.
 J. R. Ockendon. Math. Rev., 84c:34021.
 W. J. Beyn. Global bifurcation and their numerical computation. In D. Rossed, B. D. Dier, and A. Spence, editors, Bifurcation: Numerical Techniques and Applications, pages 169–181. Kluwer, Dordrecht, 1990.
 W. J. Beyn. The numerical computation of connecting orbits in dynamical systems. IMA J. Numer. Anal., 9:379–405, 1990.
 W. J. Beyn. Numerical methods for dynamical systems. In W. Light, editor, Advances in Numerical Analysis, pages 175–236. Clarendon Press, Oxford, 1992.
 M. J. Friedman and E. J. Doedel. Numerical computation and continuation of invariant manifolds connecting fixed points. SIAM J. Numer. Anal., 28:789–808, 1991.
 F. Bai, A. Spence, and A. M. Stuart. The numerical computation of heteroclinic connections in systems of gradient partial differential equations. SIAM J. Appl. Math., 53:743–769, 1993.
 L. Liu, G. Moore, and R. D. Russell. Computation and continuation of homoclinic and heteroclinic orbits with arclength parameterization. SIAM J. Sci. Comput., 18:69–93, 1997.
 H. G. Landau. Heat conduction in melting solid. Quart. Appl. Math., 8:81–94, 1950.
 R. Fazio. The Blasius problem formulated as a free boundary value problem. Acta Mech., 95:1–7, 1992.
 P. D. Ariel. Stagnation point flow  A free boundary value problem formulation. Int. J. Comput. Math., 49:123–131, 1993.
 R. Fazio. The FalknerSkan equation: numerical solutions within group invariance theory. Calcolo, 31:115–124, 1994.
 J. Zhang and B. Chen. An iterative method for solving the FalknerSkan equation. Appl. Math. Comput., 210:215–222, 2009.
 S. Zhu, Q. Wu, and X. Cheng. Numerical solution of the FalknerSkan equation based on quasilinearization. Appl. Math. Comput., 215:2472–2485, 2009.
 R. Fazio. A novel approach to the numerical solution of boundary value problems on infinite intervals. SIAM J. Numer. Anal., 33:1473–1483, 1996.
 R. Fazio. A survey on free boundary identification of the truncated boundary in numerical BVPs on infinite intervals. J. Comput. Appl. Math., 140:331–344, 2002.
 R. Fazio. A free boundary approach and Keller’s box scheme for BVPs on infinite intervals. Int. J. Computer Math., 80:1549–1560, 2003.
 C. E. Grosch and S. A. Orszag. Numerical solution of problems in unbounded regions: Coordinate transforms. J. Comput. Phys., 25:273–296, 1977.
 M. N. Koleva. Numerical solution of the heat equation in unbounded domains using quasiuniform grids. In I. Lirkov, S. Margenov, and J. Wasniewski, editors, Largescale Scientific Computing, volume 3743 of Lecture Notes in Comput. Sci., pages 509–517, 2006.
 R. Fazio and A. Jannelli. Finite difference schemes on quasiuniform grids for BVPs on infinite intervals. J. Comput. Appl. Math., 269:14–23, 2014.
 R. Fazio. A similarity approach to the numerical solution of free boundary problems. SIAM Rev., 40:616–635, 1998.
 G. W. Bluman and J. D. Cole. Similarity Methods for Differential Equations. Springer, Berlin, 1974.
 G. W. Bluman and S. Kumei. Symmetries and Differential Equations. Springer, Berlin, 1989.
 G. I. Barenblatt. Scaling, SelfSimilarity and Intermediate Asymptotics. Cambridge University Press, Cambridge, 1996.
 L. Dresner. Applications of Lie’s Theory of Ordinary and Partial Differential Equations. Institute of Physics Publishing, London, 1999.
 K. Töpfer. Bemerkung zu dem Aufsatz von H. Blasius: Grenzschichten in Flüssigkeiten mit kleiner Reibung. Z. Math. Phys., 60:397–398, 1912.
 M. S. Klamkin. On the transformation of a class of boundary value problems into initial value problems for ordinary differential equations. SIAM Rev., 4:43–47, 1962.
 T. Y. Na. Transforming boundary conditions to initial conditions for ordinary differential equations. SIAM Rev., 9:204–210, 1967.
 T. Y. Na. Further extension on transforming from boundary value to initial value problems. SIAM Rev., 20:85–87, 1968.
 R. Fazio and D. J. Evans. Similarity and numerical analysis for free boundary value problems. Int. J. Computer Math., 31:215–220, 1990. 39 : 249, 1991.
 R. Fazio. Normal variables transformation method applied to free boundary value problems. Int. J. Computer Math., 37:189–199, 1990.
 R. Fazio and S. Iacono. On the translation groups and noniterative transformation methods. In E. De Bernardis, R. Spligher, and V. Valenti, editors, Applied and Industrial Mathematics in Italy III, pages 331–340, Singapore, 2010. World Scientific.
 R. Fazio. A moving boundary hyperbolic problem for a stress impact in a bar of ratetype material. Wave Motion, 16:299–305, 1992.
 V. G. Fox, L. E. Erickson, and L. I. Fan. The laminar boundary layer on a moving continuous flat sheet in a nonnewtonian fluid. AIChE J., 15:327–333, 1969.
 G. H. Meyer. Initial Value Methods for Boundary Value Problems; Theory and Application of Invariant Imbedding. Academic Press, New York, 1973.
 T. Y. Na. Computational Methods in Engineering Boundary Value Problems. Academic Press, New York, 1979.
 P. L. Sachdev. Nonlinear Ordinary Differential Equations and their Applications. Marcel Dekker, New York, 1991.
 R. Seshadri and T. Y. Na. Group Invariance in Engineering Boundary Value Problems. Springer, New York, 1985.
 G. Birkhoff. Hydrodynamics: A Study on Logic, Fact and Similitude. Princeton University Press, 1950, Princeton, 2nd edition, 1960.
 T. Y. Na. An initial value method for the solution of a class of nonlinear equations in fluid mechanics. J. Basic Engrg. Trans. ASME, 92:503–509, 1970.
 R. Fazio. The FalknerSkan equation: numerical solutions within group invariance theory. Calcolo, 31:115–124, 1994.
 R. Fazio. Numerical transformation methods: a constructive approach. J. Comput. Appl. Math., 50:299–303, 1994.
 M. Countryman and R. Kannan. A class of nonlinear boundary value problems on the half line. Comp. Math. Appl., 28:121–130, 1994.
 U. M. Ascher and R. D. Russell. Reformulation of boundary value problems into “standard” form. SIAM Rev., 23:238–254, 1981.
 C. Runge. Ueber die numerische Auflösung von Differentialgleichungen. Mathematische Annalen, 46:167–178, 1895.
 M. W. Kutta. Beitrag zur näherungweisen Integration totaler Differentialgleichungen. Zeitschrift für Mathematik und Physik, 46:435–453, 1901.
 J. C. Butcher. Numerical Methods for Ordinary Differential Equations. Whiley, Chichester, 2003.
 R. Seshadri and T. Y. Na. Invariant solution for nonlinear viscoplastic impact. Ind. Math., 25:37–44, 1975.
 L. Dresner. Similarity Solutions of Nonlinear Partial Differential Equations, volume 88 of Research Notes in Math. Pitman, London, 1983.
 R. Fazio. Numerical length estimation for tubular flow reactors. J. Comput. Appl. Math., 41:313–321, 1992.