Moving-Water Equilibria Preserving HLL-Type Schemes for the Shallow Water Equations

We construct new HLL-type moving-water equilibria preserving upwind schemes for the one-dimensional Saint-Venant system of shallow water equations with nonflat bottom topography. The designed first-and secondorder schemes are tested on a number of numerical examples, in which we verify the well-balanced property as well as the ability of the proposed schemes to accurately capture small perturbations of moving-water steady states.     

带强迫项的高阶中立型微分方程正解的存在性

考虑一类具有强迫项的高阶中立型微分方程,通过Banach压缩映像原理,分别得到了方程存在满足(?) x(t)>0的正解x(t)的充分条件与必要条件,推广了文献中的相关结果.     

变形浅水波方程组的一个非线性变换及其应用

本文用齐次平衡方法导出了变形浅水波方程组的一个非线性变换,借助该变换得到了方程组的精确孤立波解。     

浅水方程组的数值方法

构造了浅水方程组的二阶精度的TVD格式。格式由简单的TVD Runge-Kutta型时间离散和有坡度限制的空间对称离散格式组成。数值耗散项用局部棱柱化河道流的特征变量构造。格式的主要优点是能够计算天然河道中浅水方程组的弱解并且构造简单。格式能够求出天然河道或非平底部渠道中的精确静水解。给出了渠道溃坝问题数值解与解析解的比较,验证格式精度高。实际天然河道型梯级水库溃坝的数值实验表明格式稳定,适应性强。     

解四阶抛物型方程高精度紧致差分格式

针对四阶抛物型方程周期初值问题,提出了一个两层隐式差分格式和一个三层隐式差分格式.它们的局部截断误差分别为O((Δt)2+(Δx)4)和O((Δt)2+(Δt)(Δx)2+(Δx)4),其中Δt,Δx分别为时间步长和空间步长.误差分析和数值实验均表明,本文构造的差分格式比经典的Crank-Nicolson格式和Saul’ev构造的差分格式精度更高.从精度及稳定性方面考虑,本文构造的格式也比文[5]的显式格式要好.     

一种求解Burgers方程的高精度紧致差分格式

针对Burgers方程,采用余项修正法和欧拉公式,推导了一种新的四层高精度紧致差分隐格式,其截断误差为O(τ~2+τh~2+h~4),即当τ=O(h~2)时,格式空间具有四阶精度;然后通过数值实验验证了格式的精确性和可靠性.     

A high order central-upwind scheme for hyperbolic conservation laws

A high order central-upwind scheme for approximating hyperbolic conservation laws is proposed. This construction is based on the evaluation of the local propagation speeds of the discontinuities and Peer's fourth order non-oscillatory reconstruction. The presented scheme shares the simplicity of central schemes, namely no Riemann solvers are involved. Furthermore, it avoids alternating between two staggered grids, which is particularly a challenge for problems which involve complex geometries and boundary conditions. Numerical experiments demonstrate the high resolution and non-oscillatory properties of our scheme.     

双曲守恒律的几种新数值方法的比较研究

本文就一维线性双曲方程的光滑和间断两种初值问题的求解,对双曲守恒律的三种新数值方法,即,WENO方法、间断Galerkin方法和全局复合方法,进行了数值比较实验,在精度、计算速度等方面的比较上,对这三个方法有了一个较详细的了解,得到了一些有用的结论。     

色散方程两组三阶恒稳显格式

采用组合差商法对色散方程ut=auxxx(a为常数)的初边值问题,构造了两组互为对称带参数的三层显式差分格式.它们空间宽度为4,其局部截断误差为O(τ+h3),绝对稳定.而且计算时无方向性的约束,即不管a的符号如何,每一组格式均可以计算.最后给出了数值例子,数值结果表明了理论分析的正确性.     

Finding Periodic Solutions of High Order Duffing Equations via Homotopy Method

This paper presents a detailed analysis of finding the periodic solutions for the high order Duffing equation x^（2n） ＋ g（x） = e（t） （n ≥ 1）. Firstly, we give a constructive proof for the existence of periodic solutions via the homotopy method. Then we establish an efficient and global convergence method to find periodic solutions numerically.     

