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

Thumbnail Image
Gao, Haiyang
Major Professor
Zhi J Wang
Committee Member
Journal Title
Journal ISSN
Volume Title
Research Projects
Organizational Units
Organizational Unit
Aerospace Engineering

The Department of Aerospace Engineering seeks to instruct the design, analysis, testing, and operation of vehicles which operate in air, water, or space, including studies of aerodynamics, structure mechanics, propulsion, and the like.

The Department of Aerospace Engineering was organized as the Department of Aeronautical Engineering in 1942. Its name was changed to the Department of Aerospace Engineering in 1961. In 1990, the department absorbed the Department of Engineering Science and Mechanics and became the Department of Aerospace Engineering and Engineering Mechanics. In 2003 the name was changed back to the Department of Aerospace Engineering.

Dates of Existence

Historical Names

  • Department of Aerospace Engineering and Engineering Mechanics (1990-2003)

Related Units

Journal Issue
Is Version Of

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.

Subject Categories
Sat Jan 01 00:00:00 UTC 2011