# 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],