Differential formulation of discontinuous Galerkin and related methods for compressible Euler and Navier-Stokes equations

A new approach to high-order accuracy for the numerical solution of conservation laws introduced by Huynh and extended to simplexes by the current work is renamed CPR (correction procedure or collocation penalty via reconstruction). The CPR approach employs the differential form of the equation and accounts for the jumps in flux values at the cell boundaries by a correction procedure. In addition to being simple and economical, it unifies several existing methods including discontinuous Galerkin (DG), staggered grid, spectral volume (SV), and spectral difference (SD). The approach is then extended to diffusion equation and Navier-Stokes equations. In the discretization of the diffusion terms, the BR2 (Bassi and Rebay), interior penalty, compact DG (CDG), and I-continuous approaches are used. The first three of these approaches, originally derived using the integral formulation, were recast here in the CPR framework, whereas the I-continuous scheme, originally derived for a quadrilateral mesh, was extended to a triangular mesh. The current work also includes a study of high-order curve boundaries representations. A new boundary representation based on the Bezier curve is then developed and analyzed, which is shown to have several advantages for complicated geometries. To further enhance the efficiency, the capability of h/p mesh adaptation is developed for the CPR solver. The adaptation is driven by an efficient multi-p a posteriori error estimator. P-adaptation is applied to smooth regions of the flow field while h-adaptation targets the non-smooth regions, identified by accuracy-preserving TVD marker. Several numerical tests are presented to demonstrate the capability of the technique

Hybrid Spectral Difference/Embedded Finite Volume Method for Conservation Laws

A novel hybrid spectral difference/embedded finite volume method is introduced in order to apply a discontinuous high-order method for large scale engineering applications involving discontinuities in the flows with complex geometries. In the proposed hybrid approach, the finite volume (FV) element, consisting of structured FV subcells, is embedded in the base hexahedral element containing discontinuity, and an FV based high-order shock-capturing scheme is employed to overcome the Gibbs phenomena. Thus, a discontinuity is captured at the resolution of FV subcells within an embedded FV element. In the smooth flow region, the SD element is used in the base hexahedral element. Then, the governing equations are solved by the SD method. The SD method is chosen for its low numerical dissipation and computational efficiency preserving high-order accurate solutions. The coupling between the SD element and the FV element is achieved by the globally conserved mortar method. In this paper, the 5th-order WENO scheme with the characteristic decomposition is employed as the shock-capturing scheme in the embedded FV element, and the 5th-order SD method is used in the smooth flow field. The order of accuracy study and various 1D and 2D test cases are carried out, which involve the discontinuities and vortex flows. Overall, it is shown that the proposed hybrid method results in comparable or better simulation results compared with the standalone WENO scheme when the same number of solution DOF is considered in both SD and FV elements.Comment: 27 pages, 17 figures, 2 tables, Accepted for publication in the Journal of Computational Physics, April 201

NURBS-Enhanced Finite Element Method (NEFEM)

The development of NURBS-Enhanced Finite Element Method (NEFEM) is revisited. This technique allows a seamless integration of the CAD boundary representation of the domain and the finite element method (FEM). The importance of the geometrical model in finite element simulations is addressed and the benefits and potential of NEFEM are discussed and compared with respect to other curved finite element techniques