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

#include <forcing.h>

Collaboration diagram for olb::forcing::ForcedHLBM::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 PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM, typename COLLISION>
struct olb::forcing::ForcedHLBM::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >

Definition at line 649 of file forcing.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 652 of file forcing.h.

EquilibriumF

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

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

Definition at line 651 of file forcing.h.

MomentaF

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

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

Definition at line 650 of file forcing.h.

Member Function Documentation

apply()

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

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

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

inline

Definition at line 655 of file forcing.h.

                                                                            {
      Vector<V,DESCRIPTOR::d> u{};
      MomentaF().computeU(cell, u);
 
      auto statistic = CollisionO().apply(cell, parameters);
 
      const V porosity = cell.template getField<descriptors::POROSITY>();
      const auto force = cell.template getField<descriptors::FORCE>();
 
      auto fHLBM = (1 - porosity) * (cell.template getField<descriptors::VELOCITY>() - u);
 
      // This is Kuperstokh forcing
      Vector<V,DESCRIPTOR::d> uPlus = u + force + fHLBM;
 
      V fEq[DESCRIPTOR::q] { };
      V fEq2[DESCRIPTOR::q] { };
      EquilibriumF().compute(cell, statistic.rho, u, fEq);
      EquilibriumF().compute(cell, statistic.rho, uPlus, fEq2);
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] += fEq2[iPop] - fEq[iPop];
      }
 
      if (porosity < 1) {
        return {-1,-1};
      } else {
        V rho{};
        MomentaF().computeRhoU(cell, rho, u);
        return {rho, util::normSqr<V,DESCRIPTOR::d>(u)};
      }
    };

References olb::util::normSqr().

Here is the call graph for this function:

---

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

• src/dynamics/forcing.h