{"id":6068,"date":"2023-06-02T15:31:43","date_gmt":"2023-06-02T13:31:43","guid":{"rendered":"https:\/\/www.mi.uni-koeln.de\/NumSim\/?p=6068"},"modified":"2023-06-02T15:31:43","modified_gmt":"2023-06-02T13:31:43","slug":"snapshot-subcell-limiting-results-for-dgsem-simulations","status":"publish","type":"post","link":"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/06\/02\/snapshot-subcell-limiting-results-for-dgsem-simulations\/","title":{"rendered":"Snapshot: Subcell limiting results for DGSEM simulations"},"content":{"rendered":"<p>DGSEM simulations of the two-dimensional Euler equations considering the astrophysical jet with Mach number $\\textrm{Ma} \\approx 2000$ [1].<br \/>\nThe computational domain, $\\Omega = [-0.5,0.5]^2$, is filled with a mono-atomic gas ($\\gamma = 5\/3$) at rest with<\/p>\n<p>\\begin{equation}<br \/>\n\\rho(x,y) = 0.5, \\qquad<br \/>\np(x,y) = 0.4127, \\qquad<br \/>\nv_1(x,y) = 0, \\qquad<br \/>\nv_2(x,y) = 0,<br \/>\n\\end{equation}<\/p>\n<p>and on the left boundary there is a hypersonic inflow with<\/p>\n<p>\\begin{equation}<br \/>\n\\rho(-0.5,y_B) = 5, ~<br \/>\np(-0.5,y_B) = 0.4127, ~<br \/>\nv_1(-0.5,y_B) = 800, ~<br \/>\nv_2(-0.5,y_B) = 0,<br \/>\n\\end{equation}<\/p>\n<p>for $y_B \\in [-0.05, 0.05]$, which corresponds to a Mach number of $\\textrm{Ma}=2156.91$ with respect to the speed of sound of the jet gas, and $\\textrm{Ma}=682.08$ with respect to the speed of sound of the ambient gas.<\/p>\n<p>A spatial resolution of $256\\times 256$ quadrilateral elements, a polynomial degree of $N=3$, and four different CFL numbers are used.<\/p>\n<p>The simulations are stabilized by using two different subcell limiting approaches, the monolithic convex limiting (MCL) and the FCT\/IDP limiting.<\/p>\n<p>The resulting density contours at $t=10^{-3}$:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-570x1024.jpg\" alt=\"\" width=\"570\" height=\"1024\" class=\"aligncenter size-large wp-image-6072\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-570x1024.jpg 570w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-167x300.jpg 167w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-768x1381.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet.jpg 850w\" sizes=\"(max-width: 570px) 100vw, 570px\" \/><\/a><\/p>\n<p>The dependence of the spatial discretization on the time-step size for FCT\/IDP methods causes the number of vortical structures to be highly dependent on the CFL number and the total number of time steps to be not inversely proportional to the CFL number.<\/p>\n<p>In fact, the amount of dissipation is reduced for small CFL numbers, which leads to lower minimum densities and higher maximum pressures in FCT\/IDP, as can be seen in this Table.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg\" alt=\"\" width=\"962\" height=\"378\" class=\"aligncenter size-full wp-image-6074\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg 962w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-300x118.jpg 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-768x302.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-500x196.jpg 500w\" sizes=\"(max-width: 962px) 100vw, 962px\" \/><\/a><\/p>\n<p><strong><br \/>\nReferences<\/strong><br \/>\n[1] Rueda-Ram\u00edrez, A., Bolm, B., Kuzmin, D. &#038; Gassner, G. (2023). Monolithic Convex Limiting for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods. arXiv preprint arXiv:2303.00374<\/p>\n","protected":false},"excerpt":{"rendered":"<p>DGSEM simulations of the two-dimensional Euler equations considering the astrophysical jet with Mach number $\\textrm{Ma} \\approx 2000$ [1]. The computational domain, $\\Omega = [-0.5,0.5]^2$, is filled with a mono-atomic gas ($\\gamma = 5\/3$) at rest with \\begin{equation} \\rho(x,y) = 0.5, &hellip; <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/06\/02\/snapshot-subcell-limiting-results-for-dgsem-simulations\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":14,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[49],"tags":[],"post_mailing_queue_ids":[],"_links":{"self":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/6068"}],"collection":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/users\/14"}],"replies":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/comments?post=6068"}],"version-history":[{"count":7,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/6068\/revisions"}],"predecessor-version":[{"id":6077,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/6068\/revisions\/6077"}],"wp:attachment":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/media?parent=6068"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/categories?post=6068"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/tags?post=6068"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}