Solution (source code)

= Solution

Project the steady <Euler equations for an inviscid fluid> along a streamline. With no gravity,
$$
u\,du+\frac{dp}{\rho}=0.
$$
For an isentropic fluid, $dp=c_s^2d\rho$, so
$$
\frac{d\rho}{du}=-\frac{\rho u}{c_s^2}.
$$
The scalar mass flux is $j=\rho u$, and hence
$$
\boxed{\frac{dj}{du}=\rho\left(1-\frac{u^2}{c_s^2}\right)}.
$$
It increases with speed in <subsonic flow> and decreases with speed in <supersonic flow>.