A priori error estimates for numerical methods for scalar conservation laws. Part I: The general approach

In this paper, we construct a general theory of a priori error estimates for scalar conservation laws by suitably modifying the original Kuznetsov approximation theory. As a first application of this general technique, we show that error estimates for conservation laws can be obtained without having to use explicitly any regularity properties of the approximate solution. Thus, we obtain optimal error estimates for the Engquist-Osher scheme without using the fact (i) that the solution is uniformly bounded, (ii) that the scheme is total variation diminishing, and (iii) that the discrete semigroup associated with the scheme has the -contraction property, which guarantees an upper bound for the modulus of continuity in time of the approximate solution.

     

Evolution Galerkin methods for hyperbolic systems in two space dimensions

The subject of the paper is the analysis of three new evolution Galerkin schemes for a system of hyperbolic equations, and particularly for the wave equation system. The aim is to construct methods which take into account all of the infinitely many directions of propagation of bicharacteristics. The main idea of the evolution Galerkin methods is the following: the initial function is evolved using the characteristic cone and then projected onto a finite element space. A numerical comparison is given of the new methods with already existing methods, both those based on the use of bicharacteristics as well as commonly used finite difference and finite volume methods. We discuss the stability properties of the schemes and derive error estimates.

     

非线性抛物组非均匀网格差分解的唯一性和稳定性

１．引言 １．对一维非线性抛物组,在文献山中已构造一般非均匀网格差分格式,其中差分逼近的组合系数对不同的网格点和不同的网格层可以不同,并且运用不动点原理证明了差分解的存在性和收敛性．在非均匀网格差分格式中差分逼近的组合系数为常数的情形,文献［２］证明了具有有界二阶差商的离散向量解的存在性、唯一性和稳定性．本文将对文献［１］中构造的一般非均匀网格差分格式,证明所得到的差分解的唯一性和稳定性． 考虑如下非线性抛物组其中是未知的ｍ－维向量函数是给定的矩阵函数,ｊ（ｘ,ｔ,ｕ,ｐ）。是给定的ｍ－维向量函数…     

一种高精度的修正Hermite-ENO格式

设计一种基于三单元具有六阶精度的修正Hermite-ENO格式（CHENO）,求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散,在空间层上,使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时,利用Taylor展开给出修正Hermite插值使其提高到六阶精度,并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡;在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果,是有效可行的.     

Auxiliary Principle and Iterative Algorithms for Lions-Stampacchia Variational Inequalities

In this paper, we extend the auxiliary principle (Cohen in J. Optim. Theory Appl. 49:325–333, 1988) to study a class of Lions-Stampacchia variational inequalities in Hilbert spaces. Our method consists in approximating, in the subproblems, the nonsmooth convex function by a sequence of piecewise linear and convex functions, as in the bundle method for nonsmooth optimization. This makes the subproblems more tractable. We show the existence of a solution for this Lions-Stampacchia variational inequality and explain how to build a new iterative scheme and a new stopping criterion. This iterative scheme and criterion are different from those commonly used in the special case of nonsmooth optimization. We study also the convergence of iterative sequences generated by the algorithm. This work was supported by the National Natural Science Foundation of China (10671135), the Specialized Research Fund for the Doctoral Program of Higher Education (20060610005), the National Natural Science Foundation of Sichuan Education Department of China (07ZB068) and the Open Fund (PLN0703) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University).     

An efficient splitting technique for two‐layer shallow‐water model

We consider the numerical approximation of the weak solutions of the two‐layer shallow‐water equations. The model under consideration is made of two usual one‐layer shallow‐water model coupled by nonconservative products. Because of the nonconservative products of the system, which couple both one‐layer shallow‐water subsystems, the usual numerical methods have to consider the full model. Of course, uncoupled numerical techniques, just involving finite volume schemes for the basic shallow‐water equations, are very attractive since they are very easy to implement and they are costless. Recently, a stable layer splitting technique was introduced [Bouchut and Morales de Luna, M2AN Math Model Numer Anal 42 (2008), 683–698]. In the same spirit, we exhibit new splitting technique, which is proved to be well balanced and non‐negative preserving. The main benefit issuing from the here derived uncoupled method is the ability to correctly approximate the solution of very severe benchmarks. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1396–1423, 2015     

