求解非线性对流-扩散问题的特征—差分法

近年来,已有不少文献讨论发展型方程的特征-有限元或特征-差分解法.在这一方向上,J.Douglas,Jr.及T.F Russell[1]是重要的工作,它讨论了一维对流-扩散方程C(x)?u/?t+b(x)?u/?x-?/?x(a(x)?u/?x)=f(x,t)的数值求解问题,建立了基于线性插值与基于二次插值的两种特征-差分格式,给出了误差分析,将前一种格式推广到     

求解非线性对流-扩散问题的特征—差分法

近年来,已有不少文献讨论发展型方程的特征-有限元或特征-差分解法.在这一方向上,J.Douglas,Jr.及T.F Russell[1]是重要的工作,它讨论了一维对流-扩散方程C(x)?u/?t+b(x)?u/?x-?/?x(a(x)?u/?x)=f(x,t)的数值求解问题,建立了基于线性插值与基于二次插值的两种特征-差分格式,给出了误差分析,将前一种格式推广到     

解■u/■t=i sum from p=1 to N (■~2u/■x_p~2)的绝对稳定的三层显格式

[1]给出了解 Schrodinger型方程 u_t=iu_(xx)的两个三层显格式,其稳定条件分别为.r≤1和r≤1.2071.本文对更一般的N(≥1是自然数)维方程 ?u/?t=i sum from p=1 to N (?~2u/?x_p~2) (1)建立了一个三层显格式,并证明它是绝对稳定的. 为了建立差分格式,取时间步长τ=△t,空间步长h=△x_1=△x_2=…=△x_N;并记u_(j_1j_2…j_N)~k=u(j_1△x_1,J_2△x_2,…,j_N△x_N,k△_t).     

多维抛物型方程的分支绝对稳定的显式格式

其中及R={0≤x_i≤1,j=1,2,…,p),(?)R只为区域只的边界。 对多维抛物型方程(1)的差分解法,古典显式格式的稳定性条件为r=Δt/(Δx)~2≤1/2p,十分苛刻;古典隐式格式虽是无条件稳定,却需解线性方程组。因此两者的计算量都很大,且它们的精度较低,其局部截断误差仅为O(Δt+(Δx)~2)。因此,对多维抛物型方程而言,构造显式计算、稳定性能良好且精度较高的差分格式便具有十分明显的理论意义和实用价值。本文针对上述古典显式与隐式格式所存在的问题,构造一类对任何p维空间变量的抛物型方程(1)都适用的。分支绝对稳定的显式差分格式,其局部截断误差阶为O((Δt)~2+(Δx)~2),从而避免了解线性代数方程组,大大地减少了计算工作量,且精度较高。 令Δx_k=h_k=Δx=h=1/M(k=1,2,…p)表示空间方向步长,Δt=τ=[T/N]表示时间方向步长,M、N均为正整数。 为简便计,引入下列记号     

解高阶抛物型方程的高精度显式差分格式

1 引言 1960年,Saul’ev在文中讨论了如下的高阶(2m阶)抛物型方程 μ/t=(-1)~(m 1)~2mμ/x~(2m) (1)(其中m为正整数),提出了一类含极因子α的两层差分格式。当α=0时为显式格式,其稳定性条件为,r=△t/(△x)~(2m)<1/2~(2m-1),△t,△x分别为时间及空间步长。随后,文[2],[3]利用     

非线性对流扩散问题的交替方向差分流线扩散法

1引言Peaceman,Douglas等人于1955年提出了差分格式的交替方向法。随后,Douglas,Dupont于1972年又提出了有限元格式的交替方向法[1]。其基本思想是:对两个或三个空间变量的二阶抛物型和双曲型问题,将交替方向法与Galerkin方法相结合,通过算子分裂技术,把高维问题转化为一系列低维问题,交替地沿各空间变量的方向求解。[2]、[3]和[4]给出了对更一般扩散问题(带对流项的抛物方程)的数值求解和误差分析。     

一阶三角样条插值投影算子的模

§1.引言 设区间[0,1]的分划如下: △:0=x_0相似文献   

解二维抛物型方程的高精度差分格式

本文构造了一个解二维抛物型方程的高精度三层显式差分格式,其稳定性条件为r=△t/△x2=△t/△y2≤1/4,截断误差为O(△t2+△x4).     

求解Riesz空间分数阶扩散方程的一种新的数值方法

本文利用Diethelm方法构造了一种逼近Riesz空间分数阶导数的O（△x^3-α）格式,其中1 < α < 2,△x是空间步长.进一步对一阶时间导数采用Crank-Nicolson方法离散,得到了求解Riesz空间分数阶扩散方程的一种新的有限差分格式,并用矩阵方法证明了稳定性和收敛性,其误差估计为O（△t²+△x^3-α）,其中△t为时间步长.最后,数值算例验证了差分格式的正确性和有效性.     

