一类反应扩散方程组的保持正性的差分格式

本文研究一类具有正解的反应扩散方程组的有限差分解法.构造了一个保持正性的差分格式.利用离散的最大值原理证明了差分格式解的非负性,有界性及差分格式的无条件稳定性.这些估计的证明不依赖于微分方程的解而仅仅与初边值条件有关.当微分方程的解适当光滑时,证明了差分格式的一致收敛性.最后给出了数值计算结果,并与以往方法进行了比较.计算结果说明了本文给出的方法的有效性.     

流体力学方程组的总熵增量小的守恒型差分格式

1.引言 近年来,国外许多学者对求解双曲守恒律组的高分辨率、高精度差分格式进行了深入的研究。例如MUSCL方法、TVD格式、PPM方法、各种限流的方法以及ENO格式等等。将这些方法应用于流体力学方程组,其数值实践的结果表明,在消除波后振荡、提高激波间断分辨率、提高计算精度等方面有明显的效果。在设计这些     

中心格式在Saint-Venant方程组上的应用研究

浅水波方程组对于其数值格式有较高的要求.在实际应用中作者更关心在稳态解附近的行为,特别是当计算区域出现干湿界面的时候,不但要求格式具有和谐性,而且需要保持水深恒为非负,同时又要求数值格式具有较高的精度.设计同时满足这些性质的数值格式具有一定的难度.论文的核心是总结研究了受关注的求解浅水波方程组的中心格式:KP格式,BCKN格式和T格式的各自优势以及不足之处.该文通过求解一维问题来展示各自格式在一些算例上的应用.     

求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

包含极值点和分枝点的解曲线跟踪算法的摄动解法

本文提出了采用摄动格式求解非线性方程组的解曲线跟踪算法的计算格式.文中着重讨论了解曲线上非正则点的搜索,以及从这些非正则点——转向点或分枝点——继续跟踪超临界平衡路径的计算方法.文中把这一算法应用于弹性薄壳的屈曲分析。通过柱壳和环壳的算例得到它们的整个屈曲过程的平衡路径和变形形态.     

三维热传导型半导体器件瞬态模拟问题的Crank-Nicolson差分-流线扩散有限元法及其数值分析

本文研究三维热传导型半导体器件瞬态模拟问题的数值方法.针对数学模型中各方程不同的特点,分别提出不同的有限元格式.特别针对浓度方程组是对流为主扩散问题的特点,使用Crank-Nicolson差分-流线扩散计算格式,提高了数值解的稳定性.得到的L2误差估计关于空间剖分步长是拟最优的,关于时间步长具有二阶精度.     

一类强耦合的奇异摄动对流扩散方程组的移动网格方法（英文）

本文研究一类强耦合的奇异摄动对流扩散方程组的移动网格方法.首先,利用迎风有限差分格式对方程组进行离散.然后,推出数值解的后验误差估计,并以此设计出相应的自适应网格生成算法.同时,证明数值解具有一阶一致收敛性.最后,数值实验验证了本文移动网格方法的一致收敛性.     

等熵气体动力学方程组Lax-Friedrichs格式的收敛性(Ⅰ)

§1引言 如所周知,Lax-Friedrichs格式是P.D.Lax对拟线性双曲型守恒律方程组提出的一种有限差分格式。若得到了其相应的差分逼近解的收敛性,这格式不仅提供了证明:整体广义解存在性的一种理想途径,而且能方便有效地直接用来进行整体解的数值计算。在单个守恒律方程情形,O.Oleinik,C.Conway and J.Smoller等证明了这一格式的收敛性,并得到了整体广义解的存在性。然而,对双曲型方程组,特别是气体动力学方程组,Lax-Friedrichs格式的收敛性一直没有什么结果。     

一类线性隐式A稳定的单步法

1.引言 对于Stiff方程组初值问题的数值解法,Dahlquist在[1]中引进了 A稳定的概念,并且证明了显式的线性多步法(包括显式的Runge-Kutta方法)不可能是A稳定的.现在已经有许许多多隐式A稳定或Stiff稳定的方法,但绝大多数在数值解的过程中必须解由于隐式方法所产生的非线性方程组,而非线性方程组的求解过程往往又要采用Newton-Raphson迭代方法,因此需要计算方程y’=f(x,y)的右函数f(x,y)的Jacobi矩阵以及与此有关的逆矩阵.本文的主要思想是:既然在数值解过程中要计算f(x,y)的Jacobi矩阵,那么不妨在数值公式中明显的出现f(x,y)的一阶偏导数.我们将A稳定公式     

关于预报──校正方法的注记

关于预报──校正方法的注记汪先光（苏州丝绸工学院）在讨论常微分方程初值问题数值解法时。一般针对一阶方程：建立计算格式。这些格式也适用于一阶方程组。而高阶方程可以转化为一阶方程组。所以这些格式也适用于一般高阶方程。求解（１）的数值计算格式可分为二大类，...     

双曲型守恒律方程的熵稳定格式的一些讨论

本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对,格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量,但不能唯一确定方程组的熵守恒通量,却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响.     

On entropy stable schemes for degenerate parabolic multispecies kinematic flow models

