MHD Clump Research

---

Relevant research papers and test runs.

1. Papers
• Delamarter 2001, semi-analytic test of cooling method.
• Poludnenko et al. 2002, 2D adiabatic hydrodynamic clump shocks.
• Anninos, Fragile & Murray 2004, new equilibrium cooling function.
• Fragile et al. 2004, several alternative models for radiative cooling.
• Poludnenko, Frank & Mitran 2004, eta-carinae paper (alternate parameters for bullet test).
• Fragile et al. 2005, 2D MHD clump shocks with radiative cooling.
• Patnaude & Fesen 2005, hydro shock model of the Cygnus loop.
• Nakamura et al. 2006, resolution tests.
• Yirak et al. 2007, YSO jets and clumps.
• Shin et al. 2008, 3D MHD adiabatic clump shocks.
• Inoue et al. 2009, magnetic field amplification.


2. Links
• Delamarter Cooling Test Plots
• Alexei's bullet simulations



3. Paper draft (11/18/2009, rev 58).
• To check out the paper materials using Subversion, type

  svn co https://data-store.pas.rochester.edu/svn/orda/papers/

  at a command prompt.



4. Halifax Images.
1. Bx Dense Rho With Fields.
2. By Dense Rho With Fields.
3. Comparison Dense Rho With Fields.
4. Bx Sparse Rho With Fields.
5. By Sparse Rho With Fields.
6. Comparison Sparse Rho With Fields.
7. By Sparse Rho With Fields, Hotboxes.
8. Bx Sparse Vortensity With Fields (reflected).
9. By Sparse Vortensity With Fields (reflected).


10. Bx Sparse Vortensity With Fields (.gif).
11. By Sparse Vortensity With Fields (.gif).


5. Paper-Ready Images. Bold items are awaiting approval for check-in.
1. 2D vortensity images (dense).
2. 2D vortensity images (sparse).
3. 2D vortensity images w/ fieldlines (dense).
4. 2D vortensity images w/ fieldlines (sparse).


5. Emission in shocked region vs. time (dense).
6. Emission in shocked region vs. time (sparse).


7. Vortensity histogram (volume fraction).
8. Vortensity histogram (Cumulative volume fraction).
9. Vortensity histogram (mass fraction).
10. Vortensity histogram (Cumulative mass fraction).


11. 1D y-averaged wall emission plots.
12. 1D y-averaged clump emission plots.
13. 1D y-averaged dense emission plots.
14. 1D y-averaged sparse emission plots.


15. Table of total emissions for hydro clump and wall cases.
16. Hydro wall shock emission plots.


17. Fig 1 (annotated hydro clump)


18. Fig 2a (sparse hydro case).
19. Fig 2b (sparse Bx case).
20. Fig 2c (sparse By case).
21. Fig 2d (dense hydro case).
22. Fig 2e (dense Bx case).
23. Fig 2f (dense By case).


24. Fig 4a (Front Shock Position).
25. Fig 4b (Center of Mass).
26. Fig 4c (Aft Shock Position).


27. Fig 5a (dense emission).
28. Fig 5b (sparse emission).


29. Fig 6 (Average beta vs. time).
30. Fig 6 (Average inverse beta vs. time).


31. Fig 7a (beta histogram, unweighted and non-cumulative).
32. Fig 7a (beta histogram, unweighted and cumulative).
33. Fig 7b (beta histogram, density-weighted and non-cumulative).
34. Fig 7b (beta histogram, density-weighted and cumulative).


35. Fig 8a (B/B0 histogram, unweighted and non-cumulative).
36. Fig 8a (B/B0 histogram, unweighted and cumulative).
37. Fig 8b (B/B0 histogram, density-weighted and non-cumulative).
38. Fig 8b (B/B0 histogram, density-weighted and cumulative).


39. Fig 9a (tracer/rho images, dense runs).
40. Fig 9b (tracer/rho images, sparse runs).


41. Fig 10 (Mixing Parameter).
42. Fig 11 (Average Vortensity).


43. Fig 13a (normalized kinetic energy).
44. Fig 13b (normalized thermal energy).
45. Fig 13c (normalized magnetic energy).
46. Fig 13d (normalized radiative energy).
47. Fig 13e (Initial-energy normalized thermal energy).
48. Fig 13f (Initial-energy normalized magnetic energy).


6. Beta = 2/3, Bx direction.
1. Density movie.


