| 44 | | == Tests == |
| 45 | | |
| 46 | | === Half Grid v.s Full Grid === |
| 47 | | Our simulations take advantage of the symmetry of these problems and only evolve the right half of the grid, which is then simply reflected to provide a visualization for the whole jet. We decided to run one full grid simulation to compare, just to make sure there are no boundary effects. This run is the stratified diverging wind, with the MHD equations turned on and [[latex($\beta = 100000$)]] (so essentially no magnetic fields). If we compare this to the hydro stratified diverging wind run from the previous section, we see that the symmetry of the problem is indeed preserved. [[BR]] |
| 48 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/martin/71213/fullgrid.png, 400px)]] |
| 49 | | |
| 50 | | |
| 51 | | === MHD Bubbles === |
| 52 | | Our initial runs show very distinct bubbles in the MHD runs that don't exist in the hydro case. While the MHD flows should be m = 0 unstable, there may be something else going on as well. The following images are from the same conditions, one purely hydro, one with MHD and [[latex($\beta = .8$)]] and the third MHD with a high [[latex($\beta$)]]. As we cantsee, the bubbles are very nonexistant in the hydro case, distinct in the MHD case and mild in the high [[latex($\beta$)]] case. |
| 53 | | |
| 54 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/martin/71213/hydro.png, 400px)]] |
| 55 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/martin/71213/mhd.png, 440px)]] |
| 56 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/martin/71213/highbeta.png, 500px)]] |
| 57 | | |
| 58 | | We have also tried running the weak magnetic field cases, with [[latex($\beta = 8$)]] and [[latex($\beta = 80$)]]: |
| 59 | | |
| 60 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/tests/beta/stratb8.png, 400px)]] |
| 61 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/tests/beta/stratb80.png, 400px)]] |
| 62 | | |
| 63 | | Martin, 29 July 2013. Found and fixed a bug in our module, which affected the strength of magnetic fields. Running test to see whether this gets rid of the bubble-looking features, which i think are unphysical, that we see in the mhd cases. |
| 64 | | |
| 65 | | The current runs have been updated with a solution to this issue. You can look through the old data [http://www.pas.rochester.edu/~smurugan/pn2013/mhdbubbles here]. |
| 66 | | |
| 67 | | |
| 68 | | |
| 69 | | === Torus === |
| 70 | | |
| 71 | | In our initial runs, the torus case seems to not be very different from the constant ambient case. We are trying different values of [[latex($\alpha_{t}$)]], the density contrast to see where there is a significant impact of the torus. So far, we have tried values of: |
| 72 | | [[BR]] |
| 73 | | [[latex($\alpha = .3$)]] |
| 74 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/tests/tori/a3b8.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/tests/tori/a3b8.gif movie] |
| 75 | | |
| 76 | | [[BR]] |
| 77 | | [[latex($\alpha = .9$)]] |
| 78 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/tests/tori/a9b8.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/tests/tori/a9b8.gif movie] |
| 79 | | |
| 80 | | |
| 81 | | Our current idea is that the setup of the Torus is in such a way that it will not have a great effect on the shaping of the lobe past the origin. Our new runs now reflect a mroe physical situation that combines a 1/r^2^ density dropoff along with the torus setup. You can find the data from our old setup [http://www.pas.rochester.edu/~smurugan/pn2013/torusold here]. |
| 82 | | |
| 83 | | |
| 84 | | |
| 85 | | |
| 86 | | |
| 87 | | === Hi-Res Diverging Wind Model === |
| 88 | | |
| 89 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/tests/hires/dw.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/tests/hires/dw.gif movie] |
| 90 | | |
| 91 | | |
| 92 | | === Comparison of MHD to Hydro Jets in Different Ambient Media=== |
| 93 | | Constant: |
| 94 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/comp/const_comp.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/comp/const_comp.gif movie] |
| 95 | | |
| 96 | | Stratified: |
| 97 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/comp/stratified_comp.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/comp/stratified_comp.gif movie] |
| 98 | | |
| 99 | | Toroidal: |
| 100 | | [[Image(http://www.pas.rochester.edu/~smurugan/pn2013/comp/toroidal_comp.png, 400px)]] [http://www.pas.rochester.edu/~smurugan/pn2013/comp/toroidal_comp.gif movie] |
| 101 | | |
| 102 | | |
