一类偏积分微分方程二阶差分全离散格式

本文给出了数值求解一类偏积分微分方程的二阶差分全离散格式．时间方向采用了二阶向后差分格式,积分项的离散利用了Lubich的二阶卷积求积公式,给出了稳定性的证明、误差估计及收敛性的结果,并给出了数值例子．     

一类非线性Scllr(o)dillger方程的差分/Legendre谱元法

考察一类带幂次非线性项的Schrodinger方程的Dirichlet初边值问题,提出了一个有效的计算格式,其中时间方向上应用了一种守恒的二阶差分隐格式,空间方向上采用Legendre谱元法.对于时间半离散格式,证明了该格式具有能量守恒性质,并给出了L2误差估计,对于全离散格式,应用不动点原理证明了数值解的存在唯一性,并给出了L2误差估计.最后,通过数值试验验证了结果的可信性.     

一类非线性Schrodinger方程的差分／Legendre谱元法

考察一类带幂次非线性项的Schrodinger方程的Dirichlet初边值问题，提出了一个有效的计算格式，其中时间方向上应用了一种守恒的二阶差分隐格式，空间方向上采用Legendre谱元法．对于时间半离散格式，证职了该格式具有能量守恒性质，并给出了L^2误差估计，对于全离散格式，应用不动点原理证明了数值解的存在唯一性，并给出了L^2误差估计．最后，通过数值试验验证了结果的可信性．     

带初始奇异性的多项时间分数阶扩散方程一类全离散数值方法的误差分析

本文研究了带有初始奇异性的多项时间分数阶扩散方程的一种全离散数值方法.首先,基于L1公式在渐变网格下离散多项Caputo时间分数阶导数,构造了多项时间分数阶扩散方程的时间半离散格式,证明了时间格式通过选取合适的网格参数r,时间方向的误差可以达到最优的收敛阶2-α_1,其中α_1(0 α_11)为多项时间分数阶导数阶数的最大值.然后,空间采用谱方法进行离散,得到了全离散格式,证明了全离散格式的无条件稳定性和收敛性.为了降低计算量和储存量,对多项时间分数阶扩散方程又构造了时间方向的快速算法,同时证明了该格式的收敛性.数值算例验证了算法的有效性,显示了快速算法的高效性.     

非饱和土壤水流问题的CN有限元格式

首先给出二维非饱和土壤水流问题基于Crank-Nicolson(CN)方法的具有时间二阶精度的半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN有限元格式,并给出误差估计,最后用数值例子说明全离散化CN有限元格式的优越性.这种方法可以绕开关于空间变量的半离散化格式的讨论,提高时间离散的精度,极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.     

无界域上半线性强阻尼波动方程的全离散有理谱逼近

本文运用Chebyshev有理谱方法来讨论半线性强阻尼波动方程.通过建立时间、空间方向全离散的Chebyshev有理谱格式,证明了由此格式所确定的离散算子半群存在整体吸引子,并从理论上建立了在有限时间上近似解的误差估计.     

一类非线性偏积分微分方程二阶差分全离散格式

给出了数值求解一类非线性偏积分微分方程的二阶全离散差分格式.采用了二阶向后差分格式,积分项的离散利用了Lubich的二阶卷积求积公式,给出了稳定性的证明、误差估计及收敛性的结果.     

关于九参数拟协调板元

1980年以来,唐立民等提出一种拟协调元法,用来构造椭圆型方程的离散格式.粗略地讲,该法将每个单元上的能量表达式所含导数项的面积分(假设问题二维的),用格林公式转化为单元边界上的线积分,然后采用某种数值积分,将线积分进行离散.对只含函数项的面积分,也用相应的数值积分进行离散.用此法计算单元刚度阵,比较简单、灵活.     

一类偏积分微分方程二阶差分全离散格式

本给出了数值求解一类偏积分微分方程的二阶全离散差分格式．采用了Crank-Nicolson格式；积分项的离散利用了Lubieh的二阶卷积积分公式；给出了稳定性的证明，误差估计及收敛性的结果．     

二维半线性反应扩散方程的交替方向隐格式

本文研究一类二维半线性反应扩散方程的差分方法.构造了一个二层线性化交替方向隐格式.利用离散能量估计方法证明了差分格式解的存在唯一性、差分格式在离散H~1模下的二阶收敛性和稳定性.最后给出两个数值例子验证了理论分析结果.     

The Asymptotic Behavior for Numerical Solution of a Volterra Equation

Long-time asymptotic stability and convergence properties for the numerical solution of a Volterra equation of parabolic type are studied.The methods are based on the first-second order backward difference methods.The memory term is approximated by the comvolution quadrature and the interpolant quadrature.Discretization of the spatial partial differential operators by the finite element method is also considered.     

A GENERALIZED GAUSS-TYPE QUADRATURE FORMULA AND ITS APPLICATIONS TO PSEUDOSPECTRAL METHOD

