From 8731e76a105bb03221c114c95bc4f7873ae7ba5c Mon Sep 17 00:00:00 2001
From: niikonst <niikonst@blogin1.usr.hlrn.de>
Date: Mon, 18 Sep 2023 13:01:51 +0200
Subject: [PATCH] adds wetting control depending on BC
---
.../LBM/MultiphaseScaleDistributionLBMKernel.cpp | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/src/cpu/MultiphaseFlow/LBM/MultiphaseScaleDistributionLBMKernel.cpp b/src/cpu/MultiphaseFlow/LBM/MultiphaseScaleDistributionLBMKernel.cpp
index 64a08aa86..83cb01798 100644
--- a/src/cpu/MultiphaseFlow/LBM/MultiphaseScaleDistributionLBMKernel.cpp
+++ b/src/cpu/MultiphaseFlow/LBM/MultiphaseScaleDistributionLBMKernel.cpp
@@ -5229,8 +5229,16 @@ void MultiphaseScaleDistributionLBMKernel::findNeighbors(CbArray3D<real, Indexer
//if (bcArray->getBC(x1, x2, x3)->hasVelocityBoundaryFlag(D3Q27System::INVDIR[k]))
// phi[k] = (*ph)(x1, x2, x3); // neutral wetting
//else
- phi[k] = 0.0; // unwetting
+ // phi[k] = 0.0; // unwetting
//phi[k] = (*ph)(x1, x2, x3) * 0.7;
+ if (bcArray->getBC(x1, x2, x3)->hasNoSlipBoundaryFlag(D3Q27System::INVDIR[k]))
+ {
+ if(bcArray->getBC(x1, x2, x3)->getNoSlipSecondaryOption(D3Q27System::INVDIR[k]) == 0)
+ phi[k] = (*ph)(x1, x2, x3); // neutral wetting
+ else
+ phi[k] = 0.0; // unwetting
+ }
+
}
}
}
--
GitLab