Analysis of a symplectic difference scheme for a coupled nonlinear Schrödinger system

In general, proofs of convergence and stability are difficult for symplectic schemes of nonlinear equations. In this paper, a symplectic difference scheme is proposed for an initial-boundary value problem of a coupled nonlinear Schrödinger system. An important lemma and an induction argument are used to prove the unique solvability, convergence and stability of numerical solutions. An iterative algorithm is also proposed for the symplectic scheme and its convergence is proved. Numerical examples show the efficiency of the symplectic scheme and the correction of our numerical analysis.     

Perturbation of Volterra Difference Equations

A fixed point theorem is used to investigate nonlinear Volterra difference equations that are perturbed versions of linear equations. Sufficient conditions are established to ensure that the stability properties of linear Volterra difference equations are preserved under perturbation. The existence of asymptotically periodic solutions of perturbed Volterra difference equations is also proved.     

A new conservative fourth‐order accurate difference scheme for solving a model of nonlinear dispersive equations

This article is devoted to the study of a nonlinear conservative fourth‐order difference scheme for a model of nonlinear dispersive equations that is governed by the RLW‐KdV equation. The existence of the approximate solution and the convergence of the difference scheme are proved, by using the energy method. In addition, the convergent order in maximum norm is 2 in temporal direction and 4 in spatial direction. The unconditional stability as well as uniqueness of the difference scheme is also derived. An application on the RLW and MRLW equations is discussed numerically in details. Furthermore, interaction of solitary waves with different amplitudes are shown. The 3 invariants of the motion are evaluated to determine the conservation proprieties of the system. The temporal evaluation of a Maxwellian initial pulse is then studied. Some numerical examples are given to validate the theoretical results.     

Conservative schemes for the symmetric Regularized Long Wave equations

In this paper, we study the Symmetric Regularized Long Wave (SRLW) equations by finite difference method. We design some numerical schemes which preserve the original conservative properties for the equations. The first scheme is two-level and nonlinear-implicit. Existence of its difference solutions are proved by Brouwer fixed point theorem. It is proved by the discrete energy method that the scheme is uniquely solvable, unconditionally stable and second-order convergent for U in L_∞ norm, and for N in L₂ norm on the basis of the priori estimates. The second scheme is three-level and linear-implicit. Its stability and second-order convergence are proved. Both of the two schemes are conservative so can be used for long time computation. However, they are coupled in computing so need more CPU time. Thus we propose another three-level linear scheme which is not only conservative but also uncoupled in computation, and give the numerical analysis on it. Numerical experiments demonstrate that the schemes are accurate and efficient.     

Study of a Generalized Nonlinear Euler-Poisson-Darboux System: Numerical and Bessel Based Solutions

In this paper a nonlinear Euler-Poisson-Darboux system is considered. In a first part, we proved the genericity of the hypergeometric functions in the development of exact solutions for such a systemin some special cases leading to Bessel type differential equations. Next, a finite difference scheme in two-dimensional case has been developed. The continuous system is transformed into an algebraic quasi linear discrete one leading to generalized Lyapunov-Sylvester operators. The discrete algebraic system is proved to be uniquely solvable, stable and convergent based on Lyapunov criterion of stability and Lax-Richtmyer equivalence theorem for the convergence. A numerical example has been provided at the end to illustrate the efficiency of the numerical scheme developed in section 3. The present method is thus proved to be more accurate than existing ones and lead to faster algorithms.     

Difference Schemes on Nonuniform Grids for the Two-Dimensional Convection–Diffusion Equation

New second-order accurate monotone difference schemes on nonuniform spatial grids for two-dimensional stationary and nonstationary convection–diffusion equations are proposed. The monotonicity and stability of the solutions of the computational methods with respect to the boundary conditions, the initial condition, and the right-hand side are proved. Two-sided and corresponding a priori estimates are obtained in the grid norm of C. The convergence of the proposed algorithms to the solution of the original differential problem with the second order is proved.     

