Order Accurate Finite Volume Schemes Research Articles

Modelling the effects of diffusive-viscous waves in a 3-D fluid-saturated media using two numerical approaches

SUMMARYThe accurate numerical modelling of 3-D seismic wave propagation is essential in understanding details to seismic wavefields which are, observed on regional and global scales on the Earth’s surface. The diffusive-viscous wave (DVW) equation was proposed to study the connection between fluid saturation and frequency dependence of reflections and to characterize the attenuation property of the seismic wave in a fluid-saturated medium. The attenuation of DVW is primarily described by the active attenuation parameters (AAP) in the equation. It is, therefore, imperative to acquire these parameters and to additionally specify the characteristics of the DVW. In this paper, quality factor, Q is used to obtain the AAP, and they are compared to those of the visco-acoustic wave. We further derive the 3-D numerical schemes based on a second order accurate finite-volume scheme with a second order Runge–Kutta approximation for the time discretization and a fourth order accurate finite-difference scheme with a fourth order Runge–Kutta approximation for the time discretization. We then simulate the propagation of seismic waves in a 3-D fluid-saturated medium based on the derived schemes. The numerical results indicate stronger attenuation when compared to the visco-acoustic case.

Adaptive-Mesh-Refinement for hyperbolic systems of conservation laws based on a posteriori stabilized high order polynomial reconstructions

In this paper we propose a third order accurate finite volume scheme based on a posteriori limiting of polynomial reconstructions within an Adaptive-Mesh-Refinement (AMR) simulation code for hydrodynamics equations in 2D. The a posteriori limiting is based on the detection of problematic cells on a so-called candidate solution computed at each stage of a third order Runge–Kutta scheme. Such detection may include different properties, derived from physics, such as positivity, from numerics, such as a non-oscillatory behavior, or from computer requirements such as the absence of NaN's. Troubled cell values are discarded and re-computed starting again from the previous time-step using a more dissipative scheme but only locally, close to these cells. By locally decrementing the degree of the polynomial reconstructions from 2 to 0 we switch from a third-order to a first-order accurate but more stable scheme. The entropy indicator sensor is used to refine/coarsen the mesh. This sensor is also employed in an a posteriori manner because if some refinement is needed at the end of a time step, then the current time-step is recomputed with the refined mesh, but only locally, close to the new cells. We show on a large set of numerical tests that this a posteriori limiting procedure coupled with the entropy-based AMR technology can maintain not only optimal accuracy on smooth flows but also stability on discontinuous profiles such as shock waves, contacts, interfaces, etc. Moreover numerical evidences show that this approach is at least comparable in terms of accuracy and cost to a more classical CWENO approach within the same AMR context.

High Order Accurate Direct Arbitrary-Lagrangian-Eulerian ADER-MOOD Finite Volume Schemes for Non-Conservative Hyperbolic Systems with Stiff Source Terms

AbstractIn this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-conservative products and stiff source terms. This scheme is constructed with a single stencil polynomial reconstruction operator, a one-step space-time ADER integration which is suitably designed for dealing even with stiff sources, a nodal solver with relaxation to determine the mesh motion, a path-conservative integration technique for the treatment of non-conservative products and ana posterioristabilization procedure derived from the so-called Multidimensional Optimal Order Detection (MOOD) paradigm. In this work we consider the seven equation Baer-Nunziato model of compressible multi-phase flows as a representative model involving non-conservative products as well as relaxation source terms which are allowed to become stiff. The new scheme is validated against a set of test cases on 2D/3D unstructured moving meshes on parallel machines and the high order of accuracy achieved by the method is demonstrated by performing a numerical convergence study. Classical Riemann problems and explosion problems with exact solutions are simulated in 2D and 3D. The overall numerical code is also profiled to provide an estimate of the computational cost required by each component of the whole algorithm.

Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers

