From 25b1a74030dc5fbcc4789739c3947cc69ca82a56 Mon Sep 17 00:00:00 2001
From: Ned Coltman <edward.coltman@iws.uni-stuttgart.de>
Date: Thu, 16 Sep 2021 14:24:05 +0200
Subject: [PATCH] [rans][boundarytypes] Introduce setWall boundary designation
for the rans models
---
dumux/freeflow/rans/boundarytypes.hh | 106 +++++++++++++++++++++++++++
1 file changed, 106 insertions(+)
create mode 100644 dumux/freeflow/rans/boundarytypes.hh
diff --git a/dumux/freeflow/rans/boundarytypes.hh b/dumux/freeflow/rans/boundarytypes.hh
new file mode 100644
index 0000000000..5d71503060
--- /dev/null
+++ b/dumux/freeflow/rans/boundarytypes.hh
@@ -0,0 +1,106 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 3 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file
+ * \ingroup RANSModel
+ * \copydoc Dumux::RANSBoundaryTypes
+ */
+#ifndef FREEFLOW_RANS_BOUNDARY_TYPES_HH
+#define FREEFLOW_RANS_BOUNDARY_TYPES_HH
+
+#include <dumux/common/typetraits/typetraits.hh>
+#include <dumux/freeflow/navierstokes/boundarytypes.hh>
+
+namespace Dumux {
+
+/*!
+ * \ingroup RANSModel
+ * \brief Class to specify the type of a boundary condition for the RANS extension to the Navier-Stokes model.
+ */
+template <class ModelTraits, int numEq>
+class RANSBoundaryTypes : public NavierStokesBoundaryTypes<numEq>
+{
+ using ParentType = NavierStokesBoundaryTypes<numEq>;
+ static constexpr auto dimWorld = ModelTraits::dim();
+ using Indices = typename ModelTraits::Indices;
+ static_assert(dimWorld > 1, "Wall conditions cannot be set for 1D domains.");
+public:
+ RANSBoundaryTypes()
+ {
+ for (int eqIdx=0; eqIdx < numEq; ++eqIdx)
+ resetEq(eqIdx);
+ }
+
+ void setWall()
+ {
+ isWall_ = true;
+ for (int eqIdx=0; eqIdx < numEq; ++eqIdx)
+ {
+ if constexpr (dimWorld == 3)
+ {
+ if ((eqIdx == Indices::velocityXIdx)
+ || (eqIdx == Indices::velocityYIdx)
+ || (eqIdx == Indices::velocityZIdx))
+ BoundaryTypes<numEq>::boundaryInfo_[eqIdx].isDirichlet = true;
+ }
+ else if constexpr (dimWorld == 2)
+ {
+ if ((eqIdx == Indices::velocityXIdx)
+ || (eqIdx == Indices::velocityYIdx))
+ BoundaryTypes<numEq>::boundaryInfo_[eqIdx].isDirichlet = true;
+ }
+ else
+ DUNE_THROW(Dune::NotImplemented, "1D Turbulence models are not supported");
+
+ if constexpr (numTurbulenceEq(ModelTraits::turbulenceModel()) == 1)
+ {
+ if (eqIdx == Indices::viscosityTildeIdx)
+ BoundaryTypes<numEq>::boundaryInfo_[eqIdx].isDirichlet = true;
+ }
+ else if constexpr (numTurbulenceEq(ModelTraits::turbulenceModel()) == 2)
+ {
+ if (eqIdx == Indices::turbulentKineticEnergyEqIdx || eqIdx == Indices::dissipationEqIdx)
+ BoundaryTypes<numEq>::boundaryInfo_[eqIdx].isDirichlet = true;
+ }
+ }
+ }
+
+ /*!
+ * \brief Returns true if some equation is used to specify a
+ * Wall boundary condition.
+ */
+ bool hasWall() const
+ { return isWall_; }
+
+ /*!
+ * \brief Reset the boundary types for one equation.
+ */
+ void resetEq(const int eqIdx)
+ {
+ ParentType::resetEq(eqIdx);
+ isWall_ = false;
+ }
+
+protected:
+ bool isWall_;
+};
+
+} // end namespace Dumux
+
+#endif
--
GitLab