From 89e6d3e48e2491e135c00a208fabab881d438bb5 Mon Sep 17 00:00:00 2001
From: "TESLA03\\Master" <a.wellmann@tu-bs.de>
Date: Tue, 22 Jun 2021 12:03:04 +0200
Subject: [PATCH] Use variables for the extent of the coarse grid
---
apps/gpu/LBM/MusselOyster/MusselOyster.cpp | 38 ++++++++++++----------
1 file changed, 20 insertions(+), 18 deletions(-)
diff --git a/apps/gpu/LBM/MusselOyster/MusselOyster.cpp b/apps/gpu/LBM/MusselOyster/MusselOyster.cpp
index d0e152c00..e208281d9 100644
--- a/apps/gpu/LBM/MusselOyster/MusselOyster.cpp
+++ b/apps/gpu/LBM/MusselOyster/MusselOyster.cpp
@@ -156,25 +156,26 @@ void multipleLevel(const std::string& configPath)
TriangularMesh::make("C:/Users/Master/Documents/MasterAnna/STL/" + bivalveType + "_Level1.stl");
// bounding box mussel:
- // x = -18, 58
- // y = -17, 18
- // z = -5, 13
+ const real bbxm = -18.0;
+ const real bbxp = 58.0;
+ const real bbym = -17.0;
+ const real bbyp = 18.0;
+ const real bbzm = -5.0;
+ const real bbzp = 13.0;
// bounding box oyster:
- // x = 0, 115
- // y = 0, 27
- // z = 0, 63
-
- const real xSpaceM = 40.0;
- const real xSpaceP = 300.0;
- const real ySpaceP = 60.0;
- const real zSpacePM = 30.0;
- if (bivalveType == "MUSSEL")
- gridBuilder->addCoarseGrid(-18.0 - xSpaceM, -16.5, -5.0 - zSpacePM,
- 58.0 + xSpaceP, 18.0 + ySpaceP, 13.0 + zSpacePM, dx);
- else if (bivalveType == "OYSTER")
- gridBuilder->addCoarseGrid(0.0 - xSpaceM, 0.5, 0.0 - zSpacePM,
- 115.0 + xSpaceP, 27.0 + ySpaceP, 63.0 + zSpacePM, dx);
-
+ // const real bbxm = 0.0;
+ // const real bbxp = 115.0;
+ // const real bbym = 0.0;
+ // const real bbyp = 27.0;
+ // const real bbzm = 0.0;
+ // const real bbzp = 63.0;
+
+ const real xGridMin = bbxm - 40.0;
+ const real xGridMax = bbxp + 250.0;
+ const real yGridMin = bbym + 0.5;
+ const real yGridMax = bbyp + 60.0;
+ const real zGridMin = bbzm - 30.0;
+ const real zGridMax = bbzp + 30.0;
gridBuilder->setNumberOfLayers(6, 8);
gridBuilder->addGrid(bivalveRef_1_STL, 1);
@@ -182,6 +183,7 @@ void multipleLevel(const std::string& configPath)
gridBuilder->addGeometry(bivalveSTL);
gridBuilder->setPeriodicBoundaryCondition(false, false, true);
+ gridBuilder->addCoarseGrid(xGridMin, yGridMin, zGridMin, xGridMax, yGridMax, zGridMax, dx);
gridBuilder->buildGrids(LBM, true); // buildGrids() has to be called before setting the BCs!!!!
//////////////////////////////////////////////////////////////////////////
--
GitLab