三维数值流形方法的理论研究

在二维数值流形方法的基础上,对三维数值流形进行了理论研究．研究了三维覆盖位移函数,进行了三维数值流形的力学分析,给出了三维流形单元的刚度矩阵,详细推导了三维数值流形的Hammer积分及剖分规则,系统地研究三维数值流形的理论体系与数值实现方法．作为数值算例,给出了相应的悬臂梁的计算结果,计算结果表明算法的精度和计算效益较高．     

渗流问题灰色数值模型的解法研究

灰色数值模型的求解是研究灰色数值模型的一个重要问题 .本文根据灰集合、灰数及其灰色运算规则 ,在渗流系统的基本灰色数值模型的基础上 ,分析了求解这类模型的一整套灰色数值算法 ,并对灰色数值算法、普通算法和经典数值方法的计算结果进行了全面比较 ,论证了灰色数值算法对灰信息传递的正确性和对渗流系统描述的合理性 .     

时滞积分微分方程的Rosenbrock方法

本文研究时滞积分微分方程的数值方法.通过改造现有常及离散型延迟微分方程的数值方法,并匹配以适当数值求积公式,构造了求解时滞积分微分方程的Rosenbrock方法,导出了其稳定性准则.数值例子阐明了所获方法的计算有效性.     

数值流形方法的变分原理与应用

针对线弹性体静力问题,根据数值流形方法的特点及相应的位移模式,得到了面向物理覆盖的数值流形方法的变分原理,详细推导了基于变分原理的数值流形方法的理论计算公式,建立了数值流形方法的控制方程。作为实际应用,给出了相应的数值算例,结果表明,求解精度和效益令人满意。     

Abel变换数值反演的积分算子方法

本文研究了Abel变换的数值反演问题.利用Abel变换的理论反演公式与数值求导的积分算子法相结合的方法,对反演公式中奇异积分合理处理,获得Abel变换数值反演的一种算法,并进行了理论分析与数值实验. 结果表明该算法具有计算简单、数值稳定等优点.     

具有马尔可夫调制的随机微分方程数值解的收敛性和稳定性

研究了一类具有马尔可夫调制的线性随机微分方程Euler数值解的收敛性和稳定性,建立了Euler数值解MS稳定性的定义,确定了Euler数值解MS稳定的条件.     

《高等数值分析》教学案例的建设与思考

《高等数值分析》是一门与实际联系紧密的数学公共课程,它理论深刻,应用广泛.本文结合实际应用,为《高等数值分析》中常微分方程数值解部分设计了一个教学案例,通过理论分析和数值实验向学生展示了刚性问题的概念和相关数值方法,并对《高等数值分析》课程教学案例的设计进行了思考.     

复频域卷积的离散算法

导出了适用于计算机进行计算的复频域卷积的离散算法,应用Durbin拉氏变换数值反演法对复频域卷积结果进行数值反演,可获得时域数值解.将该数值解与解析解进行比较表明,数值解具有较高的精度.     

数值分析双语教学的探讨

首先指出了数值分析双语教学的目的和意义.其次,分析了数值分析双语教学中存在的问题.最后,给出了数值分析双语教学的一些措施.     

线性传输方程的几种数值格式的比较

本文研究了线性传输方程的数值计算问题.利用Godunov格式、Entropy格式、Ultra-bee格式和Entropy-Ultra-bee格式对线性传输方程进行了数值计算,获得了相应的数值结果.数值实验结果表明Entropy-Ultra-bee格式结合了Entropy格式和Ultra-bee格式的优点,在整个计算区域都有比较高的分辨率,而且没有出现非物理振荡.     

LONG TIME ASYMPTOTIC BEHAVIOR OF SOLUTION OF IMPLICIT DIFFERENCE SCHEME FOR A SEMI-LINEAR PARABOLIC EQUATION

In this paper, the solution of back-Euler implicit difference scheme for a semi-linea rparabolic equation is proved to converge to the solution of difference scheme for the corresponding semi-linear elliptic equation as t tends to infinity. The long asymptotic behavior of its discrete solution is obtained which is analogous to that of its continuous solution. At last, a few results are also presented for Crank-Nicolson scheme.     

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

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

de Rham transforms for subdivision schemes

