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

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

非定常Navier-Stokes方程基于完全重叠型区域分解的有限元并行算法

基于完全重叠型区域分解技巧,提出三种求解非定常Navier-Stokes方程的有限元并行算法.其基本思想是首先对空间施行完全重叠区域分解,然后各个处理器使用向后Euler格式独立并行求解关于时间t的常微分方程;对于非线性的对流项,分别采用半隐格式和全隐格式进行处理.算法中每个处理器所负责的子问题是一个全局问题,它定义在整个求解区域上,但绝大部分自由度来自其所负责的子区域,从而使得算法实现简单,通信需求少.数值算例验证了算法的有效性及其良好的并行性能.     

粒子输运方程的线性间断有限元方法

将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象.     

高分辨率间断有限元方法

间断有限元方法是集高分辨率有限差分方法和有限体积方法的优点发展起来的一种数值方法,在计算流体动力学问题上显示了优良的效能.利用守恒问题给出间断有限元方法的基本概念和过程,利用简单算例给出该方法的精度分析和限制器对精度的影响,并给出浅水波问题、交通流问题和波传播问题的数值模拟结果,进一步,综合评介该方法在椭圆、抛物、对流扩散、Hamilton-Jacobi方程、Navier-Stokes方程等的实际应用进展.     

间断有限元方法在弹尾超音速喷流计算中的应用

采用间断有限元方法对超音速无粘喷流流动进行数值模拟.将二维双曲守恒方程的间断有限元方法发展到轴对称Euler方程,并就某导弹尾部超音速伴随射流进行数值计算.计算结果与实验照片反映的流动特征吻合较好,与高精度、高分辨率TVD格式的计算结果相比,间断有限元方法的计算结果在轴线反射点附近具有较高的分辨率,表明该方法对激波具有较强的捕捉能力,在激波阵面上不会产生振荡或抹平间断现象.     

一阶双曲型方程的AGE方法

采用数值边界条件的办法导出了求解一阶双曲型方程的AGE方法,并给出了其稳定性证明和计算实例。     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

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.     

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.     

Moving Finite Element Methods for a System of Semi-Linear Fractional Diffusion Equations

This paper studies a system of semi-linear fractional diffusion equations which arise in competitive predator-prey models by replacing the second-order derivatives in the spatial variables with fractional derivatives of order less than two. Moving finite element methods are proposed to solve the system of fractional diffusion equations and the convergence rates of the methods are proved. Numerical examples are carried out to confirm the theoretical findings. Some applications in anomalous diffusive Lotka-Volterra and Michaelis-Menten-Holling predator-prey models are studied.     

Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes I: Finite Element Solutions

In this paper we developed accurate finite element methods for solving 3-D Poisson-Nernst-Planck (PNP) equations with singular permanent charges for electrodiffusion in solvated biomolecular systems. The electrostatic Poisson equation was defined in the biomolecules and in the solvent, while the Nernst-Planck equation was defined only in the solvent. We applied a stable regularization scheme to remove the singular component of the electrostatic potential induced by the permanent charges inside biomolecules, and formulated regular, well-posed PNP equations. An inexact-Newton method was used to solve the coupled nonlinear elliptic equations for the steady problems; while an Adams-Bashforth-Crank-Nicolson method was devised for time integration for the unsteady electrodiffusion. We numerically investigated the conditioning of the stiffness matrices for the finite element approximations of the two formulations of the Nernst-Planck equation, and theoretically proved that the transformed formulation is always associated with an ill-conditioned stiffness matrix. We also studied the electroneutrality of the solution and its relation with the boundary conditions on the molecular surface, and concluded that a large net charge concentration is always present near the molecular surface due to the presence of multiple species of charged particles in the solution. The numerical methods are shown to be accurate and stable by various test problems, and are applicable to real large-scale biophysical electrodiffusion problems.     

自适应间断有限元方法求解三维欧拉方程

将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性.     

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.     

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.     

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.     

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)})$.     

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.     

Superconvergence of Mixed Methods for Optimal Control Problems Governed by Parabolic Equations

In this paper, we investigate the superconvergence results for optimal control problems governed by parabolic equations with semi-discrete mixed finite element approximation. We use the lowest order mixed finite element spaces to discrete the state and costate variables while use piecewise constant function to discrete the control variable. Superconvergence estimates for both the state variable and its gradient variable are obtained.     

Hyperbolic Octonion Formulation of Gravitational Field Equations

In this paper, the Maxwell-Proca type field equations of linear gravity are formulated in terms of hyperbolic octonions (split octonions). A hyperbolic octonionic gravitational wave equation with massive gravitons and gravitomagnetic monopoles is proposed. The real gravitoelectromagnetic field equations are recovered and written in compact form from an octonionic potential. In the absence of charges, this reduces to the Klein-Gordon equation of motion for the massive graviton. The analogy between massive gravitational theory and electromagnetism is shown in terms of the present formulation.     

Collocation Methods for Hyperbolic Partial Differential Equations with Singular Sources

A numerical study is given on the spectral methods and the high order WENO finite difference scheme for the solution of linear and nonlinear hyperbolic partial differential equations with stationary and non-stationary singular sources. The singular source term is represented by the $δ$-function. For the approximation of the $δ$-function, the direct projection method is used that was proposed in [6]. The $δ$-function is constructed in a consistent way to the derivative operator. Nonlinear sine-Gordon equation with a stationary singular source was solved with the Chebyshev collocation method. The $δ$-function with the spectral method is highly oscillatory but yields good results with small number of collocation points. The results are compared with those computed by the second order finite difference method. In modeling general hyperbolic equations with a non-stationary singular source, however, the solution of the linear scalar wave equation with the non-stationary singular source using the direct projection method yields non-physical oscillations for both the spectral method and the WENO scheme. The numerical artifacts arising when the non-stationary singular source term is considered on the discrete grids are explained.