B-splines methods with redefined basis functions for solving fourth order parabolic partial differential equations

In this work, we discuss two methods for solving a fourth order parabolic partial differential equation. In Method-I, we decompose the given equation into a system of second order equations and solve them by using cubic B-spline method with redefined basis functions. In Method-II, the equation is solved directly by applying quintic B-spline method with redefined basis functions. Stability of these methods have been discussed. Both methods are unconditionally stable. These methods are tested on four examples. The computed results are compared wherever possible with those already available in literature. We have developed Method-I for fourth order non homogeneous parabolic partial differential equation from which we can obtain displacement and bending moment both simultaneously, while Method-II gives only displacement. The results show that the derived methods are easily implemented and approximate the exact solution very well.     

孤立波数值方法综述(Ⅰ)

This paper is devoted to introduce the history of numerical study of solitons and the main results for K. D. V. equation, R, L. W. equation, Klein-Gordon equation, Sine-Gordon equation and so on. The technique for constructing difference schemes based on conservation and characteristic is discussed. The two main methods for analysis of stability is given too. Then several finite element methods are discuaaed. The final part of this paper is for various. spectral methods.     

孤立波数值方法综述(Ⅰ)

This paper is devoted to introduce the history of numerical study of solitons and the main results for K. D. V. equation, R. L. W. equation, Klein-Gordon equation, Sine-Gordon equation and so on.The technique for constructing difference schemes based on conservation and characteristic is discussed. The two main methods for analysis of stability is given too. Then several finite element methods are discussed. The final part of this paper is for various-spectral methods.     

Preconditioned tensor format conjugate gradient squared and biconjugate gradient stabilized methods for solving stein tensor equations

This article is concerned with solving the high order Stein tensor equation arising in control theory. The conjugate gradient squared (CGS) method and the biconjugate gradient stabilized (BiCGSTAB) method are attractive methods for solving linear systems. Compared with the large-scale matrix equation, the equivalent tensor equation needs less storage space and computational costs. Therefore, we present the tensor formats of CGS and BiCGSTAB methods for solving high order Stein tensor equations. Moreover, a nearest Kronecker product preconditioner is given and the preconditioned tensor format methods are studied. Finally, the feasibility and effectiveness of the new methods are verified by some numerical examples.     

A diagonal splitting method for solving semidiscretized parabolic partial differential equations

In this work, a diagonal splitting idea is presented for solving linear systems of ordinary differential equations. The resulting methods are specially efficient for solving systems which have arisen from semidiscretization of parabolic partial differential equations (PDEs). Unconditional stability of methods for heat equation and advection–diffusion equation is shown in maximum norm. Generalization of the methods in higher dimensions is discussed. Some illustrative examples are presented to show efficiency of the new methods.     

改进的两步高阶Euler-Chebyshev方法

给出非线性方程求根的Euler-Chebyshev方法的改进方法,证明了方法的收敛性,它们七次和九次收敛到单根.给出数值试验,且与牛顿法及其它较高阶的方程求根方法做了比较.结果表明方法具有很好的优越性,它丰富了非线性方程求根的方法,在理论上和应用上都有一定的价值.     

Multi-symplectic methods for the Ito-type coupled KdV equation

In this paper, we find that the Ito-type coupled KdV equation can be written as a multi-symplectic Hamiltonian partial differential equation (PDE). Then, multi-symplectic Fourier pseudospectral method and multi-symlpectic wavelet collocation method are constructed for this equation. In the numerical experiments, we show the effectiveness of the proposed methods. Some comparisons between the proposed methods are also made with respect to global conservation properties.     

Exact solitary wave solutions for a generalized KdV–Burgers equation

Using the special truncated expansion method, the solitary wave solutions are constructed for the compound Korteweg–de Vries–Burgers (KdVB) equation. Exact and explicit solitary wave solutions for a generalized KdVB equation are obtained by introducing a suitable ansatz equation. The generalized two-dimensional KdVB equation is discussed. Some particular cases of the generalized KdVB equation are solved by using these methods.     

高阶广义KDV方程和KP方程的数值解法

采用一种线性隐格式来解3-阶和5-阶的广义Korteweg-De Vries(KDV)方程,对这种方法做一推广,就能应用到它的二维形式广义Kadomtsew-Petviashvili(KP)方程．这种方法是无条件稳定的,且是无损耗的．数值实验描述了一个线性孤波运动的情形以及两个孤波交互的情形,从结果来看,它们满足孤立子解的两个守恒-动量守恒和能量守恒．     

Incomplete Semiiterative Methods for Solving Operator Equations in Banach Space

There are several methods for solving operator equations in a Banach space. The successive approximation methods require the spectral radius of the iterative operator be less that 1 for convergence. In this paper, we try to use the incomplete semiiterative methods to solve a linear operator equation in Banach space. Usually the special semiiterative methods are convergent even when the spectral radius of the iterative operator of an operator of an operator equation is greater than 1.     

Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations

Perturbation methods depend on a small parameter which is difficult to be found for real-life nonlinear problems. To overcome this shortcoming, two new but powerful analytical methods are introduced to solve nonlinear heat transfer problems in this article; one is He's variational iteration method (VIM) and the other is the homotopy-perturbation method (HPM). The VIM is to construct correction functionals using general Lagrange multipliers identified optimally via the variational theory, and the initial approximations can be freely chosen with unknown constants. The HPM deforms a difficult problem into a simple problem which can be easily solved. Nonlinear convective–radiative cooling equation, nonlinear heat equation (porous media equation) and nonlinear heat equation with cubic nonlinearity are used as examples to illustrate the simple solution procedures. Comparison of the applied methods with exact solutions reveals that both methods are tremendously effective.     

