AN ADAPTIVE VERSION OF GLIMM＇S SCHEME

This article describes a local error estimator for Glimm＇s scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm＇s scheme by an optimal choice. As a by-product of the local error estimator, the procedure provides a global error estimator that is shown numerically to be a very accurate estimate of the error in L1 （R） for all times. Although there is partial mathematical evidence for the error estimator proposed, at this stage the error estimator must be considered ad- hoc. Nonetheless, the error estimator is simple to compute, relatively inexpensive, without adjustable parameters and at least as accurate as other existing error estimators. Numerical experiments in 1-D for Burgers＇ equation and for Euler＇s system are performed to measure the asymptotic accuracy of the resulting scheme and of the error estimator.     

环形空腔内自然对流问题的混合有限元方法的误差估计

本文从研究环形空腔内自然对流问题的流体运动状态及温度分布规律出发,构造出Boussinesq方程组初、边值问题,给出了两种混合有限元计算格式及它们的有限元误差估计.     

A posteriori error estimates for upwind finite volume schemes for nonlinear conservation laws in multi dimensions

In this paper we shall derive a posteriori error estimates in the -norm for upwind finite volume schemes for the discretization of nonlinear conservation laws on unstructured grids in multi dimensions. This result is mainly based on some fundamental a priori error estimates published in a recent paper by C. Chainais-Hillairet. The theoretical results are confirmed by numerical experiments.

     

非结构网格上一类满足局部极值原理的三阶精度有限体积方法

对二维标量双曲型守恒律方程,发展了一类满足局部极值原理的非结构网格有限体积格式.其构造思想是,以单调数值通量为基础,通过应用基于最小二乘法的二次重构和极值限制器,使数值解满足局部极值原理.为保证数值解在光滑区域达到三阶精度,该格式可结合局部光滑探测器使用.本文从理论上分析了格式的稳定性条件,数值实验验证了格式的精度和对间断的分辨能力.     

AN ADAPTIVE VARIANT OF CGNR ALGORITHM

