1
10
10
1
1000
0.1
0.01
0.1
inlet_1
inlet_2
outlet
wall
80000
19
0.14
0.585
79
0.1
0
0
0
0.764
1
1
0
0
0
0
0
0.7
a = 1.4078*10^(-5);
b = -5.5855*10^(-3);
c = -2.8886*10^(-3);
d = 1000.4;
density =
temperature*(temperature*(a*temperature+b)+c)+d;
996.0697
p = 0.3;
a = 1.4078*10^(-5);
b = -5.5855*10^(-3);
c = -2.8886*10^(-3);
d = 1000.4;
d_f =
temperature*(temperature*(a*temperature+b)+c)+d;
d_p = 5680;
density = p*d_f + (1-p)*d_p;
p = 0.7;
a = 1.4078*10^(-5);
b = -5.5855*10^(-3);
c = -2.8886*10^(-3);
d = 1000.4;
d_f =
temperature*(temperature*(a*temperature+b)+c)+d;
d_p = 4199;
density = p*d_f + (1-p)*d_p;
p = 0.9;
a = 1.4078*10^(-5);
b = -5.5855*10^(-3);
c = -2.8886*10^(-3);
d = 1000.4;
d_f =
temperature*(temperature*(a*temperature+b)+c)+d;
d_p = 2719;
density = p*d_f + (1-p)*d_p;
a = 1.4078*10^(-5);
b = -5.5855*10^(-3);
c = -2.8886*10^(-3);
d = 1000.4;
density =
temperature*(temperature*(a*temperature+b)+c)+d;
0.01
a = -1.9296*10^(-9);
b = 4.7256*10^(-7);
c = -4.2088*10^(-5);
d = 1.6947*10^(-3);
molecular_viscosity = temperature*(temperature*(a*temperature+b)+c)+d;
0.000877
p = 0.5;
ita = 1+2.5*(1-p);
a = -1.9296*10^(-9);
b = 4.7256*10^(-7);
c = -4.2088*10^(-5);
d = 1.6947*10^(-3);
molecular_viscosity =ita*(temperature*(temperature*(a*temperature+b)+c)+d);
p = 0.7;
ita = 1+2.5*(1-p);
a = -1.9296*10^(-9);
b = 4.7256*10^(-7);
c = -4.2088*10^(-5);
d = 1.6947*10^(-3);
molecular_viscosity =(temperature*(temperature*(a*temperature+b)+c)+d);
p = 0.9;
ita = 1+2.5*(1-p);
a = -1.9296*10^(-9);
b = 4.7256*10^(-7);
c = -4.2088*10^(-5);
d = 1.6947*10^(-3);
molecular_viscosity =(temperature*(temperature*(a*temperature+b)+c)+d);
a = -1.9296*10^(-9);
b = 4.7256*10^(-7);
c = -4.2088*10^(-5);
d = 1.6947*10^(-3);
molecular_viscosity = temperature*(temperature*(a*temperature+b)+c)+d;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
4181.252
p = 0.3;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
Cp_f = temperature*(temperature*(a*temperature+b)+c)+d;
Cp_p = 455;
specific_heat = p*Cp_f + (1-p)*Cp_p;
p = 0.7;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
Cp_f = temperature*(temperature*(a*temperature+b)+c)+d;
Cp_p = 663;
specific_heat = p*Cp_f + (1-p)*Cp_p;
p = 0.3;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
Cp_f = temperature*(temperature*(a*temperature+b)+c)+d;
Cp_p = 871;
specific_heat = p*Cp_f + (1-p)*Cp_p;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374*10^(-13);
b2 = -2.1701*10^(-9);
c2 = 4.7970*10^(-7);
d2 = 1.3538*10^(-4);
thermal_conductivity = (temperature*(temperature*(a2*temperature+b2)+c2)+d2)*specific_heat;
0.613653
p = 0.3;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374*10^(-13);
b2 = -2.1701*10^(-9);
c2 = 4.7970*10^(-7);
d2 = 1.3538*10^(-4);
lamda_f = (temperature*(temperature*(a2*temperature+b2)+c2)+d2)*specific_heat;
lamda_p = 2;
thermal_conductivity = p*lamda_f + (1-p)*lamda_p;
p = 0.7;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374*10^(-13);
b2 = -2.1701*10^(-9);
c2 = 4.7970*10^(-7);
d2 = 1.3538*10^(-4);
lamda_f = (temperature*(temperature*(a2*temperature+b2)+c2)+d2)*specific_heat;
lamda_p = 100;
thermal_conductivity = p*lamda_f + (1-p)*lamda_p;
p = 0.9;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374*10^(-13);
b2 = -2.1701*10^(-9);
c2 = 4.7970*10^(-7);
d2 = 1.3538*10^(-4);
lamda_f = (temperature*(temperature*(a2*temperature+b2)+c2)+d2)*specific_heat;
lamda_p = 202.4;
thermal_conductivity = p*lamda_f + (1-p)*lamda_p;
a = -1.0224*10^(-4);
b = 2.9201*10^(-2);
c = -1.822;
d = 4209.9;
specific_heat = temperature*(temperature*(a*temperature+b)+c)+d;
a2 = -3.0374*10^(-13);
b2 = -2.1701*10^(-9);
c2 = 4.7970*10^(-7);
d2 = 1.3538*10^(-4);
thermal_conductivity = (temperature*(temperature*(a2*temperature+b2)+c2)+d2)*specific_heat;
101325
20
0
0
-9.81
0
0
0
all[]
porous_1
porous_2
porous_3
fluid
1.0
off
porosity = 0.5;
p1 = 0.5;
p3 = 0.9;
r = sqrt(y^2+z^2);
porosity = (0.068-r)/0.002*(p3-p1)+p1;
porosity = 0.9;
0
1
temperature = 19.;
temperature = 19;
79
19
1
SGDH
0.5
0
1
0
1
3
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
velocity[0] = 0;
velocity[1] = 0.;
velocity[2] = 0.;
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
2
1
1