THE TWO-LEVEL STABILIZED FINITE ELEMENT METHOD BASED ON MULTISCALE ENRICHMENT FOR THE STOKES EIGENVALUE PROBLEM

In this paper,we first propose a new stabilized finite element method for the Stokes eigenvalue problem.This new method is based on multiscale enrichment,and is derived from the Stokes eigenvalue problem itself.The convergence of this new stabilized method is proved and the optimal priori error estimates for the eigenfunctions and eigenvalues are also obtained.Moreover,we combine this new stabilized finite element method with the two-level method to give a new two-level stabilized finite element method for the Stokes eigenvalue problem.Furthermore,we have proved a priori error estimates for this new two-level stabilized method.Finally,numerical examples confirm our theoretical analysis and validate the high effectiveness of the new methods.     

解Stokes特征值问题的一种两水平稳定化有限元方法

基于局部Gauss积分,研究了解Stokes特征值问题的一种两水平稳定化有限元方法．该方法涉及在网格步长为H的粗网格上解一个Stokes特征值问题,在网格步长为h=O(H2)的细网格上解一个Stokes问题．这样使其能够仍旧保持最优的逼近精度,求得的解和一般的稳定化有限元解具有相同的收敛阶,即直接在网格步长为h的细网格上解一个Stokes特征值问题．因此,该方法能够节省大量的计算时间．数值试验验证了理论结果．     

Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered.This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h = O(H 2),which can still maintain the asymptotically optimal accuracy.It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution,which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h.Hence,the two-level stabilized finite element method can save a large amount of computational time.Moreover,numerical tests confirm the theoretical results of the present method.     

Darcy-Brinkman方程的残量型后验误差估计

有限元后验误差估计为有限元后处理技术提供了有效的理论分析.在有限元后处理方法中,基于残量型的后验误差方法在计算科学和数值模拟中占据着重要的地位,它对于局部奇异问题有着很好的逼近效果,整体上降低了数值模拟计算时的矩阵规模,使得计算资源能更合理的分配.本文设计了Darcy-Brinkman方程的残量型后验误差估计子,并给出了估计子的下界,通过数值算例验证该方法的有效性.     

不同角度Y型狭窄血管血液流动的数值模拟

根据粘性不可压Navier-Stokes方程,建立Y型分又血管中血液流动的数学模型,进而采用有限元方法研究不同分又角Y型血管动脉狭窄位置对血液流动的影响,得到了不同角度不同狭窄位置和无狭窄病变时的数值模拟结果,主要给出了各种情况下血液流动的流线图和压力图.一方面,观察流线图可知,血液流经狭窄区域时,出现流动分离,并在一定区域产生涡流,且随着狭窄位置不同,涡流位置和涡流区域面积也随之不同;另一方面,从计算的压力图中可以看到当血液流过狭窄区域时,压力发生迅速变化,且相同分叉角度下狭窄位置不同,狭窄区域压力不同;狭窄位置相同时,不同分叉角度的血管分又区域压力也有差别.