In this thesis we implement and study the effects of different convective outflow boundary conditions for the high order spectral element solver Nek5000 in the context of solving convective-dominated fluid flow problems. By numerical testing we show that the convective boundary conditions preserve the spatial and temporal convergence rates of the solver. We also study highly convective test cases such as a single vortex propagating through the outflow boundary, and the typical Kármán vortex shedding problem to analyze the accuracy and stability. A detailed comparison with the natural boundary condition that corresponds to the variational form of the incompressible Navier–Stokes equations (the Nek5000 “O” condition), and a stabilized version of it (by Dong et al. (2014)), are also presented. Our results show a major advantage of using the convective boundary conditions over the natural counterpart in solving convective problems, both according to stability and accuracy. Analytic and numerical results show that the natural condition has big stability problems for high Reynolds numbers, which make the use of stabilization methods or damping regions crucial. But, the (Dong) stabilized natural condition does not improve accuracy, and damping regions are computationally expensive. The convective conditions show very good accuracy if its convection speed is approximated accurately, and our results indicate that it can be used without damping regions efficiently. Our results also show that the magnitude of reflections significantly depends on the amplitude of the disturbances that move through the boundary. The convective boundary condition can handle large disturbances without producing significant reflections, while the natural one or a stabilized version of it in general can not. / I det här examensarbetet har vi implementerat och studerat effekterna av olika konvektiva randvillkor för spektralelement lösaren Nek5000 vid beräkningar av konvektivt dominanta flödesproblem. Med hjälp av numeriska tester bevisar vi att de nya implementationerna bevarar lösarens konvergens i både rum och tid. Vi studerar noggrannheten hos de konvektiva randvillkoren genom konvektivt dominanta testfall i form av en ensam virvel som propagerar genom utflödet, samt det klassiska Kármán-virvel-gata problemet. En detaljerad jämförelse med det naturliga randvillkoret tillhörande den svaga formuleringen av de inkompressibla Navier-Stokes ekvationerna (Nek5000 “O”) och en stabiliserad version av denna är också presenterade. Våra resultat visar tydliga fördelar med att använda de konvektiva randvillkoren mot det naturliga vid lösningar av konvektiva problem, både stabilitetsmässigt och noggrannhetsmässigt. Analytiska och numeriska resultat visar att det naturliga randvillkoret har stora stabilitetsproblem vid höga Reynolds-tal, vilket medför att specifika stabilitetsversioner eller dämpningsregioner måste användas. Men stabiliserande naturliga randvillkor (Dong) förbättrar inte noggrannheten och dämpningsregioner är dyra beräkningsmässigt. De konvektiva randvillkoren har uppvisat en väldigt god noggrannhet om konvektionshastigheten i villkoren är noggranna approximationer. Analyser av amplituden hos reflektioner har också undersökts. Våra resultat visar ett signifikant linjärt förhållande mellan storleken på störningar som genomborrar utflödes-randen och de reflektioner dessa störningar skapar. De konvektiva randvillkoren visar sig klara starka störningar bra, vilket det naturliga och den stabiliserade versionen av det naturliga randvillkoret generellt sett inte gör.
| Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:kth-164094 |
| Date | January 2015 |
| Creators | Boström, Erik |
| Publisher | KTH, Numerisk analys, NA |
| Source Sets | DiVA Archive at Upsalla University |
| Language | English |
| Detected Language | Swedish |
| Type | Student thesis, info:eu-repo/semantics/bachelorThesis, text |
| Format | application/pdf |
| Rights | info:eu-repo/semantics/openAccess |
| Relation | TRITA-MAT-E ; 2015:11 |
Page generated in 0.0014 seconds