University of Illinois Urbana-Champaign

Adjoint-based optimization for hyperbolic balance laws in the presence of discontinuities

Fikl, Alexandru

Content Files
FIKL-THESIS-2016.pdf

Permalink

https://hdl.handle.net/2142/95395

Description

Title
Adjoint-based optimization for hyperbolic balance laws in the presence of discontinuities
Author(s)
Fikl, Alexandru
Issue Date
2016-12-06
Director of Research (if dissertation) or Advisor (if thesis)
Le Chenadec, Vincent
Committee Member(s)
Sayadi, Taraneh
Department of Study
Aerospace Engineering
Discipline
Aerospace Engineering
Degree Granting Institution
University of Illinois at Urbana-Champaign
Degree Name
M.S.
Degree Level
Thesis
Date of Ingest
2017-03-01T15:49:27Z
Keyword(s)
  • Adjoint
  • Optimization
  • Hyperbolic
  • Conservation laws
  • Interface
  • Thinc
Abstract
In this thesis, we are interested in optimization in multiphase flows using discrete adjoint-based methods. The main issues we will endeavor to study are the impact on the linearized and adjoint equations of discontinuous solutions, the behavior of the THINC-family of numerical schemes under linearization and the rigorous formulation of an optimization problem when the interface is represented by a Heaviside marker function. We will study the two main discontinuous wave patterns that appear in hyperbolic balance laws: contact discontinuities and shocks. These two have very different characteristics in the forward problem and we will see that they require different treatments for the adjoint formulations as well. Specifically, shocks require additional information, in the form of the shock location, to properly define first-order variations necessary for gradient methods. Another important aspect in the numerical treatment of multiphase flows is the use of an anti-diffusive numerical scheme for the interface advection equation. The scheme we will investigate is the fairly recent THINC scheme, because it is one of the few numerical schemes in the field that is differentiable (with respect to the volume fraction). However, we will need to extend the classic formulation of the THINC scheme to force it to behave correctly in the linearized and adjoint regime. Even with these extensions, we will see that the adjoint converges everywhere except at the interface, due to the discontinuity. Finally, we define a simplified optimization problem for the volume fraction, where the velocity is given by a velocity potential instead of the Navier-Stokes equations. This allows us to specifically study the formulation of the cost functional and the viability of using a Heaviside function as a representation of the interface (as opposed to e.g. level set methods). We will show that such a formulation is indeed possible and leads to well-posed optimization problems.
Graduation Semester
2016-12
Type of Resource
text
Permalink
http://hdl.handle.net/2142/95395
Copyright and License Information
Copyright 2016 Alexandru Fikl

Owning Collections

Graduate Dissertations and Theses at Illinois PRIMARY
Graduate Theses and Dissertations at Illinois
Dissertations and Theses - Aerospace Engineering