For any subdivision scheme, we define its de Rham transform, which generalizes the de Rham and Chaikin corner cutting. The main property of the de Rham transform is that it preserves a sum rule. This allows comparison of the Hölder regularity of a given subdivision scheme with that of its de Rham transform. A graphical comparison is made for three different families of subdivision schemes, the last one being the generalized four-point scheme.     

非线性Schrdinger方程初边值问题的守恒数值格式

该文对非线性 Schrodinger方程提出了一种新的守恒差分格式 ,并证明了该格式的收敛性与稳定性 ,通过数值计算获得如下结论 ,提出的差分格式在取适当的参数后 ,精度上好于文 [7]中的格式     

A high-order compact difference scheme for 2D Laplace and Poisson equations in non-uniform grid systems

In this study, a high-order compact scheme for 2D Laplace and Poisson equations under a non-uniform grid setting is developed. Based on the optimal difference method, a nine-point compact difference scheme is generated. Difference coefficients at each grid point and source term are derived. This is accomplished through the consideration of compatibility between the partial differential equation and its difference discretization. Theoretically, the proposed scheme has third- to fourth-order accuracy; its fourth-order accuracy is achieved under uniform grid settings. Two examples are provided to examine performance of the proposed scheme. Compared with the traditional five-point difference scheme, the proposed scheme can produce more accurate results with faster convergence. Another reference scheme with the same nine-point grid stencil is derived based on the five-point scheme. The two nine-point schemes have the same coefficients for each grid points; however, their coefficients for the source term are different. The overall accuracy level of the solution resulting from the proposed scheme is higher than that of the nine-point reference scheme. It is also indicated that the smoothness of grids has significant effects on accuracy and convergence of the solutions; efforts in optimizing the grid configuration and allocation can improve solution accuracy and efficiency. Consequently, with the proposed method, solution under the non-uniform grid setting with appropriate grid allocation would be more accurate than that under the uniform-grid manipulation, with the same number of grid points.     

二维热传导方程的三层显式差分格式

对二维热传导方程构造了一个稳定的三层显式差分格式求其数值解,其背景源于高维热力学反问题迭代算法中对正问题小计算量算法的需求。首先建立一个含参数的一般差分格式去逼近微分方程,并得到了最优截断误差。然后导出了参数应满足的条件以保证差分格式的稳定性。最后给出了数值的例子并和其它算法进行比较,说明了格式在精度上的有效性和计算量上的优越性。     

Locally one-dimensional difference scheme for a fractional tracer transport equation

A locally one-dimensional scheme for a fractional tracer transport equation of order is considered. An a priori estimate is obtained for the solution of the scheme, and its convergence is proved in the uniform metric.     

A new line search inexact restoration approach for nonlinear programming

A new general scheme for Inexact Restoration methods for Nonlinear Programming is introduced. After computing an inexactly restored point, the new iterate is determined in an approximate tangent affine subspace by means of a simple line search on a penalty function. This differs from previous methods, in which the tangent phase needs both a line search based on the objective function (or its Lagrangian) and a confirmation based on a penalty function or a filter decision scheme. Besides its simplicity the new scheme enjoys some nice theoretical properties. In particular, a key condition for the inexact restoration step could be weakened. To some extent this also enables the application of the new scheme to mathematical programs with complementarity constraints.     

CONSERVATION OF THREE—POINT COMPACT SCHEMES ON SINGLE AND MULTIBLOCK PATCHED GRIDS FOR HYPERBOLIC PROBLEMS

For nonlinear hyperbloic problems,Conservation of the numerical scheme is important for convergence to the correct weak solutions.In this paper the the conservation of the well-known compact scheme up to fourth order of accuracy on a single and uniform grid is studied,and a conservative interface treatment is derived for compact schemes on patched grids .For a pure initial value problem,the compact scheme is shown to be equivalent to a scheme in the usual conservative form .For the case of a mixed initial boundary value problem,the compact scheme is conservative only if the rounding errors are small enough.For a pactched grid interface,a conservative interface condition useful for mesh fefiement and for parallel computation is derived and its order of local accuracy is analyzed.     

The Milstein Scheme for Stochastic Delay Differential Equations Without Using Anticipative Calculus

The Milstein scheme is the simplest nontrivial numerical scheme for stochastic differential equations with a strong order of convergence one. The scheme has been extended to the stochastic delay differential equations but the analysis of the convergence is technically complicated due to anticipative integrals in the remainder terms. This article employs an elementary method to derive the Milstein scheme and its first order strong rate of convergence for stochastic delay differential equations.