diff --git a/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp b/src/cpu/VirtualFluidsCore/LBM/RheologyK17LBMKernel.cpp
index 43713cbc676e2aacb6775bd8912ff0c1b59cfb77..c7bd879cbf87f97fe41a70d126e3a70c7905790c 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);