几乎各向同性的高维空间分数阶扩散方程的分块快速正则Hermite分裂预处理方法

一类空间分数阶扩散方程经过有限差分离散后所得到的离散线性方程组的系数矩阵是两个对角矩阵与Toeplitz型矩阵的乘积之和.在本文中,对于几乎各向同性的二维或三维空间分数阶扩散方程的离散线性方程组,采用预处理Krylov子空间迭代方法,我们利用其系数矩阵的特殊结构和具体性质构造了一类分块快速正则Hermite分裂预处理子.通过理论分析,我们证明了所对应的预处理矩阵的特征值大部分都聚集于1的附近.数值实验也表明,这类分块快速正则Hermite分裂预处理子可以明显地加快广义极小残量(GMRES)方法和稳定化的双共轭梯度(BiCGSTAB)方法等Krylov子空间迭代方法的收敛速度.     

A NEW STABILIZED FINITE ELEMENT METHOD FOR SOLVING THE ADVECTION-DIFFUSION EQUATIONS

1. IntroductionConsider the advection--diffusion equationin a bounded polygonal domain fl c IR2 with the boundary an, where o < K << 1 is thediffusion parameter, rr > 0 is a given positive constant, g(x) is a given vector field representingthe flow with V…     

基于分数阶扩散方程的离散线性代数方程组迭代方法研究

本文构建一类双参数拟Toeplitz分裂（TQTS）迭代方法求解变系数非定常空间分数阶扩散方程.TQTS迭代法是基于QTS迭代法引入双参技术建立而成,通过选取适当的参数使迭代矩阵谱半径变得更小,从而有效提升收敛的速度.然后对TQTS迭代法进行收敛性分析,获得相应的收敛区域,并对迭代法中涉及的参数进行讨论,获得使迭代矩阵谱半径上界达到最小的最优参数的表达式.最后通过数值仿真实验验证TQTS迭代法的有效性,实验结果表明TQTS迭代法改进效果十分突出,在迭代时间和步数上均有明显的减小.     

一类分数阶多项延迟微分方程的Jacobi谱配置方法

本文利用Jacobi谱配置方法数值求解了一类分数阶多项延迟微分方程,并证明了该方法是收敛的,通过若干数值算例验证了相应的理论结果,结果表明Jacobi谱配置方法求解这类方程是非常高效的,同时也为这类分数阶延迟微分方程的数值求解提供了新的选择,对分数阶泛函方程的数值方法的研究有一定的指导意义.     

A METHOD FOR SOLVING THE INVERSE SCATTERING PROBLEM FOR SHAPE AND IMPEDANCE

The inverse problem considered in this paper is to determine the shape and the impedance of an obstacle from a knowledge of the time-harmonic incident field and the phase and amplitude of the far field pattern of the scattered wave in two-dimension. Single-layer potential is used to approach the scattered waves. An approximation method is presented and the convergence of the proposed method is established. Numerical examples are given to show that this method is both accurate and easy to use.     

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＦｉｎｉｔｅｖｏｌｕｍｅｅｌｅｍｅｎｔｍｅｔｈｏｄｏｒＦＶＥＭｕｓｅｓａｖｏｌｕｍｅｉｎｔｅｇｒａｌｆｏｒｍｕｌａｔｉｏｎｏｆｔｈｅｄｉｆ ｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｉｔｈａｆｉｎｉｔｅｐａｒｔｉｔｉｏｎｉｎｇｓｅｔｏｆｖｏｌｕｍｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｅｑｕａｔｉｏｎ ,ｔｈｅｎｒｅ ｓｔｒｉｃｔｓｔｈｅａｄｍｉｓｓｉｂｌｅｆｕｎｃｔｉｏｎｓｔｏａｌｉｎｅａｒｆｉｎｉｔｅｅｌｅｍｅｎｔｓｐａｃｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｓｏｌｕｔｉｏｎ ([1 -4 ] ) .…     

一个求解水平线性互补问题的高阶可行内点算法的多项式复杂性

