Add Knudsen functor to locate refinement areas

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);