all[]
u = 0.43;
v = 0;
w = 0;
101325
1
1273.15
0
RR = 1.0*0.025;
xr2 = RR^2;
xcen = 0.;
ycen = 0.;
xray2 = (x-xcen)^2 + (y-ycen)^2;
mesh_vi1 = 1;
if (xray2 < xr2) mesh_vi1 = 1e10;
20
1e-08
1
1000
0.001
1017.24
0.02495
0
0
0
0
0
0
10
11
2
3
4
0
0
0
0
0
0
0
0
0
1
Sr_a = 0.16;
U0 = 0.43;
c = 0.05;
a = 0.5*c;
mass = 1.;
f = Sr_a*U0/(2*a);
omega = 2*pi*f;
zz = a*sin(omega*t);
ww = a*omega*cos(omega*t);
acc = -a*omega^2*sin(omega*t);
F = acc*mass;
disp = a*sin(omega*t);
mesh_x=0;
mesh_y=disp;
mesh_z=0;
1
1
0.01
1000
0
1