进化型方程的差分格式(Ⅰ)

本文对一些常系数进化型方程证明了除常微分方程和一阶双曲型方程特殊的差分格式之外不存在绝对相容、绝对稳定的显式格式.作为推论,我们指出Hadjidimos关于他的格式是绝对相容的结论是错误的,以及对于一般常微分方程组的显式格式,能用的步长一定要满足条件△t=o(R~(-0.5),其中R为方程组右函数的.Jacobi阵的谱半径.对于热传导方程,本文利用差分算子的分裂技巧构造了两类可以显式求解的格式.它们是绝对     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP ——In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

On a class of self-similar processes with stationary increments in higher order Wiener chaoses

We study a class of self-similar processes with stationary increments belonging to higher order Wiener chaoses which are similar to Hermite processes. We obtain an almost sure wavelet-like expansion of these processes. This allows us to compute the pointwise and local Hölder regularity of sample paths and to analyse their behaviour at infinity. We also provide some results on the Hausdorff dimension of the range and graphs of multidimensional anisotropic self-similar processes with stationary increments defined by multiple Wiener–Itô integrals.     

On a class of Riemann surfaces

It is considered the class of Riemann surfaces with dimT1 = 0, where T1 is a subclass of exact harmonic forms which is one of the factors in the orthogonal decomposition of the spaceΩH of harmonic forms of the surface, namely The surfaces in the class OHD and the class of planar surfaces satisfy dimT1 = 0. A.Pfluger posed the question whether there might exist other surfaces outside those two classes. Here it is shown that in the case of finite genus g, we should look for a surface S with dimT1 = 0 among the surfaces of the form Sg\K , where Sg is a closed surface of genus g and K a compact set of positive harmonic measure with perfect components and very irregular boundary.     

A REPRESENTATION FORMULA RELATED TO SCHR(O)DINGER OPERATORS

Schr(o)dinger operator is a central subject in the mathematical study of quantum mechanics.Consider the Schrodinger operator H = -△ V on R, where △ = d2/dx2 and the potential function V is real valued. In Fourier analysis, it is well-known that a square integrable function admits an expansion with exponentials as eigenfunctions of -△. A natural conjecture is that an L2 function admits a similar expansion in terms of "eigenfunctions" of H, a perturbation of the Laplacian (see [7], Ch. Ⅺ and the notes), under certain condition on V.     

CONTENTS OF VOLUME 30

正~~     

Instructions for contributions

正Applied Mathematics-A Journal of Chinese Universities,Series B(Appl.Math.J.Chinese Univ.,Ser.B)is a comprehensive applied mathematics journal jointly sponsored by Zhejiang University,China Society for Industrial and Applied Mathematics,and Springer-Verlag.It is a quarterly journal with     

Journal of Mathematical Research with Applications Guide for Authors

正Journal overview:Journal of Mathematical Research with Applications(JMRA),formerly Journal of Mathematical Research and Exposition(JMRE)created in 1981,one of the transactions of China Society for Industrial and Applied Mathematics,is a home for original research papers of the highest quality in all areas of mathematics with applications.The target audience comprises:pure and applied mathematicians,graduate students in broad fields of sciences and technology,scientists and engineers interested in mathematics.     

Staircase transportation problems with superadditive rewards and cumulative capacities

A cumulative-capacitated transportation problem is studied. The supply nodes and demand nodes are each chains. Shipments from a supply node to a demand node are possible only if the pair lies in a sublattice, or equivalently, in a staircase disjoint union of rectangles, of the product of the two chains. There are (lattice) superadditive upper bounds on the cumulative flows in all leading subrectangles of each rectangle. It is shown that there is a greatest cumulative flow formed by the natural generalization of the South-West Corner Rule that respects cumulative-flow capacities; it has maximum reward when the rewards are (lattice) superadditive; it is integer if the supplies, demands and capacities are integer; and it can be calculated myopically in linear time. The result is specialized to earlier work of Hoeffding (1940), Fréchet (1951), Lorentz (1953), Hoffman (1963) and Barnes and Hoffman (1985). Applications are given to extreme constrained bivariate distributions, optimal distribution with limited one-way product substitution and, generalizing results of Derman and Klein (1958), optimal sales with age-dependent rewards and capacities.To our friend, Philip Wolfe, with admiration and affection, on the occasion of his 65th birthday.Research was supported respectively by the IBM T.J. Watson and IBM Almaden Research Centers and is a minor revision of the IBM Research Report [6].     

Boundedness for commutators of multipliers on Hardy type spaces

In this paper, we study the commutators generalized by multipliers and a BMO function. Under some assumptions, we establish its boundedness properties from certain atomic Hardy space Hb^p（R^n） into the Lebesgue space L^p with p 〈 1.