Error estimation of spectral method for solving three-dimensional vorticity equation

Much work has been done for spectral scheme of P.D.E. (see [1]). Recently the author proposed a technique to prove the strict error estimation of spectral scheme for non-linear problems such as K.D.V.-Burgers' equation, two-dimensional vorticity equation and so on ([2]–[4]). In this paper we generalize this technique into three-dimensional vorticity equation. Under some conditions these error estimations imply convergence. The more smooth the solution of P.D.E., the more accurate the approximate solution.The author is     

A FOURIER-CHEBYSHEV SPECTRAL METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATIONS WITH UNILATERAL PERIODIC BOUNDARY CONDITIONS

A mixed Fourier-Chebyshev spectral method is constructed for three-dimensional vorticity equations with unilateral periodic boundary condition. The generalized stability and the convergence are analyzed. The optimal error estimation is given. The technique in this paper is also suitable for other nonlinear problems.     

The Fourier Pseudospectral Method for Two-Dimensional Vorticity Equations

In this paper we develop a Fourier pseudospectral method forsolving two-dimensional vorticity equations. We prove the generalizedstability of the schemes and give convergence estimations dependingon the smoothness of the solution of the vorticity equations. Spectral methods have been applied widely to the partial differentialequations of fluid dynamics [4–11]. Guo Ben-yu proposeda technique to estimate strictly the error of the spectral schemesfor the K.D.V.-Burgers equation, the two-dimensional vorticityequations, and the Navier-Stokes equations [5,6,8]. On the otherhand, the authors [7,10] developed a pseudospectral method byusing Riesz spherical means to get better results. In this paper,we generalize this method to two-dimensional vorticity equations.The generalized stability and the convergence are proved. Thenumerical results show the advantage of such a method.     

Spectral-finite element method for solving two-dimensional vorticity equations

In this paper, we construct a spectral-finite element scheme for solving semi-periodical two-dimensional vorticity equations. The error between the genuine solution and approximate solution is estimated strictly. The numerical results show the advantages of such a method. The technique used in this paper can be easily generalized to three-dimensional problems.     

A Spectral Approximation of the Barotropic Vorticity Equation

A spectral scheme is considered for solving the barotropic vorticity equation. The error estimates are proved strictly. The technique used in this paper is also useful for other nonlinear problems defined on a spherical surface.     

广义Benjamin-Bona-Mahony方程Fourier谱方法的大时间误差估计

§ 1　ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＢｅｎｊａｍｉｎ Ｂｏｎａ Ｍａｈｏｎｙｅｑｕａｔｉｏｎｕｔ+ｕｘ+ｕｕｘ -ｕｘｘ-ｕｘｘｔ =0 ( 1 .1 )ｉｎｃｏｒｐｏｒａｔｅｓｎｏｎｌｉｎｅａｒｄｉｓｐｅｒｓｉｖｅａｎｄｄｉｓｓｉｐａｔｉｖｅｅｆｆｅｃｔｓ ,ａｎｄｈａｓｂｅｅｎｐｒｏｐｏｓｅｄａｓａｍｏｄｅｌｆｏｒｂｏｔｈｔｈｅｂｏｒｅｐｒｏｐａｇａｔｉｏｎａｎｄｔｈｅｗａｔｅｒｗａｖｅｓ[1,2 ] .Ｔｈｅｅｘｉｓｔｅｎｃｅａｎｄｕｎｉｑｕｅｎｅｓｓｏｆｓｏｌｕｔｉｏｎｓｆｏｒｔｈｉｓｅｑｕａｔｉｏｎｈａｖｅｂｅｅ…     

The Large Time Convergence of Spectral Method for Generalized Kuramoto-Sivashinsky Equations

In this paper we use the spectral method to analyse the generalized Kuramoto-Sivashinsky equations. We prove the existence and uniqueness of global smooth solution of the equations. Finally, we obtain the error estimation between spectral approximate solution and exact solution on large time.     

双相滞热传导方程的有限元分析

本文考虑一类具有广泛应用背景的双相滞热传导方程混合边界问题．建立了其有限元和交替方向有限元的两种数值逼近格式．利用微分方程的先验估计理论与技巧,作出了数值解的L^2—范数估计结果．基于一系列的误差估计,也研究了两种逼近格式数值的稳定性和收敛性。     

Implicit a posteriori error estimates for the Maxwell equations

An implicit a posteriori error estimation technique is presented and analyzed for the numerical solution of the time-harmonic Maxwell equations using Nédélec edge elements. For this purpose we define a weak formulation for the error on each element and provide an efficient and accurate numerical solution technique to solve the error equations locally. We investigate the well-posedness of the error equations and also consider the related eigenvalue problem for cubic elements. Numerical results for both smooth and non-smooth problems, including a problem with reentrant corners, show that an accurate prediction is obtained for the local error, and in particular the error distribution, which provides essential information to control an adaptation process. The error estimation technique is also compared with existing methods and provides significantly sharper estimates for a number of reported test cases.

     

Analysis and convergence of the MAC scheme. II. Navier-Stokes equations

The MAC discretization scheme for the incompressible Navier-Stokes equations is interpreted as a covolume approximation to the equations. Using some results from earlier papers dealing with covolume error estimates for div-curl equation systems, and under certain conditions on the data and the solutions of the Navier-Stokes equations, we obtain first-order error estimates for both the vorticity and the pressure.

     

A meshless local Galerkin method for the numerical solution of Hammerstein integral equations based on the moving least squares technique

