diff --git a/prm/benchmarks/supersonic-nozzle-2d.prm b/prm/benchmarks/supersonic-nozzle-2d.prm new file mode 100644 index 000000000..084814f45 --- /dev/null +++ b/prm/benchmarks/supersonic-nozzle-2d.prm @@ -0,0 +1,48 @@ +subsection A - TimeLoop + set basename = supersonic-nozzle-2d + + set enable output full = true + + set final time = 0.01 + set timer granularity = 0.000001 +end + +subsection B - Equation + set dimension = 2 + set equation = euler + set gamma = 1.4 +end + +subsection C - Discretization + set geometry = reader + set mesh refinement = 2 + set mesh repartitioning = false + subsection reader + set filename = prm/meshes/nozzle_0.5-3mm.msh + end +end + +subsection E - InitialValues + + set configuration = contrast + set direction = 0., 1. + set position = 0., -0.424 + + subsection contrast + set primitive state right= 0.00000002323, 0.0, 2. + set primitive state left= 0.06969, 4.76, 6000000. + end +end +subsection H - TimeIntegrator + set cfl min = 0.8 + set cfl max = 0.8 + set cfl recovery strategy = none + + set time stepping scheme = erk 33 +end + + +subsection I - VTUOutput + set use mpi io = true + set schlieren beta = 10 +end \ No newline at end of file diff --git a/prm/meshes/nozzle_0.5-3mm.geo b/prm/meshes/nozzle_0.5-3mm.geo new file mode 100644 index 000000000..db9e0b3f0 --- /dev/null +++ b/prm/meshes/nozzle_0.5-3mm.geo @@ -0,0 +1,81 @@ +// Gmsh project created on Fri Jun 21 11:10:30 2024 +Point(1) = {-0.10, -0.145, 0.0, 1}; +Point(2) = {0.10, -0.145, 0.0, 1}; +Point(3) = {-0.025, -0.1, 0.0, 1}; +Point(4) = {0.025, -0.1, 0.0, 1}; +Point(5) = {-0.025, 0., 0.0, 1}; +Point(6) = {0.025, 0., 0.0, 1}; +Point(7) = {-0.15, 1., 0.0, 1}; +Point(8) = {0.15, 1., 0.0, 1}; +Point(9) = {-0.4, 1., 0.0, 1}; +Point(10) = {0.4, 1., 0.0, 1}; +Point(11) = {-0.4, 1.55, 0.0, 1}; +Point(12) = {0.4, 1.55, 0.0, 1}; +Point(13) = {-0.1, -0.425, 0., 1}; +Point(14) = {0.1, -0.425, 0., 1}; +Line(3) = {3, 1}; +//+ +Line(4) = {4, 2}; +//+ +Line(5) = {11, 9}; +//+ +Line(6) = {7, 9}; +//+ +Line(7) = {11, 12}; +//+ +Line(8) = {12, 10}; +//+ +Line(9) = {10, 8}; +//+ +Line(10) = {7, 8}; +//+ +Line(11) = {3, 4}; +//+ +Line(12) = {1, 2}; + +//+ +Line(13) = {7, 5}; +//+ +Line(14) = {8, 6}; +//+ +Line(15) = {5, 3}; +//+ +Line(16) = {6, 4}; +//+ +Line(17) = {5, 6}; +//+ +Line(18) = {1, 13}; +//+ +Line(19) = {13, 14}; +//+ +Line(20) = {14, 2}; +//+ +Curve Loop(1) = {6, -5, 7, 8, 9, -10}; +//+ +Plane Surface(1) = {1}; +//+ +Curve Loop(2) = {13, 17, -14, -10}; +//+ +Plane Surface(2) = {2}; +//+ +Curve Loop(3) = {15, 11, -16, -17}; +//+ +Plane Surface(3) = {3}; +//+ +Curve Loop(4) = {3, 12, -4, -11}; +//+ +Plane Surface(4) = {4}; +//+ +Curve Loop(5) = {18, 19, 20, -12}; +//+ +Plane Surface(5) = {5}; +//+ + +//+ +Physical Surface("mesh", 21) = {1, 2, 3, 5, 4}; +//+ +Physical Curve(3) = {5, 6, 8, 9, 14, 13, 15, 16, 3, 4, 18, 20}; +//+ +Physical Curve(4) = {19}; +//+ +Physical Curve(0) = {7}; diff --git a/prm/meshes/nozzle_0.5-3mm.msh b/prm/meshes/nozzle_0.5-3mm.msh new file mode 100644 index 000000000..75ec8203f --- /dev/null +++ b/prm/meshes/nozzle_0.5-3mm.msh @@ -0,0 +1,1579 @@ +$MeshFormat +4.1 0 8 +$EndMeshFormat +$PhysicalNames +1 +2 21 "mesh" +$EndPhysicalNames +$Entities +14 18 5 0 +1 -0.1 -0.145 0 0 +2 0.1 -0.145 0 0 +3 -0.025 -0.1 0 0 +4 0.025 -0.1 0 0 +5 -0.025 0 0 0 +6 0.025 0 0 0 +7 -0.15 1 0 0 +8 0.15 1 0 0 +9 -0.4 1 0 0 +10 0.4 1 0 0 +11 -0.4 1.55 0 0 +12 0.4 1.55 0 0 +13 -0.1 -0.425 0 0 +14 0.1 -0.425 0 0 +3 -0.1 -0.145 0 -0.025 -0.1 0 1 3 2 3 -1 +4 0.025 -0.145 0 0.1 -0.1 0 1 3 2 4 -2 +5 -0.4 1 0 -0.4 1.55 0 1 3 2 11 -9 +6 -0.4 1 0 -0.15 1 0 1 3 2 7 -9 +7 -0.4 1.55 0 0.4 1.55 0 1 0 2 11 -12 +8 0.4 1 0 0.4 1.55 0 1 3 2 12 -10 +9 0.15 1 0 0.4 1 0 1 3 2 10 -8 +10 -0.15 1 0 0.15 1 0 0 2 7 -8 +11 -0.025 -0.1 0 0.025 -0.1 0 0 2 3 -4 +12 -0.1 -0.145 0 0.1 -0.145 0 0 2 1 -2 +13 -0.15 0 0 -0.02499999999999999 1 0 1 3 2 7 -5 +14 0.02499999999999999 0 0 0.15 1 0 1 3 2 8 -6 +15 -0.025 -0.1 0 -0.025 0 0 1 3 2 5 -3 +16 0.025 -0.1 0 0.025 0 0 1 3 2 6 -4 +17 -0.025 0 0 0.025 0 0 0 2 5 -6 +18 -0.1 -0.425 0 -0.1 -0.145 0 1 3 2 1 -13 +19 -0.1 -0.425 0 0.1 -0.425 0 1 4 2 13 -14 +20 0.1 -0.425 0 0.1 -0.145 0 1 3 2 14 -2 +1 -0.4 1 0 0.4 1.55 0 1 21 6 6 -5 7 8 9 -10 +2 -0.15 0 0 0.15 1 0 1 21 4 13 17 -14 -10 +3 -0.025 -0.1 0 0.025 0 0 1 21 4 15 11 -16 -17 +4 -0.1 -0.145 0 0.1 -0.1 0 1 21 4 3 12 -4 -11 +5 -0.1 -0.425 0 0.1 -0.145 0 1 21 4 18 19 20 -12 +$EndEntities +$Nodes +37 467 1 467 +0 1 0 1 +1 +-0.1 -0.145 0 +0 2 0 1 +2 +0.1 -0.145 0 +0 3 0 1 +3 +-0.025 -0.1 0 +0 4 0 1 +4 +0.025 -0.1 0 +0 5 0 1 +5 +-0.025 0 0 +0 6 0 1 +6 +0.025 0 0 +0 7 0 1 +7 +-0.15 1 0 +0 8 0 1 +8 +0.15 1 0 +0 9 0 1 +9 +-0.4 1 0 +0 10 0 1 +10 +0.4 1 0 +0 11 0 1 +11 +-0.4 1.55 0 +0 12 0 1 +12 +0.4 1.55 0 +0 13 0 1 +13 +-0.1 -0.425 0 +0 14 0 1 +14 +0.1 -0.425 0 +1 3 0 3 +15 +16 +17 +-0.04374999999083107 -0.1112499999944986 0 +-0.06249999998164187 -0.1224999999889851 0 +-0.08124999999081799 -0.1337499999944908 0 +1 4 0 3 +18 +19 +20 +0.04374999999083107 -0.1112499999944986 0 +0.06249999998164187 -0.1224999999889851 0 +0.08124999999081799 -0.1337499999944908 0 +1 5 0 7 +21 +22 +23 +24 +25 +26 +27 +-0.4 1.48042738718979 0 +-0.4 1.410854774379371 0 +-0.4 1.341282161599911 0 +-0.4 1.271709548820521 0 +-0.4 1.203782161667184 0 +-0.4 1.135854774513998 0 +-0.4 1.06792738725692 0 +1 6 0 3 +28 +29 +30 +-0.212500000030537 1 0 +-0.2750000000611938 1 0 +-0.3375000000306223 1 0 +1 7 0 7 +31 +32 +33 +34 +35 +36 +37 +-0.2999311357117285 1.55 0 +-0.1998622714235879 1.55 0 +-0.09979340712695361 1.55 0 +0.0002754571696805819 1.55 0 +0.1002065929710035 1.55 0 +0.200137728772232 1.55 0 +0.3000688643860884 1.55 0 +1 8 0 7 +38 +39 +40 +41 +42 +43 +44 +0.4 1.48050852307638 0 +0.4 1.411017046152713 0 +0.4 1.341525569251003 0 +0.4 1.272034092349278 0 +0.4 1.204025569234956 0 +0.4 1.136017046120932 0 +0.4 1.068008523060349 0 +1 9 0 3 +45 +46 +47 +0.3375000000307191 1 0 +0.2750000000611938 1 0 +0.212500000030571 1 0 +1 10 0 3 +48 +49 +50 +-0.07500000003668431 1 0 +-7.34325378282108e-11 1 0 +0.07499999996320933 1 0 +1 11 0 3 +51 +52 +53 +-0.01250000000611455 -0.1 0 +-1.223875861766643e-11 -0.1 0 +0.01249999999388059 -0.1 0 +1 12 0 11 +54 +55 +56 +57 +58 +59 +60 +61 +62 +63 +64 +-0.08314162103424362 -0.145 0 +-0.06628324206849162 -0.145 0 +-0.04942486310611849 -0.145 0 +-0.03256648414374536 -0.145 0 +-0.01590055059858765 -0.145 0 +0.0007653829465739115 -0.145 0 +0.01743131649474372 -0.145 0 +0.0340972500428855 -0.145 0 +0.05057293755344688 -0.145 0 +0.06704862506402059 -0.145 0 +0.0835243125319933 -0.145 0 +1 13 0 19 +65 +66 +67 +68 +69 +70 +71 +72 +73 +74 +75 +76 +77 +78 +79 +80 +81 +82 +83 +-0.1404381934189695 0.9235055473517561 0 +-0.1308763868379172 0.847011094703338 0 +-0.1213145802573996 0.7705166420591969 0 +-0.1117527736768696 0.694022189414957 0 +-0.1024623257724813 0.6196986061798502 0 +-0.09317187786809183 0.5453750229447347 0 +-0.08388142997992776 0.4710514398394221 0 +-0.07459098209172726 0.3967278567338182 0 +-0.06830982519176859 0.3464786015341488 0 +-0.06202866829180992 0.2962293463344794 0 +-0.05574751137551688 0.2459800910041351 0 +-0.04946635445922384 0.1957308356737908 0 +-0.04566038193129304 0.1652830554503444 0 +-0.04185440940336224 0.134835275226898 0 +-0.03804843684444059 0.1043874947555248 0 +-0.03424246428551893 0.07393971428415147 0 +-0.03193184821723438 0.0554547857378751 0 +-0.02962123214894985 0.03696985719159884 0 +-0.02731061607447492 0.01848492859579942 0 +1 14 0 19 +84 +85 +86 +87 +88 +89 +90 +91 +92 +93 +94 +95 +96 +97 +98 +99 +100 +101 +102 +0.1404417120467 0.9235336963736003 0 +0.1308834240933762 0.84706739274701 0 +0.1213251361411772 0.7706010891294176 0 +0.1117668481889598 0.6941347855116785 0 +0.1024806589105454 0.6198452712843636 0 +0.09319446963215131 0.5455557570572105 0 +0.08390828032651966 0.4712662426121573 0 +0.07462209102084899 0.396976728166792 0 +0.06833943571409715 0.3467154857127772 0 +0.0620567804073453 0.2964542432587625 0 +0.0557741250623971 0.2461930004991768 0 +0.04949146971744887 0.195931757739591 0 +0.04568228763075501 0.1654583010460401 0 +0.04187310554406112 0.134984844352489 0 +0.03806392350447164 0.1045113880357732 0 +0.03425474146488217 0.07403793171905737 0 +0.03194105611161785 0.05552844889294284 0 +0.02962737075835353 0.03701896606682831 0 +0.02731368537917676 0.01850948303341415 0 +1 15 0 11 +103 +104 +105 +106 +107 +108 +109 +110 +111 +112 +113 +-0.025 -0.008320887085287996 0 +-0.025 -0.01664177417057405 0 +-0.025 -0.02496266126225663 0 +-0.025 -0.03328354835393922 0 +-0.025 -0.04161687535162916 0 +-0.025 -0.0499502023493191 0 +-0.025 -0.05828352932485886 0 +-0.025 -0.06661685630040591 0 +-0.025 -0.07496264221220524 0 +-0.025 -0.08330842812401168 0 +-0.025 -0.09165421406200386 0 +1 16 0 11 +114 +115 +116 +117 +118 +119 +120 +121 +122 +123 +124 +0.025 -0.00831257953634283 0 +0.025 -0.01662515907268408 0 +0.025 -0.02493773860952384 0 +0.025 -0.0332503181463636 0 +0.025 -0.04158363382427744 0 +0.025 -0.04991694950219128 0 +0.025 -0.05825026518540654 0 +0.025 -0.06658358086860489 0 +0.025 -0.07493768567063895 0 +0.025 -0.0832917904726839 0 +0.025 -0.09164589523634831 0 +1 17 0 3 +125 +126 +127 +-0.01250000000611455 0 0 +-1.223875861766643e-11 0 0 +0.01249999999388059 0 0 +1 18 0 11 +128 +129 +130 +131 +132 +133 +134 +135 +136 +137 +138 +-0.1 -0.1643542012279202 0 +-0.1 -0.1837084024558237 0 +-0.1 -0.203062603692792 0 +-0.1 -0.2224168049297332 0 +-0.1 -0.2455094337226784 0 +-0.1 -0.2686020625155881 0 +-0.1 -0.291694691317036 0 +-0.1 -0.314787320118585 0 +-0.1 -0.3423404900468721 0 +-0.1 -0.369893659975048 0 +-0.1 -0.3974468299875857 0 +1 19 0 11 +139 +140 +141 +142 +143 +144 +145 +146 +147 +148 +149 +-0.08333553653291947 -0.425 0 +-0.06667107306582606 -0.425 0 +-0.05000660960664088 -0.425 0 +-0.03334214614745569 -0.425 0 +-0.01667547955787649 -0.425 0 +-8.812968297314261e-06 -0.425 0 +0.01665785357667086 -0.425 0 +0.03332452012163059 -0.425 0 +0.04999339008511866 -0.425 0 +0.06666226004857376 -0.425 0 +0.08333113002430761 -0.425 0 +1 20 0 11 +150 +151 +152 +153 +154 +155 +156 +157 +158 +159 +160 +0.1 -0.3971640973436882 0 +0.1 -0.3693281946874231 0 +0.1 -0.3414922920270518 0 +0.1 -0.3136563893666805 0 +0.1 -0.2905961731253152 0 +0.1 -0.2675359568838517 0 +0.1 -0.2444757406716081 0 +0.1 -0.2214155244594352 0 +0.1 -0.2023116433198466 0 +0.1 -0.1832077621802492 0 +0.1 -0.1641038810901278 0 +2 1 0 67 +161 +162 +163 +164 +165 +166 +167 +168 +169 +170 +171 +172 +173 +174 +175 +176 +177 +178 +179 +180 +181 +182 +183 +184 +185 +186 +187 +188 +189 +190 +191 +192 +193 +194 +195 +196 +197 +198 +199 +200 +201 +202 +203 +204 +205 +206 +207 +208 +209 +210 +211 +212 +213 +214 +215 +216 +217 +218 +219 +220 +221 +222 +223 +224 +225 +226 +227 +-0.2132938084975145 1.395555971440597 0 +-0.2392242703847875 1.255313692128066 0 +0.2133642549921056 1.395728409138389 0 +0.2391878054093051 1.255544332934081 0 +4.514633166138177e-05 1.106415459793024 0 +0.1257560444497584 1.114201013355665 0 +-0.1257212985670526 1.114234139275252 0 +0.2600947269247008 1.12635556789028 0 +-0.2601178319250703 1.126251666564875 0 +0.0001883800818643407 1.188336797446429 0 +0.09207534215570556 1.225819478614264 0 +-0.09184390128226624 1.22563384069935 0 +0.0001371903667899556 1.359558547013693 0 +0.000126589062740913 1.462094642782964 0 +0.2051733784802625 1.471812703976191 0 +0.1164404201234327 1.383516838029624 0 +0.1046748626236101 1.465752266272805 0 +0.3073582102759602 1.404219552481517 0 +0.3031623353993226 1.476619561840829 0 +0.2259173154895813 1.323678807192205 0 +0.1619935462544438 1.242949326642699 0 +0.05954761730563889 1.283801327302675 0 +0.1391410567921438 1.30884531331389 0 +0.3192614593354991 1.264737539418085 0 +0.3130899039560679 1.333539790501081 0 +0.03565090731140585 1.204766706679941 0 +-0.03518912807737815 1.20478906026443 0 +-0.05903108078601581 1.283504466098612 0 +0.0002210199231436304 1.246415760101094 0 +0.1126027350037846 1.16995622557921 0 +0.06204045096832572 1.10868588666841 0 +0.0001036108404640036 1.154335444433642 0 +0.05321636048827538 1.159209465216122 0 +0.1362795499599922 1.057373849542891 0 +1.095076268754856e-05 1.054138109376955 0 +0.06855448318884121 1.055058898057339 0 +-0.06194793508693222 1.108841177430768 0 +-0.1124944549144779 1.169940719781085 0 +-0.05301639265544201 1.159234681894389 0 +-0.1362597390806178 1.05739549897737 0 +-0.06852281943779215 1.055097332601764 0 +0.3297102167983371 1.131406077749215 0 +0.2508150795378487 1.190149501191306 0 +0.324919476228049 1.197562126291399 0 +0.2678496627282548 1.06307379286865 0 +0.3337716332574162 1.065620630210797 0 +0.1918718042900419 1.120451473567559 0 +0.1792726043055483 1.180859001473306 0 +0.2023356812385735 1.060222740013594 0 +-0.3192338570428735 1.264526203263012 0 +-0.250846173980518 1.190008407276004 0 +-0.3297184664609963 1.131261795619289 0 +-0.3249290087651571 1.197372510806928 0 +-0.1620155988411709 1.242557413514103 0 +-0.1918718346078768 1.120398128211697 0 +-0.1792658162688962 1.180713765573937 0 +-0.2678496155468384 1.06300987275513 0 +-0.3337744806871998 1.065550549272353 0 +-0.2023307284744054 1.060200893665674 0 +-0.3072653348192199 1.404058100734719 0 +-0.2050443854664297 1.47172815625577 0 +-0.3030735003892254 1.476541242794607 0 +-0.116367977614409 1.383230486002739 0 +-0.1044412754930646 1.465674099135136 0 +-0.2258874403332681 1.323415757784205 0 +-0.313058190782124 1.333328415881823 0 +-0.1390047164781991 1.308578370809387 0 +2 2 0 57 +228 +229 +230 +231 +232 +233 +234 +235 +236 +237 +238 +239 +240 +241 +242 +243 +244 +245 +246 +247 +248 +249 +250 +251 +252 +253 +254 +255 +256 +257 +258 +259 +260 +261 +262 +263 +264 +265 +266 +267 +268 +269 +270 +271 +272 +273 +274 +275 +276 +277 +278 +279 +280 +281 +282 +283 +284 +3.256832099331371e-06 0.8475192167935834 0 +6.712770042346591e-06 0.6963566348017055 0 +7.077370760059237e-06 0.5490811345970505 0 +9.52570470289394e-06 0.4136292889405176 0 +1.013778287606553e-05 0.3021838702367814 0 +8.571720313383345e-06 0.2100386222911862 0 +7.096873312220555e-06 0.1391320977228473 0 +4.19583442297128e-06 0.08249140671291205 0 +2.946349648812055e-06 0.03860923306518894 0 +-0.06543761245599083 0.8473623540266123 0 +1.654156443267434e-06 0.9237034401579128 0 +-0.07021842329084046 0.9236409543611115 0 +0.06544430506470678 0.8473891019661306 0 +0.07022181973117082 0.9236547513475944 0 +-0.05591987872724494 0.6958305234515301 0 +4.697563226589187e-06 0.7716271079741809 0 +-0.06066253212092503 0.7712991176842749 0 +0.0559335738541732 0.6958849194980312 0 +0.06067267025120524 0.7713405343527939 0 +-0.04659472921618596 0.5479819058260884 0 +7.467369359564786e-06 0.6220217864228818 0 +-0.05123681019609492 0.6213283705940706 0 +0.04661168087828301 0.5480618848266661 0 +0.0512530439021824 0.6213988671284707 0 +-0.03780645424673212 0.4106626875582928 0 +4.333329883592692e-06 0.4789253470298152 0 +-0.04202165075286883 0.4767538996219682 0 +0.03782935755874824 0.4107994094730041 0 +0.042041130158506 0.4768594818327732 0 +-0.03102932894016014 0.3004265278621326 0 +1.104142835779264e-05 0.3552863006061386 0 +-0.03423496360446895 0.3527982941222524 0 +0.0310520888197439 0.300548284242142 0 +0.03425978559816961 0.3529172674028761 0 +-0.02515307624625421 0.2074809826672947 0 +6.046998835819179e-06 0.2536916573463482 0 +-0.02794047101197707 0.251539910034499 0 +0.02517307858197061 0.2075954256314517 0 +0.02796134779774429 0.2516473211585039 0 +-0.0209380870047601 0.137868772583783 0 +8.59633653351637e-06 0.1724345123566292 0 +-0.02289591942192308 0.1704242508148285 0 +0.02095369768516925 0.1379612786241512 0 +0.02291489946799052 0.1705151690735918 0 +-0.01737534177095312 0.08097661863552491 0 +3.572084827702885e-06 0.1092563152919942 0 +-0.01906485953927339 0.107899029768432 0 +0.01738530029378433 0.08104098511786952 0 +0.01907719604466421 0.1079651815021989 0 +-0.0148144678705717 0.03812351201225138 0 +3.820263798798426e-06 0.05938185812378294 0 +-0.0160057864003869 0.05828368903177573 0 +0.01482039316775696 0.03816437267843846 0 +0.01601405787458411 0.058325379800794 0 +1.503402345097698e-06 0.01910379546942349 0 +-0.01365531741864644 0.01891716066372254 0 +0.01365833439193384 0.01893265713667808 0 +2 3 0 33 +285 +286 +287 +288 +289 +290 +291 +292 +293 +294 +295 +296 +297 +298 +299 +300 +301 +302 +303 +304 +305 +306 +307 +308 +309 +310 +311 +312 +313 +314 +315 +316 +317 +-2.526607720786112e-12 -0.01663568459705796 0 +-7.17403783948023e-13 -0.03327862042667908 0 +-3.576248455358644e-13 -0.04993806951730116 0 +-6.710474190207982e-13 -0.06661185341451353 0 +-2.394065631503017e-12 -0.08330225084452467 0 +-0.01250000000175519 -0.01663917880164086 0 +-5.009384834026021e-12 -0.008317566790137221 0 +-0.01250000000340767 -0.008319396006339757 0 +0.01249999999827545 -0.01663089261731599 0 +0.0124999999965984 -0.008315244799248306 0 +-0.01250000000049963 -0.03328418146700401 0 +-1.268954667905942e-12 -0.02495550916947142 0 +-0.01250000000088145 -0.02496017871970053 0 +0.01249999999951554 -0.03326897900839347 0 +0.01249999999912794 -0.02494794177597526 0 +-0.01250000000024771 -0.0499450409373201 0 +-4.226826137362916e-13 -0.04160647819804649 0 +-0.0125000000002931 -0.04161292629187771 0 +0.01249999999975709 -0.04992847125985111 0 +0.0124999999997106 -0.04159653427152513 0 +-0.01250000000046526 -0.06661744564421895 0 +-4.087261128507547e-13 -0.05827305444450066 0 +-0.01250000000028127 -0.05827955205780896 0 +0.0124999999995392 -0.06660222027844988 0 +0.01249999999972041 -0.05826314823007122 0 +-0.01250000000166737 -0.08330577697526809 0 +-1.191193194430284e-12 -0.07495544331530243 0 +-0.01250000000082654 -0.07496011347550156 0 +0.01249999999834146 -0.08329750460497537 0 +0.0124999999991751 -0.07494786484582765 0 +-4.912608514429095e-12 -0.09165085207691877 0 +-0.01250000000335577 -0.09165269511977785 0 +0.0124999999966453 -0.09164854220120303 0 +2 4 0 29 +318 +319 +320 +321 +322 +323 +324 +325 +326 +327 +328 +329 +330 +331 +332 +333 +334 +335 +336 +337 +338 +339 +340 +341 +342 +343 +344 +345 +346 +0.02860204362953757 -0.1256020744033201 0 +-0.02795253395252555 -0.1255835251189625 0 +0.0003505614911618749 -0.1236638399733709 0 +0.05279325374344632 -0.1304753167726549 0 +-0.0522965501475279 -0.1304562802748023 0 +0.0001701549109202772 -0.1121182403210541 0 +0.0146266721098999 -0.1241292881943874 0 +0.0274817988165765 -0.1148942092921103 0 +0.01361071790575484 -0.1126044805386362 0 +-0.02715999286636564 -0.1148761714300451 0 +-0.01394211033368743 -0.1241446086139232 0 +-0.01327397785564842 -0.1126030165835662 0 +-0.02988252183387667 -0.1354263933118131 0 +0.0005520554603648254 -0.1345429501920538 0 +-0.01480501303631099 -0.1347797710330038 0 +0.03090013562483596 -0.1354393241819223 0 +0.01589548410939002 -0.1347808223699551 0 +0.04045345777952157 -0.1279136090471275 0 +0.05747703111939562 -0.1266406111539971 0 +0.04225827901905877 -0.1199782750297362 0 +0.05859712230980578 -0.1366199417504551 0 +0.04530083124876786 -0.136333097038436 0 +0.07034133785011429 -0.1355862883060414 0 +-0.03986304401693938 -0.1278469416628678 0 +-0.05801230907990203 -0.1365897218272945 0 +-0.04445615831442058 -0.136310919526907 0 +-0.05732454250894184 -0.126528990088689 0 +-0.04198681456452053 -0.1199406529344008 0 +-0.07005136189180408 -0.1355597396399805 0 +2 5 0 121 +347 +348 +349 +350 +351 +352 +353 +354 +355 +356 +357 +358 +359 +360 +361 +362 +363 +364 +365 +366 +367 +368 +369 +370 +371 +372 +373 +374 +375 +376 +377 +378 +379 +380 +381 +382 +383 +384 +385 +386 +387 +388 +389 +390 +391 +392 +393 +394 +395 +396 +397 +398 +399 +400 +401 +402 +403 +404 +405 +406 +407 +408 +409 +410 +411 +412 +413 +414 +415 +416 +417 +418 +419 +420 +421 +422 +423 +424 +425 +426 +427 +428 +429 +430 +431 +432 +433 +434 +435 +436 +437 +438 +439 +440 +441 +442 +443 +444 +445 +446 +447 +448 +449 +450 +451 +452 +453 +454 +455 +456 +457 +458 +459 +460 +461 +462 +463 +464 +465 +466 +467 +-0.06642393152214293 -0.1861238003696422 0 +-0.03290723459840507 -0.1875589617999193 0 +0.0004826226830231394 -0.1879576404074532 0 +0.03375461319670166 -0.1873955668864127 0 +0.06690673544544402 -0.1858001428419833 0 +-0.06652215346426048 -0.2285306606292938 0 +-0.03308248658343373 -0.2311519638676335 0 +0.0002873703681992572 -0.2318663365887817 0 +0.03358185151916766 -0.2308784392631879 0 +0.06680995636737984 -0.2279794182688292 0 +-0.06658331889165471 -0.2736722342067517 0 +-0.03318898500114319 -0.2766638220772877 0 +0.0001659717612850591 -0.2775007073925811 0 +0.03347659538458531 -0.2763361585509644 0 +0.06674946622505423 -0.2730129654102407 0 +-0.06662216009791527 -0.3215757991775318 0 +-0.03325619574945739 -0.3243768838895322 0 +8.893664065979254e-05 -0.3251460722146645 0 +0.03340998178539659 -0.3240961907696367 0 +0.06671086717267316 -0.3209792034190805 0 +-0.06664906900559732 -0.3725585477316471 0 +-0.03330289989664165 -0.3741229578801423 0 +3.523716471874128e-05 -0.37457008757634 0 +0.03336372839567159 -0.3739830790809739 0 +0.06668423613577446 -0.3722378819026758 0 +-0.0832092834975387 -0.1850039043327792 0 +-0.0663576250609987 -0.1655008876009508 0 +-0.08317701652345862 -0.1649624754461848 0 +-0.0496528153850609 -0.1869787490581385 0 +-0.03277735723241079 -0.1662158467977762 0 +-0.04955009033420098 -0.1659231199186373 0 +-0.0161969704002822 -0.1878807050525225 0 +0.0006154271758174905 -0.1664180134377203 0 +-0.01606161756346906 -0.1663789385371946 0 +0.01713363883334696 -0.1877973900116069 0 +0.03388514037466567 -0.1661318300232611 0 +0.0172693724479823 -0.166336651482423 0 +0.05034301033706197 -0.1867392062807639 0 +0.0669732195794167 -0.1653367632998641 0 +0.05044643962308279 -0.1658005905121937 0 +0.08345621872757433 -0.1845939571282734 0 +0.08348857724001332 -0.1647561389429657 0 +-0.08325876745115746 -0.226274587165453 0 +-0.06647874538780561 -0.2070593128293657 0 +-0.08323664697240381 -0.2052897772140204 0 +-0.04979515852894207 -0.2301186336124149 0 +-0.03300648477381549 -0.209163930393109 0 +-0.049733228727178 -0.2083285313288243 0 +-0.01638818233011 -0.2317302236076177 0 +0.000373544517714486 -0.2097390974866308 0 +-0.01630435164291948 -0.2096284998068811 0 +0.01694367588961568 -0.2315900557485716 0 +0.03365692436471934 -0.2089385883407913 0 +0.01702697056641164 -0.2095153052582375 0 +0.05020279173497478 -0.2297144753596893 0 +0.06685258619106546 -0.2066011137444083 0 +0.05026384960664065 -0.2079966174033169 0 +0.0834073052477963 -0.2255454056040742 0 +0.08342908868233873 -0.2046987063645196 0 +-0.08329021955512059 -0.2713263275487437 0 +-0.06655639799132125 -0.250766145972198 0 +-0.08327634528896305 -0.2484058685560031 0 +-0.04988212206144 -0.2754588175606963 0 +-0.03314233858762588 -0.2536449036605182 0 +-0.04984400572944039 -0.2524932481523223 0 +-0.01650604886350525 -0.277338422208911 0 +0.0002193526510274899 -0.2544421361405412 0 +-0.01645428924046528 -0.2542877140221463 0 +0.01682658238501596 -0.277172712617925 0 +0.03352287190221952 -0.2533405139783524 0 +0.01687810633004347 -0.2541340294085523 0 +0.0501169583836297 -0.2749736521647163 0 +0.06677614569361254 -0.2501515292734963 0 +0.05015471535783724 -0.2520417113234097 0 +0.08337615506088593 -0.2704721297227388 0 +0.08338992794247208 -0.2475945337216948 0 +-0.08331029908575791 -0.3191245754527317 0 +-0.06660480865644493 -0.2972566263543129 0 +-0.08330133157346908 -0.2947761582873036 0 +-0.04993706675921018 -0.3232778052081231 0 +-0.03322619375942255 -0.3002489617089821 0 +-0.04991255048311639 -0.299056440683734 0 +-0.01658070471070104 -0.3249955214048255 0 +0.0001233563068653257 -0.3010774530672224 0 +-0.01654737877345012 -0.3009148657497045 0 +0.01675229809757817 -0.3248550420108436 0 +0.03343971758538373 -0.2999279142191809 0 +0.01678549083825439 -0.3007569306549813 0 +0.05006249962355318 -0.3228509140321704 0 +0.06672805575495978 -0.2965944452490156 0 +0.05008679129703973 -0.2985805306254206 0 +0.08335621271896083 -0.3183277314907139 0 +0.08336511220205239 -0.2939170317074145 0 +-0.08332420955187217 -0.371326481098081 0 +-0.0666366198631131 -0.3467656697519828 0 +-0.08331778312310247 -0.3448170184596007 0 +-0.04997516611560808 -0.3734946475790252 0 +-0.0332812755407447 -0.3490443628015074 0 +-0.04995753091774554 -0.3481427218182548 0 +-0.01663267929637692 -0.374477862517811 0 +6.008743961256674e-05 -0.3496768234198649 0 +-0.01660863192827605 -0.3495481810635205 0 +0.01670061359614478 -0.374408832695171 0 +0.03338519506747753 -0.3488319296594218 0 +0.01672456294589746 -0.3494429314941181 0 +0.0500247815493562 -0.3732696557488637 0 +0.06669657805261363 -0.3462946148151626 0 +0.05004227074996934 -0.347809690928893 0 +0.08334242814664146 -0.3708895886760755 0 +0.08334880687563082 -0.3441733620202669 0 +-0.06666035053230808 -0.398710813384908 0 +-0.0833300252185426 -0.3981182382466381 0 +-0.03332262801447777 -0.3994925639522371 0 +-0.04999122269648745 -0.3991730507908433 0 +1.265304160087189e-05 -0.3997211026611473 0 +-0.01665456754387353 -0.3996726982051931 0 +0.03334403259535038 -0.3994217340295844 0 +0.01667875871795278 -0.3996382217398972 0 +0.06667297730066417 -0.39854796715408 0 +0.05000876698186856 -0.3990592264161568 0 +0.08333663529472293 -0.3978975382456735 0 +$EndNodes +$Elements +19 536 1 536 +1 3 1 4 +1 3 15 +2 15 16 +3 16 17 +4 17 1 +1 4 1 4 +5 4 18 +6 18 19 +7 19 20 +8 20 2 +1 5 1 8 +9 11 21 +10 21 22 +11 22 23 +12 23 24 +13 24 25 +14 25 26 +15 26 27 +16 27 9 +1 6 1 4 +17 7 28 +18 28 29 +19 29 30 +20 30 9 +1 7 1 8 +21 11 31 +22 31 32 +23 32 33 +24 33 34 +25 34 35 +26 35 36 +27 36 37 +28 37 12 +1 8 1 8 +29 12 38 +30 38 39 +31 39 40 +32 40 41 +33 41 42 +34 42 43 +35 43 44 +36 44 10 +1 9 1 4 +37 10 45 +38 45 46 +39 46 47 +40 47 8 +1 13 1 20 +41 7 65 +42 65 66 +43 66 67 +44 67 68 +45 68 69 +46 69 70 +47 70 71 +48 71 72 +49 72 73 +50 73 74 +51 74 75 +52 75 76 +53 76 77 +54 77 78 +55 78 79 +56 79 80 +57 80 81 +58 81 82 +59 82 83 +60 83 5 +1 14 1 20 +61 8 84 +62 84 85 +63 85 86 +64 86 87 +65 87 88 +66 88 89 +67 89 90 +68 90 91 +69 91 92 +70 92 93 +71 93 94 +72 94 95 +73 95 96 +74 96 97 +75 97 98 +76 98 99 +77 99 100 +78 100 101 +79 101 102 +80 102 6 +1 15 1 12 +81 5 103 +82 103 104 +83 104 105 +84 105 106 +85 106 107 +86 107 108 +87 108 109 +88 109 110 +89 110 111 +90 111 112 +91 112 113 +92 113 3 +1 16 1 12 +93 6 114 +94 114 115 +95 115 116 +96 116 117 +97 117 118 +98 118 119 +99 119 120 +100 120 121 +101 121 122 +102 122 123 +103 123 124 +104 124 4 +1 18 1 12 +105 1 128 +106 128 129 +107 129 130 +108 130 131 +109 131 132 +110 132 133 +111 133 134 +112 134 135 +113 135 136 +114 136 137 +115 137 138 +116 138 13 +1 19 1 12 +117 13 139 +118 139 140 +119 140 141 +120 141 142 +121 142 143 +122 143 144 +123 144 145 +124 145 146 +125 146 147 +126 147 148 +127 148 149 +128 149 14 +1 20 1 12 +129 14 150 +130 150 151 +131 151 152 +132 152 153 +133 153 154 +134 154 155 +135 155 156 +136 156 157 +137 157 158 +138 158 159 +139 159 160 +140 160 2 +2 1 3 84 +141 173 174 177 176 +142 174 34 35 177 +143 177 35 36 175 +144 176 177 175 163 +145 163 175 179 178 +146 175 36 37 179 +147 179 37 12 38 +148 178 179 38 39 +149 173 176 183 182 +150 176 163 180 183 +151 183 180 164 181 +152 182 183 181 171 +153 163 178 185 180 +154 178 39 40 185 +155 185 40 41 184 +156 180 185 184 164 +157 173 182 189 188 +158 182 171 186 189 +159 189 186 170 187 +160 188 189 187 172 +161 171 190 193 186 +162 190 166 191 193 +163 193 191 165 192 +164 186 193 192 170 +165 166 194 196 191 +166 194 8 50 196 +167 196 50 49 195 +168 191 196 195 165 +169 172 187 199 198 +170 187 170 192 199 +171 199 192 165 197 +172 198 199 197 167 +173 167 197 201 200 +174 197 165 195 201 +175 201 195 49 48 +176 200 201 48 7 +177 41 42 204 184 +178 42 43 202 204 +179 204 202 168 203 +180 184 204 203 164 +181 43 44 206 202 +182 44 10 45 206 +183 206 45 46 205 +184 202 206 205 168 +185 164 203 208 181 +186 203 168 207 208 +187 208 207 166 190 +188 181 208 190 171 +189 168 205 209 207 +190 205 46 47 209 +191 209 47 8 194 +192 207 209 194 166 +193 24 210 213 25 +194 210 162 211 213 +195 213 211 169 212 +196 25 213 212 26 +197 162 214 216 211 +198 214 172 198 216 +199 216 198 167 215 +200 211 216 215 169 +201 26 212 218 27 +202 212 169 217 218 +203 218 217 29 30 +204 27 218 30 9 +205 169 215 219 217 +206 215 167 200 219 +207 219 200 7 28 +208 217 219 28 29 +209 161 220 222 221 +210 220 22 21 222 +211 222 21 11 31 +212 221 222 31 32 +213 173 223 224 174 +214 223 161 221 224 +215 224 221 32 33 +216 174 224 33 34 +217 161 225 226 220 +218 225 162 210 226 +219 226 210 24 23 +220 220 226 23 22 +221 173 188 227 223 +222 188 172 214 227 +223 227 214 162 225 +224 223 227 225 161 +2 2 3 80 +225 7 65 239 48 +226 65 66 237 239 +227 239 237 228 238 +228 48 239 238 49 +229 49 238 241 50 +230 238 228 240 241 +231 241 240 85 84 +232 50 241 84 8 +233 66 67 244 237 +234 67 68 242 244 +235 244 242 229 243 +236 237 244 243 228 +237 228 243 246 240 +238 243 229 245 246 +239 246 245 87 86 +240 240 246 86 85 +241 68 69 249 242 +242 69 70 247 249 +243 249 247 230 248 +244 242 249 248 229 +245 229 248 251 245 +246 248 230 250 251 +247 251 250 89 88 +248 245 251 88 87 +249 70 71 254 247 +250 71 72 252 254 +251 254 252 231 253 +252 247 254 253 230 +253 230 253 256 250 +254 253 231 255 256 +255 256 255 91 90 +256 250 256 90 89 +257 72 73 259 252 +258 73 74 257 259 +259 259 257 232 258 +260 252 259 258 231 +261 231 258 261 255 +262 258 232 260 261 +263 261 260 93 92 +264 255 261 92 91 +265 74 75 264 257 +266 75 76 262 264 +267 264 262 233 263 +268 257 264 263 232 +269 232 263 266 260 +270 263 233 265 266 +271 266 265 95 94 +272 260 266 94 93 +273 76 77 269 262 +274 77 78 267 269 +275 269 267 234 268 +276 262 269 268 233 +277 233 268 271 265 +278 268 234 270 271 +279 271 270 97 96 +280 265 271 96 95 +281 78 79 274 267 +282 79 80 272 274 +283 274 272 235 273 +284 267 274 273 234 +285 234 273 276 270 +286 273 235 275 276 +287 276 275 99 98 +288 270 276 98 97 +289 80 81 279 272 +290 81 82 277 279 +291 279 277 236 278 +292 272 279 278 235 +293 235 278 281 275 +294 278 236 280 281 +295 281 280 101 100 +296 275 281 100 99 +297 82 83 283 277 +298 83 5 125 283 +299 283 125 126 282 +300 277 283 282 236 +301 236 282 284 280 +302 282 126 127 284 +303 284 127 6 102 +304 280 284 102 101 +2 3 3 48 +305 5 103 292 125 +306 103 104 290 292 +307 292 290 285 291 +308 125 292 291 126 +309 126 291 294 127 +310 291 285 293 294 +311 294 293 115 114 +312 127 294 114 6 +313 104 105 297 290 +314 105 106 295 297 +315 297 295 286 296 +316 290 297 296 285 +317 285 296 299 293 +318 296 286 298 299 +319 299 298 117 116 +320 293 299 116 115 +321 106 107 302 295 +322 107 108 300 302 +323 302 300 287 301 +324 295 302 301 286 +325 286 301 304 298 +326 301 287 303 304 +327 304 303 119 118 +328 298 304 118 117 +329 108 109 307 300 +330 109 110 305 307 +331 307 305 288 306 +332 300 307 306 287 +333 287 306 309 303 +334 306 288 308 309 +335 309 308 121 120 +336 303 309 120 119 +337 110 111 312 305 +338 111 112 310 312 +339 312 310 289 311 +340 305 312 311 288 +341 288 311 314 308 +342 311 289 313 314 +343 314 313 123 122 +344 308 314 122 121 +345 112 113 316 310 +346 113 3 51 316 +347 316 51 52 315 +348 310 316 315 289 +349 289 315 317 313 +350 315 52 53 317 +351 317 53 4 124 +352 313 317 124 123 +2 4 3 40 +353 4 53 326 325 +354 53 52 323 326 +355 326 323 320 324 +356 325 326 324 318 +357 52 51 329 323 +358 51 3 327 329 +359 329 327 319 328 +360 323 329 328 320 +361 320 328 332 331 +362 328 319 330 332 +363 332 330 57 58 +364 331 332 58 59 +365 318 324 334 333 +366 324 320 331 334 +367 334 331 59 60 +368 333 334 60 61 +369 4 325 337 18 +370 325 318 335 337 +371 337 335 321 336 +372 18 337 336 19 +373 318 333 339 335 +374 333 61 62 339 +375 339 62 63 338 +376 335 339 338 321 +377 19 336 340 20 +378 336 321 338 340 +379 340 338 63 64 +380 20 340 64 2 +381 57 330 343 56 +382 330 319 341 343 +383 343 341 322 342 +384 56 343 342 55 +385 319 327 345 341 +386 327 3 15 345 +387 345 15 16 344 +388 341 345 344 322 +389 55 342 346 54 +390 342 322 344 346 +391 346 344 16 17 +392 54 346 17 1 +2 5 3 144 +393 1 128 374 54 +394 128 129 372 374 +395 374 372 347 373 +396 54 374 373 55 +397 55 373 377 56 +398 373 347 375 377 +399 377 375 348 376 +400 56 377 376 57 +401 57 376 380 58 +402 376 348 378 380 +403 380 378 349 379 +404 58 380 379 59 +405 59 379 383 60 +406 379 349 381 383 +407 383 381 350 382 +408 60 383 382 61 +409 61 382 386 62 +410 382 350 384 386 +411 386 384 351 385 +412 62 386 385 63 +413 63 385 388 64 +414 385 351 387 388 +415 388 387 159 160 +416 64 388 160 2 +417 129 130 391 372 +418 130 131 389 391 +419 391 389 352 390 +420 372 391 390 347 +421 347 390 394 375 +422 390 352 392 394 +423 394 392 353 393 +424 375 394 393 348 +425 348 393 397 378 +426 393 353 395 397 +427 397 395 354 396 +428 378 397 396 349 +429 349 396 400 381 +430 396 354 398 400 +431 400 398 355 399 +432 381 400 399 350 +433 350 399 403 384 +434 399 355 401 403 +435 403 401 356 402 +436 384 403 402 351 +437 351 402 405 387 +438 402 356 404 405 +439 405 404 157 158 +440 387 405 158 159 +441 131 132 408 389 +442 132 133 406 408 +443 408 406 357 407 +444 389 408 407 352 +445 352 407 411 392 +446 407 357 409 411 +447 411 409 358 410 +448 392 411 410 353 +449 353 410 414 395 +450 410 358 412 414 +451 414 412 359 413 +452 395 414 413 354 +453 354 413 417 398 +454 413 359 415 417 +455 417 415 360 416 +456 398 417 416 355 +457 355 416 420 401 +458 416 360 418 420 +459 420 418 361 419 +460 401 420 419 356 +461 356 419 422 404 +462 419 361 421 422 +463 422 421 155 156 +464 404 422 156 157 +465 133 134 425 406 +466 134 135 423 425 +467 425 423 362 424 +468 406 425 424 357 +469 357 424 428 409 +470 424 362 426 428 +471 428 426 363 427 +472 409 428 427 358 +473 358 427 431 412 +474 427 363 429 431 +475 431 429 364 430 +476 412 431 430 359 +477 359 430 434 415 +478 430 364 432 434 +479 434 432 365 433 +480 415 434 433 360 +481 360 433 437 418 +482 433 365 435 437 +483 437 435 366 436 +484 418 437 436 361 +485 361 436 439 421 +486 436 366 438 439 +487 439 438 153 154 +488 421 439 154 155 +489 135 136 442 423 +490 136 137 440 442 +491 442 440 367 441 +492 423 442 441 362 +493 362 441 445 426 +494 441 367 443 445 +495 445 443 368 444 +496 426 445 444 363 +497 363 444 448 429 +498 444 368 446 448 +499 448 446 369 447 +500 429 448 447 364 +501 364 447 451 432 +502 447 369 449 451 +503 451 449 370 450 +504 432 451 450 365 +505 365 450 454 435 +506 450 370 452 454 +507 454 452 371 453 +508 435 454 453 366 +509 366 453 456 438 +510 453 371 455 456 +511 456 455 151 152 +512 438 456 152 153 +513 137 138 458 440 +514 138 13 139 458 +515 458 139 140 457 +516 440 458 457 367 +517 367 457 460 443 +518 457 140 141 460 +519 460 141 142 459 +520 443 460 459 368 +521 368 459 462 446 +522 459 142 143 462 +523 462 143 144 461 +524 446 462 461 369 +525 369 461 464 449 +526 461 144 145 464 +527 464 145 146 463 +528 449 464 463 370 +529 370 463 466 452 +530 463 146 147 466 +531 466 147 148 465 +532 452 466 465 371 +533 371 465 467 455 +534 465 148 149 467 +535 467 149 14 150 +536 455 467 150 151 +$EndElements diff --git a/prm/todo/ideal-blast.prm b/prm/todo/ideal-blast.prm index 192035980..2bf11310e 100644 --- a/prm/todo/ideal-blast.prm +++ b/prm/todo/ideal-blast.prm @@ -54,7 +54,6 @@ subsection E - InitialValues end subsection F - HyperbolicModule - set cfl with boundary dofs = false set limiter iterations = 2 set limiter newton max iterations = 2 set limiter newton tolerance = 1e-10