diff --git a/methods/conservedquantitycomputation.py b/methods/conservedquantitycomputation.py
index e6c049225396e7930be1d040f9a97cf774754d49..960e4ea169a8750351b8f1221131931931f207a3 100644
--- a/methods/conservedquantitycomputation.py
+++ b/methods/conservedquantitycomputation.py
@@ -149,10 +149,7 @@ class DensityVelocityComputation(AbstractConservedQuantityComputation):
         eqs = [sp.Eq(self._symbolOrder0, zerothOrderMoment)]
 
         firstOrderMoments = [a - b for a, b in zip(firstOrderMoments, velOffset)]
-        if self._compressible:
-            eqs += [sp.expand(sp.Eq(l, r * density)) for l, r in zip(self._symbolsOrder1, firstOrderMoments)]
-        else:
-            eqs += [sp.Eq(l, r) for l, r in zip(self._symbolsOrder1, firstOrderMoments)]
+        eqs += [sp.Eq(l, r) for l, r in zip(self._symbolsOrder1, firstOrderMoments)]
 
         return EquationCollection(eqs, [])