Differences

This shows you the differences between two versions of the page.

--- dustdevils [2025/04/01 08:02] – [vertical vorticity at the lowest level] admin
+++ dustdevils [2025/04/12 19:29] (current) – admin
@@ Line 1: / Line 1: @@
-===== Dustdevils redux =====
+===== Dust devils redux: using numerical python =====
-This page shows a revisit to {{ :dustdevils:fiedlerkanak2001.pdf | Rayleigh-BeÂnard convection as a tool for studying dust devils}}, but with the computations done with numerical python.
+This page shows a revisit to {{ :dustdevils:fiedlerkanak2001.pdf | Rayleigh-Benard convection as a tool for studying dust devils}}, but with the computations done with numerical python.
+It is a substantial upgrade to [[https://12characters.net/dustdevil/ |my old dust devil page]].
 To reproduce these simulations, you may use my ipython notebook:
-[[https://github.com/bfiedler/metr4323/blob/main/N140_convection3d.ipynb|3D convection, periodic]]
+[[https://github.com/bfiedler/metr4323/blob/main/N140_convection3d.ipynb|3D convection, periodic]].
-In cell number 28, set ''nexp'' to 91:
+To reproduce the visualizations, you will need {{https://github.com/bfiedler/dustdevils/blob/main/dustdevils.tar.gz | my companion python programs}}.
+In the notebook, in cell number 28, set ''nexp'' to 91:
 <code>
 # choose an integer, see choices for nexp below
@@ Line 48: / Line 52: @@
 Try
 <code>
-./trace.py R4096g0w4 1 400
+./hist.py R4096g0w4 1 400
 </code>
 if all went okay,
 <code>
-./trace.py R4096g0w4 1 400 > R4096g0w4.hist
+./hist.py R4096g0w4 1 40 > R4096g0w4.hist
-</code>
+</code>0
 And then
 <code>
-./plottrace.py -s R4096g0w4.hist
+./plothist.py -s R4096g0w4.hist
 </code>
 That should produce the following image, which shows a dust devil event just prior to t=30.
@@ Line 71: / Line 74: @@
 <code>
 mkdir pngs
-./hcontour.py -o pngs R4096g0w4 000001 000400
+./npy2contour.py -o pngs R4096g0w4 000001 000400
 ./maketitle.py -o titleh.png
 ./png2mp4.py pngs/0 -t titleh.png
 </code>
-That will produce ''video.mp4'', showing twin dust devils peaking at t=29.9:
+That will produce ''video.mp4'', showing twin dust devils peaking at t=29.8:
 {{ :dustdevils:wzeta.mp4 |}}
@@ Line 82: / Line 85: @@
 ==== isosurface animation ====
-The extreme event (a.k.a. dustdevil).
+The extreme event (a.k.a. dust devils) in 3-D, with isosurfaces.
+<code>
+mkdir vista
+./npy2vista.py  R4096g0w4 000270 000320 -o vista
+./maketitlevista.py
+./png2mp4.py vista/0 -t titlevista.png
+</code>
 {{ :dustdevils:pzeta.mp4 |}}
 ==== Some interactive images at $t=29.8$ ====
+<code>
+./npy2vista.py -s R4096g0w4 000298
+open scene.html
+</code>
 Use your mouse, have fun ...
 {{url>https://12characters.net/rbdd/t298zm40pm2.html}}
 $\zeta=-40$ (purple)   , $p=-2$ (green), and $s$ at the surface
+I did not finish implementing the command line options for ''npy2vista.py'',
+but the source code can be easily edited to select the desired isosurfaces.
+(You may want to copy ''npy2vista.py'' to a file of some other name, and edit that...)
+Change these lines:
+<code>
+#    plotter.add_mesh(contours_w_neg, color="lightblue", opacity=0.5, label="w = -0.5")
+#    plotter.add_mesh(contours_w_pos, color="green", opacity=0.5, label="w = 0.5")
+    plotter.add_mesh(contours_p1, color="green", opacity=0.5, label="p = -2")
+#    plotter.add_mesh(contours_p2, color="purple", opacity=1.0, label="p = -3")
+#    plotter.add_mesh(contours_p3, color="yellow", opacity=.5, label="s =  2")
+#    plotter.add_mesh(contours_p4, color="red", opacity=.5, label="c =  6")
+    plotter.add_mesh(contours_p5, color="purple", opacity=.5, label="$\zeta$ =  -40")
+    plotter.add_mesh(edges, color='black', line_width=3, label="Domain Boundaries")
+    plotter.add_mesh(origin_marker, color="yellow", label="Origin Marker")
+    plotter.add_text(time,color='pink')
+#    plotter.add_mesh(contour_2d, cmap="coolwarm", opacity=0.8, label="2D Contour (p at z=0)")
+    plotter.add_mesh(slice_plane, scalars="s", cmap="OrRd", opacity=1.0, clim=[0,3.0 ])
+</code>
+to
+<code>
+    plotter.add_mesh(contours_w_neg, color="lightblue", opacity=0.5, label="w = -0.5")
+    plotter.add_mesh(contours_w_pos, color="green", opacity=0.5, label="w = 0.5")
+#    plotter.add_mesh(contours_p1, color="green", opacity=0.5, label="p = -2")
+#    plotter.add_mesh(contours_p2, color="purple", opacity=1.0, label="p = -3")
+#    plotter.add_mesh(contours_p3, color="yellow", opacity=.5, label="s =  2")
+    plotter.add_mesh(contours_p4, color="red", opacity=.5, label="c =  6")
+    plotter.add_mesh(contours_p5, color="purple", opacity=.5, label="$\zeta$ =  -40")
+    plotter.add_mesh(edges, color='black', line_width=3, label="Domain Boundaries")
+    plotter.add_mesh(origin_marker, color="yellow", label="Origin Marker")
+    plotter.add_text(time,color='pink')
+#    plotter.add_mesh(contour_2d, cmap="coolwarm", opacity=0.8, label="2D Contour (p at z=0)")
+    plotter.add_mesh(slice_plane, scalars="s", cmap="OrRd", opacity=1.0, clim=[0,3.0 ])
+</code>
+And again:
+<code>
+./npy2vista.py -s R4096g0w4 000298
+open scene.html
+</code>
 {{url>https://12characters.net/rbdd/t298wpm5c6.html}}
 $w=.5$ (green), $w=-.5$ (grey) , $c=6$ (red)
+The suitably renamed ''scene.html'' can then be uploaded to a website, as done here.