基于对偶混合变分原理的Signorini问题的数值模拟

基于Signorini问题的对偶混合变分形式,提出了一种非协调有限元逼近格式,证明了离散的B-B条件,获得了Raviart-Thomas(k=0)有限元逼近的误差界O(h^3/4),并且Uzawa型算法对协调与非协调有限元逼近格式进行了数值求解.根据数值结果的分析和比较,表明应用非协调有限元逼近格式求解更有效.     

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.     

Optimal Convergence Analysis of Two-Level Nonconforming Finite Element Iterative Methods for 2D/3D MHD Equations

Several two-level iterative methods based on nonconforming finite element methods are applied for solving numerically the 2D/3D stationary incompressible MHD equations under different uniqueness conditions. These two-level algorithms are motivated by applying the m iterations on a coarse grid and correction once on a fine grid. A one-level Oseen iterative method on a fine mesh is further studied under a weak uniqueness condition. Moreover, the stability and error estimate are rigorously carried out, which prove that the proposed methods are stable and effective. Finally, some numerical examples corroborate the effectiveness of our theoretical analysis and the proposed methods.     

Coefficient Jump-Independent Approximation of the Conforming and Nonconforming Finite Element Solutions

A counterexample is constructed. It confirms that the error of conforming finite element solution is proportional to the coefficient jump, when solving interface elliptic equations. The Scott-Zhang operator is applied to a nonconforming finite element. It is shown that the nonconforming finite element provides the optimal order approximation in interpolation, in $L^2$-projection, and in solving elliptic differential equation, independent of the coefficient jump in the elliptic differential equation. Numerical tests confirm the theoretical finding.     

Finite element method for modeling radiative transfer in semitransparent graded index cylindrical medium

Both Galerkin finite element method (GFEM) and least squares finite element method (LSFEM) are developed and their performances are compared for solving the radiative transfer equation of graded index medium in cylindrical coordinate system (RTEGC). The angular redistribution term of the RTEGC is discretized by finite difference approach and after angular discretization the RTEGC is formulated into a discrete-ordinates form, which is then discretized based on Galerkin or least squares finite element approach. To overcome the RTEGC-led numerical singularity at the origin of cylindrical coordinate system, a pole condition is proposed as a special mathematical boundary condition. Compared with the GFEM, the LSFEM has very good numerical properties and can effectively mitigate the nonphysical oscillation appeared in the GFEM solutions. Various problems of both axisymmetry and nonaxisymmetry, and with medium of uniform refractive index distribution or graded refractive index distribution are tested. The results show that both the finite element approaches have good accuracy to predict the radiative heat transfer in semitransparent graded index cylindrical medium, while the LSFEM has better numerical stability.     

Lower Bounds for Eigenvalues of the Stokes Operator

In this paper, we propose a condition that can guarantee the lower bound property of the discrete eigenvalue produced by the finite element method for the Stokes operator. We check and prove this condition for four nonconforming methods and one conforming method. Hence they produce eigenvalues which are smaller than their exact counterparts.     

Navier-Stokes方程的最小二乘等几何方法

基于Hermite多项式的C¹型单元构造复杂,限制了最小二乘有限元法的应用.引入高阶光滑的非均匀有理B样条作为基函数简化C¹型单元构造,提出求解黏性不可压流动Navier-Stokes方程的最小二乘等几何方法.用Newton法或Picard法对Navier-Stokes方程线性化,用线性化偏微分方程的余量定义最小二乘泛函,导出最小二乘变分方程,用NURBS构造高阶光滑的有限维空间来近似速度场和压力场.计算表明:本文方法计算的二维顶盖驱动流数值解能准确描述流动状况,计算的二维通道内圆柱绕流全局质量损失由最小二乘有限元法的6%降为0.018%,该方法可用于Navier-Stokes方程的求解,并且具有较好的质量守恒性.     

Least-squares spectral element solution of incompressible Navier–Stokes equations with adaptive refinement

Least-squares spectral element solution of steady, two-dimensional, incompressible flows are obtained by approximating velocity, pressure and vorticity variable set on Gauss–Lobatto–Legendre nodes. Constrained Approximation Method is used for h- and p-type nonconforming interfaces of quadrilateral elements. Adaptive solutions are obtained using a posteriori error estimates based on least squares functional and spectral coefficient. Effective use of p-refinement to overcome poor mass conservation drawback of least-squares formulation and successful use of h- and p-refinement together to solve problems with geometric singularities are demonstrated. Capabilities and limitations of the developed code are presented using Kovasznay flow, flow past a circular cylinder in a channel and backward facing step flow.     

Finite Element $θ$-Schemes for the Acoustic Wave Equation

In this paper, we investigate the stability and convergence of a family of implicit finite difference schemes in time and Galerkin finite element methods in space for the numerical solution of the acoustic wave equation. The schemes cover the classical explicit second-order leapfrog scheme and the fourth-order accurate scheme in time obtained by the modified equation method. We derive general stability conditions for the family of implicit schemes covering some well-known CFL conditions. Optimal error estimates are obtained. For sufficiently smooth solutions, we demonstrate that the maximal error in the $L^2$-norm error over a finite time interval converges optimally as $\mathcal{O}(h^{p+1}+∆t^s)$, where $p$ denotes the polynomial degree, $s$=2 or 4, $h$ the mesh size, and $∆t$ the time step.     

