I'm reading an article from Shiyi Chen et.al named "On boundary conditions in lattice Boltzmann methods". In the first section they proceed to the well known Chapman-Enskog expansion in order to derive the Navier-Stokes equation from the Lattice Boltzmann Equation (in the D2Q9 model):
$\frac{\partial f_i}{\partial t} + e_i.\nabla f_i = \Omega_i$, where $f_i$ is the velocity discretized distribution function, $e_i$ is a unit vector and $\Omega_i$ is the BGK collision operator. After they have proceeded to the Chapman expansion, they state that:
"If, in addition, a low Mach number (Ma) assumption is invoked as the nearly incompressible limit is approached, i.e.,
$Ma=u/C_s$ (where $C_s$ is the speed of sound) and $u' \sim \rho' \sim p' \sim O(Ma^2)$ [...] Here, the primed quantities denote fluctuations".
I don't really understand how do they come to $u' \sim \rho' \sim p' \sim O(Ma^2)$. Can someone explains me why do they have this relation ?
Have a good day.