一维Maxwell方程间断解的多区域Legendre-Galerkin谱方法

研究了非一致介质一维Maxwell方程间断问题的多区间Legendre谱方法,建立了半离散和全离散Crank-Nicolson格式,设计了并行算法,并分析了方法的稳定性和收敛性.数值算例显示了多区间Legendre谱方法对于间断问题的有效性.     

一类非线性Schrodinger方程的差分／Legendre谱元法

考察一类带幂次非线性项的Schrodinger方程的Dirichlet初边值问题，提出了一个有效的计算格式，其中时间方向上应用了一种守恒的二阶差分隐格式，空间方向上采用Legendre谱元法．对于时间半离散格式，证职了该格式具有能量守恒性质，并给出了L^2误差估计，对于全离散格式，应用不动点原理证明了数值解的存在唯一性，并给出了L^2误差估计．最后，通过数值试验验证了结果的可信性．     

一维Burgers方程和KdV方程的广义有限谱方法

给出了高精度的广义有限谱方法．为使方法在时间离散方面保持高精度,采用了Adams-Bashforth 预报格式和Adams-Moulton校正格式,为了避免由Korteweg-de Vries(KdV)方程的弥散项引起的数值振荡, 给出了两种数值稳定器．以Legendre多项式、Chebyshev多项式和Hermite多项式为基函数作为例子,给出的方法与具有分析解的Burgers方程的非线性对流扩散问题和KdV方程的单孤独波和双孤独波传播问题进行了比较,结果非常吻合．     

一类非线性Scllr(o)dillger方程的差分/Legendre谱元法

考察一类带幂次非线性项的Schrodinger方程的Dirichlet初边值问题,提出了一个有效的计算格式,其中时间方向上应用了一种守恒的二阶差分隐格式,空间方向上采用Legendre谱元法.对于时间半离散格式,证明了该格式具有能量守恒性质,并给出了L2误差估计,对于全离散格式,应用不动点原理证明了数值解的存在唯一性,并给出了L2误差估计.最后,通过数值试验验证了结果的可信性.     

长短波方程的高阶紧致格式（英文）

在空间方向用高阶紧致格式离散,时间方向分别用CNI格式、Richardson格式和分裂步CNI格式离散,得到了长短波方程的一些数值格式.这些格式在时间方向是二阶收敛的,空间方向是四阶的,而用到的模版与二阶中心差分格式是一样的.数值结果表明,与中心格式相比,新提出的格式较已有格式计算效率更高.同时,从数值结果可以猜测CNI格式和分裂步CNI格式能够保持原问题的一些守恒量.     

广义对称正则长波方程的勒让德和切贝雪夫拟谱方法

本文考虑了具齐次边界条件的广义对称正则长波方程的Legendre和Chebyshev拟谱方法，构造了半离散和全离散的Legendre和Chebyshev拟谱格式，从理论上得到了这些格式对应的最优误差估计。     

二维阻尼非线性sine-Gordon方程的共形多辛Fourier拟谱格式

为二维阻尼非线性sine-Gordon方程构造了一个新的共形多辛Fourier拟谱格式.基于原系统的共形多辛哈密尔顿形式,首先在时间和空间方向上分别用辛中点和Fourier拟谱方法进行离散,得到一个全离散格式.随后证明了构造的格式保持离散的共形多辛守恒律.最后数值实验验证了格式的有效性.     

分数阶Burgers方程的有限元计算(英文)

建立了一维和二维分数阶Burgers方程的有限元格式.时间分数阶导数使用L1方法离散,空间方向使用有限元方法离散.通过选择合适的基函数,将离散后的方程转化成一个非线性代数方程组,并应用牛顿迭代方法求解.数值实验显示出了方法的有效性.     

发展型Navier-Stokes方程基于压力稳定化方法的谱逼近格式

本文采用压力稳定化方法近似模拟不可压缩条件,进而构造了发展型非周期NavierStokes方程的全离散Legendre谱逼近计算格式,严格分析了格式的广义稳定性与收敛性．本文建立的逼近结果也适用其它非周期问题．     

运动壁面槽道流动的直接数值模拟

采用谱方法,在曲线坐标系下对不可压缩Newton流体的N-S方程进行求解,采用定义在物理空间中的流动物理量以避免使用协变、逆变形式的控制方程．在计算空间采用Fourier-Chebyshev谱方法进行空间离散,时间推进采用高精度时间分裂法．为了减小时间分裂带来的误差,采用了高精度的压力边界条件．与其他求解协变、逆变形式控制方程的谱方法相比,该方法在保持谱精度的同时减小了计算量．首先通过静止波形壁面和行波壁面槽道湍流的直接数值模拟,对该数值方法进行了验证;其次,作为初步应用,利用该方法研究了槽道湍流中周期振动凹坑所产生的流动结构．     

Mixed Laguerre–Legendre Spectral Method for Incompressible Flow in an Infinite Strip

This paper concerns the mixed Laguerre–Legendre spectral approximation and its application to numerical simulation of incompressible flow in an infinite strip. Some approximation results in weighted Sobolev spaces are given. A Laguerre–Legendre spectral scheme for the stream function form of Navier–Stokes equations is constructed. The stability and the convergence of the proposed scheme are proved. The numerical experiments show the high accuracy of this method. The main techniques used in this paper are also applicable to other nonlinear partial differential equations in an infinite strip.     

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.     

