TY - JOUR

T1 - Robust numerical methods for boundary-layer equations for a model problem of flow over a symmetric curved surface

AU - Ansari, A.R.

AU - Hossain, B.

AU - Koren, B.

AU - Shishkin, G.I.

PY - 2006

Y1 - 2006

N2 - We investigate the model problem of flow of a viscous incompressible fluid past a symmetric curved surface when the flow is parallel to its axis. This problem is known to exhibit boundary layers. Also the problem does not have solutions in closed form, it is modelled by boundary-layer equations. Using a self-similar approach based on a Blasius series expansion (up to three terms), the boundary-layer equations can be reduced to a Blasius-type problem consisting of a system of eight third-order ordinary differential equations on a semi-infinite interval. Numerical methods need to be employed to attain the solutions of these equations and their derivatives, which are required for the computation of the velocity components, on a finite domain with accuracy independent of the viscosity v, which can take arbitrary values from the interval (0,1]. To construct a robust numerical method we reduce the original problem on a semi-infinite axis to a problem on the finite interval [0, K], where K = K(N) = ln N. Employing numerical experiments we justify that the constructed numerical method is parameter robust.

AB - We investigate the model problem of flow of a viscous incompressible fluid past a symmetric curved surface when the flow is parallel to its axis. This problem is known to exhibit boundary layers. Also the problem does not have solutions in closed form, it is modelled by boundary-layer equations. Using a self-similar approach based on a Blasius series expansion (up to three terms), the boundary-layer equations can be reduced to a Blasius-type problem consisting of a system of eight third-order ordinary differential equations on a semi-infinite interval. Numerical methods need to be employed to attain the solutions of these equations and their derivatives, which are required for the computation of the velocity components, on a finite domain with accuracy independent of the viscosity v, which can take arbitrary values from the interval (0,1]. To construct a robust numerical method we reduce the original problem on a semi-infinite axis to a problem on the finite interval [0, K], where K = K(N) = ln N. Employing numerical experiments we justify that the constructed numerical method is parameter robust.

U2 - 10.1080/13926292.2006.9637324

DO - 10.1080/13926292.2006.9637324

M3 - Article

VL - 11

SP - 365

EP - 378

JO - Mathematical Modelling and Analysis

JF - Mathematical Modelling and Analysis

SN - 1392-6292

IS - 4

ER -