Changes between Version 22 and Version 23 of u/smurugan/pn2013

Timestamp:

07/23/13 10:20:21 (12 years ago)

Author:

smurugan

Comment:

—

u/smurugan/pn2013

@@ -18 +18 @@
 * '''Diverging Wind: ''' This is the same setup as the jet, but instead of injecting gas directly along the y axis, it will sweep an arc of 10^o^.
 
-We also introduce toroidal magnetic fields in the MHD runs, which are initialized off the grid, according to the specifications in [http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1989ApJ...344...89L Lind et. al, 1989 ApJ]. We choose Rm = .6Rjet, or 300 AU, based on previous papers. This parameter measures how electric current flows along the jet. To determine initial field strength, we use the magnetic [[latex($\beta$)]] (not to be confused with [[latex($\beta_{t}$)]] which is associated with the shape of the torus), which is the ratio of the thermal pressure to the magnetic pressure. High values of [[latex($\beta$)]] will mean weak magnetic fields. For our initial runs, we chose [[latex($\beta = 1$)]].  
+We also introduce toroidal magnetic fields in the MHD runs, which are initialized off the grid, according to the specifications in [http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1989ApJ...344...89L Lind et. al, 1989 ApJ]. We choose Rm = .6Rjet, or 300 AU, based on previous papers. This parameter measures how electric current flows along the jet. To determine initial field strength, we use the magnetic [[latex($\beta$)]] (not to be confused with [[latex($\beta_{t}$)]] which is associated with the shape of the torus), which is the ratio of the thermal pressure to the magnetic pressure. High values of [[latex($\beta$)]] will mean weak magnetic fields. For our initial runs, we chose [[latex($\beta = 1$)]]. 
 
 
@@ -36 +36 @@
 ||= '''Jet''' =|| [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/jet.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/jet.png, 400px)]] [[BR]] ''Status: Completed.'' ||  [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/stratified_jet.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/stratified_jet.png, 400px)]] [[BR]] ''Status: Completed.''  ||  [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/toroidal_jet.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/toroidal_jet.png, 400px)]] [[BR]] ''Status: Completed.''  ||
 ||= '''Diverging  Winds''' =||  [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/diverging_wind.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/diverging_wind.png, 400px)]] [[BR]] ''Status: Completed.''   ||  [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/stratified_diverging_wind.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/stratified_diverging_wind.png, 400px)]] [[BR]] ''Status: Completed.''   ||  [http://www.pas.rochester.edu/~smurugan/pn2013/mhd/toroidal_diverging_wind.gif movie][[Image(http://www.pas.rochester.edu/~smurugan/pn2013/mhd/toroidal_diverging_wind.png, 400px)]] [[BR]] ''Status: Completed.''  ||
+
+
+
+
+== Tests ==
+
+=== Half Grid v.s Full Grid ===
+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]]
+[[Image(http://www.pas.rochester.edu/~smurugan/pn2013/martin/71213/fullgrid.png, 400px)]] 
+
+
+=== MHD Bubbles ===