def _get_fets_eval(self):
return FETS3D8H()
from mayavi import mlab
from simulator.api import \
TStepBC, XDomainFEGrid
from .mlab_decorators import decorate_figure
L = 600.0
H = 100.0
B = 50.0
L_c = 5.0
a = 5.0
w_max = 1
dgrid1 = XDomainFEGrid(dim_u=3,
coord_max=(L, H, B),
shape=(20, 5, 1),
fets=FETS3D8H())
x_x, x_y, y_z = dgrid1.mesh.geo_grid.point_x_grid
L_1 = x_x[1, 0]
d_L = L_c - L_1
x_x[1:, :] += d_L * (L - x_x[1:, :]) / (L - L_1)
a_H = a / H
n_a = int(a_H * dgrid1.shape[1])
fixed_right_bc = BCSlice(slice=dgrid1.mesh[-1, 0, -1, 0, :, :],
var='u',
dims=[1],
value=0)
fixed_x = BCSlice(slice=dgrid1.mesh[0, n_a:, :, 0, :, :],
var='u',
dims=[0],
value=0)
control_bc = BCSlice(slice=dgrid1.mesh[0, -1, :, 0, -1, :],