diff options
| -rw-r--r-- | apps/adrian/cylinder2d/cylinder2d.cpp | 12 |
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); |