High Order Padé Schemes for Nonlinear Wave Propagation Problems

Split-step Padémethod and split-step fourier method are applied to the higher- order nonlinear Schrdinger equation.It is proved that a combination of Padé scheme and spectral method is the most effective method,which has a spectral-like resolution and good stability nature.In particular,we propose an unconditional stable implicit Padé scheme to solve odd order nonlinear equations.Numerical results demonstrate the excellent performance of Padé schemes for high order nonlinear equations.     

ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS

This paper discusses the accelerating of nonlinear parabolic equations. Two iterative methods for solving the implicit scheme new nonlinear iterative methods named by the implicit-explicit quasi-Newton （IEQN） method and the derivative free implicit-explicit quasi-Newton （DFIEQN） method are introduced, in which the resulting linear equations from the linearization can preserve the parabolic characteristics of the original partial differential equations. It is proved that the iterative sequence of the iteration method can converge to the solution of the implicit scheme quadratically. Moreover, compared with the Jacobian Free Newton-Krylov （JFNK） method, the DFIEQN method has some advantages, e.g., its implementation is easy, and it gives a linear algebraic system with an explicit coefficient matrix, so that the linear （inner） iteration is not restricted to the Krylov method. Computational results by the IEQN, DFIEQN, JFNK and Picard iteration methods are presented in confirmation of the theory and comparison of the performance of these methods.     

A FINITE DIFFERENCE SCHEME FOR SOLVING THE NONLINEAR POISSON-BOLTZMANN EQUATION MODELING CHARGED SPHERES

In this work, we propose an efficient numerical method for computing the electrostaticinteraction between two like-charged spherical particles which is governed by the nonlinearPoisson-Boltzmann equation. The nonlinear problem is solved by a monotone iterativemethod which leads to a sequence of linearized equations. A modified central finite differ-ence scheme is developed to solve the linearized equations on an exterior irregular domainusing a uniform Cartesian grid. With uniform grids, the method is simple, and as aconsequence, multigrid solvers can be employed to speed up the convergence. Numericalexperiments on cases with two isolated spheres and two spheres confined in a chargedcylindrical pore are carried out using the proposed method. Our numerical schemes arefound efficient and the numerical results are found in good agreement with the previouspublished results.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

ASYMPTOTIC BEHAVIOR OF SOLUTIONS TO THE HEAT EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS

The purpose of this paper is to investigate the stability and asymptotic behav-ior of the time-dependent solutions to a linear parabolic equation with nonlinear boundarycondition in relation to their corresponding steady state solutions. Then, the above resultsare extended to a semilinear parabolic equation with nonlinear boundary condition by an-alyzing the corresponding eigenvalue problem and using the method of upper and lowersolutions.     

A MONOTONE COMPACT IMPLICIT SCHEME FOR NONLINEAR REACTION-DIFFUSION EQUATIONS

A monotone compact implicit finite difference scheme with fourth-order accuracy in space and second-order in time is proposed for solving nonlinear reaction-diffusion equations. An accelerated monotone iterative method for the resulting discrete problem is presented. The sequence of iteration converges monotonically to the unique solution of the discrete problem, and the convergence rate is either quadratic or nearly quadratic, depending on the property of the nonlinear reaction. The numerical results illustrate the high accuracy of the proposed scheme and the rapid convergence rate of.the iteration.     

New numerical methods for the coupled nonlinear Schrödinger equations

In this paper, three numerical schemes with high accuracy for the coupled Schrodinger equations are studied. The conserwtive properties of the schemes are obtained and the plane wave solution is analysised. The split step Runge-Kutta scheme is conditionally stable by linearized analyzed. The split step compact scheme and the split step spectral method are unconditionally stable. The trunction error of the schemes are discussed. The fusion of two solitions colliding with different β is shown in the figures. The numerical experments demonstrate that our algorithms are effective and reliable.     

HIGH-ORDER NONLINEAR GALERKIN METHODS

