Navier-Stokes方程的一种并行两水平有限元方法

基于区域分解技巧,提出了一种求解定常Navier-Stokes方程的并行两水平有限元方法．该方法首先在一粗网格上求解Navier-Stokes方程,然后在细网格的子区域上并行求解粗网格解的残差方程,以校正粗网格解．该方法实现简单,通信需求少．使用有限元局部误差估计,推导了并行方法所得近似解的误差界,同时通过数值算例,验证了其高效性．     

定常Navier-Stokes方程的一种高效稳定有限元方法

提出了二维定常Navier-Stokes(N-S)方程的一种两层稳定有限元方法.该方法基于局部高斯积分技术,通过不满足inf-sup条件的低次等阶有限元对N-S方程进行有限元求解.该方法在粗网格上解定常N-S方程,在细网格上只需解一个Stokes方程.误差分析和数值试验都表明:两层稳定有限元方法与直接在细网格上采用的传统有限元方法得到的解具有同阶的收敛性,但两层稳定有限元方法节省了大量的工作时间.     

非定常对流扩散方程的局部和并行有限元算法

对基于两重网格的非定常对流扩散方程的局部和并行有限元算法进行了研究.算法的理论依据是两重网格的思想,解的低频分量可以用一个整体的粗网格空间来逼近,高频分量可以用局部和并行的细网格空间来逼近.因此,这种局部和并行算法仅仅涉及一个粗网格上的整体逼近和细网格上的局部校正.得到了算法的误差估计,一些数值例子验证了算法的有效性.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

复杂边界旋转Navier-Stokes方程的几何方法及二度并行算法

本文提出了一种求解复杂边界旋转Navier-Stokes方程的微分几何方法及其二度并行算法.此方法可用于求解透平机械内部叶片间流动和飞行器外部绕流等复杂流动问题.假设流动区域可以用一系列光滑曲面■_k,k=1,2,…,K分割为一系列子区域(称作流层),通过应用微分几何的方法,三维N-S算子可以分解为两类算子之和:建立在曲面■_k切空间上"膜算子"和曲面■_k法线方向的"挠曲算子",将挠曲算子应用欧拉中心差商来逼近,由此得到建立在■_k上的"2D-3C"N-S方程.求解2D-3C N-S方程并且反复迭代直到收敛.我们得到"二度并行算法",它是2D-3C N-S方程并行算法与k方向的同时并行.这个算法的优点在于,(1)可以改进由于复杂边界造成的不规则三维网格引起的逼近解的精度;(2)为克服边界层的数值效应,在边界层内可以构造很密的流层,形成三维多尺度的网格,是一个很好的边界层算法;(3)这个方法不同于经典的区域分解算法,这里的每个子区域只需要求解一个"2D-3C"N-S方程,而经典区域分解方法要在每个子区域上求解三维问题.     

定常Navier-Stokes方程基于完全重叠型区域分解的并行稳定化有限元方法

基于完全重叠型区域分解技巧,针对低阶P_1-P_1有限元,本文提出求解二维定常不可压缩Navier-Stokes方程的并行稳定化有限元方法,其稳定项是基于两局部Gauss积分的压力投影.该方法的基本思想是,使用一局部加密的多尺度网格计算给定子区域上的局部稳定化有限元解.理论分析上借助有限元解的局部先验误差估计,推导出并行稳定化方法所得速度和压力解的误差界.选取适当的算法参数比例,该方法能取得与标准稳定化有限元方法相同的收敛阶,同时减少大量的计算时间.最后给出两类数值算例验证并行稳定化方法的高效性.     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

二阶椭圆问题带单位分解技巧的两重网格方法

标准的两重网格方法是一种求解二阶椭圆问题的局部并行方法,其计算所得数值解在整个求解区域上并不连续使用单位分解技术,将各个子区域上的局部解粘合在一起,从而得到全局连续解,并证明此解在H1范数意义下最优.更进一步,可以证明通过在粗网格上修正,能够改善其L2误差.数值例子验证了理论的正确性.     

粘性Cahn-Hilliard方程的时间双层网格混合有限元方法

主要针对在求解粘性Cahn-Hilliard方程时非线性项引起的时间耗时问题,提出了时间双层网格混合有限元方法.在空间上采用混合有限元方法进行离散,时间上采用Crank-Nicolson格式.首先在时间粗网格上,通过非线性牛顿迭代方法求解非线性混合有限元系统.其次基于初始迭代数值解和拉格朗日插值公式在时间细网格上求解线性混合有限元系统,然后证明了该方法的稳定性和误差估计,并通过数值算例对理论部分进行验证.结果表明,理论与数值算例相一致.     

