[rans][oneeq] add density in all terms of spalart-amaras model

dumux/freeflow/rans/oneeq/staggered/fluxvariables.hh

@@ -96,7 +96,7 @@ public:
         // calculate advective flux
         auto upwindTermK = [](const auto& volVars)
         {
-            return volVars.viscosityTilde();
+            return volVars.viscosityTilde() * volVars.density();
         };
 
         flux[viscosityTildeEqIdx]
@@ -109,9 +109,9 @@ public:
         const auto& outsideVolVars = elemVolVars[scvf.outsideScvIdx()];
 
         // effective diffusion coefficients
-        Scalar insideCoeff = (insideVolVars.kinematicViscosity() + insideVolVars.viscosityTilde())
+        Scalar insideCoeff = (insideVolVars.viscosity() + insideVolVars.viscosityTilde() * insideVolVars.density())
                              / insideVolVars.sigma();
-        Scalar outsideCoeff = (outsideVolVars.kinematicViscosity() + outsideVolVars.viscosityTilde())
+        Scalar outsideCoeff = (outsideVolVars.viscosity() + outsideVolVars.viscosityTilde() * outsideVolVars.density())
                               / outsideVolVars.sigma();
 
         // scale by extrusion factor

dumux/freeflow/rans/oneeq/staggered/localresidual.hh

@@ -81,7 +81,7 @@ public:
                                                            const VolumeVariables& volVars) const
     {
         CellCenterPrimaryVariables storage = ParentType::computeStorageForCellCenter(problem, scv, volVars);
-        storage[viscosityTildeEqIdx] = volVars.viscosityTilde();
+        storage[viscosityTildeEqIdx] = volVars.viscosityTilde() * volVars.density();
         return storage;
     }
 
@@ -99,18 +99,19 @@ public:
 
         source[viscosityTildeEqIdx] += volVars.cb1() * (1.0 - volVars.ft2())
                                        * volVars.stressTensorScalarProductTilde()
-                                       * volVars.viscosityTilde();
+                                       * volVars.viscosityTilde() * volVars.density();
 
         source[viscosityTildeEqIdx] -= (volVars.cw1() * volVars.fW()
                                         - volVars.cb1() * volVars.ft2() / problem.karmanConstant() / problem.karmanConstant())
                                        * volVars.viscosityTilde() * volVars.viscosityTilde()
-                                       / volVars.wallDistance() / volVars.wallDistance();
+                                       / volVars.wallDistance() / volVars.wallDistance() * volVars.density();;
 
         for (unsigned int dimIdx = 0; dimIdx < ModelTraits::dim(); ++dimIdx)
         {
             source[viscosityTildeEqIdx] += volVars.cb2() / volVars.sigma()
                                            * volVars.storedViscosityTildeGradient()[dimIdx]
-                                           * volVars.storedViscosityTildeGradient()[dimIdx];
+                                           * volVars.storedViscosityTildeGradient()[dimIdx]
+                                           * volVars.density();
         }
 
         return source;