Parallel spectral-element direction splitting method for incompressible Navier–Stokes equations

An efficient parallel algorithm for the time dependent incompressible Navier–Stokes equations is developed in this paper. The time discretization is based on a direction splitting method which only requires solving a sequence of one-dimensional Poisson type equations at each time step. Then, a spectral-element method is used to approximate these one-dimensional problems. A Schur-complement approach is used to decouple the computation of interface nodes from that of interior nodes, allowing an efficient parallel implementation. The unconditional stability of the full discretized scheme is rigorously proved for the two-dimensional case. Numerical results are presented to show that this algorithm retains the same order of accuracy as a usual spectral-element projection type schemes but it is much more efficient, particularly on massively parallel computers.     

Spectral element solution of steady incompressible viscous free-surface flows

In this paper we present a new approach for the solution of the steady incompressible Navier-Stokes equations in a domain bounded in part by a free surface. In the spatial discretization procedure, a Legendre spectral element method is used to generate the discrete equations. For effective solution of the set of algebraic equations, the geometry is decoupled from the fluid velocity and pressure. In addition, two different algorithms are proposed depending on the importance of surface tension effects. Numerical results are presented to demonstrate the effectiveness of the proposed algorithms.     

Klein-Gordon-Zakharov方程初边值问题的Legendre谱方法

研究Klein-Gordon-Zakharov方程初边值问题的Legendre谱方法.在先验估计的基础上,证明了该格式的稳定性和收敛性,并得到最优阶误差估计.另外,还设计了一个半隐格式,并给出数值例子.在文章的后面给出了多区域谱格式,数值结果表明精度要高于单区域.     

基于Picard迭代的PN×PN-2谱元法求解定常不可压缩Navier-Stokes方程

该文给出了一种求解二维定常不可压缩Navier-Stokes方程的基于Picard线性化迭代的P_N×P_N-2谱元法.通过Picard线性化将不可压缩Navier-Stokes方程的求解转化为一系列线性的Stokes-type方程,再利用非交错网格的P_N×P_N-2谱元法计算每个迭代步的Stokes-type方程.为了消除伪压力模,压力离散比速度离散低两阶,非交错网格的应用使得方程的离散方便且不会带来相应的插值误差,从而保证了谱精度.通过此方法数值计算了有精确解的Stokes流动、Kovasznay流动和方腔顶盖驱动流,结果表明,迭代收敛非常快,误差收敛达到了谱精度收敛,并且避免了压力震荡的出现,表明了该文方法准确可靠.     

Error analysis of a fully discrete projection method for magnetohydrodynamic system

In this paper, we develop and analyze a finite element projection method for magnetohydrodynamics equations in Lipschitz domain. A fully discrete scheme based on Euler semi-implicit method is proposed, in which continuous elements are used to approximate the Navier–Stokes equations and H ( curl ) conforming Nédélec edge elements are used to approximate the magnetic equation. One key point of the projection method is to be compatible with two different spaces for calculating velocity, which leads one to obtain the pressure by solving a Poisson equation. The results show that the proposed projection scheme meets a discrete energy stability. In addition, with the help of a proper regularity hypothesis for the exact solution, this paper provides a rigorous optimal error analysis of velocity, pressure and magnetic induction. Finally, several numerical examples are performed to demonstrate both accuracy and efficiency of our proposed scheme.     

Analysis of the Yosida method for the incompressible Navier–Stokes equations

The Yosida method was introduced in (Quarteroni et al., to appear) for the numerical approximation of the incompressible unsteady Navier–Stokes equations. From the algebraic viewpoint, it can be regarded as an inexact factorization of the matrix arising from the space and time discretization of the problem. However, its differential interpretation resides on an elliptic stabilization of the continuity equation through the Yosida regularization of the Laplacian (see (Brezis, 1983, Ciarlet and Lions, 1991)). The motivation of this method as well as an extensive numerical validation were given in (Quarteroni et al., to appear).In this paper we carry out the analysis of this scheme. In particular, we consider a first-order time advancing unsplit method. In the case of the Stokes problem, we prove unconditional stability and moreover that the splitting error introduced by the Yosida scheme does not affect the overall accuracy of the solution, which remains linear with respect to the time step. Some numerical experiments, for both the Stokes and Navier–Stokes equations, are presented in order to substantiate our theoretical results.     

DFT modal analysis of spectral element methods for the 2D elastic wave equation

The DFT modal analysis is a dispersion analysis technique that transforms the equations of a numerical scheme to the discrete Fourier transform domain sampled in the mesh nodes. This technique provides a natural matching of exact and approximate modes of propagation. We extend this technique to spectral element methods for the 2D isotropic elastic wave equation, by using a Rayleigh quotient approximation of the eigenvalue problem that characterizes the dispersion relation, taking full advantage of the tensor product representation of the spectral element matrices. Numerical experiments illustrate the dependence of dispersion errors on the grid resolution, polynomial degree, and discretization in time. We consider spectral element methods with Chebyshev and Legendre collocation points.