hp-version interior penalty DGFEMs for the biharmonic equation

We construct hp-version interior penalty discontinuous Galerkin finite element methods (DGFEMs) for the biharmonic equation, including symmetric and nonsymmetric interior penalty discontinuous Galerkin methods and their combinations: semisymmetric methods. Our main concern is to establish the stability and to develop the a priori error analysis of these methods. We establish error bounds that are optimal in h and slightly suboptimal in p. The theoretical results are confirmed by numerical experiments

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

International audienceWe design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters

A C 0 finite element method for the biharmonic problem without extrinsic penalization

A symmetric C0 finite element method for the biharmonic problem is constructed and analyzed. In our approach, we introduce one-sided discrete second-order derivatives and Hessian matrices to formulate our scheme. We show that the method is stable and converge with optimal order in a variety of norms. A distinctive feature of the method is that the results hold without extrinsic penalization of the gradient across interelement boundaries. Numerical experiments are given that support the theoretical results, and the extension to Kirchhoff plates is also discusse

Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment

This is the peer reviewed version of the following article: [García-Blanco, R., Borzacchiello, D., Chinesta, F., and Diez, P. (2017) Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment. Int. J. Numer. Meth. Engng, 111: 529–552. doi: 10.1002/nme.5470], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5470/full. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The parametric analysis of electric grids requires carrying out a large number of Power Flow computations. The different parameters describe loading conditions and grid properties. In this framework, the Proper Generalized Decomposition (PGD) provides a numerical solution explicitly accounting for the parametric dependence. Once the PGD solution is available, exploring the multidimensional parametric space is computationally inexpensive. The aim of this paper is to provide tools to monitor the error associated with this significant computational gain and to guarantee the quality of the PGD solution. In this case, the PGD algorithm consists in three nested loops that correspond to 1) iterating algebraic solver, 2) number of terms in the separable greedy expansion and 3) the alternated directions for each term. In the proposed approach, the three loops are controlled by stopping criteria based on residual goal-oriented error estimates. This allows one for using only the computational resources necessary to achieve the accuracy prescribed by the end- user. The paper discusses how to compute the goal-oriented error estimates. This requires linearizing the error equation and the Quantity of Interest to derive an efficient error representation based on an adjoint problem. The efficiency of the proposed approach is demonstrated on benchmark problems.Peer ReviewedPostprint (author's final draft

hp-version interior penalty DGFEMs for the biharmonic equation

We construct hp-version interior penalty discontinuous Galerkin finite element methods (DGFEMs) for the biharmonic equation, including symmetric and nonsymmetric interior penalty discontinuous Galerkin methods and their combinations: semisymmetric methods. Our main concern is to establish the stability and to develop the a priori error analysis of these methods. We establish error bounds that are optimal in h and slightly suboptimal in p. The theoretical results are confirmed by numerical experiments. © 2006 Elsevier B.V. All rights reserved

Goal-oriented error estimation based on equilibrated-flux reconstruction for finite element approximations of elliptic problems

We propose an approach for goal-oriented error estimation in finite element approximations of second-order elliptic problems that combines the dual-weighted residual method and equilibrated-flux reconstruction methods for the primal and dual problems. The objective is to be able to consider discretization schemes for the dual solution that may be different from those used for the primal solution. It is only assumed here that the discretization methods come with a priori error estimates and an equilibrated-flux reconstruction algorithm. A high-order discontinuous Galerkin (dG) method is actually the preferred choice for the approximation of the dual solution thanks to its flexibility and straightforward construction of equilibrated fluxes. One contribution of the paper is to show how the order of the dG method for asymptotic exactness of the proposed estimator can be chosen in the cases where a conforming finite element method, a dG method, or a mixed Raviart- Thomas method are used for the solution of the primal problem. Numerical experiments are also presented to illustrate the performance and convergence of the error estimates in quantities of interest with respect to the mesh size

A robust goal-oriented estimator based on the construction of equilibrated fluxes for discontinuous Galerkin finite element approximations of convection-diffusion problems

We propose an a posteriori error estimator with respect to quantities of interest for dis-continuous Galerkin approximations of convection-diffusion boundary-value problems. The error estimator is based on the construction of equilibrated fluxes in Raviart-Thomas finite element spaces and on the solution of the dual problem. We show that it is asymptotically exact in both the elliptic and hyperbolic regimes if the dual problem is approximated by a discontinuous Galerkin method of order one greater than that of the primal problem. We show in this case that the effectivity index behaves as (1+Pe 1/2)o(h), where Pe is the Péclet number and h the mesh diameter. It follows that the quality of the effectivity index may deteriorate for large values of Pe, but we put in evidence that it suffices to increase the approximation order of the dual problem to keep the effectivity index close to unity even on coarse meshes. Two-dimensional numerical examples demonstrate the robustness of the error estimator in both the diffusion and advection regimes

hp-version interior penalty DGFEMs for the biharmonic equation

We construct hp-version interior penalty discontinuous Galerkin finite element methods (DGFEMs) for the biharmonic equation, including symmetric and nonsymmetric interior penalty discontinuous Galerkin methods and their combinations: semisymmetric methods. Our main concern is to establish the stability and to develop the a priori error analysis of these methods. We establish error bounds that are optimal in h and slightly suboptimal in p. The theoretical results are confirmed by numerical experiments. © 2006 Elsevier B.V. All rights reserved