Error estimate for finite volume scheme

We study the convergence of a finite volume scheme for the linear advection equation with a Lipschitz divergence-free speed in R ^d . We prove a h ^1/2-error estimate in the L ^∞(0,t;L ¹)-norm for BV data. This result was expected from numerical experiments and is optimal.     

Convergence of higher order upwind finite volume schemes on unstructured grids for scalar conservation laws in several space dimensions

Summary. We prove convergence of a class of higher order upwind finite volume schemes on unstructured grids for scalar conservation laws in several space dimensions. The result is applied to the discontinuous Galerkin method due to Cockburn, Hou and Shu. Received April 15, 1993 / Revised version received March 13, 1995     

A posteriori error estimate for finite volume approximations to singularly perturbed nonlinear convection-diffusion equations

Summary. This paper is devoted to the study of a posteriori and a priori error estimates for the scalar nonlinear convection diffusion equation . The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the -norm in the situation, where the diffusion parameter is smaller or comparable to the mesh size. Numerical experiments underline the theoretical results. Received February 25, 1999 / Revised version received July 6, 1999 / Published online August 2, 2000     

A combined finite volume–nonconforming/mixed-hybrid finite element scheme for degenerate parabolic problems

We propose and analyze a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the diffusion term, which generally involves an inhomogeneous and anisotropic diffusion tensor, over an unstructured simplicial mesh of the space domain by means of the piecewise linear nonconforming (Crouzeix–Raviart) finite element method, or using the stiffness matrix of the hybridization of the lowest-order Raviart–Thomas mixed finite element method. The other terms are discretized by means of a cell-centered finite volume scheme on a dual mesh, where the dual volumes are constructed around the sides of the original mesh. Checking the local Péclet number, we set up the exact necessary amount of upstream weighting to avoid spurious oscillations in the convection-dominated case. This technique also ensures the validity of the discrete maximum principle under some conditions on the mesh and the diffusion tensor. We prove the convergence of the scheme, only supposing the shape regularity condition for the original mesh. We use a priori estimates and the Kolmogorov relative compactness theorem for this purpose. The proposed scheme is robust, only 5-point (7-point in space dimension three), locally conservative, efficient, and stable, which is confirmed by numerical experiments.This work was supported by the GdR MoMaS, CNRS-2439, ANDRA, BRGM, CEA, EdF, France.     

A two-grid method for finite volume element approximations of second-order nonlinear hyperbolic equations

The two-grid method is studied for solving a two-dimensional second-order nonlinear hyperbolic equation using finite volume element method. The method is based on two different finite element spaces defined on one coarse grid with grid size H and one fine grid with grid size h, respectively. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. A prior error estimate in the H¹-norm is proved to be O(h+H³|lnH|) for the two-grid semidiscrete finite volume element method. With these proposed techniques, solving such a large class of second-order nonlinear hyperbolic equations will not be much more difficult than solving one single linearized equation. Finally, a numerical example is presented to validate the usefulness and efficiency of the method.     

Error estimates for a mixed finite element discretization of some degenerate parabolic equations

We consider a numerical scheme for a class of degenerate parabolic equations, including both slow and fast diffusion cases. A particular example in this sense is the Richards equation modeling the flow in porous media. The numerical scheme is based on the mixed finite element method (MFEM) in space, and is of one step implicit in time. The lowest order Raviart–Thomas elements are used. Here we extend the results in Radu et al. (SIAM J Numer Anal 42:1452–1478, 2004), Schneid et al. (Numer Math 98:353–370, 2004) to a more general framework, by allowing for both types of degeneracies. We derive error estimates in terms of the discretization parameters and show the convergence of the scheme. The features of the MFEM, especially of the lowest order Raviart–Thomas elements, are now fully exploited in the proof of the convergence. The paper is concluded by numerical examples.     

Second order blended multidimensional upwind residual distribution scheme for steady and unsteady computations

