Discretized Newton-relaxation solution of high resolution flux-difference split schemes

A method is presented for the implicit solution of high resolution flux-difference split schemes for which it is impractical to obtain the Jacobian matrix analytically. The method is used frequently in numerical schemes for solving nonlinear systems of equations, but apparently has not been used to date for approximating the Jacobian of the Roe numerical flux vectors. The method consists of obtaining a discretized Jacobian by using simple finite-differences of the flux vectors, and then using this discretized Jacobian in a Newton-relaxation solution method. The method is applied to both compressible and incompressible flows. Numerical results are presented that demonstrate the method works for both inviscid and viscous flows.