3-D composite velocity solutions for subsonic/transonic flows

A composite velocity procedure for the three-dimensional reduced Navier-Stokes equations is developed. In the spirit of matched asymptotic expansions, the velocity components are written as a combination multiplicative and additive composite of viscouslike velocities and pseudopotential or inviscid velocities. The solution procedure is then consistent with both asymptotic inviscid flow and boundary layer theory. For transonic flow cases, the Enquist-Osher flux biasing scheme developed for the full potential equation is used. A quasi-conservation form of the governing equations is used in the shock region to capture the correct rotational shock with the standard nonconservation form of the equations used in nonshock regions. The consistent coupled strongly implicit procedure coupled with a plane relaxation procedure is used to solve the discretized equations.