Pseudospectral solutions to some singular nonlinear BVPs

We solve by Laguerre collocation method two second order genuinely nonlinear singularly perturbed boundary value problems and a fourth order nonlinear one formulated on the half-line. First, we are concerned with solutions to a model equation for low Reynolds number flow, the so-called Cohen, Fokas, Lagerstrom model. Then we compute monotonously increasing solutions, the so-called “bubble type solutions”, to density profile equation for the description of the formation of microscopical bubbles in a nonhomogeneous fluid. We compare the results obtained in the latter case with those carried out by shooting or polynomial collocation both coupled with domain truncation. A fourth order nonlinear boundary value problem supplied with mixed boundary conditions in origin as well as at infinity is also solved. The solution and its first three derivatives are obtained and compared with their counterparts obtained by Keller’s box scheme. The collocation based on Laguerre functions avoids the domain truncation, exactly impose any kind of boundary conditions at infinity and resolve with reliable accuracy the sharp interior or exterior layer of the solutions. The method is fairly easy implementable and efficient in treating various types of nonlinearities. The accuracy of the method is also tested on the classical Blasius problem (a third order one). We solve with the same method the eigenvalue problem obtained by linearization around the constant solution, which corresponds to the case of a homogeneous fluid (without bubbles), and observe that this solution is stable.

[1]  George C. Hsiao Singular Perturbations for a Nonlinear Differential Equation with a Small Parameter , 1973 .

[2]  Tao Tang,et al.  Boundary Layer Resolving Pseudospectral Methods for Singular Perturbation Problems , 1996, SIAM J. Sci. Comput..

[3]  Riccardo Fazio,et al.  A survey on free boundary identification of the truncated boundary in numerical BVPs on infinite intervals , 2002 .

[4]  M. L. Morgado,et al.  Bubbles and droplets in nonlinear physics models: Analysis and numerical simulation of singular nonlinear boundary value problem , 2008 .

[5]  P. M. Limaa,et al.  Analytical – numerical investigation of bubble-type solutions of nonlinear singular problems , 2006 .

[6]  Gennady Y. Kulikov,et al.  Analysis and numerical approximation of singular boundary value problems with the p-Laplacian in fluid mechanics , 2014, J. Comput. Appl. Math..

[7]  P. Henrici Bounds for iterates, inverses, spectral variation and fields of values of non-normal matrices , 1962 .

[8]  Françoise Chaitin-Chatelin,et al.  Lectures on finite precision computations , 1996, Software, environments, tools.

[9]  Riccardo Fazio,et al.  A free boundary approach and keller's box scheme for BVPs on infinite intervals , 2003, Int. J. Comput. Math..

[10]  J. Boyd,et al.  Pseudospectral methods on a semi-infinite interval with application to the Hydrogen atom: a comparison of the mapped Fourier-sine method with Laguerre series and rational Chebyshev expansions , 2003 .

[11]  Joost Rommes,et al.  Application of the Jacobi–Davidson method to accurate analysis of singular linear hydrodynamic stability problems , 2013 .

[12]  Athanassios S. Fokas,et al.  Proof of some asymptotic results for a model equation for low Reynolds number flow , 1978 .

[13]  Othmar Koch,et al.  Efficient Numerical Solution of the Density Profile Equation in Hydrodynamics , 2007, J. Sci. Comput..

[14]  Satish C. Reddy,et al.  A MATLAB differentiation matrix suite , 2000, TOMS.

[15]  Peter A. Markowich Analysis of Boundary Value Problems on Infinite Intervals , 1983 .

[16]  Calin-Ioan Gheorghiu Laguerre collocation solutions to boundary layer type problems , 2012, Numerical Algorithms.

[17]  H. B. Keller,et al.  Boundary Value Problems on Semi-Infinite Intervals and Their Numerical Solution , 1980 .

[18]  Steven Pruess,et al.  Mathematical software for Sturm-Liouville problems , 1993, TOMS.

[19]  Lloyd N. Trefethen,et al.  Computation of pseudospectra , 1999, Acta Numerica.

[20]  Donal O'Regan Solvability of some singular boundary value problems on the semi-infinite interval , 1996 .

[21]  L. A. Rubel An estimate of the error due to the truncated boundary in the numerical solution of the Blasius equation , 1955 .