44 OstreamManager clout(std::cout,
"setInterpolatedVelocityBoundary");
46 bool includeOuterCells =
false;
47 if (indicator->getSuperGeometry().getOverlap() == 1) {
48 includeOuterCells =
true;
49 clout <<
"WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
52 setInterpolatedVelocityBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice.
getBlock(iC),omega, indicator->getBlockIndicatorF(iC),includeOuterCells);
56 addPoints2CommBC(sLattice,std::forward<
decltype(indicator)>(indicator), _overlap);
64 using namespace boundaryhelper;
66 const int margin = includeOuterCells ? 0 : 1;
67 std::vector<int> discreteNormal(4,0);
68 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
69 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
70 && indicator(iX, iY, iZ)) {
72 discreteNormal = indicator.
getBlockGeometry().getStatistics().getType(iX, iY, iZ);
73 if (discreteNormal[0] == 0) {
74 if (discreteNormal[1] != 0 && discreteNormal[1] == -1) {
75 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
82 else if (discreteNormal[1] != 0 && discreteNormal[1] == 1) {
83 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
90 else if (discreteNormal[2] != 0 && discreteNormal[2] == -1) {
91 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
98 else if (discreteNormal[2] != 0 && discreteNormal[2] == 1) {
99 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
106 else if (discreteNormal[3] != 0 && discreteNormal[3] == -1) {
107 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
114 else if (discreteNormal[3] != 0 && discreteNormal[3] == 1) {
115 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
124 else if (discreteNormal[0] == 1) {
126 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
131 promisePostProcessorForNormal<T,DESCRIPTOR,OuterVelocityCornerProcessor3D>(
135 else if (discreteNormal[0] == 2) {
137 dynamics = _block.
getDynamics(PlainMixinDynamicsForNormalMomenta<T,DESCRIPTOR,
142 else if (discreteNormal[0] == 3) {
143 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
148 promisePostProcessorForNormalSpecial<T,DESCRIPTOR,OuterVelocityEdgeProcessor3D>(
152 else if (discreteNormal[0] == 4) {
153 dynamics = _block.
getDynamics(PlainMixinDynamicsForNormalSpecialMomenta<T,DESCRIPTOR,
159 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.
virtual Dynamics< T, DESCRIPTOR > * getDynamics(DynamicsPromise< T, DESCRIPTOR > &&)=0
Return pointer to dynamics yielded by promise.
void setInterpolatedVelocityBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material)
Initialising the setInterpolatedVelocityBoundary function on the superLattice domain Interpolated Bou...
void addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.
void setBoundary(BlockLattice< T, DESCRIPTOR > &block, int iX, int iY, Dynamics< T, DESCRIPTOR > *dynamics, PostProcessorGenerator2D< T, DESCRIPTOR > *postProcessor)
virtual AbstractParameters< T, DESCRIPTOR > & getParameters(BlockLattice< T, DESCRIPTOR > &block)=0
Parameters access for legacy post processors.
Communication after propagation.