高维抛物方程迎风区域分裂差分方法

结合迎风方法和区域分裂思想,采用一阶迎风、二阶修正迎风法逼近高维抛物方程的对流项.内边界处和子区域分别对应区域分裂显隐格式;并运用极值原理和嵌入定理给出了收敛性分析,最后给出数值试验,说明其实际意义.     

时间分数阶期权定价模型的一类有效差分方法

时间分数阶期权定价模型(时间分数阶Black-Scholes方程)数值解法的研究具有重要的理论意义和实际应用价值.对时间分数阶Black-Scholes方程构造了显-隐格式和隐-显差分格式,讨论了两类格式解的存在唯一性,稳定性和收敛性.理论分析证实,显-隐格式和隐-显格式均为无条件稳定和收敛的,两种格式具有相同的计算量.数值试验表明:显-隐和隐-显格式的计算精度与经典Crank-Nicolson(C-N)格式的计算精度相当,其计算效率(计算时间)比C-N格式提高30%.数值试验验证了理论分析,表明本文的显-隐和隐-显差分方法对求解时间分数阶期权定价模型是高效的,证实了时间分数阶Black-Scholes方程更符合实际金融市场.     

油藏数值模拟中动边值问题的迎风差分格式

考虑了三维油藏数值模拟中的动边值问题,对压力方程,给出中心差分格式;对饱和度方程给出隐式迎风差分格式及修正的迎风差分格式,并证明了格式的收敛性。数值算例与理论结果是一致的。     

对流扩散方程区域分裂特征差分方法

1引言对流扩散方程是许多物理问题的数学模型,研究其稳定的数值解法具有重要的应用价值．而标准的差分法和有限元法通常会失效,出现数值振荡．80年代,Douglas和Russel提出了特征线方法,在一定程度上克服了数值振荡,保证了数值的稳定,尤其对“对流占优”问题,更能突出特征法的优越性,并有了大量的理论成果[1,2,3]．区域分裂是一种解决大规模的科学与工程计算问题的有效方法,Dawson,Du和Dupont对热传导方程给出了非重叠区域分裂格式及分析,由于内边界的显格式,需要一定的稳定性条件Δt≤CH2;而Du等在[5]给出了抛物方程的几种区域分裂格式,对区域分裂法的     

非线性波动方程的交替显-隐差分方法

1．引言众所周知，非线性波动方程在自然科学领域有广泛的物理背景，诸如物理、化学反应方程，机械动力学方程，地球物理与大气海洋方程等．差分方法求解非线性波动方程已有研究，如［1］和IZ］就给出了非线性波动方程组的显式和隐式差分格式以及收敛性分析．虽然古典的显式差分格式易于并行计算，但是它的稳定性条件差（条件稳定）；古典的隐式差分格式稳定性条件好（绝对稳定）；但对非线性问题，一般需要线性化，然后求解一个线性代数方程组，并行计算能力差．本文正是在这样一种前题下，给出了一维问题的一种交替分段显一隐差分格式，…     

常微分方程的半隐Runge-Kutta法

针对某些非线性常微分方程,提出一种算子分裂半隐Runge-Kutta方法,对于非线性部分采用显式计算,对于刚性强的线性部分采用隐式处理.给出了格式的推导,分析了绝对稳定性,并证明了半隐二阶格式的收敛性.相比于显式Runge-Kutta法,半隐格式计算量相近,但改进了稳定性,数值结果显示了方法的合理性和有效性.最后,将算子分裂半隐Runge-Kutta方法应用于数值求解Zakharov偏微分方程组.     

关于U_t=aU_(xx)差分格式的进一步研究(英文)

对于热传导方程构造了两个高阶精度的差分格式，一个是三层七点显格式，另一个是三层九点隐格式．证明了差分格式的收敛性和稳定性，最后给出数值计算结果．     

关于Ut=aUxx差分格式的进一步研究

对于热传导方程构造了两个高阶精度的差分格式，一个是三层七点显格式．另一个是三层九点隐格式．证明了差分格式的收敛性和稳定性，最后给出数值计算结果。     

一个求解Euler方程的特殊矩阵分裂格式

