diff --git a/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh b/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
index ba1034ac59b52522c8afd881d9f91f90e701d5ae..6bb51f95c925b701589dc830959c1343549ae90e 100644
--- a/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
+++ b/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
@@ -4,12 +4,15 @@
 #include <cuda_runtime.h>
 
 #include <basics/DataTypes.h>
+#include <basics/constants/NumericConstants.h>
 
 #include <lbm/collision/CollisionParameter.h>
 #include <lbm/collision/TurbulentViscosity.h>
 
 #include "LBM/GPUHelperFunctions/KernelUtilities.h"
 
+using namespace vf::basics::constant;
+
 namespace vf::gpu
 {
 
@@ -58,9 +61,9 @@ __global__ void runCollision(CollisionFunctor collision, GPUCollisionParameter c
         para.forceZ = (collisionParameter.forces[2] + collisionParameter.bodyForceZ[k_000]) * c1o2 * collisionParameter.forceFactor;
 
         // Reset body force. To be used when not using round-off correction.
-        collisionParameter.bodyForceX[k_000] = 0.0f;
-        collisionParameter.bodyForceX[k_000] = 0.0f;
-        collisionParameter.bodyForceX[k_000] = 0.0f;
+        collisionParameter.bodyForceX[k_000] = c0o1;
+        collisionParameter.bodyForceY[k_000] = c0o1;
+        collisionParameter.bodyForceZ[k_000] = c0o1;
 
         ////////////////////////////////////////////////////////////////////////////////////
         //!> Round-off correction