Hi Mentz! Note that ##\xi^{a}## is hypersurface orthogonal, which by Frobenius' theorem implies ##\xi_{[a}\nabla_{b}\xi_{c]} = 0## which one can then use to show that the twist ##\omega_{ab}= 0## (geometrically, the lack of twisting is what allows us to use this foliation). Hence ##B_{ab} = B_{(ab)}## where ##B_{ab} = \nabla_{b}\xi_{a}## is the usual quantity from the kinematical decomposition. If the shear vanishes then ##B_{(ab)} = \frac{1}{3}h_{ab}\theta## where ##\theta = \nabla_{a}\xi^{a}## is the expansion. Since ##K_{ab} = B_{ba} = B_{(ab)}##, ##K_{ab}## vanishes identically iff ##\nabla_{a}\xi^{a} = 0## identically. This need not vanish in general.
As a simple counter example, consider the flat Friedman universe given in synchronous coordinates by ##ds^{2} = -d\tau^{2} + a^{2}(\tau)\{dx^{2} + dy^{2} + dz^{2}\}## where the function ##a(\tau)## to be solved for is called the scale factor. The geodesic congruence here corresponds to a family of locally inertial, isotropic observers comoving with the Hubble flow (the parameter ##\tau## is the proper time as read on a clock carried by anyone of these observers) and is given by their 4-velocity field ##u^{a}##. When represented in synchronous coordinates, this is just ##u^{\mu} = \delta^{\mu}_{\tau}## and ##\nabla_{\mu}u^{\mu} = \Gamma ^{\mu}_{\mu\tau} = \partial_{\tau}\ln\sqrt{\left | g \right |}\neq 0## in general.