CRL 618

4 feb '12.

Three runs from the table below have completed. Bin is in the process of producing movies and plots:

​https://clover.pas.rochester.edu/trac/astrobear/blog/blin02052012

Four of the new runs (see table below) are waiting in bluegene's queue. We should have them -along with plots, diagrams, etc.- in a week time, before Adam's visit to Bruce. Additionally, Bin has 2 of these sims in the standard queue at bluehive. Note that the grid is slightly larger than before so we can see the edges of the object. Also, I have reduced 1 amr level because it adds time to the runs and post processing and doesn't give us significant info about the object's dynamics. Bruce and I have discussed about the separation of the ambient rings. The setup in runs 4, 5 and 8 (see table) seems to be consistent with the obs.

Telecon with Bruce, Adam and Jason on 26 jan '12.

  • We don't need high resolution (256x128x128+2amr) to compare with the observations, which is the main objective of the paper. We will drop 1, 2 may be, refinement levels which should decrease the simulation production time.
  • We don't need a constant ambient density model
  • We'll do velocity and density line-outs along the jet-axis at y={-2/4,-¼,0,¼,2/4} Rclump
  • We'll do PV diagrams with a slit displaced Rclump/2 from the symmetry axis
  • Separation of the ambient rings form observations: ~1e16 cm = 0.9kpc*3e21cm/kpc/206265. We're using a separatin that is twice as long in the models for some experimentation showed the rings to be too close otherwise.

New simulations:

Model nameRunning orderRunning schedulesNo. clumpsJetDens contrast with ambientAmbientResolution
a127 Jan. Completed10200 (as the previous runs)r-2+rings128x80x80+2amrhttp://www.pas.rochester.edu/~blin/feb2012/clump-strat0110.jpeg
b227 Jan. Completed01200r-2+ringsdittohttp://www.pas.rochester.edu/~blin/feb2012/jet-strat0110.jpeg
c328 Jan. Completed1050r-2150x100x100+2amrhttp://www.pas.rochester.edu/~martinhe/2011/5feb12-515.png
d428 Jan Completed1050r-2+ringsdittohttp://www.pas.rochester.edu/~martinhe/2011/4feb1420.png
e530 Jan Running0150r-2ditto
f631 Jan Running0150r-2+ringsditto
g731 Jan Queued2; one at the initial condition and another after the 1st one has left the grid050,50r-2ditto
h81 Feb Queued2; one at the initial condition and another after the 1st one has left the grid050,50r-2+ringsditto

New synthetic diagrams along the correct slit position angle (turns out you cannot turn the slit in shape. The PV diagrams that I posted on Saturday correspond to a slit position angle of 0 degrees, but we want this to be 90deg. Thus I've rotated the objects with respect to the slit and had used more data from the simulations in order to see bow shocks).

Notes:

  • Clump/jet are shown going up because shape's slit is like that (there's no way around this in shape. I could rotate the synthetic images later on (with gimp) for the paper versions).
  • Gray scales in the left column (density maps) are logarithmic and have the same limits for all images.
  • Gray scales in the emission columns (2 and 4) are logarithmic and do not have the same limits.
  • Gray scales in the pv columns (3 and 5) are linear and do not have the same limits.

Clump, stratified ambient density:

Log(dens) [cm-3] Emission (dens2) 90oPV 90oEmission (dens2) 60oPV 60o
http://www.pas.rochester.edu/~martinhe/c-s-dens0001.jpeghttp://www.pas.rochester.edu/~martinhe/c-s-emiss-90.jpghttp://www.pas.rochester.edu/~martinhe/c-s-pv-90.jpghttp://www.pas.rochester.edu/~martinhe/c-s-emiss-60.jpghttp://www.pas.rochester.edu/~martinhe/c-s-pv-60.jpg

Clump, constant ambient density:

Log(dens) [cm-3] Emission (dens2) 90oPV 90oEmission (dens2) 60oPV 60o
http://www.pas.rochester.edu/~martinhe/c-n-dens0000.jpeghttp://www.pas.rochester.edu/~martinhe/c-n-emiss-90.jpghttp://www.pas.rochester.edu/~martinhe/c-n-pv-90.jpghttp://www.pas.rochester.edu/~martinhe/c-n-emiss-60.jpghttp://www.pas.rochester.edu/~martinhe/c-n-pv-60.jpg

Jet, stratified ambient density:

Log(dens) [cm-3] Emission (dens2) 90oPV 90oEmission (dens2) 60oPV 60o
http://www.pas.rochester.edu/~martinhe/j-s-dens0000.jpeghttp://www.pas.rochester.edu/~martinhe/j-s-emiss-90.jpghttp://www.pas.rochester.edu/~martinhe/j-s-pv-90.jpghttp://www.pas.rochester.edu/~martinhe/j-s-emiss-60.jpghttp://www.pas.rochester.edu/~martinhe/j-s-pv-60.jpg

Jet, constant ambient density:

Log(dens) [cm-3] Emission (dens2) 90oPV 90oEmission (dens2) 60oPV 60o
http://www.pas.rochester.edu/~martinhe/j-n-dens0000.jpeghttp://www.pas.rochester.edu/~martinhe/j-n-emiss-90.jpghttp://www.pas.rochester.edu/~martinhe/j-n-pv-90.jpghttp://www.pas.rochester.edu/~martinhe/j-n-emiss-60.jpghttp://www.pas.rochester.edu/~martinhe/j-n-pv-60.jpg

Some notes comparing Dennis et al. 2008 with our sims.

-His clump/jet vels are 100m/s. Ours are 400km/s.

-His emission and pv diagrams are for 2.5D so that's why they have more resolution, but his 3d emission maps have a little less resolution than our shape ones.

-His Fig8 is the average vx. Our vx lineouts are not averaged; they are along a single line n the x direction.

See:

​https://clover.pas.rochester.edu/trac/astrobear/blog/blin01162012

and subsequent comments.

Comments

1. martinhe -- 13 years ago

Bruce's comments on the 26 Jan '12 teleconf: ​http://www.pas.rochester.edu/~martinhe/balick26jan12.pdf