27#ifndef SET_SLIP_BOUNDARY_HH
28#define SET_SLIP_BOUNDARY_HH
35template <
typename,
typename,
int NX,
int NY,
int NZ>
44 template <
typename CELL,
typename V =
typename CELL::value_t>
46 using DESCRIPTOR =
typename CELL::descriptor_t;
47 int reflectionPop[DESCRIPTOR::q];
51 int mult = 2 / (NX*NX + NY*NY + NZ*NZ);
53 for (
int iPop = 1; iPop < DESCRIPTOR::q; iPop++) {
54 reflectionPop[iPop] = 0;
56 int scalarProduct = descriptors::c<DESCRIPTOR>(iPop,0)*NX + descriptors::c<DESCRIPTOR>(iPop,1)*NY + descriptors::c<DESCRIPTOR>(iPop,2)*NZ;
57 if ( scalarProduct < 0) {
60 mirrorDirection0 = -descriptors::c<DESCRIPTOR>(iPop,0);
61 mirrorDirection1 = -descriptors::c<DESCRIPTOR>(iPop,1);
62 mirrorDirection2 = -descriptors::c<DESCRIPTOR>(iPop,2);
65 mirrorDirection0 = descriptors::c<DESCRIPTOR>(iPop,0) - mult*scalarProduct*NX;
66 mirrorDirection1 = descriptors::c<DESCRIPTOR>(iPop,1) - mult*scalarProduct*NY;
67 mirrorDirection2 = descriptors::c<DESCRIPTOR>(iPop,2) - mult*scalarProduct*NZ;
71 for (
int i = 1; i < DESCRIPTOR::q; i++) {
72 if (descriptors::c<DESCRIPTOR>(i,0)==mirrorDirection0
73 && descriptors::c<DESCRIPTOR>(i,1)==mirrorDirection1
74 && descriptors::c<DESCRIPTOR>(i,2)==mirrorDirection2) {
75 reflectionPop[iPop] = i;
81 for (
int iPop = 1; iPop < DESCRIPTOR::q ; ++iPop) {
82 if (reflectionPop[iPop]!=0) {
84 x_b[iPop] = x_b[reflectionPop[iPop]];
92template<
typename T,
typename DESCRIPTOR>
99template<
typename T,
typename DESCRIPTOR>
102 bool includeOuterCells =
false;
105 if (indicator->getSuperGeometry().getOverlap() == 1) {
106 includeOuterCells =
true;
107 clout <<
"WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
109 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
113 addPoints2CommBC(sLattice, std::forward<
decltype(indicator)>(indicator), _overlap);
117template<
typename T,
typename DESCRIPTOR>
122 const int margin = includeOuterCells ? 0 : 1;
123 std::vector<int> discreteNormal(4, 0);
124 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
125 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
126 && indicator(iX, iY, iZ)) {
127 discreteNormal = blockGeometryStructure.getStatistics().getType(iX, iY, iZ);
128 if (discreteNormal[1]!=0 || discreteNormal[2]!=0 || discreteNormal[3]!=0) {
130 bool _output =
false;
132 clout <<
"setSlipBoundary<" << discreteNormal[1] <<
","<< discreteNormal[2] <<
","<< discreteNormal[3] <<
">(" << iX <<
", "<< iX <<
", " << iY <<
", " << iY <<
", " << iZ <<
", " << iZ <<
" )" << std::endl;
136 boundaryhelper::promisePostProcessorForNormal<T, DESCRIPTOR, FullSlipBoundaryPostProcessor3D>(
137 Vector <int,3> (discreteNormal.data()+1)
142 clout <<
"Warning: Could not setSlipBoundary (" << iX <<
", " << iY <<
", " << iZ <<
"), discreteNormal=(" << discreteNormal[0] <<
","<< discreteNormal[1] <<
","<< discreteNormal[2] <<
","<< discreteNormal[3] <<
"), set to bounceBack" << std::endl;
143 block.template defineDynamics<BounceBack>({iX, iY, iZ});
Base block indicator functor.
BlockGeometry< T, 3 > & getBlockGeometry()
Get underlying block geometry structure.
Platform-abstracted block lattice for external access and inter-block interaction.
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.
Smart pointer for managing the various ways of passing functors around.
class for marking output with some text
Representation of a statistic for a parallel 2D geometry.
std::unique_ptr< SuperIndicatorF< T, D > > getMaterialIndicator(std::vector< int > &&materials)
Returns a material indicator using the given vector of materials.
Base indicator functor (discrete)
Super class maintaining block lattices for a cuboid decomposition.
BlockLattice< T, DESCRIPTOR > & getBlock(int locC)
Return BlockLattice with local index locC.
LoadBalancer< T > & getLoadBalancer()
Read and write access to the load balancer.
Top level namespace for all of OpenLB.
void addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.
void setSlipBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 2 > &superGeometry, int material)
Initialising the SlipBoundary on the superLattice domain.
OperatorScope
Block-wide operator application scopes.
@ PerCell
Per-cell application, i.e. OPERATOR::apply is passed a CELL concept implementation.
void apply(CELL &x_b) any_platform
static constexpr OperatorScope scope
Communication after propagation.