Error Analysis Of Summation By Parts Formulations PDF Download

In this thesis, we consider the numerical solution of initial boundary value problems (IBVPs). Boundary and interface conditions are derived such that the IBVP under consideration is well-posed. We also study the dual problem and the related dual boundary/interface conditions. Once the continuous problem is analyzed, we use finite difference operators with the Summation- By-Parts property (SBP) and a weak boundary/interface treatment using the Simultaneous-Approximation-Terms (SAT) technique to construct high-order accurate numerical schemes. We focus in particular on stability, conservation and dual consistency. The energy method is used as our main analysis tool for both the continuous and numerical problems. The contributions of this thesis can be divided into two parts. The first part focuses on the coupling of different IBVPs. Interface conditions are derived such that the continuous problem satisfy an energy estimate and such that the discrete problem is stable. In the first paper, two hyperbolic systems of different size posed on two domains are considered. We derive the dual problem and dual interface conditions. It is also shown that a specific choice of penalty matrices leads to dual consistency. As an application, we study the coupling of the Euler and wave equations. In the fourth paper, we examine how to couple the compressible and incompressible Navier-Stokes equations. In order to obtain a sufficient number of interface conditions, the decoupled heat equation is added to the incompressible equations. The interface conditions include mass and momentum balance and two variants of heat transfer. The typical application in this case is the atmosphere-ocean coupling. The second part of the thesis focuses on the relation between the primal and dual problem and the relation between dual consistency and conservation. In the second and third paper, we show that dual consistency and conservation are equivalent concepts for linear hyperbolic conservation laws. We also show that these concepts are equivalent for symmetric or symmetrizable parabolic problems in the fifth contribution. The relation between the primal and dual boundary conditions for linear hyperbolic systems of equations is investigated in the sixth and last paper. It is shown that for given well-posed primal/dual boundary conditions, the corresponding well-posed dual/primal boundary conditions can be obtained by a simple scaling operation. It is also shown how one can proceed directly from the well-posed weak primal problem to the well-posed weak dual problem. Den här avhandlingen handlar om numeriska metoder för att lösa initial och randvärdes problem. Studien fokuserar på härledningen av rand/kopplingsvillkor som garanterar välställdhet. Det duala problemet och dess duala rand/kopplingsvillkor studeras också. Dessa problem diskretiseras genom att använda noggranna finita differensscheman på SBP-form (eng. summation-by-parts), kombinerat med en svag randbehandling benämnd SAT (eng. simultaneous approximation term). Vi fokuserar särskilt på stabilitet, konservation och dualkonsistens. Det främsta analysverktyget för både det kontinuerliga och diskreta problemet är energimetoden. Den första delen av avhandlingen behandlar välställdhet och stabilitet för koppling av olika system av ekvationer. Kopplingsvillkoren är härledda så att det kontinuerliga problemet uppfyller en energiuppskattning och så att det diskreta problemet är stabilt. I den första artikeln görs analysen för koppling av två olika hyperboliska system på första ordningens form. Som tillämpning kopplar vi Euler och vågekvationerna. Koppling mellan kompressibla och inkompressibla Navier-Stokes ekvationer studeras i den fjärde artikeln. För att få rätt antal kopplingsvillkor lägger vi till värmeledningsekvationen till de inkompressibla ekvationerna. Kopplingsvillkoren innefattar massans och rörelsemängdens bevarande samt två varianter av värmeöverföring. Den typiska tillämpningen är koppling mellan atmosfär och hav. Den andra delen av avhandlingen fokuserar på relationen mellan det primära och duala problemet och relationen mellan dualkonsistens och konservation. I den andra och tredje artikeln visar vi att dualkonsistens och konservation är ekvivalenta koncept för linjära hyperboliska konserveringslagar. I den femte artikeln, visas att dessa koncept är ekvivalenta även för paraboliska problem. Relationen mellan de primära och duala randvilkoren för linjära hyperboliska system av ekvationer i två dimensioner studeras i den sista artikeln. Vi visar att primära/duala välställda randvilkor ger duala/primära välställda randvilkor genom en enkel skalningsoperation. Det visas också att man kan gå direkt från det välställda svaga primära problemet till det välställda svaga duala problemet.