olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL > Struct Template Reference

#include <shanChenForcedSingleComponentCoupling.h>

Collaboration diagram for olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >:

Classes
struct	G
struct	OMEGA
struct	RHO0

Public Types
using	parameters = typename POTENTIAL::parameters::template include<G,RHO0,OMEGA>

Public Member Functions
int	getPriority () const
template<typename CELL , typename PARAMETERS >
void	apply (CELL &cell, PARAMETERS &parameters) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCellWithParameters

Detailed Description

template<typename T, typename DESCRIPTOR, typename POTENTIAL>
struct olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >

Definition at line 33 of file shanChenForcedSingleComponentCoupling.h.

Member Typedef Documentation

parameters

template<typename T , typename DESCRIPTOR , typename POTENTIAL >

using olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >::parameters = typename POTENTIAL::parameters::template include<G,RHO0,OMEGA>

Definition at line 40 of file shanChenForcedSingleComponentCoupling.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , typename POTENTIAL >

template<typename CELL , typename PARAMETERS >

void olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 47 of file shanChenForcedSingleComponentCoupling.h.

  {
    using V = typename CELL::value_t;
    //using DESCRIPTOR = typename CELL::descriptor_t;
 
    auto j = cell.template getField<descriptors::VELOCITY>();
    lbm<DESCRIPTOR>::computeJ(cell, j);
    cell.template setField<descriptors::VELOCITY>(j);
 
    const V g = parameters.template get<G>();
    const V rho0 = parameters.template get<RHO0>();
    const V omega = parameters.template get<OMEGA>();
 
    const V rho = cell.template getFieldComponent<descriptors::STATISTIC>(0) * rho0;
    // Computation of the common velocity, shared among the two populations
    const V rhoTot = rho * omega;
 
    Vector<V,DESCRIPTOR::d> uTot{};
    auto u = cell.template getField<descriptors::VELOCITY>();
    uTot = (u*rho0*omega) / rhoTot;
 
    // Computation of the interaction potential
    Vector<V,DESCRIPTOR::d> rhoContribution;
    const V psi = POTENTIAL().compute(rho, parameters);
 
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      const V rho_ = cell.neighbor(descriptors::c<DESCRIPTOR>(iPop))
                         .template getFieldComponent<descriptors::STATISTIC>(0)
                   * rho0;
      const V psi_ = POTENTIAL().compute(rho_, parameters);
      rhoContribution += psi * psi_
                       * descriptors::c<DESCRIPTOR>(iPop)
                       * descriptors::t<V,DESCRIPTOR>(iPop);
    }
 
    // Computation and storage of the final velocity, consisting
    //   of u and the momentum difference due to interaction
    //   potential plus external force
    auto force = cell.template getField<descriptors::EXTERNAL_FORCE>();
    cell.template setField<descriptors::VELOCITY>(uTot);
    cell.template setField<descriptors::FORCE>(force - g*rhoContribution/rho);
  }

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

Here is the call graph for this function:

getPriority()

template<typename T , typename DESCRIPTOR , typename POTENTIAL >

int olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >::getPriority ( ) const

inline

Definition at line 42 of file shanChenForcedSingleComponentCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , typename POTENTIAL >

constexpr OperatorScope olb::ShanChenForcedSingleComponentPostProcessor< T, DESCRIPTOR, POTENTIAL >::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 34 of file shanChenForcedSingleComponentCoupling.h.

---

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

• src/dynamics/shanChenForcedSingleComponentCoupling.h