Massachusetts Institute of Technology
Viscosity stabilized adjoint method for unsteady compressible Navier-Stokes equations
Abstract
dc:description.abstractDesign optimization methods are a popular tool in computational fluid dynamics for designing components or finalizing the flow parameters of a system. The adjoint method accelerates the design process by providing gradients of the design objective with respect to the system parameters. But, typically, adjoint-based design optimization methods have used low fidelity simulations like Reynolds Averaged Navier-Stokes (RANS). To reliably capture the complex flow phenomena like turbulent boundary layers, turbulent wakes and fluid separation involved in high Reynolds number flows, high fidelity simulations like large eddy simulation (LES) are required. Unfortunately, due to the chaotic dynamics of turbulence, the adjoint method for LES diverges and produces incorrect gradients. In this thesis, the adjoint method for unsteady flow equations is modified by adding artificial viscosity to the adjoint equations. The additional viscosity stabilizes the adjoint solution and maintains reasonable accuracy of the gradients obtained from it. The accuracy of the method is assessed on multiple turbulent flow problems, including subsonic flow over a cylinder and transonic flow over a gas turbine vane. The utility of the method is then tested in performing shape optimization of the trailing edge of a transonic turbine vane. The optimal design, found using a modified gradient-based Bayesian optimization algorithm, shows approximately 15% better aero-thermal performance than the baseline design. Such design optimizations are possible due to the availability of massively parallel supercomputers. Designing high performance fluid flow solvers for the next generation supercomputers is a challenging task. In this thesis, a two-level computational graph method for writing optimized distributed flow solvers on heterogeneous architectures is presented. A checkpoint-based automatic differentiation method is used to derive the corresponding adjoint flow solver in this framework.
Degree
thesis:*- Department dc:contributor.department
- Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
- Grantor dc:publisher
- Massachusetts Institute of Technology
- Year dc:date.issued
- 2018
Author and committee
dc:creator, dc:contributor.*- Author dc:creator
-
- Talnikar, Chaitanya Anil.
- Advisor dc:contributor.advisor
-
- Qiqi Wang.
Subjects
dc:subject × 1Rights
dc:rights- Statement dc:rights
-
- MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission.
- Licence dc:rights.uri
- Language dc:language.iso
- eng
Identifiers
dc:identifier.*- Handle dc:identifier.uri
- http://hdl.handle.net/1721.1/120425
- OAI identifier oai:identifier
- oai:dspace.mit.edu:1721.1/120425