From 94d3e79a8617f88dc0219cfdeedfa3147833719d Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Mon, 24 Jun 2019 14:43:36 +0200
Subject: Initialize at openlb-1-3

---
 src/boundary/rtlbmBoundaryDynamics.h | 177 +++++++++++++++++++++++++++++++++++
 1 file changed, 177 insertions(+)
 create mode 100644 src/boundary/rtlbmBoundaryDynamics.h

(limited to 'src/boundary/rtlbmBoundaryDynamics.h')

diff --git a/src/boundary/rtlbmBoundaryDynamics.h b/src/boundary/rtlbmBoundaryDynamics.h
new file mode 100644
index 0000000..576551b
--- /dev/null
+++ b/src/boundary/rtlbmBoundaryDynamics.h
@@ -0,0 +1,177 @@
+/*  This file is part of the OpenLB library
+ *
+ *  Copyright (C) 2017 Albert Mink
+ *  E-mail contact: info@openlb.net
+ *  The most recent release of OpenLB can be downloaded at
+ *  <http://www.openlb.net/>
+ *
+ *  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 2
+ *  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, write to the Free
+ *  Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
+ *  Boston, MA  02110-1301, USA.
+*/
+
+#ifndef RTLBM_BOUNDARY_DYNAMICS_H
+#define RTLBM_BOUNDARY_DYNAMICS_H
+
+
+#include "dynamics/dynamics.h"
+
+namespace olb {
+
+/** Defines incoming (axis parallel) directions on flat walls.
+ *
+ * \param directions   takes values 0, 1, 2; which corresponds to x, y, z plane normal
+ * \param orientation  takes values -1, 1; which is the sign of plane normal
+ *
+ *  Note:
+ *  Dynamics computes the arriving density on walls and redistributes it equally to
+ *  all inward pointing directions.
+ */
+template<typename T, typename DESCRIPTOR, int direction, int orientation>
+class RtlbmDiffuseBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for flat boundary
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+template<typename T, typename DESCRIPTOR, int plane, int normal1, int normal2>
+class RtlbmDiffuseEdgeBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseEdgeBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for flat boundary
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+template<typename T, typename DESCRIPTOR, int xNormal, int yNormal, int zNormal>
+class RtlbmDiffuseCornerBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseCornerBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for corner
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+
+/** Defines incoming directions on flat walls.
+ *  Emposed dirichlet density is distributed equivalently to all incoming/unkown directions.
+ *
+ * \param directions   takes values 0, 1, 2; which corresponds to x, y, z plane normal
+ * \param orientation  takes values -1, 1; which is the sign of plane normal
+ *
+ *  Note:
+ *  Dynamics ensure that cell has constant 0th moment.
+ *  To set the constant 0th moment defineRho() must be called.
+ */
+template<typename T, typename DESCRIPTOR, int direction, int orientation>
+class RtlbmDiffuseConstBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseConstBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for flat boundary
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+
+
+/** Defines incoming directions on edge boundaries.
+*   Emposed dirichlet density is distributed equivalently to all incoming/unkown directions.
+*
+* \param plane    takes values 0, 1, 2; which denotes the alignment of edge either along x, y or z axis
+* \param normal1  takes values -1 or 1; normal orientation of a surface side 'right'
+* \param normal2  takes values -1 or 1; normal orientation of a surface side 'left'
+*/
+template<typename T, typename DESCRIPTOR, int plane, int normal1, int normal2>
+class RtlbmDiffuseConstEdgeBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseConstEdgeBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for edges
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+
+/** Defines incoming directions on corner boundaries.
+*   Emposed dirichlet density is distributed equivalently to all incoming/unkown directions.
+*
+* \param xNormal  takes value -1 or 1
+* \param yNormal  takes value -1 or 1
+* \param zNormal  takes value -1 or 1
+*/
+template<typename T, typename DESCRIPTOR, int xNormal, int yNormal, int zNormal>
+class RtlbmDiffuseConstCornerBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDiffuseConstCornerBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for corner
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+
+
+
+template<typename T, typename DESCRIPTOR, int direction, int orientation>
+class RtlbmDirectedBoundaryDynamics : public BasicDynamics<T,DESCRIPTOR> {
+public:
+  /// Constructor
+  RtlbmDirectedBoundaryDynamics(T omega_, Momenta<T,DESCRIPTOR>& momenta_);
+  /// Compute equilibrium distribution function
+  T computeEquilibrium(int iPop, T rho, const T u[DESCRIPTOR::d], T uSqr) const override;
+  /// Collision step for directed boundary walls
+  void collide(Cell<T,DESCRIPTOR>& cell, LatticeStatistics<T>& statistics) override;
+  /// place holder
+  T getOmega() const override;
+  /// place holder
+  void setOmega(T omega_) override;
+};
+}  // namespace olb
+
+#endif
-- 
cgit v1.2.3