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

#include <freeEnergyDynamics.h>

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

Public Types
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) any_platform

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::FreeEnergy::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 90 of file freeEnergyDynamics.h.

Member Typedef Documentation

EquilibriumF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

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

Definition at line 91 of file freeEnergyDynamics.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::FreeEnergy::type< DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 94 of file freeEnergyDynamics.h.

                                                                            {
      V fEq[DESCRIPTOR::q] { };
      const auto statistic = EquilibriumF().compute(cell, parameters, fEq);
      const V omega = parameters.template get<descriptors::OMEGA>();
      const V gamma = parameters.template get<GAMMA>();
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] *= V{1} - omega;
        cell[iPop] += omega * fEq[iPop];
      }
      const V tmp = gamma * descriptors::invCs2<V,DESCRIPTOR>()
                  * cell.template getField<descriptors::CHEM_POTENTIAL>();
      for (int iPop=1; iPop < DESCRIPTOR::q; ++iPop) {
        cell[iPop] -= omega * descriptors::t<V,DESCRIPTOR>(iPop) * (statistic.rho - tmp);
      }
      cell[0] += omega * (V{1} - descriptors::t<V,DESCRIPTOR>(0)) * (statistic.rho - tmp);
      return statistic;
    };

---

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

• src/dynamics/freeEnergyDynamics.h