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

#include <forcing.h>

Collaboration diagram for olb::forcing::LiangTRT< Forced >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >:

Public Types
using	MomentaF = typename Forced<MOMENTA>::template type<DESCRIPTOR>
using	CollisionO = typename COLLISION::template type<DESCRIPTOR,Forced<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

Static Public Attributes
static constexpr bool	is_vectorizable

Detailed Description

template<template< typename > typename Forced = momenta::Forced>
template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM, typename COLLISION>
struct olb::forcing::LiangTRT< Forced >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >

Definition at line 205 of file forcing.h.

Member Typedef Documentation

CollisionO

template<template< typename > typename Forced = momenta::Forced>

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

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

Definition at line 207 of file forcing.h.

MomentaF

template<template< typename > typename Forced = momenta::Forced>

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

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

Definition at line 206 of file forcing.h.

Member Function Documentation

apply()

template<template< typename > typename Forced = momenta::Forced>

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

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

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

inline

Definition at line 217 of file forcing.h.

                                                                            {
      V u[DESCRIPTOR::d];
      MomentaF().computeU(cell, u);
      CollisionO().apply(cell, parameters);
      const V omega = parameters.template get<descriptors::OMEGA>();
      const auto force = cell.template getField<descriptors::FORCE>();
      const auto rho = cell.template getField<descriptors::RHO>();
      const auto gradRho = cell.template getField<descriptors::NABLARHO>();
 
      V fTerm[DESCRIPTOR::q], fTermPlus[DESCRIPTOR::q], fTermMinus[DESCRIPTOR::q];
 
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        V c_u{};
        for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
          c_u += descriptors::c<DESCRIPTOR>(iPop,iD)*u[iD];
        }
        fTerm[iPop] = 0.;
        for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
          fTerm[iPop] += descriptors::c<DESCRIPTOR>(iPop,iD) * force[iD] * rho +c_u * descriptors::c<DESCRIPTOR>(iPop,iD) * gradRho[iD];
        }
        fTerm[iPop] *= descriptors::t<V,DESCRIPTOR>(iPop)*descriptors::invCs2<V,DESCRIPTOR>();
      }
 
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        fTermPlus[iPop]  = V{0.5} * (fTerm[iPop] + fTerm[descriptors::opposite<DESCRIPTOR>(iPop)]);
        fTermMinus[iPop] = V{0.5} * (fTerm[iPop] - fTerm[descriptors::opposite<DESCRIPTOR>(iPop)]);
      }
 
      V tau = V{1} / omega;
      const V tau2 = parameters.template get<collision::ITRT::TAU_MINUS>();
      const auto nablaRho = cell.template getField<descriptors::NABLARHO>();
      V ratio = util::norm<DESCRIPTOR::d>(nablaRho) / parameters.template get<collision::ITRT::MAXNABLARHO>();
      const V omega2 = V{1} / (tau + ratio * (tau2 - tau));
 
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] += ( V{1} - omega * V{0.5} ) * fTermPlus[iPop]
                    + ( V{1} - omega2 * V{0.5} ) * fTermMinus[iPop];
      }
 
      return {rho, util::normSqr<V,DESCRIPTOR::d>(u)};
    };

References olb::descriptors::c(), olb::descriptors::invCs2(), olb::util::norm(), olb::util::normSqr(), olb::descriptors::opposite(), and olb::descriptors::t().

Here is the call graph for this function:

Member Data Documentation

is_vectorizable

template<template< typename > typename Forced = momenta::Forced>

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

bool olb::forcing::LiangTRT< Forced >::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::is_vectorizable

staticconstexpr

Initial value:

= dynamics::is_vectorizable_v<CollisionO>
      
      && !std::is_same_v<MOMENTA, momenta::IncBulkTuple<momenta::ForcedMomentum<momenta::IncompressibleBulkMomentum>>>

Definition at line 212 of file forcing.h.

---

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

• src/dynamics/forcing.h