Buggy MacCormack Explicit Scheme for Compressible NS Equations
Outputs horizontal velocity contours. Can generate "wave-like" behavior before scheme fails. Not sure if the oscillation near walls is from bad BC's or buried somewhere in Predictor/Corrector stepping. Designed using Pletcher 3rd edition, reworked for general 2D coordinates. Code available on request (Fortran 90 w/OpenMP)
Ideas?