y0 = 0.;
y1 = 2.34;
x = 5.1;
y = y1*s + y0*(1.-s);
z = 0.;
-1
end
200
y0 = 0.;
y1 = 2.34;
x = 6.2;
y = y1*s + y0*(1.-s);
z = 0.;
-1
end
200
10
0.7
450000
0
0.001
cable
top
bottom
outlet
inlet
#exemple :
uref2 = 0.5;
dh = 0.13;
rho0 = 1025;
mu0 = 0.001;
re = sqrt(uref2)*dh*rho0/mu0;
if (re < 2000){
# in this case u*^2 is directly calculated to not have a problem with
# xlmbda=64/Re when Re->0
ustar2 = 8.*mu0*sqrt(uref2)/rho0/dh;}
else if (re<4000){
xlmbda = 0.021377 + 5.3115e-6*re;
ustar2 = uref2*xlmbda/8.;}
else {
xlmbda = 1/( 1.8*log(re)/log(10.)-1.64)^2;
ustar2 = uref2*xlmbda/8.;}
cmu = 0.09;
kappa = 0.42;
k = ustar2/sqrt(cmu);
epsilon = ustar2^1.5/(kappa*dh*0.1);
d2s3 = 2/3;
r11 = d2s3*k;
r22 = d2s3*k;
r33 = d2s3*k;
r12 = 0;
r13 = 0;
r23 = 0;
1
0.5
0
0
0
mesh_x=0.5*Um*T/pi*sin(2*pi*t/T);
mesh_y=0;
mesh_z=0;
0
0
0
0
0
0
0
0
0
4.3e-06
0
0
0
0.0001
standard
1
1
0
1025
molecular_viscosity = -1.;
0.001
1017.24
0.02495
0
0
0
0
0
0
0.1
25
front or back
0
0
0.5
1
1
all[]
0.5
0
1
PR=1.0*0.025;
xr2 = PR*2;
xcen = 0.0;
ycen = 0.;
xray2 = (x-xcen)^2 + (y-ycen)^2;
mesh_viscosity_1 = 1;
if (xray2 < xr2) mesh_viscosity_1 = 1e10;
1e-08
20
2
101325
293.15
0.12345
cmu = 0.09;
#example
trii = (0.02*0.1)^2;
cmu = 0.09;
r11 = trii;
r22 = trii;
r33 = trii;
r12 = 0.;
r13 = 0.;
r23 = 0.;
k = 0.5*(r11+r22+r33);
epsilon = k^1.5*cmu/0.1;
0
1
0
1
4
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
10
1
10