Two-Scale Picard Stabilized Finite Volume Method for the Incompressible Flow

In this paper, we consider a two-scale stabilized finite volume method for the two-dimensional stationary incompressible flow approximated by the lowest equal-order element pair $P_1-P_1$ which does not satisfy the inf-sup condition. The two-scale method consists of solving a small non-linear system on the coarse mesh and then solving a linear Stokes equations on the fine mesh. Convergence of the optimal order in the $H^1$-norm for velocity and the $L^2$-norm for pressure is obtained. The error analysis shows there is the same convergence rate between the two-scale stabilized finite volume solution and the usual stabilized finite volume solution on a fine mesh with relation $h =\mathcal{O}(H^2)$. Numerical experiments completely confirm theoretic results. Therefore, this method presented in this paper is of practical importance in scientific computation.     

A Stabilized Finite Element Method for Non-Stationary Conduction-Convection Problems

This paper is concerned with a stabilized finite element method based on two local Gauss integrations for the two-dimensional non-stationary conduction-convection equations by using the lowest equal-order pairs of finite elements. This method only offsets the discrete pressure space by the residual of the simple and symmetry term at element level in order to circumvent the inf-sup condition. The stability of the discrete scheme is derived under some regularity assumptions. Optimal error estimates are obtained by applying the standard Galerkin techniques. Finally, the numerical illustrations agree completely with the theoretical expectations.     

离散元与有限元结合的多尺度方法及其应用

在深入研究复杂结构和非均质材料冲击响应和破坏机理的过程中,往往遇到多尺度计算问题.提出并建立起离散元与有限元结合的多尺度方法,该方法采用离散元对感兴趣的局部进行细观尺度的模拟,利用有限元进行宏观的模拟,从而节约了计算时间.采用一种特殊的过渡层衔接离散元区和有限元区.将这一方法应用于激光辐照下预应力铝板的破坏响应,并将得到的模拟结果与实验进行了比较.     

Two-Level Finite Element Iterative Algorithm Based on Stabilized Method for the Stationary Incompressible Magnetohydrodynamics

In this paper, based on the stabilization technique, the Oseen iterative method and the two-level finite element algorithm are combined to numerically solve the stationary incompressible magnetohydrodynamic (MHD) equations. For the low regularity of the magnetic field, when dealing with the magnetic field sub-problem, the Lagrange multiplier technique is used. The stabilized method is applied to approximate the flow field sub-problem to circumvent the inf-sup condition restrictions. One- and two-level stabilized finite element algorithms are presented, and their stability and convergence analysis is given. The two-level method uses the Oseen iteration to solve the nonlinear MHD equations on a coarse grid of size H, and then employs the linearized correction on a fine grid with grid size h. The error analysis shows that when the grid sizes satisfy h=O(H2), the two-level stabilization method has the same convergence order as the one-level one. However, the former saves more computational cost than the latter one. Finally, through some numerical experiments, it has been verified that our proposed method is effective. The two-level stabilized method takes less than half the time of the one-level one when using the second class Nédélec element to approximate magnetic field, and even takes almost a third of the computing time of the one-level one when adopting the first class Nédélec element.     

Finite Element Iterative Methods for the 3D Steady Navier–Stokes Equations

In this work, a finite element (FE) method is discussed for the 3D steady Navier–Stokes equations by using the finite element pair Xh×Mh. The method consists of transmitting the finite element solution (uh,ph) of the 3D steady Navier–Stokes equations into the finite element solution pairs (uhn,phn) based on the finite element space pair Xh×Mh of the 3D steady linearized Navier–Stokes equations by using the Stokes, Newton and Oseen iterative methods, where the finite element space pair Xh×Mh satisfies the discrete inf-sup condition in a 3D domain Ω. Here, we present the weak formulations of the FE method for solving the 3D steady Stokes, Newton and Oseen iterative equations, provide the existence and uniqueness of the FE solution (uhn,phn) of the 3D steady Stokes, Newton and Oseen iterative equations, and deduce the convergence with respect to (σ,h) of the FE solution (uhn,phn) to the exact solution (u,p) of the 3D steady Navier–Stokes equations in the H1−L2 norm. Finally, we also give the convergence order with respect to (σ,h) of the FE velocity uhn to the exact velocity u of the 3D steady Navier–Stokes equations in the L2 norm.     

An Error Analysis for the Finite Element Approximation to the Steady-State Poisson-Nernst-Planck Equations

