24#ifndef LATTICE_PHYS_BOUNDARY_FORCE_2D_HH
25#define LATTICE_PHYS_BOUNDARY_FORCE_2D_HH
44template<
typename T,
typename DESCRIPTOR>
50 _indicatorF(std::move(indicatorF))
52 this->
getName() =
"physBoundaryForce";
53 for (
int iC = 0; iC < this->
_sLattice.getLoadBalancer().size(); ++iC) {
57 _indicatorF->getBlockIndicatorF(iC),
62template<
typename T,
typename DESCRIPTOR>
68 superGeometry.getMaterialIndicator(material),
72template <
typename T,
typename DESCRIPTOR>
78 _indicatorF(indicatorF),
79 _blockGeometry(indicatorF.getBlockGeometry())
81 this->
getName() =
"physBoundaryForce";
84template <
typename T,
typename DESCRIPTOR>
87 for (
int i = 0; i < this->getTargetDim(); ++i) {
91 if (_indicatorF(input)) {
92 for (
int iPop = 1; iPop < DESCRIPTOR::q; ++iPop) {
96 if (_blockGeometry.get({input[0] + descriptors::c<DESCRIPTOR >(iPop,0), input[1] + descriptors::c<DESCRIPTOR >(iPop,1)}) == 1) {
98 T f = this->_blockLattice.get(input[0] + descriptors::c<DESCRIPTOR >(iPop,0), input[1] + descriptors::c<DESCRIPTOR >(iPop,1))[iPop];
101 f += this->_blockLattice.get(input[0], input[1])[descriptors::opposite<DESCRIPTOR >(iPop)];
103 for (
int i = 0; i < this->getTargetDim(); ++i) {
104 output[i] -= descriptors::c<DESCRIPTOR >(iPop,i) * f;
108 for (
int i = 0; i < this->getTargetDim(); ++i) {
109 output[i] = this->_converter.getPhysForce(output[i]);
Base block indicator functor (discrete)
BlockLatticePhysBoundaryForce2D returns pointwise phys force acting on a boundary.
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysBoundaryForce2D(BlockLattice< T, DESCRIPTOR > &blockLattice, BlockIndicatorF2D< T > &indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Platform-abstracted block lattice for external access and inter-block interaction.
Smart pointer for managing the various ways of passing functors around.
std::string & getName()
read and write access to name
std::vector< std::unique_ptr< BlockF2D< T > > > _blockF
Super functors may consist of several BlockF2D<W> derived functors.
Representation of a statistic for a parallel 2D geometry.
SuperLattice< T, DESCRIPTOR > & _sLattice
functor to get pointwise phys force acting on a boundary with a given material on local lattice
SuperLatticePhysBoundaryForce2D(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
The description of a generic interface for all functor classes – header file.
This file contains indicator functions.
Wrapper functions that simplify the use of MPI.
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.