def METHOD_NAME(self):
    # Setup walls
    for i in range(3):
        n = np.zeros(3)
        n[i] = 1
        self.lbf.add_boundary_from_shape(espressomd.shapes.Wall(
            normal=-n, dist=-(self.system.box_l[i] - AGRID)))
        self.lbf.add_boundary_from_shape(
            espressomd.shapes.Wall(normal=n, dist=AGRID))
    # setup sphere without slip in the middle
    sphere_shape = espressomd.shapes.Sphere(
        radius=RADIUS, center=self.system.box_l / 2, direction=1)
    self.lbf.add_boundary_from_shape(sphere_shape)
    sphere_volume = 4. / 3. * np.pi * RADIUS**3
    # Equilibration
    last_force = np.inf * np.ones(3)
    self.system.integrator.run(100)
    while True:
        self.system.integrator.run(10)
        force = np.linalg.norm(self.lbf.boundary['sphere'].get_force())
        if np.linalg.norm(force - last_force) < 0.01:
            break
        last_force = force
    # Check force balance
    boundary_force = np.zeros(3)
    for b in self.lbf.boundary:
        boundary_force += b.get_force()
    fluid_nodes = np.sum(np.logical_not(
        self.lbf[:, :, :].is_boundary).astype(int))
    fluid_volume = fluid_nodes * AGRID**3
    applied_force = fluid_volume * np.array(LB_PARAMS['ext_force_density'])
    np.testing.assert_allclose(
        boundary_force,
        applied_force,
        atol=0.08 * np.linalg.norm(applied_force))
    # Check buoyancy force on the sphere
    expected_force = np.array(
        [0, -sphere_volume * DENS * G, 0])
    np.testing.assert_allclose(
        np.copy(self.lbf.boundary['sphere'].get_force()), expected_force,
        atol=np.linalg.norm(expected_force) * 0.02)