A stabilized finite element method based on local polynomial pressure projection for the stationary Navier–Stokes equations

This article considers a stabilized finite element approximation for the branch of nonsingular solutions of the stationary Navier–Stokes equations based on local polynomial pressure projection by using the lowest equal-order elements. The proposed stabilized method has a number of attractive computational properties. Firstly, it is free from stabilization parameters. Secondly, it only requires the simple and efficient calculation of Gauss integral residual terms. Thirdly, it can be implemented at the element level. The optimal error estimate is obtained by the standard finite element technique. Finally, comparison with other methods, through a series of numerical experiments, shows that this method has better stability and accuracy.     

Patch‐wise local projection stabilized finite element methods for convection–diffusion–reaction problems

In this article, we develop patch‐wise local projection‐stabilized conforming and nonconforming finite element methods for the convection–diffusion–reaction problems. It is a composition of the standard Galerkin finite element method, the patch‐wise local projection stabilization, and weakly imposed Dirichlet boundary conditions on the discrete solution. In this paper, a priori error analysis is established with respect to a patch‐wise local projection norm for the conforming and the nonconforming finite element methods. The numerical experiments confirm the efficiency of the proposed stabilization technique and validate the theoretical convergence rates.     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

The Cauchy–Riemann Equations: Discretization by Finite Elements,Fast Solution of the Second Variable,and a Posteriori Error Estimation

In this paper we will concentrate on the numerical solution of the Cauchy–Riemann equations. First we show that these equations bring together the finite element discretizations for the Laplace equation by standard finite elements on the one hand, and by mixed finite element methods on the other. As a consequence, methods for a posteriori error estimation for both finite element methods can derive their validity from each other. Moreover, we show that given a finite element approximation of one of the vectorfields, the missing can be accurately computed in optimal complexity.     

Vibration of a Reissner–Mindlin–Timoshenko plate–beam system

In this paper, we consider a plate–beam system in which the Reissner–Mindlin plate model is combined with the Timoshenko beam model. Natural frequencies and vibration modes for the system are calculated using the finite element method. The interface conditions at the contact between the plate and beams are discussed in some detail. The impact of regularity on the enforcement of certain interface conditions is an important feature of the paper.     

A posteriori error estimation for the Stokes–Darcy coupled problem on anisotropic discretization

This paper presents an a posteriori error analysis for the stationary Stokes–Darcy coupled problem approximated by finite element methods on anisotropic meshes in or 3. Korn's inequality for piecewise linear vector fields on anisotropic meshes is established and is applied to non‐conforming finite element method. Then the existence and uniqueness of the approximation solution are deduced for non‐conforming case. With the obtained finite element solutions, the error estimators are constructed and based on the residual of model equations plus the stabilization terms. The lower error bound is proved by means of bubble functions and the corresponding anisotropic inverse inequalities. In order to prove the upper error bound, it is vital that an anisotropic mesh corresponds to the anisotropic function under consideration. To measure this correspondence, a so‐called matching function is defined, and its discussion shows it to be useful tool. With its help, the upper error bound is shown by means of the corresponding anisotropic interpolation estimates and a special Helmholtz decomposition in both media. Copyright © 2016 John Wiley & Sons, Ltd.     

Asymptotic Gauss–Jacobi quadrature error estimation for Schwarz–Christoffel integrals

Numerical conformal mapping packages based on the Schwarz–Christoffel formula have been in existence for a number of years. Various authors, for good reasons of practical efficiency, have chosen to use composite n-point Gauss–Jacobi rules for the estimation of the Schwarz–Christoffel path integrals. These implementations rely on an ad hoc, but experimentally well-founded, heuristic for selecting the spacing of the integration end-points relative to the position of the nearby integrand singularities. In the present paper we derive an explicitly computable estimate, asymptotic as n→∞, for the relevant Gauss–Jacobi quadrature error. A numerical example illustrates the potential accuracy of the estimate even at low values of n. It is apparent that the error estimate will allow the adaptive construction of composite rules in a manner that is more efficient than has been possible hitherto.     

Error analysis of a fully discrete finite element variational multiscale method for time‐dependent incompressible Navier–Stokes equations

A finite element variational multiscale method based on two local Gauss integrations is applied to solve numerically the time‐dependent incompressible Navier–Stokes equations. A significant feature of the method is that the definition of the stabilization term is derived via two local Guass integrations at element level, making it more efficient than the usual projection‐based variational multiscale methods. It is computationally cheap and gives an accurate approximation to the quantities sought. Based on backward Euler and Crank–Nicolson schemes for temporal discretization, we derive error bounds of the fully discrete solution which are first and second order in time, respectively. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

A posteriori error estimate for the Stokes–Darcy system

We consider the a posteriori error estimates for finite element approximations of the Stokes–Darcy system. The finite element spaces adopted are the Hood–Taylor element for the velocity and the pressure in fluid region and conforming piecewise quadratic element for the pressure in porous media region. The a posteriori error estimate is based on a suitable evaluation on the residual of the finite element solution. It is proven that the a posteriori error estimate provided in this paper is both reliable and efficient. Copyright © 2011 John Wiley & Sons, Ltd.     

A nonconforming rectangular finite element pair for the Darcy–Stokes–Brinkman model

