diff options
author | Adrian Kummerlaender | 2018-10-11 21:12:22 +0200 |
---|---|---|
committer | Adrian Kummerlaender | 2018-10-11 21:12:22 +0200 |
commit | 4e8921427d4dd6717ff4ef262e7c04bbbd21a906 (patch) | |
tree | 25b6ca1f98c2c9f483739a881ce28366980acaa5 /lid_driven_cavity.cc | |
parent | 52c2bf1526abbbf17f5de44310137f77a60f5e35 (diff) | |
download | boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar.gz boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar.bz2 boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar.lz boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar.xz boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.tar.zst boltzbub-4e8921427d4dd6717ff4ef262e7c04bbbd21a906.zip |
Rename examples
Diffstat (limited to 'lid_driven_cavity.cc')
-rw-r--r-- | lid_driven_cavity.cc | 90 |
1 files changed, 90 insertions, 0 deletions
diff --git a/lid_driven_cavity.cc b/lid_driven_cavity.cc new file mode 100644 index 0000000..ca4b030 --- /dev/null +++ b/lid_driven_cavity.cc @@ -0,0 +1,90 @@ +#include <iostream> +#include <vector> +#include <string> +#include <algorithm> + +#include "lbm.h" +#include "boundary_conditions.h" + +constexpr std::size_t dimX = 100; +constexpr std::size_t dimY = dimX; + +constexpr double uLid = 0.4; +constexpr double reynolds = 1000; + +constexpr double tau = 3. * uLid * (dimX-1) / reynolds + 0.5; +constexpr double omega = 1. / tau; + +DataCellBuffer pop(dimX, dimY); +FluidBuffer fluid(dimX, dimY); + +void init() { + for ( std::size_t x = 0; x < dimX; ++x ) { + for ( std::size_t y = 0; y < dimY; ++y ) { + fluid.density(x,y) = 1.0; + fluid.velocity(x,y) = { 0.0, 0.0 }; + + static_cast<Cell&>(pop.curr(x,y)).equilibrize( + fluid.density(x,y), fluid.velocity(x,y)); + static_cast<Cell&>(pop.prev(x,y)).equilibrize( + fluid.density(x,y), fluid.velocity(x,y)); + } + } +} + +void computeLbmStep() { + pop.swap(); + + for ( std::size_t x = 0; x < dimX; ++x ) { + for ( std::size_t y = 0; y < dimY; ++y ) { + streamFluidCell(pop, x, y); + } + } + + // straight wall cell bounce back + for ( std::size_t x = 0; x < dimX; ++x ) { + computeZouHeVelocityWallCell(pop, {x, dimY-1}, { 0,-1}, uLid); + } + for ( std::size_t x = 1; x < dimX-1; ++x ) { + computeWallCell(pop, {x, 0}, { 0, 1}); + } + for ( std::size_t y = 1; y < dimY-1; ++y ) { + computeWallCell(pop, {0, y}, { 1, 0}); + computeWallCell(pop, {dimX-1, y}, {-1, 0}); + } + + // edge wall cell bounce back + computeWallCell(pop, {0, 0 }, { 1, 1}); + computeWallCell(pop, {dimX-1, 0 }, {-1, 1}); + + for ( std::size_t x = 0; x < dimX; ++x ) { + for ( std::size_t y = 0; y < dimY; ++y ) { + Cell& cell = static_cast<Cell&>(pop.curr(x,y)); + fluid.density(x,y) = cell.sum(); + fluid.velocity(x,y) = cell.velocity(fluid.density(x,y)); + + collideFluidCell(omega, pop, fluid, x, y); + } + } +} + +int main() { + init(); + + 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 <= 10000; ++t ) { + computeLbmStep(); + + if ( t % 100 == 0 ) { + std::cout << "."; + std::cout.flush(); + fluid.writeAsVTK("result/data_t" + std::to_string(t) + ".vtk"); + } + } + + std::cout << std::endl; +} |