7. Vortensity plots (overall vortensity).
1. Vortensity Composite Plot.
2. Vortensity Variance Composite Plot.
3. Vortensity (All Runs).
4. Vortensity Variance (All Runs).


8. Vortensity plots (regions with > 10% tracer).
1. Mixed-Region Vortensity Composite Plot.
2. Mixed-Region Vortensity (All Runs).
3. Mixed-Region Fractional Vortensity (All Runs).


9. Plots of Single Clump Mixing Parameter (number of cells w/ > 10% clump material).
1. Composite Plot.
2. All Runs.


10. Single-Clump Animations.
1. Bx comparison.
2. By comparison.
3. dense comparison.
4. sparse comparison.


11. Plots of Mixing Parameter (number of cells w/ > 10% clump material).
1. Composite Plot.
2. All Runs.


12. Plots of x-position of Center of Mass
1. Composite Plot.
2. All Runs.
3. Quadratic curve fits (all runs).


13. Plots of Front and Rear Shock Position
1. Dense Run.
2. Sparse Run.
3. Bx Run.
4. By Run.


14. Plots of Clump Material Ratios within Shock Region
1. Dense Run.
2. Sparse Run.
3. Bx Run.
4. By Run.
5. all runs.


15. MHD Bx-By Comparisons (Animated)
1. dense rho plot.
2. dense Bmag plot.
3. dense beta plot.
4. dense schlieren plot.
5. dense vortensity plot.
6. dense emission plot.


7. sparse rho plot.
8. sparse Bmag plot.
9. sparse beta plot.
10. sparse schlieren plot.
11. sparse vortensity plot.
12. sparse emission plot.


16. MHD Dense-Sparse Comparisons (Animated)
1. Bx rho plot.
2. Bx Bmag plot.
3. Bx vortensity plot.
4. Bx schlieren plot.
5. Bx beta plot (with field lines).
6. Bx emission plot .


7. By rho plot.
8. By Bmag plot.
9. By linear schlieren plot.
10. By vortensity plot.
11. By beta plot (with field lines).
12. By emission plot .




17. MHD Dense-Sparse Comparisons (Plots)
1. Bx Center of mass (x-coordinate) comparison plot.
2. Bx Mixing parameter comparison plot (area-based).
3. Bx Single-Clump Mixing param comparison plot (cell-based).
4. Bx Single clump tracer.
5. Bx Energy comparison plots.
6. Bx Comparison plot of total vortensity.
7. Bx Horizontal lineout of vorticity (animated).
8. Bx Vertical lineout of vorticity (animated).
9. Bx Energy comparison plots.


10. By Center of mass (x-coordinate) comparison plot.
11. By Mixing parameter comparison plot (area-based).
12. Single-Clump Mixing parameter comparison plot (cell-based).
13. Single clump tracer.
14. Bx Comparison plot of total vortensity.
15. Bx Horizontal lineout of vorticity (animated).
16. Bx Vertical lineout of vorticity (animated).
17. By Energy comparison plots.


18. MHD Sparse Clump Runs (Incomplete)
1. By density plot (large)
2. By Bmag plot (large).
3. By density plot (dense and sparse runs)
4. By B-field plot (dense and sparse runs)


19. Pairwise Swap Demonstrations
1. Advected Loop Final Frame (serial).
2. Advected loop Final Frame (pairwise).
3. Magnetic Pressure Comparison Plot (Field Loop Advection).


1. Comparisons of dense and sparse clump runs.
1. Hydro density.
2. Hydro schlieren.
3. Hydro vortensity.


2. Magnetic Field Property Plots.
1. Comparison plot of maximum field strengths (2R runs).


3. Sparsely-populated MHD clump initial frames.
1. 4 * clump radius
2. 6 * clump radius


4. Plots of Tracer-Derived Variables.
1. Mixing Parameter Plot.
2. Center of Mass Plots.


5. Comparison Plots of Normalized Energy.
1. Normalized Energy Comparison Plots.
2. Normalized Fluxless Energy Comparison Plots.
3. Magnetic Energy Comparison Plots.
4. Analytic (Uncooled) vs. Actual (Cooled) Comparison Plots.
5. Total Energy Breakdown Comparison Plots.


6. MHD Wind Clump Plots: Density, B-field, and Beta.
1. Density comparison plots.
2. Magnetic field comparison plots.
3. beta^-1 comparison plots.


