From 16c83ab939ba0f0a01651162efa3099e90a53fcf Mon Sep 17 00:00:00 2001
From: Anna Wellmann <a.wellmann@tu-bs.de>
Date: Thu, 19 Jan 2023 10:50:37 +0000
Subject: [PATCH] In DrivenCavity: calculate dt from velocityLB
---
apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
index 4ef9786cd..5a9e94c82 100644
--- a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
+++ b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
@@ -85,7 +85,7 @@ int main()
const real L = 1.0;
const real Re = 1000.0;
const real velocity = 1.0;
- const real dt = (real)0.5e-3;
+ const real velocityLB = 0.05; // LB units
const uint nx = 64;
const uint timeStepOut = 1000;
@@ -128,7 +128,7 @@ int main()
// compute parameters in lattice units
//////////////////////////////////////////////////////////////////////////
- const real velocityLB = velocity * dt / dx; // LB units
+ const real dt = velocityLB / velocity * dx;
const real vxLB = velocityLB / sqrt(2.0); // LB units
const real vyLB = velocityLB / sqrt(2.0); // LB units
--
GitLab