A Quasi-Parabolic technique for computation of three-dimensional viscous flows

A computational technique is presented for obtaining flowfield solutions to a parabolic form of the Navier-Stokes equations. The point of departure is the General Interpolant Method (GIM) which provides a discretization for partial differential equations on arbitrary three-dimensional geometries. The new scheme, termed Quasi-Parabolic, treats the parabolized equations but with 'time-like' terms appended. Addition of these extra terms, which are relaxed by iteration, avoid many of the singularities inherent in classical parabolic Navier-Stokes methods. Streamwise derivatives are approximated by three-point backward differences and the cross plane operators use an alternating forward-backward sweep. A two-step sequence is used to implement the difference scheme in the spatial dimensions and a time-like relaxation converges the quasi-parabolic procedure at each plane. Solutions are presented for flows in two and three dimensions. Inviscid flows are solved for internal and external applications and viscous flows in boundary layers and free shear layers are also computed with the GIM/Quasi-Parabolic scheme.