A NEW FLUX-CONSERVING NEWTON SCHEME FOR THE NAVIER-STOKES EQUATIONS

Publisher Website
Google Scholar
Abstract
A new numerical method is developed for the two-dimensional, steady Navier-Stokes equations. Using local polynomial expansions to represent the discrete primitive variables on each cell, we construct a scheme which has the following properties: First, the local discrete primitive variables are functional solutions of both the integral and differential forms of the Navier-Stokes equations. Second, fluxes are balanced across cell interfaces using exact functional expressions (to the order of accuracy of the local expansions). No interpolation, flux models, or flux limiters are required. Third, local and global conservation of mass, momentum, and energy are explicitly provided for. Finally, the discrete primitive variables and their derivatives are treated in a unified and consistent manner. All are treated as unknowns to be solved together for simulating the local and global flux conservation. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented. As a special case, the formulation is applied to incompressible flow, and a Newton's method scheme is developed for the solution of laminar channel flow. H is shown that, at Reynolds numbers of 100, 1000, and 2000, the developing channel flow boundary layer can be accurately resolved using as few as six to ten cells per channel width.