{"id":5291,"date":"2021-04-01T18:21:00","date_gmt":"2021-04-01T16:21:00","guid":{"rendered":"https:\/\/www.mi.uni-koeln.de\/NumSim\/?p=5291"},"modified":"2021-04-01T18:36:59","modified_gmt":"2021-04-01T16:36:59","slug":"snapshot-simulation-of-a-kelvin-helmholtz-instability-using-second-order-finite-volume-schemes-and-fourth-order-discontinuous-galerkin-methods","status":"publish","type":"post","link":"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/04\/01\/snapshot-simulation-of-a-kelvin-helmholtz-instability-using-second-order-finite-volume-schemes-and-fourth-order-discontinuous-galerkin-methods\/","title":{"rendered":"Snapshot: Simulation of a Kelvin-Helmholtz instability using second order Finite Volume schemes and fourth order Discontinuous Galerkin methods"},"content":{"rendered":"<p>We present in-viscid and viscous simulations of a Kelvin-Helmholtz instability using second a order accurate monotoniced-central finite volume (FV) method and a fourth order accurate discontinuous Galerkin (DG) method. The initial condition is given by [1]:<\/p>\n<p>$$\\rho (t=0) = \\frac{1}{2}<br \/>\n+ \\frac{3}{4} B,<br \/>\n~~~~~~~~~<br \/>\np (t=0) = 1,~~~~~~~~~~<br \/>\n$$<\/p>\n<p>$$<br \/>\nv_1 (t=0) = \\frac{1}{2} \\left( B-1 \\right),<br \/>\n~~~~~~~<br \/>\nv_2 (t=0) = \\frac{1}{10} \\sin(2 \\pi x),<br \/>\n$$<\/p>\n<p>with $$B=\\tanh \\left( 15 y + 7.5 \\right) &#8211; \\tanh(15y-7.5).$$<\/p>\n<p>We first present the FV results at end time $t=3.7$, which are computed using a monotoniced-central second order discretization of the Euler equations of gas dynamics on uniform grids.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5307\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-1022x1024.png 1022w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-768x769.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density.png 1363w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5308\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-768x767.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density.png 1367w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>The next results use a fourth order DG discretization of the Navier-Stokes equations on uniform grids using $Re=320.000$ at end time $t=3.7$. The highest resolution (4096\u00b2 DOFs) is a direct numerical simulation (DNS) of the problem, where all scales are resolved.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5309\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-768x768.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale.png 1368w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5310\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-768x769.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density.png 1367w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>It is remarkable that the numerical dissipation of the second order FV scheme causes the in-viscid simulation with 2048\u00b2 DOFs to look very similar to the viscous DNS solution at $Re=320.000$.<\/p>\n<p>[1] A.M. Rueda-Ram\u00edrez, G.J Gassner (2021). <em>A Subcell Finite Volume Positivity-Preserving Limiter for DGSEM Discretizations of the Euler Equations<\/em>. <a href=\"https:\/\/arxiv.org\/pdf\/2102.06017.pdf\">https:\/\/arxiv.org\/pdf\/2102.06017.pdf<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>We present in-viscid and viscous simulations of a Kelvin-Helmholtz instability using second a order accurate monotoniced-central finite volume (FV) method and a fourth order accurate discontinuous Galerkin (DG) method. The initial condition is given by [1]: $$\\rho (t=0) = \\frac{1}{2} &hellip; <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/04\/01\/snapshot-simulation-of-a-kelvin-helmholtz-instability-using-second-order-finite-volume-schemes-and-fourth-order-discontinuous-galerkin-methods\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":13,"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\/5291"}],"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\/13"}],"replies":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/comments?post=5291"}],"version-history":[{"count":20,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5291\/revisions"}],"predecessor-version":[{"id":5327,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5291\/revisions\/5327"}],"wp:attachment":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/media?parent=5291"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/categories?post=5291"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/tags?post=5291"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}