Several advances have been reported in the recent literature on divergence-free finite volume schemes for Magnetohydrodynamics (MHD). Almost all of these advances are restricted to structured meshes. To retain full geometric versatility, however, it is also very important to make analogous advances in divergence-free schemes for MHD on unstructured meshes. Such schemes utilize a staggered Yee-type mesh, where all hydrodynamic quantities (mass, momentum and energy density) are cell-centered, while the magnetic fields are face-centered and the electric fields, which are so useful for the time update of the magnetic field, are centered at the edges.Three important advances are brought together in this paper in order to make it possible to have high order accurate finite volume schemes for the MHD equations on unstructured meshes. First, it is shown that a divergence-free WENO reconstruction of the magnetic field can be developed for unstructured meshes in two and three space dimensions using a classical cell-centered WENO algorithm, without the need to do a WENO reconstruction for the magnetic field on the faces. This is achieved via a novel constrained L2-projection operator that is used in each time step as a postprocessor of the cell-centered WENO reconstruction so that the magnetic field becomes locally and globally divergence free. Second, it is shown that recently-developed genuinely multidimensional Riemann solvers (called MuSIC Riemann solvers) can be used on unstructured meshes to obtain a multidimensionally upwinded representation of the electric field at each edge. Third, the above two innovations work well together with a high order accurate one-step ADER time stepping strategy, which requires the divergence-free nonlinear WENO reconstruction procedure to be carried out only once per time step.The resulting divergence-free ADER–WENO schemes with MuSIC Riemann solvers give us an efficient and easily-implemented strategy for divergence-free MHD on unstructured meshes. Several stringent two- and three-dimensional problems are shown to work well with the methods presented here.

Modeling Interactions and Shoaling of Solitary Waves Using a Hybrid Finite Volume and Finite Difference Solver

This paper presents a mixed finite volume and finite difference solver with results showing the solitary wave interactions and shoaling process by solving a set of conservative forms of Boussinesq equations. A second order accurate finite volume scheme is applied to the conservative terms of the governing equations while up to the second order finite difference formulations are used to discretize the dispersive source terms with higher order derivatives. The limiters and surface gradient method are implemented in the model to remove the unwanted spurious oscillations and preserve the still water condition without introducing errors at the interfaces. The performance of the present numerical solver is tested with results of head on collisions and shoaling of solitary waves compared against those from finite element models that were developed based on fully nonlinear weakly dispersive and weakly nonlinear weakly dispersive forms of the Boussinesq equations as well as analytical solutions and experimental observations.

Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms

In [16,17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions), it is more difficult to design high order schemes which do not produce negative density or pressure. In this paper, we first show that our framework to construct positivity-preserving high order schemes in [16,17] can also be applied to Euler equations with a general equation of state. Then we discuss an extension to Euler equations with source terms. Numerical tests of the third order Runge–Kutta DG (RKDG) method for Euler equations with different types of source terms are reported.

A Genuinely High Order Total Variation Diminishing Scheme for One-Dimensional Scalar Conservation Laws

It is well known that finite difference or finite volume total variation diminishing (TVD) schemes solving one-dimensional scalar conservation laws degenerate to first order accuracy at smooth extrema [S. Osher and S. Chakravarthy, SIAM J. Numer. Anal., 21 (1984), pp. 955–984], thus TVD schemes are at most second order accurate in the $L^1$ norm for general smooth and nonmonotone solutions. However, Sanders [Math. Comp., 51 (1988), pp. 535–558] introduced a third order accurate finite volume scheme which is TVD, where the total variation is defined by measuring the variation of the reconstructed polynomials rather than the traditional way of measuring the variation of the grid values. By adopting the definition of the total variation for the numerical solutions as in [R. Sanders, Math. Comp., 51 (1988), pp. 535–558], it is possible to design genuinely high order accurate TVD schemes. In this paper, we construct a finite volume scheme which is TVD in this sense with high order accuracy (up to sixth order) in the $L^1$ norm. Numerical tests for a fifth order accurate TVD scheme will be reported, which include test cases from traffic flow models.

Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier–Stokes equations

In this paper, we propose a new unified family of arbitrary high order accurate explicit one-step finite volume and discontinuous Galerkin schemes on unstructured triangular and tetrahedral meshes for the solution of the compressible Navier–Stokes equations. This new family of numerical methods has first been proposed in [16] for purely hyperbolic systems and has been called P N P M schemes, where N indicates the polynomial degree of the test functions and M is the degree of the polynomials used for flux and source computation. A particular feature of the general P N P M schemes is that they contain classical high order accurate finite volume schemes ( N = 0 ) as well as standard discontinuous Galerkin methods ( M = N ) just as special cases, which therefore allows for a direct efficiency comparison. In the application section of this paper we first show numerical convergence results on unstructured meshes obtained for the compressible Navier–Stokes equations with Sutherland’s viscosity law, comparing all third to sixth order accurate P N P M schemes with each other. In order to validate the method also in practice we show several classical steady and unsteady CFD applications, such as the laminar boundary layer flow over a flat plate at high Reynolds numbers, flow past a NACA0012 airfoil, the unsteady flows past a circular cylinder and a sphere, the unsteady flows of a compressible mixing layer in two space dimensions and finally we also show applications to supersonic flows with shock Mach numbers up to M s = 10 .

