64 using namespace boundaryhelper;
67 const int margin = includeOuterCells ? 0 : 1;
68 std::vector<int> discreteNormal(4,0);
69 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
70 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
71 && indicator(iX, iY, iZ)) {
73 discreteNormal = indicator.
getBlockGeometry().getStatistics().getType(iX, iY, iZ);
75 if (discreteNormal[1]!=0 || discreteNormal[2]!=0 || discreteNormal[3]!=0) {
78 clout <<
"setSlipBoundary<" << discreteNormal[1] <<
","<< discreteNormal[2] <<
","<< discreteNormal[3] <<
">(" << iX <<
", "<< iX <<
", " << iY <<
", " << iY <<
", " << iZ <<
", " << iZ <<
" )" << std::endl;
80 PostProcessorGenerator3D<T, DESCRIPTOR>* postProcessor =
new SlipBoundaryProcessorGenerator3D<T, DESCRIPTOR>(iX, iX, iY, iY, iZ, iZ, discreteNormal[1], discreteNormal[2], discreteNormal[3]);
86 clout <<
"Warning: Could not setSlipBoundary (" << iX <<
", " << iY <<
", " << iZ <<
"), discreteNormal=(" << discreteNormal[0] <<
","<< discreteNormal[1] <<
","<< discreteNormal[2] <<
","<< discreteNormal[3] <<
"), set to bounceBack" << std::endl;
87 block.template defineDynamics<BounceBack>({iX, iY, iZ});
90 if (discreteNormal[0] == 0) {
91 if (discreteNormal[1] != 0 && discreteNormal[1] == -1) {
92 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
96 else if (discreteNormal[1] != 0 && discreteNormal[1] == 1) {
97 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
101 else if (discreteNormal[2] != 0 && discreteNormal[2] == -1) {
102 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
106 else if (discreteNormal[2] != 0 && discreteNormal[2] == 1) {
107 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
111 else if (discreteNormal[3] != 0 && discreteNormal[3] == -1) {
112 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
116 else if (discreteNormal[3] != 0 && discreteNormal[3] == 1) {
117 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
122 else if (discreteNormal[0] == 1) {
123 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
127 else if (discreteNormal[0] == 2) {
128 dynamics = block.
getDynamics(PlainMixinDynamicsForNormalMomenta<T,DESCRIPTOR,
133 else if (discreteNormal[0] == 3) {
134 dynamics = block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
139 else if (discreteNormal[0] == 4) {
140 dynamics = block.
getDynamics(PlainMixinDynamicsForNormalSpecialMomenta<T,DESCRIPTOR,
144 dynamics->
getParameters(block).template set<descriptors::OMEGA>(omega);
virtual void addPostProcessor(std::type_index stage, LatticeR< DESCRIPTOR::d > latticeR, PostProcessorPromise< T, DESCRIPTOR > &&promise)=0
Schedule post processor for application to latticeR in stage.
void addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.
void setSlipBoundaryWithDynamics(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material)
Initialising the setSlipBoundaryWithDynamics function on the superLattice domain Interpolated Boundar...
void setBoundary(BlockLattice< T, DESCRIPTOR > &block, int iX, int iY, Dynamics< T, DESCRIPTOR > *dynamics, PostProcessorGenerator2D< T, DESCRIPTOR > *postProcessor)