A shock-capturing discontinuous Galerkin method for the numerical solution of inviscid compressible flow
A Sobolev gradient method for treating the steady-state incompressible Navier-Stokes equations
The velocity-vorticity-pressure formulation of the steady-state incompressible Navier-Stokes equations in two dimensions is cast as a nonlinear least squares problem in which the functional is a weighted sum of squared residuals. A finite element discretization of the functional is minimized by a trust-region method in which the trustregion radius is defined by a Sobolev norm and the trust-region subproblems are solved by a dogleg method. Numerical test results show the method to be effective.
A special finite element method based on component mode synthesis
The goal of our paper is to introduce basis functions for the finite element discretization of a second order linear elliptic operator with rough or highly oscillating coefficients. The proposed basis functions are inspired by the classic idea of component mode synthesis and exploit an orthogonal decomposition of the trial subspace to minimize the energy. Numerical experiments illustrate the effectiveness of the proposed basis functions.
A stabilized finite element scheme for the Navier-Stokes equations on quadrilateral anisotropic meshes
It is well known that the classical local projection method as well as residual-based stabilization techniques, as for instance streamline upwind Petrov-Galerkin (SUPG), are optimal on isotropic meshes. Here we extend the local projection stabilization for the Navier-Stokes system to anisotropic quadrilateral meshes in two spatial dimensions. We describe the new method and prove an a priori error estimate. This method leads on anisotropic meshes to qualitatively better convergence behavior...
A Stabilized Lagrange Multiplier Method for the Finite Element Approximation of Frictional Contact Problems in Elastostatics
In this work we consider a stabilized Lagrange multiplier method in order to approximate the Coulomb frictional contact model in linear elastostatics. The particularity of the method is that no discrete inf-sup condition is needed. We study the existence and the uniqueness of solution of the discrete problem.
A stable and optimal complexity solution method for mixed finite element discretizations
We outline a solution method for mixed finite element discretizations based on dissecting the problem into three separate steps. The first handles the inhomogeneous constraint, the second solves the flux variable from the homogeneous problem, whereas the third step, adjoint to the first, finally gives the Lagrangian multiplier. We concentrate on aspects involved in the first and third step mainly, and advertise a multi-level method that allows for a stable computation of the intermediate and final...
A stable and optimal complexity solution method for mixed finite element discretizations
A stable mixed finite element method on truncated exterior domains
A Static condensation Reduced Basis Element method : approximation and a posteriori error estimation
We propose a new reduced basis element-cum-component mode synthesis approach for parametrized elliptic coercive partial differential equations. In the Offline stage we construct a Library of interoperable parametrized reference components relevant to some family of problems; in the Online stage we instantiate and connect reference components (at ports) to rapidly form and query parametric systems. The method is based on static condensation at the interdomain level, a conforming eigenfunction “port”...
A strongly nonlinear problem arising in glaciology
A strongly nonlinear problem arising in glaciology
The computation of glacier movements leads to a system of nonlinear partial differential equations. The existence and uniqueness of a weak solution is established by using the calculus of variations. A discretization by the finite element method is done. The solution of the discrete problem is proved to be convergent to the exact solution. A first simple numerical algorithm is proposed and its convergence numerically studied.
A subspace correction method for discontinuous Galerkin discretizations of linear elasticity equations
We study preconditioning techniques for discontinuous Galerkin discretizations of isotropic linear elasticity problems in primal (displacement) formulation. We propose subspace correction methods based on a splitting of the vector valued piecewise linear discontinuous finite element space, that are optimal with respect to the mesh size and the Lamé parameters. The pure displacement, the mixed and the traction free problems are discussed in detail. We present a convergence analysis of the proposed...
A Superconvergence result for mixed finite element approximations of the eigenvalue problem
In this paper, we present a superconvergence result for the mixed finite element approximations of general second order elliptic eigenvalue problems. It is known that a superconvergence result has been given by Durán et al. [Math. Models Methods Appl. Sci. 9 (1999) 1165–1178] and Gardini [ESAIM: M2AN 43 (2009) 853–865] for the lowest order Raviart-Thomas approximation of Laplace eigenvalue problems. In this work, we introduce a new way to derive the superconvergence of general second order elliptic...
A Superconvergence result for mixed finite element approximations of the eigenvalue problem∗
In this paper, we present a superconvergence result for the mixed finite element approximations of general second order elliptic eigenvalue problems. It is known that a superconvergence result has been given by Durán et al. [Math. Models Methods Appl. Sci. 9 (1999) 1165–1178] and Gardini [ESAIM: M2AN 43 (2009) 853–865] for the lowest order Raviart-Thomas approximation of Laplace eigenvalue problems. In this work, we introduce a new way to derive the superconvergence of general second order elliptic...
A technique of upstream type applied to a linear nonconforming finite element approximation of convective diffusion equations
A tensor approximation method based on ideal minimal residual formulations for the solution of high-dimensional problems
In this paper, we propose a method for the approximation of the solution of high-dimensional weakly coercive problems formulated in tensor spaces using low-rank approximation formats. The method can be seen as a perturbation of a minimal residual method with a measure of the residual corresponding to the error in a specified solution norm. The residual norm can be designed such that the resulting low-rank approximations are optimal with respect to particular norms of interest, thus allowing to take...
A Three-Dimensional Quadratic Nonconforming Element.
A three-field augmented Lagrangian formulation of unilateral contact problems with cohesive forces
We investigate unilateral contact problems with cohesive forces, leading to the constrained minimization of a possibly nonconvex functional. We analyze the mathematical structure of the minimization problem. The problem is reformulated in terms of a three-field augmented Lagrangian, and sufficient conditions for the existence of a local saddle-point are derived. Then, we derive and analyze mixed finite element approximations to the stationarity conditions of the three-field augmented Lagrangian....
A triangular mixed finite element method for the stationary semiconductor device equations
A two parameter iterative method for solving algebraic systems of domain decomposition type
An iterative procedure containing two parameters for linear algebraic systems originating from the domain decomposition technique is proposed. The optimization of the parameters is investigated. A numeric example is given as an illustration.