26#ifndef SET_FREE_ENERGY_OUTLET_BOUNDARY_3D_HH
27#define SET_FREE_ENERGY_OUTLET_BOUNDARY_3D_HH
37template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
40 setFreeEnergyOutletBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice, omega, superGeometry.
getMaterialIndicator(material), type, latticeNumber);
46template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
51 bool includeOuterCells =
false;
52 if (indicator->getSuperGeometry().getOverlap() == 1) {
53 includeOuterCells =
true;
54 clout <<
"WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
56 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
57 setFreeEnergyOutletBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice.
getBlock(iCloc),
58 indicator->getBlockIndicatorF(iCloc), omega, type, latticeNumber, includeOuterCells);
61 addPoints2CommBC(sLattice, std::forward<
decltype(indicator)>(indicator), _overlap);
67template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
69 int latticeNumber,
bool includeOuterCells)
75 setFreeEnergyInletBoundary<T,DESCRIPTOR,MixinDynamics>(_block, omega, indicator, type, latticeNumber,includeOuterCells);
76 const int margin = includeOuterCells ? 0 : 1;
77 std::vector<int> discreteNormal(4, 0);
78 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
79 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
80 && indicator(iX, iY, iZ)) {
81 discreteNormal = blockGeometryStructure.getStatistics().getType(iX, iY, iZ,
true);
82 if (discreteNormal[0] == 0) {
85 olb::boundaryhelper::promisePostProcessorForNormal<T,DESCRIPTOR, FreeEnergyConvectiveProcessor3D>(
88 clout <<
"setFreeEnergyOutletBoundary<" <<
"," <<
">(" << iX <<
", "<< iY <<
", " << iZ <<
")" << std::endl;
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 setFreeEnergyOutletBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, T omega, SuperGeometry< T, 2 > &superGeometry, int material, std::string type, int latticeNumber)
Implementation of a outlet boundary condition for the partner lattices of the binary or ternary free ...
Communication after propagation.