In this article, we consider the Darcy–Stokes–Brinkman model that can be sorted into three problems: the Darcy problems, the Stokes–Brinkman interface problems and the coupled Darcy–Stokes problems. We study finite element approximation of the model with Dirichlet boundary conditions and make a unified analysis of the three problems based on nonconforming element. Optimal error estimates for the fluid velocity and pressure are derived. Finally, we present some numerical examples verifying the theoretical predictions. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Iterative methods in penalty finite element discretizations for the Steady Navier–Stokes equations

Three penalty finite element methods are designed to solve numerically the steady Navier–Stokes equations, where the Stokes, Newton, and Oseen iteration methods are used, respectively. Moreover, the stability analysis and error estimate for these nine algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms. Meanwhile, the numerical investigations are provided to show that the proposed methods are efficient for solving the steady Navier–Stokes equations with the different viscosity. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 74‐94, 2014     

Decoupled modified characteristics finite element method for the time dependent Navier–Stokes/Darcy problem

In this paper, a modified characteristics finite element method for the time dependent Navier–Stokes/Darcy problem with the Beavers–Joseph–Saffman interface condition is presented. In this method, the Navier–Stokes/Darcy equation is decoupled into two equations, one is the Navier–Stokes equation, the other is the Darcy equation, and the Navier–Stokes equation is solved by the modified characteristics finite element method. The theory analysis shows that this method has a good convergence property. In order to show the effect of our method, some numerical results was presented. The numerical results show that this method is highly efficient. Copyright © 2013 John Wiley & Sons, Ltd.     

A Priori Error Estimates of Crank–Nicolson Finite Volume Element Method for a Hyperbolic Optimal Control Problem

In this article, a Crank–Nicolson linear finite volume element scheme is developed to solve a hyperbolic optimal control problem. We use the variational discretization technique for the approximation of the control variable. The optimal convergent order O(h² + k²) is proved for the numerical solution of the control, state and adjoint‐state in a discrete L²‐norm. To derive this result, we also get the error estimate (convergent order O(h² + k²)) of Crank–Nicolson finite volume element approximation for the second‐order hyperbolic initial boundary value problem. Numerical experiments are presented to verify the theoretical results.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1331–1356, 2016     

Composite generalized Laguerre–Legendre pseudospectral method for Fokker–Planck equation in an infinite channel

In this paper, we propose a composite generalized Laguerre–Legendre pseudospectral method for the Fokker–Planck equation in an infinite channel, which behaves like a parabolic equation in one direction, and behaves like a hyperbolic equation in other direction. We establish some approximation results on the composite generalized Laguerre–Legendre–Gauss–Radau interpolation, with which the convergence of proposed composite scheme follows. An efficient implementation is provided. Numerical results show the spectral accuracy in space of this approach and coincide well with theoretical analysis. The approximation results and techniques developed in this paper are also very appropriate for many other problems on multiple-dimensional unbounded domains, which are not of standard types.     

Ciarlet–Raviart mixed finite element approximation for an optimal control problem governed by the first bi-harmonic equation

The Ciarlet–Raviart mixed finite element approximation is constructed to solve the constrained optimal control problem governed by the first bi-harmonic equation. The optimality conditions consisting of the state and the co-state equations is derived. Also, the a priori error estimates are analyzed. In the analysis of the a priori error estimates, the improved convergent rate of the higher order than existed results is proved. Some numerical experiments are performed to confirm the theoretical analysis for the a priori error estimate.     

A two-grid method based on Newton iteration for the Navier–Stokes equations

In this paper, we consider a two-grid method for resolving the nonlinearity in finite element approximations of the equilibrium Navier–Stokes equations. We prove the convergence rate of the approximation obtained by this method. The two-grid method involves solving one small, nonlinear coarse mesh system and two linear problems on the fine mesh which have the same stiffness matrix with only different right-hand side. The algorithm we study produces an approximate solution with the optimal asymptotic in h and accuracy for any Reynolds number. Numerical example is given to show the convergence of the method.     

A posteriori error analysis of a cell-centered finite volume method for semilinear elliptic problems

In this paper, we conduct a goal-oriented a posteriori analysis for the error in a quantity of interest computed from a cell-centered finite volume scheme for a semilinear elliptic problem. The a posteriori error analysis is based on variational analysis, residual errors and the adjoint problem. To carry out the analysis, we use an equivalence between the cell-centered finite volume scheme and a mixed finite element method with special choice of quadrature.     

Bubble—函数稳定化混合有限元分析

本文针对混合结构抽象问题，基于「９」的非标准稳定化有限元方法的一般框架研究了ｂｕｂｂｌｅ－函数稳定化方法，该逼近代格式使得Ｂａｂｕｓｋａ－Ｂｒｅｚｚｉ条件是不必要的。     

An adaptive algorithm for the Crank–Nicolson scheme applied to a time-dependent convection–diffusion problem

An a posteriori upper bound is derived for the nonstationary convection–diffusion problem using the Crank–Nicolson scheme and continuous, piecewise linear stabilized finite elements with large aspect ratio. Following Lozinski et al. (2009) [13], a quadratic time reconstruction is used.A space and time adaptive algorithm is developed to ensure the control of the relative error in the L²(H¹) norm. Numerical experiments illustrating the efficiency of this approach are reported; it is shown that the error indicator is of optimal order with respect to both the mesh size and the time step, even in the convection dominated regime and in the presence of boundary layers.