The non-oscillatory central difference scheme of Nessyahu and Tadmor, in which the resolution of Riemann problems at the cell interfaces is by-passed thanks to the use of the staggered Lax-Friedrichs scheme, is extended here to a two-step, two-dimensional non-oscillatory centered scheme in finite volume formulation. The construction of the scheme rests on a finite volume extension of the Lax-Friedrichs scheme, in which the finite volume cells are the barycentric cells constructed around the nodes of an FEM triangulation, for odd time steps, and some quadrilateral cells associated with this triangulation, for even time steps.
Piecewise linear cell interpolants using least-squares gradients combined with a van Leer-type slope limiting allow for an oscillation-free second-order resolution.
The method is illustrated by an air flow calculation around a 2-D double ellipse; and airfoil NACA0012. For more detail see [1], [2], [3].
The authors will now focus their efforts on applications to problems with the Navier-Stokes equations in conservation form for compressible flows [4] and three-dimensional transonic and supersonic flow problems (see Figure 2).
Figure 1
Mesh (4318 vertices) and Mach contours (Mach = 0.85 and 1° of attack)
Figure 3
Euler flow around a double ellipse. Original grid, barycentric cells Ci and quadrilateral cells Lij
Figure 4
Double ellipse: Initial mesh (1558 vertices) and solution (pressure and Mach contours)
Figure 5
Double ellipse: Final mesh (5055 vertices) and solution (pressure and Mach contours)
Table des matières de la petite galerie du CRM
28 January 1998, webmaster@CRM.UMontreal.CA