“…Then, motivated by the work in [6], we show that an improvement on the integrability of density can ensure the effectiveness and convergence of our approximation scheme. More specifically, we show that the uniform bound of ρ γ ln(1 + ρ) in L 1 , rather than the uniform bound of ρ γ+θ in L 1 for some θ > 0 as used in [10,9,21], ensures the vanishing of artificial pressure and the strong convergence of the density. To overcome the difficulty arising from the possible large oscillations of the density ρ, we adopt the method in Lions [21] and Feireisl [9] which is based on the celebrated weak continuity of the effective viscous flux p − (λ + 2µ)divu (see also Hoff [13]).…”