def METHOD_NAME(self):
    """
    Integrate the LB fluid until steady state is reached within a certain
    accuracy.
    """
    # disable periodicity except in the flow direction
    self.system.periodicity = np.logical_not(self.params['axis'])
    local_lb_params = LB_PARAMS.copy()
    local_lb_params['ext_force_density'] = np.array(
        self.params['axis']) * EXT_FORCE
    self.lbf = self.lb_class(**local_lb_params, **self.lb_params)
    self.system.lb = self.lbf
    cylinder_shape = espressomd.shapes.Cylinder(
        center=self.system.box_l / 2.0, axis=self.params['axis'],
        direction=-1, radius=EFFECTIVE_RADIUS, length=BOX_L * 1.5)
    self.lbf.add_boundary_from_shape(cylinder_shape)
    # simulate until profile converges
    mid_indices = 3 * [int((BOX_L / AGRID) / 2)]
    diff = float("inf")
    old_val = self.lbf[mid_indices].velocity[2]
    while diff > 1E-5:
        self.system.integrator.run(5)
        new_val = self.lbf[mid_indices].velocity[
            np.nonzero(self.params['axis'])[0]]
        diff = abs(new_val - old_val)
        old_val = new_val