非线性薛定谔方程的平均向量场方法

利用平均向量场方法(AVF)对非线性薛定谔方程进行求解, 在理论上得到了一个保非线性薛定谔方程描述的系统能量守恒的AVF格式, 再分别用非线性薛定谔方程的AVF格式和辛格式数值模拟孤立波的演化行为, 并比较两个格式是否保系统能量守恒特性. 数值结果表明, AVF格式也能很好地模拟孤立波的演化行为,并且比辛格式更能保持系统的能量守恒.     

A local iteration scheme for nonlinear two-dimensional steady-state heat conduction: a BEM approach

An efficient algorithm is proposed to solve the steady-state nonlinear heat conduction equation using the boundary element method (BEM). Nonlinearity of the heat conduction equation arises from nonlinear boundary conditions and temperature dependence of thermal conductivity. Using Kirchhoff's transformation, the case of temperature dependence of thermal conductivity can be transformed to the nonlinear boundary conditions case. Applying the BEM technique, the resulting matrix equation becomes nonlinear. The nonlinearity, however, only involves the boundary nodes that have nonlinearboundary conditions. The proposed local iterative scheme reduces the entire BEM matrix equation to a smaller matrix equation whose rank is the same as the number of boundary nodes with nonlinear boundary conditions. The Newton-Raphson iteration scheme is used to solve the reduced nonlinear matrix equation. The local iterative scheme is first applied to two one-dimensional problems (analytical solutions are possible) with different nonlinear boundary conditions. It is then applied to a two-region problem. Finally, the local iterative scheme is applied to two cavity problems in which radiation plays a role in the heat transfer.     

The Conservation and Convergence of Two Finite Difference Schemes for KdV Equations with Initial and Boundary Value Conditions

Korteweg-de Vries equation is a nonlinear evolutionary partial differential equation that is of third order in space. For the approximation to this equation with the initial and boundary value conditions using the finite difference method, the difficulty is how to construct matched finite difference schemes at all the inner grid points. In this paper, two finite difference schemes are constructed for the problem. The accuracy is second-order in time and first-order in space. The first scheme is a two-level nonlinear implicit finite difference scheme and the second one is a three-level linearized finite difference scheme. The Browder fixed point theorem is used to prove the existence of the nonlinear implicit finite difference scheme. The conservation, boundedness, stability, convergence of these schemes are discussed and analyzed by the energy method together with other techniques. The two-level nonlinear finite difference scheme is proved to be unconditionally convergent and the three-level linearized one is proved to be conditionally convergent. Some numerical examples illustrate the efficiency of the proposed finite difference schemes.     

On local convergence of a symmetric semi-discrete scheme for an abstract analogue of the Kirchhoff equation

The present work considers a nonlinear abstract hyperbolic equation with a self-adjoint positive definite operator, which represents a generalization of the Kirchhoff string equation. A symmetric three-layer semi-discrete scheme is constructed for an approximate solution of a Cauchy problem for this equation. Value of the gradient in the nonlinear term of the scheme is taken at the middle point. It makes possible to find an approximate solution at each time step by inverting the linear operator. Local convergence of the constructed scheme is proved. Numerical calculations for different model problems are carried out using this scheme.     

广义Boussinesq方程的多辛方法

广义Boussinesq方程作为一类重要的非线性方程有着许多有趣的性质,基于Hamilton空间体系的多辛理论研究了广义Boussinesq方程的数值解法,构造了一种等价于多辛Box格式的新隐式多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律．对广义Boussinesq方程孤子解的数值模拟结果表明,该多辛离散格式具有较好的长时间数值稳定性．     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

Convergence of a splitting method for the nonlinear Boltzmann equation

Convergence of a splitting method scheme for the nonlinear Boltzmann equation is considered. Using the splitting method scheme, boundedness of the positive solutions in a space of continuous functions is obtained. By means of the solution boundedness and some a priori estimates, convergence of the splitting method scheme and uniqueness of the limiting element are proved. The limiting element satisfies an equivalent integral Boltzmann equation. Thereby global in time solvability of the nonlinear Boltzmann equation is shown.     

关于广义KPP方程的数值解

广义KPP(Kolmogorov-Petrovskii-Piskunov)方程是一个积分微分方程.为了要研究其数值解,我们首先将该方程转化为一个非线性双曲型方程,然后构造了一个线性化的差分格式,得到了差分格式解的存在唯一性,利用能量不等式证明了差分格式二阶收敛性和关于初值的无条件稳定性,数值结果验证了本文提出的方法.     

Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions

The nonlinear Klein–Gordon equation is used to model many nonlinear phenomena. In this paper, we propose a numerical scheme to solve the one-dimensional nonlinear Klein–Gordon equation with quadratic and cubic nonlinearity. Our scheme uses the collocation points and approximates the solution using Thin Plate Splines (TPS) radial basis functions (RBF). The implementation of the method is simple as finite difference methods. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme.     

一个解KdV方程的满足两个守恒律的差分格式

Korteweg-de Vries(KdV)方程是人们在研究一些物理问题时得到的非线性波 动方程,其解满足无穷多个守恒律．本文为该方程设计了一种差分格式,其采用的是有限 体积法．但与传统的有限体积法不同的是,它的数值解同时满足两个相关的守恒律．这样 可以更好地保持解的物理上的守恒性质．数值例子表明这一算法是有效的．     

A Nonlinear Difference Scheme and Inverse Scattering

A nonlinear partial difference equation is obtained and solved by the method of inverse scattering. In a certain continuum limit it is shown how this equation approximates the nonlinear Schrodinger equation and a related nonlinear differential-difference equation. At all times the solutions can be compared, and the scheme is shown to be convergent. These ideas apply to other nonlinear evolution equations as well.     

