diff options
Diffstat (limited to 'apps/adrian/poiseuille2d/poiseuille2d.cpp')
| -rw-r--r-- | apps/adrian/poiseuille2d/poiseuille2d.cpp | 12 |
1 files changed, 4 insertions, 8 deletions
diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp index 6d0c86f..b2fb48b 100644 --- a/apps/adrian/poiseuille2d/poiseuille2d.cpp +++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp @@ -105,13 +105,13 @@ void prepareLattice(Grid2D<T,DESCRIPTOR>& grid, sLattice.defineDynamics(sGeometry, 0, &instances::getNoDynamics<T, DESCRIPTOR>()); sLattice.defineDynamics(sGeometry, 1, &bulkDynamics); // bulk - sLattice.defineDynamics(sGeometry, 2, &instances::getNoDynamics<T, DESCRIPTOR>()); + sLattice.defineDynamics(sGeometry, 2, &bulkDynamics); // walls sLattice.defineDynamics(sGeometry, 3, &bulkDynamics); // inflow sLattice.defineDynamics(sGeometry, 4, &bulkDynamics); // outflow - sBoundaryCondition.addVelocityBoundary(sGeometry, 2, omega); - sBoundaryCondition.addVelocityBoundary(sGeometry, 3, omega); - sBoundaryCondition.addPressureBoundary(sGeometry, 4, omega); + sBoundaryCondition.addVelocityBoundary(sGeometry, 2, omega); // 0-velocity walls + sBoundaryCondition.addVelocityBoundary(sGeometry, 3, omega); // velocity inflow + sBoundaryCondition.addPressureBoundary(sGeometry, 4, omega); // pressure outflow const T Lx = converter.getLatticeLength(lx); const T Ly = converter.getLatticeLength(ly); @@ -201,10 +201,6 @@ int main(int argc, char* argv[]) prepareGeometry(fineGrid); - auto refinedOverlap = fineGrid.getRefinedOverlap(); - coarseGrid.getSuperGeometry().rename(1,0,*refinedOverlap); - coarseGrid.getSuperGeometry().rename(2,0,*refinedOverlap); - BGKdynamics<T, DESCRIPTOR> coarseBulkDynamics( coarseGrid.getConverter().getLatticeRelaxationFrequency(), instances::getBulkMomenta<T, DESCRIPTOR>()); |