TY - JOUR

T1 - An algorithm for stabilizing hybridizable discontinuous Galerkin methods for nonlinear elasticity

AU - Cockburn, Bernardo

AU - Shen, Jiguang

N1 - Funding Information:
Partially supported by NSF grant DMS 1522657.
Funding Information:
Partially supported by NSF grant DMS 1522657 .
Publisher Copyright:
© 2019 The Author(s)

PY - 2019/6

Y1 - 2019/6

N2 - It is now a well known fact that hybridizable discontinuous Galerkin for nonlinear elasticity may not converge to the exact solution if their inter-element jumps are not properly penalized or, equivalently, if the values of their stabilization function are not suitably chosen. For example, if their stabilization function is of order one, the method generates spurious oscillations in the region in which the elastic moduli tensor is indefinite. This phenomenon disappears when the values of the stabilization function are increased. Mixed methods display the same problem, as their stabilization function is identically zero. Here, we find explicit formulas for the lower bound of the values of the stabilization function which allow us to avoid this unpleasant phenomenon. Our numerical experiments show that, when we use polynomials of degree k>0 for the approximate gradient, first-order Piola-Kirchhoff and displacement, all the approximations converge with the optimal order, k+1. They also show that, if we increase the polynomial degree of the local spaces by one and update the values of the stabilization function only for the elements for which the elasticity tensor is indefinite, we obtain that all approximate solutions converge with order k+1 and that a local post processing of the displacement converges with order k+2.

AB - It is now a well known fact that hybridizable discontinuous Galerkin for nonlinear elasticity may not converge to the exact solution if their inter-element jumps are not properly penalized or, equivalently, if the values of their stabilization function are not suitably chosen. For example, if their stabilization function is of order one, the method generates spurious oscillations in the region in which the elastic moduli tensor is indefinite. This phenomenon disappears when the values of the stabilization function are increased. Mixed methods display the same problem, as their stabilization function is identically zero. Here, we find explicit formulas for the lower bound of the values of the stabilization function which allow us to avoid this unpleasant phenomenon. Our numerical experiments show that, when we use polynomials of degree k>0 for the approximate gradient, first-order Piola-Kirchhoff and displacement, all the approximations converge with the optimal order, k+1. They also show that, if we increase the polynomial degree of the local spaces by one and update the values of the stabilization function only for the elements for which the elasticity tensor is indefinite, we obtain that all approximate solutions converge with order k+1 and that a local post processing of the displacement converges with order k+2.

UR - http://www.scopus.com/inward/record.url?scp=85065210513&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85065210513&partnerID=8YFLogxK

U2 - 10.1016/j.rinam.2019.01.001

DO - 10.1016/j.rinam.2019.01.001

M3 - Article

AN - SCOPUS:85065210513

VL - 1

JO - Results in Applied Mathematics

JF - Results in Applied Mathematics

SN - 2590-0374

M1 - 100001

ER -