§1.引言 自[1]提出矢通量分裂格式以来,在求解气动方程方面得到广泛应用。矢通量分裂格式是一种求解守恒型双曲方程组的方法,它将方程中代表质量、动量和能量的矢通量按照矢通量Jacobian矩阵正负特征值分裂为两个亚矢通量项,目的在于改进显式格式和隐式格式的计算效率和提高求解时的稳定性。在求解方法上,对于二维问题,需要求解以4×4块矩阵为矩阵元的上三角矩阵和下三角矩阵,比中心差分格式需要求解两个块三     

时间分数阶慢扩散方程的一类有效差分方法

对时间分数阶慢扩散方程提出一类数值差分方法:显-隐(Explicit-Implicit, E-I)和隐-显(Implicit-Explicit, I-E)差分方法.它是将古典显式格式与古典隐式格式相结合构造出的一类有效差分格式.理论证明了格式解的存在唯一性,用傅里叶方法证明了格式的稳定性和收敛性.数值试验验证了理论分析,表明E-I格式和I-E格式在具有良好的精度且无条件稳定的情况下,计算速度比隐式格式提高了75%.从而用此格式解决分数阶慢扩散方程是可行的.     

Upwind finite difference schemes for linear conservation law with memory

In this article we consider upwind finite difference schemes for a class of linear conservation laws with memory. Assuming the positivity of the kernel, it is proved by using the energy estimates that the upwind finite difference scheme, explicit or implicit, is stable and convergent to the real solution. The numerical results of some examples, including Burger's equation with memory, are reported; the effect of memory is also discussed based on the numerical results. © 1994 John Wiley & Sons, Inc.     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

Discretization of a convection-diffusion equation

The unsteady convection-diffusion equation with constant coefficientsadmits an exact solution in the form of a convolution integral,which provides an explicit representation of the evolution operatorthrough one time step. This is used to unify many numericalschemes for the equation, showing interrelationships betweenfinite difference and finite element schemes and presentinga general framework for detailed error analysis. In particular,the upwind scheme, Lax Wendroff, QUICKEST, the ECG schemes andCrank Nicolson are all members of a family that includes powerfulnew schemes. Fourier analysis is used to obtain practical stabilityregions and some insights into accuracy; and the Peano kerneltheorem is also used to derive rigorous error bounds that canbe generalized to irregular meshes.     

Domain Decomposition–Finite Difference Approximate Inverse Preconditioned Schemes for Solving Fourth-Order Equations

A new class of explicit approximate inverse preconditioning is introduced for solving fourth-order equations, based on the coupled equation approach, by the domain decomposition method in conjunction with various finite difference approximation schemes. Explicit approximate inverse arrow-type matrix techniques, based on the concept of sparse L U-type factorization procedures, are introduced for computing a class of approximate inverses. Explicit preconditioned conjugate gradient-type schemes are presented for the efficient solution of linear systems. Applications of the method to a biharmonic problem are discussed and numerical results are given.     

抛物型方程的一种高精度区域分解有限差分算法

1引言 近年来,区域分解算法以可以将大型问题分解为一系列小型问题以减少计算规模及算法可高度并行实现等特点受到了人们的广泛关注.前人也做了很多很好的工作:参考文献[1]中C.N.Dawson等人提出了显一隐格式的区域分解算法,在时间层不分层的内边界点采用大步长向前-中心差分显格式及在内点采用古典隐格式,取得的精度为O(△t+h2+H3).参考文献[2]中给出了[1]中区域分解算法对于内边界点为等距分布的多子区域时的新的误差估计,使含H3误差项的系数比[1]中缩小了一倍.还将采用大步长日的saul'yev的非对称差分格式应用于内边界点,并给出了两个子区域和多个子区域情形下差分解的先验误差估计.     

THE UPWIND OPERATOR SPLITTING FINITE DIFFERENCE METHOD FOR COMPRESSIBLE TWO-PHASE DISPLACEMENT PROBLEM AND ANALYSIS

For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation rule of difference operators, decompo-sition of high order difference operators and prior estimates are adopted. Optimal order estimates in L2 norm are derived to determine the error in the approximate solution.     

The upwind finite difference fractional steps method for combinatorial system of dynamics of fluids in porous media and its application

For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L ² norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.     

激波捕捉差分方法研究

在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试．该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式．用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试．通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率．