一类四阶牛顿变形方法

给出非线性方程求根的一类四阶方法,也是牛顿法的变形方法.证明了方法收敛性,它们至少四次收敛到单根,线性收敛到重根.文末给出数值试验,且与牛顿法及其它牛顿变形法做了比较.结果表明方法具有很好的优越性,它丰富了非线性方程求根的方法,在理论上和应用上都有一定的价值.     

扩充偏微分方程(组)守恒律和对称的辅助方程方法及微分形式吴方法的应用

首先,我们给出了引入伴随方程(组)扩充原方程(组)的策略使给定偏微分方程(组)的扩充方程组具有对应泛瓯即,成为Lagrange系统的方法,以此为基础提出了作为偏微分方程(组)传统守恒律和对称概念的一种推广-偏微分方程(组)扩充守恒律和扩充对称的概念;其次,以得到的Lagrange系统为基础给定了确定原方程(组)扩充守恒律和扩充对称的方法,从而达到扩充给定偏微分方程(组)的首恒律和对称的目的;第三,提出了适用于一般形式微分方程(组)的计算固有守恒律的方法;第四,实现以上算法过程中,我们先把计算(扩充)守恒律和对称问题均归结为求解超定线性齐次偏微分方程组(确定方程组)的问题.然后,对此关键问题我们提出了用微分形式吴方法处理的有效算法;最后,作为方法的应用我们计算确定了非线性电报方程组在内的五个发展方程(组)的新守恒律和对称,同时也说明了方法的有效性.     

Higher order curvature information and its application in a modified diagonal Secant method

ABSTRACT

A secant equation (quasi-Newton) has one of the most important rule to find an optimal solution in nonlinear optimization. Curvature information must satisfy the usual secant equation to ensure positive definiteness of the Hessian approximation. In this work, we present a new diagonal updating to improve the Hessian approximation with a modifying weak secant equation for the diagonal quasi-Newton (DQN) method. The gradient and function evaluation are utilized to obtain a new weak secant equation and achieve a higher order accuracy in curvature information in the proposed method. Modified DQN methods based on the modified weak secant equation are globally convergent. Extended numerical results indicate the advantages of modified DQN methods over the usual ones and some classical conjugate gradient methods.     

互补问题算法的新进展

互补问题是一类重要的优化问题,在最近３０多年的时间里,人们为求解它而提出了许多算法,该文主要介绍１９９０－１９９７年之间出现的某些新算法,它们大致可归类为：（１）光滑方程法;（２）非光滑方程法;（３）可微无约束优化法;（４）ＧＬＰ投影法;（５）内点法;（６）磨光与非内点连续法,文中对每类算法及相应的收敛性结果做了描述与评论,并列出有关文献。     

极小解法求解最大-最小型关系方程的系统研究

最大-最小型关系方程已有许多学者进行过研究并给出一些解法,这些解法在应用中仍显得过于复杂.本文在文献[1]的基础上对最大-最小型关系方程的极小解的性质和特点进行了系统的研究,并对文献[1]提出的求解这类方程的图法和分枝法在理论上进行了补充和完善.研究表明,图法和分枝法是求解最大-最小型关系方程的有效方法.     

Finite volume and finite element methods for solving a one-dimensional space-fractional Boussinesq equation

In this paper, a new one-dimensional space-fractional Boussinesq equation is proposed. Two novel numerical methods with a nonlocal operator (using nodal basis functions) for the space-fractional Boussinesq equation are derived. These methods are based on the finite volume and finite element methods, respectively. Finally, some numerical results using fractional Boussinesq equation with the maximally positive skewness and the maximally negative skewness are given to demonstrate the strong potential of these approaches. The novel simulation techniques provide excellent tools for practical problems. These new numerical models can be extended to two- and three-dimensional fractional space-fractional Boussinesq equations in future research where we plan to apply these new numerical models for simulating the tidal water table fluctuations in a coastal aquifer.     

Sobolev-Volterra投影与积分微分方程有限元数值分析

本文提出一类称之为Sobolev-Volterra投影的有限元投影,研究了有关性质并将之应用于伪抛物型积分微分方程有限元方法、伪双曲型积分微分方程有限元方法以及三维伪双曲型积分微分方程交替方向有限元方法的数值分析.     

一类组合KdV方程的数值方法

1引言从Scott-Russell[1]提出在平静的水面上孤波运动的情形以后,大量的目光开始关注它的存在性、其中的性质和动态的交互情形[2],这是因为许多非线性动态的物理现象可以被描述成一个孤立子的模型[2,3].诸如,Korteweg-de Vries(KdV)方程[4,5,6],正弦     

Symplectic and multi-symplectic wavelet collocation methods for two-dimensional Schrödinger equations

In this paper, we develop symplectic and multi-symplectic wavelet collocation methods to solve the two-dimensional nonlinear Schrödinger equation in wave propagation problems and the two-dimensional time-dependent linear Schrödinger equation in quantum physics. The Hamiltonian and the multi-symplectic formulations of each equation are considered. For both formulations, wavelet collocation method based on the autocorrelation function of Daubechies scaling functions is applied for spatial discretization and symplectic method is used for time integration. The conservation of energy and total norm is investigated. Combined with splitting scheme, splitting symplectic and multi-symplectic wavelet collocation methods are also constructed. Numerical experiments show the effectiveness of the proposed methods.