In Fig. REF  we show the decay of the flow field perpendicular (\(v_\perp \) ) and parallel (\(v_{||}\) ) to the swimming direction,
which we compare to the decay in bulk using the analytic expressions [Eq. (REF )].
Due to the leading force dipole field, they decay with \(\sim r^{-2}\) .
For intermediate distances away from the squirmer the flow fields obtained from our simulations in strong confinement show a similar decay.
Indeed it has been shown analytically that for very strong confinement flow fields of microswimmers decay in general with \(\sim r^{-2}\)  [1]}, [2]}.
However, in our simulations the decay of the flow fields far away from the squirmer is faster, which probably stems from the periodic boundary conditions [3]}, [4]}.
Also, the fact that the strength of the flow field in front of the swimmer is weaker than at the back may be a signature of the periodic boundary conditions or due to the fact that MPCD does not model a perfect incompressible Stokes flow, but finite Reynolds and Schmidt numbers which account for finite fluid inertia and compressibility, respectively [5]}, [6]}.
<FIGURE><FIGURE>