Monolithic multiphysics simulation of hypersonic aerothermoelasticity using a hybridized discontinuous Galerkin method

This work presents implementation of a hybridized discontinuous Galerkin (DG) method for robust simulation of the hypersonic aerothermoelastic multiphysics system. Simulation of hypersonic vehicles requires accurate resolution of complex multiphysics interactions including the effects of high-speed turbulent flow, extreme heating, and vehicle deformation due to considerable pressure loads and thermal stresses. However, the state-of-the-art procedures for hypersonic aerothermoelasticity are comprised of low-fidelity approaches and partitioned coupling schemes. These approaches preclude robust design and analysis of hypersonic vehicles for a number of reasons. First, low-fidelity approaches limit their application to simple geometries and lack the ability to capture small scale flow features (e.g. turbulence, shocks, and boundary layers) which greatly degrades modeling robustness and solution accuracy. Second, partitioned coupling approaches can introduce considerable temporal and spatial inaccuracies which are not trivially remedied. In light of these barriers, we propose development of a monolithically-coupled hybridized DG approach to enable robust design and analysis of hypersonic vehicles with arbitrary geometries. Monolithic coupling methods implement a coupled multiphysics system as a single, or monolithic, equation system to be resolved by a single simulation approach. Further, monolithic approaches are free from the physical inaccuracies and instabilities imposed by partitioned approaches and enable time-accurate evolution of the coupled physics system. In this work, a DG method is considered due to its ability to accurately resolve second-order partial differential equations (PDEs) of all classes. We note that the hypersonic aerothermoelastic system is composed of PDEs of all three classes. Hybridized DG methods are specifically considered due to their exceptional computational efficiency compared to traditional DG methods. It is expected that our monolithic hybridized DG implementation of the hypersonic aerothermoelastic system will 1) provide the physical accuracy necessary to capture complex physical features, 2) be free from any spatial and temporal inaccuracies or instabilities inherent to partitioned coupling procedures, 3) represent a transition to high-fidelity simulation methods for hypersonic aerothermoelasticity, and 4) enable efficient analysis of hypersonic aerothermoelastic effects on arbitrary geometries

An advection-robust Hybrid High-Order method for the Oseen problem

In this work, we study advection-robust Hybrid High-Order discretizations of the Oseen equations. For a given integer $k\ge 0$, the discrete velocity unknowns are vector-valued polynomials of total degree $\le k$ on mesh elements and faces, while the pressure unknowns are discontinuous polynomials of total degree $\le k$ on the mesh. From the discrete unknowns, three relevant quantities are reconstructed inside each element: a velocity of total degree $\le(k+1)$, a discrete advective derivative, and a discrete divergence. These reconstructions are used to formulate the discretizations of the viscous, advective, and velocity-pressure coupling terms, respectively. Well-posedness is ensured through appropriate high-order stabilization terms. We prove energy error estimates that are advection-robust for the velocity, and show that each mesh element $T$ of diameter $h_T$ contributes to the discretization error with an $\mathcal{O}(h_T^{k+1})$-term in the diffusion-dominated regime, an $\mathcal{O}(h_T^{k+\frac12})$-term in the advection-dominated regime, and scales with intermediate powers of $h_T$ in between. Numerical results complete the exposition

HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems

A physics-based shock capturing method for unsteady laminar and turbulent flows

We present a shock capturing method for unsteady laminar and turbulent flows. The proposed approach relies on physical principles to increase selected transport coefficients and resolve unstable sharp features, such as shock waves and strong thermal and shear gradients, over the smallest distance allowed by the discretization. In particular, we devise various sensors to detect when the shear viscosity, bulk viscosity and thermal conductivity of the fluid do not suffice to stabilize the numerical solution. In such cases, the transport coefficients are increased as necessary to optimally resolve these features with the available resolution. The performance of the method is illustrated through numerical simulation of external and internal flows in transonic, supersonic, and hypersonic regimes.United States. Air Force. Office of Scientific Research (FA9550-16-1-0214)Pratt & Whitney Aircraft CompanyFundación Obra Social de La CaixaMassachusetts Institute of Technology. Office of the Dean for Graduate Education (Zakhartchenko Fellowship

A Divergence-Free and H(div)H(div)-Conforming Embedded-Hybridized DG Method for the Incompressible Resistive MHD equations

We proposed a divergence-free and $H(div)$-conforming embedded-hybridized discontinuous Galerkin (E-HDG) method for solving stationary incompressible viso-resistive magnetohydrodynamic (MHD) equations. In particular, the E-HDG method is computationally far more advantageous over the hybridized discontinuous Galerkin (HDG) counterpart in general. The benefit is even significant in the three-dimensional/high-order/fine mesh scenario. On a simplicial mesh, our method with a specific choice of the approximation spaces is proved to be well-posed for the linear case. Additionally, the velocity and magnetic fields are divergence-free and $H(div)$-conforming for both linear and nonlinear cases. Moreover, the results of well-posedness analysis, divergence-free property, and $H(div)$-conformity can be directly applied to the HDG version of the proposed approach. The HDG or E-HDG method for the linearized MHD equations can be incorporated into the fixed point Picard iteration to solve the nonlinear MHD equations in an iterative manner. We examine the accuracy and convergence of our E-HDG method for both linear and nonlinear cases through various numerical experiments including two- and three-dimensional problems with smooth and singular solutions. For smooth problems, the results indicate that convergence rates in the $L^2$ norm for the velocity and magnetic fields are optimal in the regime of low Reynolds number and magnetic Reynolds number. Furthermore, the divergence error is machine zero for both smooth and singular problems. Finally, we numerically demonstrated that our proposed method is pressure-robust