1　IntroductionConsider large linear systemAx =b, ( 1 .1 )where A and b are given,A isa nonsingularn×n real matrix,and b is an n-vector.When A issymmetric positive definite ( s.p.d.) ,the classical conjugate gradient method ( CG) [7] used inconjunction with incomplete factorization preconditioning technique[1 0 ] is a succesful method inview of the high robustness and efficiency.But for general non-s.p.d.systems,the situationis less satisfactory.The generalized minimal residual algorithm( …     

An adaptive wavelet viscosity method for hyperbolic conservation laws

We extend the multiscale finite element viscosity method for hyperbolic conservation laws developed in terms of hierarchical finite element bases to a (pre‐orthogonal spline‐)wavelet basis. Depending on an appropriate error criterion, the multiscale framework allows for a controlled adaptive resolution of discontinuities of the solution. The nonlinearity in the weak form is treated by solving a least‐squares data fitting problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

AN EXPLICIT MULTI-CONSERVATION FINITE-DIFFERENCE SCHEME FOR SHALLOW-WATER-WAVE EQUATION

An explicit multi-conservation finite-difference scheme for solving the spherical shallow-water-wave equation set of barotropic atmosphere has been proposed. The numerical scheme is based on a special semi-discrete form of the equations that conserves four basic physical integrals including the total energy, total mass, total potential vorticity and total enstrophy. Numerical tests show that the new scheme performs closely like but is much more time-saving than the implicit multi-conservation scheme.     

铁磁链Landau-Lifshitz方程的显式差分法

正如在研究流体动力学时,Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程起着十分重要的作用一样,在对于非平衡态磁学的研究中,描述连续铁磁体自旋场发展过程的 Ｌａｎｄａｎ－Ｌｉｆｓｈｉｔｚ方程［１］起着十分重要的作用［２］．一九九三年,美国和印度签署了一个大约 ２８０万美元的合作研究计划,在三年的时间里,对Ｌａｎｄａｕ－Ｌｉｆｓｈｉｔｚ铁磁链方程进行研究．在无阻尼的情况下,它为一完全可积的孤立子系统［３,４,５］。很多物理学家研究了它的孤立子解的存在性、逆散射方法以及相互碰撞［３,４,５］．关于解的存在性, Ａｌｏｎ…     

Klein-Gordon方程初达值问题的一个新的守恒差分格式

本文对非线性Klein－Gordon（NKG）方程的初边值问题提出了一种新的差分格式,它保持了NKG方程初边值问题的能量守恒．证明了该格式的收敛性和稳定性．特别地,由于该格式是完全隐式的,故对求长时解有着重要的作用．数值计算结果表明该方法计算速度快,精度好．     

非线性对流扩散方程沿特征线的多步有限体积元格式

对于二维非线性对流扩散方程构造了沿特征线方向的多步有限体积元格式．关于空间采用二次有限体积元方法离散,关于时间采用多步法进行离散,获得了O(Δt^2 h^2)形式的误差估计．本文最后给出的数值算例表明了方法的有效性．     

特征值问题的自适应反迭代有限元算法

1引言设Ω∈R~2为Lipschitz单连通的有界闭区域,X为定义在Ω的Sobolev空间,a(·,·)和b(·,·)为X×X→C的有界双线性或半双线性泛函,考虑变分特征值问题:求(λ,u≠0)∈C×X使得a(u,v)=λb(u,u),(?)u∈X,其中a(·,·)满足X上的"V-强制性"条件或者连续的inf-sup条件,设M_h为Q区域上的正则三角形剖分,X_h∈X为定义在M_h有限元子空间,上述变分问题对应的有限元离散问题为:求(λ_h,u_h)∈R×X,u_h≠0使得     

ADAPTIVE FINITE ELEMENT APPROXIMATION FOR A CLASS OF PARAMETER ESTIMATION PROBLEMS

In this paper, we study adaptive finite element discretisation schemes for a class of parameter estimation problem. We propose to efficient algorithms for the estimation problem use adaptive multi-meshes in developing We derive equivalent a posteriori error estimators for both the state and the control approximation, which particularly suit an adaptive multi-mesh finite element scheme. The error estimators are then implemented and tested with promising numerical results.     

对流扩散方程迎风有限元的自适应方法

本文对二维发展型对流扩散方程的迎风有限元格式给出了显式后验误差估计,证明了真实误差被后验误差估计器上下界定;并通过误差估计器建立了相应的自适应算法,数值例子表明了方法的有效性.     

自适应迎风格式求解奇异摄动两点边值问题的高精度算法(二)

0引言 考虑与文[1]相同的奇异摄动两点边值问题的数值解法: Tu(x):=-εu″(x)-p(x)u′(x)=f(x),x∈(0,1); (1) u(0)=0,u(1)=1. (2) 其中ε是一个常数,0<ε≤1,f∈C2[0,1].假定P∈C3[0,1]且存在常数β和-β使得0<β≤p(x)≤-β,|p′(x)|≤-β,(V)x∈[0,1] (3) 成立.     

一个求解多维守恒律方程组的二阶显式有限元格式

１．引言 近年来,在非结构网格上求解双曲型守恒律的数值方法引起了较为广泛的关注,出现了有限体积方法［１］,间断 Ｇａｌｅｒｋｉｎ方法 ［２］,流线扩散方法［３］,以及 ＮＮＤ格式 ［４］等．我们在［６,７］中提出了一种求解双曲型守恒律方程式的有限元方法,它是在一个求解对流扩散问题的有限元方法 ［５］的基础上发展起来的．它是一个显式有限元方法,因此计算量很小．在这个方法中,我们将任意维的问题归结为在单元棱边上的一维计算,引入了积分因子,因此在单元内部可以容纳边界层．这样,它特别适合于对流占优问题以及双曲…     

On the optimization of flux limiter schemes for hyperbolic conservation laws

A classic strategy to obtain high‐quality discretizations of hyperbolic partial differential equations is to use flux limiter (FL) functions for blending two types of approximations: a monotone first‐order scheme that deals with discontinuous solution features and a higher order method for approximating smooth solution parts. In this article, we study a new approach to FL methods. Relying on a classification of input data with respect to smoothness, we associate specific basis functions with the individual smoothness notions. Then, we construct a limiter as a linear combination of the members of parameter‐dependent families of basis functions, and we explore the possibility to optimize the parameters in interesting model situations to find a corresponding optimal limiter. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

具有跳扩散的美式期权二叉树计算格式的收敛速率

American put option with jump-diffusion can be modelled as a vari- ational inequality problem with an integral term.Under the stability condition (σ~2Δt)/(Δx~2)≤1,whereΔx=ln(S_n 1)/(S_n),the convergence rate O((Δx)~(2/3) (Δt)~(1/3))of the explicit finite scheme for this problem is obtained by using penalization technique. The binomial tree scheme of this model,which is equivalent to the explicit scheme, is convergent by the same rate.     

非线性Gerdjikov-Ivanov方程的守恒差分格式

本文考察了一类非线性Gerdjikov-Ivanov方程的周期初值问题，提出了一种守恒的差分格式，对其差分解作了先验估计．证明了格式的收敛性与稳定性，最后，通过数值计算检验了格式的可信性。     

用基于水平集方法的自适应运动网格变形方法解Hamilton-Jacobi方程组

A numerical scheme is presented for solving the Hamilton-Jacobi equation by applying adaptive moving grid methods of level-set-based deformation methods. Two numerical examples are given, which demonstrate the accuracy and efficiency of computing“extreme”and“spikes”of solutions to the Hamilton-Jacobi equation.