Improving shock-capturing robustness for higher-order finite element solvers

Simulation of high speed flows where shock waves play a significant role is still an area of development in computational fluid dynamics. Numerical simulation of discontinuities such as shock waves often suffer from nonphysical oscillations which can pollute the solution accuracy. Grid adaptation, along with shock-capturing methods such as artificial viscosity, can be used to resolve the shock by targeting the key flow features for grid refinement. This is a powerful tool, but cannot proceed without first converging on an initially coarse, unrefined mesh. These coarse meshes suffer the most from nonphysical oscillations, and many algorithms abort the solve process when detecting nonphysical values. In order to improve the robustness of grid adaptation on initially coarse meshes, this thesis presents methods to converge solutions in the presence of nonphysical oscillations. A high order discontinuous Galerkin (DG) framework is used to discretize Burgers' equation and the Euler equations. Dissipation-based globalization methods are investigated using both a pre-defined continuation schedule and a variable continuation schedule based on homotopy methods, and Burgers' equation is used as a test bed for comparing these continuation methods. For the Euler equations, a set of surrogate variables based on the primitive variables (density, velocity, and temperature) are developed to allow the convergence of solutions with nonphysical oscillations. The surrogate variables are applied to a flow with a strong shock feature, with and without continuation methods, to demonstrate their robustness in comparison to the primitive variables using physicality checks and pseudo-time continuation.

Department dc:contributor.department

Massachusetts Institute of Technology. Department of Aeronautics and Astronautics.

Grantor dc:publisher

Massachusetts Institute of Technology

Year dc:date.issued

2015

Author dc:creator

• Wagner, Carlee F

Advisor dc:contributor.advisor

• David L. Darmofal.

• Aeronautics and Astronautics.

Statement dc:rights

• M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.

Licence dc:rights.uri

• http://dspace.mit.edu/handle/1721.1/7582

Language dc:language.iso

eng

Handle dc:identifier.uri

http://hdl.handle.net/1721.1/101498

OAI identifier oai:identifier

oai:dspace.mit.edu:1721.1/101498