Poisson-Nernst-Planck equations are a coupled system of nonlinear partial differential equations consisting of the Nernst-Planck equation and the electrostatic Poisson equation with delta distribution sources, which describe the electrodiffusion of ions in a solvated biomolecular system. In this paper, some error bounds for a piecewise finite element approximation to this problem are derived. Several numerical examples including biomolecular problems are shown to support our analysis.     

A Two-Level Method for Pressure Projection Stabilized P1 Nonconforming Approximation of the Semi-Linear Elliptic Equations

In this paper, we study a new stabilized method based on the local pressure projection to solve the semi-linear elliptic equation. The proposed scheme combines nonconforming finite element pairs NCP₁−P₁triangle element and two-level method, which has a number of attractive computational properties: parameter-free, avoiding higher-order derivatives or edge-based data structures, but have more favorable stability and less support sets. Stability analysis and error estimates have been done. Finally, numerical experiments to check estimates are presented.     

一种二阶混合有限体元格式的GAMG预条件子

针对一种含跳系数椭圆问题的二阶混合有限体元格式,讨论求解相应离散系统PGMRES法的预条件子构造问题.通过严格的理论分析,建立分层基下该二阶混合有限体元刚度矩阵和二次有限元刚度矩阵的谱等价关系,并利用关于二次有限元刚度矩阵的一种基于分层思想的GAMG预条件子,为二阶混合有限体元刚度矩阵设计一种高效GAMG预条件子.数值结果验证理论分析的正确性和新预条件子的高效性与稳定性.     

间断有限元方法求解一维非平衡辐射扩散方程

研究一维非平衡辐射扩散方程的数值方法.通过求解间断系数热传导方程的广义黎曼问题,得到一种带加权数值流量,基于该数值流量构造了一类新型的间断有限元方法.在时间离散上采用向后Euler方法,形成的非线性方程组采用Picard迭代求解.数值试验表明该方法具有捕捉大梯度的能力,而且能适应扩散系数间断的情形.     

耦合FE/WB法在声分析中的应用

简要描叙FE法(finite element method)和WB法(wave based method)的理论背景以及耦合FE/WB法的数学基础.耦合FE/WB法利用两者的优势——FE法的广泛应用和WB法的高收敛特性,将FE模型中较大且几何简单的部分采用WB法代替.耦合模型具有相对较少的自由度.对于较高的频率还可以进行细分得到更高的计算精度,并利用模态缩减法进一步减少自由度数.数值算例结果表明,该耦合方法有潜力覆盖中频段的声分析.     

Two-Level Stabilized Finite Volume Methods for the Stationary Navier-Stokes Equations

In this work, two-level stabilized finite volume formulations for the 2D steady Navier-Stokes equations are considered. These methods are based on the local Gauss integration technique and the lowest equal-order finite element pair. Moreover, the two-level stabilized finite volume methods involve solving one small Navier-Stokes problem on a coarse mesh with mesh size $H$, a large general Stokes problem for the Simple and Oseen two-level stabilized finite volume methods on the fine mesh with mesh size $h$=$\mathcal{O}(H^2)$ or a large general Stokes equations for the Newton two-level stabilized finite volume method on a fine mesh with mesh size $h$=$\mathcal{O}(|\log h|^{1/2}H^3)$. These methods we studied provide an approximate solution $(\widetilde{u}_h^v,\widetilde{p}_h^v)$ with the convergence rate of same order as the standard stabilized finite volume method, which involve solving one large nonlinear problem on a fine mesh with mesh size $h$. Hence, our methods can save a large amount of computational time.     

Lagrange坐标系下二维气动方程组的RKDG有限元方法

构造Lagrange坐标系下二维可压缩气动方程组的RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法.将流体力学方程组和几何守恒律统-求解,所有计算都在固定的网格上进行,计算过程中不需要网格节点的速度信息.对几个数值算例进行数值模拟,得到较好的数值模拟结果.     

Numerical Analysis and Comparison of Four Stabilized Finite Element Methods for the Steady Micropolar Equations

In this paper, four stabilized methods based on the lowest equal-order finite element pair for the steady micropolar Navier–Stokes equations (MNSE) are presented, which are penalty, regular, multiscale enrichment, and local Gauss integration methods. A priori properties, existence, uniqueness, stability, and error estimation based on Fem approximation of all the methods are proven for the physical variables. Finally, some numerical examples are displayed to show the numerical characteristics of these methods.     

