olb::momenta::BulkPressure< MOMENTUM > Struct Template Reference

Computation of hydrodynamic pressure as zeroth moment of the population. More...

#include <elements.h>

Collaboration diagram for olb::momenta::BulkPressure< MOMENTUM >:

Public Member Functions
template<typename TYPE , typename CELL , typename PRESSURE , typename DESCRIPTOR = typename CELL::descriptor_t, typename V = typename CELL::value_t>
void	compute (CELL &cell, PRESSURE &p) any_platform
template<typename TYPE , typename CELL , typename RHO >
void	define (CELL &cell, const RHO &rho) any_platform
template<typename TYPE , typename CELL >
void	initialize (CELL &cell) any_platform
template<typename TYPE , typename CELL , typename RHO >
void	inverseShift (CELL &cell, RHO &rho) any_platform

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

Detailed Description

template<typename MOMENTUM>
struct olb::momenta::BulkPressure< MOMENTUM >

Computation of hydrodynamic pressure as zeroth moment of the population.

Definition at line 173 of file elements.h.

Member Function Documentation

compute()

template<typename MOMENTUM >

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

void olb::momenta::BulkPressure< MOMENTUM >::compute ( CELL & cell, PRESSURE & p )

inline

Definition at line 175 of file elements.h.

  {
    V sum = 0;
    for (int iPop=1; iPop < DESCRIPTOR::q; ++iPop) {
      sum += cell[iPop];
    }
    auto rho = cell.template getField<descriptors::RHO>();
    auto gradRho = cell.template getField<descriptors::NABLARHO>();
    Vector<V,DESCRIPTOR::d> u{};
    V uGradRho = 0;
    MOMENTUM().template computeU<TYPE>(cell, u);
    for (unsigned iD=0; iD < DESCRIPTOR::d; ++iD) {
      uGradRho += u[iD] * gradRho[iD];
    }
    V uSqr = util::normSqr<V,DESCRIPTOR::d>(u);
    const V s0 = -descriptors::t<V,DESCRIPTOR>(0) * descriptors::invCs2<V,DESCRIPTOR>()/V{2} * uSqr;
    p = 1/ descriptors::invCs2<V,DESCRIPTOR>() / (1-descriptors::t<V,DESCRIPTOR>(0)) * (sum + 0.5*uGradRho + rho*s0);
  }

References olb::descriptors::invCs2(), olb::util::normSqr(), and olb::descriptors::t().

Here is the call graph for this function:

define()

template<typename MOMENTUM >

template<typename TYPE , typename CELL , typename RHO >

void olb::momenta::BulkPressure< MOMENTUM >::define ( CELL & cell, const RHO & rho )

inline

Definition at line 195 of file elements.h.

{};

getName()

template<typename MOMENTUM >

static std::string olb::momenta::BulkPressure< MOMENTUM >::getName ( )

inlinestatic

Definition at line 203 of file elements.h.

                            {
    return "BulkPressure";
  }

initialize()

template<typename MOMENTUM >

template<typename TYPE , typename CELL >

void olb::momenta::BulkPressure< MOMENTUM >::initialize ( CELL & cell )

inline

Definition at line 198 of file elements.h.

{};

inverseShift()

template<typename MOMENTUM >

template<typename TYPE , typename CELL , typename RHO >

void olb::momenta::BulkPressure< MOMENTUM >::inverseShift ( CELL & cell, RHO & rho )

inline

Definition at line 201 of file elements.h.

{};

---

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

• src/dynamics/momenta/elements.h