最近,Zhao和Sun提出了一个求解sufficient线性互补问题的高阶不可行内点算法.不需要严格互补解条件,他们的算法获得了高阶局部收敛率,但他们的文章没有报告多项式复杂性结果.本文我们考虑他们所给算法的一个简化版本,即考虑求解单调水平线性互补问题的一个高阶可行内点算法.我们证明了算法的迭代复杂性是     

求解非光滑方程组的Newton-GMRES方法的全局收敛性分析

1引言对于非光滑方程组F(y)=0(1.1)的求解,这里F:D(?)R~n→R~n是一个非光滑映射,目前主要有两种求解方法.一种方法是由Pang提出的.     

求解流函数涡度方程的特征线—差分方法

本文讨论了用特征线方法与有限差分方法相结合的数值方法(特征线—差分方法)求解流函数涡度形式的Navier-Stokes方程的问题.证明了该方法的收敛性,给出了数值例子.     

GENUINE-OPTIMAL CIRCULANT PRECONDITIONERS FOR WIENER-HOPF EQUATIONS

1. IntroductionWienerHopf equations are integral equations defined on the haif line:where rr > 0, a(.) C L1(ro and g(.) E L2(at). Here R = (--oo,oo) and ty [0,oo). Inou-r discussions, we assume that a(.) is colljugate symmetric, i.e. a(--t) = a(t). WienerHop f equations arise in a variety of practical aPplicatiolls in mathematics and ellgineering, forinstance, in the linear prediction problems fOr stationary stochastic processes [8, pp.145--146],diffuSion problems and scattering problems […     

带有对流项的拟线性Sobolev方程的差分-流线扩散法

1　引　　言流线扩散法 (Stream line- Diffusion method,简称 SD方法 )是由 Hughes和 Brooks在文献 [1]中提出 ,并由 Johnson等人 (文献 [2 ]- [4])发展起来的求解对流占优扩散问题 (包括双曲型问题 )的一种有效的数值方法 .SD方法兼具有良好的数值稳定性和高阶收敛精度已广泛地应用于计算流体等诸多问题 .然而 ,传统的 SD方法均采用时空有限元求解发展型问题 .这样 ,虽可使时间和空间方向上的精度很好的协调起来 ,却增大了实际计算的复杂程度 ;对非线性问题也不便进行线性化处理 .孙澈等人在文献 [5 ]中对对流扩散问题提出了仅对空间…     

对流占优扩散方程的特征线楔形基无网格法

本文研究了一维对流占优扩散方程的初边值问题.利用特征线法与楔形基无网格法,获得了特征线楔形基无网格显格式与隐格式算法.数值实验表明算法具有精度高、计算简单等优点.     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

一种新的求非线性方程组的数值延拓法

针对迭代过程中的Jacobi奇异问题,本文提出了一种新的数值延拓法.通过构造双参数同伦算子,采用可控条件和适当选取参数的方式克服Jacobi奇异性,并分析了方法的收敛性.最后,通过数值实验对比,验证了方法的可行性和优越性.特别是具有可调控越过Jacobi奇异(点、线、面)的优势,从而也在某种程度上解决了数值延拓法严重依赖于初值的问题.     

解第一类边界积分方程的高精度机械求积法与外推

０．引言使用单层位势理论把Ｄｉｒｉｃｈｌｅｔ问题：转化为具有对数核的边界积分方程：这里Г假设为简单光滑闭曲线．熟知,若Г的容度Ｃｒ≠１,（０．２）有唯一解存在［１］．借助参数变换这里的数值解法有Ｇａｌｅｒｋｉｎ法［２］,配置法［３］,和谱方法~［４］,这些方法有一个共同缺点就是矩阵元素的生成要计算反常积分,由于离散方程的系数矩阵是满阵,使矩阵生成的工作量很庞大,甚至超过了解方程组的工作量．显然,如能找到适当求积公式离散（０．２）,则可节省大量计算．使用求积公式法解（０．２）的文献不多,［５］中提…     

DIRECTIONAL DO-NOTHING CONDITION FOR THE NAVIER-STOKES EQUATIONS

The numerical solution of flow problems usually requires bounded domains although the physical problem may take place in an unbounded or substantially larger domain. In this case, artificial boundaries are necessary. A well established artificial boundary condition for the Navier-Stokes equations diseretized by finite elements is the “do-nothing” condition. The reason for this is the fact that this condition appears automatically in the variational formulation after partial integration of the viscous term and the pressure gradient. This condition is one of the most established outflow conditions for Navier-Stokes but there are very few analytical insight into this boundary condition. We address the question of existence and stability of weak solutions for the Navier-Stokes equations with a “directional do-nothing” condition. In contrast to the usual “do-nothing” condition this boundary condition has enhanced stability properties. In the case of pure outflow, the condition is equivalent to the original one, whereas in the case of inflow a dissipative effect appears. We show existence of weak solutions and illustrate the effect of this boundary condition by computation of steady and non-steady flows.     

CONVERGENCE ANALYSIS OF THE JACOBI SPECTRAL-COLLOCATION METHOD FOR FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method,which shows that the errors of the approximate solution decay exponentially in L∞norm and weighted L2-norm. The numerical examples are given to illustrate the theoretical results.     

WEAK TRAVELLING WAVE FRONT SOLUTIONS OF GENERALIZED DIFFUSION EQUATIONS WITH REACTION

The author demonstrate that the two-point boundary value problem {p′(s)=f′(s)-λp^β(s)for s∈(0,1);β∈(0,1),p(0)=p(1)=0,p(s)&gt;0 if s∈(0,1),has a solution(λ^-，p^-（s）),where |λ^-| is the smallest parameter,under the minimal stringent restrictions on f(s), by applying the shooting and regularization methods. In a classic paper, Kohmogorov et.al.studied in 1937 a problem which can be converted into a special case of the above problem. The author also use the solution(λ^-,p^-(s)) to construct a weak travelling wave front solution u(x,t)=y(ξ),ξ=x-Ct,C=λ^-N/(N+1),of the generalized diffusion equation with reaction δ/δx（k（u）|δu/δx|^n-1 δu/δx)-δu/δt=g（u），where N&gt;0,k(s)&gt;0 a.e.on(0,1),and f(a):=n+1/N∫0ag(t)k^1/N(t)dt is absolutely continuous ou[0,1],while y(ξ) is increasing and absolutely continuous on (-∞,+∞) and (k(y(ξ))|y′(ξ)|^N)′=g(y(ξ))-Cy′(ξ)a.e.on(-∞，+∞）,y(-∞)=0,y(+∞)=1.     

关于非定常不可压Navier-Stokes方程的时间高精度隐式差分方法

The incompressible Navier-Stokes equations,upon spatial discretization,become a system of differential algebraic equations,formally of index2.But due to the special forms of the discrete gradient and disrete divergence,its index can be regarded as 1.Thus,in this paper,a systematic approach following the ODE theory and methods is presented for the construction of high-order time-accurate implicit schemes for the incompressible Navier-Stokes equations,with projection methods for efficiency of numerical solution.The 3rd order 3-step BDF with componentconsistent pressure-correction projection method is a first attempt in this direction;the related iterative solution of the auxiliary velocyty,the boundary conditions and the stability of the algorithm are discussed.Results of numerical tests on the incompressible Navier-Stokes equations with an exact solution are presented,confirming the accureacy,stability and component-consistency of the proposed method.     

ANALYSIS ON A NUMERICAL SCHEME WITH SECOND-ORDER TIME ACCURACY FOR NONLINEAR DIFFUSION EQUATIONS

A nonlinear fully implicit finite difference scheme with second-order time evolution for nonlinear diffusion problem is studied.The scheme is constructed with two-layer coupled discretization (TLCD) at each time step.It does not stir numerical oscillation,while per-mits large time step length,and produces more accurate numerical solutions than the other two well-known second-order time evolution nonlinear schemes,the Crank-Nicolson (CN)scheme and the backward difference formula second-order (BDF2) scheme.By developing a new reasoning technique,we overcome the difficulties caused by the coupled nonlinear discrete diffusion operators at different time layers,and prove rigorously the TLCD scheme is uniquely solvable,unconditionally stable,and has second-order convergence in both s-pace and time.Numerical tests verify the theoretical results,and illustrate its superiority over the CN and BDF2 schemes.