最小二乘有限元法求解非定常应力的Navier-Stokes方程

为精确求解非定常层流问题,发展一种非定常速度-应力-压力的方法.采用牛顿法对非线性对流项进行线性化处理和预处理共轭梯度法,实现了非定常应力形式Navier-Stokes方程的求解.方腔层流流动比较发现,非定常应力形式比涡量形式与试验结果更加吻合,精度更高.该方法有效地解决亚格子应力项的问题,实现基于最小二乘有限元法的湍流求解.比较方腔湍流流动的试验与仿真结果,证明本文的方法具有可行性,为湍流大涡模拟计算打下基础.     

A New Finite Volume Element Formulation for the Non-Stationary Navier-Stokes Equations

A semi-discrete scheme about time for the non-stationary Navier-Stokes equations is presented firstly, then a new fully discrete finite volume element (FVE) formulation based on macroelement is directly established from the semi-discrete scheme about time. And the error estimates for the fully discrete FVE solutions are derived by means of the technique of the standard finite element method. It is shown by numerical experiments that the numerical results are consistent with theoretical conclusions. Moreover, it is shown that the FVE method is feasible and efficient for finding the numerical solutions of the non-stationary Navier-Stokes equations and it is one of the most effective numerical methods among the FVE formulation, the finite element formulation, and the finite difference scheme.     

A Finite Volume Method Based on the Constrained Nonconforming Rotated Q1-Constant Element for the Stokes Problem

We construct a finite volume element method based on the constrained nonconforming rotated Q₁-constant element (CNRQ₁-P₀) for the Stokes problem. Two meshes are needed, which are the primal mesh and the dual mesh. We approximate the velocity by CNRQ₁ elements and the pressure by piecewise constants. The errors for the velocity in the H¹ norm and for the pressure in the L² norm are O(h) and the error for the velocity in the L² norm is O(h²). Numerical experiments are presented to support our theoretical results.     

二维多介质可压缩流的RKDG有限元方法

应用RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法、Level Set方法和Ghost Fluid方法数值模拟二维多介质可压缩流,其中Euler方程组、Level Set方程和重新初始化方程的空间离散采用DG(Discontinuous Galerkin)有限元方法,时间离散采用Runge-Kutta方法.对二维的气-气和气-液两相流进行了数值计算,得到了分辨率较高的计算结果.     

Finite Element Analysis of Maxwell's Equations in Dispersive Lossy Bi-Isotropic Media

In this paper, the time-dependent Maxwell's equations used to modeling wave propagation in dispersive lossy bi-isotropic media are investigated. Existence and uniqueness of the modeling equations are proved. Two fully discrete finite element schemes are proposed, and their practical implementation and stability are discussed.     

RKDG有限元法求解一维拉格朗日形式的Euler方程

描述一种新的求解Euler方程的拉格朗日格式,该格式用Runge-Kutta Discontinuous Galerkin(RKDG)方法在拉格朗日坐标系求解Euler方程,剖分网格随流体运动.新格式不仅保证流体的质量、动量和能量守恒,而且能够在时间和空间上同时达到二阶精度.数值算例表明在一维情况,随着拉氏网格的移动和改变,格式在时间和空间上仍保持二阶精度,并且没有数值震荡.     

An Iterative Two-Grid Method of a Finite Element PML Approximation for the Two Dimensional Maxwell Problem

In this paper, we propose an iterative two-grid method for the edge finite element discretizations (a saddle-point system) of Perfectly Matched Layer (PML) equations to the Maxwell scattering problem in two dimensions. Firstly, we use a fine space to solve a discrete saddle-point system of $H(grad)$ variational problems, denoted by auxiliary system 1. Secondly, we use a coarse space to solve the original saddle-point system. Then, we use a fine space again to solve a discrete$\boldsymbol{H}(curl)$-elliptic variational problems, denoted by auxiliary system 2. Furthermore, we develop a regularization diagonal block preconditioner for auxiliary system 1 and use $H$-$X$ preconditioner for auxiliary system 2. Hence we essentially transform the original problem in a fine space to a corresponding (but much smaller) problem on a coarse space, due to the fact that the above two preconditioners are efficient and stable. Compared with some existing iterative methods for solving saddle-point systems, such as PMinres, numerical experiments show the competitive performance of our iterative two-grid method.