On the Ma–Trudinger–Wang curvature on surfaces

We investigate the properties of the Ma–Trudinger–Wang nonlocal curvature tensor in the case of surfaces. In particular, we prove that a strict form of the Ma–Trudinger– Wang condition is stable under C ⁴ perturbation if the nonfocal domains are uniformly convex; and we present new examples of positively curved surfaces which do not satisfy the Ma–Trudinger–Wang condition. As a corollary of our results, optimal transport maps on a “sufficiently flat” ellipsoid are in general nonsmooth.     

Local and pointwise error estimates of the local discontinuous Galerkin method applied to the Stokes problem

We prove local and pointwise error estimates for the local discontinuous Galerkin method applied to the Stokes problem in two and three dimensions. By using techniques originally developed by A. Schatz [Math. Comp., 67 (1998), 877-899] to prove pointwise estimates for the Laplace equation, we prove optimal weighted pointwise estimates for both the velocity and the pressure for domains with smooth boundaries.

     

Finite element approximations to surfaces of prescribed variable mean curvature

We give an algorithm for finding finite element approximations to surfaces of prescribed variable mean curvature, which span a given boundary curve. We work in the parametric setting and prove optimal estimates in the H¹ norm. The estimates are verified computationally.     

On a mixed finite element formulation of a second‐order quasilinear problem in the plane

We analyze a mixed finite element discretization of a second‐order quasilinear problem based on the Raviart‐Thomas space. We prove that the discrete problem is solvable and provide a local uniqueness result for the solution. We also obtain optimal order L²‐error estimates for both the scalar variable and the associated flux. The main feature of our method is that it is free from the boundness conditions required in previous works on the coefficients of the quasilinear operator. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 90–103, 2004.     

Finite element convergence for the Joule heating problem with mixed boundary conditions

We prove strong convergence of conforming finite element approximations to the stationary Joule heating problem with mixed boundary conditions on Lipschitz domains in three spatial dimensions. We show optimal global regularity estimates on creased domains and prove a priori and a posteriori bounds for shape regular meshes.     

A posteriori error estimate for discontinuous Galerkin finite element method on polytopal mesh

In this work, we derive a posteriori error estimates for discontinuous Galerkin finite element method on polytopal mesh. We construct a reliable and efficient a posteriori error estimator on general polygonal or polyhedral meshes. An adaptive algorithm based on the error estimator and DG method is proposed to solve a variety of test problems. Numerical experiments are performed to illustrate the effectiveness of the algorithm.     

Refined mixed finite element method for the elasticity problem in a polygonal domain

The purpose of this article is to study a mixed formulation of the elasticity problem in plane polygonal domains and its numerical approximation. In this mixed formulation the strain tensor is introduced as a new unknown and its symmetry is relaxed by a Lagrange multiplier, which is nothing else than the rotation. Because of the corner points, the displacement field is not regular in general in the vicinity of the vertices but belongs to some weighted Sobolev space. Using this information, appropriate refinement rules are imposed on the family of triangulations in order to recapture optimal error estimates. Moreover, uniform error estimates in the Lamé coefficient λ are obtained for λ large. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 323–339, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10009     

Finite element approximation of the elasticity spectral problem on curved domains

We analyze the finite element approximation of the spectral problem for the linear elasticity equation with mixed boundary conditions on a curved non-convex domain. In the framework of the abstract spectral approximation theory, we obtain optimal order error estimates for the approximation of eigenvalues and eigenvectors. Two kinds of problems are considered: the discrete domain does not coincide with the real one and mixed boundary conditions are imposed. Some numerical results are presented.     

Error estimates for finite volume element methods for general second‐order elliptic problems

We treat the finite volume element method (FVE) for solving general second order elliptic problems as a perturbation of the linear finite element method (FEM), and obtain the optimal H¹ error estimate, H¹ superconvergence and L^p (1 < p ≤ ∞) error estimates between the solution of the FVE and that of the FEM. In particular, the superconvergence result does not require any extra assumptions on the mesh except quasi‐uniform. Thus the error estimates of the FVE can be derived by the standard error estimates of the FEM. Moreover we consider the effects of numerical integration and prove that the use of barycenter quadrature rule does not decrease the convergence orders of the FVE. The results of this article reveal that the FVE is in close relationship with the FEM. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 693–708, 2003.     

