THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

LOW ORDER NONCONFORMING RECTANGULAR FINITE ELEMENT METHODS FOR DARCY-STOKES PROBLEMS

In this paper, we consider lower order rectangular finite element methods for the singularly perturbed Stokes problem. The model problem reduces to a linear Stokes problem when the perturbation parameter is large and degenerates to a mixed formulation of Poisson＇s equation as the perturbation parameter tends to zero. We propose two 2D and two 3D nonconforming rectangular finite elements, and derive robust discretization error estimates. Numerical experiments are carried out to verify the theoretical results.     

混合有限元非嵌套网格型预处理法

1 引言 众所周知,二阶椭圆型问题混合有限元离散以后的矩阵是不定的,所以对混合法很难形成一种有效的区域分解法,在文[9]、[10]、[11]中提出了一些混合有限元方法的区域分解法,但在实际计算中有很多局限性。最近Chen对混合有限元法提出一种全新的解释并把它应用到多重网格法中,他的基本思想是混合有限元离散的代数系统实际上等价于某个非协调有限元离散的代数系统,这样可把一个不定问题转化为一个正定问题,本文将基于这种思想考虑混合有限元的区域分解法。 若按传统的Dryia-widlund两水平加性Schwarz方法,要求两层网格间具有嵌套关系,这样在应用中将带来很大的不便。本文将不要求粗网格嵌入细网格中,减少两层网格间的     

SPECTRAL AND SPECTRAL ELEMENT METHODS FOR HIGH ORDER PROBLEMS WITH MIXED BOUNDARY CONDITIONS

In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and prove their spectral accuracy by using the recent results on the Jacobi quasi-orthogonal approximation. Numerical results demonstrate the high accuracy of suggested algorithm, which also works well even for oscillating solutions.     

PRECONDITIONING HIGHER ORDER FINITE ELEMENT SYSTEMS BY ALGEBRAIC MULTIGRID METHOD OF LINEAR ELEMENTS

We present and analyze a robust preconditioned conjugate gradient method for the higher order Lagrangian finite element systems of a class of elliptic problems. An auxiliary linear element stiffness matrix is chosen to be the preconditioner for higher order finite elements. Then an algebraic multigrid method of linear finite element is applied for solving the preconditioner. The optimal condition number which is independent of the mesh size is obtained. Numerical experiments confirm the efficiency of the algorithm.     

MIXED FINITE ELEMENT METHODS FOR A STRONGLY NONLINEAR PARABOLIC PROBLEM

1.IntroductionFOrsecondorderellipticproblems,themixedmethodwasdescribedandanalyzedbymanyauthors[1--3]inthecaseoflinearequationsindivergenceform,aswellasin[4,5]forquasilinearornonlinearproblemsindivergenceform.JohnsonandThorn6e[6]consideredalternativeproofSofthepreviouslyknownerrorestimatesforsuchmethodsintheellipticcase.Theyalsoanalyzedthemixedfiniteelementmethodforthepaxabolicequationgivenbypt--ho~f.Garcia[7]studiedtheconvergenceofmixedfiniteelementapproximationstoquasilinearparabolicequati…     

MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

A PRECONDITIONER FOR COUPLING SYSTEM OF NATURAL BOUNDARY ELEMENT AND COMPOSITE GRID FINITE ELEMENT

1. IntroductionThe coupling of boundary elemellts and finite elements is of great imPortance for the nu-mercal treatment of boundary value problems posed on unbounded domains. It permits us tocombine the advanages of boundary elements for treating domains extended to infinity withthose of finite elemenis in treating the comP1icated bounded domains.The standard procedure of coupling the boundary elemeni and finite elemeni methods isdescribed as follows. First, the (unbounded) domain is divided…     

UNIFORMLY-STABLE FINITE ELEMENT METHODS FOR DARCY-STOKES-BRINKMAN MODELS

In this paper, we consider 2D and 3D Darcy-Stokes interface problems. These equations are related to Brinkman model that treats both Darcy＇s law and Stokes equations in a single form of PDE but with strongly discontinuous viscosity coefficient and zerothorder term coefficient. We present three different methods to construct uniformly stable finite element approximations. The first two methods are based on the original weak formulations of Darcy-Stokes-Brinkman equations. In the first method we consider the existing Stokes elements. We show that a stable Stokes element is also uniformly stable with respect to the coefficients and the jumps of Darcy-Stokes-Brinkman equations if and only if the discretely divergence-free velocity implies almost everywhere divergence-free one. In the second method we construct uniformly stable elements by modifying some well-known H（div）-conforming elements. We give some new 2D and 3D elements in a unified way. In the last method we modify the original weak formulation of Darcy-Stokes- Brinkman equations with a stabilization term. We show that all traditional stable Stokes elements are uniformly stable with respect to the coefficients and their jumps under this new formulation.     

