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

#include <porousBGKdynamics.h>

Collaboration diagram for olb::collision::KrauseH< 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::KrauseH< COLLISION >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 346 of file porousBGKdynamics.h.

Member Typedef Documentation

CollisionO

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 349 of file porousBGKdynamics.h.

EquilibriumF

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 348 of file porousBGKdynamics.h.

MomentaF

template<typename COLLISION >

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 347 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::KrauseH< COLLISION >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Molecular realaxation time

Turbulent realaxation time

Effective realaxation time

Definition at line 352 of file porousBGKdynamics.h.

                                                               {
      V rho, u[DESCRIPTOR::d];
      V newOmega[DESCRIPTOR::d];
      V _fieldTmp[4];
      _fieldTmp[0] = V{1};
      _fieldTmp[1] = V{};
      _fieldTmp[2] = V{};
      _fieldTmp[3] = V{};
 
      const V smagoConst = parameters.template get<collision::LES::Smagorinsky>();
      V preFactor = smagoConst*smagoConst
                    * descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                    * 2*util::sqrt(2);
 
      MomentaF().computeRhoU(cell, rho, u);
 
      // compute newOmega
      V omega = parameters.template get<descriptors::OMEGA>();
      V uSqr = u[0]*u[0];
      for (int iDim=0; iDim<DESCRIPTOR::d; iDim++) {
        uSqr += u[iDim]*u[iDim];
      }
      V tau_mol = 1./omega;
      for (int iPop=0; iPop<DESCRIPTOR::q; iPop++) {
        V fNeq = util::fabs(cell[iPop] - EquilibriumF().compute(iPop, rho, u));
        V tau_turb = 0.5*(util::sqrt(tau_mol*tau_mol+(preFactor*fNeq))-tau_mol);
        V tau_eff = tau_mol + tau_turb;
        newOmega[iPop] = 1./tau_eff;
      }
      parameters.template set<descriptors::OMEGA>(newOmega);
 
 
      V vel_denom = cell.template getField<descriptors::VELOCITY_DENOMINATOR>();
      if (vel_denom > std::numeric_limits<V>::epsilon()) {
        V porosity = cell.template getField<descriptors::POROSITY>(); // prod(1-smoothInd)
        auto vel_num = cell.template getFieldPointer<descriptors::VELOCITY_NUMERATOR>();
        porosity = 1.-porosity; // 1-prod(1-smoothInd)
        for (int i=0; i<DESCRIPTOR::d; i++)  {
          u[i] += porosity * (vel_num[i] / vel_denom - u[i]);
        }
      }
      uSqr = CollisionO().apply(cell, parameters);
 
      cell.template setField<descriptors::POROSITY>(_fieldTmp[0]);
      cell.template setField<descriptors::VELOCITY_NUMERATOR>({_fieldTmp[1], _fieldTmp[2]});
      cell.template setField<descriptors::VELOCITY_DENOMINATOR>(_fieldTmp[3]);
 
      return {rho, uSqr};
    }

References olb::util::fabs(), and olb::util::sqrt().

Here is the call graph for this function:

---

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

• src/dynamics/porousBGKdynamics.h