Entropy stable schemes for the numerical solution of initial value problems of nonlinear, possibly strongly degenerate systems of convection–diffusion equations were recently proposed in Jerez and Parés's study. These schemes extend the theoretical framework of Tadmor's study to convection–diffusion systems. They arise from entropy conservative schemes by adding a small amount of viscosity to avoid spurious oscillations. The main condition for feasibility of entropy conservative or stable schemes for a given model is that the corresponding first‐order system of conservation laws possesses a convex entropy function and corresponding entropy flux, and that the diffusion matrix multiplied by the inverse of the Hessian of the entropy is positive semidefinite. As a new contribution, it is demonstrated in the present work, first, that these schemes can naturally be extended to initial‐boundary value problems with zero‐flux boundary conditions in one space dimension, including an explicit bound on the growth of the total entropy. Second, it is shown that these assumptions are satisfied by certain diffusively corrected multiclass kinematic flow models of arbitrary size that describe traffic flow or the settling of dispersions and emulsions, where the latter application gives rise to zero‐flux boundary conditions. Numerical examples illustrate the behavior and accuracy of entropy stable schemes for these applications.     

对流扩散方程一类改进的特征线修正有限元方法

１引言在地下水污染,地下渗流驱动,核污染,半导体等问题的数值模拟中,均涉及抛物型对流扩散方程（或方程组）的数值求解问题．这些对流扩散型偏微分方程（或方程组）具有共同的特点：对流的影响远大于扩散的影响,即对流占优性,对流占优性给问题的数值求解带来许多困难,因此对流占优问题的有效数值解法一直是计算数学中重要的研究内容．用通常的差分法或有限元法进行数值求解将出现数值振荡．为了克服数值振荡,提出各种迎风方法和修正的特征方法并在这些问题上得到成功的实际应用、８０年代,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［２］等…     

High Resolution Schemes and Discrete Entropy Conditions for 2-D Linear Conservation Laws

In this paper, fully discrete entropy conditions of a class of high resolution schemes with the MmB property are discussed by using the theory of proper discrete entropy flux for the linear scalar conservation laws in two dimensions. The theoretical resluts show that the high resolution schemes satisfying fully discrete entropy conditions with proper discrete entropy flux cannot preserve second order accuracy in the case of two dimensions.     

MUSCL Type Schemes and Discrete Entropy Conditions

1.IntroductionConsiderthehyperbolicconservationlaws:Theresearchofnumericalmethodsforequations(1.1)hasbeendevelopedrapidlyinthisdecade.SincetheappearanceoftheconceptofTVD(totalvariationdiminishing)schemes,varioushighresolutionschemeshavebeenproposedl1,2,3l4]andsuccessfullyappliedtocomputationalfluiddynamics-ItiswellknownthattheconvergenceofthenumericalmethodsforhyperbolicconservationlawsdependsontheentroPyconditionofthenumericalsolutionsl5].Previouslytheconstructionofdifferenceschemeswasalway…     

Entropy flux-splittings for hyperbolic conservation laws part I: General framework

A general framework is proposed for the derivation and analysis of flux-splittings and the corresponding flux-splitting schemes for systems of conservation laws endowed with a strictly convex entropy. The approach leads to several new properties of the existing flux-splittings and to a method for the construction of entropy flux-splittings for general situations. A large family of genuine entropy flux-splittings is derived for several significant examples: the scalar conservation laws, the p-system, and the Euler system of isentropic gas dynamics. In particular, for the isentropic Euler system, we obtain a family of splittings that satisfy the entropy inequality associated with the mechanical energy. For this system, it is proved that there exists a unique genuine entropy flux-splitting that satisfies all of the entropy inequalities, which is also the unique diagonalizable splitting. This splitting can be also derived by the so-called kinetic formulation. Simple and useful difference schemes are derived from the flux-splittings for hyperbolic systems. Such entropy flux-splitting schemes are shown to satisfy a discrete cell entropy inequality. For the diagonalizable splitting schemes, an a priori L^∞ estimate is provided by applying the principle of bounded invariant regions. The convergence of entropy flux-splitting schemes is proved for the 2 × 2 systems of conservation laws and the isentropic Euler system. ©1995 John Wiley & Sons, Inc.     

一类TVD格式的熵强迫函数及熵条件

1.引言 在拟线性双曲型方程差分方法的研究中,数值解的收敛性是一个重要的问题,因为这一性质决定了数值解能否近似地反映真实的物理现象。数值解的收敛性实际上包含了三方面的内容: 1°当网格步长趋于零时,数值解序列包含一个按某种范数收敛到函数u的子序列。     

A minimum entropy principle of high order schemes for gas dynamics equations

The entropy solutions of the compressible Euler equations satisfy a minimum principle for the specific entropy (Tadmor in Appl Numer Math 2:211–219, 1986). First order schemes such as Godunov-type and Lax-Friedrichs schemes and the second order kinetic schemes (Khobalatte and Perthame in Math Comput 62:119–131, 1994) also satisfy a discrete minimum entropy principle. In this paper, we show an extension of the positivity-preserving high order schemes for the compressible Euler equations in Zhang and Shu (J Comput Phys 229:8918–8934, 2010) and Zhang et?al. (J Scientific Comput, in press), to enforce the minimum entropy principle for high order finite volume and discontinuous Galerkin (DG) schemes.     

Discrete shock profiles for systems of conservation laws

The existence of discrete shock profiles for difference schemes approximating a system of conservation laws is the major topic studied in this paper. The basic theorem established here applies to first-order accurate difference schemes; for weak shocks, this theorem provides necessary and sufficient conditions involving the truncation error of the linearized scheme which guarantee entropy satisfying or entropy violating discrete shock profiles. Several explicit difference schemes are used as examples illustrating the interplay between the entropy condition, monotonicity, and linearized stability. Entropy violating stationary shocks for second-order accurate Lax-Wendroff schemes approximating systems are also constructed. The only tools used in the proofs are local analysis and the center manifold theorem.