Convergence of the mixed finite element method for Maxwell's equations with nonlinear conductivity

In this paper, we study a numerical scheme to solve coupled Maxwell's equations with a nonlinear conductivity. This model plays an important role in the study of type‐II superconductors. The approximation scheme is based on backward Euler discretization in time and mixed conforming finite elements in space. We will prove convergence of this scheme to the unique weak solution of the problem and develop the corresponding error estimates. As a next step, we study the stability of the scheme in the quasi‐static limit ? → 0 and present the corresponding convergence rate. Finally, we support the theory by several numerical experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

Local a posteriori estimates for pointwise gradient errors in finite element methods for elliptic problems

We prove local a posteriori error estimates for pointwise gradient errors in finite element methods for a second-order linear elliptic model problem. First we split the local gradient error into a computable local residual term and a weaker global norm of the finite element error (the ``pollution term'). Using a mesh-dependent weight, the residual term is bounded in a sharply localized fashion. In specific situations the pollution term may also be bounded by computable residual estimators. On nonconvex polygonal and polyhedral domains in two and three space dimensions, we may choose estimators for the pollution term which do not employ specific knowledge of corner singularities and which are valid on domains with cracks. The finite element mesh is only required to be simplicial and shape-regular, so that highly graded and unstructured meshes are allowed.

     

On a finite element method for three‐dimensional unsteady compressible viscous flows

In this article we analyze a finite element method for three‐dimensional unsteady compressible Navier‐Stokes equations. We prove the existence and uniqueness of the numerical solution, and obtain a priori error estimates uniform in time. Numerical computations are carried out to test the orders of accuracy in the error estimates. Blend function interpolations are applied in the calculation of numerical integrations. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 432–449, 2004.     

Multiphysics mixed finite element method with Nitsche's technique for Stokes-poroelasticity problem

In this article, we propose a multiphysics mixed finite element method with Nitsche's technique for Stokes-poroelasticity problem. Firstly, we reformulate the poroelasticity part of the original problem by introducing two pseudo-pressures to into a “fluid–fluid” coupled problem so that we can use the classical stable finite element pairs to deal with this problem conveniently. Then, we prove the existence and uniqueness of weak solution of the reformulated problem. And we use Nitsche's technique to approximate the coupling condition at the interface to propose a loosely-coupled time-stepping method to solve three subproblems at each time step–a Stokes problem, a generalized Stokes problem and a mixed diffusion problem. And the proposed method does not require any restriction on the choice of the discrete approximation spaces on each side of the interface provided that appropriate quadrature methods are adopted. Also, we give the stability analysis and error estimates of the loosely-coupled time-stepping method. Finally, we give the numerical tests to show that the proposed numerical method has a good stability and no “locking” phenomenon.     

A hybrid high-order formulation for a Neumann problem on polytopal meshes

In this work, we study a hybrid high-order (HHO) method for an elliptic diffusion problem with Neumann boundary condition. The proposed method has several features, such as: (a) the support of arbitrary approximation order polynomial at mesh elements and faces on polytopal meshes, (b) the design of a local (element-wise) potential reconstruction operator and a local stabilization term, that weakly enforces the matching between local element- and face-based on degrees of freedom, and (c) cheap computational cost, thanks to static condensation and compact stencil. We prove the well-posedness of our HHO formulation, and obtain the optimal error estimates, according to previous study. Implementation aspects are thoroughly discussed. Finally, some numerical examples are provided, which are in agreement with our theoretical results.     

Stability estimates for h-p spectral element methods for elliptic problems

In a series of papers of which this is the first we study how to solve elliptic problems on polygonal domains using spectral methods on parallel computers. To overcome the singularities that arise in a neighborhood of the corners we use a geometrical mesh. With this mesh we seek a solution which minimizes a weighted squared norm of the residuals in the partial differential equation and a fractional Sobolev norm of the residuals in the boundary conditions and enforce continuity by adding a term which measures the jump in the function and its derivatives at inter-element boundaries, in an appropriate fractional Sobolev norm, to the functional being minimized. Since the second derivatives of the actual solution are not square integrable in a neighborhood of the corners we have to multiply the residuals in the partial differential equation by an appropriate power of r_k, where r_k measures the distance between the pointP and the vertexA _k in a sectoral neighborhood of each of these vertices. In each of these sectoral neighborhoods we use a local coordinate system (τ_k, θ_k) where τ_k = lnr_k and (r_k, θ_k) are polar coordinates with origin at A_k, as first proposed by Kondratiev. We then derive differentiability estimates with respect to these new variables and a stability estimate for the functional we minimize. In [6] we will show that we can use the stability estimate to obtain parallel preconditioners and error estimates for the solution of the minimization problem which are nearly optimal as the condition number of the preconditioned system is polylogarithmic inN, the number of processors and the number of degrees of freedom in each variable on each element. Moreover if the data is analytic then the error is exponentially small inN.     

Finite element approximation of a semilinear elliptic problem with a singular nonlinearity

Given , we consider the following problem: find , such that where or 3, and in . We prove and error bounds for the standard continuous piecewise linear Galerkin finite element approximation with a (weakly) acute triangulation. Our bounds are nearly optimal. In addition, for d = 1 and 2 and we analyze a more practical scheme involving numerical integration on the nonlinear term. We obtain nearly optimal and error bounds for d = 1. For this case we also present some numerical results. Received July 4, 1996 / Revised version received December 18, 1997     

Error analysis of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation

This article discusses a priori and a posteriori error estimates of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation. We use variational discretization concept to discretize the control variable and discontinuous piecewise linear finite elements to approximate the state and costate variable. Based on the error estimates of discontinuous Galerkin finite element method for the transport equation, we get a priori and a posteriori error estimates for the transport equation optimal control problem. Finally, two numerical experiments are carried out to confirm the theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1493–1512, 2017     

The heterogeneous multiscale finite element method for elliptic homogenization problems in perforated domains

In this contribution we analyze a generalization of the heterogeneous multiscale finite element method for elliptic homogenization problems in perforated domains. The method was originally introduced by E and Engquist (Commun Math Sci 1(1):87–132, 2003) for homogenization problems in fixed domains. It is based on a standard finite element approach on the macroscale, where the stiffness matrix is computed by solving local cell problems on the microscale. A-posteriori error estimates are derived in L ²(Ω) by reformulating the problem into a discrete two-scale formulation (see also, Ohlberger in Multiscale Model Simul 4(1):88–114, 2005) and using duality methods afterwards. Numerical experiments are given in order to numerically evaluate the efficiency of the error estimate.     

Finite element analysis of the vibration problem of a plate coupled with a fluid

Summary. We consider the approximation of the vibration modes of an elastic plate in contact with a compressible fluid. The plate is modelled by Reissner-Mindlin equations while the fluid is described in terms of displacement variables. This formulation leads to a symmetric eigenvalue problem. Reissner-Mindlin equations are discretized by a mixed method, the equations for the fluid with Raviart-Thomas elements and a non conforming coupling is used on the interface. In order to prove that the method is locking free we consider a family of problems, one for each thickness , and introduce appropriate scalings for the physical parameters so that these problems attain a limit when . We prove that spurious eigenvalues do not arise with this discretization and we obtain optimal order error estimates for the eigenvalues and eigenvectors valid uniformly on the thickness parameter t. Finally we present numerical results confirming the good performance of the method. Received February 4, 1998 / Revised version received May 26, 1999 / Published online June 21, 2000     

Local and parallel finite element method for solving the biharmonic eigenvalue problem of plate vibration

In this paper, we establish a new local and parallel finite element discrete scheme based on the shifted‐inverse power method for solving the biharmonic eigenvalue problem of plate vibration. We prove the local error estimation of finite element solution for the biharmonic equation/eigenvalue problem and prove the error estimation of approximate solution obtained by the local and parallel scheme. When the diameters of three grids satisfy H⁴ = ?(w²) = ?(h), the approximate solutions obtained by our schemes can achieve the asymptotically optimal accuracy. The numerical experiments show that the computational schemes proposed in this paper are effective to solve the biharmonic eigenvalue problem of plate vibration.