{"id":5772,"date":"2022-11-08T14:45:47","date_gmt":"2022-11-08T13:45:47","guid":{"rendered":"https:\/\/www.mi.uni-koeln.de\/NumSim\/?p=5772"},"modified":"2022-11-09T15:27:25","modified_gmt":"2022-11-09T14:27:25","slug":"snapshot-low-pressure-kelvin-helmholtz-instability-with-dg-and-subcell-positivity-limiter","status":"publish","type":"post","link":"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/11\/08\/snapshot-low-pressure-kelvin-helmholtz-instability-with-dg-and-subcell-positivity-limiter\/","title":{"rendered":"Snapshot: High Mach number Kelvin-Helmholtz instability with DG and subcell positivity limiter"},"content":{"rendered":"<p>We solve the compressible Euler equations of gas dynamics with a fourth-order accurate entropy-stable discontinuous Galerkin (DG) method and combine it with a first-order accurate finite volume (FV) method at the node level to impose positivity of density and pressure [1,2].<br \/>\nThe initial condition of this problem is given by:<br \/>\n\\[<br \/>\n\\begin{array}{rlrl}<br \/>\n\\rho (t=0) &#038;= \\frac{1}{2}<br \/>\n+ \\frac{3}{4} B,<br \/>\n&#038;<br \/>\np (t=0) &#038;= 0.1, \\\\<br \/>\nv_1 (t=0) &#038;= \\frac{1}{2} \\left( B-1 \\right),<br \/>\n&#038;<br \/>\nv_2 (t=0) &#038;= \\frac{1}{10} \\sin(2 \\pi x),<br \/>\n\\end{array}<br \/>\n\\]<br \/>\nwith $B=\\tanh \\left( 15 y + 7.5 \\right) &#8211; \\tanh(15y-7.5)$, where $\\rho$ is the density, $\\vec{v}=(v_1,v_2)$ is the velocity, and $p$ is the pressure.<\/p>\n<p><iframe loading=\"lazy\" title=\"Low-pressure Kelvin-Helmholtz instability with DG and subcell positivity limiter\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/sbgMK9YkusI?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<p>The video shows the evolution of the density in time using $1024^2$ and degrees of freedom. We use the entropy-conservative and kinetic energy preserving flux of Ranocha [3] for the volume fluxes and the Rusanov solver for the surface fluxes of the DG and FV methods. During this simulation, the FV method acts on average in 0.000086% of the computational domain, and never more than 0.002% of the computational domain at a specific time.<\/p>\n<p>References<br \/>\n[1] A. M. Rueda-Ram\u00edrez, G. J. Gassner, A Subcell Finite Volume Positivity-Preserving Limiter for DGSEM Discretizations of the Euler Equations, WCCM-ECCOMAS2020, pp. 1\u201312.<br \/>\n[2] A. M. Rueda-Ram\u00edrez, W. Pazner, G. J. Gassner, Subcell limiting strategies for discontinuous galerkin spectral element methods, Computers &#038; Fluids 247 (2022) 105627.<br \/>\n[3] H. Ranocha, Generalised summation-by-parts operators and entropy stability of numerical methods for hyperbolic balance laws, Cuvillier Verlag, 2018.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>We solve the compressible Euler equations of gas dynamics with a fourth-order accurate entropy-stable discontinuous Galerkin (DG) method and combine it with a first-order accurate finite volume (FV) method at the node level to impose positivity of density and pressure &hellip; <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/11\/08\/snapshot-low-pressure-kelvin-helmholtz-instability-with-dg-and-subcell-positivity-limiter\/\">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\/5772"}],"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=5772"}],"version-history":[{"count":9,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5772\/revisions"}],"predecessor-version":[{"id":5784,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5772\/revisions\/5784"}],"wp:attachment":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/media?parent=5772"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/categories?post=5772"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/tags?post=5772"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}