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

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

两个带变时间步的两重网格算法的稳定性分析

对用于求解非线性发展方程的两个带变时间步的两重网格算法，对空间变量用有限元离散，对时间变量分别用一阶精度Euler显式和二阶精度半隐式差分格式离散，然后构造两重网格算法，通过深入的稳定性分析，得出本文的算法优于标准全离散有限元算法。     

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 ANALYSIS OF NONLINEAR GALERKIN ＆ GALERKIN METHODS FOR NONLINEAR EVOLUTION EQUATIONS

Abstract. Ogr object in this artlcle is to describe tbe Galerkln scheme and nonlin-eax Galerkin scheme for the approximation of nonlinear evolution equations, and tostudy the stability of these schemes. Spatial discretizatlon can be pedormed by eitherGalerkln spectral method or nonlinear Galerldn spectral method; time discretizatlort isdone hy Euler sin.heine wklch is explicit or implicit in the nonlinear terms. According tothe stability analysis of the above schemes, the stability of nonllneex Galerkln methodis better than that of Galexkln method.     

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

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

巴拿哈空间中随机年龄结构种群方程的数值分析

本文研究年龄结构随机种群方程的离散误差,在空间离散中用到Galerkin公式,时间离散中用到显式欧拉公式.     

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.     

Local discontinuous Galerkin approximation of non-Fickian diffusion model in viscoelastic polymers

In this paper, we investigate a fully discrete local discontinuous Galerkin approximation of a non-linear non-Fickian diffusion model in viscoelastic polymers. For the spatial discretization, we adopt local discontinuous Galerkin finite element method and for the time discretization we use backward Euler method. We derive the stability estimate and a priori error estimate for the discrete scheme. Numerical examples are given to verify the theoretical findings.     

Numerical analysis of semilinear stochastic evolution equations in Banach spaces

The solution of stochastic evolution equations generally relies on numerical computation. Here, usually the main idea is to discretize the SPDE spatially obtaining a system of SDEs that can be solved by e.g., the Euler scheme. In this paper, we investigate the discretization error of semilinear stochastic evolution equations in L^p-spaces, resp. Banach spaces. The space discretization may be done by Galerkin approximation, for the time discretization we consider the implicit Euler, the explicit Euler scheme and the Crank–Nicholson scheme. In the last section, we give some examples, i.e., we consider an SPDEs driven by nuclear Wiener noise approximated by wavelets and delay equation approximated by finite differences.     

Approximation for Semilinear Stochastic Evolution Equations

We investigate the approximation by space and time discretization of quasi linear evolution equations driven by nuclear or space time white noise. An error bound for the implicit Euler, the explicit Euler, and the Crank–Nicholson scheme is given and the stability of the schemes are considered. Lastly we give some examples of different space approximation, i.e., we consider approximation by eigenfunction, finite differences and wavelets.     

Wavelet-Taylor Galerkin Method for the Burgers Equation

In this paper, we propose a wavelet-Taylor Galerkin method for the numerical solution of the Burgers equation. In deriving the computational scheme, Taylor-generalized Euler time discretization is performed prior to wavelet-based Galerkin spatial approximation. The linear system of equations obtained in the process are solved by approximate-factorization-based simple explicit schemes, and the resulting solution is compared with that from regular methods. To deal with transient advection-diffusion situations that evolve toward a convective steady state, a splitting-up strategy is known to be very effective. So the Burgers equation is also solved by a splitting-up method using a wavelet-Taylor Galerkin approach. Here, the advection and diffusion terms in the Burgers equation are separated, and the solution is computed in two phases by appropriate wavelet-Taylor Galerkin schemes. Asymptotic stability of all the proposed schemes is verified, and the L_∞ errors relative to the analytical solution together with the numerical solution are reported. AMS subject classification (2000) 65M70     

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

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

First order decoupled method of the primitive equations of the ocean I: Time discretization

In this article, we study the time discretization of the first order decoupled semi-implicit scheme of the 3D primitive equations of the ocean in the case of the non-slip and non-heat flux boundary conditions on the side. The almost unconditional stability of the scheme and the optimal error estimates of the time discrete velocity, pressure and density are provided.     

A partitioned finite element scheme based on Gauge-Uzawa method for time-dependent MHD equations

