@@ -15,10 +15,11 @@ The image below shows the KH problem initialized with McNally's test. It ran on
1515 ./pyro.py compressible kh inputs.kh kh.vbulk=0
1616 ./pyro.py compressible_rk kh inputs.kh kh.vbulk=0
1717 ./pyro.py compressible_fv4 kh inputs.kh kh.vbulk=0
18+ ./pyro.py compressible_sdc kh inputs.kh kh.vbulk=0
1819
1920 .. image :: ./solver_comparisons/kh.png
2021 :align: center
21-
22+
2223
2324We vary the velocity in the positive y-direction (vbulk) to see how effective the solvers are at preserving the initial shape.
2425
@@ -28,20 +29,24 @@ Sedov
2829The Sedov problem ran on a 128 x 128 grid, with gamma = 1.4, and until t = 0.1, which can be run as:
2930
3031.. code-block :: none
31-
32+
3233 ./pyro.py compressible sedov inputs.sedov
3334 ./pyro.py compressible_rk sedov inputs.sedov
3435 ./pyro.py compressible_fv4 sedov inputs.sedov
36+ ./pyro.py compressible_sdc sedov inputs.sedov
3537
3638 .. image :: ./solver_comparisons/sedov.png
3739 :align: center
38-
40+
3941.. image :: ./solver_comparisons/sedov_rk.png
4042 :align: center
41-
43+
4244.. image :: ./solver_comparisons/sedov_fv4.png
4345 :align: center
44-
46+
47+ .. image :: ./solver_comparisons/sedov_sdc.png
48+ :align: center
49+
4550Quad
4651^^^^
4752The quad problem ran on a 256 x 256 grid until t = 0.8, which can be run as:
@@ -51,17 +56,21 @@ The quad problem ran on a 256 x 256 grid until t = 0.8, which can be run as:
5156 ./pyro.py compressible quad inputs.quad
5257 ./pyro.py compressible_rk quad inputs.quad
5358 ./pyro.py compressible_fv4 quad inputs.quad
59+ ./pyro.py compressible_sdc quad inputs.quad
5460
5561 .. image :: ./solver_comparisons/quad.png
5662 :align: center
57-
63+
5864.. image :: ./solver_comparisons/quad_rk.png
5965 :align: center
6066
6167.. image :: ./solver_comparisons/quad_fv4.png
6268 :align: center
63-
64-
69+
70+ .. image :: ./solver_comparisons/quad_sdc.png
71+ :align: center
72+
73+
6574Bubble
6675^^^^^^
6776The bubble problem ran on a 128 x 256 grid until t = 3.0, which can be run as:
@@ -71,6 +80,7 @@ The bubble problem ran on a 128 x 256 grid until t = 3.0, which can be run as:
7180 ./pyro.py compressible bubble inputs.bubble
7281 ./pyro.py compressible_rk bubble inputs.bubble
7382 ./pyro.py compressible_fv4 bubble inputs.bubble
83+ ./pyro.py compressible_sdc bubble inputs.bubble
7484
7585 .. image :: ./solver_comparisons/bubble.png
7686 :align: center
@@ -81,6 +91,9 @@ The bubble problem ran on a 128 x 256 grid until t = 3.0, which can be run as:
8191.. image :: ./solver_comparisons/bubble_fv4.png
8292 :align: center
8393
94+ .. image :: ./solver_comparisons/bubble_sdc.png
95+ :align: center
96+
8497
8598Rayleigh-Taylor
8699^^^^^^^^^^^^^^^
@@ -91,12 +104,16 @@ The Rayleigh-Taylor problem ran on a 64 x 192 grid until t = 3.0, which can be r
91104 ./pyro.py compressible rt inputs.rt
92105 ./pyro.py compressible_rk rt inputs.rt
93106 ./pyro.py compressible_fv4 rt inputs.rt
107+ ./pyro.py compressible_sdc rt inputs.rt
94108
95109 .. image :: ./solver_comparisons/rt.png
96110 :align: center
97111
98112.. image :: ./solver_comparisons/rt_rk.png
99113 :align: center
100-
114+
101115.. image :: ./solver_comparisons/rt_fv4.png
102116 :align: center
117+
118+ .. image :: ./solver_comparisons/rt_sdc.png
119+ :align: center
0 commit comments