# creating boundaries and obstacles in the channel
# OutputVtk writes a file
# lbboundaries created boundaries for fluid
# constraints created boundaries for the cells
boundaries = []
# bottom of the channel
bottom_shape = espressomd.shapes.Rhomboid(corner=[0.0, 0.0, 0.0],
a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0],
c=[0.0, 0.0, 1.0],
direction=1)
boundaries.append(bottom_shape)
output_vtk_rhomboid(bottom_shape,
out_file=os.path.join(output_path, "wallBottom.vtk"))
# top of the channel
top_shape = espressomd.shapes.Rhomboid(corner=[0.0, 0.0, boxZ - 1],
a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0],
c=[0.0, 0.0, 1.0],
direction=1)
boundaries.append(top_shape)
output_vtk_rhomboid(top_shape,
out_file=os.path.join(output_path, "wallTop.vtk"))
# front wall of the channel
front_shape = espressomd.shapes.Rhomboid(corner=[0.0, 0.0, 0.0],
a=[boxX, 0.0, 0.0],
b=[0.0, 1.0, 0.0],
system.actors.add(lbf)
system.thermostat.set_lb(LB_fluid=lbf, gamma=1.5)
# creating boundaries and obstacles in the channel
# OutputVtk writes a file
# lbboundaries created boundaries for fluid
# constraints created boundaries for the cells
boundaries = []
# bottom of the channel
bottom_shape = shapes.Rhomboid(corner=[0.0, 0.0, 0.0], a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0], c=[0.0, 0.0, 1.0],
direction=1)
boundaries.append(bottom_shape)
output_vtk_rhomboid(
bottom_shape, out_file=output_path + "/wallBottom.vtk")
# top of the channel
top_shape = shapes.Rhomboid(corner=[0.0, 0.0, boxZ - 1], a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0], c=[0.0, 0.0, 1.0], direction=1)
boundaries.append(top_shape)
output_vtk_rhomboid(
top_shape, out_file=output_path + "/wallTop.vtk")
# front wall of the channel
front_shape = shapes.Rhomboid(corner=[0.0, 0.0, 0.0], a=[boxX, 0.0, 0.0],
b=[0.0, 1.0, 0.0], c=[0.0, 0.0, boxZ], direction=1)
boundaries.append(front_shape)
output_vtk_rhomboid(
front_shape, out_file=output_path + "/wallFront.vtk")
# creating boundaries and obstacles in the channel
# OutputVtk writes a file
# lbboundaries created boundaries for fluid
# constraints created boundaries for the cells
boundaries = []
# bottom of the channel
bottom_shape = shapes.Rhomboid(corner=[0.0, 0.0, 0.0],
a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0],
c=[0.0, 0.0, 1.0],
direction=1)
boundaries.append(bottom_shape)
output_vtk_rhomboid(bottom_shape,
out_file="output/sim" + str(simNo) + "/wallBottom.vtk")
# top of the channel
top_shape = shapes.Rhomboid(corner=[0.0, 0.0, boxZ - 1],
a=[boxX, 0.0, 0.0],
b=[0.0, boxY, 0.0],
c=[0.0, 0.0, 1.0],
direction=1)
boundaries.append(top_shape)
output_vtk_rhomboid(top_shape,
out_file="output/sim" + str(simNo) + "/wallTop.vtk")
# front wall of the channel
front_shape = shapes.Rhomboid(corner=[0.0, 0.0, 0.0],
a=[boxX, 0.0, 0.0],
b=[0.0, 1.0, 0.0],
