Newton's method solver for high-speed viscous separated flowfields

Abstract

A new method for calculating the 2-D, laminar Navier-Stokes equations is presented. The method uses Newton's method for nonlinear systems of equations to find steady-state solutions. The Navier-Stokes equations are approximated by finite differences using Roe's flux difference splitting. Second-order accuracy is attained by using Spekreijse's interpolation with Van Albada's limiter. The exact Newton's method Jacobian matrix is inverted by using recent sparse matrix routines. The symbolic manipulation package MACSYMA is used to develop and write the FORTRAN code