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/superOffBoundaryCondition2D.h | 163 +++++++++++++++++++++++++++++
1 file changed, 163 insertions(+)
create mode 100644 src/boundary/superOffBoundaryCondition2D.h
(limited to 'src/boundary/superOffBoundaryCondition2D.h')
diff --git a/src/boundary/superOffBoundaryCondition2D.h b/src/boundary/superOffBoundaryCondition2D.h
new file mode 100644
index 0000000..6fb1f56
--- /dev/null
+++ b/src/boundary/superOffBoundaryCondition2D.h
@@ -0,0 +1,163 @@
+/* This file is part of the OpenLB library
+ *
+ * Copyright (C) 2012, 2016 Jonas Kratzke, Mathias J. Krause
+ * E-mail contact: info@openlb.net
+ * The most recent release of OpenLB can be downloaded at
+ *
+ *
+ * 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.
+*/
+
+/** \file
+ * A helper for initialising 2D boundaries -- header file.
+ */
+
+
+#ifndef SUPER_OFF_BOUNDARY_CONDITION_2D_H
+#define SUPER_OFF_BOUNDARY_CONDITION_2D_H
+
+#include
+#include
+
+#include "offBoundaryCondition2D.h"
+#include "geometry/superGeometry2D.h"
+#include "core/superLattice2D.h"
+#include "io/ostreamManager.h"
+#include "functors/analytical/analyticalF.h"
+#include "utilities/functorPtr.h"
+
+/// All OpenLB code is contained in this namespace.
+namespace olb {
+
+/// A helper for initialising 2D boundaries for super lattices.
+/** Here we have methods that initializes for a given global
+ * point or global range the local postprocessors and the
+ * communicator (_commBC in SuperLattice) for boundary conditions.
+ *
+ * This class is not intended to be derived from.
+ */
+template
+class sOffLatticeBoundaryCondition2D {
+public:
+ /// Constructor
+ sOffLatticeBoundaryCondition2D(SuperLattice2D& sLattice, T epsFraction_ = 0.0001);
+ /// Copy construction
+ sOffLatticeBoundaryCondition2D(sOffLatticeBoundaryCondition2D const& rhs);
+ /// Copy assignment
+ sOffLatticeBoundaryCondition2D operator=(sOffLatticeBoundaryCondition2D rhs);
+ /// Destructor
+ ~sOffLatticeBoundaryCondition2D();
+
+ /// Set offDynamics with boundary links and post processors using indicators
+ /**
+ * Add offDynamics with initialisation of boundary links and the corresponding
+ * post processors
+ * Note: Uses information of the second neighbours of the cell (x,y)
+ * Add post processors. Ensure that offDynamics are defined!
+ *
+ * \param boundaryIndicator Indicator describing boundary cells
+ * \param bulkIndicator Indicator describing bulk cells
+ * \param geometryIndicator Indicator describing the geometry to be bounded
+ **/
+ void addVelocityBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator,
+ IndicatorF2D& geometryIndicator);
+ void addVelocityBoundary(SuperGeometry2D& superGeometry, int material,
+ IndicatorF2D& geometryIndicator,
+ std::vector bulkMaterials = std::vector(1,1));
+ void addVelocityBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator);
+ void addVelocityBoundary(SuperGeometry2D& superGeometry, int material,
+ std::vector bulkMaterials = std::vector(1,1));
+
+ void addZeroVelocityBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator,
+ IndicatorF2D& geometryIndicator);
+ void addZeroVelocityBoundary(SuperGeometry2D& superGeometry, int material,
+ IndicatorF2D& geometryIndicator,
+ std::vector bulkMaterials = std::vector(1,1));
+ void addZeroVelocityBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator);
+ void addZeroVelocityBoundary(SuperGeometry2D& superGeometry, int material,
+ std::vector bulkMaterials = std::vector(1,1));
+
+ void addPressureBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator,
+ IndicatorF2D& geometryIndicator);
+ void addPressureBoundary(SuperGeometry2D& superGeometry, int material,
+ IndicatorF2D& geometryIndicator,
+ std::vector bulkMaterials = std::vector(1,1));
+ void addPressureBoundary(FunctorPtr>&& boundaryIndicator,
+ FunctorPtr>&& bulkIndicator);
+ void addPressureBoundary(SuperGeometry2D& superGeometry, int material,
+ std::vector bulkMaterials = std::vector(1,1));
+
+ void defineU(FunctorPtr>&& indicator,
+ FunctorPtr>&& bulkIndicator,
+ AnalyticalF2D& u);
+ void defineU(SuperGeometry2D& superGeometry, int material,
+ AnalyticalF2D& u,
+ std::vector bulkMaterials = std::vector(1,1));
+
+ void defineRho(FunctorPtr>&& indicator,
+ FunctorPtr>&& bulkIndicator,
+ AnalyticalF2D& rho);
+ void defineRho(SuperGeometry2D& superGeometry, int material,
+ AnalyticalF2D& rho,
+ std::vector bulkMaterials = std::vector(1,1));
+
+ /// Adds needed Cells to the Communicator _commBC in SuperLattice
+ void addPoints2CommBC(FunctorPtr>&& indicator);
+ void addPoints2CommBC(SuperGeometry2D& superGeometry, int material);
+
+ SuperLattice2D& getSuperLattice()
+ {
+ return _sLattice;
+ };
+ std::vector* >& getBlockBCs()
+ {
+ return _blockBCs;
+ };
+ int getOverlap()
+ {
+ return _overlap;
+ };
+ void setOverlap(int overlap)
+ {
+ _overlap = overlap;
+ };
+
+ void outputOn();
+ void outputOff();
+
+private:
+ mutable OstreamManager clout;
+ SuperLattice2D& _sLattice;
+ std::vector* > _blockBCs;
+ T _epsFraction;
+ int _overlap;
+ bool _output;
+};
+
+template >
+void createBouzidiBoundaryCondition2D(sOffLatticeBoundaryCondition2D& sBC);
+
+template
+void createBounceBackBoundaryCondition2D(sOffLatticeBoundaryCondition2D& sBC);
+
+} // namespace olb
+
+#endif
--
cgit v1.2.3