1
traceur = 0;
traceur = 0;
1e+12
-1e+12
1.83e-05
1
SGDH
x = 27.5+(55-27.5)*s;
y = 0.5;
z = 8.25;
-1
end
200
x = 27.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 28.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 28.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 29.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 29.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 30.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 30.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 31.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 31.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 32.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 32.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 33.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 33.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 34.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 34.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 35.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 35.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 36.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 36.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 37.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 37.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 38.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 38.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 39.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 39.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 40.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 40.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 41.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 41.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 42.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 42.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 43.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 43.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 44.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 44.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 45.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 45.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 46.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 46.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 47.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 47.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 48.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 48.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 49.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 49.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 50.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 50.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 51.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 51.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 52.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 52.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 53.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 53.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 54.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 54.75;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 0.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 27.25;
y = 0.5;
z = 25.5*s;
-1
end
200
x = 27.5+(55-27.5)*s;
y = 0.5;
z = 8.75;
-1
end
200
3000
1.5
10
2
1000
0.1
0.5
0.1
Top or inlet
bottom
outlet
roof
wall
293.15
0
ustar = 0.41 * Uref /log((zref-H+z0)/z0);
k = ((ustar)^2)/(0.09^0.5);
epsilon = ((ustar)^3)/(0.41*(z-H+z0));
1
dir_x = sin(phi*pi/180);
dir_y = cos(phi*pi/180);
dir_z = 0;
ustar = k * Uref /log((zref-H+z0)/z0);
u_norm = ustar / k * log((z-H+z0)/z0);
0
1
1
0
0
0
0
0
0
0
0
0.1
0
0
0
0
0
1
0
0
0
0
0
0.1
1
1
0
1.204
0.01
1.83e-05
1017.24
0.02495
101325
293.15
0
0
-9.81
0
0
0
0.1
25
Symmetry
0
1
0
1
1
all[]
z>8.5
z > 0 and z < 8.5 and x > 27.5 and x < 55.0
0
S = 1 / volume;
dS = 0;
1
temperature = 293.15;
temperature = 293.15;
1e+12
0
1
SGDH
ustar = k * Uref /(log((zref-H+z0)/z0));
k = ((ustar)^2)/(0.09^0.5);
epsilon = ((ustar)^3)/(k*(z-H+z0));
1
0
1
0
1
3
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
velocity[0] = sin(phi*pi/180) * Uref * log((z-H+z0)/ z0) / log((zref - H+z0)/z0);
velocity[1] = cos(phi*pi/180) * Uref * log((z-H+z0)/ z0) / log((zref - H+z0)/z0);
velocity[2] = 0.;
2
1
1