Turbulent flows occur in various fields and are a central, yet an extremely complex, topic in fluid dynamics. Understanding the behaviour of fluids is vital for multiple research areas including, but not limited to: biological models, weather prediction, and engineering design models for automobiles and aircraft. In this thesis, we study a number of fundamental problems that arise in 2D turbulent flows, using the 2D Navier-Stokes system. Introducing optimization techniques for systems described by partial differential equations (PDE), we frame these problems such that they can be solved using computational methods. We utilize adjoint calculus to build the computational framework to be implemented in an iterative gradient flow procedure, using the "optimize-then-discretize" approach. Pseudospectral methods are employed for solving PDEs in a numerically efficient manner. The use of optimization methods together with computational mathematics in this work provides an illuminating perspective on fluid mechanics.
We first apply these techniques to better understand enstrophy dissipation in 2D Navier-Stokes flows, in the limit of vanishing viscosity. By defining an optimization problem to determine optimal initial conditions, multiple branches of local maximizers were obtained each corresponding to a different mechanism producing maximum enstrophy dissipation. Viewing this quantity as a function of viscosity revealed quantitative agreement with an analytic bound, demonstrating the sharpness of this bound. We also introduce an extension of this problem, where enstrophy dissipation is maximized in the context of kinetic theory using the Boltzmann equation.
Secondly, these PDE-constrained optimization techniques were used to probe the fundamental limitations on the performance of the Leith eddy-viscosity closure model for 2D Large-Eddy Simulations of the Navier-Stokes system. Obtained by solving an optimization problem with a non-standard structure, the results demonstrate the optimal eddy viscosities do not converge to a well-defined limit as regularization and discretization parameters are refined, hence the problem of determining an optimal eddy viscosity is ill-posed.
Further extending the problem of finding optimal eddy-viscosity closures, we consider imposing an additional nonlinear constraint on the control variable in the problem, in the form of requiring the time-averaged enstrophy be preserved. To address this problem in a novel way, we employ adjoint calculus to characterize a subspace tangent to the constraint manifold, which allows one to approximately enforce the constraint. Not only do we demonstrate that this produces better results when compared to the case without constraints, but this also provides a flexible computational framework for approximate enforcement of general nonlinear constraints. Lastly in this thesis, we introduce an optimization problem to study the Kolmogorov-Richardson energy cascade, where a pathway towards solutions is outlined. / Thesis / Doctor of Philosophy (PhD)
| Identifer | oai:union.ndltd.org:mcmaster.ca/oai:macsphere.mcmaster.ca:11375/27996 |
| Date | January 2022 |
| Creators | Matharu, Pritpal |
| Contributors | Protas, Bartosz, Mathematics and Statistics |
| Source Sets | McMaster University |
| Language | English |
| Detected Language | English |
| Type | Thesis |
Page generated in 0.0022 seconds