summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-01-21 20:39:38 +0100
committerAdrian Kummerlaender2019-06-24 15:17:09 +0200
commit6372c796c62d60d2642ffc81365ba9e7a87fdc68 (patch)
tree308614c69e69bf3ef3eb28f2b361c3a44fc10c3f
parent80e945199941d2c80ef8ff6016f4406f3b37869c (diff)
downloadgrid_refinement_openlb-6372c796c62d60d2642ffc81365ba9e7a87fdc68.tar
grid_refinement_openlb-6372c796c62d60d2642ffc81365ba9e7a87fdc68.tar.gz
grid_refinement_openlb-6372c796c62d60d2642ffc81365ba9e7a87fdc68.tar.bz2
grid_refinement_openlb-6372c796c62d60d2642ffc81365ba9e7a87fdc68.tar.xz
grid_refinement_openlb-6372c796c62d60d2642ffc81365ba9e7a87fdc68.zip
Add Knudsen functor to locate refinement areas
-rw-r--r--apps/adrian/cylinder2d/cylinder2d.cpp12
1 files changed, 7 insertions, 5 deletions
diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index 3baa4ef..2e18ebd 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -171,12 +171,14 @@ void getResults(const std::string& prefix,
auto& sGeometry = grid.getSuperGeometry();
SuperVTMwriter2D<T> vtmWriter(prefix + "cylinder2d");
- SuperLatticePhysVelocity2D<T, DESCRIPTOR> velocity(sLattice, converter);
- SuperLatticePhysPressure2D<T, DESCRIPTOR> pressure(sLattice, converter);
- SuperLatticeGeometry2D<T, DESCRIPTOR> geometry(sLattice, sGeometry);
+ SuperLatticePhysVelocity2D<T,DESCRIPTOR> velocity(sLattice, converter);
+ SuperLatticePhysPressure2D<T,DESCRIPTOR> pressure(sLattice, converter);
+ SuperLatticeGeometry2D<T,DESCRIPTOR> geometry(sLattice, sGeometry);
+ SuperLatticeKnudsen2D<T,DESCRIPTOR> knudsen(sLattice);
vtmWriter.addFunctor(geometry);
vtmWriter.addFunctor(velocity);
vtmWriter.addFunctor(pressure);
+ vtmWriter.addFunctor(knudsen);
const int statIter = converter.getLatticeTime(maxPhysT/10.);
@@ -235,8 +237,8 @@ int main(int argc, char* argv[])
sOnLatticeBoundaryCondition2D<T, DESCRIPTOR> fineSBoundaryCondition(fineGrid.getSuperLattice());
createLocalBoundaryCondition2D<T, DESCRIPTOR>(fineSBoundaryCondition);
- const Vector<T,2> fineExtend2 {5*cylinderD, 1.5*cylinderD};
- const Vector<T,2> fineOrigin2 {1.25*cylinderD, (ly-fineExtend2[1])/2};
+ const Vector<T,2> fineExtend2 {1.75*cylinderD, 1.75*cylinderD};
+ const Vector<T,2> fineOrigin2 {2*cylinderD-fineExtend2[0]/2, 2*cylinderD-fineExtend2[1]/2};
auto& fineGrid2 = fineGrid.refine(fineOrigin2, fineExtend2);
prepareGeometry(fineGrid2);