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

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

Recovering the local volatility in Black–Scholes model by numerical differentiation

In this article, a numerical method for recovering the local volatility in Black–Scholes model is proposed based on the Dupire formula in which the numerical derivatives are used. By Tikhonov regularization, a new numerical differentiation method in two-dimensional (2-D) case is presented. The convergent analysis and numerical examples are also given. It shows that our method is efficient and stable.     

The numerical method of successive interpolations for Fredholm functional integral equations

A new numerical method for Fredholm functional integral equations is proposed. The method combines the fixed point technique with numerical integration and cubic spline interpolation. The convergence and the numerical stability of the method are proved and tested on some numerical examples.     

用线方法求解带混合边值条件的抛物方程

研制了分别用显式Euler法、隐式Euler法、Crank-Nicolson格式(梯形方法)求解带第一、第二及混合边值条件的抛物问题的应用软件,通过求解若干抛物问题对该软件作了测试,获得了预期的数值结果,讨论了时间和空间步长的变化对格式计算结果的影响,得到了三种方法的稳定性、收敛精度和计算量.     

Implicit numerical approximation scheme for the fractional Fokker-Planck equation

In this paper, we consider an anomalous subdiffusion process, governed by fractional Fokker-Planck equation. An effective numerical method for approximating Fokker-Planck equation in a bounded domain is presented. The stability and convergence of the numerical method are analyzed. Some numerical examples are presented to show the application of the present technique. The numerical results exhibit the good performance of our theoretical analysis.     

A method for verifying the accuracy of numerical solutions of symmetric saddle point linear systems

A fast numerical verification method is proposed for evaluating the accuracy of numerical solutions for symmetric saddle point linear systems whose diagonal blocks of the coefficient matrix are semidefinite matrices. The method is based on results of an algebraic analysis of a block diagonal preconditioning. Some numerical experiments are present to illustrate the usefulness of the method.     

Numerical solution of an inverse initial boundary-value problem for the full time-dependent Maxwell's equations in the presence of imperfections of small volume

We consider the numerical solution, in a three-dimensional bounded domain, of the inverse problem for identifying the location of small electromagnetic imperfections in a medium with homogeneous background. Our numerical algorithm is based on the coupling of a discontinuous Galerkin method for the time-dependent Maxwell's equations, on the exact controllability method and on a Fourier inversion. Several numerical results are given with one and two imperfections and the robustness and accuracy of the numerical method used for the dynamic detection problem are shown.     

A numerical approximation method for solving a three-dimensional space Galilei invariant fractional advection-diffusion equation

In this paper, a numerical approximation method for solving a three-dimensional space Galilei invariant fractional advection diffusion equation is presented. The convergence and stability of the numerical approximation method are discussed by a new technique of Fourier analysis. The solvability of the numerical approximation method also is analyzed. Finally, applying Richardson extrapolation technique, a high-accuracy algorithm is structured and the numerical example demonstrated the theoretical results.     

Numerical analysis for a variable-order nonlinear cable equation

In this paper, a variable-order nonlinear cable equation is considered. A numerical method with first-order temporal accuracy and fourth-order spatial accuracy is proposed. The convergence and stability of the numerical method are analyzed by Fourier analysis. We also propose an improved numerical method with second-order temporal accuracy and fourth-order spatial accuracy. Finally, the results of a numerical example support the theoretical analysis.     

Dimensional splitting with front tracking and adaptive grid refinement

Front tracking in combination with dimensional splitting is analyzed as a numerical method for scalar conservation laws in two space dimensions. An analytic error bound is derived, and convergence rates based on numerical experiments are presented. Numerical experiments indicate that large CFL numbers can be used without loss of accuracy for a wide range of problems. A new method for grid refinement is introduced. The method easily allows for dynamical changes in the grid, using, for instance, the total variation in each grid cell as a criterion for introducing new or removing existing refinements. Several numerical examples are included, highlighting the features of the numerical method. A comparison with a high-resolution method confirms that dimensional splitting with front tracking is a highly viable numerical method for practical computations. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 627–648, 1998     

Numerical analysis for a variable-order nonlinear cable equation

