[rans][iofields] add writeFlatWallBoundedFields to the rans IO fields

dumux/common/parameters.hh

@@ -333,7 +333,6 @@ private:
         // parameters in the mpfa group
         defaultParams["MPFA.Q"] = "0.0";
 
-
         // merge the global default tree but do not overwrite if the parameter already exists
         mergeTree_(params, defaultParams, false);
     }

dumux/freeflow/rans/iofields.hh

@@ -40,7 +40,8 @@ struct RANSIOFields
     {
         NavierStokesIOFields::initOutputModule(out);
 
-        static const bool isFlatWallBounded = getParamFromGroup<bool>(out.paramGroup(), "RANS.IsFlatWallBounded");
+        static const bool isFlatWallBounded = getParamFromGroup<bool>(out.paramGroup(), "RANS.IsFlatWallBounded", false);
+        static const bool writeFlatWallBoundedFields = getParamFromGroup<bool>(out.paramGroup(), "RANS.WriteFlatWallBoundedFields", isFlatWallBounded);
 
         static constexpr auto dim = decltype(std::declval<typename OutputModule::VolumeVariables>().ccVelocityVector())::dimension;
 
@@ -54,7 +55,7 @@ struct RANSIOFields
         out.addVolumeVariable([](const auto& v){ return v.viscosity() / v.density(); }, "nu");
         out.addVolumeVariable([](const auto& v){ return v.kinematicEddyViscosity(); }, "nu_t");
         out.addVolumeVariable([](const auto& v){ return v.wallDistance(); }, "l_w");
-        if (isFlatWallBounded)
+        if (writeFlatWallBoundedFields)
         {
             out.addVolumeVariable([](const auto& v){ return v.yPlus(); }, "y^+");
             out.addVolumeVariable([](const auto& v){ return v.uPlus(); }, "u^+");

test/freeflow/rans/CMakeLists.txt

@@ -68,7 +68,9 @@ dumux_add_test(NAME test_ff_rans_lauferpipe_komega_nochannel
                         --files ${CMAKE_SOURCE_DIR}/test/references/test_ff_rans_lauferpipe_komega_nochannel-reference.vtu
                                 ${CMAKE_CURRENT_BINARY_DIR}/test_ff_rans_lauferpipe_komega_nochannel-00047.vtu
                         --command "${CMAKE_CURRENT_BINARY_DIR}/test_ff_rans_lauferpipe_komega params.input
-                        -Problem.Name test_ff_rans_lauferpipe_komega_nochannel")
+                        -Problem.Name test_ff_rans_lauferpipe_komega_nochannel
+                        -RANS.IsFlatWallBounded False
+                        -RANS.WriteFlatWallBoundedFields False")
 
 dumux_add_test(NAME test_ff_rans_lauferpipe_oneeq
               SOURCES main.cc

test/freeflow/rans/params_nonisothermal.input

@@ -23,6 +23,7 @@ EddyViscosityModel = "baldwinLomax"
 TurbulentPrandtlNumber = 0.85
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [Assembly]
 NumericDifferenceMethod = 0

test/freeflow/rans/params_nonisothermal_wallfunction.input

@@ -23,6 +23,7 @@ YPlusThreshold = 10. # should be small (10-30) for coarse grids
 TurbulentPrandtlNumber = 0.85
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [Assembly]
 NumericDifferenceMethod = 0

test/freeflow/rans/params_wallfunction.input

@@ -23,6 +23,7 @@ YPlusThreshold = 10. # should be small (10-30) for coarse grids
 
 [RANS]
 WallNormalAxis = 1
+WriteFlatWallBoundedFields = True
 
 [Assembly]
 NumericDifferenceMethod = 0

test/freeflow/rans/scripts/params_verbose.input

@@ -21,6 +21,7 @@ SandGrainRoughness = 0.0 # [m] # not implemented for EddyViscosityModel = 3
 EddyViscosityModel = "baldwinLomax"
 UseStoredEddyViscosity = false
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [KEpsilon]
 YPlusThreshold = 60. # should be large (30-60) for fine grids

test/freeflow/ransnc/params.input

@@ -23,6 +23,7 @@ EddyViscosityModel = "prandtl"
 TurbulentSchmidtNumber = 0.7
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [Assembly]
 NumericDifferenceMethod = 0

test/freeflow/ransnc/params_nonisothermal.input

@@ -21,6 +21,7 @@ TurbulentSchmidtNumber = 0.7
 TurbulentPrandtlNumber = 0.85
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [Assembly]
 NumericDifferenceMethod = 0

test/freeflow/ransnc/params_nonisothermal_wallfunction.input

@@ -19,6 +19,7 @@ TurbulentSchmidtNumber = 0.7
 TurbulentPrandtlNumber = 0.85
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [KEpsilon]
 YPlusThreshold = 10. # should be small (10-30) for coarse grids

test/freeflow/ransnc/params_wallfunction.input

@@ -22,6 +22,7 @@ TurbulentSchmidtNumber = 0.7
 TurbulentPrandtlNumber = 0.85
 WallNormalAxis = 1
 IsFlatWallBounded = True
+WriteFlatWallBoundedFields = True
 
 [KEpsilon]
 YPlusThreshold = 10. # should be small (10-30) for coarse grids