1　ＩｎｔｒｏｄｕｃｔｉｏｎＳｐｅｃｔｒａｌｍｅｔｈｏｄｅｍｐｌｏｙｓｇｌｏｂａｌｐｏｌｙｎｏｍｉａｌｓａｓｔｒｉａｌｆｕｎｃｔｉｏｎｓｉｎｓｐａｔｉａｌｄｉｓｃｒｅｔｉｚａｔｉｏｎｏｆｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓ.Ｉｔｏｆｔｅｎｐｒｏｖｉｄｅｓａｃｃｕｒａｔｅｎｕｍｅｒｉｃａｌｒｅｓｕｌｔｓｗｉｔｈｒｅｌａｔｉｖｅｌｙｌｅｓｓｄｅｇｒｅｅｏｆｆｒｅｅｄｏｍｓ,ａｎｄｔｈｕｓｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｕｓｅｄｆｏｒｃｏｍｐｕｔａｔｉｏｎｓｉｎｓｃｉｅｎｃｅａｎｄｅｎｇｉｎ…     

The Uniform <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> Behavior for Time Discretization of an Evolution Equation

The uniform L^2 stability and convergence properties for the time discretization of an evolution equation with a memory term are studied.The methods are based on the second-order backward difference methods.The memory term is approximated by the second-order convolution quadrature and interpolant quadrature.     

非线性积分方程迭代配置法的渐近展开及其外推

非线性积分方程迭代配置法的渐近展开及其外推韩国强（华南理工大学计算机工程与科学系）ＡＳＹＭＰＴＯＴＩＣＥＲＲＯＲＥＸＭＮＳＩＯＮＳＡＮＤＥＸＴＲＡＰＯＬＡＴＩＯＮＦＯＲＴＨＥＩＴＥＲＡＴＥＤＣＯＬＬＯＣＡＴＩＯＮＭＥＴＨＯＤＳＯＦＮＯＮＬＩＮＥＡＲＩ...     

The numerical evaluation of the error term in Gaussian quadrature rules

A method for the numerical evaluation of the error term in Gaussian quadrature rules is derived by means of Chebyshev polynomials of the first kind.     

LOCAL GAUSSIAN-COLLOCATION SCHEME TO APPROXIMATE THE SOLUTION OF NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS USING VOLTERRA INTEGRAL EQUATIONS

This work describes an accurate and effective method for numerically solving a class of nonlinear fractional differential equations.To start the method,we equivalently convert these types of differential equations to nonlinear fractional Volterra integral equations of the second kind by integrating from both sides of them.Afterward,the solution of the mentioned Volterra integral equations can be estimated using the collocation method based on locally supported Gaussian functions.The local Gaussian-collocation scheme estimates the unknown function utilizing a small set of data instead of all points in the solution domain,so the proposed method uses much less computer memory and volume computing in comparison with global cases.We apply the composite non-uniform Gauss-Legendre quadrature formula to estimate singular-fractional integrals in the method.Because of the fact that the proposed scheme requires no cell structures on the domain,it is a meshless method.Furthermore,we obtain the error analysis of the proposed method and demon-strate that the convergence rate of the approach is arbitrarily high.Illustrative examples clearly show the reliability and efficiency of the new technique and confirm the theoretical error estimates.     

High-order ADI orthogonal spline collocation method for a new 2D fractional integro-differential problem

We use the generalized L1 approximation for the Caputo fractional derivative, the second-order fractional quadrature rule approximation for the integral term, and a classical Crank-Nicolson alternating direction implicit (ADI) scheme for the time discretization of a new two-dimensional (2D) fractional integro-differential equation, in combination with a space discretization by an arbitrary-order orthogonal spline collocation (OSC) method. The stability of a Crank-Nicolson ADI OSC scheme is rigourously established, and error estimate is also derived. Finally, some numerical tests are given.     

The backward euler anisotropic a posteriori error analysis for parabolic integro‐differential equations

We derive two optimal a posteriori error estimators for an implicit fully discrete approximation to the solutions of linear integro‐differential equations of the parabolic type. A continuous, piecewise linear finite element space is used for the space discretization and the time discretization is based on an implicit backward Euler method. The a posteriori error indicator corresponding to space discretization is derived using the anisotropic interpolation estimates in conjunction with a Zienkiewicz‐Zhu error estimator to approach the error gradient. The error due to time discretization is derived using continuous, piecewise linear polynomial in time. We use the linear approximation of the Volterra integral term to estimate the quadrature error in the second estimator. Numerical experiments are performed on the isotropic mesh to validate the derived results.© 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1309–1330, 2016     

Error estimates for Filon quadrature formulae

In this paper, upper bounds for the error of (generalized) Filon quadrature formulae are stated. Furthermore, the main term of this error is derived, yielding simple modified quadrature rules of higher asymptotical precision.     

About some quadrature formulas

In this paper, we studied a class of quadrature formulas obtained by using the connection between the monospline functions and the quadrature formulas. For this class we obtain the optimal quadrature formula with regard to the error and we give some inequalities for the remainder term of this optimal quadrature formula.