Skip to content

Commit 5a1aa3c

Browse files
committed
updated comparison page and pngs
1 parent ac57eea commit 5a1aa3c

4 files changed

Lines changed: 68 additions & 12 deletions

File tree

docs/source/compressible_compare.rst

Lines changed: 68 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -3,11 +3,12 @@ Compressible solver comparisons
33

44
We run various problems run with the different compressible solvers in pyro (standard Riemann, Runge-Kutta, fourth order).
55

6+
67
Kelvin-Helmholtz
78
^^^^^^^^^^^^^^^^
89
The McNally Kelvin-Helmholtz problem sets up a heavier fluid moving in the negative x-direction sandwiched between regions of lighter fluid moving in the positive x-direction.
910

10-
The image below shows the KH problem initialized with McNally's test. It had a resolution of 128 in the x and y-directions, gamma = 1.7, and ran until t = 2.0. This is run with:
11+
The image below shows the KH problem initialized with McNally's test. It ran on a 128 x 128 grid, with gamma = 1.7, and ran until t = 2.0. This is run with:
1112

1213
.. code-block:: none
1314
@@ -16,28 +17,83 @@ The image below shows the KH problem initialized with McNally's test. It had a r
1617
./pyro.py compressible_fv4 kh inputs.kh kh.vbulk=0
1718
1819
.. image:: ./solver_comparisons/kh.png
20+
:align: center
21+
1922

2023
We vary the velocity in the positive y-direction (vbulk) to see how effective the solvers are at preserving the initial shape.
2124

25+
2226
Sedov
2327
^^^^^
24-
25-
The Sedov problem ran with a resolution of 128 in the x and y-directions, gamma = 1.4, and until t = 0.1, which can be run as:
28+
The Sedov problem ran on a 128 x 128 grid, with gamma = 1.4, and until t = 0.1, which can be run as:
2629

2730
.. code-block:: none
2831
29-
./pyro.py compressible sedov inputs.sedov
30-
31-
.. image:: ./solver_comparisons/sedov_standard.png
32+
./pyro.py compressible sedov inputs.sedov
33+
./pyro.py compressible_rk sedov inputs.sedov
34+
./pyro.py compressible_fv4 sedov inputs.sedov
3235
33-
.. code-block:: none
34-
35-
./pyro.py compressible_rk sedov inputs.sedov
36+
.. image:: ./solver_comparisons/sedov.png
3637

3738
.. image:: ./solver_comparisons/sedov_rk.png
3839

40+
.. image:: ./solver_comparisons/sedov_fv4.png
41+
42+
Quad
43+
^^^^
44+
The quad problem ran on a 256 x 256 grid until t = 0.8, which can be run as:
45+
3946
.. code-block:: none
40-
41-
./pyro.py compressible_fv4 sedov inputs.sedov
4247
43-
.. image:: ./solver_comparisons/sedov_fourth.png
48+
./pyro.py compressible quad inputs.quad
49+
./pyro.py compressible_rk quad inputs.quad
50+
./pyro.py compressible_fv4 quad inputs.quad
51+
52+
.. image:: ./solver_comparisons/quad.png
53+
:width: 30%
54+
55+
.. image:: ./solver_comparisons/quad_rk.png
56+
:width: 30%
57+
58+
.. image:: ./solver_comparisons/quad_fv4.png
59+
:width: 30%
60+
61+
62+
Bubble
63+
^^^^^^
64+
The bubble problem ran on a 128 x 256 grid until t = 3.0, which can be run as:
65+
66+
.. code-block:: none
67+
68+
./pyro.py compressible bubble inputs.bubble
69+
./pyro.py compressible_rk bubble inputs.bubble
70+
./pyro.py compressible_fv4 bubble inputs.bubble
71+
72+
.. image:: ./solver_comparisons/bubble.png
73+
:width: 30%
74+
75+
.. image:: ./solver_comparisons/bubble_rk.png
76+
:width: 30%
77+
78+
.. image:: ./solver_comparisons/bubble_fv4.png
79+
:width: 30%
80+
81+
82+
Rayleigh-Taylor
83+
^^^^^^^^^^^^^^^
84+
The Rayleigh-Taylor problem ran on a 64 x 192 grid until t = 3.0, which can be run as:
85+
86+
.. code-block:: none
87+
88+
./pyro.py compressible rt inputs.rt
89+
./pyro.py compressible_rk rt inputs.rt
90+
./pyro.py compressible_fv4 rt inputs.rt
91+
92+
.. image:: ./solver_comparisons/rt.png
93+
:align: center
94+
95+
.. image:: ./solver_comparisons/rt_rk.png
96+
:align: center
97+
98+
.. image:: ./solver_comparisons/rt_fv4.png
99+
:align: center
-57.6 KB
Binary file not shown.
-11.6 KB
Loading
-56.9 KB
Binary file not shown.

0 commit comments

Comments
 (0)