Superconvergent estimates of conforming finite element method for nonlinear time‐dependent Joule heating equations

In this article, we study the superconvergence analysis of conforming bilinear finite element method (FEM) for nonlinear Joule heating equations. Based on the rigorous estimates together with high accuracy analysis of this element, mean value technique and interpolation postprocessing approach, the superclose and superconvergent estimates about the related variables in H¹‐norm are derived for semidiscrete and a linearized backward Euler fully discrete schemes, which extends the results of optimal estimates obtained for conforming FEMs in the previous literature. At last, a numerical experiment is performed to verify the theoretical analysis.     

Automorphisms of Hilbert schemes of points on a generic projective K3 surface

We study automorphisms of the Hilbert scheme of n points on a generic projective K3 surface S, for any . We show that is either trivial or generated by a non‐symplectic involution and we determine numerical and divisorial conditions which allow us to distinguish between the two cases. As an application of these results we prove that, for any , there exist infinitely many admissible degrees for the polarization of the K3 surface S such that admits a non‐natural involution. This provides a generalization of the results of [7] for .     

Very-high-order weno schemes

We study weno(2r − 1) reconstruction [D.S. Balsara, C.W. Shu, Monotonicity prserving weno schemes with increasingly high-order of accuracy, J. Comput. Phys. 160 (2000) 405–452], with the mapping (wenom) procedure of the nonlinear weights [A.K. Henrick, T.D. Aslam, J.M. Powers, Mapped weighted-essentially-non-oscillatory schemes: achieving optimal order near critical points, J. Comput. Phys. 207 (2005) 542–567], which we extend up to weno17 (r=9)$(r=9)$. We find by numerical experiment that these procedures are essentially nonoscillatory without any stringent cfl limitation (cfl∈[0.8,1])$(\text{<small-caps>cfl</small-caps>}\in [0.8\text{,}1])$, for scalar hyperbolic problems (both linear and scalar conservation laws), provided that the exponent p_β$p_{\beta }$ in the definition of the Jiang–Shu [G.S. Jiang, C.W. Shu, Efficient implementation of weighted eno schemes, J. Comput. Phys. 126 (1996) 202–228] nonlinear weights be taken as p_β=r$p_{\beta }=r$, as originally proposed by Liu et al. [X.D. Liu, S. Osher, T. Chan, Weighted essentially nonoscillatory schemes, J. Comput. Phys. 115 (1994) 200–212], instead of p_β=2$p_{\beta }=2$ (this is valid both for weno and wenom reconstructions), although the optimal value of the exponent is probably p_β(r)∈[2,r]$p_{\beta }(r)\in [2\text{,}r]$. Then, we apply the family of very-high-order wenom_pβ=r${\text{<small-caps>wenom</small-caps>}}_{p_{\beta }=r}$ reconstructions to the Euler equations of gasdynamics, by combining local characteristic decomposition [A. Harten, B. Engquist, S. Osher, S.R. Chakravarthy, Uniformly high-order accurate essentially nonoscillatory schemes iii, J. Comput. Phys. 71 (1987) 231–303], with recursive-order-reduction (ror) aiming at aleviating the problems induced by the nonlinear interactions of characteristic fields within the stencil. The proposed ror algorithm, which generalizes the algorithm of Titarev and Toro [V.A. Titarev, E.F. Toro, Finite-volume weno schemes for 3-D conservation laws, J. Comput. Phys. 201 (2004) 238–260], is free of adjustable parameters, and the corresponding rorwenom_pβ=r${\text{<small-caps>rorwenom</small-caps>}}_{p_{\beta }=r}$ schemes are essentially nonoscillatory, as Δx→0$\mathrm{\Delta }x\to 0$, up to r=9$r=9$, for all of the test-cases studied. Finally, the unsplit linewise 2-D extension of the schemes is evaluated for several test-cases.