12.14.2020 New setup runs
See bliu-2020/12/14 for details..
Lx = 6cm
dx=6/120/2^5=15.625e-4 cm
d_wire=0.025 cm.
Size_ratio=16
12.7.2020 New setup runs
11.20.2020 Meeting summary
bliu-2020/11/18
11.10.2020 Simulations Summary
1. Setup 1
Box: 160x160 mm
radius_wire = 2 mm
distance_between_centers = 9 mm
MagField_direction = 0, 0, -1
rhoWire = 2.86e19 g/cc
tempWire = 1.39e3 K
WindMaterial = 27
rhoWind = 2.86e+17 1/cc
velWind = 6e1 km/s
BWind ! varies see table
TempWind = 12 ev or 1.39e5 K
rhoAmb = 2.86e+17 1/cc
tempAmb = 1.39e5 K
velAmb = 6e1 km/s
dx: 0.3125mm
runTime: 4.92 ms or 1.8 domain crossing time or 32.8 wire distance crossing time..
imomentumProtect=1
| Runs | Results | diffusion parameter
|
| hydro, no cooling | rho; | diff_Alpha2=9d4
|
| hydro, Al cooling | rho; Temp; Mach; | diff_Alpha=0.1
|
| Bz=0.1T, beta=138, no cooling | rho; Mach; mag pressure; | diff_Alpha2=9d4
|
| Bz=1T, beta=1.38, no cooling | rho; Mach; mag pressure; | diff_Alpha2=9d4
|
| Bz=1T, beta=1.38, no cooling | rhoScaled; Mach; mag pressure; | diff_Alpha=0.1
|
| Bz=5T, beta=0.055, no cooling | rhoScaled; Temperature; Mach; mag pressure; | diff_Alpha=0.1
|
| Bz=5T, beta=0.055, Al cooling | rhoScaled; Temperature; Mach; mag pressure; | diff_Alpha=0.1
|
| Bz=10T, beta=1.38e-2, no cooling | rhoScaled; Mach; mag pressure; | diff_Alpha=0.1
|
2. Setup 2
Box: 20x20 mm
radius_wire = 0.25 mm
distance_between_centers = 7.5 mm
MagField_direction = 0, -1, 0
rhoWire = 2.86e19 g/cc
WindMaterial = 27
rhoWind = 2.86E+17 1/cc
velWind = 6e1 km/s
BWind !varies see table
TempWind = 12 ev or 1.39e5 K
dx: 0.03125mm
runTime: 0.924 ms or 3.6 domain crossing time, or 9.6 wire distance crossing time
imomentumProtect=1
| Runs | Results
|
| hydro, no cooling | rho scaled; Mach;
|
| By=0.1T, beta=138, no cooling | 
|
3. Setup 3
Box: 20x20 mm
radius_wire = 0.25 mm
distance_between_centers = 3.0 mm
MagField_direction = 0, 0, -1
rhoWire = 2.86e19 g/cc
WindMaterial = 27
rhoWind = 2.86E+17 1/cc
velWind = 6e6 cm/s
BWind ! varies see table
TempWind = 12 ev or 1.39e5 K
rhoAmb = 2.86e+17 1/cc
tempAmb = 1.39e5 K
velAmb = 6e1 km/s
dx: 0.3125mm
runTime: 0.924 ms or 3.6 domain crossing time, or 9.6 wire distance crossing time
imomentumProtect=1
11.2.2020 fix high-density tail issue, non-cooling
Updated boundary condition in the reflect() surface function in sweep_scheme.f90 and fixed the high-density tail issue behind the wires
! Now the reflective surface is just to the right of ul
! Need to set ur based on desired boundary properties
!!if (ul(ivx) > 0) then
!! ur=ul
!! ur(ivx)=ur(ivx)*vfact
!!else
!! ! ur=ul !! vn = 1 for flow towards the boundary always
!! continue
!! !leave ur alone
!!end if
!if (ul(ivx) > 0) then
ur=ul
ur(ivx)=ur(ivx)*vfact
!else
! ur=ul !! vn = 1 for flow towards the boundary always
! continue
!leave ur alone
!end if
10.26.2020 non-cooling Larger box
set vel_amb = 0 doesn't help the high-density tail issue behind the clump for low mach case
high mach case also has a high-density tail for a short time.
Exp set up: B=10^4 Tesla, beta=1.3e-8 ?
10.19.2020 non-cooling Larger box
nWire=100*nWind
vAmb = vWind =6e1 km/s
problem.f90
10.12.2020 experiment runs
With rhoWire = 8.73 g/cc and keeping the pressure equilibrium, the temperature of wire will be <1K!
nWind = 2.86e17 1/cc
tempWind = 1.39e5 K
rhoWire = 8.73 g/cc
copper atom mass: 62.93*amu
nWire ~ 8.73e22 1/cc
tempWire ~ 4.55e-1 K
TimeScale 1.23E-08
Lbox 2 cm
Lscale 1
vWind 6.00E+05 cm/s
round 3
tFinal 8.13E+02 c.u
10.5.2020 Setup
set up 1
Box: 20x20 mm
radius_wire = 0.2 mm
distance_between_centers = 0.9 mm
MagField_direction = 0, 0, -1
rhoWire = 8.73 g/cc
WindMaterial = 27
rhoWind = 2.86E+17 1/cc
velWind = 6e6 cm/s
BWind = 5e4 Gauss
TempWind = 12 ev or 1.39e5 K
set up 2
Box: 20x20 mm (?)
radius_wire = 0.25 mm
distance_between_centers = 7.5 mm
MagField_direction = 0, -1, 0
rhoWire = 8.73 g/cc
WindMaterial = 27
rhoWind = 2.86E+17 1/cc
velWind = 6e6 cm/s
BWind = 5e4 Gauss
TempWind = 12 ev or 1.39e5 K
set up 3
Box: 20x20 mm (?)
radius_wire = 0.25 mm
distance_between_centers = 3.0 mm
MagField_direction = 0, 0, -1
rhoWire = 8.73 g/cc
WindMaterial = 27
rhoWind = 2.86E+17 1/cc
velWind = 6e6 cm/s
BWind = 5e4 Gauss
TempWind = 12 ev or 1.39e5 K
HeFrac = .084d0 ! Ratio of Helium Nuclei to Hydrogen Nuclei
xH = .01d0 ! Ionization fraction of Hydrogenn
xHe = 0d0 ! Ionization fraction of Helium (singly ionized)
xHeII = 0d0 ! Ionization fraction of Helium (doubly ionized)
ZFrac = .00126d0 ! Ratio of Metal Nuclei to Hydrogen Nuclei
xZ = .001d0 ! Ionization fraction of metals
Zmu = 16d0 ! Mean atomic mass of metals - (including electrons from ionized metals)
Xmu = 0d0 ! Mean atomic mass of everything - (including all free electrons and nuclei)
9.2.2020
9.1.2020 from JC
8.24.2020
OutflowClumps
rho_outflow 0.5
rho_ambient 0.1
rho_clump 1
mach =10
beta=10
Documents
- MHD votex
