One of the key open problems in the field of theoretical fluid mechanics concerns the possibility of the singularity formation in solutions of the 3D Navier-Stokes system in finite time. This phenomenon is associated with the behaviour of the enstrophy, which is an L2 norm of the vorticity and must become unbounded if such a singularity occurs. Although there is no blow-up in the 2D Navier-Stokes equation, we would like to investigate how much enstrophy can a planar incompressible flow in a bounded domain produce given certain initial enstrophy. We address this issue by formulating an optimization problem in which the time derivative of the enstrophy serves as the objective functional and solve it using tools of the optimization theory and calculus of variations. We propose an efficient computational approach which is based on the iterative steepest-ascent procedure. In addition, we introduce an easy-to-implement method of computing the gradient of the objective functional. Finally, we present computational results addressing the key question of this project and provide numerical evidence that the maximum enstrophy growth exhibits the scaling dE/dt ~ C*E*E for C>0 and very small E. All computations are performed using the Chebyshev spectral method. / Thesis / Master of Science (MSc) / For many decades, scientists have been investigating fundamental aspects of the Navier-Stokes equation, a central mathematical model arising in fluid mechanics. Although the equation is widely used by engineers to describe numerous flow phenomena, it is still an open question whether the Navier-Stokes system always admits physically meaningful solutions. To address this issue, we want to explore its mathematical aspects deeper by analyzing the behaviour of the enstrophy, which is a quantity associated with the vorticity of the flow and a convenient measure of the regularity of the solution. In this study, we consider a planar and incompressible flow bounded by solid walls. Using basic tools of mathematical analysis and optimization theory, we propose a computational method enabling us to find out how much enstrophy can such a flow produce instantaneously. We present numerical evidence that this instantaneous growth of enstrophy has a well-defined asymptotic behavior, which is consistent with physical assumptions.
| Identifer | oai:union.ndltd.org:mcmaster.ca/oai:macsphere.mcmaster.ca:11375/23087 |
| Date | January 2017 |
| Creators | Sliwiak, Adam |
| Contributors | Protas, Bartosz, Computational Engineering and Science |
| Source Sets | McMaster University |
| Language | English |
| Detected Language | English |
| Type | Thesis |
Page generated in 0.0026 seconds