scholarly journals The convergence of $\alpha$ schemes for conservation laws II: Fully-discrete case

2014 ◽  
Vol 21 (2) ◽  
pp. 201-220
Author(s):  
Nan Jiang
Author(s):  
Elias S. W. Shiu ◽  
Xiaoyi Xiong

AbstractFor a general fully continuous life insurance model, the variance of the loss-at-issue random variable is the expectation of the square of the discounted value of the net amount at risk at the moment of death. In 1964 Jim Hickman gave an elementary and elegant derivation of this result by the method of integration by parts. One might expect that the method of summation by parts could be used to treat the fully discrete case. However, there are two difficulties. The summation-by-parts formula involves shifting an index, making it somewhat unwieldy. In the fully discrete case, the variance of the loss-at-issue random variable is more complicated; it is the expectation of the square of the discounted value of the net amount at risk at the end of the year of death times a survival probability factor. The purpose of this note is to show that one can indeed use the method of summation by parts to find the variance of the loss-at-issue random variable for a fully discrete life insurance policy.


Acta Numerica ◽  
2003 ◽  
Vol 12 ◽  
pp. 451-512 ◽  
Author(s):  
Eitan Tadmor

We study the entropy stability of difference approximations to nonlinear hyperbolic conservation laws, and related time-dependent problems governed by additional dissipative and dispersive forcing terms. We employ a comparison principle as the main tool for entropy stability analysis, comparing the entropy production of a given scheme against properly chosen entropy-conservative schemes.To this end, we introduce general families of entropy-conservative schemes, interesting in their own right. The present treatment of such schemes extends our earlier recipe for construction of entropy-conservative schemes, introduced in Tadmor (1987b). The new families of entropy-conservative schemes offer two main advantages, namely, (i) their numerical fluxes admit an explicit, closed-form expression, and (ii) by a proper choice of their path of integration in phase space, we can distinguish between different families of waves within the same computational cell; in particular, entropy stability can be enforced on rarefactions while keeping the sharp resolution of shock discontinuities.A comparison with the numerical viscosities associated with entropy-conservative schemes provides a useful framework for the construction and analysis of entropy-stable schemes. We employ this framework for a detailed study of entropy stability for a host of first- and second-order accurate schemes. The comparison approach yields a precise characterization of the entropy stability of semi-discrete schemes for both scalar problems and systems of equations.We extend these results to fully discrete schemes. Here, spatial entropy dissipation is balanced by the entropy production due to time discretization with a suffciently small time-step, satisfying a suitable CFL condition. Finally, we revisit the question of entropy stability for fully discrete schemes using a different approach based on homotopy arguments. We prove entropy stability under optimal CFL conditions.


2021 ◽  
Vol 155 (18) ◽  
pp. 184101
Author(s):  
Nathan A. Seifert ◽  
Kirill Prozument ◽  
Michael J. Davis

2019 ◽  
Vol 53 (1) ◽  
pp. 105-144 ◽  
Author(s):  
Lingling Zhou ◽  
Yinhua Xia ◽  
Chi-Wang Shu

In this paper, we discuss the stability and error estimates of the fully discrete schemes for linear conservation laws, which consists of an arbitrary Lagrangian–Eulerian discontinuous Galerkin method in space and explicit total variation diminishing Runge–Kutta (TVD-RK) methods up to third order accuracy in time. The scaling arguments and the standard energy analysis are the key techniques used in our work. We present a rigorous proof to obtain stability for the three fully discrete schemes under suitable CFL conditions. With the help of the reference cell, the error equations are easy to establish and we derive the quasi-optimal error estimates in space and optimal convergence rates in time. For the Euler-forward scheme with piecewise constant elements, the second order TVD-RK method with piecewise linear elements and the third order TVD-RK scheme with polynomials of any order, the usual CFL condition is required, while for other cases, stronger time step restrictions are needed for the results to hold true. More precisely, the Euler-forward scheme needs τ ≤ ρh2 and the second order TVD-RK scheme needs $ \tau \le \rho {h}^{\frac{4}{3}}$ for higher order polynomials in space, where τ and h are the time and maximum space step, respectively, and ρ is a positive constant independent of τ and h.


2004 ◽  
Vol 01 (03) ◽  
pp. 493-520
Author(s):  
BRUNO DESPRÉS

We prove the convergence of numerical approximations of compressive solutions for scalar conservation laws with convex flux. This new proof of convergence is fully discrete and does not use Kuznetsov's approach. We recover the well-known rate of convergence in O(Δx½). With the same fully discrete approach, we also prove a rate of convergence in O(Δx) uniformly in time, if the initial data is a shock, or asymptotically after the compression of the initial profile. Numerical experiments confirm the theoretical analysis.


Sign in / Sign up

Export Citation Format

Share Document