1
0.1
allface
1
300
0
300
0.1
3
300
0
0
0
mass = (0.1 * 44 + 0.2 * 18 + 0.7 * 28);
V = 8.314 * temperature / 1e5;
density = mass * 1e-3 / V;
1.17862
0.01
1.83e-05
Y1 = 0.1 * 44 / (0.1 * 44 + 0.2 * 18 + 0.7 * 28);
Y2 = 0.2 * 18 / (0.1 * 44 + 0.2 * 18 + 0.7 * 28);
Y3 = 0.7 * 28 / (0.1 * 44 + 0.2 * 18 + 0.7 * 28);
CP1 = 1000 * (7.70 + 0.005390 * temperature - 0.830 * 1e-6 * temperature * temperature);
CP2 = 1000 * (6.89 + 0.003284 * temperature - 0.343 * 1e-6 * temperature * temperature);
CP1 = 1000 * (6.30 + 0.001819 * temperature - 0.345 * 1e-6 * temperature * temperature);
1017.24
0.02495
101325
300
0
0
0
0
0
0
all[]
0.2
3
1
3
2
2
1
te=800.0;
r=sqrt(x*x+y*y);
if(r<1){
if(z<0.375){
tc=400+(1800-400)/0.375*z;
temperature=(tc-te)*(1-3*r*r+2*r*r*r)+te;
}
else{
tc=1800+(800-1800)/(4-0.375)*(z-0.375);
temperature=(tc-te)*(1-3*r*r+2*r*r*r)+te;
}
}
else{
temperature = 800.0;
}
1e+12
1
SGDH
1
3
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
2
1
1