Multidimensional upwind residual distribution (RD) schemes have become an appealing alternative to more widespread finite volume and finite element methods (FEM) for solving compressible fluid flows. The RD approach allows to construct nonlinear second order and non-oscillatory methods at the same time. They are routinely used for steady state calculations of the complex flow problems, e.g., 3D turbulent transonic industrial-type simulations [H. Deconinck, K. Sermeus, R. Abgrall, Status of multidimensional upwind residual distribution schemes and applications in aeronautics, AAIA Paper 2000-2328, AIAA, 2000; K. Sermeus, H. Deconinck, Drag prediction validation of a multi-dimensional upwind solver, CFD-based aircraft drag prediction and reduction, VKI Lecture Series 2003-02, Von Karman Institute for Fluid Dynamics, Chausée do Waterloo 72, B-1640 Rhode Saint Genèse, Belgium, 2003].     

Two-grid methods for finite volume element approximations of nonlinear parabolic equations

Two-grid methods are studied for solving a two dimensional nonlinear parabolic equation using finite volume element method. The methods are based on one coarse-grid space and one fine-grid space. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine-grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy h=O(H³|lnH|)$h=O\left(H^3|lnH|\right)$. As a result, solving such a large class of nonlinear parabolic equations will not be much more difficult than solving one single linearized equation.     

Finite volume schemes for nonhomogeneous scalar conservation laws: error estimate

Summary. In this paper, we study finite volume schemes for the nonhomogeneous scalar conservation law with initial condition . The source term may be either stiff or nonstiff. In both cases, we prove error estimates between the approximate solution given by a finite volume scheme (the scheme is totally explicit in the nonstiff case, semi-implicit in the stiff case) and the entropy solution. The order of these estimates is in space-time -norm (h denotes the size of the mesh). Furthermore, the error estimate does not depend on the stiffness of the source term in the stiff case. Received October 21, 1999 / Published online February 5, 2001     

A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A finite volume scheme for the global shallow water model on the cubed-sphere mesh is proposed and studied in this paper. The new cell-centered scheme is based on Osher’s Riemann solver together with a high-order spatial reconstruction. On each patch interface of the cubed-sphere only one layer of ghost cells is needed in the scheme and the numerical flux is calculated symmetrically across the interface to ensure the numerical conservation of total mass. The discretization of the topographic term in the equation is properly modified in a well-balanced manner to suppress spurious oscillations when the bottom topography is non-smooth. Numerical results for several test cases including a steady-state nonlinear geostrophic flow and a zonal flow over an isolated mountain are provided to show the flexibility of the scheme. Some parallel implementation details as well as some performance results on a parallel supercomputer with more than one thousand processor cores are also provided.     

A finite volume scheme for nonlinear parabolic equations derived from one-dimensional local dirichlet problems

In this paper, we propose a new method to compute the numerical flux of a finite volume scheme, used for the approximation of the solution of the nonlinear partial differential equation u_t+div(qf(u))−ΔΦ(u)=0 in a 1D, 2D or 3D domain. The function Φ is supposed to be strictly increasing, but some values s such that Φ′(s)=0 can exist. The method is based on the solution, at each interface between two control volumes, of the nonlinear elliptic two point boundary value problem (qf(υ)+(Φ(υ))′)′=0 with Dirichlet boundary conditions given by the values of the discrete approximation in both control volumes. We prove the existence of a solution to this two point boundary value problem. We show that the expression for the numerical flux can be yielded without referring to this solution. Furthermore, we prove that the so designed finite volume scheme has the expected stability properties and that its solution converges to the weak solution of the continuous problem. Numerical results show the increase of accuracy due to the use of this scheme, compared to some other schemes.     

Error and extrapolation of a compact LOD method for parabolic differential equations

This paper is concerned with a compact locally one-dimensional (LOD) finite difference method for solving two-dimensional nonhomogeneous parabolic differential equations. An explicit error estimate for the finite difference solution is given in the discrete infinity norm. It is shown that the method has the accuracy of the second-order in time and the fourth-order in space with respect to the discrete infinity norm. A Richardson extrapolation algorithm is developed to make the final computed solution fourth-order accurate in both time and space when the time step equals the spatial mesh size. Numerical results demonstrate the accuracy and the high efficiency of the extrapolation algorithm.     