In this paper, we mainly introduce a partitioned scheme based on Gauge-Uzawa finite element method for the 2D time-dependent incompressible magnetohydrodynamics (MHD) equations. It is a fully decoupled projection method which combines the Gauge and Uzawa methods within a variational formulation. Firstly, the temporal discretization is based on backward Euler technique for the linear term and semi-implicit scheme for the nonlinear term. Secondly, the spatial approximation of fluid velocity, hydrodynamic pressure, and magnetic field apply the mixed element method. Finally, the validity, reliability, and accuracy of the proposed algorithms are supported by numerical experiments.     

The composite Milstein methods for the numerical solution of Ito stochastic differential equations

In this paper, we present the composite Milstein methods for the strong solution of Ito stochastic differential equations. These methods are a combination of semi-implicit and implicit Milstein methods. We give a criterion for choosing either the implicit or the semi-implicit scheme at each step of our numerical solution. The stability and convergence properties are investigated and discussed for the linear test equation. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. The stability properties of the composite Milstein methods are found to be more superior compared to those of the Milstein, the Euler and even better than the composite Euler method. This superiority in stability makes the methods a better candidate for the solution of stiff SDEs.     

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     

A theory of discretization for nonlinear evolution inequalities applied to parabolic Signorini problems

We present a discretization theory for a class of nonlinear evolution inequalities that encompasses time dependent monotone operator equations and parabolic variational inequalities. This discretization theory combines a backward Euler scheme for time discretization and the Galerkin method for space discretization. We include set convergence of convex subsets in the sense of Glowinski-Mosco-Stummel to allow a nonconforming approximation of unilateral constraints. As an application we treat parabolic Signorini problems involving the p-Laplacian, where we use standard piecewise polynomial finite elements for space discretization. Without imposing any regularity assumption for the solution we establish various norm convergence results for piecewise linear as well piecewise quadratic trial functions, which in the latter case leads to a nonconforming approximation scheme. Entrata in Redazione il 16 marzo 1998, in versione riveduta il 15 febbraio 1999.     

High order accurate semi-implicit WENO schemes for hyperbolic balance laws

In this paper we propose a family of well-balanced semi-implicit numerical schemes for hyperbolic conservation and balance laws. The basic idea of the proposed schemes lies in the combination of the finite volume WENO discretization with Roe’s solver and the strong stability preserving (SSP) time integration methods, which ensure the stability properties of the considered schemes [S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001) 89-112]. While standard WENO schemes typically use explicit time integration methods, in this paper we are combining WENO spatial discretization with optimal SSP singly diagonally implicit (SDIRK) methods developed in [L. Ferracina, M.N. Spijker, Strong stability of singly diagonally implicit Runge-Kutta methods, Appl. Numer. Math. 58 (2008) 1675-1686]. In this way the implicit WENO numerical schemes are obtained. In order to reduce the computational effort, the implicit part of the numerical scheme is linearized in time by taking into account the complete WENO reconstruction procedure. With the proposed linearization the new semi-implicit finite volume WENO schemes are designed.A detailed numerical investigation of the proposed numerical schemes is presented in the paper. More precisely, schemes are tested on one-dimensional linear scalar equation and on non-linear conservation law systems. Furthermore, well-balanced semi-implicit WENO schemes for balance laws with geometrical source terms are defined. Such schemes are then applied to the open channel flow equations. We prove that the defined numerical schemes maintain steady state solution of still water. The application of the new schemes to different open channel flow examples is shown.     

A stabilized semi-implicit Galerkin scheme for Navier-Stokes equations

By introducing a time relaxation term for the time derivative of higher frequency components, we proposed a stabilized semi-implicit Galerkin scheme for evolutionary Navier-Stokes equations in this paper. Analysis shows that such a scheme has weaker stability conditions than that of a classical semi-implicit Galerkin scheme and, when a suitable relaxation parameter σ is chosen, it generates an approximate solution with the same accuracy as the classical one. That means the proposed scheme might use a larger time step to generate a bounded approximate solution. Thus it is more suitable for long time simulations.     

Fractional stochastic Burgers‐type equation in Hölder space—Wellposedness and approximations

In this work, we use the spectral Galerkin method to prove the existence of a pathwise unique mild solution of a fractional stochastic partial differential equation of Burgers type in a Hölder space. We get the temporal regularity, and using a combination of Galerkin and exponential‐Euler methods, we obtain a full discretization scheme of the solution. Moreover, we calculate the rates of convergence for both approximations (Galerkin and full discretization) with respect to time and to space.