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

The velocity is stored in the external field U, except for the component "direction", which is computed by means of the population and the pressure. More...

#include <elements.h>

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

Classes
struct	VELOCITY

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 , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>
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::FixedPressureMomentum< direction, orientation >

The velocity is stored in the external field U, except for the component "direction", which is computed by means of the population and the pressure.

Definition at line 463 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::FixedPressureMomentum< direction, orientation >::compute ( CELL & cell, J & j )

inline

Definition at line 468 of file elements.h.

  {
    computeU<TYPE>(cell, j);
    const V rho = TYPE().computeRho(cell);
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      j[iD] *= rho;
    }
  }

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::FixedPressureMomentum< direction, orientation >::computeU ( CELL & cell, U & u )

inline

Definition at line 479 of file elements.h.

  {
    auto values = cell.template getField<VELOCITY>();
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      u[iD] = values[iD];
    }
    const V rho = TYPE().computeRho(cell);
 
    constexpr auto onWallIndices = util::populationsContributingToVelocity<DESCRIPTOR,direction,0>();
    constexpr auto normalIndices = util::populationsContributingToVelocity<DESCRIPTOR,direction,orientation>();
 
    V rhoOnWall = V{};
    for (auto e : onWallIndices) {
      rhoOnWall += cell[e];
    }
 
    V rhoNormal = V{};
    for (auto e : normalIndices) {
      rhoNormal += cell[e];
    }
 
    u[direction] = orientation * ((V{2}*rhoNormal+rhoOnWall+V{1}) / rho-V{1});
  }

define()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename U >

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

inline

Definition at line 505 of file elements.h.

  {
    cell.template setField<VELOCITY>(u);
  }

getName()

template<int direction, int orientation>

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

inlinestatic

Definition at line 520 of file elements.h.

                             {
    return "FixedPressureMomentum";
  }

initialize()

template<int direction, int orientation>

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

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

inline

Definition at line 511 of file elements.h.

  {
    V u[DESCRIPTOR::d] { V{} };
    cell.template setField<VELOCITY>(u);
  }

inverseShift()

template<int direction, int orientation>

template<typename TYPE , typename CELL , typename U >

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

inline

Definition at line 518 of file elements.h.

{};

---

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

• src/dynamics/momenta/elements.h