4. Schlieren comparison plots.
5. Vortensity comparison plots.


6. Vorticity horizontal sum plot (all).
7. Vorticity vertical sum plot (all).


8. Vorticity total sum plot.
9. sqrt(vortensity) total sum plot.


7. 2D MHD Wind Clump Separation Plots, 200 frames, Scale = 400 AU, beta = 10, rho_amb = 100 cm^-3
1. Hydrodynamic density plot.


2. Bx density plot.
3. Bx magnetic field plot.
4. Bx beta^-1 plot with field lines.
5. Energy Plot (Bx)
6. Fractional Energy Plot (Bx)


7. By density plot.
8. By magnetic field plot.
9. By beta^-1 plot with field lines.
10. Energy Plot (By)
11. Fractional Energy Plot (By)


12. Magnetic Energy Plot (Bx, By)


8. 2D MHD Wind Clump Separation Plots, scale = 400 AU, random perturbation placement, beta = 10, ambient density = 100 cm^-3
1. Radius = 200 AU, Minimum Separation = 2R, Bx (density plot).
2. Radius = 200 AU, Minimum Separation = 2R, Bx (fixed density plot).
3. Radius = 200 AU, Minimum Separation = 2R, Bx (B-field plot).
4. Radius = 200 AU, Minimum Separation = 2R, Bx (beta^-1 plot).


5. Energy Plots of Bx simulation.
6. 
   
7. Radius = 200 AU, Minimum Separation = 2R, By (density plot).
8. Radius = 200 AU, Minimum Separation = 2R, By (flat scale B-field plot).
9. Radius = 200 AU, Minimum Separation = 2R, By (deep scale B-field plot).


9. 2D MHD Wind Clump Separation Plots, scale = 400 AU, beta = 10, ambient density = 100 cm^-3
1. Radius = 200 AU, Separation = 400 AU, Bx
2. Radius = 200 AU, Separation = 400 AU, Bx (fixed)


3. sparse run (28 clumps)
4. dense run (54 clumps)


10. 2D MHD Wind Clump Separation Plots, scale = 100 AU, beta = 10, ambient density = 100 cm^-3
1. Radius = 200 AU, Separation = 400 AU, Bx
2. Radius = 200 AU, Separation = 400 AU, By


11. 2D MHD Wind Clump Plots, scale = 100 AU, equal wind pressure, beta = 10, ambient density = 100 cm^-3
1. 1 clump, Mach 10, Bx, density plot (with field lines, boxes).
2. 1 clump, Mach 10, By, density plot (with field lines, boxes).


3. 5 clumps, Mach 10, Bx, density plot (with field lines, boxes).
4. 5 clumps, Mach 10, By, density plot (with field lines, boxes).


5. 14 clumps, Mach 10, Bx, density plot (with field lines, boxes).
6. 14 clumps, Mach 10, By, density plot (with field lines, boxes).


7. 14 clumps, Mach 10, beta = 4, By, density plot (no field lines).


12. 2D MHD Wind Clump Plots, scale = 400 AU, equal wind pressure, ambient density = 100 cm^-3
1. Mach 10,beta = 10, 64 cells/radius, density plot (with field lines).
2. Mach 10,beta = 10, 64 cells/radius, density plot (no field lines).
3. Mach 10,beta = 10, 64 cells/radius, B-field plot (magnitude).


4. Mach 10,beta = 10, 64 cells/radius, CFL = .25, density plot.
5. Mach 10,beta = 10, 64 cells/radius, cooling off, density plot.


13. 2D Wind Clump Plots, scale = 400 AU, equal wind pressure, ambient density = 100 cm^-3
1. hydro, Mach 10, 64 cells/radius, density plot.
2. hydro, Mach 10, 64 cells/radius, density plot (small).


3. hydro, Mach 25, 64 cells/radius, density plot.


14. 2D Wind Clump Plots, scale = 400 AU, ambient density = 100 cm^-3
• hydro, Mach 10, 64 cells/clump, density plot.
• hydro, Mach 10, 128 cells/clump, density plot.
• hydro, Mach 25, 64 cells/clump, density plot.


15. 2D Shock Clump Plots, scale = 400 AU, Mach = 50, ambient density = 100 cm^-3
• hydro, 128 cells/clump, density plot.


