Numerical integration of partial differential equations using cubic splines

This paper presents techniques for the numerical solution of partial differential equations using cubic spline collocation.

The main spline relations are presented and incorporated into solution procedures for partial differential equations. The computational algorithm in every case is a tridiagonal matrix system amenable to efficient inversion methods. Truncation errors and stability are briefly discussed. Finally, some examples of their application to parabolic and hyperbolic systems with mixed boundary conditions are presented.

The results obtained are encouraging and justify further research in this field.     

Polynomial solutions of certain differential equations

In this paper, the Chebyshev matrix method is applied generalisations of the Hermite, Laguerre, Legendre and Chebyshev differential equations which have polynomial solution. The method is based on taking the truncated Chebyshev series expansions of the functions in equation, and then substituting their matrix forms into the result equation. Thereby the given equation reduces to a matrix equation, which corresponds to a system of linear algebraic equations with unknown Chebyshev coefficients.     

Parallelizable restarted iterative methods for nonsymmetric linear systems. part I: Theory

Large sparse nonsymmetric problems of the form A u = b are frequently solved using restarted conjugate gradient-type algorithms such as the popular GCR and GMRES algorithms. In this study we define a new class of algorithms which generate the same iterates as the standard GMRES algorithm but require as little as half of the computational expense. This performance improvement is obtained by using short economical three-term recurrences to replace the long recurrence used by GMRES. The new algorithms are shown to have good numerical properties in typical cases, and the new algorithms may be easily modified to be as numerically safe as standard GMRES. Numerical experiments with these algorithms are given in Part II, in which we demonstrate the improved performance of the new schemes on different computer architectures.     

Third-order temporal discrete scheme for the non-stationary Navier–Stokes equations

In this article, we suggest a new third-order time discrete scheme for the two-dimensional non-stationary Navier–Stokes equations. After presenting the Galerkin finite element approximation for the spatial discretization, we consider an implicit/explicit time discrete scheme for the problem, which is based on the two-step Adams–Moulton scheme (implicit scheme) for the linear term and the three-step Adams–Bashforth scheme (explicit scheme) for the nonlinear term. In this method, we only need to solve a linearized discrete system at each time step, so the scheme can converge fast and the computational cost can be reduced. Moreover, under some assumptions, we deduce the stability and optimal error estimate for the velocity in L ²-norm.     

分数阶微积分在图像处理中的研究综述*

综述了关于分数阶微积分理论在数字图像底层处理中的应用研究,具体包括:分数阶微积分和分数阶偏微分方程的基本理论及分数阶傅里叶变换的基本性质;分数阶微分滤波器的构造及在图像增强中的应用研究;分数阶积分滤波器的构造及在图像去噪中的应用研究;分数阶偏微分方程在图像去噪中的应用研究。最后,总结了分数阶微积分理论在图像处理中已取得的研究成果,并结合已有的基于分数阶微积分理论的图像底层处理模型,展望了分数阶微积分理论在图像处理中的应用前景。     

Nodal DG-FEM solution of high-order Boussinesq-type equations

A discontinuous Galerkin finite-element method (DG-FEM) solution to a set of high-order Boussinesq-type equations for modelling highly nonlinear and dispersive water waves in one horizontal dimension is presented. The continuous equations are discretized using nodal polynomial basis functions of arbitrary order in space on each element of an unstructured computational domain. A fourth-order explicit Runge-Kutta scheme is used to advance the solution in time. Methods for introducing artificial damping to control mild nonlinear instabilities are also discussed. The accuracy and convergence of the model with both h (grid size) and p (order) refinement are confirmed for the linearized equations, and calculations are provided for two nonlinear test cases in one horizontal dimension: harmonic generation over a submerged bar, and reflection of a steep solitary wave from a vertical wall. Test cases for two horizontal dimensions will be considered in future work.     

一次同余方程组的一种矩阵解法

同余方程组是数论中一个极为有用的课题。文中考虑的是具有一个未知数而具有不同的模的一次同余方程组。一般可以用中国剩余定理或者递推算法等方法给出一次同余方程组的解法。利用矩阵的初等变换和矩阵的Sm ith标准形给出了求解一次同余方程组的一种矩阵解法。     

The generalized Riemann problem method for the shallow water equations with bottom topography

This paper extends the generalized Riemann problem method (GRP) to the system of shallow water equations with bottom topography. The main contribution is that the generalized Riemann problem method (J. Comput. Phys. 1984; 55 (1):1–32) is used to evaluate the midpoint values of solutions at each cell interface so that the bottom topography effect is included in numerical fluxes, and at the same step the source term is discretized with an interface method in which only mid‐point values are plugged in. This scheme is well balanced between the flux gradient and bottom topography when incorporating the surface gradient method (SGM) (J. Comput. Phys. 2001; 168 (1):1–25) into data reconstruction step, and it is also suitable for both steady and unsteady flow simulations. We illustrate the accuracy of this scheme by several 1‐D and 2‐D numerical experiments. Copyright © 2005 John Wiley & Sons, Ltd.     

Influence of Implicit Integration Scheme on Prediction of Shear Band Formation

Implicit integration schemes for elastoplastic constitutive equations have been developed in recent years as an alternative to explicit schemes. The consistent tangent constitutive matrix Dcon that results from implicit schemes makes the global stiffness matrix consistent with the implicit integration procedure and differs from the traditional continuum tangent constitutive matrix Dep that results from explicit schemes. Onset of strain localization and shear banding has been traditionally predicted using the continuum tangent constitutive matrix. It is shown that different criteria for onset of shear-band formation are obtained depending on whether Dcon or Dep is used. It is shown that shear band prediction using Dcon is step-size dependent, and that the use of Dcon influences the predicted onset of strain localization in frictional materials. An analytical equation for prediction of the onset of shear-band formation using Dcon for the Mohr–Coulomb model is developed, and a numerical example is presented.     

微振动对空间实验室微重力环境的影响研究

针对大型刚柔组合体空间结构(空间实验室)的特点,建立系统的动力学方程,并对机柜内设备、太阳阵的转动、宇航员活动产生的激扰力进行了简化。利用计算机建立组合体的简化模型,对空间组合体在各种微振动激扰下的加速度响应进行了仿真研究,并将这些仿真结果用于空间实验室微重力环境的预测。