{"id":5483,"date":"2021-10-19T18:23:51","date_gmt":"2021-10-19T16:23:51","guid":{"rendered":"https:\/\/www.mi.uni-koeln.de\/NumSim\/?p=5483"},"modified":"2021-10-19T18:23:51","modified_gmt":"2021-10-19T16:23:51","slug":"snapshot-simulation-of-a-stiff-multi-species-reacting-flow","status":"publish","type":"post","link":"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/10\/19\/snapshot-simulation-of-a-stiff-multi-species-reacting-flow\/","title":{"rendered":"Snapshot: Simulation of a stiff multi-species reacting flow"},"content":{"rendered":"

Simulation of an asymmetric two-dimensional stiff multi-species reacting flow for the compressible Euler equations. The computation domain [0,6]x[0,2]<\/code> of this detonation problem is split into three zones. Zone A [0,0.5]x[0,2]<\/code> as well as zone B (0.5,6]x[1.2,2]<\/code> is filled with burnt gas, namely O2,OH<\/code> and H2O<\/code>. The last zone C (0.5,6]x[0,1.2)<\/code> contains unburnt gas, namely H2<\/code> and O2<\/code>. <\/p>\n

The reacting model consists of two reactions H2+O2->2OH<\/code>, 2OH+H2->2H2O<\/code> and five species H2,O2,OH,H2O,N2<\/code> whereas reaction 1 has a smaller ignition temperature as well as a much faster chemical rate than reaction 2. Species N2<\/code> acts as a dilute catalyst here. <\/p>\n

This example has been calculated using Trixi.jl with a P4est mesh using 3072x1024=3.145.728<\/code> elements, Dirichlet boundary conditions in x-direction as well as solid wall boundaries in y-direction. The final time is T = 0.1. The convection part is solved with FV whereas the chemical network is solved with KROME. (Note: Not all of these features are currently available in the main code of Trixi.jl.)<\/em><\/p>\n

Pressure distribution at final time T=0.1:
\n\"\"<\/a><\/p>\n

Temperature distribution at final time T=0.1:
\n\"\"<\/a><\/p>\n

Density distribution at final time T=0.1:
\n\"\"<\/a><\/p>\n

OH-fraction distribution at final time T=0.1:
\n\"\"<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"

Simulation of an asymmetric two-dimensional stiff multi-species reacting flow for the compressible Euler equations. The computation domain [0,6]x[0,2] of this detonation problem is split into three zones. Zone A [0,0.5]x[0,2] as well as zone B (0.5,6]x[1.2,2] is filled with burnt … Continue reading →<\/span><\/a><\/p>\n","protected":false},"author":12,"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\/5483"}],"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\/12"}],"replies":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/comments?post=5483"}],"version-history":[{"count":24,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5483\/revisions"}],"predecessor-version":[{"id":5515,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/posts\/5483\/revisions\/5515"}],"wp:attachment":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/media?parent=5483"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/categories?post=5483"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/tags?post=5483"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}