High Res MHD Colliding Flows

Here is the page that consists of our diagnostics and figures for the paper: ​CollidingFlowsFigures

Initial conditions:

Beta	10, everywhere
Number density	1cm-3, everywhere
Flow diameter	40 pc
Mach of flows	1.5
Velocity	136 (cu) inside flows, 0 outside
B field	Field runs parallel to flows, everywhere in box
B pressure (gauss)	(B2/2)*Bscale = 20 \mu ~esu = 20 \mu G cm^2
B-mag	Bx = 31.40 (cu), B constant everywhere in box
Sound Speed	90.66 (cu), constant everywhere in box
Pressure	4931 (cu), ~constant everywhere in box
Temp	4931 (cu), n=1

Questions: Field strength is B*Bscale?

Scales:

velscale	8091.24708
timescale	!381360279481384.188
bscale	0.!416530471377859306E-07
nscale	1
lscale	0.!308568024999999590E+19 = 1pc
scalegrav	0.!204668643805330057E-01
pscale	0.!138065029999999951E-15
tempscale	1

velscale = lScale/TimeScale

TimeScale=lscale/sqrt(Boltzmann*TempScale/Xmu/amu)

Bscale=SQRT(4d0 * Pi * pScale)

rScale=Xmu*amu*nScale

ScaleGrav=G*rscale*TimeScale**2

pscale = nscale(cm^{-3})*K_B(cgs)*tempscale(K)

5 Production Runs:

Beta	Shear Angle	Purpose	Final time	Final Frame
10	0	Compare with Jonathan's work	27 Myr	273
10	15	Compare with Christina's work	" "	"
10	30	" "	" "	"
10	60	" "	complete out to 32.5 Myr	328
1	60	Compare a strong field	20 Myr	200

For orientation, here is an image that shows the flows and tilted interface,

TableView of Column Density Maps

	Beta 1, Shear 60	Beta 10, Shear 60	Shear 30	Shear 15	No Shear
down-the-barrel (sight along the flow)
movies	b1s60cdm1.gif​	b10s60cdm1.gif​	b10s30cdm1.gif​	b10s15cdm1.gif​	b10noshearcdm1.gif​
down-x (normal angle to flow, colliding from top and bottom)
movies	b1s60cdm2.gif​	b10s60cdm2.gif​	b10s30cdm2.gif​	b10s15cdm2.gif​	b10noshearcdm2.gif​
down-y (normal angle to flow, colliding from left and rt.)
movies	b1s60cdm3.gif​	b10s60cdm3.gif​	b10s30cdm3.gif​	b10s15cdm3.gif​	b10noshearcdm3.gif​

CDM Movies by run

B10, Shear 0
movies	B10_0_vistacdm.gif​
B10, Shear 15
movies	b10s15_galaxy_cdm.gif​
B10, Shear 30
movies	b10s30_galaxy_cdm.gif​
B10, Shear 60
movies	b10s60_galaxy_cdm.gif​
B1, Shear 60
movies	b1s60_galaxy_cdm.gif​

Sink Formation

Sink Formation

Comparing sink formation for different runs, and for the high vs low resolution runs. As expected, low res runs produce more sinks. Both resolution runs don't produce sinks in the strong field/shear case. The no shear case wasn't run in low resolution also, possibly misleading in plot.

Interesting how the high and low resolution runs take on opposite trends here,

Same high res data in tabular form:

	Beta 10, Shear 60	Shear 30	Shear 15	No Shear
No. sinks by 27 Myr:	0	1	1	4 (1 more sink between 20&27 Myr)
First sink forms at frame:	—	138	188	111
Time (Myr):	—	13.62	18	11

For two levels of refinement:

	Beta 1, Shear 60	Beta 10, Shear 60	Shear 30	Shear 15	No Shear
No. sinks by 20 Myr:	0	0	0	4	4
First sink forms at frame:	—	—	—	124	122
Time (Myr):	—	—	—	12.238	12.09

Here​ is a movie showing the S30 case down the barrel and down y (left 2 columns) and S15 down y in right column, since I found it interesting that the low res S30 case didn't make a sink.

Resolution Studies

Resolution Studies:

Was thinking of doing a convergence study for the Shear 60 case so I could see if there was a convergence effect for the different field strengths?

Something like 4 little resolution tests out to frame 50:

Beta	Shear Angle	maxlevel	Status
10	60	3	complete
10	60	4	gives error
10	60	6	complete
1	60	3	complete
1	60	4	gives error
1	60	6	complete

Given the maxlevel = 4 is giving an error, we ran our lower maxlevel out to 3. Also ran 2 convergence runs on the Beta = 10, Shear 15 case to compensate:

Beta	Shear Angle	maxlevel	Status
10	15	4	complete
10	15	6	complete