Thoughts

User Tools

Site Tools

Trace:
dustdevils

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
dustdevils [2025/04/01 09:11] – [isosurface animation] admindustdevils [2025/04/01 11:56] (current) – admin
Line 1: Line 1:
-===== Dustdevils redux =====+===== Dust devils redux =====
  
  
-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: 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> <code>
 # choose an integer, see choices for nexp below # choose an integer, see choices for nexp below
Line 48: Line 52:
 Try  Try 
 <code> <code>
-./trace.py R4096g0w4 1 400+./hist.py R4096g0w4 1 400
 </code> </code>
 if all went okay, if all went okay,
 <code> <code>
-./trace.py R4096g0w4 1 400 > R4096g0w4.hist +./hist.py R4096g0w4 1 40 > R4096g0w4.hist 
-</code>+</code>0
  
 And then  And then 
 <code> <code>
-./plottrace.py -s R4096g0w4.hist +./plothist.py -s R4096g0w4.hist 
 </code>  </code> 
  
 That should produce the following image, which shows a dust devil event just prior to t=30. That should produce the following image, which shows a dust devil event just prior to t=30.
- 
  
  
Line 71: Line 74:
 <code> <code>
 mkdir pngs mkdir pngs
-./hcontour.py -o pngs R4096g0w4 000001 000400+./npy2contour.py -o pngs R4096g0w4 000001 000400
 ./maketitle.py -o titleh.png ./maketitle.py -o titleh.png
 ./png2mp4.py pngs/0 -t titleh.png ./png2mp4.py pngs/0 -t titleh.png
Line 86: Line 89:
 <code> <code>
 mkdir vista mkdir vista
-./looknpy.py  R4096g0w4 000270 000320 -o vista+./npy2vista.py  R4096g0w4 000270 000320 -o vista
 ./maketitlevista.py ./maketitlevista.py
 ./png2mp4.py vista/0 -t titlevista.png ./png2mp4.py vista/0 -t titlevista.png
Line 94: Line 97:
  
 ==== Some interactive images at $t=29.8$ ==== ==== Some interactive images at $t=29.8$ ====
 +
 +<code>
 +./npy2vista.py -s R4096g0w4 000298
 +open scene.html
 +</code>
 +
 Use your mouse, have fun ... Use your mouse, have fun ...
  
 {{url>https://12characters.net/rbdd/t298zm40pm2.html}} {{url>https://12characters.net/rbdd/t298zm40pm2.html}}
 $\zeta=-40$ (purple)   , $p=-2$ (green), and $s$ at the surface $\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}} {{url>https://12characters.net/rbdd/t298wpm5c6.html}}
 $w=.5$ (green), $w=-.5$ (grey) , $c=6$ (red) $w=.5$ (green), $w=-.5$ (grey) , $c=6$ (red)
  
 +The suitably renamed ''scene.html'' can then be uploaded to a website, as done here.
dustdevils.1743516702.txt.gz · Last modified: 2025/04/01 09:11 by admin
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki