From accb459243367089905aa79306770be3fd0fb221 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Mon, 2 Jul 2018 16:47:13 +0200
Subject: [PATCH] [rans][twoeq] Implement turbulent dilation term

---
 .../twoeq/kepsilon/staggered/fluxvariables.hh | 18 ++++
 .../twoeq/komega/staggered/fluxvariables.hh   | 19 +++++
 .../lowrekepsilon/staggered/fluxvariables.hh  | 19 +++++
 test/references/pipe_laufer_kepsilon.vtu      | 44 +++++-----
 .../pipe_laufer_kepsilon_wallfunction.vtu     | 84 +++++++++----------
 test/references/pipe_laufer_komega.vtu        | 44 +++++-----
 test/references/pipe_laufer_lowrekepsilon.vtu | 44 +++++-----
 .../test_flatplate_kepsilon2cni.vtu           | 34 ++++----
 ...st_flatplate_kepsilon2cni_wallfunction.vtu | 10 +--
 test/references/test_flatplate_komega2cni.vtu | 34 ++++----
 .../test_flatplate_lowrekepsilon2cni.vtu      | 34 ++++----
 11 files changed, 220 insertions(+), 164 deletions(-)

diff --git a/dumux/freeflow/rans/twoeq/kepsilon/staggered/fluxvariables.hh b/dumux/freeflow/rans/twoeq/kepsilon/staggered/fluxvariables.hh
index 320161ba2f..1ee5672027 100644
--- a/dumux/freeflow/rans/twoeq/kepsilon/staggered/fluxvariables.hh
+++ b/dumux/freeflow/rans/twoeq/kepsilon/staggered/fluxvariables.hh
@@ -176,6 +176,24 @@ public:
         }
         return flux;
     }
+
+    /*!
+    * \brief Returns the momentum flux over all staggered faces.
+    */
+    FacePrimaryVariables computeMomentumFlux(const Problem& problem,
+                                             const Element& element,
+                                             const SubControlVolumeFace& scvf,
+                                             const FVElementGeometry& fvGeometry,
+                                             const ElementVolumeVariables& elemVolVars,
+                                             const ElementFaceVariables& elemFaceVars)
+    {
+        const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
+
+        return ParentType::computeFrontalMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + ParentType::computeLateralMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + 2.0 / ModelTraits::dim() * insideVolVars.density() * insideVolVars.turbulentKineticEnergy()
+                 * scvf.area() * scvf.directionSign() * insideVolVars.extrusionFactor();
+    }
 };
 
 } // end namespace
diff --git a/dumux/freeflow/rans/twoeq/komega/staggered/fluxvariables.hh b/dumux/freeflow/rans/twoeq/komega/staggered/fluxvariables.hh
index fd4645f064..2c723f196d 100644
--- a/dumux/freeflow/rans/twoeq/komega/staggered/fluxvariables.hh
+++ b/dumux/freeflow/rans/twoeq/komega/staggered/fluxvariables.hh
@@ -70,6 +70,7 @@ class KOmegaFluxVariablesImpl<TypeTag, BaseFluxVariables, DiscretizationMethod::
     using Indices = typename GET_PROP_TYPE(TypeTag, ModelTraits)::Indices;
     using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace;
     using CellCenterPrimaryVariables = typename GET_PROP_TYPE(TypeTag, CellCenterPrimaryVariables);
+    using FacePrimaryVariables = typename GET_PROP_TYPE(TypeTag, FacePrimaryVariables);
 
     static constexpr int turbulentKineticEnergyEqIdx = Indices::turbulentKineticEnergyEqIdx - ModelTraits::dim();
     static constexpr int dissipationEqIdx = Indices::dissipationEqIdx - ModelTraits::dim();
@@ -166,6 +167,24 @@ public:
         }
         return flux;
     }
+
+    /*!
+    * \brief Returns the momentum flux over all staggered faces.
+    */
+    FacePrimaryVariables computeMomentumFlux(const Problem& problem,
+                                             const Element& element,
+                                             const SubControlVolumeFace& scvf,
+                                             const FVElementGeometry& fvGeometry,
+                                             const ElementVolumeVariables& elemVolVars,
+                                             const ElementFaceVariables& elemFaceVars)
+    {
+        const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
+
+        return ParentType::computeFrontalMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + ParentType::computeLateralMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + 2.0 / ModelTraits::dim() * insideVolVars.density() * insideVolVars.turbulentKineticEnergy()
+                 * scvf.area() * scvf.directionSign() * insideVolVars.extrusionFactor();
+    }
 };
 
 } // end namespace
diff --git a/dumux/freeflow/rans/twoeq/lowrekepsilon/staggered/fluxvariables.hh b/dumux/freeflow/rans/twoeq/lowrekepsilon/staggered/fluxvariables.hh
index 489bab2140..6d8e276b92 100644
--- a/dumux/freeflow/rans/twoeq/lowrekepsilon/staggered/fluxvariables.hh
+++ b/dumux/freeflow/rans/twoeq/lowrekepsilon/staggered/fluxvariables.hh
@@ -70,6 +70,7 @@ class LowReKEpsilonFluxVariablesImpl<TypeTag, BaseFluxVariables, DiscretizationM
     using Indices = typename GET_PROP_TYPE(TypeTag, ModelTraits)::Indices;
     using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace;
     using CellCenterPrimaryVariables = typename GET_PROP_TYPE(TypeTag, CellCenterPrimaryVariables);
+    using FacePrimaryVariables = typename GET_PROP_TYPE(TypeTag, FacePrimaryVariables);
 
     static constexpr int turbulentKineticEnergyEqIdx = Indices::turbulentKineticEnergyEqIdx - ModelTraits::dim();
     static constexpr int dissipationEqIdx = Indices::dissipationEqIdx - ModelTraits::dim();
@@ -167,6 +168,24 @@ public:
         }
         return flux;
     }
+
+    /*!
+    * \brief Returns the momentum flux over all staggered faces.
+    */
+    FacePrimaryVariables computeMomentumFlux(const Problem& problem,
+                                             const Element& element,
+                                             const SubControlVolumeFace& scvf,
+                                             const FVElementGeometry& fvGeometry,
+                                             const ElementVolumeVariables& elemVolVars,
+                                             const ElementFaceVariables& elemFaceVars)
+    {
+        const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
+
+        return ParentType::computeFrontalMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + ParentType::computeLateralMomentumFlux(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars)
+               + 2.0 / ModelTraits::dim() * insideVolVars.density() * insideVolVars.turbulentKineticEnergy()
+                 * scvf.area() * scvf.directionSign() * insideVolVars.extrusionFactor();
+    }
 };
 
 } // end namespace
diff --git a/test/references/pipe_laufer_kepsilon.vtu b/test/references/pipe_laufer_kepsilon.vtu
index c9011c3cab..7852291bd2 100644
--- a/test/references/pipe_laufer_kepsilon.vtu
+++ b/test/references/pipe_laufer_kepsilon.vtu
@@ -100,28 +100,28 @@
           0.394357 0.396274 0.398377 0.399961
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          2.39118 2.149 1.97741 1.80969 1.64451 1.48447 1.33036 1.18165 1.03735 0.896392 0.757719 0.620388
-          0.483555 0.346533 0.208838 0.0699234 2.39131 2.14926 1.97756 1.80978 1.64457 1.48451 1.33038 1.18166
-          1.03736 0.896392 0.757716 0.620384 0.48355 0.346528 0.208833 0.0699209 2.3911 2.14952 1.97771 1.80989
-          1.64464 1.48456 1.33041 1.18168 1.03737 0.896397 0.757716 0.620381 0.483546 0.346522 0.208828 0.0699177
-          2.39062 2.14971 1.9778 1.80998 1.64472 1.48462 1.33046 1.18171 1.03739 0.896411 0.757723 0.620383
-          0.483544 0.346518 0.208822 0.0699132 2.3899 2.14991 1.97784 1.81006 1.64482 1.48472 1.33054 1.18177
-          1.03744 0.896444 0.757745 0.620395 0.483549 0.346518 0.20882 0.069907 2.38904 2.15013 1.97786 1.8101
-          1.6449 1.48482 1.33064 1.18187 1.03752 0.896508 0.757794 0.62043 0.483572 0.346531 0.208825 0.0698985
-          2.38814 2.15037 1.97786 1.81009 1.64493 1.48488 1.33073 1.18197 1.03762 0.896605 0.757878 0.620497
-          0.483622 0.346566 0.208849 0.0698875 2.38743 2.15059 1.97787 1.81008 1.64492 1.48489 1.33075 1.18201
-          1.03768 0.896675 0.757954 0.620573 0.483691 0.346625 0.208895 0.0698781 2.38743 2.15059 1.97787 1.81008
-          1.64492 1.48489 1.33075 1.18201 1.03768 0.896675 0.757954 0.620573 0.483691 0.346625 0.208895 0.0698781
-          2.38814 2.15037 1.97786 1.81009 1.64493 1.48488 1.33073 1.18197 1.03762 0.896605 0.757878 0.620497
-          0.483622 0.346566 0.208849 0.0698875 2.38904 2.15013 1.97786 1.8101 1.6449 1.48482 1.33064 1.18187
-          1.03752 0.896508 0.757794 0.62043 0.483572 0.346531 0.208825 0.0698985 2.3899 2.14991 1.97784 1.81006
-          1.64482 1.48472 1.33054 1.18177 1.03744 0.896444 0.757745 0.620395 0.483549 0.346518 0.20882 0.069907
-          2.39062 2.14971 1.9778 1.80998 1.64472 1.48462 1.33046 1.18171 1.03739 0.896411 0.757723 0.620383
-          0.483544 0.346518 0.208822 0.0699132 2.3911 2.14952 1.97771 1.80989 1.64464 1.48456 1.33041 1.18168
-          1.03737 0.896397 0.757716 0.620381 0.483546 0.346522 0.208828 0.0699177 2.39131 2.14926 1.97756 1.80978
-          1.64457 1.48451 1.33038 1.18166 1.03736 0.896392 0.757716 0.620384 0.48355 0.346528 0.208833 0.0699209
-          2.39118 2.149 1.97741 1.80969 1.64451 1.48447 1.33036 1.18165 1.03735 0.896392 0.757719 0.620388
-          0.483555 0.346533 0.208838 0.0699234
+          2.40991 2.1777 2.02219 1.86754 1.71053 1.55527 1.40402 1.25719 1.11427 0.974409 0.836659 0.700126
+          0.563974 0.427512 0.290253 0.151657 2.41004 2.17797 2.02241 1.86784 1.71096 1.55588 1.40483 1.2582
+          1.11545 0.975712 0.838047 0.701557 0.565415 0.428938 0.291648 0.153024 2.40905 2.18168 2.02678 1.87093
+          1.71297 1.55735 1.40612 1.25948 1.11678 0.97709 0.839456 0.702971 0.56681 0.430294 0.292953 0.154287
+          2.40736 2.17214 2.01217 1.85687 1.70113 1.5475 1.39772 1.25206 1.11003 0.970799 0.833468 0.697181
+          0.56114 0.424686 0.287361 0.148706 2.40614 2.16627 1.9986 1.83899 1.6833 1.53165 1.38399 1.24007
+          1.09935 0.961066 0.824407 0.688575 0.552827 0.416538 0.279284 0.140667 2.40522 2.16416 1.99179 1.82621
+          1.66566 1.5121 1.36503 1.22275 1.08374 0.946897 0.811373 0.676375 0.541199 0.40527 0.268213 0.129708
+          2.40439 2.16383 1.98963 1.8212 1.65646 1.4981 1.34695 1.20235 1.06282 0.926654 0.792283 0.658442
+          0.524182 0.38891 0.252306 0.114102 2.40372 2.16399 1.98923 1.82001 1.65392 1.49342 1.33933 1.19118
+          1.04802 0.908754 0.772211 0.637248 0.502792 0.367951 0.232118 0.0946715 2.40372 2.16399 1.98923 1.82001
+          1.65392 1.49342 1.33933 1.19118 1.04802 0.908754 0.772211 0.637248 0.502792 0.367951 0.232118 0.0946715
+          2.40439 2.16383 1.98963 1.8212 1.65646 1.4981 1.34695 1.20235 1.06282 0.926654 0.792283 0.658442
+          0.524182 0.38891 0.252306 0.114102 2.40522 2.16416 1.99179 1.82621 1.66566 1.5121 1.36503 1.22275
+          1.08374 0.946897 0.811373 0.676375 0.541199 0.40527 0.268213 0.129708 2.40614 2.16627 1.9986 1.83899
+          1.6833 1.53165 1.38399 1.24007 1.09935 0.961066 0.824407 0.688575 0.552827 0.416538 0.279284 0.140667
+          2.40736 2.17214 2.01217 1.85687 1.70113 1.5475 1.39772 1.25206 1.11003 0.970799 0.833468 0.697181
+          0.56114 0.424686 0.287361 0.148706 2.40905 2.18168 2.02678 1.87093 1.71297 1.55735 1.40612 1.25948
+          1.11678 0.97709 0.839456 0.702971 0.56681 0.430294 0.292953 0.154287 2.41004 2.17797 2.02241 1.86784
+          1.71096 1.55588 1.40483 1.2582 1.11545 0.975712 0.838047 0.701557 0.565415 0.428938 0.291648 0.153024
+          2.40991 2.1777 2.02219 1.86754 1.71053 1.55527 1.40402 1.25719 1.11427 0.974409 0.836659 0.700126
+          0.563974 0.427512 0.290253 0.151657
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.43725e-05 1.43725e-05 1.43726e-05 1.43726e-05 1.43726e-05 1.43726e-05 1.43726e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05
diff --git a/test/references/pipe_laufer_kepsilon_wallfunction.vtu b/test/references/pipe_laufer_kepsilon_wallfunction.vtu
index 944d369012..3a511558eb 100644
--- a/test/references/pipe_laufer_kepsilon_wallfunction.vtu
+++ b/test/references/pipe_laufer_kepsilon_wallfunction.vtu
@@ -180,48 +180,48 @@
           0.683226 0.678068 0.673622 0.669819 0.666566 0.663733 0.66118 0.659363
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          1.31886 1.26628 1.21589 1.16659 1.11807 1.06974 1.021 0.971335 0.920488 0.868392 0.815126 0.760846
-          0.705738 0.64999 0.59377 0.537223 0.480469 0.423605 0.366708 0.309839 0.253053 0.196404 0.13996 0.0837914
-          0.0278886 1.31855 1.26633 1.21595 1.16665 1.11813 1.06981 1.02106 0.971398 0.920548 0.868448 0.815177
-          0.760893 0.705782 0.65003 0.593806 0.537255 0.480498 0.423631 0.366731 0.30986 0.253072 0.196424 0.139982
-          0.083816 0.0278978 1.31815 1.26634 1.21596 1.16666 1.11814 1.06982 1.02107 0.97141 0.920561 0.868462
-          0.815194 0.760912 0.705802 0.65005 0.593827 0.537276 0.480518 0.42365 0.366749 0.309877 0.253089 0.196441
-          0.14 0.0838363 0.0279062 1.3178 1.26633 1.21595 1.16666 1.11814 1.06982 1.02107 0.971407 0.920561
-          0.868468 0.815203 0.760925 0.705818 0.650069 0.593846 0.537296 0.480539 0.423671 0.366769 0.309897 0.253108
-          0.19646 0.140019 0.0838542 0.0279154 1.31751 1.26632 1.21595 1.16665 1.11813 1.0698 1.02106 0.971394
-          0.920552 0.868463 0.815203 0.760929 0.705826 0.650081 0.593862 0.537314 0.480558 0.42369 0.366789 0.309917
-          0.253128 0.196479 0.140038 0.0838718 0.027927 1.31727 1.26632 1.21595 1.16665 1.11812 1.06979 1.02104
-          0.97138 0.920538 0.868451 0.815195 0.760925 0.705827 0.650085 0.59387 0.537325 0.480572 0.423706 0.366806
-          0.309934 0.253146 0.196497 0.140056 0.0838889 0.0279405 1.31708 1.26632 1.21595 1.16665 1.11812 1.06979
-          1.02103 0.971368 0.920525 0.868438 0.815183 0.760916 0.70582 0.650082 0.59387 0.537329 0.480579 0.423716
-          0.366818 0.309948 0.253161 0.196513 0.140072 0.0839045 0.0279551 1.31693 1.26633 1.21596 1.16665 1.11812
-          1.06978 1.02103 0.97136 0.920516 0.868428 0.815173 0.760906 0.705812 0.650075 0.593866 0.537327 0.480579
-          0.423719 0.366824 0.309956 0.253171 0.196525 0.140084 0.0839172 0.0279691 1.31684 1.26633 1.21596 1.16665
-          1.11812 1.06978 1.02102 0.971356 0.920511 0.868423 0.815167 0.760899 0.705805 0.650069 0.59386 0.537322
-          0.480576 0.423718 0.366824 0.309959 0.253175 0.196531 0.140092 0.0839256 0.0279801 1.31679 1.26633 1.21596
-          1.16665 1.11812 1.06978 1.02102 0.971354 0.920509 0.86842 0.815164 0.760897 0.705802 0.650065 0.593856
-          0.537318 0.480572 0.423715 0.366822 0.309958 0.253175 0.196532 0.140095 0.083929 0.0279859 1.31679 1.26633
-          1.21596 1.16665 1.11812 1.06978 1.02102 0.971354 0.920509 0.86842 0.815164 0.760897 0.705802 0.650065
-          0.593856 0.537318 0.480572 0.423715 0.366822 0.309958 0.253175 0.196532 0.140095 0.083929 0.0279859 1.31684
-          1.26633 1.21596 1.16665 1.11812 1.06978 1.02102 0.971356 0.920511 0.868423 0.815167 0.760899 0.705805
-          0.650069 0.59386 0.537322 0.480576 0.423718 0.366824 0.309959 0.253175 0.196531 0.140092 0.0839256 0.0279801
-          1.31693 1.26633 1.21596 1.16665 1.11812 1.06978 1.02103 0.97136 0.920516 0.868428 0.815173 0.760906
-          0.705812 0.650075 0.593866 0.537327 0.480579 0.423719 0.366824 0.309956 0.253171 0.196525 0.140084 0.0839172
-          0.0279691 1.31708 1.26632 1.21595 1.16665 1.11812 1.06979 1.02103 0.971368 0.920525 0.868438 0.815183
-          0.760916 0.70582 0.650082 0.59387 0.537329 0.480579 0.423716 0.366818 0.309948 0.253161 0.196513 0.140072
-          0.0839045 0.0279551 1.31727 1.26632 1.21595 1.16665 1.11812 1.06979 1.02104 0.97138 0.920538 0.868451
-          0.815195 0.760925 0.705827 0.650085 0.59387 0.537325 0.480572 0.423706 0.366806 0.309934 0.253146 0.196497
-          0.140056 0.0838889 0.0279405 1.31751 1.26632 1.21595 1.16665 1.11813 1.0698 1.02106 0.971394 0.920552
-          0.868463 0.815203 0.760929 0.705826 0.650081 0.593862 0.537314 0.480558 0.42369 0.366789 0.309917 0.253128
-          0.196479 0.140038 0.0838718 0.027927 1.3178 1.26633 1.21595 1.16666 1.11814 1.06982 1.02107 0.971407
-          0.920561 0.868468 0.815203 0.760925 0.705818 0.650069 0.593846 0.537296 0.480539 0.423671 0.366769 0.309897
-          0.253108 0.19646 0.140019 0.0838542 0.0279154 1.31815 1.26634 1.21596 1.16666 1.11814 1.06982 1.02107
-          0.97141 0.920561 0.868462 0.815194 0.760912 0.705802 0.65005 0.593827 0.537276 0.480518 0.42365 0.366749
-          0.309877 0.253089 0.196441 0.14 0.0838363 0.0279062 1.31855 1.26633 1.21595 1.16665 1.11813 1.06981
-          1.02106 0.971398 0.920548 0.868448 0.815177 0.760893 0.705782 0.65003 0.593806 0.537255 0.480498 0.423631
-          0.366731 0.30986 0.253072 0.196424 0.139982 0.083816 0.0278978 1.31886 1.26628 1.21589 1.16659 1.11807
-          1.06974 1.021 0.971335 0.920488 0.868392 0.815126 0.760846 0.705738 0.64999 0.59377 0.537223 0.480469
-          0.423605 0.366708 0.309839 0.253053 0.196404 0.13996 0.0837914 0.0278886
+          1.33828 1.28455 1.23331 1.18365 1.13539 1.08798 1.04073 0.992993 0.944348 0.894579 0.843655 0.791658
+          0.73873 0.685032 0.630722 0.575942 0.520816 0.465448 0.409927 0.354331 0.298739 0.243244 0.187969 0.133065
+          0.0785534 1.33523 1.28118 1.2307 1.18261 1.13621 1.0905 1.04452 0.997555 0.949236 0.899472 0.848349
+          0.796045 0.742762 0.688699 0.634036 0.578926 0.523498 0.467859 0.412101 0.356306 0.300563 0.24498 0.189707
+          0.13491 0.0805312 1.33511 1.28078 1.22894 1.17913 1.13087 1.08346 1.03612 0.988209 0.939322 0.889279
+          0.838075 0.785815 0.732653 0.678758 0.624291 0.569392 0.51418 0.458755 0.403197 0.347579 0.291971 0.236452
+          0.181128 0.126113 0.0713737 1.33502 1.28099 1.22877 1.17827 1.12909 1.0806 1.03211 0.983086 0.933225
+          0.882403 0.830631 0.777997 0.724621 0.670634 0.616164 0.561323 0.506212 0.450912 0.395493 0.340013 0.28453
+          0.229106 0.173821 0.118761 0.0638905 1.33487 1.2812 1.2289 1.17817 1.12862 1.07963 1.03053 0.980803
+          0.930175 0.878569 0.826051 0.772756 0.718837 0.664437 0.609675 0.554651 0.499442 0.44411 0.388707 0.333276
+          0.277859 0.222504 0.167276 0.112244 0.0573576 1.3347 1.28133 1.22905 1.17825 1.12855 1.07935 1.02996
+          0.979886 0.928832 0.876728 0.823653 0.769764 0.715248 0.660284 0.605023 0.549584 0.494052 0.438485 0.382923
+          0.327396 0.271931 0.216561 0.161339 0.106318 0.0514325 1.33454 1.28141 1.22915 1.17834 1.1286 1.0793
+          1.02978 0.979536 0.928269 0.875904 0.822513 0.76825 0.713307 0.657874 0.602126 0.546206 0.490227 0.434267
+          0.378377 0.322591 0.266932 0.211425 0.156109 0.10103 0.0461037 1.33441 1.28145 1.22922 1.17841 1.12865
+          1.07932 1.02974 0.979413 0.92804 0.875543 0.821988 0.767526 0.712341 0.656621 0.60054 0.544249 0.487874
+          0.431513 0.375237 0.319098 0.263133 0.207372 0.151858 0.0966296 0.0415992 1.33432 1.28148 1.22926 1.17846
+          1.1287 1.07935 1.02974 0.979377 0.927953 0.87539 0.821754 0.767192 0.711885 0.656017 0.599756 0.543252
+          0.486629 0.42999 0.373415 0.31697 0.260706 0.204672 0.148919 0.093496 0.0383225 1.33428 1.28148 1.22927
+          1.17848 1.12872 1.07936 1.02975 0.97937 0.927925 0.875335 0.821665 0.767062 0.711704 0.655774 0.599438
+          0.542843 0.486111 0.429343 0.372622 0.316015 0.259581 0.203378 0.147468 0.0919084 0.03663 1.33428 1.28148
+          1.22927 1.17848 1.12872 1.07936 1.02975 0.97937 0.927925 0.875335 0.821665 0.767062 0.711704 0.655774
+          0.599438 0.542843 0.486111 0.429343 0.372622 0.316015 0.259581 0.203378 0.147468 0.0919084 0.03663 1.33432
+          1.28148 1.22926 1.17846 1.1287 1.07935 1.02974 0.979377 0.927953 0.87539 0.821754 0.767192 0.711885
+          0.656017 0.599756 0.543252 0.486629 0.42999 0.373415 0.31697 0.260706 0.204672 0.148919 0.093496 0.0383225
+          1.33441 1.28145 1.22922 1.17841 1.12865 1.07932 1.02974 0.979413 0.92804 0.875543 0.821988 0.767526
+          0.712341 0.656621 0.60054 0.544249 0.487874 0.431513 0.375237 0.319098 0.263133 0.207372 0.151858 0.0966296
+          0.0415992 1.33454 1.28141 1.22915 1.17834 1.1286 1.0793 1.02978 0.979536 0.928269 0.875904 0.822513
+          0.76825 0.713307 0.657874 0.602126 0.546206 0.490227 0.434267 0.378377 0.322591 0.266932 0.211425 0.156109
+          0.10103 0.0461037 1.3347 1.28133 1.22905 1.17825 1.12855 1.07935 1.02996 0.979886 0.928832 0.876728
+          0.823653 0.769764 0.715248 0.660284 0.605023 0.549584 0.494052 0.438485 0.382923 0.327396 0.271931 0.216561
+          0.161339 0.106318 0.0514325 1.33487 1.2812 1.2289 1.17817 1.12862 1.07963 1.03053 0.980803 0.930175
+          0.878569 0.826051 0.772756 0.718837 0.664437 0.609675 0.554651 0.499442 0.44411 0.388707 0.333276 0.277859
+          0.222504 0.167276 0.112244 0.0573576 1.33502 1.28099 1.22877 1.17827 1.12909 1.0806 1.03211 0.983086
+          0.933225 0.882403 0.830631 0.777997 0.724621 0.670634 0.616164 0.561323 0.506212 0.450912 0.395493 0.340013
+          0.28453 0.229106 0.173821 0.118761 0.0638905 1.33511 1.28078 1.22894 1.17913 1.13087 1.08346 1.03612
+          0.988209 0.939322 0.889279 0.838075 0.785815 0.732653 0.678758 0.624291 0.569392 0.51418 0.458755 0.403197
+          0.347579 0.291971 0.236452 0.181128 0.126113 0.0713737 1.33523 1.28118 1.2307 1.18261 1.13621 1.0905
+          1.04452 0.997555 0.949236 0.899472 0.848349 0.796045 0.742762 0.688699 0.634036 0.578926 0.523498 0.467859
+          0.412101 0.356306 0.300563 0.24498 0.189707 0.13491 0.0805312 1.33828 1.28455 1.23331 1.18365 1.13539
+          1.08798 1.04073 0.992993 0.944348 0.894579 0.843655 0.791658 0.73873 0.685032 0.630722 0.575942 0.520816
+          0.465448 0.409927 0.354331 0.298739 0.243244 0.187969 0.133065 0.0785534
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.43726e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05 1.43727e-05
diff --git a/test/references/pipe_laufer_komega.vtu b/test/references/pipe_laufer_komega.vtu
index cb6911f932..ecd36ff582 100644
--- a/test/references/pipe_laufer_komega.vtu
+++ b/test/references/pipe_laufer_komega.vtu
@@ -100,28 +100,28 @@
           0.160624 0.160739 0.160808 0.160845
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          10.3546 9.31494 8.55387 7.89439 7.27651 6.66926 6.05838 5.43919 4.81182 4.1781 3.53997 2.89896
-          2.25615 1.61223 0.967644 0.32264 10.3565 9.31572 8.55405 7.89442 7.27651 6.66925 6.05837 5.43918
-          4.81181 4.1781 3.53997 2.89896 2.25615 1.61223 0.967643 0.32264 10.3555 9.31679 8.55435 7.8945
-          7.27652 6.66924 6.05837 5.43917 4.8118 4.17809 3.53996 2.89896 2.25615 1.61223 0.967643 0.32264
-          10.3539 9.31802 8.55482 7.89468 7.27659 6.66926 6.05836 5.43916 4.8118 4.17808 3.53996 2.89896
-          2.25615 1.61223 0.967642 0.322639 10.3516 9.31944 8.55554 7.89503 7.27674 6.66932 6.05838 5.43916
-          4.81179 4.17808 3.53995 2.89895 2.25615 1.61223 0.967641 0.322639 10.3486 9.321 8.55646 7.89558
-          7.27705 6.66947 6.05845 5.43919 4.8118 4.17808 3.53995 2.89895 2.25614 1.61223 0.96764 0.322638
-          10.3452 9.32263 8.55747 7.89629 7.27751 6.66975 6.0586 5.43927 4.81184 4.1781 3.53996 2.89895
-          2.25615 1.61223 0.96764 0.322637 10.3425 9.32402 8.55838 7.897 7.27805 6.67013 6.05884 5.43941
-          4.81192 4.17815 3.53999 2.89897 2.25616 1.61224 0.967645 0.322637 10.3425 9.32402 8.55838 7.897
-          7.27805 6.67013 6.05884 5.43941 4.81192 4.17815 3.53999 2.89897 2.25616 1.61224 0.967645 0.322637
-          10.3452 9.32263 8.55747 7.89629 7.27751 6.66975 6.0586 5.43927 4.81184 4.1781 3.53996 2.89895
-          2.25615 1.61223 0.96764 0.322637 10.3486 9.321 8.55646 7.89558 7.27705 6.66947 6.05845 5.43919
-          4.8118 4.17808 3.53995 2.89895 2.25614 1.61223 0.96764 0.322638 10.3516 9.31944 8.55554 7.89503
-          7.27674 6.66932 6.05838 5.43916 4.81179 4.17808 3.53995 2.89895 2.25615 1.61223 0.967641 0.322639
-          10.3539 9.31802 8.55482 7.89468 7.27659 6.66926 6.05836 5.43916 4.8118 4.17808 3.53996 2.89896
-          2.25615 1.61223 0.967642 0.322639 10.3555 9.31679 8.55435 7.8945 7.27652 6.66924 6.05837 5.43917
-          4.8118 4.17809 3.53996 2.89896 2.25615 1.61223 0.967643 0.32264 10.3565 9.31572 8.55405 7.89442
-          7.27651 6.66925 6.05837 5.43918 4.81181 4.1781 3.53997 2.89896 2.25615 1.61223 0.967643 0.32264
-          10.3546 9.31494 8.55387 7.89439 7.27651 6.66926 6.05838 5.43919 4.81182 4.1781 3.53997 2.89896
-          2.25615 1.61223 0.967644 0.32264
+          10.6453 9.55531 8.77428 8.10983 7.49347 6.88963 6.28215 5.66564 5.04017 4.40769 3.77035 3.12982
+          2.48729 1.84354 1.19905 0.554104 10.5361 9.66225 8.9024 8.23953 7.62493 7.02338 6.41801 5.80312
+          5.17877 4.54703 3.91014 3.26989 2.62753 1.98388 1.33945 0.694535 10.4618 9.60949 8.88364 8.23139
+          7.61981 7.01914 6.41395 5.799 5.17454 4.54269 3.90571 3.26541 2.62301 1.97933 1.33489 0.689963
+          10.4181 9.53176 8.83575 8.19932 7.59321 6.99407 6.38901 5.77375 5.14889 4.51671 3.87949 3.23902
+          2.59652 1.95278 1.30829 0.663341 10.3907 9.45505 8.77015 8.15046 7.55299 6.95719 6.35303 5.7377
+          5.11248 4.47991 3.84237 3.20167 2.55901 1.91516 1.27061 0.625627 10.3736 9.39272 8.69341 8.08255
+          7.49458 6.90409 6.30211 5.68737 5.06207 4.42922 3.7914 3.15045 2.50761 1.86364 1.21901 0.573977
+          10.3639 9.35465 8.62294 8.00229 7.41658 6.83034 6.23138 5.61835 4.99394 4.36151 3.72387 3.08299
+          2.44016 1.79618 1.15152 0.506485 10.3594 9.3409 8.58174 7.93463 7.33432 6.74432 6.1472 5.53762
+          4.91645 4.28653 3.65065 3.01094 2.36887 1.72536 1.081 0.436131 10.3594 9.3409 8.58174 7.93463
+          7.33432 6.74432 6.1472 5.53762 4.91645 4.28653 3.65065 3.01094 2.36887 1.72536 1.081 0.436131
+          10.3639 9.35465 8.62294 8.00229 7.41658 6.83034 6.23138 5.61835 4.99394 4.36151 3.72387 3.08299
+          2.44016 1.79618 1.15152 0.506485 10.3736 9.39272 8.69341 8.08255 7.49458 6.90409 6.30211 5.68737
+          5.06207 4.42922 3.7914 3.15045 2.50761 1.86364 1.21901 0.573977 10.3907 9.45505 8.77015 8.15046
+          7.55299 6.95719 6.35303 5.7377 5.11248 4.47991 3.84237 3.20167 2.55901 1.91516 1.27061 0.625627
+          10.4181 9.53176 8.83575 8.19932 7.59321 6.99407 6.38901 5.77375 5.14889 4.51671 3.87949 3.23902
+          2.59652 1.95278 1.30829 0.663341 10.4618 9.60949 8.88364 8.23139 7.61981 7.01914 6.41395 5.799
+          5.17454 4.54269 3.90571 3.26541 2.62301 1.97933 1.33489 0.689963 10.5361 9.66225 8.9024 8.23953
+          7.62493 7.02338 6.41801 5.80312 5.17877 4.54703 3.91014 3.26989 2.62753 1.98388 1.33945 0.694535
+          10.6453 9.55531 8.77428 8.10983 7.49347 6.88963 6.28215 5.66564 5.04017 4.40769 3.77035 3.12982
+          2.48729 1.84354 1.19905 0.554104
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.43713e-05 1.43715e-05 1.43716e-05 1.43717e-05 1.43718e-05 1.43719e-05 1.4372e-05 1.43721e-05 1.43721e-05 1.43722e-05 1.43723e-05 1.43724e-05
diff --git a/test/references/pipe_laufer_lowrekepsilon.vtu b/test/references/pipe_laufer_lowrekepsilon.vtu
index 9d9cc6aeef..e41e2f05ef 100644
--- a/test/references/pipe_laufer_lowrekepsilon.vtu
+++ b/test/references/pipe_laufer_lowrekepsilon.vtu
@@ -100,28 +100,28 @@
           0.302786 0.304296 0.305596 0.306453
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          3.63626 3.43035 3.22516 3.00849 2.7864 2.5609 2.33281 2.10242 1.86953 1.63356 1.39373 1.14937
-          0.900345 0.647057 0.390256 0.130578 3.63616 3.43063 3.22537 3.00861 2.78647 2.56094 2.33284 2.10244
-          1.86954 1.63357 1.39373 1.14937 0.900344 0.647056 0.390254 0.130577 3.63594 3.4308 3.22561 3.00881
-          2.7866 2.56102 2.33289 2.10247 1.86956 1.63358 1.39373 1.14937 0.90034 0.647051 0.390249 0.130574
-          3.63551 3.43089 3.22579 3.00904 2.78679 2.56116 2.33298 2.10253 1.86959 1.63359 1.39374 1.14937
-          0.900336 0.647044 0.390242 0.13057 3.63487 3.43092 3.22586 3.00923 2.78702 2.56135 2.33313 2.10263
-          1.86966 1.63364 1.39376 1.14938 0.900336 0.647039 0.390235 0.130563 3.63412 3.43088 3.22585 3.00932
-          2.78719 2.56157 2.33333 2.1028 1.86979 1.63373 1.39382 1.14941 0.900351 0.647044 0.390234 0.130554
-          3.63335 3.43086 3.22578 3.00931 2.78727 2.56172 2.33354 2.10303 1.86999 1.63389 1.39393 1.14949
-          0.900401 0.647074 0.390252 0.130543 3.63273 3.43087 3.22575 3.00928 2.78725 2.56175 2.33363 2.10318
-          1.87018 1.63409 1.39412 1.14963 0.900511 0.647152 0.390306 0.130535 3.63273 3.43087 3.22575 3.00928
-          2.78725 2.56175 2.33363 2.10318 1.87018 1.63409 1.39412 1.14963 0.900511 0.647152 0.390306 0.130535
-          3.63335 3.43086 3.22578 3.00931 2.78727 2.56172 2.33354 2.10303 1.86999 1.63389 1.39393 1.14949
-          0.900401 0.647074 0.390252 0.130543 3.63412 3.43088 3.22585 3.00932 2.78719 2.56157 2.33333 2.1028
-          1.86979 1.63373 1.39382 1.14941 0.900351 0.647044 0.390234 0.130554 3.63487 3.43092 3.22586 3.00923
-          2.78702 2.56135 2.33313 2.10263 1.86966 1.63364 1.39376 1.14938 0.900336 0.647039 0.390235 0.130563
-          3.63551 3.43089 3.22579 3.00904 2.78679 2.56116 2.33298 2.10253 1.86959 1.63359 1.39374 1.14937
-          0.900336 0.647044 0.390242 0.13057 3.63594 3.4308 3.22561 3.00881 2.7866 2.56102 2.33289 2.10247
-          1.86956 1.63358 1.39373 1.14937 0.90034 0.647051 0.390249 0.130574 3.63616 3.43063 3.22537 3.00861
-          2.78647 2.56094 2.33284 2.10244 1.86954 1.63357 1.39373 1.14937 0.900344 0.647056 0.390254 0.130577
-          3.63626 3.43035 3.22516 3.00849 2.7864 2.5609 2.33281 2.10242 1.86953 1.63356 1.39373 1.14937
-          0.900345 0.647057 0.390256 0.130578
+          3.68129 3.48919 3.32115 3.1244 2.91495 2.69903 2.47891 2.25556 2.02906 1.79891 1.5642 1.32415
+          1.07851 0.827727 0.572661 0.31402 3.65246 3.50504 3.36615 3.18091 2.97611 2.76276 2.54456 2.32287
+          2.09791 1.86916 1.6357 1.3967 1.15189 0.90173 0.647091 0.388718 3.65097 3.46424 3.31018 3.13393
+          2.93588 2.72592 2.50933 2.28834 2.06357 1.83471 1.60098 1.36162 1.11645 0.865944 0.611014 0.352461
+          3.65076 3.44884 3.26728 3.08696 2.89345 2.68718 2.47268 2.25272 2.02834 1.79948 1.56551 1.32579
+          1.08022 0.829334 0.574074 0.315312 3.65051 3.44482 3.24527 3.04788 2.85193 2.64861 2.43695 2.21884
+          1.99549 1.76709 1.53319 1.29331 1.04744 0.796226 0.540654 0.281686 3.65011 3.44399 3.23847 3.02629
+          2.81653 2.60927 2.39928 2.18357 1.96204 1.73482 1.50151 1.26177 1.0158 0.764366 0.508554 0.249428
+          3.64954 3.44402 3.23701 3.02 2.79976 2.57984 2.36169 2.144 1.92378 1.69839 1.4663 1.22719
+          0.981523 0.730269 0.474605 0.215629 3.64901 3.44414 3.23687 3.0189 2.79602 2.57056 2.34378 2.11644
+          1.88847 1.65883 1.42579 1.18775 0.943999 0.694769 0.440897 0.183235 3.64901 3.44414 3.23687 3.0189
+          2.79602 2.57056 2.34378 2.11644 1.88847 1.65883 1.42579 1.18775 0.943999 0.694769 0.440897 0.183235
+          3.64954 3.44402 3.23701 3.02 2.79976 2.57984 2.36169 2.144 1.92378 1.69839 1.4663 1.22719
+          0.981523 0.730269 0.474605 0.215629 3.65011 3.44399 3.23847 3.02629 2.81653 2.60927 2.39928 2.18357
+          1.96204 1.73482 1.50151 1.26177 1.0158 0.764366 0.508554 0.249428 3.65051 3.44482 3.24527 3.04788
+          2.85193 2.64861 2.43695 2.21884 1.99549 1.76709 1.53319 1.29331 1.04744 0.796226 0.540654 0.281686
+          3.65076 3.44884 3.26728 3.08696 2.89345 2.68718 2.47268 2.25272 2.02834 1.79948 1.56551 1.32579
+          1.08022 0.829334 0.574074 0.315312 3.65097 3.46424 3.31018 3.13393 2.93588 2.72592 2.50933 2.28834
+          2.06357 1.83471 1.60098 1.36162 1.11645 0.865944 0.611014 0.352461 3.65246 3.50504 3.36615 3.18091
+          2.97611 2.76276 2.54456 2.32287 2.09791 1.86916 1.6357 1.3967 1.15189 0.90173 0.647091 0.388718
+          3.68129 3.48919 3.32115 3.1244 2.91495 2.69903 2.47891 2.25556 2.02906 1.79891 1.5642 1.32415
+          1.07851 0.827727 0.572661 0.31402
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.43723e-05 1.43723e-05 1.43724e-05 1.43724e-05 1.43724e-05 1.43725e-05 1.43725e-05 1.43725e-05 1.43726e-05 1.43726e-05 1.43726e-05 1.43727e-05
diff --git a/test/references/test_flatplate_kepsilon2cni.vtu b/test/references/test_flatplate_kepsilon2cni.vtu
index 26b353e6ae..c7de836616 100644
--- a/test/references/test_flatplate_kepsilon2cni.vtu
+++ b/test/references/test_flatplate_kepsilon2cni.vtu
@@ -80,23 +80,23 @@
           0.99875 0.999139 0.999624 1 1 1 1 1
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          0.00153595 0.00100403 0.000823987 0.000678554 0.000550298 0.000434428 0.000328971 0.000230522 0.00013714 4.68887e-05 0.00153583 0.00100406
-          0.000823994 0.000678558 0.0005503 0.000434429 0.000328972 0.000230522 0.00013714 4.68887e-05 0.00153565 0.00100411 0.000824004 0.000678563
-          0.000550304 0.000434431 0.000328974 0.000230523 0.00013714 4.68888e-05 0.00153541 0.00100417 0.000824018 0.00067857 0.000550308 0.000434434
-          0.000328975 0.000230524 0.000137141 4.6889e-05 0.0015351 0.00100426 0.000824039 0.00067858 0.000550315 0.000434438 0.000328978 0.000230526
-          0.000137142 4.68892e-05 0.00153472 0.00100438 0.000824067 0.000678594 0.000550324 0.000434444 0.000328982 0.000230528 0.000137143 4.68895e-05
-          0.00153435 0.00100455 0.000824108 0.000678613 0.000550338 0.000434452 0.000328987 0.000230531 0.000137145 4.68901e-05 0.0015341 0.00100478
-          0.000824167 0.000678642 0.000550356 0.000434463 0.000328994 0.000230536 0.000137148 4.6891e-05 0.00153378 0.00100514 0.000824259 0.000678684
-          0.000550383 0.000434481 0.000329006 0.000230544 0.000137153 4.6893e-05 0.00153287 0.00100577 0.000824424 0.000678755 0.000550422 0.00043451
-          0.000329025 0.000230558 0.000137163 4.68978e-05 0.00153033 0.00100698 0.00082473 0.000678885 0.000550484 0.000434563 0.000329064 0.000230588
-          0.000137186 4.69105e-05 0.00152652 0.00100905 0.000825272 0.000679126 0.000550593 0.000434667 0.000329144 0.000230652 0.00013724 4.69428e-05
-          0.00152109 0.0010121 0.000826159 0.000679536 0.000550781 0.00043486 0.000329295 0.00023078 0.000137351 4.7017e-05 0.00151341 0.00101619
-          0.000827526 0.000680171 0.00055107 0.000435169 0.00032955 0.000231004 0.00013755 4.717e-05 0.00150283 0.00102139 0.000829145 0.00068096
-          0.000551705 0.000435638 0.000329956 0.00023137 0.000137888 4.74707e-05 0.00148903 0.00102804 0.00083124 0.00068194 0.000552355 0.000436097
-          0.000330336 0.000231706 0.000138205 4.79116e-05 0.00147198 0.00103619 0.000833971 0.000683087 0.000553083 0.000436655 0.000330757 0.000232054
-          0.00013854 4.86628e-05 0.00145211 0.00104551 0.000837393 0.000684503 0.000553936 0.000437291 0.000331215 0.000232411 0.000138904 4.99257e-05
-          0.00143129 0.0010549 0.000841212 0.000686112 0.000554897 0.000437997 0.000331714 0.000232792 0.000139323 5.16819e-05 0.00141483 0.00106202
-          0.000844403 0.000687479 0.000555714 0.000438597 0.000332142 0.000233121 0.000139708 5.32684e-05
+          0.00159478 0.00104976 0.000869155 0.000725998 0.000596037 0.000479388 0.000374821 0.000278437 0.000187944 0.000101031 0.0015948 0.00104985
+          0.000869203 0.000726039 0.000596072 0.00047942 0.000374852 0.000278468 0.000187975 0.000101064 0.00159482 0.00104997 0.000869271 0.000726098
+          0.000596123 0.000479465 0.000374895 0.000278511 0.000188019 0.00010111 0.00159485 0.00105014 0.000869366 0.00072618 0.000596193 0.000479529
+          0.000374956 0.000278572 0.000188082 0.000101175 0.00159488 0.00105037 0.000869499 0.000726296 0.000596292 0.000479617 0.000375041 0.000278657
+          0.000188169 0.000101266 0.00159495 0.0010507 0.000869687 0.000726458 0.00059643 0.000479742 0.00037516 0.000278776 0.000188292 0.000101394
+          0.00159516 0.00105115 0.000869951 0.000726686 0.000596625 0.000479917 0.000375328 0.000278944 0.000188464 0.000101573 0.00159555 0.00105176
+          0.000870318 0.000727004 0.000596897 0.000480162 0.000375562 0.000279178 0.000188705 0.000101824 0.00159589 0.0010526 0.000870823 0.000727438
+          0.00059727 0.0004805 0.000375886 0.000279502 0.000189037 0.00010217 0.00159555 0.0010538 0.000871497 0.000728005 0.000597756 0.000480947
+          0.000376316 0.000279932 0.000189478 0.000102626 0.00159406 0.00105557 0.000872362 0.0007287 0.000598344 0.000481499 0.000376849 0.000280465
+          0.000190024 0.000103184 0.00158684 0.00105815 0.000873446 0.000729507 0.000598997 0.000482131 0.000377455 0.00028107 0.000190638 0.000103802
+          0.00157727 0.00106164 0.00087481 0.000730427 0.00059969 0.00048282 0.000378109 0.000281714 0.000191286 0.000104438 0.00156675 0.00106609
+          0.000876577 0.000731498 0.000600415 0.000483564 0.000378807 0.000282396 0.000191964 0.000105091 0.00155494 0.00106956 0.000878234 0.000734346
+          0.000601418 0.0004844 0.000379588 0.000283153 0.000192716 0.000105834 0.00154099 0.00107404 0.000875119 0.00072705 0.000598749 0.000484412
+          0.000381335 0.000286034 0.000196315 0.000109901 0.00152416 0.00108154 0.000874765 0.000721497 0.000590704 0.000474777 0.000370452 0.000274199
+          0.000183798 9.72585e-05 0.00150445 0.00109094 0.000877506 0.000720615 0.000587157 0.000468596 0.000361486 0.000262458 0.000169451 8.14457e-05
+          0.00148371 0.00110043 0.000881341 0.00072193 0.000587235 0.00046752 0.000358976 0.000258275 0.000163458 7.4815e-05 0.00146726 0.00110756
+          0.000884549 0.000723295 0.000587988 0.00046793 0.000359006 0.000257895 0.000162705 7.47477e-05
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.71759e-05 1.71926e-05 1.71956e-05 1.71969e-05 1.71978e-05 1.71984e-05 1.71988e-05 1.7199e-05 1.71991e-05 1.71992e-05 1.71024e-05 1.71577e-05
diff --git a/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu b/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
index 27b0b3c411..a919fe6013 100644
--- a/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
+++ b/test/references/test_flatplate_kepsilon2cni_wallfunction.vtu
@@ -32,11 +32,11 @@
           1 1
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          0.00053604 0.000470416 0.000412942 0.000356245 0.000300092 0.000244574 0.000189766 0.000135716 8.23746e-05 2.86046e-05 0.000530008 0.000470289
-          0.000413008 0.000356306 0.000300144 0.000244621 0.000189815 0.000135772 8.24722e-05 2.96917e-05 0.000525566 0.000470018 0.000413027 0.00035635
-          0.000300193 0.000244673 0.000189869 0.000135835 8.26262e-05 3.09723e-05 0.000522663 0.0004698 0.000413057 0.000356402 0.000300248 0.00024473
-          0.000189926 0.000135896 8.27676e-05 3.19309e-05 0.000521235 0.000469677 0.000413073 0.000356432 0.000300281 0.000244764 0.00018996 0.000135934
-          8.28511e-05 3.24247e-05
+          0.000596665 0.000529826 0.000470942 0.000412733 0.000355025 0.000297952 0.000241617 0.000186087 0.000131325 7.62009e-05 0.000582596 0.000516241
+          0.00045397 0.00039349 0.000334466 0.000276793 0.0002204 0.000165229 0.000111174 5.78648e-05 0.000577914 0.000515347 0.000452855 0.00039179
+          0.000332085 0.000273669 0.000216487 0.000160493 0.000105668 5.26627e-05 0.000575024 0.000515127 0.000452855 0.000391765 0.000331992 0.000273475
+          0.000216159 0.000160002 0.000105059 5.26604e-05 0.000573598 0.000515007 0.000452876 0.000391799 0.000332027 0.000273505 0.000216177 0.000160005
+          0.00010508 5.30607e-05
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.44385e-05 1.44992e-05 1.45555e-05 1.46078e-05 1.46565e-05 1.47021e-05 1.47447e-05 1.47846e-05 1.48221e-05 1.48571e-05 1.4378e-05 1.43872e-05
diff --git a/test/references/test_flatplate_komega2cni.vtu b/test/references/test_flatplate_komega2cni.vtu
index 84f136219f..17c93bf34e 100644
--- a/test/references/test_flatplate_komega2cni.vtu
+++ b/test/references/test_flatplate_komega2cni.vtu
@@ -80,23 +80,23 @@
           0.997714 0.998724 0.999082 0.99934 0.999649 1 1 1
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          0.00140474 0.000879006 0.000707804 0.000578441 0.000470591 0.000374938 0.000286234 0.000201767 0.000120188 4.03348e-05 0.00140463 0.00087904
-          0.000707811 0.000578443 0.000470592 0.000374939 0.000286234 0.000201767 0.000120188 4.03348e-05 0.00140448 0.000879085 0.000707819 0.000578446
-          0.000470593 0.00037494 0.000286234 0.000201767 0.000120188 4.03349e-05 0.00140426 0.000879147 0.000707832 0.000578451 0.000470596 0.000374941
-          0.000286235 0.000201768 0.000120188 4.0335e-05 0.00140394 0.000879229 0.000707849 0.000578458 0.000470599 0.000374942 0.000286236 0.000201769
-          0.000120189 4.03352e-05 0.00140347 0.00087934 0.000707874 0.000578467 0.000470603 0.000374945 0.000286238 0.00020177 0.00012019 4.03355e-05
-          0.00140279 0.000879493 0.000707909 0.000578481 0.000470609 0.000374948 0.00028624 0.000201771 0.000120191 4.03359e-05 0.00140179 0.000879714
-          0.00070796 0.0005785 0.000470618 0.000374953 0.000286242 0.000201773 0.000120192 4.03365e-05 0.00140035 0.000880057 0.000708036 0.000578529
-          0.000470631 0.00037496 0.000286247 0.000201776 0.000120194 4.03376e-05 0.0013983 0.000880637 0.000708168 0.000578582 0.000470658 0.000374975
-          0.000286256 0.000201783 0.000120199 4.034e-05 0.00139544 0.00088166 0.000708431 0.000578709 0.000470732 0.00037502 0.000286285 0.000201803
-          0.000120214 4.0348e-05 0.00139148 0.000883392 0.000708945 0.000578995 0.000470922 0.000375148 0.000286373 0.000201864 0.00012026 4.03736e-05
-          0.00138595 0.000886082 0.000709783 0.000579475 0.00047127 0.000375407 0.000286564 0.000202006 0.00012037 4.04408e-05 0.00137824 0.000889947
-          0.000710997 0.000580108 0.000471728 0.000375777 0.000286867 0.000202253 0.000120577 4.05903e-05 0.0013677 0.000895157 0.000712677 0.00058089
-          0.000472223 0.000376164 0.000287206 0.000202563 0.000120871 4.08869e-05 0.00135406 0.000901833 0.00071492 0.00058187 0.000472761 0.000376519
-          0.000287491 0.000202833 0.000121166 4.14348e-05 0.00133729 0.000910043 0.000717859 0.000583115 0.000473389 0.000376874 0.000287728 0.000203032
-          0.000121416 4.239e-05 0.00131781 0.00091946 0.000721549 0.000584697 0.000474176 0.000377296 0.000287975 0.000203213 0.000121677 4.38069e-05
-          0.00129742 0.000928966 0.000725657 0.000586497 0.000475093 0.000377804 0.000288283 0.000203438 0.000122007 4.54844e-05 0.00128131 0.000936192
-          0.000729076 0.000588014 0.000475873 0.000378247 0.000288562 0.000203652 0.000122325 4.68857e-05
+          0.00140474 0.000879006 0.000707804 0.000578441 0.000470591 0.000374938 0.000286234 0.000201767 0.000120188 4.03348e-05 0.00140464 0.00087904
+          0.000707811 0.000578443 0.000470592 0.000374939 0.000286234 0.000201767 0.000120188 4.03349e-05 0.00140452 0.000879085 0.00070782 0.000578447
+          0.000470594 0.00037494 0.000286235 0.000201768 0.000120189 4.03355e-05 0.00140438 0.000879148 0.000707834 0.000578454 0.000470598 0.000374944
+          0.000286239 0.000201772 0.000120192 4.0339e-05 0.00140425 0.000879239 0.000707859 0.00057847 0.000470614 0.00037496 0.000286255 0.000201788
+          0.000120209 4.03559e-05 0.0014042 0.000879382 0.000707915 0.000578521 0.000470669 0.00037502 0.00028632 0.000201855 0.000120277 4.04255e-05
+          0.00140435 0.000879653 0.000708069 0.00057869 0.000470867 0.000375243 0.00028656 0.000202105 0.000120533 4.06867e-05 0.00140493 0.000880277
+          0.00070854 0.000579257 0.000471553 0.000376024 0.000287402 0.000202985 0.000121432 4.16074e-05 0.00140634 0.000881875 0.000710011 0.000581116
+          0.000473825 0.000378616 0.000290205 0.000205911 0.000124424 4.46726e-05 0.00140901 0.000885783 0.000714218 0.000586622 0.000480609 0.000386369
+          0.000298585 0.000214656 0.000133369 5.38335e-05 0.00141298 0.000893281 0.000723146 0.000598759 0.000495747 0.000403696 0.000317308 0.000234198
+          0.00015337 7.43116e-05 0.00141716 0.000902954 0.000733566 0.00061298 0.00051408 0.000425209 0.000340934 0.000259141 0.000179119 0.000100783
+          0.00141975 0.000912137 0.000738939 0.000617966 0.000521161 0.000435375 0.000353949 0.000274352 0.000195907 0.000118747 0.00141905 0.000921206
+          0.000740743 0.000614203 0.000514224 0.000428365 0.000349017 0.000272227 0.000196444 0.000121289 0.00141389 0.000931443 0.00074388 0.000610992
+          0.000504682 0.000413827 0.000331948 0.000255079 0.000180848 0.000107506 0.00140379 0.00094261 0.000749361 0.00061226 0.000501174 0.000404905
+          0.000317668 0.000236371 0.000159237 8.44216e-05 0.00138888 0.00095384 0.000755524 0.000616005 0.000502673 0.000403614 0.000312928 0.000227667
+          0.00014641 6.8352e-05 0.00137005 0.000964596 0.000761025 0.000619581 0.000505296 0.000405325 0.000313488 0.00022672 0.000143645 6.46147e-05
+          0.00134982 0.000974446 0.000765681 0.000622111 0.000507075 0.000406749 0.00031466 0.000227626 0.00014432 6.61873e-05 0.00133372 0.000981707
+          0.000769171 0.000623742 0.000508016 0.000407398 0.000315182 0.000228109 0.000144914 6.78487e-05
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.71798e-05 1.71959e-05 1.71985e-05 1.71995e-05 1.71999e-05 1.72001e-05 1.72003e-05 1.72004e-05 1.72005e-05 1.72005e-05 1.71158e-05 1.71689e-05
diff --git a/test/references/test_flatplate_lowrekepsilon2cni.vtu b/test/references/test_flatplate_lowrekepsilon2cni.vtu
index f37ba72a96..27504a7361 100644
--- a/test/references/test_flatplate_lowrekepsilon2cni.vtu
+++ b/test/references/test_flatplate_lowrekepsilon2cni.vtu
@@ -80,23 +80,23 @@
           0.997699 0.998529 0.999005 0.99918 0.999373 0.999587 0.999791 1
         </DataArray>
         <DataArray type="Float32" Name="p_rel" NumberOfComponents="1" format="ascii">
-          0.00143114 0.00090716 0.000721685 0.000580979 0.000463752 0.000361613 0.000270003 0.000186157 0.000108239 3.49194e-05 0.00143102 0.000907196
-          0.000721693 0.000580983 0.000463754 0.000361615 0.000270003 0.000186158 0.000108239 3.49196e-05 0.00143086 0.000907245 0.000721703 0.000580988
-          0.000463757 0.000361616 0.000270005 0.000186158 0.00010824 3.49199e-05 0.00143063 0.000907312 0.000721719 0.000580994 0.000463761 0.000361619
-          0.000270006 0.00018616 0.000108241 3.49202e-05 0.00143031 0.000907403 0.00072174 0.000581004 0.000463766 0.000361623 0.000270009 0.000186161
-          0.000108242 3.49207e-05 0.00142986 0.000907525 0.000721771 0.000581018 0.000463774 0.000361628 0.000270012 0.000186164 0.000108243 3.49214e-05
-          0.00142925 0.000907692 0.000721814 0.000581037 0.000463785 0.000361635 0.000270017 0.000186167 0.000108246 3.49223e-05 0.00142846 0.000907924
-          0.000721875 0.000581065 0.000463801 0.000361645 0.000270024 0.000186172 0.000108249 3.49237e-05 0.00142744 0.000908251 0.000721966 0.000581105
-          0.000463823 0.000361659 0.000270034 0.000186179 0.000108253 3.49258e-05 0.00142596 0.000908745 0.000722104 0.000581165 0.000463858 0.000361681
-          0.000270049 0.00018619 0.000108261 3.49296e-05 0.0014236 0.000909554 0.000722337 0.000581268 0.000463917 0.00036172 0.000270077 0.00018621
-          0.000108277 3.49386e-05 0.00141985 0.000910898 0.00072275 0.000581457 0.00046403 0.000361799 0.000270136 0.000186257 0.000108316 3.49633e-05
-          0.00141424 0.00091305 0.000723438 0.000581783 0.000464236 0.00036195 0.000270257 0.000186358 0.000108404 3.50284e-05 0.0014063 0.000916325
-          0.000724504 0.00058228 0.000464547 0.000362185 0.000270452 0.00018653 0.000108564 3.51809e-05 0.00139562 0.00092101 0.00072608 0.000583
-          0.00046497 0.000362485 0.000270697 0.00018675 0.000108785 3.55025e-05 0.00138196 0.000927283 0.000728293 0.000584007 0.000465546 0.00036287
-          0.000270984 0.00018699 0.000109036 3.61241e-05 0.00136513 0.000935119 0.000731247 0.000585352 0.00046631 0.000363371 0.000271345 0.000187276
-          0.000109353 3.71821e-05 0.00134551 0.000944102 0.000734931 0.000587042 0.000467269 0.000363998 0.000271792 0.000187629 0.000109769 3.86489e-05
-          0.00132495 0.000953151 0.000739006 0.000588927 0.00046834 0.000364701 0.000272297 0.000188032 0.000110264 4.03017e-05 0.00130868 0.000960024
-          0.00074239 0.000590506 0.000469235 0.000365288 0.00027272 0.000188372 0.000110694 4.16435e-05
+          0.00143115 0.000907161 0.000721686 0.00058098 0.000463752 0.000361614 0.000270003 0.000186158 0.00010824 3.49203e-05 0.00143116 0.000907217
+          0.000721705 0.000580992 0.000463763 0.000361624 0.000270013 0.000186167 0.00010825 3.49312e-05 0.00143134 0.000907329 0.00072175 0.000581025
+          0.000463792 0.000361652 0.000270041 0.000186197 0.000108281 3.49652e-05 0.00143184 0.000907551 0.000721853 0.000581103 0.000463861 0.00036172
+          0.000270111 0.00018627 0.000108359 3.50497e-05 0.00143296 0.000907995 0.000722074 0.000581272 0.000464016 0.000361872 0.000270268 0.000186435
+          0.000108535 3.52409e-05 0.00143522 0.000908886 0.000722538 0.000581633 0.000464347 0.0003622 0.000270605 0.000186791 0.000108915 3.56541e-05
+          0.0014395 0.000910683 0.000723509 0.000582395 0.000465047 0.000362894 0.000271321 0.000187545 0.00010972 3.65296e-05 0.00144685 0.000914271
+          0.00072555 0.000584009 0.000466532 0.000364365 0.000272836 0.00018914 0.000111421 3.83774e-05 0.00145711 0.000921098 0.000729826 0.000587453
+          0.00046971 0.00036751 0.000276069 0.000192535 0.000115034 4.22953e-05 0.00146695 0.00093236 0.000738326 0.000594667 0.000476463 0.000374212
+          0.00028295 0.000199741 0.000122674 5.0544e-05 0.0014718 0.000945916 0.000752134 0.00060796 0.000489521 0.000387408 0.000296566 0.000213979
+          0.000137693 6.66329e-05 0.0014705 0.000955518 0.000765545 0.000624148 0.000507473 0.000406688 0.000317047 0.000235645 0.000160568 9.09417e-05
+          0.0014648 0.000958893 0.000769735 0.000631248 0.000517861 0.000420035 0.000332889 0.000253568 0.000180233 0.000112093 0.00145645 0.000960301
+          0.000766887 0.000626435 0.000512656 0.000415605 0.000329944 0.000252419 0.000180896 0.000114086 0.00144595 0.000963841 0.000764611 0.000619874
+          0.000502277 0.000401912 0.000313651 0.000234401 0.000162045 9.47715e-05 0.00143292 0.00097043 0.000765782 0.000617613 0.00049674 0.00039289
+          0.000300911 0.000217805 0.000141646 7.10536e-05 0.00141675 0.000979166 0.000769398 0.000618893 0.000496272 0.000390602 0.000296561 0.000211103
+          0.000132338 5.96853e-05 0.00139752 0.000988841 0.000773917 0.000621393 0.000497838 0.000391455 0.000296678 0.000210388 0.000130774 5.81909e-05
+          0.00137715 0.000998219 0.000778445 0.000623822 0.000499502 0.000392753 0.000297728 0.000211233 0.000131553 5.99404e-05 0.00136094 0.00100521
+          0.000781994 0.000625612 0.000500649 0.000393625 0.00029846 0.000211898 0.000132312 6.16052e-05
         </DataArray>
         <DataArray type="Float32" Name="nu" NumberOfComponents="1" format="ascii">
           1.71781e-05 1.71954e-05 1.71986e-05 1.71999e-05 1.72007e-05 1.72012e-05 1.72015e-05 1.72017e-05 1.72019e-05 1.7202e-05 1.71101e-05 1.71671e-05
-- 
GitLab