Remarks on the numerical treatment of Prandtl's boundary-layer equations (Q1096483)

The Crocco transformation for the boundary-layer equations very successfully used by Oleinik and Nickel for theoretical purposes (existence of 2-D, uniqueness of 2-D and 3-D solutions) will be used here for the construction of numerical procedures. An approximate discretization of the transformed equations leads to a nicely structured system of nonlinear equations that can be devided into small parts to be solved one after the other (in the 3-D case even partially by parallel processing). The Jacobians of these parts are often M-matrices such that SOR iterations work. The inverse transformation can numerically be realized very simply.