From 6980149cdce3282b63fcc00cd1f888e719219dfc Mon Sep 17 00:00:00 2001
From: kutscher <kutscher@irmb.tu-bs.de>
Date: Fri, 25 Sep 2020 14:32:37 +0200
Subject: [PATCH] fix calculation of dxux, dyuy, dzuz, Dxy, Dxz, Dyz in
RheologyK17LBMKernel
---
src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp b/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp
index 43713cbc6..c7bd879cb 100644
--- a/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp
+++ b/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp
@@ -621,6 +621,15 @@ void RheologyK17LBMKernel::calculate(int step)
omega = Thixotropy::getHerschelBulkleyCollFactor(omega, shearRate, drho);
//omega = Thixotropy::getBinghamCollFactor(omega, shearRate, drho);
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+
+ dxux = c1o2 * (-omega) * (mxxMyy + mxxMzz);// +c1o2 * OxxPyyPzz * (mfaaa - mxxPyyPzz);
+ dyuy = dxux + omega * c3o2 * mxxMyy;
+ dzuz = dxux + omega * c3o2 * mxxMzz;
+
+ Dxy = -three * omega * mfbba;
+ Dxz = -three * omega * mfbab;
+ Dyz = -three * omega * mfabb;
+
mxxMyy += omega * (-mxxMyy) - three * (one + c1o2 * (-omega)) * (vx2 * dxux - vy2 * dyuy);
mxxMzz += omega * (-mxxMzz) - three * (one + c1o2 * (-omega)) * (vx2 * dxux - vz2 * dzuz);
--
GitLab