椭圆型偏微分方程的弱有限元方法 献给林群教授80华诞

弱有限元方法 (weak Galerkin finite element methods),简称WG方法,是求解偏微分方程的一种全新且高效的数值方法.WG有限元方法的主要思想是,对间断函数引入广义弱微分算子,并将其应用于通常的变分形式中以对相应的偏微分方程进行数值求解,数值解的连续性则通过稳定子以弱形式来实现.本文以二阶椭圆问题为例,详细介绍弱有限元方法的原理和基础,并给出相关的理论分析.此外,本文简单介绍其他椭圆问题的弱有限元方法的离散格式.弱有限元方法的最大特点是,(1)有限元剖分允许任意多边形或多面体;(2)总体刚度矩阵可通过单元刚度矩阵叠加而得;(3)逼近函数构造简单,且极易满足相应的稳定性条件;(4)格式可做杂交处理以并行消去单元内部自由度.     

《数值计算与计算机应用》2021年第42卷第4期摘要

潘佳佳,李会元,二阶椭圆问题的弱迦辽金四边形谱元方法[J].数值计算与计算机应用,2021,42(4):303-322.摘要:本文对二阶椭圆方程特征值问题的弱伽辽金谱元方法开展相关数值研究.与弱有限元方法类似,弱伽辽金谱元方法的逼近函数空间包括各个单元上的独立内部分量、并辅以各单元边界分量作为单元与单元间的联系.本文聚焦任意凸四边形网格剖分下的弱伽辽金四边形谱元方法,弱逼近函数中的各内部分量与边界分量分别由参考正方形单元与参考单元边界上的正交多项式通过双线性变换来构造;而弱梯度逼近空间则由参考正方形上的正交多项式通过Piola变换构造.在此基础上,本文提出了二阶椭圆方程特征值问题的弱伽辽金四边形谱元方法逼近格式和实现算法,并通过对离散弱梯度核空间的系统研究。     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

拟线性对流占优扩散方程的扩展特征混合有限元方法数值模拟

讨论了拟线性对流占优扩散问题的数值模拟.对对流部分采用特征线格式进行离散,以消除流动锋线前沿的数值弥散现象,保证格式的稳定性;而对扩散部分采用扩展混合有限元方法,同时逼近未知函数,未知函数的梯度及伴随向量函数.理论分析和数值算例表明,此方法是稳定的,具有最优L2逼近精度.     

关于线弹性问题混合格式稳定化思想的探讨

首先从混合有限元理论出发,探讨线弹性问题混合变分格式所满足的稳定性条件,从而保证解的存在唯一性.使用连续的分片线性函数和分片常数来分别逼近应力和位移,详细分析了混合格式下稳定化的必要性,有助于更加深入地了解稳定化的基本思想.然后,通过在混合格式中引入位移的跳跃惩罚项,展示了一个无闭锁稳定化混合有限元方法,并证明了此方法是稳定的且是线性收敛的.     

半线性抛物最优控制问题全离散插值系数有限元方法的收敛性分析

本文考虑全离散插值系数有限元方法求解半线性抛物最优控制问题,其中控制变量用分片常数函数逼近,状态变量和对偶状态变量用分片线性函数逼近.对于方程中的半线性项,先用插值系数技巧处理,再用牛顿迭代法求解.通过引入一些辅助变量和投影算子,并利用有限元空间的逼近性质,得到半线性抛物最优控制问题插值系数有限元方法的收敛性结果;数值算例结果验证了理论结果的正确性.     

弱有限元方法简论

本文简述弱有限元方法(weak Galerkin finite element met,hods)的数学基本原理和计算机实现.弱有限元方法对间断函数引入广义弱微分,并将其应用于偏微分方程相应的变分形式进行数值求解,而数值解的弱连续性则通过稳定子或光滑子来实现.弱有限元方法针对广义函数而构建,是经典有限元方法的一种自然拓广,且能够弥补经典有限元方法的某些缺憾,也因此在科学与工程计算领域具有广泛的应用前景.     

