Numerical approximations using high order finite differences on summation-byparts (SBP) form are investigated for discontinuous and fully nonlinear systems of partial differential equations. Stability and conservation properties of the approximations are obtained through a weak imposition of interface and boundary conditions with the simultaneous-approximation-term (SAT) technique. The SBP-SAT approximations replicate the continuous integration by parts rule. From this property, well-posedness and integral properties of the continuous problem are mimicked, and energy estimates leading to stability are obtained. The first part of the thesis focuses on the simulations of discontinuous linear advection problems. An artificial interface is introduced, separating parts of the spatial domain characterized by different wave speeds. A set of flexible stability conditions at the interface are derived, which can be adapted to yield conservative or non-conservative approximations. This model can be interpreted as a simplified version of nonlinear problems involving jumps at shocks, or as a prototypical of wave propagation through different materials. In the second part of the thesis, the vorticity/stream function formulation of the nonlinear momentum equation for an incompressible inviscid fluid is considered. SBP operators are used to derive a new Arakawa-like Jacobian with mimetic properties by combining different consistent approximations of the convection terms. Energy and enstrophy conservation is obtained for periodic problems using schemes with arbitrarily high order of accuracy. These properties are crucial for long-term numerical calculations in climate and weather forecasts or ocean circulation predictions. The third and final contribution of the thesis is dedicated to the incompressible Navier-Stokes problem. First, different completely general formulations of energy bounding boundary conditions are derived for the nonlinear equations. The boundary conditions can be used at both far field and solid wall boundaries. The discretisation in time and space with weakly imposed initial and boundary conditions using the SBP-SAT framework is proved to be stable and the divergence free condition is approximated with the design order of the scheme. Next, the same formulations are considered in a linearised setting, whereupon the spectra associated with the initial boundary value problem and its SBP-SAT discretisation are derived using the Laplace-Fourier technique. The influence of different boundary conditions on the spectrum and in particular the convergence to steady state is studied. / Numeriska approximationer av ekvationer som styr fysikaliska lagar är avgörande i många tillämpningar. Förutom en matematisk modell som kan fånga huvuddragen i ett verkligt problem är det nödvändigt att kunna utföra tillförlitliga simuleringar. Denna avhandling behandlar numeriska approximationer som med hög noggrannhet bevarar både rent matematiska aspekter av ekvationerna så väl som viktiga egenskaper hos modellen. Dessutom ges särskild uppmärksamhet åt modeller med diskontinuiteter och icke-linjära beteenden. Den första delen av avhandlingen handlar om diskontinuerliga problem. Det fysiska rummet kan ha olika egenskaper i olika regioner, något som kan resultera i instabila lösningar. Tillvägagångssättet består av att införa artificiella gränssnitt som skiljer dessa regioneråt. På detta sätt kan varje region behandlas separat, men på liknande sätt. Exempel på naturliga tillämningsområden är vågutbredning genom olika material och jordbävningssimuleringar. I den andra delen av avhandlingen visar vi att om den numeriska approximationen imiterar partiell integration, då följer också de väsentliga egenskaperna hos modellen på ett naturligt sätt. Att fysikaliska egenskaper bevaras är nödvändigt för att bibehålla stabilitet under långa simuleringstider för bland annat geofysiska problem. Den sista delen av avhandlingen är ägnasåt en av de mest använda modellerna inom strömningsmekanik, nämligen Navier-Stokes ekvationer. Studien fokuserar på härledningen av randvillkor som garanterar att lösningen inte växer på ett oförutsett och okontrollerat vis. Slutligen visas att de härledda randvillkoren på ett korrekt och noggrant sätt återskapar den dissipativa mekanism som ger upphov till jämviktstillstånd.
| Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:liu-140497 |
| Date | January 2017 |
| Creators | La Cognata, Cristina |
| Publisher | Linköpings universitet, Beräkningsmatematik, Linköpings universitet, Tekniska fakulteten, Linköping |
| Source Sets | DiVA Archive at Upsalla University |
| Language | English |
| Detected Language | English |
| Type | Doctoral thesis, comprehensive summary, info:eu-repo/semantics/doctoralThesis, text |
| Format | application/pdf |
| Rights | info:eu-repo/semantics/openAccess |
| Relation | Linköping Studies in Science and Technology. Dissertations, 0345-7524 ; 1880 |
Page generated in 0.002 seconds