olb::momenta::GuoZhaoForcedStress< STRESS > Struct Template Reference

#include <elements.h>

Collaboration diagram for olb::momenta::GuoZhaoForcedStress< STRESS >:

Public Member Functions
template<typename TYPE , typename CELL , typename RHO , typename U , typename PI , typename V = typename CELL::value_t, typename DESCRIPTOR = typename CELL::descriptor_t>
void	compute (CELL &cell, const RHO &rho, const U &u, PI &pi) any_platform

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

Detailed Description

template<typename STRESS>
struct olb::momenta::GuoZhaoForcedStress< STRESS >

Definition at line 1401 of file elements.h.

Member Function Documentation

compute()

template<typename STRESS >

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

void olb::momenta::GuoZhaoForcedStress< STRESS >::compute ( CELL & cell, const RHO & rho, const U & u, PI & pi )

inline

Definition at line 1403 of file elements.h.

  {
    V uNew[DESCRIPTOR::d] { };
    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 (unsigned iD=0; iD < DESCRIPTOR::d; ++iD) {
      uNew[iD] = u[iD] * (c_0 + util::sqrt(c_0*c_0 + c_1*uMag)) - V{0.5} * bodyF[iD] * epsilon;
    }
    STRESS().template compute<TYPE>(cell, rho, uNew, pi);
    V forceTensor[util::TensorVal<DESCRIPTOR>::n];
    // Creation of body force tensor (rho/2.)*(G_alpha*U_beta + U_alpha*G_Beta)
    int iPi = 0;
    for (int iAlpha=0; iAlpha < DESCRIPTOR::d; ++iAlpha) {
      for (int iBeta=iAlpha; iBeta < DESCRIPTOR::d; ++iBeta) {
        forceTensor[iPi] = V{0.5} * rho * (bodyF[iAlpha]*uNew[iBeta] + uNew[iAlpha]*bodyF[iBeta]);
        ++iPi;
      }
    }
    // Creation of second-order moment off-equilibrium tensor
    for (int iPi=0; iPi < util::TensorVal<DESCRIPTOR>::n; ++iPi) {
      pi[iPi] += forceTensor[iPi];
    }
  }

References olb::util::sqrt().

Here is the call graph for this function:

getName()

template<typename STRESS >

static std::string olb::momenta::GuoZhaoForcedStress< STRESS >::getName ( )

inlinestatic

Definition at line 1435 of file elements.h.

                             {
    return "GuoZhaoForcedStress<" + STRESS().getName() + ">";
  }

---

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

• src/dynamics/momenta/elements.h