43 OstreamManager clout(std::cout,
"setInterpolatedConvectionBoundary");
45 bool includeOuterCells =
false;
46 if (indicator->getSuperGeometry().getOverlap() == 1) {
47 includeOuterCells =
true;
48 clout <<
"WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
50 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
51 setInterpolatedConvectionBoundary<T,DESCRIPTOR>(sLattice.
getBlock(iCloc), omega,
52 indicator->getBlockIndicatorF(iCloc), uAv, includeOuterCells);
55 addPoints2CommBC<T, DESCRIPTOR>(sLattice, std::forward<
decltype(indicator)>(indicator), _overlap);
61 T* uAv,
bool includeOuterCells)
64 const int margin = includeOuterCells ? 0 : 1;
65 OstreamManager clout(std::cout,
"setInterpolatedConvectionBoundary");
67 std::vector<int> discreteNormal(3, 0);
68 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY) {
69 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY}) >= margin
70 && indicator(iX, iY)) {
72 discreteNormal = indicator.
getBlockGeometry().getStatistics().getType(iX, iY);
73 if (discreteNormal[0] == 0) {
74 if (discreteNormal[1] == -1) {
76 clout <<
"setInterpolatedConvectionBoundary<" << 0 <<
","<< -1 <<
">(" << iX <<
", "<< iY <<
", " << omega <<
" )" << std::endl;
79 < T,DESCRIPTOR,0,-1 > (iX, iX, iY, iY, uAv);
81 else if (discreteNormal[1] == 1) {
83 clout <<
"setInterpolatedConvectionBoundary<" << 0 <<
"," << 1 <<
">(" << iX <<
", "<< iY <<
", " << omega <<
" )" << std::endl;
86 < T,DESCRIPTOR,0,1 > (iX, iX, iY, iY, uAv);
88 else if (discreteNormal[2] == -1) {
90 clout <<
"setInterpolatedConvectionBoundary<" << 1 <<
"," << -1 <<
">(" << iX <<
", "<< iY <<
", " << omega <<
" )" << std::endl;
93 < T,DESCRIPTOR,1,-1 > (iX, iX, iY, iY, uAv);
95 else if (discreteNormal[2] == 1) {
97 clout <<
"setInterpolatedConvectionBoundary<" << 1 <<
"," << 1 <<
">(" << iX <<
", "<< iY <<
", " << omega <<
" )" << std::endl;
100 < T,DESCRIPTOR,1,1 > (iX, iX, iY, iY, uAv);
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 setInterpolatedConvectionBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material, T *uAv)
Initialising the InterpolatedConvectionBoundary on the superLattice domain This is an interpolated bo...