The nonlinear Galerkin methods are numerical schemes well adapted to the long-term integration of nonlinear evolution partial differential equations. In this paper, a class of high-order nonlinear Galerkin methods are provided. Moreover, convergence results with high-order spectral accuracy are derived for the schemes introduced.     

AD GALERKIN ANALYSIS FOR NONLINEAR PSEUDO—HYPERBOLIC EQUATIONS

AD(Alternating direction)Galerkin schemes for d-dimensional nonlinear pseudo-hyperbolic equations are studied.By using patch approximation technique,AD procedure is realized,and calculation,work is simplified.By using Galerkin approach,highly computational accuracy is kept.By using various priori estimate techniques for differential equations,difficulty coming form non-linearity is treated,and optimal H^1 and L^2 convergence prop-erties are demonstrated.Moreover,although all the existed AD Galerkin schemes using patch approximation are limited to have only one order accuracy in time increment,yet the schemes formulated in this paper have second order accuracy in it.This implies an essential advancement in AD Galerkin aualysis.     

A FOURIER SPECTRAL SCHEME FOR SOLVING NONLINEAR KLEIN-GORDON EQUATION

A Fourier spectral scheme is proposed for solving the periodic problem of nonlinear Klein-Gordon equation. Its stability and convergence are investigated. Numerical results are also presented.     

Convergence and stability of finite element nonlinear Galerkin method for the Navier-Stokes equations

In this article, we present a new fully discrete finite element nonlinear Galerkin method, which are well suited to the long time integration of the Navier-Stokes equations. Spatial discretization is based on two-grid finite element technique; time discretization is based on Euler explicit scheme with variable time step size. Moreover, we analyse the boundedness, convergence and stability condition of the finite element nonlinear Galerkin method. Our discussion shows that the time step constraints of the method depend only on the coarse grid parameter and the time step constraints of the finite element Galerkin method depend on the fine grid parameter under the same convergence accuracy. Received February 2, 1994 / Revised version received December 6, 1996     

Stability of Galerkin and Inertial Algorithms with variable time step size

This paper provides Galerkin and Inertial Algorithms for solving a class of nonlinear evolution equations. Spatial discretization can be performed by either spectral or finite element methods; time discretization is done by Euler explicit or Euler semi-implicit difference schemes with variable time step size. Moreover, the boundedness and stability of these algorithms are studied. By comparison, we find that the boundedness and stability of Inertial Algorithm are superior to the ones of Galerkin Algorithm in the case of explicit scheme and the boundedness and stability of two algorithms are same in the case of semi-implicit scheme.     

Stability and convergence of optimum spectral non‐linear Galerkin methods

Our objective in this article is to present some numerical schemes for the approximation of the 2‐D Navier–Stokes equations with periodic boundary conditions, and to study the stability and convergence of the schemes. Spatial discretization can be performed by either the spectral Galerkin method or the optimum spectral non‐linear Galerkin method; time discretization is done by the Euler scheme and a two‐step scheme. Our results show that under the same convergence rate the optimum spectral non‐linear Galerkin method is superior to the usual Galerkin methods. Finally, numerical example is provided and supports our results. Copyright © 2001 John Wiley & Sons, Ltd.     

非定常 N-S 方程的非线性 Galerkin 算法及其误差估计

本文给出了二维非定常N-S方程的三种数值格式,其中空间变量用谱非线性Galerkin算法进行离散,时间变量用有限差分离散,并研究了这些格式数值解的逼近精度．最后,给出了部分数值计算结果．     

两类变时间步长的非线性Galerkin算法的稳定性

1．引言近年来,随着计算机的飞速发展,人们越来越关心非线性发展方程解的渐进行为．为了较精确地描述解在时间t→∞时的渐进行为,人们发展了一类惯性算法,即非线性Galerkin算法．该算法是将来解空间分解为低维部分和高维部分,相应的方程可以分别投影到它们上面,它的解也相应地分解为两部分,大涡分量和小涡分量;然后核算法给出大涡分量和小涡分量之间依赖关系的一种近似,以便容易求出相应的近似解．许多研究表明,非线性Galerkin算法比通常的Galerkin算法节省可观的计算量．当数值求解微分方程时,计算机只能对已知数据进行有限位…     