基于时间有限元方法的旋转柔性叶片动力学响应分析

将时间有限元方法引入到柔性多体系统的数值计算中,研究了旋转柔性叶片系统的刚-柔耦合响应问题.首先,基于非线性梁理论,建立了旋转柔性叶片系统的中心刚体-柔性梁模型,构造柔性叶片系统考虑一次近似耦合的Lagrange函数;其次,采用假设模态方法对空间坐标进行离散,建立系统的时间有限元格式;最后,通过数值实验,分析了柔性叶片的动力学响应.该方法直接构造了系统的离散积分格式,并自动保证了该格式是保辛的,因而具有较高的数值精度和稳定性.数值结果表明:时间有限元可以有效地求解旋转柔性叶片系统内低频大范围运动与高频弹性振动之间的刚-柔耦合问题.     

具有弱对称应力的线弹性方程的稳定化格式

利用稳定化方法讨论拉格朗日乘子法得到的具有弱对称应力的线弹性问题. 用线性元和分片常数分别逼近变分问题的应力和位移. 并通过添加稳定项$G_1(\cdot,\cdot)$, $G_2(\cdot,\cdot)$和$G_3(\cdot,\cdot)$ 使相应混合离散变分问题满足弱BB条件. 接着详细研究了变分问题的解与稳定混合有限元解之间的误差估计,最后用两个数值算例验证理论分析的有效性.     

双调和方程混合元的一种新格式

本文介绍了双调和方程混合元的一种新格式,用双二次多项式逼近流函数,双一次多项式逼近涡函数.在拟一致矩形剖分的条件下,证明了此格式具有与C-R格式中分别用双二次多项式逼近相同的收敛阶.     

随机平面线弹性问题的一类弱Galerkin方法

本文针对带有随机杨氏模量和荷载的平面线弹性问题,提出了一类随机弱Galerkin有限元方法.先利用Karhunen-Loève展开把随机项参数化,将方程转化为一个确定性问题;再采用弱Galerkin有限元法和k-/p-型方法分别离散空间区域和随机场.在弱Galerkin离散中,用分片s(s≥1)和s+1次多项式逼近单元...     

A discrete maximum principle for the weak Galerkin finite element method on nonuniform rectangular partitions

This article establishes a discrete maximum principle (DMP) for the approximate solution of convection–diffusion–reaction problems obtained from the weak Galerkin (WG) finite element method on nonuniform rectangular partitions. The DMP analysis is based on a simplified formulation of the WG involving only the approximating functions defined on the boundary of each element. The simplified weak Galerkin (SWG) method has a reduced computational complexity over the usual WG, and indeed provides a discretization scheme different from the WG when the reaction terms are present. An application of the SWG on uniform rectangular partitions yields some 5- and 7-point finite difference schemes for the second order elliptic equation. Numerical experiments are presented to verify the DMP and the accuracy of the scheme, particularly the finite difference scheme.     

A hybridized weak Galerkin finite element scheme for the Stokes equations

In this paper a hybridized weak Galerkin(HWG) finite element method for solving the Stokes equations in the primary velocity-pressure formulation is introduced.The WG method uses weak functions and their weak derivatives which are defined as distributions.Weak functions and weak derivatives can be approximated by piecewise polynomials with various degrees.Different combination of polynomial spaces leads to different WG finite element methods,which makes WG methods highly flexible and efficient in practical computation.A Lagrange multiplier is introduced to provide a numerical approximation for certain derivatives of the exact solution.With this new feature,the HWG method can be used to deal with jumps of the functions and their flux easily.Optimal order error estimates are established for the corresponding HWG finite element approximations for both primal variables and the Lagrange multiplier.A Schur complement formulation of the HWG method is derived for implementation purpose.The validity of the theoretical results is demonstrated in numerical tests.     

