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

#include <porousBGKdynamics.h>

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

Detailed Description

template<typename COLLISION>
template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::PSM< COLLISION >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 224 of file porousBGKdynamics.h.

Member Typedef Documentation

CollisionO

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 227 of file porousBGKdynamics.h.

EquilibriumF

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 226 of file porousBGKdynamics.h.

MomentaF

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 225 of file porousBGKdynamics.h.

Member Function Documentation

apply()

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

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

inline

Definition at line 230 of file porousBGKdynamics.h.

                                                               {
      V rho, u[DESCRIPTOR::d];
      const V epsilon = 1. - cell.template getField<descriptors::POROSITY>();
      const V omega = parameters.template get<descriptors::OMEGA>();
      const V paramA = V{1.} / omega - V{0.5};
      MomentaF().computeRhoU(cell, rho, u);
      // velocity at the boundary
      auto u_s = cell.template getField<descriptors::VELOCITY_SOLID>();
      if (epsilon < 1e-5) {
        return CollisionO().apply(cell, parameters);
      }
      else {
        // speed up paramB
        V paramB = (epsilon * paramA) / ((1. - epsilon) + paramA);
        // speed up paramC
        V paramC = (1. - paramB);
        V omega_s[DESCRIPTOR::q];
        V cell_tmp[DESCRIPTOR::q];
 
        V fEq[DESCRIPTOR::q] { };
        V fEq_s[DESCRIPTOR::q] { };
        EquilibriumF().compute(cell, rho, u, fEq);
        EquilibriumF().compute(cell, rho, u_s, fEq_s);
        for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
          cell_tmp[iPop] = cell[iPop];
          //if constexpr (MODE == 2) {
            omega_s[iPop] = (fEq_s[iPop] - cell[iPop]) + (V{1} - omega) * (cell[iPop] - fEq[iPop]);
          //} else if constexpr (MODE == 3) {
          //  omega_s[iPop] = (cell[descriptors::opposite<DESCRIPTOR>(iPop)] - fEq_s[descriptors::opposite<DESCRIPTOR>(iPop)])
          //                  - (cell[iPop] - fEq_s[iPop]);
          //}
        }
        CollisionO().apply(cell, parameters);
        for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
          cell[iPop] = cell_tmp[iPop] + paramC * (cell[iPop] - cell_tmp[iPop]);
          cell[iPop] += paramB * omega_s[iPop];
        }
        for (int iVel=0; iVel<DESCRIPTOR::d; ++iVel) {
          u[iVel] = paramC * u[iVel] + paramB * u_s[iVel];
        }
      }
 
      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/porousBGKdynamics.h