Crank–Nicolson Finite Element Scheme and Modified Reduced-Order Scheme for Fractional Sobolev Equation

Abstract

In this paper, a Crank–Nicolson finite difference/finite element method is considered to obtain the numerical solution for a time fractional Sobolev equation. Firstly, the classical finite element method is presented. Stability and error estimation for the fully discrete scheme are rigorously established. However, the amount of calculation and computing time are too large due to many degrees of freedom of classical finite element scheme and nonlocality of fractional differential operator. And then the modified reduced-order finite element scheme with low dimensions and sufficiently high accuracy, which is based on proper orthogonal decomposition technique, is provided. Stability and convergence for the reduced-order scheme are also studied. At last, numerical examples show that the results of numerical computation are consistent with previous theoretical conclusions.     

A stabilized finite element method for the time-dependent Stokes equations based on Crank–Nicolson Scheme

A stabilized finite element method for the time-dependent Stokes equations based on Crank–Nicolson scheme is considered in this paper. The method combines the Crank–Nicolson scheme with a stabilized finite element method which uses the lowest equal-order element pair, i.e., the stabilized finite element method is applied for the spatial approximation and the time discretization is based on the Crank–Nicolson scheme. Moreover, we present optimal error estimates and prove that the scheme is unconditionally stable and convergent. Finally, numerical tests confirm the theoretical results of the presented method.     

Low order Crouzeix-Raviart type nonconforming finite element methods for the 3D time-dependent Maxwell’s equations

Two Crouzeix-Raviart type nonconforming elements are used in a finite element scheme as well in a mixed finite element scheme for time-dependent Maxwell’s equations in three dimensions. The error estimates are obtained under anisotropic meshes, which are the same as those for conforming elements under regular meshes.     

A low order anisotropic nonconforming characteristic finite element method for a convection-dominated transport problem

A low order anisotropic nonconforming rectangular finite element method for the convection-diffusion problem with a modified characteristic finite element scheme is studied in this paper. The O(h²) order error estimate in L²-norm with respect to the space, one order higher than the expanded characteristic-mixed finite element scheme with order O(h), and the same as the conforming case for a modified characteristic finite element scheme under regular meshes, is obtained by use of some distinct properties of the interpolation operator and the mean value technique, instead of the so-called elliptic projection, which is an indispensable tool in the convergence analysis of the previous literature. Lastly, some numerical results of the element are provided to verify our theoretical analysis.     

Finite volume element approximation and analysis for a kind of semiconductor device simulation

In this paper, we present a finite volume element scheme for a kind of two dimensional semiconductor device simulation. A general framework is developed for finite volume element approximation of the semiconductor problems. We construct a fully discrete finite volume element scheme based on triangulations with a piecewise linear finite element space and a general type of control volume. Optimal-order convergence in H ¹-norm is derived.     

A new combined finite element‐upwind finite volume method for convection‐dominated diffusion problems

In this article, we develop a combined finite element‐weighted upwind finite volume method for convection‐dominated diffusion problems in two dimensions, which discretizes the diffusion term with the standard finite element scheme, and the convection and source terms with the weighted upwind finite volume scheme. The developed method leads to a totally new scheme for convection‐dominated problems, which overcomes numerical oscillation, avoids numerical dispersion, and has high‐order accuracy. Stability analyses of the scheme are given for the problems with constant coefficients. Numerical experiments are presented to illustrate the stability and optimal convergence of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 799–818, 2016     

A Multilevel Correction Method for Steklov Eigenvalue Problem by Nonconforming Finite Element Methods

In this paper, a multilevel correction scheme is proposed to solve the Steklov eigenvalue problem by nonconforming finite element methods. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an Steklov eigenvalue problem on the coarsest finite element space. This correction scheme can increase the overall efficiency of solving eigenvalue problems by the nonconforming finite element method. Furthermore, as same as the direct eigenvalue solving by nonconforming finite element methods, this multilevel correction method can also produce the lower-bound approximations of the eigenvalues.     

