Imperial College London
Dr Matthew Eaton

DrMatthewEaton

Faculty of Engineering, Department of Mechanical Engineering

Reader
 
 
 
//

Contact

 

+44 (0)20 7594 7053m.eaton

 
 
//

Location

 

657City and Guilds BuildingSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Owens:2017:10.1016/j.jcp.2017.09.020,
author = {Owens, A and Kophazi, J and Eaton, M},
doi = {10.1016/j.jcp.2017.09.020},
journal = {Journal of Computational Physics},
pages = {847--870},
title = {Optimal trace inequality constants for interior penalty discontinuous Galerkin discretisations of elliptic operators using arbitrary elements with non-constant Jacobians},
url = {http://dx.doi.org/10.1016/j.jcp.2017.09.020},
volume = {350},
year = {2017}
}
Download

RIS format (EndNote, RefMan)

TY  - JOUR
AB  - In this paper, a new method to numerically calculate the trace inequality constants, which arise in the calculation of penalty parameters for interior penalty discretisations of elliptic operators, is presented. These constants are provably optimal for the inequality of interest. As their calculation is based on the solution of a generalised eigenvalue problem involving the volumetric and face stiffness matrices, the method is applicable to any element type for which these matrices can be calculated, including standard finite elements and the non-uniform rational B-splines of isogeometric analysis. In particular, the presented method does not require the Jacobian of the element to be constant, and so can be applied to a much wider variety of element shapes than are currently available in the literature. Numerical results are presented for a variety of finite element and isogeometric cases. When the Jacobian is constant, it is demonstrated that the new method produces lower penalty parameters than existing methods in the literature in all cases, which translates directly into savings in the solution time of the resulting linear system. When the Jacobian is not constant, it is shown that the naive application of existing approaches can result in penalty parameters that do not guarantee coercivity of the bilinear form, and by extension, the stability of the solution. The method of manufactured solutions is applied to a model reaction-diffusion equation with a range of parameters, and it is found that using penalty parameters based on the new trace inequality constants result in better conditioned linear systems, which can be solved approximately 11% faster than those produced by the methods from the literature.
AU  - Owens,A
AU  - Kophazi,J
AU  - Eaton,M
DO  - 10.1016/j.jcp.2017.09.020
EP  - 870
PY  - 2017///
SN  - 0021-9991
SP  - 847
TI  - Optimal trace inequality constants for interior penalty discontinuous Galerkin discretisations of elliptic operators using arbitrary elements with non-constant Jacobians
T2  - Journal of Computational Physics
UR  - http://dx.doi.org/10.1016/j.jcp.2017.09.020
UR  - https://www.sciencedirect.com/science/article/pii/S0021999117306794
UR  - http://hdl.handle.net/10044/1/50698
VL  - 350
ER  -
Download