[lowrekepsilon] Cleanup, decrease verbosity of output

dumux/freeflow/rans/twoeq/lowrekepsilon/staggered/fluxvariables.hh

@@ -109,8 +109,6 @@ public:
             = ParentType::advectiveFluxForCellCenter(elemVolVars, elemFaceVars, scvf, upwindTermK, isOutflowK);
         flux[dissipationEqIdx - ModelTraits::dim()]
             = ParentType::advectiveFluxForCellCenter(elemVolVars, elemFaceVars, scvf, upwindTermEpsilon, isOutflowEpsilon);
-        Dune::dverb << " k_adv " << ParentType::advectiveFluxForCellCenter(elemVolVars, elemFaceVars, scvf, upwindTermK, isOutflowK);
-        Dune::dverb << " e_adv " << ParentType::advectiveFluxForCellCenter(elemVolVars, elemFaceVars, scvf, upwindTermEpsilon, isOutflowEpsilon);
 
         // calculate diffusive flux
         const auto& insideScv = fvGeometry.scv(scvf.insideScvIdx());
@@ -157,9 +155,6 @@ public:
                 += coeff_k / distance
                    * (insideVolVars.turbulentKineticEnergy() - outsideVolVars.turbulentKineticEnergy())
                    * scvf.area();
-            Dune::dverb << " k_diff " << coeff_k / distance
-                                                 * (insideVolVars.turbulentKineticEnergy() - outsideVolVars.turbulentKineticEnergy())
-                                                 * scvf.area();
         }
         if (!isOutflowEpsilon)
         {
@@ -167,9 +162,6 @@ public:
                 += coeff_e / distance
                    * (insideVolVars.dissipationTilde() - outsideVolVars.dissipationTilde())
                    * scvf.area();
-            Dune::dverb << " e_diff " << coeff_e / distance
-                                      * (insideVolVars.dissipationTilde() - outsideVolVars.dissipationTilde())
-                                      * scvf.area();
         }
         return flux;
     }

dumux/freeflow/rans/twoeq/lowrekepsilon/vtkoutputfields.hh

@@ -51,22 +51,8 @@ public:
     template <class VtkOutputModule>
     static void add(VtkOutputModule& vtk)
     {
-        vtk.addVolumeVariable([](const auto& v){ return v.turbulentKineticEnergy(); }, "turbulentKineticEnergy");
-        vtk.addVolumeVariable([](const auto& v){ return v.dissipationTilde(); }, "dissipation");
-        vtk.addVolumeVariable([](const auto& v){ return v.stressTensorScalarProduct(); }, "stressTensorScalarProduct");
-        vtk.addVolumeVariable([](const auto& v){ return v.kinematicEddyViscosity(); }, "eddyViscosity");
-        vtk.addVolumeVariable([](const auto& v){ return 2.0 * v.kinematicEddyViscosity()
-                                               * v.stressTensorScalarProduct(); }, "production_k");
-        vtk.addVolumeVariable([](const auto& volVars){ return volVars.cOneEpsilon() * volVars.fOne()
-                                    * volVars.dissipationTilde() / volVars.turbulentKineticEnergy()
-                                    * 2.0 * volVars.kinematicEddyViscosity()
-                                    * volVars.stressTensorScalarProduct(); }, "production_eps");
-        vtk.addVolumeVariable([](const auto& volVars){ return volVars.dissipationTilde(); }, "destruction_k");
-        vtk.addVolumeVariable([](const auto& volVars){ return volVars.cTwoEpsilon() * volVars.fTwo()
-                                    * volVars.dissipationTilde() * volVars.dissipationTilde()
-                                    / volVars.turbulentKineticEnergy(); }, "destruction_eps");
-        vtk.addVolumeVariable([](const auto& v){ return v.dValue(); }, "dValue");
-        vtk.addVolumeVariable([](const auto& v){ return v.eValue(); }, "eValue");
+        vtk.addVolumeVariable([](const auto& v){ return v.turbulentKineticEnergy(); }, "k");
+        vtk.addVolumeVariable([](const auto& v){ return v.dissipationTilde(); }, "epsilon");
     }
 };