wiki:u/bliu/MHDJets

Context Navigation

Version 43 (modified by Baowei Liu, 5 years ago) ( diff )
—

Runs following cooling curve 2.24.2020

1 ev ~ 11604.5 K

points picked on the curve:
xp=[0.3 0.5 0.7 0.9 1.2]
yp=xp+log10(11604.5)

T_ps = 10^yp = 1.0e+05*[0.2315    0.3670    0.5816    0.9218    1.8392]

T_ps=0.75*M_H/K_B*Vjet^2 

M_H = 1.6733 * 1e-24;
KB = 1.380658 * 1e-16;

Vjet = sqrt(4*KB*T_ps/M_H/3);

Vjet =   1.0e+06 *[1.5960    2.0093    2.5295    3.1845    4.4982] cm/s
     = [15.96 20.093 25.295 31.845 44.982] km/s
mach=[1.1459    1.4426    1.8161    2.2863    3.2295]

Runs to do

Common parameters:

rho_amb	1e18 1/cm3
T_amb	8.4e4 K
rho_jet	6e18 1/cm3
T_jet	1.4e4 K

Runs:

	T_ps (K)	Vjet (km/s)	Mach
Run1	23150	15.96	1.1459
Run2	36700	20.093	1.4426
Run3	58160	25.295	1.8161
Run4	92180	31.845	2.2863
Run5	183920	44.982	3.2295

Fiducial Runs 2.6.2020

runs	njet	Tjet	vjet	njet*vjet	njetvjetvjet	Pjet	namb	Tamb	Pamb
1	n1=6e18 1/cm3	T1=1.4e4 K	V1=70 km/s	n1*V1	n1*V1²	n1*T1	n1/6	6*T1	n1*T1
2	n1/36	36*T1	36*V1	n1*V1	36n1V1*V1	n1*T1	n1/6	6*T1	n1*T1
3	n1/36	36*T1	6*V1	n1*V1/6	n1V1V1	n1*T1	n1/6	6*T1	n1*T1

Mach 5 Heavy Jets

rho_amb	1e18 1/cm3
T_amb	8.4e4 K
rho_jet	6e18 1/cm3
T_jet	1.4e4 K
vel_jet	70 km/s

no cooling	heavy jets no cooling
Al cooling	heavy jets Al cooling short run

Mach 30 light Jets

rho_amb	1e18 1/cm3
T_amb	8.4e4 K
rho_jet	0.167e18 1/cm3
T_jet	50.4e4 K
vel_jet	2520 km/s

no cooling	mach 30 light jets no cooling
Al cooling	mach 30 light jets Al cooling
	compare

Mach 5 light Jets

rho_amb	1e18 1/cm3
T_amb	8.4e4 K
rho_jet	0.167e18 1/cm3
T_jet	50.4e4 K
vel_jet	420 km/s

no cooling	mach 5 light jets no cooling
Al cooling	mach 5 light jets Al cooling
	compare

Al cooling function from Eddie

cooling.f90; AlCooling.f90; Alcooling.tab; new format Al Cooling table;

cooling function plot

http://www.pas.rochester.edu/~bliu/MHDCollidingFlows/FromEddie/AlCooling.png

Lab runs 2.5.2020

parameters

rho_amb	1e18 1/cm3
rho_jet	6e18 1/cm3
vel_jet 1	60 km/s
vel_jet 2	80 km/s

problem.f90; global.data; physics.data; problem.data; outflow objects

No Cooling, B=0, Tjet=100k, open_angle=15	no cooling 100k
No Cooling, B=0, Tjet=100k, open_angle=0	no cooling 100k thin
No Cooling, B=0, Tjet=100k, open_angle=0	no cooling 100k wide
No Cooling, B=0, Tjet=1000k	no cooling 1000k

Lab runs before 2.5.2020

parameters

rho_amb		1e18 1/cm3
rho_jet		6e18 1/cm3
vel_jet 1		60 km/s
vel_jet 2		80 km/s

problem.f90; global.data; physics.data; problem.data;

No Cooling, B=0, Tjet=100k	no cooling 100k	wrong time scale, should be 100 larger
No Cooling, B=0, Tjet=10000k	no cooling 10000K	wrong time scale, should be 100 larger
Analytic Cooling, B=0, Tjet=10000k	ana cooling 10000k	wrong time scale, should be 100 larger
Z Cooling, B=0, Tjet=10000k	z cooling 10000k	wrong time scale, should be 100 larger

Test Runs before 1.29.2020

3D MHD Colliding Jets

rho_amb	0.01	1/cm3
rho_jet	70000	1/cm3
vel_jet 1	80	km/s
vel_jet 2	60	km/s

No Cooling, B=0, Tjet=10k		rho
No Cooling, B=0, Tjet=1000k		rho
DM Cooling, B=0, Tjet=10k		rho
DM Cooling, B=0, Tjet=1000k		rho
No Cooling, toroidal B=1e-5, Tjet=10k	;;
DM Cooling, toroidal B=1e-5, Tjet=1000k	;	rho;
No Cooling, toroidal B=1e-6, Tjet=10k	;;
DM Cooling, toroidal B=1e-5, Tjet=10k	;;	rho;Bx;By