12000
0
0.03
walls
0
0
0
density = rho0;
5907
0.01
molecular_viscosity = mu0;
0.002037
370
148
101325
50
0
0
0
0
0
0
all[]
loop
pump
off
0.1
0
1
1.0
off
0
## mhd momentum
I=.1; # amps
B=1; # teslay (Ns/C/m
width=.004; # in mm
height=.0016; # cannehl height
A=width*height; # area of the current cross section m^2
f=I/A*B; # N/m^3, kg/m^2/s^2 volumetric force
Su = 0;
Sv = f;
Sw = 0;
dSudu = 0;
dSudv = 0;
dSudw = 0;
dSvdu = 0;
dSvdv = 0;
dSvdw = 0;
dSwdu = 0;
dSwdv = 0;
dSwdw = 0;
1
1e+12
-273.15
1
10
1
0
velocity[0] = .001;
velocity[1] = .002;
velocity[2] = .003;
2
1
1