Modified two-grid method for solving coupled Navier-Stokes/Darcy model based on Newton iteration

A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the err...     

Guaranteed and robust discontinuous Galerkin a posteriori error estimates for convection-diffusion-reaction problems

We propose and study a posteriori error estimates for convection-diffusion-reaction problems with inhomogeneous and anisotropic diffusion approximated by weighted interior-penalty discontinuous Galerkin methods. Our twofold objective is to derive estimates without undetermined constants and to analyze carefully the robustness of the estimates in singularly perturbed regimes due to dominant convection or reaction. We first derive locally computable estimates for the error measured in the energy (semi)norm. These estimates are evaluated using -conforming diffusive and convective flux reconstructions, thereby extending the previous work on pure diffusion problems. The resulting estimates are semi-robust in the sense that local lower error bounds can be derived using suitable cutoff functions of the local Péclet and Damköhler numbers. Fully robust estimates are obtained for the error measured in an augmented norm consisting of the energy (semi)norm, a dual norm of the skew-symmetric part of the differential operator, and a suitable contribution of the interelement jumps of the discrete solution. Numerical experiments are presented to illustrate the theoretical results.     

A posteriori error estimates for mixed finite element approximations of elliptic problems

We derive residual based a posteriori error estimates of the flux in L ²-norm for a general class of mixed methods for elliptic problems. The estimate is applicable to standard mixed methods such as the Raviart–Thomas–Nedelec and Brezzi–Douglas–Marini elements, as well as stabilized methods such as the Galerkin-Least squares method. The element residual in the estimate employs an elementwise computable postprocessed approximation of the displacement which gives optimal order.     

On computing the pressure by thep version of the finite element method for Stokes problem

Summary This paper introduces and analyzes two ways of extracting the hydrostatic pressure when solving Stokes problem using thep version of the finite element method. When one uses a localH ¹ projection, we show that optimal rates of convergence for the pressure approximation is achieved. When the pressure is not inH ¹. or the value of the pressure is only needed at a few points, one may extract the pressure pointwise using e.g. a single layer potential recovery. Negative, zero, and higher norm estimates for the Stokes velocity are derived within the framework of thep version of the F.E.M.Partially supported by ONR grants N00014-87-K-0427 and N00014-90-J-1238     

Local and parallel finite element algorithms for the stokes problem

Based on two-grid discretizations, some local and parallel finite element algorithms for the Stokes problem are proposed and analyzed in this paper. These algorithms are motivated by the observation that for a solution to the Stokes problem, low frequency components can be approximated well by a relatively coarse grid and high frequency components can be computed on a fine grid by some local and parallel procedure. One technical tool for the analysis is some local a priori estimates that are also obtained in this paper for the finite element solutions on general shape-regular grids. Y. He was partially subsidized by the NSF of China 10671154 and the National Basic Research Program under the grant 2005CB321703; A. Zhou was partially supported by the National Science Foundation of China under the grant 10425105 and the National Basic Research Program under the grant 2005CB321704; J. Li was partially supported by the NSF of China under the grant 10701001. J. Xu was partially supported by Alexander von Humboldt Research Award for Senior US Scientists, NSF DMS-0609727 and NSFC-10528102.     

A penalty finite element method based on the Euler implicit/explicit scheme for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations, where the time discretization is based on the Euler implicit/explicit scheme with some implicit linear terms and an explicit nonlinear term, and the finite element spatial discretization is based on the P₁b-P₁ element pair, which satisfies the discrete inf-sup condition. This method allows us to separate the computation of the velocity from the computation of the pressure with a larger time-step size Δt, so that the numerical velocity and the pressure are easily computed. An optimal error estimate of the numerical velocity and the pressure is provided for the fully discrete penalty finite element method when the penalty parameter ?, the time-step size Δt and the mesh size h satisfy the following stability conditions: ?c₁≤1, Δtκ₁≤1 and h²≤β₁Δt, respectively, for some positive constants c₁, κ₁ and β₁. Finally, some numerical tests to confirm the theoretical results of the penalty finite element method are provided.     

A framework for obtaining guaranteed error bounds for finite element approximations

We give an overview of our recent progress in developing a framework for the derivation of fully computable guaranteed posteriori error bounds for finite element approximation including conforming, non-conforming, mixed and discontinuous finite element schemes. Whilst the details of the actual estimator are rather different for each particular scheme, there is nonetheless a common underlying structure at work in all cases. We aim to illustrate this structure by treating conforming, non-conforming and discontinuous finite element schemes in a single framework. In taking a rather general viewpoint, some of the finer details of the analysis that rely on the specific properties of each particular scheme are obscured but, in return, we hope to allow the reader to ‘see the wood despite the trees’.     

Analysis of a pressure-stabilized finite element approximation of the stationary Navier-Stokes equations

Summary. The purpose of this paper is to analyze a finite element approximation of the stationary Navier-Stokes equations that allows the use of equal velocity-pressure interpolation. The idea is to introduce as unknown of the discrete problem the projection of the pressure gradient (multiplied by suitable algorithmic parameters) onto the space of continuous vector fields. The difference between these two vectors (pressure gradient and projection) is introduced in the continuity equation. The resulting formulation is shown to be stable and optimally convergent, both in a norm associated to the problem and in the norm for both velocities and pressure. This is proved first for the Stokes problem, and then it is extended to the nonlinear case. All the analysis relies on an inf-sup condition that is much weaker than for the standard Galerkin approximation, in spite of the fact that the present method is only a minor modification of this. Received May 25, 1998 / Revised version received August 31, 1999 / Published online July 12, 2000     

