A pseudospectral based method of lines for solving integro-differential boundary-layer equations. Application to the mixed convection over a heated horizontal plate

Abstract The method of lines is well suited for solving numerically parabolic boundary-layer equations because it avoids the numerical difficulties associated to the integration of the continuity equation, which is subsumed into the momentum equations as an integral of the main velocity component. To deal with these integrals, as well as with any other integral operator entering the boundary layer equations in some particular problems, it is very efficient to discretize the transversal coordinate using pseudospectral methods. The resulting ordinary differential equations (ODEs) can be then written in a very compact form, suitable for general-purpose methods and software developed for the numerical integration of ODEs. We present here such a numerical method applied to the boundary-layer equations governing the mixed convection over a heated horizontal plate. These parabolic equation can be written in such a way that the natural convection appears as an integro-differential term in the usual horizontal momentum equation, so that the discretization by pseudospectral methods of the vertical coordinate derivative is very appropriate. Several Matlab based solvers are compared to integrate the resulting ODEs. To validate the numerical results they are compared with analytic solutions valid near the leading edge of the boundary-layer.

[1]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[2]  L. Richardson.  The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[3]  Elsayed M. E. Elbarbary.  Chebyshev finite difference method for the solution of boundary-layer equations , 2005, Appl. Math. Comput..

[4]  A Collocation Finite Element Method for Integration of the Boundary Layer Equations , 1984 .

[5]  R. Vaillancourt,et al.  Behind and beyond the Matlab ODE suite , 2000 .

[6]  F. S. V. Bazán.  Chebyshev pseudospectral method for computing numerical solution of convection-diffusion equation , 2008, Appl. Math. Comput..

[7]  H. Steinrück.  Mixed convection over a cooled horizontal plate: non-uniqueness and numerical instabilities of the boundary-layer equations , 1994, Journal of Fluid Mechanics.

[8]  U. Ghia,et al.  Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications , 2007 .

[9]  N. Afzal,et al.  Mixed Convection Over a Horizontal Plate , 1984 .

[10]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[11]  W. Schneider,et al.  A similarity solution for combined forced and free convection flow over a horizontal plate , 1979 .

[12]  Vanessa Hertzog,et al.  Computational Fluid Mechanics And Heat Transfer , 2016 .

[13]  J. Schetz.  Boundary Layer Analysis , 1992 .

[14]  W. Gill,et al.  A theoretical investigation of natural convection effects in forced horizontal flows , 1962 .

[15]  H. Fernholz.  Boundary Layer Theory , 2001 .

[16]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[17]  Robert L. Ash,et al.  Application of spectral collocation techniques to the stability of swirling flows , 1989 .

[18]  Ioan Pop,et al.  Similarity solutions for mixed convection boundary layer flow over a permeable horizontal flat plate , 2010, Appl. Math. Comput..

[19]  Ephraim M Sparrow,et al.  Buoyancy effects on horizontal boundary-layer flow and heat transfer , 1962 .

[20]  L. Trefethen.  Spectral Methods in MATLAB , 2000 .

[21]  L. Trefethen,et al.  Stability of the method of lines , 1992, Spectra and Pseudospectra.