Least‐squares mixed finite element methods for the incompressible Navier‐Stokes equations

Least‐squares mixed finite element schemes are formulated to solve the evolutionary Navier‐Stokes equations and the convergence is analyzed. We recast the Navier‐Stokes equations as a first‐order system by introducing a vorticity flux variable, and show that a least‐squares principle based on L² norms applied to this system yields optimal discretization error estimates. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 441–453, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10015     

Error estimates for mixed finite element approximations of the viscoelasticity wave equation

This paper studies mixed finite element approximations to the solution of the viscoelasticity wave equation. Two new transformations are introduced and a corresponding system of first‐order differential‐integral equations is derived. The semi‐discrete and full‐discrete mixed finite element methods are then proposed for the problem based on the Raviart–Thomas–Nedelec spaces. The optimal error estimates in L²‐norm are obtained for the semi‐discrete and full‐discrete mixed approximations of the general viscoelasticity wave equation. Copyright © 2004 John Wiley & Sons, Ltd.     

Symmetric mixed covolume methods for parabolic problems

We present a mixed covolume method for a system of first order partial differential equations resulting from the mixed formulation of the general self‐adjoint parabolic problem with a variable nondiagonal diffusion tensor. The lowest order Raviart‐Thomas mixed element space on rectangles is used. We prove the first order optimal rate of convergence for approximate pressure as well as for approximate velocity. We also prove the second order superconvergence both for approximate velocity and pressure in certain discrete norms. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 561–583, 2002     

Mixed covolume method for parabolic problems on triangular grids

In this paper, we present a mixed covolume method for parabolic equations on triangular grids. This method use the lowest order Raviart–Thomas (R–T) mixed finite element space as the trial space. We prove the optimal order of convergence for the approximate pressure and velocity in L²-norm. Furthermore, we obtain the quasi-optimal error estimates for the approximate pressure in L^∞-norm.     

A splitting positive definite mixed element method for miscible displacement of compressible flow in porous media

A miscible displacement of one compressible fluid by another in a porous medium is governed by a nonlinear parabolic system. A new mixed finite element method, in which the mixed element system is symmetric positive definite and the flux equation is separated from pressure equation, is introduced to solve the pressure equation of parabolic type, and a standard Galerkin method is used to treat the convection‐diffusion equation of concentration of one of the fluids. The convergence of the approximate solution with an optimal accuracy in L²‐norm is proved. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 229–249, 2001     

An implicit–explicit residual error estimator for the coupling of dual‐mixed finite elements and boundary elements in elastostatics

We consider the coupling of dual‐mixed finite elements and boundary elements to solve a mixed Dirichlet–Neumann problem of plane elasticity. We derive an a‐posteriori error estimate that is based on the solution of local Dirichlet problems and on a residual term defined on the coupling interface. The general error estimate does not make use of any special finite element or boundary element spaces. Here the residual term is given in a negative order Sobolev norm. In practical applications, where a certain boundary element subspace is used, this norm can be estimated by weighted local L²‐norms. Copyright © 2001 John Wiley & Sons, Ltd.     

An optimal‐order error estimate on an H1‐Galerkin mixed method for a nonlinear parabolic equation in porous medium flow

We present an H¹‐Galerkin mixed finite element method for a nonlinear parabolic equation, which models a compressible fluid flow process in subsurface porous media. The method possesses the advantages of mixed finite element methods while avoiding directly inverting the permeability tensor, which is important especially in a low permeability zone. We conducted theoretical analysis to study the existence and uniqueness of the numerical solutions of the scheme and prove an optimal‐order error estimate for the method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Superconvergence of H 1-Galerkin mixed finite element methods for parabolic problems

In this article, we study the semidiscrete H ¹-Galerkin mixed finite element method for parabolic problems over rectangular partitions. The well-known optimal order error estimate in the L ²-norm for the flux is of order 𝒪(h ^k+1) (SIAM J. Numer. Anal. 35 (2), (1998), pp. 712–727), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas element. We derive a superconvergence estimate of order 𝒪(h ^k+3) between the H ¹-Galerkin mixed finite element approximation and an appropriately defined local projection of the flux variable when k ≥ 1. A the new approximate solution for the flux with superconvergence of order 𝒪(h ^k+3) is realized via a postprocessing technique using local projection methods.     

Error estimates of expanded mixed methods for optimal control problems governed by hyperbolic integro‐differential equations

In this article, we investigate the L^∞(L²) ‐error estimates of the semidiscrete expanded mixed finite element methods for quadratic optimal control problems governed by hyperbolic integrodifferential equations. The state and the costate are discretized by the order k Raviart‐Thomas mixed finite element spaces, and the control is approximated by piecewise polynomials of order k(k ≥ 0). We derive error estimates for both the state and the control approximation. Numerical experiments are presented to test the theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Convergence analysis of an upwind mixed element method for advection diffusion problems

We consider a upwinding mixed element method for a system of first order partial differential equations resulting from the mixed formulation of a general advection diffusion problem. The system can be used to model the transport of a contaminant carried by a flow. We use the lowest order Raviart-Thomas mixed finite element space. We show the first order convergence both for concentration and concentration flux in L²(Ω).     

An approximation of incompressible miscible displacement in porous media by mixed finite elements and symmetric finite volume element method of characteristics

A nonlinear system of two coupled partial differential equations models miscible displacement of one incompressible fluid by another in a porous medium. A sequential implicit time‐stepping procedure is defined, in which the pressure and Darcy velocity of the mixture are approximated by a mixed finite element method and the concentration is approximated by a combination of a modified symmetric finite volume element method and the method of characteristics. Optimal order convergence in H¹ and in L² are proved for full discrete schemes. Finally, some numerical experiments are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A combined hybrid mixed element method for incompressible miscible displacement problem with local discontinuous Galerkin procedure

