diff --git a/dumux/discretization/cellcentered/tpfa/darcyslaw.hh b/dumux/discretization/cellcentered/tpfa/darcyslaw.hh
index 75867b75f72bdbbf1f0be95d152bf7cdb48f18ef..4f1945674864105861e03000064944a614496dfd 100644
--- a/dumux/discretization/cellcentered/tpfa/darcyslaw.hh
+++ b/dumux/discretization/cellcentered/tpfa/darcyslaw.hh
@@ -130,9 +130,9 @@ public:
             // do averaging for the density over all neighboring elements
             const auto rho = [&]()
             {
-                // boundaries
+                // use "Dirichlet" density on boundaries
                 if (scvf.boundary())
-                    return insideVolVars.density(phaseIdx);
+                    return outsideVolVars.density(phaseIdx);
 
                 // inner faces with two neighboring elements
                 else if (scvf.numOutsideScvs() == 1)