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.     

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.     

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.     

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

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

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.     

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.     

Moving Finite Element Simulations for Reaction-Diffusion Systems

This work is concerned with the numerical simulations for two reaction-diffusion systems, i.e., the Brusselator model and the Gray-Scott model. The numerical algorithm is based upon a moving finite element method which helps to resolve large solution gradients. High quality meshes are obtained for both the spot replication and the moving wave along boundaries by using proper monitor functions. Unlike [33], this work finds out the importance of the boundary grid redistribution which is particularly important for a class of problems for the Brusselator model. Several ways for verifying the quality of the numerical solutions are also proposed, which may be of important use for comparisons.     

Error Analysis and Adaptive Methods of Least Squares Nonconforming Finite Element for the Transport Equations

In this paper, we consider a least squares nonconforming finite element of low order for solving the transport equations. We give a detailed overview on the stability and the convergence properties of our considered methods in the stability norm. Moreover, we derive residual type a posteriori error estimates for the least squares nonconforming finite element methods under $H^{−1}$-norm, which can be used as the error indicators to guide the mesh refinement procedure in the adaptive finite element method. The theoretical results are supported by a series of numerical experiments.     

一类各项异性半线性椭圆方程自然边界元与有限元耦合法

将冯康和余德浩提出的自然边界归化方法用于研究一类半线性椭圆方程外区域问题,提出一种自然边界元与有限元的耦合算法、针对某一类半线性椭圆方程,应用变分原理,研究其弱解性及Galerkin逼近,得到有限元解的误差估计及收敛阶O(h^n),最后给出相应数值例子。     

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.     

Finite Difference/Element Method for a Two-Dimensional Modified Fractional Diffusion Equation

We present the finite difference/element method for a two-dimensional modified fractional diffusion equation. The analysis is carried out first for the time semi-discrete scheme, and then for the full discrete scheme. The time discretization is based on the $L1$-approximation for the fractional derivative terms and the second-order backward differentiation formula for the classical first order derivative term. We use finite element method for the spatial approximation in full discrete scheme. We show that both the semi-discrete and full discrete schemes are unconditionally stable and convergent. Moreover, the optimal convergence rate is obtained. Finally, some numerical examples are tested in the case of one and two space dimensions and the numerical results confirm our theoretical analysis.     

An Inf-Sup Stabilized Finite Element Method by Multiscale Functions for the Stokes Equations

In the paper, an inf-sup stabilized finite element method by multiscale functions for the Stokes equations is discussed. The key idea is to use a Petrov-Galerkin approach based on the enrichment of the standard polynomial space for the velocity component with multiscale functions. The inf-sup condition for $P_1$-$P_0$ triangular element (or $Q_1$-$P_0$ quadrilateral element) is established. And the optimal error estimates of the stabilized finite element method for the Stokes equations are obtained.     

Finite Element Analysis of the Dominant Mode Variation in Groove Guide

In this paper, the dominant mode in groove guide is analyzed by finite element method. For the guide with a shallow groove, the electric field lines (abbreviated as E lines) are perpendicular to the parallel plate, which is different from the literatures, and the dominant mode distribution is irrelevant to the groove width; for the guide with a narrow but deeper groove, the E lines are parallel to the plate region, which represents the characteristics of lower attenuation, for the guide with a deeper and wider groove, a complicated E lines are obtained. These results reveal that the dominant mode distribution in groove guide is varied, which replenish our understanding of groove guide, and the results have important values in design, manufacture, and application of groove guide.     

Two-Grid Discretization Scheme for Nonlinear Reaction Diffusion Equation by Mixed Finite Element Methods

In this paper, we study an efficient scheme for nonlinear reaction-diffusion equations discretized by mixed finite element methods. We mainly concern the case when pressure coefficients and source terms are nonlinear. To linearize the nonlinear mixed equations, we use the two-grid algorithm. We first solve the nonlinear equations on the coarse grid, then, on the fine mesh, we solve a linearized problem using Newton iteration once. It is shown that the algorithm can achieve asymptotically optimal approximation as long as the mesh sizes satisfy $H =\mathcal{O}(h^{\frac{1}{2}})$. As a result, solving such a large class of nonlinear equations will not be much more difficult than getting solutions of one linearized system.     

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.     

Finite Element Analysis of the Field Patterns in Elliptical Waveguide

In this paper, finite element method is first used to discuss the field patterns in a hollow perfectly conducting pipe with elliptical cross section. The electrical field lines of the dominant TE mode as well as the magnetic field lines of the lowest TM mode are presented. The results are the same with literatures, and also give a correction to the field pattern of the lowest TM mode in elliptical waveguide.     

有限体积元数值方法在大气污染模式中的应用

运用有限体积元方法分析求解大气污染模型问题,分别选取试探函数空间和检验函数空间为一次元函数空间和分片常数函数空间,并且给出L²估计和H¹估计,通过数值实验与有限差分方法进行分析与比较,说明其有效性.为改善大气污染问题的模拟提供实用有效的方法.     

Analysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods

In this paper, we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method. We use two Newton iterations on the fine grid in our methods. Firstly, we solve an original nonlinear problem on the coarse nonlinear grid, then we use Newton iterations on the fine grid twice. The two-grid idea is from Xu$'$s work [SIAM J. Numer. Anal., 33 (1996), pp. 1759-1777] on standard finite method. We also obtain the error estimates for the algorithms of the two-grid method. It is shown that the algorithm achieves asymptotically optimal approximation rate with the two-grid methods as long as the mesh sizes satisfy $h=\mathcal{O}(H^{(4k+1)/(k+1)})$.     

Uniform Convergence Analysis of a Higher Order Hybrid Stress Quadrilateral Finite Element Method for Linear Elasticity Problems

This paper derives a higher order hybrid stress finite element method on quadrilateral meshes for linear plane elasticity problems. The method employs continuous piecewise bi-quadratic functions in local coordinates to approximate the displacement vector and a piecewise-independent 15-parameter mode to approximate the stress tensor. Error estimation shows that the method is free from Poisson-locking and has second-order accuracy in the energy norm. Numerical experiments confirm the theoretical results.     

Discrete Maximum Principle Based on Repair Technique for Finite Element Scheme of Anisotropic Diffusion Problems

In this paper, we construct a global repair technique for the finite element scheme of anisotropic diffusion equations to enforce the repaired solutions satisfying the discrete maximum principle. It is an extension of the existing local repair technique. Both of the repair techniques preserve the total energy and are easy to be implemented. The numerical experiments show that these repair techniques do not destroy the accuracy of the finite element scheme, and the computational cost of the global repair technique is lower than the local repair technique when the diffusion tensors are highly anisotropic.