Energy and entropy stable numerical methods with injected boundary conditions

Abstract

I denne avhandlingen studerer vi de kompressible Navier-Stokes-likningene formulert med både adiabatiske veggrandvilkår og fjernfeltvilkår. Selv om det er ukjent om disse likningene er velformulerte er de av stor interesse, og de er mye brukt innen numerisk fluiddynamikk. Et resultat av Strang (1964) sier at for ikke-lineære problem diskretisert ved hjelp av en differansemetode som er lineærstabil, er denne metoden konvergent for glatte løsninger. Altså finnes det teori vi kan bruke i analysen av Navier-Stokes-likningene. Derfor studerer vi her teori for velformulerte lineære problem, og stabilitet for numeriske metoder. Dette gjøres både for de kompressible Navier-Stokes-likningene, men også for lineære partielle differensiallikninger som modellproblem. Videre utleder vi entropiestimat for de ikke-lineære Navier-Stokes-likningene, et estimat som virker som et kriterium for den svake løsningen vi leter etter; den skal i tillegg til likningene tilfredsstille termodynamikkens andre lov.

Hovedfokuset ved dette arbeidet er stabil håndtering av de adiabatiske veggrandvilkårene og fjernfeltvilkår for Navier-Stokes-likningene. Vi beviser at heftelsesvilkåret (eng.: no- slip condition) kan bli implementert eksakt og fremdeles resultere i et entropiestimat når teknikken brukes i kombinasjon med delvissummasjonsoperatorer (SBP-operatorer) som har diagonale normmatriser og randmatriser. Vi introduserer også en ny metodikk for å sette fjernfeltvilkår, og beviser at den fører til et entropistabilt skjema for de kompressible Navier-Stokes-likningene. Teknikken er i tillegg lineært velformulert. Gjennom hele arbeidet bruker vi SBP-operatorer på grunn av deres gode stabilitetsegenskaper. Vi beviser også at en litt endret versjon av SBP-operatoren som tilnærmer den andrederiverte ved hjelp av endelig-volummetoden gitt av Chandrashekar (2016) er (svakt) konsistent, noe som gjør den egnet til å diskretisere de viskøse leddene i Navier-Stokes-likningene på ustrukturerte gitter.

The compressible Navier-Stokes equation subject to both adiabatic wall boundary conditions and far-field boundary conditions are studied in this thesis. Although the well- posedness of these equations is generally unknown, they are of wide interest and are extensively used in computational fluid dynamics. A result by Strang (1964) states that if a non-linear problem is discretised using a difference method that is linearly stable, then this method is convergent for smooth solutions. That is, there exists theory we can use in the analysis of the Navier-Stokes equations. Thus, we study linear well-posedness and stability of numerical schemes both in the context of the compressible Navier-Stokes equations, but also linear partial differential equations as model problems. Furthermore, entropy estimates are derived for the fully non-linear Navier-Stokes equations, which pose as an admissibility criterion for the relevant weak solution we seek; it should additionally satisfy the second law of thermodynamics.

The main focus of this work is the stable imposition of the adiabatic wall and far-field boundary conditions for the Navier-Stokes equations. In particular, we prove that the no-slip condition can be imposed strongly and still yield an entropy estimate when used in combination with diagonal-norm summation-by-parts (SBP) operators with diagonal boundary operators. Furthermore, we introduce a new methodology for setting far- field boundary conditions, and prove that it leads to an entropy stable scheme for the compressible Navier-Stokes equations. The procedure is additionally linearly well-posed. Throughout, we employ SBP operators due to their remarkable stability properties. We also prove that a slightly modified version of the finite-volume SBP approximation of the second-derivative given by Chandrashekar (2016) is (weakly) consistent, thus making it suitable for discretising the viscous terms of the Navier-Stokes equations on unstructured grids.