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

#include <elements.h>

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

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 , typename DESCRIPTOR = typename CELL::descriptor_t>
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 , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>
void	inverseShift (CELL &cell, U &u) any_platform

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

Detailed Description

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

Definition at line 1125 of file elements.h.

Member Function Documentation

compute()

template<typename MOMENTUM >

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

void olb::momenta::GuoZhaoForcedMomentum< MOMENTUM >::compute ( CELL & cell, J & j )

inline

Definition at line 1127 of file elements.h.

  {
    MOMENTUM().template compute<TYPE>(cell, j);
  }

computeU()

template<typename MOMENTUM >

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

void olb::momenta::GuoZhaoForcedMomentum< MOMENTUM >::computeU ( CELL & cell, U & u )

inline

Definition at line 1133 of file elements.h.

  {
    MOMENTUM().template computeU<TYPE>(cell, u);
    const V epsilon = cell.template getField<descriptors::EPSILON>();
    const V nu      = cell.template getField<descriptors::NU>();
    const V k       = cell.template getField<descriptors::K>();
    const auto bodyF = cell.template getFieldPointer<descriptors::BODY_FORCE>();
 
    const V uMag = util::sqrt( util::normSqr<V,DESCRIPTOR::d>(u) );
    const V Fe = 0.;//1.75/util::sqrt(150.*util::pow(epsilon,3));
    const V c_0 = 0.5*(1 + 0.5*epsilon*nu/k);
    const V c_1 = 0.5*epsilon*Fe/util::sqrt(k);
 
    for (int iVel=0; iVel < DESCRIPTOR::d; ++iVel) {
      u[iVel] += bodyF[iVel] * V(0.5) * epsilon;
      u[iVel] /= (c_0 + util::sqrt(c_0*c_0 + c_1*uMag));
    }
  }

References olb::util::sqrt().

Here is the call graph for this function:

define()

template<typename MOMENTUM >

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

void olb::momenta::GuoZhaoForcedMomentum< MOMENTUM >::define ( CELL & cell, const U & u )

inline

Definition at line 1153 of file elements.h.

  {
    MOMENTUM().template define<TYPE>(cell, u);
  }

getName()

template<typename MOMENTUM >

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

inlinestatic

Definition at line 1182 of file elements.h.

                             {
    return "GuoZhaoForcedMomentum<" + MOMENTUM().getName() + ">";
  }

initialize()

template<typename MOMENTUM >

template<typename TYPE , typename CELL >

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

inline

Definition at line 1159 of file elements.h.

  {
    MOMENTUM().template initialize<TYPE>(cell);
  }

inverseShift()

template<typename MOMENTUM >

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

void olb::momenta::GuoZhaoForcedMomentum< MOMENTUM >::inverseShift ( CELL & cell, U & u )

inline

Definition at line 1165 of file elements.h.

                                                   {
    const V epsilon = cell.template getField<descriptors::EPSILON>();
    const V nu      = cell.template getField<descriptors::NU>();
    const V k       = cell.template getField<descriptors::K>();
    const auto bodyF = cell.template getFieldPointer<descriptors::BODY_FORCE>();
 
    const V uMag = util::sqrt( util::normSqr<V,DESCRIPTOR::d>(u) );
    const V Fe = 0.;//1.75/util::sqrt(150.*util::pow(epsilon,3));
    const V c_0 = 0.5*(1 + 0.5*epsilon*nu/k);
    const V c_1 = 0.5*epsilon*Fe/util::sqrt(k);
 
    for (int iVel=0; iVel < DESCRIPTOR::d; ++iVel) {
      u[iVel] *= (c_0 + util::sqrt(c_0*c_0 + c_1*uMag));
      u[iVel] -= bodyF[iVel] * V(0.5) * epsilon;
    }
  }

References olb::util::sqrt().

Here is the call graph for this function:

---

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

• src/dynamics/momenta/elements.h