44 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;
51 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
52 setFreeEnergyInletBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice.
getBlock(iCloc),omega,
53 indicator->getBlockIndicatorF(iCloc), type, latticeNumber, includeOuterCells);
56 addPoints2CommBC(sLattice, std::forward<
decltype(indicator)>(indicator), _overlap);
63 int latticeNumber,
bool includeOuterCells)
66 const int margin = includeOuterCells ? 0 : 1;
69 std::vector<int> discreteNormal(4, 0);
70 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
71 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
72 && indicator(iX, iY, iZ)) {
73 discreteNormal = blockGeometryStructure.getStatistics().getType(iX, iY, iZ,
true);
75 if (discreteNormal[0] == 0) {
77 if (discreteNormal[1] != 0 && discreteNormal[1] == -1) {
78 if (latticeNumber == 1) {
79 if (type ==
"density") {
84 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
94 else if (discreteNormal[1] != 0 && discreteNormal[1] == 1) {
95 if (latticeNumber == 1) {
96 if (type ==
"density") {
101 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
107 dynamics = _block.template getDynamics<FreeEnergyInletOutletDynamics<T,DESCRIPTOR,0,1>>();
111 else if (discreteNormal[2] != 0 && discreteNormal[2] == -1) {
112 if (latticeNumber == 1) {
113 if (type ==
"density") {
118 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
128 else if (discreteNormal[2] != 0 && discreteNormal[2] == 1) {
129 if (latticeNumber == 1) {
130 if (type ==
"density") {
135 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
141 dynamics = _block.template getDynamics<FreeEnergyInletOutletDynamics<T,DESCRIPTOR,1,1>>();
145 else if (discreteNormal[3] != 0 && discreteNormal[3] == -1) {
146 if (latticeNumber == 1) {
147 if (type ==
"density") {
152 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
162 else if (discreteNormal[3] != 0 && discreteNormal[3] == 1) {
163 if (latticeNumber == 1) {
164 if (type ==
"density") {
169 dynamics = _block.template getDynamics<
typename MixinDynamics::template exchange_momenta<
175 dynamics = _block.template getDynamics<FreeEnergyInletOutletDynamics<T,DESCRIPTOR,2,1>>();
178 dynamics->
getParameters(_block).template set<descriptors::OMEGA>(omega);
179 if (latticeNumber != 1) {
181 auto cell = _block.
get(iX,iY,iZ);
186 if (latticeNumber == 1){
189 olb::boundaryhelper::promisePostProcessorForNormal<T,DESCRIPTOR, FreeEnergyChemPotBoundaryProcessor3DA>(
194 olb::boundaryhelper::promisePostProcessorForNormal<T,DESCRIPTOR, FreeEnergyChemPotBoundaryProcessor3DB>(
199 clout <<
"setFreeEnergyInletBoundary<" <<
"," <<
">(" << iX <<
", "<< iY <<
", " << iZ <<
")" << std::endl;
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 setFreeEnergyInletBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material, std::string type, int latticeNumber)
Implementation of a inlet boundary condition for the partner lattices of the binary or ternary free e...
Communication after propagation.