一维非线性Schrödinger 方程的两个无条件收敛的守恒紧致差分格式

本文对一维非线性 Schrödinger 方程给出两个紧致差分格式, 运用能量方法和两个新的分析技 巧证明格式关于离散质量和离散能量守恒, 而且在最大模意义下无条件收敛. 对非线性紧格式构造了 一个新的迭代算法, 证明了算法的收敛性, 并在此基础上给出一个新的线性化紧格式. 数值算例验证 了理论分析的正确性, 并通过外推进一步提高了数值解的精度.     

ON LINEARIZED FINITE DIFFERENCE SIMULATION FOR THE MODEL OF NUCLEAR REACTOR DYNAMICS

1 hoeductIOuThe dynamics models of one--dimensional continuous medium nuclear reactor are the foelowing initial--boundary value problem of the formsubject to the innal conditionsand the boundary conditionsIn (1. 1), x denotes position along the reactor, which is regarded as a rod of length L, t denotes the time, u(t) the logarithm of the loud reactor POwer, v(x,t) the deviation of the temperature from equilibrium, a(x) the ratio of the temperature coefficient of reactivity to theynean life of…     

Numerical approximation of solution for the coupled nonlinear Schrödinger equations

In this article, a compact finite difference scheme for the coupled nonlinear Schrödinger equations is studied. The scheme is proved to conserve the original conservative properties. Unconditional stability and convergence in maximum norm with order O(τ² + h⁴) are also proved by the discrete energy method. Finally, numerical results are provided to verify the theoretical analysis.     

A finite difference scheme for solving the undamped Timoshenko beam equations with both ends free

In this paper, a boundary feedback system of a class of non-uniform undamped Timoshenko beam with both ends free is considered. A linearized three-level difference scheme for the Timoshenko beam equations is derived by the method of reduction of order on uniform meshes. The unique solvability, unconditional stability and convergence of the difference scheme are proved by the discrete energy method. The convergence order in maximum norm is of order two in both space and time. The validity of this theoretical analysis is verified experimentally.     

河床流体模型方程的稳定化有限差分法及非线性稳定性分析

1 引言本文考虑如下问题: μ(x+2π,t)=μ(x,t), x∈R,t∈[0,τ], (1.2) μ(x,0) =μ_0(x) β,ε,σ∈R,ε,σ>0. (1.3) 该模型描述河床流体流动,其中μ(x,t)为实值函数,它代表河床流体中微粒沉淀(concen—tration)在空间方向上的周期小扰动。G.H.Ganser和D.A.Drew用摄动法对该问题进行了分析,认为该问题是非线性不稳定的。 数值研究表明,对该问题,采用通常的差分方法和Galerkin有限元是不稳定的。文     

Estimates for the solutions of systems of linear algebraic equations and their application to the investigation of the convergence of the method of finite differences

One proves theorems regarding the estimates of the solutions of systems of linear algebraic equations and inequalities. On the basis of these theorems one suggests a method for estimating the norm of the inverse matrix of a system of difference equations which approximates a boundary-value problem for an integrodifferential equation. The method allows us to eliminate the restrictions which are usually imposed on the coefficients of the integrodifferential equation in order to ensure the diagonal dominance in the system of the difference equations. One considers an application to nonlinear problems.     

General difference schemes with intrinsic parallelism for nonlinear parabolic systems

The boundary value problem for nonlinear parabolic system is solved by the finite difference method with intrinsic parallelism. The existence of the discrete vector solution for the general finite difference schemes with intrinsic parallelism is proved by the fixed-point technique in finite-dimensional Euclidean space. The convergence and stability theorems of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. The limitation vector function is just the unique generalized solution of the original problem for the parabolic system.     

一类变系数微分方程差分格式的收敛性(英文)

本文首先对一类变系数微分方程建立有限差分格式.然后利用矩阵的特征值和范数理论,讨论该格式解的收敛性和唯一性.通过数值算例,说明该格式既有效又便于模拟.并且文中所用方法还能用于高阶微分方程和某些非线性微分方程问题的研究.     

