Why nonconservative schemes converge to wrong solutions: error analysis

1994/01/01 by Thomas Y. Hou, Philippe G. LeFloch · 7 citations
Engineering · Mathematics · #Computational Fluid Dynamics and Aerodynamics #Navier-Stokes equation solutions #Advanced Mathematical Physics Problems

paper · pdf · doi:10.1090/s0025-5718-1994-1201068-0

Abstract

This paper attempts to give a qualitative and quantitative description of the numerical error introduced by using finite difference schemes in nonconservative form for scalar conservation laws. We show that these schemes converge strongly in <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper L Subscript l o c Superscript 1"> <mml:semantics> <mml:msubsup> <mml:mi>L</mml:mi> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:mi>loc</mml:mi> </mml:mrow> <mml:mn>1</mml:mn> </mml:msubsup> <mml:annotation encoding="application/x-tex">L_\operatorname loc¹</mml:annotation> </mml:semantics> </mml:math> </inline-formula> norm to the solution of an inhomogeneous conservation law containing a Borel measure source term. Moreover, we analyze the properties of this Borel measure, and derive a sharp estimate for the <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper L Superscript 1"> <mml:semantics> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:msup> <mml:mi>L</mml:mi> <mml:mn>1</mml:mn> </mml:msup> </mml:mrow> <mml:annotation encoding="application/x-tex">L¹</mml:annotation> </mml:semantics> </mml:math> </inline-formula> error between the limit function given by the scheme and the correct solution. In general, the measure source term is of the order of the entropy dissipation measure associated with the scheme. In certain cases, the error can be small for short times, which makes it difficult to detect numerically. But generically, such an error will grow in time, and this would lead to a large error for large-time calculations. Finally, we show that a local correction of any high-order accurate scheme in nonconservative form is sufficient to ensure its convergence to the correct solution.

Cited by

• Numerical methods for nonconservative hyperbolic systems: a theoretical framework.
• Convergence analysis of upwind type schemes for the aggregation equation\n with pointy potential
• Godunov-type methods for free-surface shallow flows: A review
• Modelling of stellar convection. [europepmc]
• High-resolution shock-capturing numerical simulations of three-phase immiscible fluids from the unsaturated to the saturated zone. [europepmc]
• Investigating the migration of immiscible contaminant fluid flow in homogeneous and heterogeneous aquifers with high-precision numerical simulations. [europepmc]
• A Combination of Residual Distribution and the Active Flux Formulations or a New Class of Schemes That Can Combine Several Writings of the Same Hyperbolic Problem: Application to the 1D Euler Equations. [europepmc]

Related

• The finite precision computation and the nonconvergence of difference scheme
• Error Analysis of Explicit Partitioned Runge-Kutta Schemes for Conservation Laws
• Exact Total Variation Minimizers as Non-Oscillatory Limiters in High-Order Methods for Conservation Laws
• Conservative iterative methods for implicit discretizations of conservation laws
• Accurate numerical schemes for approximating initial-boundary value problems for systems of conservation laws