两点混合边值问题的紧有限体积格式

本文针对常系数和变系数两点混合边值问题提出一种紧有限体积格式,该格式形成的线性代数方程组具有三对角性质,可以使用追赶法求解.证明格式按照H1半范数具有四阶收敛精度.利用节点计算值,给出单元中点值和一阶导数值的高精度后处理计算公式,这两个公式同样具有四阶精度.数值算例验证了理论分析的正确性,并说明了格式的有效性.     

A noniterative algebraic solution for Riccati equations satisfying two-point boundary-value problems

A noniterative algebraic method is presented for solving differential Riccati equations which satisfy two-point boundary-value problems. This class of numerical problems arises in quadratic optimization problems where the cost functionals are composed of both continuous and discrete state penalties, leading to piecewise periodic feedback gains. The necessary condition defining the solution for the two-point boundary value problem is cast in the form of a discrete-time algebraic Riccati equation, by using a formal representation for the solution of the differential Riccati equation. A numerical example is presented which demonstrates the validity of the approach.The authors would like to thank Dr. Fernando Incertis, IBM Madrid Scientific Center, who reviewed this paper and pointed out that the two-point boundary-value necessary condition could be manipulated into the form of a discrete-time Riccati equation. His novel approach proved to be superior to the authors' previously proposed iterative continuation method.     

An Explicit Pseudo-Spectral Scheme with Almost Unconditional Stability for the Cahn-Hilliard Equation

In this paper, an explicit fully discrete three-level pseudo-spectral scheme with almost unconditional stability for the Cahn-Hilliard equation is proposed. Stability and convergence of the scheme are proved by Sobolev's inequalities and the bounded extensive method of the nonlinear function (B.N. Lu (1995)). The scheme possesses the almost same stable condition and convergent accuracy as the Creak-Nicloson scheme but it is an explicit scheme. Thus the iterative method to solve nonlinear algebraic system is avoided. Moreover, the linear stability of the critical point $u_0$ is investigated and the linear dispersive relation is obtained. Finally, the numerical results are supplied, which check the theoretical results.     

Analysis of new conservative difference scheme for two-dimensional Rosenau-RLW equation

In this article, numerical solution for the Rosenau-RLW equation in 2D is considered and a conservative Crank–Nicolson finite difference scheme is proposed. Existence of the numerical solutions for the difference scheme has been shown by Browder fixed point theorem. A priori bound and uniqueness as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability and a second-order accuracy on both space and time of the difference scheme are proved. Numerical experiments are given to support our theoretical results.     

A Generalized Lyapunov-Sylvester Computational Method for Numerical Solutions of NLS Equation with Singular Potential

In the present paper a numerical method is developed to approximate the solution of two-dimensional Nonlinear Schrdinger equation in the presence of a singular potential. The method leads to generalized Lyapunov-Sylvester algebraic operators that are shown to be invertible using original topological and differential calculus issued methods. The numerical scheme is proved to be consistent, convergent and stable using the Lyapunov criterion, lax equivalence theorem and the properties of the generalized Lyapunov-Sylvester operators.     

Fourier spectral approximation to long-time behavior of three dimensional Ginzburg-Landau type equation*

In this paper, we consider a three dimensional Ginzburg–Landau type equation with a periodic initial value condition. A fully discrete Galerkin–Fourier spectral approximation scheme is constructed, and then the dynamical properties of the discrete system are analyzed. First, the existence and convergence of global attractors of the discrete system are proved by a priori estimates and error estimates of the discrete solution, and the numerical stability and convergence of the discrete scheme are proved. Furthermore, the long-time convergence and stability of the discrete scheme are proved. *This work was supported by the National Natural Science Foundation of China (No.: 10432010 and 10571010)     

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation (PME)is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13–32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on $f$log $f$ as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete $W$^1,∞norm and a refined estimate is applied to derive the optimal error order.     

Klein-Gordon方程初达值问题的一个新的守恒差分格式

本文对非线性Klein－Gordon（NKG）方程的初边值问题提出了一种新的差分格式,它保持了NKG方程初边值问题的能量守恒．证明了该格式的收敛性和稳定性．特别地,由于该格式是完全隐式的,故对求长时解有着重要的作用．数值计算结果表明该方法计算速度快,精度好．     

Combined higher order finite volume and finite element scheme for double porosity and non-linear adsorption of transport problem in porous media

In this work, a dual porosity model of reactive solute transport in porous media is presented. This model consists of a nonlinear-degenerate advection-diffusion equation including equilibrium adsorption to the reaction combined with a first-order equation for the non-equilibrium adsorption interaction processes. The numerical scheme for solving this model involves a combined high order finite volume and finite element scheme for approximation of the advection-diffusion part and relaxation-regularized algorithm for nonlinearity-degeneracy. The combined finite volume-finite element scheme is based on a new formulation developed by Eymard et al. (2010) [10]. This formulation treats the advection and diffusion separately. The advection is approximated by a second-order local maximum principle preserving cell-vertex finite volume scheme that has been recently proposed whereas the diffusion is approximated by a finite element method. The result is a conservative, accurate and very flexible algorithm which allows the use of different mesh types such as unstructured meshes and is able to solve difficult problems. Robustness and accuracy of the method have been evaluated, particularly error analysis and the rate of convergence, by comparing the analytical and numerical solutions for first and second order upwind approaches. We also illustrate the performance of the discretization scheme through a variety of practical numerical examples. The discrete maximum principle has been proved.     

Finite difference discretization of the Rosenau-RLW equation

In this paper, numerical solutions of the Rosenau-RLW equation are considered using Crank–Nicolson type finite difference method. Existence of the numerical solutions is derived by Brouwer fixed point theorem. A priori bound and the error estimates as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability, second-order convergence and uniqueness of the scheme are proved using discrete energy method. Some numerical experiments have been conducted in order to verify the theoretical results.     

