diff options
| author | Adrian Kummerlaender | 2018-10-11 21:03:59 +0200 |
|---|---|---|
| committer | Adrian Kummerlaender | 2018-10-11 21:04:29 +0200 |
| commit | 52c2bf1526abbbf17f5de44310137f77a60f5e35 (patch) | |
| tree | 605ef5a744e6ea0148b52c85f50e017beaae8bbe | |
| parent | a2476312c64a565736beaf791df1a464e7e2ef68 (diff) | |
| download | boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar.gz boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar.bz2 boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar.lz boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar.xz boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.tar.zst boltzbub-52c2bf1526abbbf17f5de44310137f77a60f5e35.zip | |
Adapt lid driven cavity example to Re 1000
I suspect something is not quite right with my Zou He velocity boundary impl.
| -rw-r--r-- | cavity2d.cc | 17 |
1 files 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 ) { |