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

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

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

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

求解复系数椭圆方程的两层网格算法

<正>1引言两层网格方法是用来求解非对称不定问题和非线性问题的一种非常有效的数值方法[1,2].其主要思想是,借助于两层网格空间,将细网格上的复杂问题转化为求解一个细网格空间的简单问题和一个粗网格上的问题.由于粗网格空间相对于细网格空间很小,所以减少了计算代价,并且仍能得到原问题的最优解.因此,两层网格算法被广泛研究并被用于求解多种问题,例如,求解非对称和非线性椭圆方程[1,2,3,4],非线性弹性方程[5],Navier-Stokes方程[6,7,8]及特征值问题[9,10].HSS迭代方法是求解大规模稀疏非埃尔米特正定     

Navier-Stokes方程流函数形式两重网格算法的误差分析

对定常Navier-Stokes方程流函数形式两重网格有限元算法进行了误差分析。此方法包括在粗网格上求解一个非线性问题，在细网格上求解一个线性问题，然后再在粗网格上求解一个线性校正问题。分析了包括校正项和不包括校正项两种方法的误差，得出对于任意固定的Beynolds数，能达到最优逼近阶。     

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

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

大Reynolds数Navier-Stokes方程带回溯技巧的两水平方法

本文考虑了一种求解大Reynolds数定常Navier-Stokes方程带回溯(backtracking)技巧的两水平有限元方法.其基本思想是,首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上求解一个亚格子模型稳定化的线性Newton问题,最后又回到粗网格上求解线性化的校正问题.通过适当的稳定化参数和粗细网格尺寸的选取,本文的算法能取得最优渐近收敛阶.数值实验检验了理论分析的正确性和算法的有效性.     

非线性整规划中的精确光滑罚函数

本文提出了几个非线性整规划 的全局精确光滑罚函数，每个罚函数有两个参数，并且给出了每个罚函数的精确罚参数的估计值，最后，我们举例说明了所提出的罚方法在具有整系数多项式目标函数以约束函数的整数规划中的应用。     

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

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

带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

基于Picard迭代的PN×PN-2谱元法求解定常不可压缩Navier-Stokes方程

该文给出了一种求解二维定常不可压缩Navier-Stokes方程的基于Picard线性化迭代的P_N×P_N-2谱元法.通过Picard线性化将不可压缩Navier-Stokes方程的求解转化为一系列线性的Stokes-type方程,再利用非交错网格的P_N×P_N-2谱元法计算每个迭代步的Stokes-type方程.为了消除伪压力模,压力离散比速度离散低两阶,非交错网格的应用使得方程的离散方便且不会带来相应的插值误差,从而保证了谱精度.通过此方法数值计算了有精确解的Stokes流动、Kovasznay流动和方腔顶盖驱动流,结果表明,迭代收敛非常快,误差收敛达到了谱精度收敛,并且避免了压力震荡的出现,表明了该文方法准确可靠.     

A Two‐Parameter Stabilized Finite Element Method for Incompressible Flows

Based on two‐grid discretizations, a two‐parameter stabilized finite element method for the steady incompressible Navier–Stokes equations at high Reynolds numbers is presented and studied. In this method, a stabilized Navier–Stokes problem is first solved on a coarse grid, and then a correction is calculated on a fine grid by solving a stabilized linear problem. The stabilization term for the nonlinear Navier–Stokes equations on the coarse grid is based on an elliptic projection, which projects higher‐order finite element interpolants of the velocity into a lower‐order finite element interpolation space. For the linear problem on the fine grid, either the same stabilization approach (with a different stabilization parameter) as that for the coarse grid problem or a completely different stabilization approach could be employed. Error bounds for the discrete solutions are estimated. Algorithmic parameter scalings of the method are also derived. The theoretical results show that, with suitable scalings of the algorithmic parameters, this method can yield an optimal convergence rate. Numerical results are provided to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 425–444, 2017     

Two‐level Newton iterative method for the 2D/3D steady Navier‐Stokes equations

