olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM > Struct Template Reference

#include <porousBGKdynamics.h>

Collaboration diagram for olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >:

Public Types
using	MomentaF = typename MOMENTA::template type<DESCRIPTOR>
using	EquilibriumF = typename EQUILIBRIUM::template type<DESCRIPTOR,MOMENTA>
using	CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Public Member Functions
template<typename CELL , typename pVELOCITY , typename V = typename CELL::value_t>
void	calculate (CELL &cell, pVELOCITY &pVelocity)
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters)

Static Public Attributes
static constexpr bool	is_vectorizable = false

Detailed Description

template<typename COLLISION, bool isStatic = false>
template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 110 of file porousBGKdynamics.h.

Member Typedef Documentation

CollisionO

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Definition at line 113 of file porousBGKdynamics.h.

EquilibriumF

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::EquilibriumF = typename EQUILIBRIUM::template type<DESCRIPTOR,MOMENTA>

Definition at line 112 of file porousBGKdynamics.h.

MomentaF

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 111 of file porousBGKdynamics.h.

Member Function Documentation

apply()

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>

CellStatistic< V > olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 133 of file porousBGKdynamics.h.

                                                               {
      V rho, u[DESCRIPTOR::d];
      MomentaF().computeRhoU(cell, rho, u);
 
      auto statistic = CollisionO().apply(cell, parameters);
 
      V velDenominator = cell.template getFieldComponent<descriptors::VELOCITY_DENOMINATOR>(0);
 
      // use Kuperstokh forcing by default
      V uPlus[DESCRIPTOR::d]{ };
      V diff[DESCRIPTOR::q]{ };
      V fEqPlus[DESCRIPTOR::q]{ };
      V fEq[DESCRIPTOR::q]{ };
 
      if (velDenominator > std::numeric_limits<V>::epsilon()) {
        for (int iDim=0; iDim<DESCRIPTOR::d; ++iDim) {
          uPlus[iDim] = u[iDim];
        }
        calculate(cell, uPlus);
        if constexpr (!isStatic) {
          particles::resetParticleRelatedFields<DESCRIPTOR,CELL,V>(cell);
        }
 
        EquilibriumF().compute(cell, rho, uPlus, fEqPlus);
        EquilibriumF().compute(cell, rho, u, fEq);
        for (int tmp_iPop=0; tmp_iPop < DESCRIPTOR::q; tmp_iPop++) {
          diff[tmp_iPop] += fEqPlus[tmp_iPop] - fEq[tmp_iPop];
          cell[tmp_iPop] += diff[tmp_iPop];
        }
      }
 
      particles::resetParticleContactRelatedFields<DESCRIPTOR,CELL,V>(cell);
 
      return statistic;
    }

References olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::calculate(), olb::particles::resetParticleContactRelatedFields(), and olb::particles::resetParticleRelatedFields().

Here is the call graph for this function:

calculate()

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

template<typename CELL , typename pVELOCITY , typename V = typename CELL::value_t>

void olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::calculate ( CELL & cell, pVELOCITY & pVelocity )

inline

Definition at line 118 of file porousBGKdynamics.h.

                                                      {
      if constexpr (isStatic) {
        for (int i=0; i<DESCRIPTOR::d; i++)  {
          pVelocity[i] -= (1.-(cell.template getField<descriptors::POROSITY>())) * pVelocity[i];
        }
      } else {
        for (int i=0; i<DESCRIPTOR::d; i++)  {
          pVelocity[i] += (1.-cell.template getField<descriptors::POROSITY>())
                        * (cell.template getFieldComponent<descriptors::VELOCITY_NUMERATOR>(i)
                        / cell.template getField<descriptors::VELOCITY_DENOMINATOR>() - pVelocity[i]);
        }
      }
    }

Here is the caller graph for this function:

Member Data Documentation

is_vectorizable

template<typename COLLISION , bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

bool olb::collision::PorousParticle< COLLISION, isStatic >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::is_vectorizable = false

staticconstexpr

Definition at line 115 of file porousBGKdynamics.h.

---

The documentation for this struct was generated from the following file:

• src/dynamics/porousBGKdynamics.h