In this paper, a variable-order nonlinear cable equation is considered. A numerical method with first-order temporal accuracy and fourth-order spatial accuracy is proposed. The convergence and stability of the numerical method are analyzed by Fourier analysis. We also propose an improved numerical method with second-order temporal accuracy and fourth-order spatial accuracy. Finally, the results of a numerical example support the theoretical analysis.     

Numerical method with high order accuracy for solving a anomalous subdiffusion equation

In this paper, a numerical method with second order temporal accuracy and fourth order spatial accuracy is developed to solve a anomalous subdiffusion equation; by Fourier analysis, the convergence, stability and solvability of the numerical method are analyzed; the theoretical results are strongly supported by the numerical experiment.     

Numerical model for the shallow water equations on a curvilinear grid with the preservation of the Bernoulli integral

Methodological aspects concerning the construction of a two-dimensional numerical model for reservoir flows based on the shallow water equations are considered. A numerical scheme is constructed by applying the control volume method on staggered grids in combination with the Bernoulli integral, which is used to interpolate the desired fields inside a grid cell. The implementation of the method yields a monotone numerical scheme. The results of numerical integration are compared with the exact solution.     

A semi‐discrete defect correction finite element method for unsteady incompressible magnetohydrodynamics equations

In this report, we give a semi‐discrete defect correction finite element method for the unsteady incompressible magnetohydrodynamics equations. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear magnetohydrodynamics equations is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect‐correction technique. Then, we give the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. In order to show the effect of our method, some numerical results are shown. Copyright © 2017 John Wiley & Sons, Ltd.     

一类强色散DGH方程的多辛普雷斯曼格式

DGH方程作为一类重要的非线性水波方程有着许多广泛的应用前景.基于Hamilton系统的多辛理论研究了一类强色散DGH方程的数值解法,利用多辛普雷斯曼方法构造了一种典型的半隐式的多辛格式.分析了该格式的局部能量和动量守恒律误差,并给出了数值算例.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

An efficient algorithm for solving generalized pantograph equations with linear functional argument

A numerical method for solving the generalized (retarded or advanced) pantograph equation under initial value conditions is presented. To display the validity and applicability of the numerical method four illustrative examples are presented. The results reveal that this method is very effective and highly promising when compared with other numerical methods, such as Adomian decomposition method, spline methods and Taylor method.     

A finite volume element formulation and error analysis for the non-stationary conduction–convection problem

The non-stationary conduction–convection problem including the velocity vector field and the pressure field as well as the temperature field is studied with a finite volume element (FVE) method. A fully discrete FVE formulation and the error estimates between the fully discrete FVE solutions and the accuracy solution are provided. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the FVE method is feasible and efficient for finding the numerical solutions of the non-stationary conduction–convection problem and is one of the most effective numerical methods by comparing the results of the numerical simulations of the FVE formulation with those of the numerical simulations of the finite element method and the finite difference scheme for the non-stationary conduction–convection problem.     

单物种人口模型指数隐式Euler方法的振动性

本文研究了用以描述单物种人口模型的延迟Logistic方程的数值振动性.对方程应用隐式Euler方法进行求解,针对离散格式定义了指数隐式Euler方法,证明了该方法的收敛阶为1.根据线性振动性理论获得了数值解振动的充分条件.进而还对非振动数值解的性质作了讨论.最后用数值算例对理论结果进行了验证.     

Convergence and Mean-Square Stability of Exponential Euler Method for Semi-Linear Stochastic Delay Integro-Differential Equations

In this paper, the numerical methods for semi-linear stochastic delay integro-differential equations are studied. The uniqueness, existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained. Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied. It is proved that the exponential Euler method is convergent with strong order $\frac{1}{2}$ and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size. In addition, numerical experiments are presented to confirm the theoretical results.     

Numerical differentiation with noisy signal

The aim of this paper is to present a new numerical method, which ables one to filter and compute numerical derivatives of a function whose values are known in some points from experimental measurements, inducing noisy data. We use a piecewise cubic spline interpolation to generate a function whose Fourier coefficients give an approximation of the numerical derivatives we are looking for. Error and stability analysis of this numerical algorithm are provided. Numerical results are presented for data smoothing and for the first and second derivatives computed from noisy data. They show that this method gives good numerical results. Comparison with other methods is done.