A posteriori error estimates for a finite element approximation of transmission problems with sign changing coefficients

We perform the a posteriori error analysis of residual type of transmission problem with sign changing coefficients. According to Bonnet-BenDhia et al. (2010) [9], if the contrast is large enough, the continuous problem can be transformed into a coercive one. We further show that a similar property holds for the discrete problem for any regular meshes, extending the framework from Bonnet-BenDhia et al. [9]. The reliability and efficiency of the proposed estimator are confirmed by some numerical tests.     

Tuning the mesh of a mixed method for the stream function Vorticity formulation of the Navier-Stokes equations

Summary In this article we study a new mixed method for the Stokes and Navier-Stokes equations. The method uses two meshes, one very fine for and a coarser one for . Error estimates show that boundary layers do not require to refine the mesh for the stream function as much as for the vorticity when the Reynolds number is large. We prove estimates and study implementation problems.     

Error estimates for a mixed finite element discretization of some degenerate parabolic equations

We consider a numerical scheme for a class of degenerate parabolic equations, including both slow and fast diffusion cases. A particular example in this sense is the Richards equation modeling the flow in porous media. The numerical scheme is based on the mixed finite element method (MFEM) in space, and is of one step implicit in time. The lowest order Raviart–Thomas elements are used. Here we extend the results in Radu et al. (SIAM J Numer Anal 42:1452–1478, 2004), Schneid et al. (Numer Math 98:353–370, 2004) to a more general framework, by allowing for both types of degeneracies. We derive error estimates in terms of the discretization parameters and show the convergence of the scheme. The features of the MFEM, especially of the lowest order Raviart–Thomas elements, are now fully exploited in the proof of the convergence. The paper is concluded by numerical examples.     

A finite volume method for solving Navier-Stokes problems

We develop a finite volume method for solving the Navier-Stokes equations on a triangular mesh. We prove that the unique solution of the finite volume method converges to the true solution with optimal order for velocity and for pressure in discrete H¹ norm and L² norm respectively.     

A stabilized mixed finite element method for the biharmonic equation based on biorthogonal systems

We propose a stabilized finite element method for the approximation of the biharmonic equation with a clamped boundary condition. The mixed formulation of the biharmonic equation is obtained by introducing the gradient of the solution and a Lagrange multiplier as new unknowns. Working with a pair of bases forming a biorthogonal system, we can easily eliminate the gradient of the solution and the Lagrange multiplier from the saddle point system leading to a positive definite formulation. Using a superconvergence property of a gradient recovery operator, we prove an optimal a priori estimate for the finite element discretization for a class of meshes.     

A stabilized finite element method based on two local Gauss integrations for the Stokes equations

This paper considers a stabilized method based on the difference between a consistent and an under-integrated mass matrix of the pressure for the Stokes equations approximated by the lowest equal-order finite element pairs (i.e., the _P1$P_1$–_P1$P_1$ and _Q1$Q_1$–_Q1$Q_1$ pairs). 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. Optimal error estimates are obtained by applying the standard Galerkin technique. Finally, the numerical illustrations agree completely with the theoretical expectations.     

Pointwise error estimates for a streamline diffusion finite element scheme

Summary Pointwise error estimates for a streamline diffusion scheme for solving a model convection-dominated singularly perturbed convection-diffusion problem are given. These estimates improve pointwise error estimates obtained by Johnson et al.[5].     

Optimal error estimate of the penalty finite element method for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations. The time discretization of the penalty Navier-Stokes equations is based on the backward Euler scheme; the spatial discretization of the time discretized penalty Navier-Stokes equations is based on a finite element space pair which satisfies some approximate assumption. An optimal error estimate of the numerical velocity and pressure is provided for the fully discrete penalty finite element method when the parameters and are sufficiently small.

     

Stabilized discontinuous finite element approximations for Stokes equations

In this paper, we derive two stabilized discontinuous finite element formulations, symmetric and nonsymmetric, for the Stokes equations and the equations of the linear elasticity for almost incompressible materials. These methods are derived via stabilization of a saddle point system where the continuity of the normal and tangential components of the velocity/displacements are imposed in a weak sense via Lagrange multipliers. For both methods, almost all reasonable pair of discontinuous finite element spaces can be used to approximate the velocity and the pressure. Optimal error estimate for the approximation of both the velocity of the symmetric formulation and pressure in L²$L^2$ norm are obtained, as well as one in a mesh-dependent norm for the velocity in both symmetric and nonsymmetric formulations.     

The L.P.D.E.M., a mixed finite difference method for elliptic problems

Summary A mixed finite difference method is analyzed for solving certain elliptic problems. This method, called L.P.D.E.M. (Locally exact Partial Differential Equation Method) was initially proposed in the frame of hydrodynamic lubrication. Convergence is obtained. Relations between this scheme and homogenization theory are also discussed. For a one-dimensional elliptic equation with no zero-order term and in conservative form, this method is an exact one. Some numerical results will also be given.     

A posteriori error estimates for mixed finite element solutions of convex optimal control problems

In this paper, we present an a posteriori error analysis for mixed finite element approximation of convex optimal control problems. We derive a posteriori error estimates for the coupled state and control approximations under some assumptions which hold in many applications. Such estimates can be used to construct reliable adaptive mixed finite elements for the control problems.