A Legendre–Galerkin Chebyshev collocation method for the Burgers equation with a random perturbation on boundary condition

In the paper, we apply the generalized polynomial chaos expansion and spectral methods to the Burgers equation with a random perturbation on its left boundary condition. Firstly, the stochastic Galerkin method combined with the Legendre–Galerkin Chebyshev collocation scheme is adopted, which means that the original equation is transformed to the deterministic nonlinear equations by the stochastic Galerkin method and the Legendre–Galerkin Chebyshev collocation scheme is used to deal with the resulting nonlinear equations. Secondly, the stochastic Legendre–Galerkin Chebyshev collocation scheme is developed for solving the stochastic Burgers equation; that is, the stochastic Legendre–Galerkin method is used to discrete the random variable meanwhile the nonlinear term is interpolated through the Chebyshev–Gauss points. Then a set of deterministic linear equations can be obtained, which is in contrast to the other existing methods for the stochastic Burgers equation. The mean square convergence of the former method is analyzed. Numerical experiments are performed to show the effectiveness of our two methods. Both methods provide alternative approaches to deal with the stochastic differential equations with nonlinear terms.     

PARALLEL COMPUTING METHOD OF PURE ALTERNATIVE SEGMENT EXPLICIT-IMPLICIT DIFFERENCE SCHEME FOR NONLINEAR LELAND EQUATION

The research on the numerical solution of the nonlinear Leland equation has important theoretical significance and practical value. To solve nonlinear Leland equation, this paper offers a class of difference schemes with parallel nature which are pure alternative segment explicit-implicit(PASE-I) and implicit-explicit(PASI-E) schemes. It also gives the existence and uniqueness,the stability and the error estimate of numerical solutions for the parallel difference schemes. Theoretical analysis demonstrates that PASE-I and PASI-E schemes have obvious parallelism, unconditionally stability and second-order convergence in both space and time. The numerical experiments verify that the calculation accuracy of PASE-I and PASI-E schemes are better than that of the existing alternating segment Crank-Nicolson scheme, alternating segment explicit-implicit and implicit-explicit schemes. The speedup of PASE-I scheme is 9.89, compared to classical Crank-Nicolson scheme. Thus the schemes given by this paper are high efficient and practical for solving the nonlinear Leland equation.     

A numerical scheme based on a solution of nonlinear advection–diffusion equations

A mathematical formulation of the two-dimensional Cole–Hopf transformation is investigated in detail. By making use of the Cole–Hopf transformation, a nonlinear two-dimensional unsteady advection–diffusion equation is transformed into a linear equation, and the transformed equation is solved by the spectral method previously proposed by one of the authors. Thus a solution to initial value problems of nonlinear two-dimensional unsteady advection–diffusion equations is derived. On the base of the solution, a numerical scheme explicit with respect to time is presented for nonlinear advection–diffusion equations. Numerical experiments show that the present scheme possesses the total variation diminishing properties and gives solutions with good quality.     

A complex tanh-function method applied to nonlinear equations of Schrödinger type

A complex tanh-function method is introduced for constructing exact travelling wave solutions of nonlinear partial differential equations with complex phases and solutions. The scheme is implemented for obtaining multiple soliton solutions to the nonlinear cubic Schrödinger equation and a generalized Schrödinger-like equation. In additon. an ansätze is proposed to obtain stationary soliton solutions of the cubic Schrödinger equation.     

非线性Schrdinger方程的高精度守恒差分格式

本文首先分析线性Schrodinger方程一种高阶差分格式的构造方法,得到方程的耗散项.在此基础上对三次非线性Schrodinger方程,提出了一种精度为O(r2 h2)的差分格式,证明了该格式保持了连续方程的两个守恒量,且是收敛的与稳定的.并通过数值例子与已有隐格式进行了比较,结果表明,本文格式在计算量类似的情况下,提高了数值精度.     

A second order numerical scheme for the solution of the one-dimensional Boussinesq equation

A predictor–corrector (P-C) scheme is applied successfully to a nonlinear method arising from the use of rational approximants to the matrix-exponential term in a three-time level recurrence relation. The resulting nonlinear finite-difference scheme, which is analyzed for local truncation error and stability, is solved using a P-C scheme, in which the predictor and the corrector are explicit schemes of order 2. This scheme is accelerated by using a modification (MPC) in which the already evaluated values are used for the corrector. The behaviour of the P-C/MPC schemes is tested numerically on the Boussinesq equation already known from the bibliography free of boundary conditions. The numerical results are derived for both the bad and the good Boussinesq equation and conclusions from the relevant known results are derived.      

一类多维非线性反应扩散方程有限差分格式的稳定性问题(英文)

本文研究了一类多维线性反应扩散方程差分格式的稳定性.利用量未知元方法,建立了具有增量未知元的有限差分格式;然后利用非线性Galerkin方法,得到该差分格式的稳定性条件.通过对该格式的稳定性分析,说明和经典的差分格式的稳定性相比较,带有增量未知元的有限差分格式的稳定性得到了提高.     

守恒型扩散方程非线性离散格式的性质分析和快速求解

对于守恒型扩散方程,研究其二阶时间精度非线性全隐有限差分离散格式的性质,证明了其解的存在唯一性.研究了二阶时间精度的Picard-Newton迭代格式,证明了迭代解对原问题真解的二阶时间和空间收敛性,以及对非线性离散解的二次收敛速度,实现了非线性问题的快速求解.本文中方法也适用于一阶时间精度格式的分析,并可推广至对流扩散问题.数值实验验证了二阶时间精度Picard-Newton迭代格式的高精度和高效率.     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.