Prox-regularization and solution of ill-posed elliptic variational inequalities

In this paper new methods for solving elliptic variational inequalities with weakly coercive operators are considered. The use of the iterative prox-regularization coupled with a successive discretization of the variational inequality by means of a finite element method ensures well-posedness of the auxiliary problems and strong convergence of their approximate solutions to a solution of the original problem.In particular, regularization on the kernel of the differential operator and regularization with respect to a weak norm of the space are studied. These approaches are illustrated by two nonlinear problems in elasticity theory.     

A coupled finite element-differential quadrature element method and its accuracy for moving load problem

In this paper a mixed method, which combines the finite element method and the differential quadrature element method (DQEM), is presented for solving the time dependent problems. In this study, the finite element method is first used to discretize the spatial domain. The DQEM is then employed as a step-by-step DQM in time domain to solve the resulting initial value problem. The resulting algebraic equations can be solved by either direct or iterative methods. Two general formulations using the DQM are also presented for solving a system of linear second-order ordinary differential equations in time. The application of the formulation is then shown by solving a sample moving load problem. Numerical results show that the present mixed method is very efficient and reliable.     

A weak form quadrature element method for plane elasticity problems

A weak form quadrature element method is proposed and applied to analysis of plane elasticity problems. A variational formulation of plane elasticity problems is established and the differential quadrature analog of the derivatives in the functional is introduced. Several typical plane elasticity problems are studied to verify the proposed method. Results show that the method is highly efficient and promising. It is applied to the analysis of nearly incompressible materials and shown to be robust against volumetric locking. Similarities and dissimilarities, advantages and disadvantages as compared with other numerical methods, typically the p-version finite element method are discussed.     

A new type of full multigrid method for the elasticity eigenvalue problem

This paper introduces a new type of full multigrid method for the elasticity eigenvalue problem. The main idea is to avoid solving large scale elasticity eigenvalue problem directly by transforming the solution of the elasticity eigenvalue problem into a series of solutions of linear boundary value problems defined on a multilevel finite element space sequence and some small scale elasticity eigenvalue problems defined on the coarsest correction space. The involved linear boundary value problems will be solved by performing some multigrid iterations. Besides, some efficient techniques such as parallel computing and adaptive mesh refinement can also be absorbed in our algorithm. The efficiency and validity of the multigrid methods are verified by several numerical experiments.     

Numerical and theoretical study of weak Galerkin finite element solutions of Turing patterns in reaction–diffusion systems

In this paper, we introduce numerical schemes and their analysis based on weak Galerkin finite element framework for solving 2‐D reaction–diffusion systems. Weak Galerkin finite element method (WGFEM) for partial differential equations relies on the concept of weak functions and weak gradients, in which differential operators are approximated by weak forms through the Green's theorem. This method allows the use of totally discontinuous functions in the approximation space. In the current work, the WGFEM solves reaction–diffusion systems to find unknown concentrations (u, v) in element interiors and boundaries in the weak Galerkin finite element space WG(P₀, P₀, RT₀) . The WGFEM is used to approximate the spatial variables and the time discretization is made by the backward Euler method. For reaction–diffusion systems, stability analysis and error bounds for semi‐discrete and fully discrete schemes are proved. Accuracy and efficiency of the proposed method successfully tested on several numerical examples and obtained results satisfy the well‐known result that for small values of diffusion coefficient, the steady state solution converges to equilibrium point. Acquired numerical results asserted the efficiency of the proposed scheme.     

Hybridized weak Galerkin finite element method for linear elasticity problem in mixed form

A hybridization technique is applied to the weak Galerkin finite element method (WGFEM) for solving the linear elasticity problem in mixed form. An auxiliary function, the Lagrange multiplier defined on the boundary of elements, is introduced in this method. Consequently, the computational costs are much lower than the standard WGFEM. Optimal order error estimates are presented for the approximation scheme. Numerical results are provided to verify the theoretical results.