Premium
Numerical techniques for flow problems with singularities
Author(s) -
Farrell Paul A.,
Hegarty Alan F.,
Miller John J. H.,
O'Riordan Eugene,
Shishkin Grigorii I.
Publication year - 2003
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.536
Subject(s) - mathematics , boundary layer , polygon mesh , mathematical analysis , boundary value problem , numerical analysis , prandtl number , uniform convergence , geometry , physics , mechanics , computer science , heat transfer , computer network , bandwidth (computing)
Abstract This paper deals with grid approximations to Prandtl's boundary value problem for boundary layer equations on a flat plate in a region including the boundary layer, but outside a neighbourhood of its leading edge. The perturbation parameter ε= Re‐1 takes any values from the half‐interval (0,1] ; here Re is the Reynolds number. To demonstrate our numerical techniques we consider the case of the self‐similar solution. By using piecewise uniform meshes, which are refined in a neighbourhood of the parabolic boundary layer, we construct a finite difference scheme that converges ε‐uniformly. We present the technique of experimental substantiation of ε‐uniform convergence for both the numerical solution and its normalized (scaled) difference derivatives, outside a neighbourhood of the leading edge of the plate. By numerical experiments we demonstrate the efficiency of numerical techniques based on the fitted mesh method. We discuss also the applicability of fitted operator methods for the numerical approximation of the Prandtl problem. It is shown that the use of meshes refined in the parabolic boundary layer region is necessary for achieving ε‐uniform convergence. Copyright © 2003 John Wiley & Sons, Ltd.