16. Plots of 2D Bullets Using Cylindrical Symmetry and Radius = 200 AU
• log density plot, rho_amb = 10^2, mach = 20.
• log density plot, rho_amb = 10^2, mach = 20, no viscosity.
• log density plot, rho_amb = 10^2, minTemp = 9000 K, mach = 20.
• log density plot, rho_amb = 10^3, mach = 20.


17. Plots of Clumps in Pressure Balance (Mach 1) using Radius 200 AU
• Bullet in pressure balance (rho = 10^2).
• Bullet in pressure balance (rho = 10^2, tanh smoothing compensation).
• Bullet in pressure balance (rho = 10^2, mach = 1, corrected.


18. Plots of 2D Bullets Using Mach 20, Radius 50 AU, ambient density 10^2 cm^-3
• Density plot (greyscale).
• Density plot (Schlieren).


19. Plots of 2D Clump Shocks Using Mach 50, ambient density 10 cm^-3
• Density plot (greyscale).


• Temperature plot (VisIt).
• Temperature plot (greyscale).


20. Plots of 2D Clump Shocks Using Mach 50, ambient density 1 cm^-3
• Density plot (VisIt).
• Density plot (greyscale).


• Temperature plot (VisIt).
• Temperature plot (greyscale).


21. Temperature Plots of 2D Clump Shocks Using Mach 50, ambient density 0.1 cm^-3
• 1D plot, ambient density 0.1, adiabatic.
• 2D plot, ambient density 0.1, adiabatic.
• 1D plot, ambient density 0.1, radiative cooling.
• 2D plot, ambient density 0.1, radiative cooling.


• 1D plot, rho_amb = 0.1 (tall), radiative cooling.
• 2D plot, rho_amb = 0.1 (tall), radiative cooling.
• 2D greyscale plot, rho_amb = 0.1 (tall), radiative cooling.
• 2D density plot, rho_amb = 0.1 (tall), radiative cooling.
• 2D greyscale density plot, rho_amb = 0.1 (tall), radiative cooling.


• 1D plot, B-field in x, MHD + cooling.
• 2D field-line plot, B-field in x, MHD + cooling.
• 2D plot, B-field in x, MHD + cooling.
• 1D plot, B-field in y, MHD + cooling.
• 2D field-line plot, B-field in y, MHD + cooling.
• 2D plot, B-field in y, MHD + cooling.


22. Temperature Plots of 1D Clump Shocks Using Mach 25, minimum temp 9000 K
• ambient density 0.1, ambient temp = 100 K.
• ambient density 1, ambient temp = 100 K.
• ambient density 10, ambient temp = 100 K.
• ambient density 100, ambient temp = 100 K.


• cell plot, ambient density 10, ambient temp = 100 K.


23. Temperature Plots of 1D Clump Shocks Using Non-Fragile Parameters
• adiabatic, ambient density 0.1, Mach 50, ambient temp = 100 K.
• cooling, ambient density 0.1, Mach 50, ambient temp = 100 K.
• cooling (hi-res), ambient density 0.1, Mach 50, ambient temp = 100 K.
• cooling, ambient density 0.1, Mach 10, ambient temp = 100 K.


24. Temperature Plots of 1D Clump Shocks Using Fragile Parameters
• adiabatic, ambient density 0,1, chi 100, rclump = 1 pc.
• cooling, ambient density 0.1, chi 100, rclump = 1 pc.
• adiabatic, ambient density 10, chi 1000, rclump = 1 pc.
• cooling, ambient density 10, chi 1000, rclump = 1 pc.


25. Temperature Plots of Steady Shocks
• Temperature Plot, density = 1.
• Temperature Plot, density = 5.
• Temperature Plot, density = 10.
• Temperature Plot, density = 50.
• Temperature Plot, density = 100.


• Cooling Length Plots.


26. 14 Randomly-Distributed Clumps
• Density Plot, Fixed-Grid, Bx, Wide-Ranging Densities
• Schlieren Plot, Fixed-Grid, Bx, Wide-Ranging Densities
• Density Plot, 2 Extra Levels AMR, Bx, Constrained Density Range


27. 3 Clumps, Beta = 1
    • Density Plot, Fixed-Grid, Bx
    • Density Plot, 2 Extra Levels AMR, By (61 frames on Calhoun)
    • Density Plot, 2 Extra Levels AMR, Bx (Full 100 frames on Nova)
    
    
28. 3 Clumps, Beta = 4
    • Density Plot, 2 Extra Levels AMR, By