Density–velocity equations with bulk modulus for computational hydro-acoustics