Error estimates for a mixed finite volume method for the p-Laplacian problem

In this work we propose and analyze a mixed finite volume method for the p-Laplacian problem which is based on the lowest order Raviart–Thomas element for the vector variable and the P1 nonconforming element for the scalar variable. It is shown that this method can be reduced to a P1 nonconforming finite element method for the scalar variable only. One can then recover the vector approximation from the computed scalar approximation in a virtually cost-free manner. Optimal a priori error estimates are proved for both approximations by the quasi-norm techniques. We also derive an implicit error estimator of Bank–Weiser type which is based on the local Neumann problems.This work was supported by the Post-doctoral Fellowship Program of Korea Science & Engineering Foundation (KOSEF).     

Error estimate of a first-order time discretization scheme for the geodynamo equations

In this paper, we analyze a first-order time discretization scheme for a nonlinear geodynamo model and carry out the convergence analysis of this numerical scheme. It is concluded that our numerical scheme converges with first-order accuracy in the sense of L²$L^2$-norm with respect to the velocity field u$\mathit{u}$ and the magnetic field B$\mathit{B}$ and with half-order accuracy in time for the total kinematic pressure P.     

Error bounds for the large time step Glimm scheme applied to scalar conservation laws

Summary. In this paper we derive an error bound for the large time step, i.e. large Courant number, version of the Glimm scheme when used for the approximation of solutions to a genuinely nonlinear, i.e. convex, scalar conservation law for a generic class of piecewise constant data. We show that the error is bounded by for Courant numbers up to 1. The order of the error is the same as that given by Hoff and Smoller [5] in 1985 for the Glimm scheme under the restriction of Courant numbers up to 1/2. Received April 10, 2000 / Revised version received January 16, 2001 / Published online September 19, 2001     

Convergence of a finite volume scheme for nonlinear degenerate parabolic equations

Summary. One approximates the entropy weak solution u of a nonlinear parabolic degenerate equation by a piecewise constant function using a discretization in space and time and a finite volume scheme. The convergence of to u is shown as the size of the space and time steps tend to zero. In a first step, estimates on are used to prove the convergence, up to a subsequence, of to a measure valued entropy solution (called here an entropy process solution). A result of uniqueness of the entropy process solution is proved, yielding the strong convergence of to{\it u}. Some on a model equation are shown. Received September 27, 2000 / Published online October 17, 2001     

Superconvergent biquadratic finite volume element method for two-dimensional Poisson’s equations

In this paper, a kind of biquadratic finite volume element method is presented for two-dimensional Poisson’s equations by restricting the optimal stress points of biquadratic interpolation as the vertices of control volumes. The method can be effectively implemented by alternating direction technique. It is proved that the method has optimal energy norm error estimates. The superconvergence of numerical gradients at optimal stress points is discussed and it is proved that the method has also superconvergence displacement at nodal points by a modified dual argument technique. Finally, a numerical example verifies the theoretical results and illustrates the effectiveness of the method.     

Two-grid finite volume element method for linear and nonlinear elliptic problems

Two-grid finite volume element discretization techniques, based on two linear conforming finite element spaces on one coarse and one fine grid, are presented for the two-dimensional second-order non-selfadjoint and indefinite linear elliptic problems and the two-dimensional second-order nonlinear elliptic problems. With the proposed techniques, solving the non-selfadjoint and indefinite elliptic problem on the fine space is reduced into solving a symmetric and positive definite elliptic problem on the fine space and solving the non-selfadjoint and indefinite elliptic problem on a much smaller space; solving a nonlinear elliptic problem on the fine space is reduced into solving a linear problem on the fine space and solving the nonlinear elliptic problem on a much smaller space. Convergence estimates are derived to justify the efficiency of the proposed two-grid algorithms. A set of numerical examples are presented to confirm the estimates. The work is supported by the National Natural Science Foundation of China (Grant No: 10601045).