Essentially non-oscillatory Residual Distribution schemes for hyperbolic problems

The Residual Distribution (RD) schemes are an alternative to standard high order accurate finite volume schemes. They have several advantages: a better accuracy, a much more compact stencil, easy parallelization. However, they face several problems, at least for steady problems which are the only cases considered here. The solution is obtained via an iterative method. The iterative convergence must be good in order to get spatially accurate solutions, as suggested by the few theoretical results available for the RD schemes. In many cases, especially for systems, the iterative convergence is not sufficient to guaranty the theoretical accuracy. In fact, up to our knowledge, the iterative convergence is correct in only two cases: for first order monotone schemes and the (scalar) Struij’s PSI scheme which is a multidimensional upwind scheme. Up to our knowledge, the iterative convergence is poor for systems, except for the blended scheme of Deconinck et al. [Á. Csík, M. Ricchiuto, H. Deconinck, A conservative formulation of the multidimensional upwind residual distribution schemes for general nonlinear conservation laws, J. Comput. Phys. 179(2) (2002) 286–312] and Abgrall [R. Abgrall, Toward the ultimate conservative scheme: following the quest, J. Comput. Phys. 167(2) (2001) 277–315] which are also a genuinely multidimensional upwind scheme. A second drawback is that their construction relies, up to now, on a single first order scheme: the N scheme. However, it is known that standard first order finite volume schemes can be rephrased into a Residual Distribution framework. Unfortunately, the standard way of upgrading the order of accuracy to second order leads to very unsatisfactory results but clearly the construction of good schemes based on a wider class of first order schemes would be interesting. In this paper, we analyze these two problems, and show they are linked. We propose a fix and demonstrate its efficiency on several test cases that cover a wide range of applications. Our solution extends considerably the number of working RD schemes.

3D adaptive central schemes: Part I. Algorithms for assembling the dual mesh

Central schemes are frequently used for incompressible and compressible flow calculations. The present paper is the first in a forthcoming series where a new approach to a 2nd order accurate Finite Volume scheme operating on Cartesian grids is discussed. Here we start with an adaptively refined Cartesian primal grid in 3D and present a construction technique for the staggered dual grid based on L ∞ -Voronoi cells. The local refinement constellation on the primal grid leads to a finite number of uniquely defined local patterns on a primal cell. Assembling adjacent local patterns forms the dual grid. All local patterns can be analysed in advance. Later, running the numerical scheme on staggered grids, all necessary geometric information can instantly be retrieved from lookup-tables. The new scheme is compared to established ones in terms of algorithmic complexity and computational effort.

A fourth order accurate finite volume scheme for numerical simulations of turbulent Rayleigh–Bénard convection in cylindrical containers

A finite volume scheme, which is based on fourth order accurate central differences in spatial directions and on a hybrid explicit/semi-implicit time stepping scheme, was developed to solve the incompressible Navier–Stokes and energy equations on cylindrical staggered grids. This includes a new fourth order accurate discretization of the velocity and temperature fields at the singularity of the cylindrical coordinate system and a new stability condition [J. Appl. Numer. Anal. Comput. Math. 1 (2004) 315–326]. The method was applied in direct numerical simulations of turbulent Rayleigh–Bénard convection for different Rayleigh numbers Ra = 10 γ , γ = 5 , … , 8 , in wide cylinders with the aspect ratios a ≡ H / R = 0.2 and a = 0.4 (where R denotes the radius and H – the height of the cylinder). To cite this article: O. Shishkina, C. Wagner, C. R. Mecanique 333 (2005).

Large Eddy Simulation of Constant Heat Flux Turbulent Channel Flow With Property Variations: Quasi-Developed Model and Mean Flow Results

Turbulent planar channel flow has been computed for uniform wall heating and cooling fluxes strong enough to cause significant property variations using large eddy simulation. Channels with both walls either heated or cooled were considered, with wall-to-bulk temperature ratios as high as 1.5 for the heated case, and as low as 0.56 for the cooled case. An implicit, second order accurate finite volume scheme was used to solve the time dependent filtered set of equations to determine the large eddy motion, while a dynamic subgrid-scale model was used to account for the subgrid scale effects. Step-periodicity was used based on a quasi-developed assumption. The effects of strong heating and cooling on the flow were investigated and compared with the results obtained under low heating conditions.