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
---
.../advectionDiffusionBoundaryPostProcessor3D.h | 159 +++++++++++++++++++++
1 file changed, 159 insertions(+)
create mode 100644 src/boundary/advectionDiffusionBoundaryPostProcessor3D.h
(limited to 'src/boundary/advectionDiffusionBoundaryPostProcessor3D.h')
diff --git a/src/boundary/advectionDiffusionBoundaryPostProcessor3D.h b/src/boundary/advectionDiffusionBoundaryPostProcessor3D.h
new file mode 100644
index 0000000..9961040
--- /dev/null
+++ b/src/boundary/advectionDiffusionBoundaryPostProcessor3D.h
@@ -0,0 +1,159 @@
+/* This file is part of the OpenLB library
+ *
+ * Copyright (C) 2006, 2016 Robin Trunk
+ * E-mail contact: info@openlb.net
+ * The most recent release of OpenLB can be downloaded at
+ *
+ *
+ * Generic version of the collision, which modifies the particle
+ * distribution functions, by Orestis Malaspinas.
+ *
+ * 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.
+*/
+
+#include "core/postProcessing.h"
+#include "core/blockLattice3D.h"
+
+namespace olb {
+
+/**
+* This class interpolates missing f_i from values
+* near the boundary to get a more stable outflow
+* condition for the density. It is assumed that the
+* next two cells in direction of this f_i have
+* viable values.
+*/
+template
+class ConvectionBoundaryProcessor3D : public LocalPostProcessor3D {
+public:
+ ConvectionBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
+ int z1_, int discreteNormalX_,
+ int discreteNormalY_, int discreteNormalZ_);
+ int extent() const override
+ {
+ return 0;
+ }
+ int extent(int whichDirection) const override
+ {
+ return 0;
+ }
+ void process(BlockLattice3D& blockLattice) override;
+ void processSubDomain ( BlockLattice3D& blockLattice,
+ int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
+private:
+ int interpolationPop[DESCRIPTOR::q];
+ int x0, x1, y0, y1, z0, z1;
+};
+
+template
+class ConvectionBoundaryProcessorGenerator3D : public PostProcessorGenerator3D {
+public:
+ ConvectionBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_,
+ int z0_, int z1_, int discreteNormalX_,
+ int discreteNormalY_, int discreteNormalZ_);
+ PostProcessor3D* generate() const override;
+ PostProcessorGenerator3D* clone() const override;
+private:
+ int discreteNormalX;
+ int discreteNormalY;
+ int discreteNormalZ;
+};
+
+/**
+* This class copies missing values in the
+* external field from the neighbour in normal direction.
+* Therefore it is assumed this neighbour is a fluid cell.
+*/
+template
+class ExtFieldBoundaryProcessor3D : public LocalPostProcessor3D {
+public:
+ ExtFieldBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
+ int z1_, int discreteNormalX_, int discreteNormalY_,
+ int discreteNormalZ_, int offset_);
+ int extent() const override
+ {
+ return 0;
+ }
+ int extent(int whichDirection) const override
+ {
+ return 0;
+ }
+ void process(BlockLattice3D& blockLattice) override;
+ void processSubDomain ( BlockLattice3D& blockLattice,
+ int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
+private:
+ int x0, x1, y0, y1, z0, z1;
+ int discreteNormalX, discreteNormalY, discreteNormalZ;
+ int offset;
+ bool par;
+};
+
+template
+class ExtFieldBoundaryProcessorGenerator3D : public PostProcessorGenerator3D {
+public:
+ ExtFieldBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_, int z0_,
+ int z1_, int discreteNormalX_, int discreteNormalY_,
+ int discreteNormalZ_, int offset_);
+ PostProcessor3D* generate() const override;
+ PostProcessorGenerator3D* clone() const override;
+private:
+ int discreteNormalX, discreteNormalY, discreteNormalZ;
+ int offset;
+};
+
+/**
+* This class resets some values of the distribution
+* on the boundary that can have arbitrary values
+* to be zero and thus ensures a correct computation
+* of the density that is about to leave the domain.
+*/
+template
+class ZeroDistributionBoundaryProcessor3D : public LocalPostProcessor3D {
+public:
+ ZeroDistributionBoundaryProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_,
+ int z1_, int discreteNormalX_, int discreteNormalY_,
+ int discreteNormalZ_);
+ int extent() const override
+ {
+ return 0;
+ }
+ int extent(int whichDirection) const override
+ {
+ return 0;
+ }
+ void process(BlockLattice3D& blockLattice) override;
+ void processSubDomain ( BlockLattice3D& blockLattice,
+ int x0_, int x1_, int y0_, int y1_ , int z0_, int z1_) override;
+private:
+ int resetPop[DESCRIPTOR::q];
+ int x0, x1, y0, y1, z0, z1;
+};
+
+template
+class ZeroDistributionBoundaryProcessorGenerator3D : public PostProcessorGenerator3D {
+public:
+ ZeroDistributionBoundaryProcessorGenerator3D(int x0_, int x1_, int y0_, int y1_,
+ int z0_, int z1_, int discreteNormalX_,
+ int discreteNormalY_, int discreteNormalZ_);
+ PostProcessor3D* generate() const override;
+ PostProcessorGenerator3D* clone() const override;
+private:
+ int discreteNormalX;
+ int discreteNormalY;
+ int discreteNormalZ;
+};
+}
+
--
cgit v1.2.3