A combination method of the Newton iteration and two‐level finite element algorithm is applied for solving numerically the steady Navier‐Stokes equations under the strong uniqueness condition. This algorithm is motivated by applying the m Newton iterations for solving the Navier‐Stokes problem on a coarse grid and computing the Stokes problem on a fine grid. Then, the uniform stability and convergence with respect to ν of the two‐level Newton iterative solution are analyzed for the large m and small H and h << H. Finally, some numerical tests are made to demonstrate the effectiveness of the method. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Error estimates for two‐level penalty finite volume method for the stationary Navier–Stokes equations

Two‐level penalty finite volume method for the stationary Navier–Stokes equations based on the P₁ ? P₀ element is considered in this paper. The method involves solving one small penalty Navier–Stokes problem on a coarse mesh with mesh size H = ?^{1 / 4}h^{1 / 2}, a large penalty Stokes problem on a fine mesh with mesh size h, where 0 < ? < 1 is a penalty parameter. The method we study provides an approximate solution with the convergence rate of same order as the penalty finite volume solution (u_?h,p_?h), which involves solving one large penalty Navier–Stokes problem on a fine mesh with the same mesh size h. However, our method can save a large amount of computational time. Copyright © 2013 John Wiley & Sons, Ltd.     

Two-level iteration penalty and variational multiscale method for steady incompressible flows

In this paper, we study two-level iteration penalty and variational multiscale method for the approximation of steady Navier-Stokes equations at high Reynolds number. Comparing with classical penalty method, this new method does not require very small penalty parameter $\varepsilon$. Moreover, two-level mesh method can save a large amount of CPU time. The error estimates in $H^1$ norm for velocity and in $L^2$ norm for pressure are derived. Finally, two numerical experiments are shown to support the efficiency of this new method.     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

Iterative methods in penalty finite element discretizations for the Steady Navier–Stokes equations

Three penalty finite element methods are designed to solve numerically the steady Navier–Stokes equations, where the Stokes, Newton, and Oseen iteration methods are used, respectively. Moreover, the stability analysis and error estimate for these nine algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms. Meanwhile, the numerical investigations are provided to show that the proposed methods are efficient for solving the steady Navier–Stokes equations with the different viscosity. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 74‐94, 2014     

A penalty‐FEM for navier‐stokes type variational inequality with nonlinear damping term

In this article, we consider a penalty finite element (FE) method for incompressible Navier‐Stokes type variational inequality with nonlinear damping term. First, we establish penalty variational formulation and prove the well‐posedness and convergence of this problem. Then we show the penalty FE scheme and derive some error estimates. Finally, we give some numerical results to verify the theoretical rate of convergence. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 918–940, 2017     

Navier-Stokes方程带Backtracking技巧的两重网格算法

1 引 言考虑二维不可压 Navier-Stokes方程:     

Alternating Direction Method with Self-Adaptive Penalty Parameters for Monotone Variational Inequalities

The alternating direction method is one of the attractive approaches for solving linearly constrained separate monotone variational inequalities. Experience on applications has shown that the number of iterations depends significantly on the penalty parameter for the system of linear constraint equations. While the penalty parameter is a constant in the original method, in this paper we present a modified alternating direction method that adjusts the penalty parameter per iteration based on the iterate message. Preliminary numerical tests show that the self-adaptive adjustment technique is effective in practice.     

A three-step Oseen-linearized finite element method for incompressible flows with damping

This paper introduces a three-step Oseen-linearized finite element method for the 2D/3D steady incompressible Navier–Stokes equations with nonlinear damping term. Within this method, we solve a nonlinear problem over a coarse grid followed by solving two Oseen-linearized problems over a fine grid, which possess the same stiffness matrices with only various right-hand sides. We theoretically analyze the stability of the present method, and derive optimal error estimates of the finite element solutions. We conduct a series of numerical experiments which support the theoretical analysis and test the effectiveness of the proposed method. We demonstrate numerically that there is a significant improvement in the accuracy of the approximate solutions over those for the standard two-level method.