An Economical Finite Element Scheme for Navier-Stokes Equations

In this paper, a new finite element scheme for Navier-Stokes equations is proposed, in which three different partitions (in the two dimensional case) are used to construct finite element subspaces of the velocity field and the pressure. The error estimate of the finite element than approximation is given. The precision of this new scheme has the same order as the scheme $Q_2/P_0$, but it is more economical that the scheme $Q_2/P_0$.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

A monotone scheme for Hamilton–Jacobi equations via the nonstandard finite difference method

A usual way of approximating Hamilton–Jacobi equations is to couple space finite element discretization with time finite difference discretization. This classical approach leads to a severe restriction on the time step size for the scheme to be monotone. In this paper, we couple the finite element method with the nonstandard finite difference method, which is based on Mickens' rule of nonlocal approximation. The scheme obtained in this way is unconditionally monotone. The convergence of the new method is discussed and numerical results that support the theory are provided. Copyright © 2009 John Wiley & Sons, Ltd.     

An unfitted interface penalty method for the numerical approximation of contrast problems

We aim to approximate contrast problems by means of a numerical scheme which does not require that the computational mesh conforms with the discontinuity between coefficients. We focus on the approximation of diffusion-reaction equations in the framework of finite elements. In order to improve the unsatisfactory behavior of Lagrangian elements for this particular problem, we resort to an enriched approximation space, which involves elements cut by the interface. Firstly, we analyze the H¹-stability of the finite element space with respect to the position of the interface. This analysis, applied to the conditioning of the discrete system of equations, shows that the scheme may be ill posed for some configurations of the interface. Secondly, we propose a stabilization strategy, based on a scaling technique, which restores the standard properties of a Lagrangian finite element space and results to be very easily implemented. We also address the behavior of the scheme with respect to large contrast problems ending up with a choice of Nitsche?s penalty terms such that the extended finite element scheme with penalty is robust for the worst case among small sub-elements and large contrast problems. The theoretical results are finally illustrated by means of numerical experiments.     

Numerical studies on a novel split-step quadratic B-spline finite element method for the coupled Schrödinger–KdV equations

We propose a novel split-step quadratic B-spline finite element method for solving the initial-boundary value problem of the coupled Schrödinger–KdV equations. A full-discrete finite element scheme is constructed. The conserved properties of the full-discrete scheme are proved. Detailed numerical results show the efficiency of our scheme.     

A fast and robust numerical scheme for solving models of charge carrier transport and ion vacancy motion in perovskite solar cells

Drift-diffusion models that account for the motion of ion vacancies and electronic charge carriers are important tools for explaining the behaviour, and guiding the development, of metal halide perovskite solar cells. Computing numerical solutions to such models in realistic parameter regimes, where the short Debye lengths give rise to boundary layers in which the solution varies extremely rapidly, is challenging. Two suitable numerical methods, that can effectively cope with the spatial stiffness inherent to such problems, are presented and contrasted (a finite element scheme and a finite difference scheme). Both schemes are based on an appropriate choice of non-uniform spatial grid that allows the solution to be computed accurately in the boundary layers. An adaptive time step is employed in order to combat a second source of stiffness, due to the disparity in timescales between the motion of the ion vacancies and electronic charge carriers. It is found that the finite element scheme provides significantly higher accuracy, in a given compute time, than both the finite difference scheme and some previously used alternatives (Chebfun and pdepe). An example transient sweep of a current-voltage curve for realistic parameter values can be computed using this finite element scheme in only a few seconds on a standard desktop computer.     

Finite-element method for time-dependent euler equation

This paper presents finite element methods to approximate inviscid incompressible flow problems. First we emphasize the conservation properties of these problems, and we show that finite element methods appear as a very natural way to find conservative schemes such as Arakawa's scheme. We give convergence theorems and an error analysis of finite element discretization schemes. We turn then to the time differencing problem. We derive stability and convergence results for a second-order semi-implicit scheme and for the leap-frog scheme.     

板的低阶间断与连续有限元法统一分析

