From 52c2bf1526abbbf17f5de44310137f77a60f5e35 Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Thu, 11 Oct 2018 21:03:59 +0200
Subject: Adapt lid driven cavity example to Re 1000

I suspect something is not quite right with my Zou He velocity boundary impl.
---
 cavity2d.cc | 17 +++++++++--------
 1 file changed, 9 insertions(+), 8 deletions(-)

diff --git a/cavity2d.cc b/cavity2d.cc
index 0006252..ca4b030 100644
--- a/cavity2d.cc
+++ b/cavity2d.cc
@@ -6,13 +6,13 @@
 #include "lbm.h"
 #include "boundary_conditions.h"
 
-constexpr std::size_t dimX = 64;
+constexpr std::size_t dimX = 100;
 constexpr std::size_t dimY = dimX;
 
-constexpr double uLid     = 0.1;
-constexpr double reynolds = 100;
+constexpr double uLid     = 0.4;
+constexpr double reynolds = 1000;
 
-constexpr double tau   = 3 * uLid * (dimX-1) / reynolds + 0.5;
+constexpr double tau   = 3. * uLid * (dimX-1) / reynolds + 0.5;
 constexpr double omega = 1. / tau;
 
 DataCellBuffer pop(dimX, dimY);
@@ -71,11 +71,12 @@ void computeLbmStep() {
 int main() {
 	init();
 
-	std::cout << "Re:   " << reynolds << std::endl;
-	std::cout << "uLid: " << uLid     << std::endl;
-	std::cout << "tau:  " << tau      << std::endl;
+	std::cout << "Re:    " << reynolds << std::endl;
+	std::cout << "uLid:  " << uLid     << std::endl;
+	std::cout << "tau:   " << tau      << std::endl;
+	std::cout << "omega: " << omega    << std::endl;
 
-	for ( std::size_t t = 0; t <= 6000; ++t ) {
+	for ( std::size_t t = 0; t <= 10000; ++t ) {
 		computeLbmStep();
 
 		if ( t % 100 == 0 ) {
-- 
cgit v1.2.3