Analysis of the Three-Dimensional Model of Diffusion of Minority Charge Carriers Generated by an Electron Probe in a Heterogeneous Semiconductor Material by Means of Projection Methods

Algorithms for using the Galerkin projection method and the projection least squares method to analyze the three-dimensional model of the diffusion of minority charge carriers generated by an electron probe in a semiconductor material are presented. The results obtained using these methods are compared with the analytical solution. An estimate of the error is given, and the condition for the computation stability of the projection least squares method in the form of the limiting relation is obtained.     

High-order explicit-implicit numerical methods for nonlinear anomalous diffusion equations

In this paper, the high-order finite difference/element methods for the nonlinear anomalous diffusion equations of subdiffusion and superdiffusion are developed, where the high-order finite difference methods are used to approximate the time-fractional derivatives and the finite element methods are used in the spatial domain. The stability and error estimates are proved for both cases of superdiffusion and subdiffusion. Numerical examples are provided to confirm the theoretical analysis.     

A deficiency problem of the least squares finite element method for solving radiative transfer in strongly inhomogeneous media

The accuracy and stability of the least squares finite element method (LSFEM) and the Galerkin finite element method (GFEM) for solving radiative transfer in homogeneous and inhomogeneous media are studied theoretically via a frequency domain technique. The theoretical result confirms the traditional understanding of the superior stability of the LSFEM as compared to the GFEM. However, it is demonstrated numerically and proved theoretically that the LSFEM will suffer a deficiency problem for solving radiative transfer in media with strong inhomogeneity. This deficiency problem of the LSFEM will cause a severe accuracy degradation, which compromises the performance of the LSFEM too much and makes it not a good choice to solve radiative transfer in strongly inhomogeneous media. It is also theoretically proved that the LSFEM using the one dimensional linear element is equivalent to a second order form of radiative transfer equation discretized by the central difference scheme.     

Numerical Analysis and Comparison of Three Iterative Methods Based on Finite Element for the 2D/3D Stationary Micropolar Fluid Equations

In this paper, three iterative methods (Stokes, Newton and Oseen iterative methods) based on finite element discretization for the stationary micropolar fluid equations are proposed, analyzed and compared. The stability and error estimation for the Stokes and Newton iterative methods are obtained under the strong uniqueness conditions. In addition, the stability and error estimation for the Oseen iterative method are derived under the uniqueness condition of the weak solution. Finally, numerical examples test the applicability and the effectiveness of the three iterative methods.     

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.     

A Posteriori Error Estimates of Triangular Mixed Finite Element Methods for Semilinear Optimal Control Problems

In this paper, we present an a posteriori error estimates of semilinear quadratic constrained optimal control problems using triangular mixed finite element methods. The state and co-state are approximated by the order $k\leq 1$ Raviart- Thomas mixed finite element spaces and the control is approximated by piecewise constant element. We derive a posteriori error estimates for the coupled state and control approximations. A numerical example is presented in confirmation of the theory.     

自适应有限元方法和后验误差估计

自适应有限元方法是科学研究和工程设计领域中非常有效的一种求解偏微分方程的数值计算方法。这种方法是为了以尽可能小的代价取得尽可能好的计算效果。后验误差估计是实现自适应有限元计算的关键性手段。文章综合介绍了自适应有限元方法和后验误差估计在求解椭圆型方程、抛物型方程和双曲型方程方面所取得的比较新的成就。     

Convergence Analysis of Carey Nonconforming Finite Element for the Second-Order Elliptic Problem with the Lowest Regularity

Russian Physics Journal - In this paper, the Carey nonconforming finite element method (NFEM) for the second order elliptic problem is discussed. By means of the different techniques from the...     

两种偏最小二乘分光光度法同时测定三组分混合物

本文研究两种偏最小二乘法（经典编最小二乘法（ＣＰＬＳ）和基于核心矩阵的偏最小二乘法（ＫＰＬＳ）同时测定三组分混合物，根据数学原理编制三个程序（ＳＰＧＲＡＦＡ，ＳＰＧＲＰＬＳ和ＳＰＧＲＫＰＬＳ）执行这些计算，八个误差函数用以推断因子数目，因为核心矩阵维数小于原始数据矩阵，所以ＫＰＬＳ法适于计算具有较多光谱数和较少样品数的数据矩阵，实验结果显示对相互重叠的光谱用这两种方法均能获得令人满意且十分吻合的结     

Error Estimates and Superconvergence of Mixed Finite Element Methods for Optimal Control Problems with Low Regularity

In this paper, we investigate the error estimates and superconvergence property of mixed finite element methods for elliptic optimal control problems. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. We derive $L^2$ and $L^\infty$-error estimates for the control variable. Moreover, using a recovery operator, we also derive some superconvergence results for the control variable. Finally, a numerical example is given to demonstrate the theoretical results.     

空间望远镜分块式主镜建模与面形控制方法

空间望远镜分块式主镜的面形是由其背后布置的若干致动器控制的,是一个复杂的控制系统。应用BP神经网络的方法建立了以致动器作用力作为输入、镜面形变的Zernike多项式拟合系数作为输出的镜面形变模型。利用镜面有限元分析的大量数据对该模型进行了离线训练,并在最小二乘法的基础上,设计了加入单纯形修正算法的主镜面形静态控制器。仿真结果表明,应用该控制器对空间望远镜进行在线控制,控制精度优于最小二乘控制法。