基于Reissner-Mindlin板问题的间断Galerkin有限元逼近, 建立了一个对挠度空间和角位移空间取连续或间断元都适用的低阶有限元离散格式. 取剪切力空间为分片常数元, 挠度空间和角位移空间无论取间断元还是连续元, 格式都是一致稳定的, 并给出了H¹范数估计及L²范数估计. 作为应用,对几类低阶有限元空间讨论. 结果表明, 该格式对常见的低阶有限元空间都适用, 并且若至少有一个元连续时, 该格式需要的空间比[1,2]中的都要简单.       

Node-pair finite volume/finite element schemes for the Euler equation in cylindrical and spherical coordinates

A numerical scheme is presented for the solution of the compressible Euler equations in both cylindrical and spherical coordinates. The unstructured grid solver is based on a mixed finite volume/finite element approach. Equivalence conditions linking the node-centered finite volume and the linear Lagrangian finite element scheme over unstructured grids are reported and used to devise a common framework for solving the discrete Euler equations in both the cylindrical and the spherical reference systems. Numerical simulations are presented for the explosion and implosion problems with spherical symmetry, which are solved in both the axial–radial cylindrical coordinates and the radial–azimuthal spherical coordinates. Numerical results are found to be in good agreement with one-dimensional simulations over a fine mesh.     

Combined higher order finite volume and finite element scheme for double porosity and non-linear adsorption of transport problem in porous media

In this work, a dual porosity model of reactive solute transport in porous media is presented. This model consists of a nonlinear-degenerate advection-diffusion equation including equilibrium adsorption to the reaction combined with a first-order equation for the non-equilibrium adsorption interaction processes. The numerical scheme for solving this model involves a combined high order finite volume and finite element scheme for approximation of the advection-diffusion part and relaxation-regularized algorithm for nonlinearity-degeneracy. The combined finite volume-finite element scheme is based on a new formulation developed by Eymard et al. (2010) [10]. This formulation treats the advection and diffusion separately. The advection is approximated by a second-order local maximum principle preserving cell-vertex finite volume scheme that has been recently proposed whereas the diffusion is approximated by a finite element method. The result is a conservative, accurate and very flexible algorithm which allows the use of different mesh types such as unstructured meshes and is able to solve difficult problems. Robustness and accuracy of the method have been evaluated, particularly error analysis and the rate of convergence, by comparing the analytical and numerical solutions for first and second order upwind approaches. We also illustrate the performance of the discretization scheme through a variety of practical numerical examples. The discrete maximum principle has been proved.     

Solution of double nonlinear problems in porous media by a combined finite volume–finite element algorithm

The combined finite volume–finite element scheme for a double nonlinear parabolic convection-dominated diffusion equation which models the variably saturated flow and contaminant transport problems in porous media is extended. Whereas the convection is approximated by a finite volume method (Multi-Point Flux Approximation), the diffusion is approximated by a finite element method. The scheme is fully implicit and involves a relaxation-regularized algorithm. Due to monotonicity and conservation properties of the approximated scheme and in view of the compactness theorem we show the convergence of the numerical scheme to the weak solution. Our scheme is applied for computing two dimensional examples with different degrees of complexity. The numerical results demonstrate that the proposed scheme gives good performance in convergence and accuracy.     

两个带变时间步的两重网格算法的稳定性分析

对用于求解非线性发展方程的两个带变时间步的两重网格算法，对空间变量用有限元离散，对时间变量分别用一阶精度Euler显式和二阶精度半隐式差分格式离散，然后构造两重网格算法，通过深入的稳定性分析，得出本文的算法优于标准全离散有限元算法。     

自然对流换热问题基于混合元法的差分格式及其数值模拟

研究自然对流换热问题,通过对于空间变量采用有限元离散而对于时间变量用差分离散,导出一种基于混合有限元法的最低阶的差分格式,这种格式可以同时求出流体的速度、温度和压力的数值解,并给出了模拟方腔流的自然换热的数值例子。