The spectral method for the flow between two concentric rotating spheres, part II: The full discrete nonlinear Legendre-Galerkin method

The nonlinear Galerkin methods are numerical schemes for evolutionary partial differential equations based on the theory of inertial manifolds and approximate inertial manifolds. In this paper, we consider the flow between two concentric rotating spheres, and combine the Legendre-Galerkin spectral methods in Part I together with the nonlinear Galerkin method, then construct the full discrete nonlinear Legendre-Galerkin spectral scheme, and derive the stability conditions and its error estimate.     

Multi-level spectral galerkin method for the navier-stokes problem I : spatial discretization

A multi-level spectral Galerkin method for the two-dimensional non-stationary Navier-Stokes equations is presented. The method proposed here is a multiscale method in which the fully nonlinear Navier-Stokes equations are solved only on a low-dimensional space subsequent approximations are generated on a succession of higher-dimensional spaces j=2, . . . ,J, by solving a linearized Navier-Stokes problem around the solution on the previous level. Error estimates depending on the kinematic viscosity 0<ν<1 are also presented for the J-level spectral Galerkin method. The optimal accuracy is achieved when We demonstrate theoretically that the J-level spectral Galerkin method is much more efficient than the standard one-level spectral Galerkin method on the highest-dimensional space . The work of this author was supported in part by the NSF of China 10371095, City University of Hong Kong Research Project 7001093 Hong Kong and the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 1084/02P)     

Uniform stability of spectral nonlinear Galerkin methods

This article provides a stability analysis for the backward Euler schemes of time discretization applied to the spatially discrete spectral standard and nonlinear Galerkin approximations of the nonstationary Navier‐Stokes equations with some appropriate assumption of the data (λ, u₀, f). If the backward Euler scheme with the semi‐implicit nonlinear terms is used, the spectral standard and nonlinear Galerkin methods are uniform stable under the time step constraint Δt ≤ (2/λλ₁). Moreover, if the backward Euler scheme with the explicit nonlinear terms is used, the spectral standard and nonlinear Galerkin methods are uniform stable under the time step constraints Δt = O(λ) and Δt = O(λ), respectively, where λ ≤ λ, which shows that the restriction on the time step of the spectral nonlinear Galerkin method is less than that of the spectral standard Galerkin method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

非线性Galerkin算法的收敛性和复杂性

本文给出了数值求解非线性发展方程的全离散非线性Ｇａｌｅｒｋｉｎ算法,即将空间离散时的谱非线性Ｇａｌｅｒｋｉｎ算法和时间离散的Ｅｕｌｅｒ差分格式相结合,得到了显式和隐式两种全离散数值格式,相应地也考虑了显式和隐式的Ｇａｌｅｒｋｉｎ全离散格式,并分别分析了上述四种全离散格式的收敛性和复杂性,经过比较得出结论;在某些约束条件下,非线性Ｇａｌｅｒｋｉｎ算法和Ｇａｌｅｒｋｉｎ算法具有相同阶的收敛速度,然而前     

Numerical analysis of a modified finite element nonlinear Galerkin method

Summary. A fully discrete modified finite element nonlinear Galerkin method is presented for the two-dimensional equation of Navier-Stokes type. The spatial discretization is based on two finite element spaces X_H and X_h defined on a coarse grid with grid size H and a fine grid with grid size h << H, respectively; the time discretization is based on the Euler explicit scheme with respect to the nonlinear term. We analyze the stability and convergence rate of the method. Comparing with the standard finite element Galerkin method and the nonlinear Galerkin method, this method can admit a larger time step under the same convergence rate of same order. Hence this method can save a large amount of computational time. Finally, we provide some numerical tests on this method, the standard finite element Galerkin method, and the nonlinear Galerkin method, which are in a good agreement with the theoretical analysis.Mathematics Subject Classification (2000): 35Q30, 65M60, 65N30, 76D05