olb::collision::FreeEnergyInletOutlet< direction, orientation >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM > Struct Template Reference

#include <freeEnergyDynamics.h>

Collaboration diagram for olb::collision::FreeEnergyInletOutlet< direction, orientation >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >:

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

Public Member Functions
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<int direction, int orientation>
template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::FreeEnergyInletOutlet< direction, orientation >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 125 of file freeEnergyDynamics.h.

Member Typedef Documentation

MomentaF

template<int direction, int orientation>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::FreeEnergyInletOutlet< direction, orientation >::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 126 of file freeEnergyDynamics.h.

Member Function Documentation

apply()

template<int direction, int orientation>

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

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

inline

Definition at line 129 of file freeEnergyDynamics.h.

                                                                            {
      // Do a standard collision neglecting the chemical potential term.
      V rho, u[DESCRIPTOR::d];
      MomentaF().computeRhoU(cell, rho, u);
      const V omega = parameters.template get<descriptors::OMEGA>();
      V uSqr = lbm<DESCRIPTOR>::bgkCollision(cell, rho, u, omega);
      for (int iPop=1; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] -= omega * descriptors::t<V,DESCRIPTOR>(iPop) * rho;
      }
      cell[0] += omega * (V{1} - descriptors::t<V,DESCRIPTOR>(0)) * rho;
      // Distribute the missing density to the unknown distribution functions.
      constexpr auto missingIndices = util::subIndexOutgoing<DESCRIPTOR,direction,orientation>();
      V missingRho = rho - V{1};
      V missingWeightSum = 0;
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
 
        // std::find can conflict with CUDA
        /*if (std::find(missingIndices.begin(), missingIndices.end(), iPop) != missingIndices.end()) {
          missingWeightSum += descriptors::t<V,DESCRIPTOR>(iPop);
        } else {
          missingRho -= cell[iPop];
        }*/
 
        bool contains = false;
        for(long unsigned int i=0; i < missingIndices.size(); i++){
          if(missingIndices[i] == (long unsigned int)iPop){
            contains = true;
          }
        }
 
        if(contains){
          missingWeightSum += descriptors::t<V,DESCRIPTOR>(iPop);
        } else {
          missingRho -= cell[iPop];
        }
 
      }
 
      for (unsigned iPop=0; iPop < missingIndices.size(); ++iPop) {
        cell[missingIndices[iPop]] = missingRho * descriptors::t<V,DESCRIPTOR>(missingIndices[iPop]) / missingWeightSum;
      }
      return {rho, uSqr};
    };

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/freeEnergyDynamics.h