Analysis of finite element method for one-dimensional time-dependent Schrödinger equation on unbounded domain

This paper addresses the theoretical analysis of a fully discrete scheme for the one-dimensional time-dependent Schrödinger equation on unbounded domain. We first reduce the original problem into an initial-boundary value problem in a bounded domain by introducing a transparent boundary condition, then fully discretize this reduced problem by applying Crank–Nicolson scheme in time and linear or quadratic finite element approximation in space. By a rigorous analysis, this scheme has been proved to be unconditionally stable and convergent, its convergence order has also be obtained. Finally, two numerical examples are performed to show the accuracy of the scheme.     

Rosenbrock schemes with complex coefficients for stiff and differential algebraic systems

Many applied problems are described by differential algebraic systems. Complex Rosenbrock schemes are proposed for the numerical integration of differential algebraic systems by the ?-embedding method. The method is proved to converge quadratically. The scheme is shown to be applicable even to superstiff systems. The method makes it possible to perform computations with a guaranteed accuracy. An equation is derived that describes the leading term of the error in the method as a function of time. An algorithm extending the method to systems of differential equations for complex-valued functions is proposed. Examples of numerical computations are given.     

广义混合非线性Schr？dinger方程的拟谱方法

本文讨论了广义混合非线性Schrodinger方程的周期初值问题,构造了守恒的半离散Fourier拟谱格式,对其近似解进行了先验估计,并证明了格式的收敛性.证明了该方程存在孤立子解,并给出其孤立子解的精确表达式.研究了线性化方程的稳定性问题,即在初值有扰动的情况下,该方程只有振荡解和鞍点.最后,通过数值例子验证了格式的可信性,数值计算表明,本格式时间方向可取大步长且是长时间稳定的,我们还计算了孤立子解,并绘出了在初值有扰动的情况下,相空间的轨线图.     

A numerical method for solving linear integral equations of the second kind on the non-rectangular domains based on the meshless method

The main purpose of this article is to describe a numerical scheme for solving two-dimensional linear Fredholm integral equations of the second kind on a non-rectangular domain. The method approximates the solution by the discrete collocation method based on radial basis functions (RBFs) constructed on a set of disordered data. The proposed method does not require any background mesh or cell structures, so it is meshless and consequently independent of the geometry of domain. This approach reduces the solution of the two-dimensional integral equation to the solution of a linear system of algebraic equations. The error analysis of the method is provided. The proposed scheme is also extended to linear mixed Volterra–Fredholm integral equations. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the new technique.     

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

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

广义对称正则长波方程的守恒差分格式

文章考虑了具有齐次边界条件的广义对称正则长波方程的有限差分格式.提出了一个守恒并且线性非耦合的三层有限差分格式,由于格式在计算中只需要解三对角线性方程组,从而避免了其中的迭代计算.文中先讨论了一个离散守恒量,然后我们利用离散泛函分析方法证明了格式的收敛性和稳定性,从理论上得到了收敛阶为O(h~2+τ~2).通过数值试验表明,所提的方法是可靠有效的.     

Two differential equation systems for equality-constrained optimization

This paper presents two differential systems, involving first and second order derivatives of problem functions, respectively, for solving equality-constrained optimization problems. Local minimizers to the optimization problems are proved to be asymptotically stable equilibrium points of the two differential systems. First, the Euler discrete schemes with constant stepsizes for the two differential systems are presented and their convergence theorems are demonstrated. Second, we construct algorithms in which directions are computed by these two systems and the stepsizes are generated by Armijo line search to solve the original equality-constrained optimization problem. The constructed algorithms and the Runge–Kutta method are employed to solve the Euler discrete schemes and the differential equation systems, respectively. We prove that the discrete scheme based on the differential equation system with the second order information has the locally quadratic convergence rate under the local Lipschitz condition. The numerical results given here show that Runge–Kutta method has better stability and higher precision and the numerical method based on the differential equation system with the second information is faster than the other one.     

一维Allen-Cahn方程紧差分格式的离散最大化原则和能量稳定性研究

本文主要研究相场模拟中的Allen-Cahn模型,考虑一维Allen-Cahn方程紧差分方法的数值逼近.建立具有O(∫τ2+h4)精度的全离散紧差分格式,证明在合理的步长比和时间步长的约束下,其数值解满足离散最大化原则,在此基础上,研究了全离散格式的能量稳定性.最后给出数值算例.     

First-Order and Second-Order,Chaos-Free,Finite Difference Schemes for Fisher Equation

A new class of finite difference schemes is constructed for Fisher partial differential equation i.e. the reaction-diffusion equation with stiff source term: $au(1-u)$. These schemes have the properties that they reduce to high fidelity algorithms in the diffusion-free case namely in which the numerical solutions preserve the properties in the exact solutions for arbitrary time step-size and reaction coefficient α＞0 and all nonphysical spurious solutions including bifurcations and chaos that normally appear in the standard discrete models of Fisher partial differential equation will not occur. The implicit schemes so developed obtain the numerical solutions by solving a single linear algebraic system at each step. The boundness and asymptotic behaviour of numerical solutions obtained by all these schemes are given. The approach constructing the above schemes can be extended to reaction-diffusion equations with other stiff source terms.     

A three‐level linearized finite difference scheme for the camassa–holm equation

The Camassa–Holm (CH) system is a strong nonlinear third‐order evolution equation. So far, the numerical methods for solving this problem are only a few. This article deals with the finite difference solution to the CH equation. A three‐level linearized finite difference scheme is derived. The scheme is proved to be conservative, uniquely solvable, and conditionally second‐order convergent in both time and space in the discrete L_∞ norm. Several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 451–471, 2014