Implementation of Vigneron's streamwise pressure gradient approximation in the PNS equations

For single sweep parabolized Navier-Stokes solvers, the streamwise pressure gradient must be modified in the subsonic region to eliminate numerical instabilities. The accuracy of this modification on the solution of the parabolized Navier-Stokes equations with Vigneron's technique is shown to depend on how the numerical approximation of the pressure gradient is formed. A simple test case of supersonic laminar flow over a flat plate is computed with two different numerical methods for solving the PNS equations. Significant errors in the temperature profile and skin friction coefficient are demonstrated using a fully conservative differencing treatment of Vigneron's splitting for the pressure gradient typically used in parabolized Navier-Stokes solvers. The physical reason for this error is discussed. An alternate formulation is demonstrated which minimizes these errors.