A numerical method for solving the equations of compressible viscous flow

Although much progress has already been made in solving problems in aerodynamic design, many new developments are still needed before the equations for unsteady compressible viscous flow can be solved routinely. This paper describes one such development. A method for solving these equations has been devised that: (1) is second-order accurate in space and time; (2) is unconditionally stable; (3) preserves conservation form; (4) requires no block or scalar tridiagonal inversions; (5) is simple and straightforward to program (estimated 10% modification for the update of many existing programs); (6) is more efficient than present methods; and (7) should easily adapt to current and future computer architectures. Computational results for laminar and turbulent flows at Reynolds numbers from 300,000 to 3 x 10 to the 7th and at CFL numbers as high as 1000 are compared with theory and experiment.