求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

A parallel finite element variational multiscale method based on fully overlapping domain decomposition for incompressible flows

Based on fully overlapping domain decomposition and a recent variational multiscale method, a parallel finite element variational multiscale method for convection dominated incompressible flows is proposed and analyzed. In this method, each processor computes a local finite element solution in its own subdomain using a global mesh that is locally refined around its own subdomain, where a stabilization term based on two local Gauss integrations is adopted to stabilize the numerical form of the Navier–Stokes equations. Using the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 856–875, 2015     

A second-order finite element variational multiscale scheme for the fully discrete unsteady Navier–Stokes equations

In this report, we present and study a fully discrete finite element variational multiscale scheme for the unsteady incompressible Navier–Stokes equations where high Reynolds numbers are allowed. The scheme uses conforming finite element pairs for spatial discretization and a three-point difference formula for temporal discretization which is of second-order, where a stabilization term based on two local Gauss integrations is employed to stabilize the numerical scheme. We prove stability of the scheme, derive a priori error estimates for the fully discrete solution, and finally, give some numerical results to verify the theoretical predictions and demonstrate the effectiveness of the proposed numerical scheme.     

A parallel two-level finite element variational multiscale method for the Navier–Stokes equations

A combination method of two-grid discretization approach with a recent finite element variational multiscale algorithm for simulation of the incompressible Navier–Stokes equations is proposed and analyzed. The method consists of a global small-scale nonlinear Navier–Stokes problem on a coarse grid and local linearized residual problems in overlapped fine grid subdomains, where the numerical form of the Navier–Stokes equations on the coarse grid is stabilized by a stabilization term based on two local Gauss integrations at element level and defined by the difference between a consistent and an under-integrated matrix involving the gradient of velocity. By the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method.     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

Local and parallel finite element algorithms based on two-grid discretizations for the transient Stokes equations

Based on two-grid discretizations, some local and parallel finite element algorithms for the d-dimensional (d = 2,3) transient Stokes equations are proposed and analyzed. Both semi- and fully discrete schemes are considered. With backward Euler scheme for the temporal discretization, the basic idea of the fully discrete finite element algorithms is to approximate the generalized Stokes equations using a coarse grid on the entire domain, then correct the resulted residue using a finer grid on overlapped subdomains by some local and parallel procedures at each time step. By the technical tool of local a priori estimate for the fully discrete finite element solution, errors of the corresponding solutions from these algorithms are estimated. Some numerical results are also given which show that the algorithms are highly efficient.     

加罚N-S方程的有限元非线性Galerkin方法

非线性Galerkin方法是对耗散型非线性发展方程的一种数值解法,其空间变量不象一般Galerkin方法那样在线性空间上离散,而是在非线性流形上离散,所得逼近解在时间变量增大时可以更快地逼近其精确解.精细的理论分析可见[1],[2]等,在有限元逼近基础上将此方法应用到Navier-Stokes方程上的工作可参见[3],[4],这些文章主要针对速度与压力同时求解的混合元情形做了讨论.本文在[4]的基础上对加罚Navier-Stokes方程的一种非线性Galerkin方法的半离散和全离散有限元逼近格式分别进行了误差估     

Two-Grid Algorithm of $H^{1}$-Galerkin Mixed Finite Element Methods for Semilinear Parabolic Integro-Differential Equations

In this paper, we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by $H^{1}$-Galerkin mixed finite element methods. We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization, and backward Euler scheme for temporal discretization. Firstly, a priori error estimates and some superclose properties are derived. Secondly, a two-grid scheme is presented and its convergence is discussed. In the proposed two-grid scheme, the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy. Finally, a numerical experiment is implemented to verify theoretical results of the proposed scheme. The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice $h=H^2$.     

A finite element modified method of characteristics for convective heat transport

We propose a finite element modified method of characteristics for numerical solution of convective heat transport. The flow equations are the incompressible Navier‐Stokes equations including density variation through the Boussinesq approximation. The solution procedure consists of combining an essentially non‐oscillatory modified method of characteristics for time discretization with finite element method for space discretization. These numerical techniques associate the geometrical flexibility of the finite elements with the ability offered by modified method of characteristics to solve convection‐dominated flows using time steps larger than its Eulerian counterparts. Numerical results are shown for natural convection in a squared cavity and heat transport in the strait of Gibraltar. Performance and accuracy of the method are compared to other published data. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008