olb::momenta::FixedTemperatureMomentum< direction, orientation > Struct Template Reference

The conduction is computed from density and population. More...

#include <elements.h>

Collaboration diagram for olb::momenta::FixedTemperatureMomentum< direction, orientation >:

Public Member Functions
template<typename TYPE , typename CELL , typename J , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>
void	compute (CELL &cell, J &j) any_platform
template<typename TYPE , typename CELL , typename U , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>
void	computeU (CELL &cell, U &u) any_platform
template<typename TYPE , typename CELL , typename U >
void	define (CELL &cell, const U &u) any_platform
template<typename TYPE , typename CELL >
void	initialize (CELL &cell) any_platform
template<typename TYPE , typename CELL , typename U >
void	inverseShift (CELL &cell, U &u) any_platform

Static Public Member Functions
static std::string	getName ()

Detailed Description

template<int direction, int orientation>
struct olb::momenta::FixedTemperatureMomentum< direction, orientation >

The conduction is computed from density and population.

Be careful that computeU(...) computes the conduction whereas compute(...) computes the transport (= conduction + convective transport)

Definition at line 575 of file elements.h.

Member Function Documentation

compute()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename J , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>

void olb::momenta::FixedTemperatureMomentum< direction, orientation >::compute ( CELL & cell, J & j )

inline

Definition at line 578 of file elements.h.

  {
    const V temp = TYPE().computeRho(cell);
    const auto uNS = cell.template getFieldPointer<descriptors::VELOCITY>();
    computeU<TYPE>(cell, j);
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      j[iD] += temp * uNS[iD];
    }
  }

computeU()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename U , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>

void olb::momenta::FixedTemperatureMomentum< direction, orientation >::computeU ( CELL & cell, U & u )

inline

Definition at line 590 of file elements.h.

  {
    constexpr auto onWallIndices = util::populationsContributingToVelocity<DESCRIPTOR,direction,0>();
    constexpr auto normalIndices = util::populationsContributingToVelocity<DESCRIPTOR,direction,orientation>();
 
    const V temp = TYPE().computeRho(cell);
    const auto uNS = cell.template getField<descriptors::VELOCITY>();
 
    V uOnWall[DESCRIPTOR::d] = { V{} };
    for (unsigned fIndex=0; fIndex<onWallIndices.size(); ++fIndex) {
      for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
        uOnWall[iD] += (cell[onWallIndices[fIndex]]
                          - equilibrium<DESCRIPTOR>::firstOrder(onWallIndices[fIndex],temp,uNS.data()))
                          * descriptors::c<DESCRIPTOR>(onWallIndices[fIndex],iD);
      }
    }
 
    V uNormal[DESCRIPTOR::d] = { V{} };
    for (unsigned fIndex=0; fIndex<normalIndices.size(); ++fIndex) {
      for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
        uNormal[iD] += (cell[normalIndices[fIndex]]
                          - equilibrium<DESCRIPTOR>::firstOrder(normalIndices[fIndex],temp,uNS.data()))
                          * descriptors::c<DESCRIPTOR>(normalIndices[fIndex],iD);
      }
    }
 
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      u[iD] = uOnWall[iD] + V(2) * uNormal[iD];
    }
  }

References olb::equilibrium< DESCRIPTOR >::firstOrder().

Here is the call graph for this function:

define()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename U >

void olb::momenta::FixedTemperatureMomentum< direction, orientation >::define ( CELL & cell, const U & u )

inline

Definition at line 623 of file elements.h.

{};

getName()

template<int direction, int orientation>

static std::string olb::momenta::FixedTemperatureMomentum< direction, orientation >::getName ( )

inlinestatic

Definition at line 631 of file elements.h.

                             {
    return "FixedTemperatureMomentum";
  }

initialize()

template<int direction, int orientation>

template<typename TYPE , typename CELL >

void olb::momenta::FixedTemperatureMomentum< direction, orientation >::initialize ( CELL & cell )

inline

Definition at line 626 of file elements.h.

{};

inverseShift()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename U >

void olb::momenta::FixedTemperatureMomentum< direction, orientation >::inverseShift ( CELL & cell, U & u )

inline

Definition at line 629 of file elements.h.

{};

---

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

• src/dynamics/momenta/elements.h