A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation ...A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.展开更多
采用Mortar有限单元法(mortar finite element method,MFEM)能够得到正定、对称的系数矩阵,而且刚度矩阵是分块对称的,这种特点适合于并行迭代求解。阐述了非重叠Mortar有限单元法(non-overlapping MFEM,NO-MFEM)的基本原理,介绍了适合...采用Mortar有限单元法(mortar finite element method,MFEM)能够得到正定、对称的系数矩阵,而且刚度矩阵是分块对称的,这种特点适合于并行迭代求解。阐述了非重叠Mortar有限单元法(non-overlapping MFEM,NO-MFEM)的基本原理,介绍了适合于NO-MFEM并行计算的区域分解策略以及并行求解的基本流程。针对简单2维静电场问题,使用NO-MFEM进行了并行计算,并与理论值和串行计算结果进行对比,验证了所提方法的有效性。同时,对于非协调网格造成的计算误差进行了分析。NO-MFEM法的并行计算为工程应用中优化设计问题的区域分解和并行求解提供了一种新的选择。展开更多
Mortar元法(mortar element method,MEM)是一种新型区域分解算法,它允许将求解区域分解为多个子域,在各个区域以最适合子域特征的方式离散。在各个区域的交界面上,边界节点不要求逐点匹配,而是通过建立加权积分形式的Mortar条件使得交...Mortar元法(mortar element method,MEM)是一种新型区域分解算法,它允许将求解区域分解为多个子域,在各个区域以最适合子域特征的方式离散。在各个区域的交界面上,边界节点不要求逐点匹配,而是通过建立加权积分形式的Mortar条件使得交界面上的传递条件在分布意义上满足。Mortar有限元法(mortar finite element method,MFEM)将MEM和有限元法(finite element method,FEM)相结合,在各区域中分别使用FEM网格离散,区域的交界面上通过施加Mortar条件实现区域间的自由度连续。该文阐述了非重叠Mortar有限单元法(non-overlapping MFEM,NO-MFEM)的基本原理,介绍了NO-MFEM的程序实现过程,使用NO-MFEM对2维静磁场问题和3维静电场问题进行了计算,并与FEM模型结果进行对比,验证了该文方法的有效性。将NO-MFEM应用于电磁分析,丰富了电磁场数值计算理论,为运动涡流问题和大规模问题的分析提供了新的选择。展开更多
The lowest order H1-Galerkin mixed finite element method （for short MFEM） is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas elem...The lowest order H1-Galerkin mixed finite element method （for short MFEM） is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O（h2）/O（h2 ＋ r2） in Hi-norm and H（div; Ω）-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.展开更多
In this paper, a fully discrete scheme based on the LI approximation in temporal direction for the fractional derivative of order in (0,1) and nonconforming mixed finite element method (MFEM) in spatial direction is e...In this paper, a fully discrete scheme based on the LI approximation in temporal direction for the fractional derivative of order in (0,1) and nonconforming mixed finite element method (MFEM) in spatial direction is established. First, we prove a novel result of the consistency error estimate with order O(h^2)of EQ1^rot element (see Lemma 2.3). Then, by using the proved character of EQ1^rot element, we present the superconvergent estimates for the original variable u in the broken H^1-norm and the flux →p =△u in the (L^2)^2-norm under a weaker regularity of the exact solution. Finally, numerical results are provided to confirm the theoretical analysis.展开更多
An H^1-Galerkin mixed finite element method (MFEM) is discussed for the two-dimensional Ginzburg-Landau equation with the bilinear element and zero order Raviart-Thomas elemen t (Q11+Q10×Qo01). A linearized Crank...An H^1-Galerkin mixed finite element method (MFEM) is discussed for the two-dimensional Ginzburg-Landau equation with the bilinear element and zero order Raviart-Thomas elemen t (Q11+Q10×Qo01). A linearized Crank-Nicolson fully-discrete scheme is developed and a time-discrete system is introduced to split the error into two parts which are called the temporal error and the spatial error, respectively. On one hand, the regularity of the time-discrete system is deduced through the temporal error estimation. On the other hand, the superconvergent estimates of u in H^1-norm and →q in H(div;Ω)-norm with order 0(h^2+τ^2) are obtained unconditionally based on the achievement of the spatial result. At last, a numerical experiment is included to illustrate the feasibility of the proposed method. Here, h is the subdivision parame ter and τ is the time step.展开更多
In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,t...In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,the superclose error estimates for both the velocity in broken H-norm and the pressure in L2-norm are first obtained,which play a key role to bound the numerical solution in Lx-norm.Then the corresponding global superconvergence results are derived through a suitable interpolation postprocessing approach.Finally,some numerical results are provided to demonstrated the theoretical analysis.展开更多
In this paper,we consider the mixed finite element method(MFEM)of the elasticity problem in two and three dimensions(2D and 3D).We develop a new residual based stabilization method to overcome the inf-sup difficulty,a...In this paper,we consider the mixed finite element method(MFEM)of the elasticity problem in two and three dimensions(2D and 3D).We develop a new residual based stabilization method to overcome the inf-sup difficulty,and use Langrange elements to approximate the stress and displacement.The new method is unconditionally stable,and its stability can be obtained directly from C´ea’s lemma.Optimal error estimates for the H1-norm of the displacement and H(div)-norm of the stress can be obtained at the same time.Numerical results show the excellent stability and accuracy of the new method.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.10971203,11271340,and 11101381)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20094101110006)
文摘A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.
文摘采用Mortar有限单元法(mortar finite element method,MFEM)能够得到正定、对称的系数矩阵,而且刚度矩阵是分块对称的,这种特点适合于并行迭代求解。阐述了非重叠Mortar有限单元法(non-overlapping MFEM,NO-MFEM)的基本原理,介绍了适合于NO-MFEM并行计算的区域分解策略以及并行求解的基本流程。针对简单2维静电场问题,使用NO-MFEM进行了并行计算,并与理论值和串行计算结果进行对比,验证了所提方法的有效性。同时,对于非协调网格造成的计算误差进行了分析。NO-MFEM法的并行计算为工程应用中优化设计问题的区域分解和并行求解提供了一种新的选择。
文摘Mortar元法(mortar element method,MEM)是一种新型区域分解算法,它允许将求解区域分解为多个子域,在各个区域以最适合子域特征的方式离散。在各个区域的交界面上,边界节点不要求逐点匹配,而是通过建立加权积分形式的Mortar条件使得交界面上的传递条件在分布意义上满足。Mortar有限元法(mortar finite element method,MFEM)将MEM和有限元法(finite element method,FEM)相结合,在各区域中分别使用FEM网格离散,区域的交界面上通过施加Mortar条件实现区域间的自由度连续。该文阐述了非重叠Mortar有限单元法(non-overlapping MFEM,NO-MFEM)的基本原理,介绍了NO-MFEM的程序实现过程,使用NO-MFEM对2维静磁场问题和3维静电场问题进行了计算,并与FEM模型结果进行对比,验证了该文方法的有效性。将NO-MFEM应用于电磁分析,丰富了电磁场数值计算理论,为运动涡流问题和大规模问题的分析提供了新的选择。
基金Supported in part by the National Natural Science Foundation of China under Grant Nos.11671369,11271340the Natural Science Foundation of the Education Department of Henan Province under Grant Nos.14A110009,16A110022
文摘The lowest order H1-Galerkin mixed finite element method （for short MFEM） is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart- Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O（h2）/O（h2 ＋ r2） in Hi-norm and H（div; Ω）-norm axe deduced for the semi-discrete and the fully-discrete schemes, where h, r- denote the mesh size and the time step, respectively, which improve the results in the previous literature.
基金the National Natural Science Foundation of China (No. 1167136911271340).
文摘In this paper, a fully discrete scheme based on the LI approximation in temporal direction for the fractional derivative of order in (0,1) and nonconforming mixed finite element method (MFEM) in spatial direction is established. First, we prove a novel result of the consistency error estimate with order O(h^2)of EQ1^rot element (see Lemma 2.3). Then, by using the proved character of EQ1^rot element, we present the superconvergent estimates for the original variable u in the broken H^1-norm and the flux →p =△u in the (L^2)^2-norm under a weaker regularity of the exact solution. Finally, numerical results are provided to confirm the theoretical analysis.
基金Natural Science Foundation of China (Grant Nos. 11671369, 11271340).
文摘An H^1-Galerkin mixed finite element method (MFEM) is discussed for the two-dimensional Ginzburg-Landau equation with the bilinear element and zero order Raviart-Thomas elemen t (Q11+Q10×Qo01). A linearized Crank-Nicolson fully-discrete scheme is developed and a time-discrete system is introduced to split the error into two parts which are called the temporal error and the spatial error, respectively. On one hand, the regularity of the time-discrete system is deduced through the temporal error estimation. On the other hand, the superconvergent estimates of u in H^1-norm and →q in H(div;Ω)-norm with order 0(h^2+τ^2) are obtained unconditionally based on the achievement of the spatial result. At last, a numerical experiment is included to illustrate the feasibility of the proposed method. Here, h is the subdivision parame ter and τ is the time step.
基金This work is supported by National Natural Science Foundation of China(Nos.11671369,11271340).
文摘In this paper,the superconvergence properties of the time-dependent Navier-Stokes equations are investigated by a low order nonconforming mixed finite element method(MFEM).In terms of the integral identity technique,the superclose error estimates for both the velocity in broken H-norm and the pressure in L2-norm are first obtained,which play a key role to bound the numerical solution in Lx-norm.Then the corresponding global superconvergence results are derived through a suitable interpolation postprocessing approach.Finally,some numerical results are provided to demonstrated the theoretical analysis.
基金This work was supported by High-Level Personal Foundation of Henan University of Technology(No.2015BS018)National Natural Science Fund of China(Nos.11671369,11271340 and 11601124).
文摘In this paper,we consider the mixed finite element method(MFEM)of the elasticity problem in two and three dimensions(2D and 3D).We develop a new residual based stabilization method to overcome the inf-sup difficulty,and use Langrange elements to approximate the stress and displacement.The new method is unconditionally stable,and its stability can be obtained directly from C´ea’s lemma.Optimal error estimates for the H1-norm of the displacement and H(div)-norm of the stress can be obtained at the same time.Numerical results show the excellent stability and accuracy of the new method.