Entropy dissipative one‐leg multistep time approximations of nonlinear diffusive equations

New one‐leg multistep time discretizations of nonlinear evolution equations are investigated. The main features of the scheme are the preservation of the non‐negativity and the entropy dissipation structure of the diffusive equations. The key ideas are to combine Dahlquist's G‐stability theory with entropy dissipation methods and to introduce a nonlinear transformation of variables, which provides a quadratic structure in the equations. It is shown that G‐stability of the one‐leg scheme is sufficient to derive discrete entropy dissipation estimates. The general result is applied to a cross‐diffusion system from population dynamics and a nonlinear fourth‐order quantum diffusion model, for which the existence of semidiscrete weak solutions is proved. Under some assumptions on the operator of the evolution equation, the second‐order convergence of solutions is shown. Moreover, some numerical experiments for the population model are presented, which underline the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1119–1149, 2015     

Moving finite element methods for time fractional partial differential equations

With the aim of simulating the blow-up solutions, a moving finite element method, based on nonuniform meshes both in time and in space, is proposed in this paper to solve time fractional partial differential equations (FPDEs). The unconditional stability and convergence rates of 2-α for time and r for space are proved when the method is used for the linear time FPDEs with α-th order time derivatives. Numerical examples are provided to support the theoretical findings, and the blow-up solutions for the nonlinear FPDEs are simulated by the method.     

The discontinuous Galerkin method for fractional degenerate convection-diffusion equations

We propose and study discontinuous Galerkin methods for strongly degenerate convection-diffusion equations perturbed by a fractional diffusion (Lévy) operator. We prove various stability estimates along with convergence results toward properly defined (entropy) solutions of linear and nonlinear equations. Finally, the qualitative behavior of solutions of such equations are illustrated through numerical experiments.     

THE UNCONDITIONAL CONVERGENT DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR QUASILINEAR PARABOLIC SYSTEMS WITH TWO DIMENSIONS

In the present work we are going to solve the boundary value problem for the quasilinear parabolic systems of partial differential equations with two space dimensions by the finite difference method with intrinsic parallelism.Some fundamental behaviors of general finite difference schemes with intrinsic parallelism for the mentioned problems are studied.By the method of a priori estimation of the discrete solutions of the nonlinear difference systems,and the interpolation formulas of the various norms of the discrete functions and the fixed-point technique in finite dimensional Euclidean space,the existennce of the discrete vector solutions of the nonliear difference system with intrinsic parallelism are proved .Moreover the convergence of the discrete vector solutions of these difference schemes to the unique generalizd solution of the original quasilinear parabolic problem is proved.     

On a decoupled linear FEM integrator for eddy-current-LLG

We propose a numerical integrator for the coupled system of the eddy-current equation with the nonlinear Landau–Lifshitz–Gilbert equation. The considered effective field contains a general field contribution, and we particularly cover exchange, anisotropy, applied field and magnetic field (stemming from the eddy-current equation). Even though the considered problem is nonlinear, our scheme requires only the solution of two linear systems per time-step. Moreover, our algorithm decouples both equations so that in each time-step, one linear system is solved for the magnetization, and afterwards one linear system is solved for the magnetic field. Unconditional convergence – at least of a subsequence – towards a weak solution is proved, and our analysis even provides existence of such weak solutions. Numerical experiments with micromagnetic benchmark problems underline the performance and the stability of the proposed algorithm.     

An Unconditionally Stable Difference Approximation for a Class of Nonlinear Dispersive Equations

An unconditionally stable leapfrog finite difference scheme for a class of nonlinear dispersive equations is presented and analyzed. The solvability of the difference equation which is a tridiagonal circular linear system is discussed. Moreover, the convergence and stability of the difference scheme are also investigated by a standard argument so that more difficult priori estimations are avoided. Finally, numerical examples are given.