From d900d8c794bb9d50f528bc8e72ceb594fbc292c8 Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Thu, 17 Jan 2019 13:37:26 +0100
Subject: Add named types to Grid2D constructor

This allows for readable differentiation between constructor overloads.
---
 apps/adrian/cylinder2d/cylinder2d.cpp     | 12 ++++++------
 apps/adrian/poiseuille2d/poiseuille2d.cpp |  6 +++---
 2 files changed, 9 insertions(+), 9 deletions(-)

(limited to 'apps/adrian')

diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index cb0cfc6..739c02c 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -1,8 +1,8 @@
-/*  Lattice Boltzmann sample, written in C++, using the OpenLB
- *  library
+/*
+ *  Lattice Boltzmann grid refinement sample, written in C++,
+ *  using the OpenLB library
  *
- *  Copyright (C) 2007, 2012 Jonas Latt, Mathias J. Krause
- *  Vojtech Cvrcek, Peter Weisbrod
+ *  Copyright (C) 2019 Adrian Kummerländer
  *  E-mail contact: info@openlb.net
  *  The most recent release of OpenLB can be downloaded at
  *  <http://www.openlb.net/>
@@ -170,7 +170,7 @@ void getResults(const std::string& prefix,
     vtmWriter.createMasterFile();
   }
 
-  if (iT%20==0) {
+  if (iT%100==0) {
     vtmWriter.write(iT);
   }
 
@@ -192,7 +192,7 @@ int main(int argc, char* argv[])
   const Vector<T,2> coarseExtend {lx, ly};
   IndicatorCuboid2D<T> coarseCuboid(coarseExtend, coarseOrigin);
 
-  Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, baseTau, Re);
+  Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, RelaxationTime<T>(baseTau), Re);
   prepareGeometry(coarseGrid);
 
   const Vector<T,2> fineExtend {3.0, 1.5};
diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index d92934b..19dc4d7 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -38,9 +38,9 @@ typedef double T;
 
 const T lx  = 4.0;            // length of the channel
 const T ly  = 1.0;            // height of the channel
-const int N = 30;             // resolution of the model
+const int N = 50;             // resolution of the model
 const T Re  = 100.;           // Reynolds number
-const T baseTau = 0.8;        // Relaxation time of coarsest grid
+const T uMax = 0.02;          // Max lattice speed
 const T maxPhysT = 60.;       // max. simulation time in s, SI unit
 const T physInterval = 0.25;  // interval for the convergence check in s
 const T residuum = 1e-5;      // residuum for the convergence check
@@ -184,7 +184,7 @@ int main(int argc, char* argv[])
   const Vector<T,2> coarseExtend {lx/2, ly};
   IndicatorCuboid2D<T> coarseCuboid(coarseExtend, coarseOrigin);
 
-  Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, baseTau, Re);
+  Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, LatticeVelocity<T>(uMax), Re);
   prepareGeometry(coarseGrid);
 
   const T coarseDeltaX = coarseGrid.getConverter().getPhysDeltaX();
-- 
cgit v1.2.3