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.