olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION > Struct Template Reference

#include <porousForcedBGKDynamics.h>

Collaboration diagram for olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >:

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

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

Static Public Attributes
static constexpr bool	is_vectorizable = false

Detailed Description

template<bool isStatic = false>
template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM, typename COLLISION>
struct olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >

Definition at line 108 of file porousForcedBGKDynamics.h.

Member Typedef Documentation

CollisionO

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Definition at line 111 of file porousForcedBGKDynamics.h.

EquilibriumF

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::EquilibriumF = combined_equilibrium<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Definition at line 110 of file porousForcedBGKDynamics.h.

MomentaF

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 109 of file porousForcedBGKDynamics.h.

Member Function Documentation

apply()

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

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

CellStatistic< V > olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 131 of file porousForcedBGKDynamics.h.

                                                                            {
      V u[DESCRIPTOR::d];
      MomentaF().computeU(cell, u);
      const auto velDenominator = cell.template getField<descriptors::VELOCITY_DENOMINATOR>();
      const auto statistic = CollisionO().apply(cell, parameters);
      const auto force = cell.template getField<descriptors::FORCE>();
      V uPlusDeltaU[DESCRIPTOR::d];
      for (int iVel=0; iVel < DESCRIPTOR::d; ++iVel) {
        uPlusDeltaU[iVel] = u[iVel] + force[iVel];
      }
      if (velDenominator > std::numeric_limits<V>::epsilon()) {
        calculate(cell, uPlusDeltaU);
        if constexpr (!isStatic) {
          // reset external field for next timestep
          cell.template setField<descriptors::POROSITY>(1.);
          cell.template setField<descriptors::VELOCITY_DENOMINATOR>(0.);
          cell.template setField<descriptors::VELOCITY_NUMERATOR>({0.,0.,0.});
        }
      }
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] += EquilibriumF().compute(iPop, statistic.rho, uPlusDeltaU)
                    - EquilibriumF().compute(iPop, statistic.rho, u);
      }
 
      // Reset contact helper if utilized
      if constexpr (DESCRIPTOR::template provides<descriptors::CONTACT_DETECTION>()) {
        cell.template setField<descriptors::CONTACT_DETECTION>(0);
      }
 
      return statistic;
    };

References olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::calculate().

Here is the call graph for this function:

calculate()

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

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

void olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::calculate ( CELL & cell, VELOCITY & u )

inline

Definition at line 116 of file porousForcedBGKDynamics.h.

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

Here is the caller graph for this function:

Member Data Documentation

is_vectorizable

template<bool isStatic = false>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

constexpr bool olb::forcing::PorousParticleKupershtokh< isStatic >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::is_vectorizable = false

staticconstexpr

Definition at line 113 of file porousForcedBGKDynamics.h.

---

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

• src/dynamics/porousForcedBGKDynamics.h