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

#include <porousBGKdynamics.h>

Collaboration diagram for olb::collision::DBBParticleBGK::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >:

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

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 DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::DBBParticleBGK::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 283 of file porousBGKdynamics.h.

Member Typedef Documentation

EquilibriumF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 285 of file porousBGKdynamics.h.

MomentaF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 284 of file porousBGKdynamics.h.

Member Function Documentation

apply()

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

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

inline

Definition at line 288 of file porousBGKdynamics.h.

                                                               {
      V rho, u[DESCRIPTOR::d],eta[DESCRIPTOR::q],uPlus[DESCRIPTOR::d],tmp_cell[(DESCRIPTOR::q+1)/2];
      MomentaF().computeRhoU(cell, rho, u);
      const V omega = parameters.template get<descriptors::OMEGA>();
      V tmpMomentumLoss[DESCRIPTOR::d]{ };
      auto velNumerator   = cell.template getFieldPointer<descriptors::VELOCITY_NUMERATOR>();
      auto zeta           = cell.template getFieldPointer<descriptors::ZETA>();
      auto velDenominator = cell.template getField<descriptors::VELOCITY_DENOMINATOR>();
 
      if (velDenominator > 1) {
        rho /= velDenominator;
      }
      for (int tmp_iPop=1; 2*tmp_iPop<DESCRIPTOR::q; tmp_iPop++) {
        eta[tmp_iPop]=6.*descriptors::t<V,DESCRIPTOR>(tmp_iPop)*rho*(descriptors::c<DESCRIPTOR>(tmp_iPop,0)*(velNumerator[0])+descriptors::c<DESCRIPTOR>(tmp_iPop,1)*(velNumerator[1]));
        tmp_cell[tmp_iPop]=(zeta[tmp_iPop]*(-cell[tmp_iPop]+cell[descriptors::opposite<DESCRIPTOR>(tmp_iPop)]+eta[tmp_iPop]));
        cell[tmp_iPop]+=tmp_cell[tmp_iPop]/(1.+2.*(zeta[tmp_iPop]));
        cell[descriptors::opposite<DESCRIPTOR>(tmp_iPop)]-=tmp_cell[tmp_iPop]/(1.+2.*(zeta[tmp_iPop]));
        zeta[tmp_iPop] = 0.;
        zeta[descriptors::opposite<DESCRIPTOR>(tmp_iPop)] = 0.;
      }
 
      cell.template setField<descriptors::POROSITY>(1.);
      cell.template setField<descriptors::VELOCITY_DENOMINATOR>(0);
 
      MomentaF().computeRhoU(cell, rho, uPlus);
      V uPlusSqr = lbm<DESCRIPTOR>::bgkCollision(cell, rho, uPlus, omega);
 
      V diff[DESCRIPTOR::q] = {};
      for (int tmp_iPop=0; tmp_iPop<DESCRIPTOR::q; tmp_iPop++) {
        diff[tmp_iPop] += EquilibriumF().compute(tmp_iPop, rho, uPlus)
                          - EquilibriumF().compute(tmp_iPop, rho, u);
 
        for (int iDim=0; iDim<DESCRIPTOR::d; iDim++) {
          tmpMomentumLoss[iDim] -= descriptors::c<DESCRIPTOR>(tmp_iPop,iDim) * diff[tmp_iPop];
        }
      }
 
      for (int i_dim=0; i_dim<DESCRIPTOR::d; i_dim++) {
        velNumerator[i_dim] = tmpMomentumLoss[i_dim];
      }
 
      return {rho, util::normSqr<V,DESCRIPTOR::d>(u)};
    }

References olb::lbm< DESCRIPTOR >::bgkCollision().

Here is the call graph for this function:

---

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

• src/dynamics/porousBGKdynamics.h