b=[0.0, 1.0, 0.0],
c=[0.0, 0.0, boxZ],
direction=1)
boundaries.append(back_shape)
output_vtk_rhomboid(back_shape,
out_file=os.path.join(output_path, "wallBack.vtk"))
# obstacle - cylinder A
cylA_shape = espressomd.shapes.Cylinder(center=[11.0, 2.0, boxZ / 2.],
axis=[0.0, 0.0, 1.0],
length=boxZ,
radius=2.0,
direction=1)
boundaries.append(cylA_shape)
output_vtk_cylinder(cylA_shape,
n=20,
out_file=os.path.join(output_path, "cylinderA.vtk"))
# obstacle - cylinder B
cylB_shape = espressomd.shapes.Cylinder(center=[16.0, 8.0, boxZ / 2.],
axis=[0.0, 0.0, 1.0],
length=boxZ,
radius=2.0,
direction=1)
boundaries.append(cylB_shape)
output_vtk_cylinder(cylB_shape,
n=20,
out_file=os.path.join(output_path, "cylinderB.vtk"))
# obstacle - cylinder C
cylC_shape = espressomd.shapes.Cylinder(center=[11.0, 12.0, boxZ / 2.],
boundaries.append(front_shape)
output_vtk_rhomboid(
front_shape, out_file=output_path + "/wallFront.vtk")
# back wall of the channel
back_shape = shapes.Rhomboid(corner=[0.0, boxY - 1.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(back_shape)
output_vtk_rhomboid(
back_shape, out_file=output_path + "/wallBack.vtk")
# obstacle - cylinder A
cylA_shape = shapes.Cylinder(center=[11.0, 2.0, boxZ / 2.], axis=[0.0, 0.0, 1.0],
length=boxZ, radius=2.0, direction=1)
boundaries.append(cylA_shape)
output_vtk_cylinder(
cylA_shape, n=20, out_file=output_path + "/cylinderA.vtk")
# obstacle - cylinder B
cylB_shape = shapes.Cylinder(center=[16.0, 8.0, boxZ / 2.], axis=[0.0, 0.0, 1.0],
length=boxZ, radius=2.0, direction=1)
boundaries.append(cylB_shape)
output_vtk_cylinder(
cylB_shape, n=20, out_file=output_path + "/cylinderB.vtk")
# obstacle - cylinder C
cylC_shape = shapes.Cylinder(center=[11.0, 12.0, boxZ / 2.], axis=[0.0, 0.0, 1.0],
length=boxZ, radius=2.0, direction=1)
boundaries.append(cylC_shape)
output_vtk_cylinder(
cylC_shape, n=20, out_file=output_path + "/cylinderC.vtk")
b=[0.0, 1.0, 0.0],
c=[0.0, 0.0, boxZ],
direction=1)
boundaries.append(back_shape)
output_vtk_rhomboid(back_shape,
out_file="output/sim" + str(simNo) + "/wallBack.vtk")
# obstacle - cylinder A
cylA_shape = shapes.Cylinder(center=[11.0, 2.0, 7.0],
axis=[0.0, 0.0, 1.0],
length=7.0,
radius=2.0,
direction=1)
boundaries.append(cylA_shape)
output_vtk_cylinder(cylA_shape,
n=20,
out_file="output/sim" + str(simNo) + "/cylinderA.vtk")
# obstacle - cylinder B
cylB_shape = shapes.Cylinder(center=[16.0, 8.0, 7.0],
axis=[0.0, 0.0, 1.0],
length=7.0,
radius=2.0,
direction=1)
boundaries.append(cylB_shape)
output_vtk_cylinder(cylB_shape,
n=20,
out_file="output/sim" + str(simNo) + "/cylinderB.vtk")
# obstacle - cylinder C
cylC_shape = shapes.Cylinder(center=[11.0, 12.0, 7.0],