A New Discontinuous Galerkin Method for Parabolic Equations with Discontinuous Coefficient

In this paper, a new discontinuous Galerkin method is developed for the parabolic equation with jump coefficients satisfying the continuous flow condition. Theoretical analysis shows that this method is $L^2$ stable. When the finite element space consists of interpolative polynomials of degrees $k$, the convergent rate of the semi-discrete discontinuous Galerkin scheme has an order of$\mathcal{O}(h^k)$. Numerical examples for both 1-dimensional and 2-dimensional problems demonstrate the validity of the new method.     

抛物型方程的一个新的高精度恒稳定的隐式差分格式

本文用待定参数法对一维抛物型方程构造出一个截断误差为 0 (△ t3+△ x6)的隐式差分格式 ,格式绝对稳定且可用追赶法求解 .     

A Finite Volume Implicit Euler Scheme for the Linearized Shallow Water Equations: Stability and Convergence

The linearized shallow water equations are discretized in space by a finite volume method and in time by an implicit Euler scheme. Stability and convergence of the scheme are proved.     

Local Discontinuous Galerkin Methods for Three Classes of Nonlinear Wave Equations

In this paper, we further develop the local discontinuous Galerkin method to solve three classes of nonlinear wave equations formulated by the general KdV-Burgers type equations, the general fifth-order KdV type equations and the fully nonlinear K(n, n, n) equations, and prove their stability for these general classes of nonlinear equations. The schemes we present extend the previous work of Yan and Shu [30, 31] and of Levy, Shu and Yan [24] on local discontinuous Galerkin method solving partial differential equations with higher spatial derivatives. Numerical examples for nonlinear problems are shown to illustrate the accuracy and capability of the methods. The numerical experiments include stationary solitons, soliton interactions and oscillatory solitary wave solutions.The numerical experiments also include the compacton solutions of a generalized fifthorder KdV equation in which the highest order derivative term is nonlinear and the fully nonlinear K (n, n, n) equations.     

Development and Comparison of Numerical Fluxes for LWDG Methods

The discontinuous Galerkin （DO） or local discontinuous Galerkin （LDG） method is a spatial discretization procedure for convection-diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters. The Lax- Wendroff time discretization procedure is an altemative method for time discretization to the popular total variation diminishing （TVD） Runge-Kutta time discretizations. In this paper, we develop fluxes for the method of DG with Lax-Wendroff time discretization procedure （LWDG） based on different numerical fluxes for finite volume or finite difference schemes, including the first-order monotone fluxes such as the Lax-Friedfichs flux, Godunov flux, the Engquist-Osher flux etc. and the second-order TVD fluxes. We systematically investigate the performance of the LWDG methods based on these different numerical fluxes for convection terms with the objective of obtaining better performance by choosing suitable numerical fluxes. The detailed numerical study is mainly performed for the one-dimensional system case, addressing the issues of CPU cost, accuracy, non-oscillatory property, and resolution of discontinuities. Numerical tests are also performed for two dimensional systems.     

流动控制方程中对流项离散格式的比对分析

首先在有限体积法的基础上,针对流体流动控制方程中一阶对流项的离散问题,通过选用不同的控制节点来产生、分析已有的插值函数,从而形成不同的离散格式;其次,通过应用一定的数值算例来对各离散格式进行了相应的数值比对、分析,得出了影响问题求解的一些因素,选取出了一种相对比较稳定、高效的对流项离散格式.     

求解一维非定常对流扩散反应方程的高精度紧致差分格式

提出了数值求解一维非定常对流扩散反应方程的一种高精度紧致隐式差分格式,其截断误差为O(τ~4+τ~2h~2+h~4),即格式整体具有四阶精度.差分方程在每一时间层上只用到了三个网格节点,所形成的代数方程组为三对角型,可采用追赶法进行求解,最后通过数值算例验证了格式的精确性和可靠性.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.