In this paper, a computational scheme is proposed to estimate the solution of one- and two-dimensional Fredholm-Hammerstein integral equations of the second kind. The method approximates the solution using the discrete Galerkin method based on the moving least squares (MLS) approach as a locally weighted least squares polynomial fitting. The discrete Galerkin technique for integral equations results from the numerical integration of all integrals in the system corresponding to the Galerkin method. Since the proposed method is constructed on a set of scattered points, it does not require any background meshes and so we can call it as the meshless local discrete Galerkin method. The implication of the scheme for solving two-dimensional integral equations is independent of the geometry of the domain. The new method is simple, efficient and more flexible for most classes of nonlinear integral equations. The error analysis of the method is provided. The convergence accuracy of the new technique is tested over several Hammerstein integral equations and obtained results confirm the theoretical error estimates.     

非线性KdV-Schr(o)dinger方程Fourier谱逼近的大时间性态

近几年来,对具弱阻尼的非线性发展方程的研究越来越受到人们的关注.大部分情况下,由于精确解无法得到,我们只有通过求数值解来研究方程解的性质.本文讨论具弱阻尼的非线性KdV-Schroedinger方程Fourier谱逼近的大时间性态问题.我们构造了方程的Fourier近似谱格式,并对方程的近似解作了相应的先验估计及方程近似解与精确解之间的误差估计.最后,证明了近似吸引子AN的存在性及其弱上半连续性dω(AN,A)→0.     

Uniformly convergent novel finite difference methods for singularly perturbed reaction–diffusion equations

In this paper, we construct a kind of novel finite difference (NFD) method for solving singularly perturbed reaction–diffusion problems. Different from directly truncating the high‐order derivative terms of the Taylor's series in the traditional finite difference method, we rearrange the Taylor's expansion in a more elaborate way based on the original equation to develop the NFD scheme for 1D problems. It is proved that this approach not only can highly improve the calculation accuracy but also is uniformly convergent. Then, applying alternating direction implicit technique, the newly deduced schemes are extended to 2D equations, and the uniform error estimation based on Shishkin mesh is derived, too. Finally, numerical experiments are presented to verify the high computational accuracy and theoretical prediction.     

一类半线性抛物型方程全离散Chebyshev拟谱逼近的大时间性态

In the paper, the nonperidic initial value problem for a class of semilinear parabolic equations is considered. We construct the full discrete Chebyshev pseudospectral scheme and analyze the error of approximate solution for it. We obtain the error estimation on large time using the local continuation method and the existence of approximate global attractor.     

Global Error Estimation with Runge--Kutta Methods

An analysis of global error estimation for Runge—Kuttasolutions of ordinary differential equations is presented. Thebasic technique is that of Zadunaisky in which the global erroris computed from a numerical solution of a neighbouring problemrelated to the main problem by some method of interpolation.It is shown that Runge—Kutta formulae which permit validglobal error estimation using low-degree interpolation can bedeveloped, thus leading to more accurate and computationallyconvenient algorithms than was hitherto expected. Some specialRunge—Kutta processes up to order 4 are presented togetherwith numerical results.     

二维不可压缩流的谱方法

本文以二维涡度方程为模型,介绍了谱方法和拟谱方法以及它们与差分方法和有限元法相结合的混合解法．这些方法可推广应用于其它一些类似的非线性问题．本文还给出了这些方法的某些数值例子和误差估计结果     

The Crank–Nicolson finite spectral element method and numerical simulations for 2D non-stationary Navier–Stokes equations

In this paper, we first build a semi-discretized Crank–Nicolson (CN) model about time for the two-dimensional (2D) non-stationary Navier–Stokes equations about vorticity–stream functions and discuss the existence, stability, and convergence of the time semi-discretized CN solutions. And then, we build a fully discretized finite spectral element CN (FSECN) model based on the bilinear trigonometric basic functions on quadrilateral elements for the 2D non-stationary Navier–Stokes equations about the vorticity–stream functions and discuss the existence, stability, and convergence of the FSECN solutions. Finally, we utilize two sets of numerical experiments to check out the correctness of theoretical consequences.     

On the application of mosaic-skeleton approximations of matrices for the acceleration of computations in the vortex method for the three-dimensional Euler equations

We study the possibility of using fast matrix multiplication methods for the approximation of the velocity field when solving the system of differential equations describing the vorticity transport in an ideal incompressible fluid in Lagrangian coordinates. We suggest a numerical scheme that permits effectively using the fast matrix multiplication (the method of mosaic-skeleton approximations). We show that the functions used for the computation of the velocity field and moving grids appearing in the solution of the problem permit one to use the above-mentioned method. We prove the convergence of the resulting numerical solution to the exact solution with regard of the error contributed by the use of the algorithm for approximate fast multiplication of matrices by vectors.     

Convergence of Fourier spectral method for resonant long-short nonlinear wave interaction

In this paper, the evolution equations with nonlinear term describing the resonance interaction between the long wave and the short wave are studied. The semi-discrete and fully discrete Crank-Nicholson Fourier spectral schemes are given. An energy estimation method is used to obtain error estimates for the approximate solutions. The numerical results obtained are compared with exact solution and found to be in good agreement.     

On using the Lagrange coefficients for a posteriori error estimation

A posteriori error estimation of the objective functional is considered by means of a differential presentation of a finite-difference scheme and adjoint equations. The local approximation error is presented as a Taylor series remainder in the Lagrangian form. The field of the Lagrange coefficients is determined by a high-accuracy finite-difference template affecting the computation results. The feasibility of using the Lagrange coefficients for refining the solution and for estimating its uncertainty is considered.