In this article, we propose a combined hybrid discontinuous mixed finite element method for miscible displacement problem with local discontinuous Galerkin method. Here, to obtain more accurate approximation and deal with the discontinuous case, we use the hybrid mixed element method to approximate the pressure and velocity, and use the local discontinuous Galerkin finite element method for the concentration. Compared with other combined methods, this method can improve the efficiency of computation, deal with the discontinuous problem well and keep local mass balance. We study the convergence of this method and give the corresponding optimal error estimates in L^∞(L²) for velocity and concentration and the super convergence in L^∞(H¹) for pressure. Finally, we also present some numerical examples to confirm our theoretical analysis.     

Two-grid methods for –P1 mixed finite element approximation of general elliptic optimal control problems with low regularity

In this paper, we present a two-grid mixed finite element scheme for distributed optimal control governed by general elliptic equations. –P₁ mixed finite elements are used for the discretization of the state and co-state variables, whereas piecewise constant function is used to approximate the control variable. We first use a new approach to obtain the superclose property between the centroid interpolation and the numerical solution of the optimal control u with order h² under the low regularity. Based on the superclose property, we derive the optimal a priori error estimates. Then, using a postprocessing projection operator, we get a second-order superconvergent result for the control u. Next, we construct a two-grid mixed finite element scheme and analyze a priori error estimates. In the two-grid scheme, the solution of the elliptic optimal control problem on a fine grid is reduced to the solution of the elliptic optimal control problem on a much coarser grid and the solution of a linear algebraic system on the fine grid and the resulting solution still maintains an asymptotically optimal accuracy. Finally, a numerical example is presented to verify the theoretical results.     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

Unconditional superconvergent analysis of a new mixed finite element method for Ginzburg–Landau equation

In this article, unconditional superconvergent analysis of a linearized fully discrete mixed finite element method is presented for a class of Ginzburg–Landau equation based on the bilinear element and zero‐order Nédélec's element pair (Q₁₁/Q₀₁ × Q₁₀). First, a time‐discrete system is introduced to split the error into temporal error and spatial error, and the corresponding error estimates are deduced rigorously. Second, the unconditional superclose and optimal estimate of order O(h² + τ) for u in H¹‐norm and p = ?u in L²‐norm are derived respectively without the restrictions on the ratio between h and τ, where h is the subdivision parameter and τ, the time step. Third, the global superconvergent results are obtained by interpolated postprocessing technique. Finally, some numerical results are carried out to confirm the theoretical analysis.     

Superconvergence of the split least‐squares method for second‐order hyperbolic equations

This article studies superconvergence phenomena of the split least‐squares mixed finite element method for second‐order hyperbolic equations. By selecting the least‐squares functional properly, the procedure can be split into two independent symmetric positive definite subprocedures, one of which is for the primitive unknown and the other is for the flux. Based on interpolation operators and an auxiliary projection, superconvergent H¹ error estimates for the primary variable u and L² error estimates for the introduced flux variable σ are obtained under the standard quasiuniform assumptions on finite element partition. A numerical example is given to show the performance of the introduced scheme. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 222‐238, 2014     

Monte Carlo approximate tensor moment simulations

An algorithm to generate samples with approximate first‐order, second‐order, third‐order, and fourth‐order moments is presented by extending the Cholesky matrix decomposition to a Cholesky tensor decomposition of an arbitrary order. The tensor decomposition of the first‐order, second‐order, third‐order, and fourth‐order objective moments generates a non‐linear system of equations. The algorithm solves these equations by numerical methods. The results show that the optimization algorithm delivers samples with an approximate residual error of less than 10¹⁶ between the components of the objective and the sample moments. The algorithm is extended for a n‐th‐order approximate tensor moment version, and simulations of non‐normal samples replicated from distributions with asymmetries and heavy tails are presented. An application for sensitivity analysis of portfolio risk assessment with Value‐at‐Risk (VaR) is provided. A comparison with previous methods available in the literature suggests that the methodology proposed reduces the error of the objective moments in the generated samples. ? ? JEL Classification: C14, C15, G32.
Copyright © 2016 John Wiley & Sons, Ltd.     

Convergence analysis of an approximation to miscible fluid flows in porous media by combining mixed finite element and finite volume methods

This article deals with development and analysis of a numerical method for a coupled system describing miscible displacement of one incompressible fluid by another through heterogeneous porous media. A mixed finite element (MFE) method is employed to discretize the Darcy flow equation combined with a conservative finite volume (FV) method on unstructured grids for the concentration equation. It is shown that the FV scheme satisfies a discrete maximum principle. We derive L^∞ and BV estimates under an appropriate CFL condition. Then we prove convergence of the approximate solutions to a weak solution of the coupled system. Numerical results are presented to see the performance of the method in two space dimensions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Analysis of a combined mixed finite element and discontinuous Galerkin method for incompressible two‐phase flow in porous media

We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two‐phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L²(H¹) for saturation and in L ^∞ (H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.     

Two new least-squares mixed finite element procedures for convection-dominated Sobolev equations

Two new least-squares mixed finite element procedures are formulated for solving convection-dominated Sobolev equations. Optimal H(div;Ω)×H ¹(Ω) norms error estimates are derived under the standard mixed finite spaces. Moreover, these two schemes provide the approximate solutions with first-order and second-order accuracy in time increment, respectively.