A Locking-Free hp DPG Method for Linear Elasticity with Symmetric Stresses

We present two new methods for linear elasticity that simultaneously yield stress and displacement approximations of optimal accuracy in both the mesh size h and polynomial degree p. This is achieved within the recently developed discontinuous Petrov- Galerkin (DPG) framework. In this framework, both the stress and the displacement ap- proximations are discontinuous across element interfaces. We study locking-free convergence properties and the interrelationships between the two DPG methods

Computational Engineering

The focus of this Computational Engineering Workshop was on the mathematical foundation of state-of-the-art and emerging finite element methods in engineering analysis. The 52 participants included mathematicians and engineers with shared interest on discontinuous Galerkin or Petrov-Galerkin methods and other generalized nonconforming or mixed finite element methods

An analysis of the practical DPG method

In this work we give a complete error analysis of the Discontinuous Petrov Galerkin (DPG) method, accounting for all the approximations made in its practical implementation. Specifically, we consider the DPG method that uses a trial space consisting of polynomials of degree $p$ on each mesh element. Earlier works showed that there is a "trial-to-test" operator $T$, which when applied to the trial space, defines a test space that guarantees stability. In DPG formulations, this operator $T$ is local: it can be applied element-by-element. However, an infinite dimensional problem on each mesh element needed to be solved to apply $T$. In practical computations, $T$ is approximated using polynomials of some degree $r > p$ on each mesh element. We show that this approximation maintains optimal convergence rates, provided that $r\ge p+N$, where $N$ is the space dimension (two or more), for the Laplace equation. We also prove a similar result for the DPG method for linear elasticity. Remarks on the conditioning of the stiffness matrix in DPG methods are also included.Comment: Mathematics of Computation, 201

Various applications of discontinuous Petrov-Galerkin (DPG) finite element methods

Discontinuous Petrov-Galerkin (DPG) finite element methods have garnered significant attention since they were originally introduced. They discretize variational formulations with broken (discontinuous) test spaces and are crafted to be numerically stable by implicitly computing a near-optimal discrete test space as a function of a discrete trial space. Moreover, they are completely general in the sense that they can be applied to a variety of variational formulations, including non-conventional ones that involve non-symmetric functional settings, such as ultraweak variational formulations. In most cases, these properties have been harnessed to develop numerical methods that provide robust control of relevant equation parameters, like in convection-diffusion problems and other singularly perturbed problems. In this work, other features of DPG methods are systematically exploited and applied to different problems. More specifically, the versatility of DPG methods is elucidated by utilizing the underlying methodology to discretize four distinct variational formulations of the equations of linear elasticity. By taking advantage of interface variables inherent to DPG discretizations, an approach to coupling different variational formulations within the same domain is described and used to solve interesting problems. Moreover, the convenient algebraic structure in DPG methods is harnessed to develop a new family of numerical methods called discrete least-squares (DLS) finite element methods. These involve solving, with improved conditioning properties, a discrete least-squares problem associated with an overdetermined rectangular system of equations, instead of directly solving the usual square systems. Their utility is demonstrated with illustrative examples. Additionally, high-order polygonal DPG (PolyDPG) methods are devised by using the intrinsic discontinuities present in ultraweak formulations. The resulting methods can handle heavily distorted non-convex polygonal elements and discontinuous material properties. A polygonal adaptive strategy was also proposed and compared with standard techniques. Lastly, the natural high-order residual-based a posteriori error estimator ingrained within DPG methods was further applied to problems of physical relevance, like the validation of dynamic mechanical analysis (DMA) calibration experiments of viscoelastic materials, and the modeling of form-wound medium-voltage stator coils sitting inside large electric machinery.Computational Science, Engineering, and Mathematic

An Analysis of the Practical DPG Method

We give a complete error analysis of the Discontinuous Petrov Galerkin (DPG) method, accounting for all the approximations made in its practical implementation. Specifically, we consider the DPG method that uses a trial space consisting of polynomials of degree p on each mesh element. Earlier works showed that there is a trial-to-test operator T, which when applied to the trial space, defines a test space that guarantees stability. In DPG formulations, this operator T is local: it can be applied element-by-element. However, an infinite dimensional problem on each mesh element needed to be solved to apply T. In practical computations, T is approximated using polynomials of some degree r \u3e p on each mesh element. We show that this approximation maintains optimal convergence rates, provided that r p + N, where N is the space dimension (two or more), for the Laplace equation. We also prove a similar result for the DPG method for linear elasticity. Remarks on the conditioning of the stiffness matrix in DPG methods are also included

A superconvergent hybridisable discontinuous Galerkin method for linear elasticity

The first superconvergent hybridisable discontinuous Galerkin method for linear elastic problems capable of using the same degree of approximation for both the primal and mixed variables is presented. The key feature of the method is the strong imposition of the symmetry of the stress tensor by means of the well known and extensively used Voigt notation, circumventing the use of complex mathematical concepts to enforce the symmetry of the stress tensor either weakly or strongly. A novel procedure to construct element by element a superconvergent postprocessed displacement is proposed. Contrary to other hybridisable discontinuous Galerkin formulations, the methodology proposed here is able to produce a superconvergent displacement field for low-order approximations. The resulting method is robust and locking-free in the nearly incompressible limit. An extensive set of numerical examples is utilised to provide evidence of the optimality of the method and its superconvergent properties in two and three dimensions and for different element type