summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-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);