二阶椭圆问题的一类广义有限元法

本文提出了求解二阶椭圆问题的一类广义有限元方法,分析了广义有限元方法的优越性,证明了二阶椭圆问题的广义有限元方法具有比标准的Galerkin有限元方法更高阶的收敛速度,根据插值算子的性质,进一步证明了有限元解的亏量迭代校正收敛到广义有限元解,并用数值例子说明广义有限元方法是有效的.     

COUPLING OF FINITE ELEMENT AND BOUNDARY ELEMENT METHODS FOR THE SCATTERING BY PERIODIC CHIRAL STRUCTURES

Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.     

Stokes方程的稳定化间断有限元法

本文对定常的Stokes方程提出了一种新的间断有限元法,通过对通常的间断Galerkin有限元法应用稳定化思想,建立了一个相容的稳定间断有限元格式,对速度和压力的任意分片多项式空间Pl(K),Pm(K)的间断有限元逼近证明了解的存在唯一性,给出了关于速度和压力的L2 范数的最优误差估计．     

求解三维高次拉格朗日有限元方程的代数多重网格法

本文针对带有间断系数的三维椭圆问题,讨论任意四面体剖分下的二次拉格朗日有限元方程的代数多重网格法．通过分析线性和高次有限元空间之间的关系,我们给出了一种新的网格粗化算法和构造提升算子的代数途径．进一步,我们还对新的代数多重网格法给出了收敛性分析．数值实验表明这种代数多重网格法对求解二次拉格朗日有限元方程是健壮和有效的。     

一类四阶抛物型积分-微分方程的混合间断时空有限元法

构造四阶抛物型积分-微分方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明离散解的稳定性,存在唯一性和收敛性.     

哈密尔顿系统的有限元法

利用常微分方程的连续有限元法,结合函数的M-型展开,对非线性哈密尔顿系统证明了连续一、二次有限元分在3阶量、5阶量意义下近似保辛,且保持能量守恒.在数值实验中结合庞加莱截面,哈密尔顿混沌数值试验结果与理论相吻合.     

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＦｉｎｉｔｅｖｏｌｕｍｅｅｌｅｍｅｎｔｍｅｔｈｏｄｏｒＦＶＥＭｕｓｅｓａｖｏｌｕｍｅｉｎｔｅｇｒａｌｆｏｒｍｕｌａｔｉｏｎｏｆｔｈｅｄｉｆ ｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｉｔｈａｆｉｎｉｔｅｐａｒｔｉｔｉｏｎｉｎｇｓｅｔｏｆｖｏｌｕｍｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｅｑｕａｔｉｏｎ ,ｔｈｅｎｒｅ ｓｔｒｉｃｔｓｔｈｅａｄｍｉｓｓｉｂｌｅｆｕｎｃｔｉｏｎｓｔｏａｌｉｎｅａｒｆｉｎｉｔｅｅｌｅｍｅｎｔｓｐａｃｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｓｏｌｕｔｉｏｎ ([1 -4 ] ) .…     

SEMI-DISCRETE AND FULLY DISCRETE PARTIALPROJECTION FINITE ELEMENT METHODS FOR THEVIBRATING TIMOSHENKO BEAM

1.IntroductionTheTimoshenkobeammodelisgivenbywherethebeamisconsidereddamped,drepresentsthebeamthicknessandI~[0,1].000istherotationofverticalfibersinthebeamandw(x)istheverticaldisplacementofthebeam'scenterline(underaverticalloadgivenbyg(x)).Analogoust...     

HIGH ORDER FINITE DIFFERENCE/SPECTRAL METHODS TO A WATER WAVE MODEL WITH NONLOCAL VISCOSITY

In this paper, efficient numerical schemes are proposed for solving the water wave model with nonlocal viscous term that describe the propagation of surface water wave. By using the Caputo fractional derivative definition to approximate the nonlocal fractional operator, finite difference method in time and spectral method in space are constructed for the considered model. The proposed method employs known 5/2 order scheme for fractional derivative and a mixed linearization for the nonlinear term. The analysis shows that the proposed numerical scheme is unconditionally stable and error estimates are provided to predict that the second order backward differentiation plus 5/2 order scheme converges with order 2 in time, and spectral accuracy in space. Several numerical results are provided to verify the efficiency and accuracy of our theoretical claims. Finally, the decay rate of solutions is investigated.