olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY > Struct Template Reference

Compute and update cell-wise OMEGA using Oswald-de-waele model. More...

#include <powerLawBGKdynamics.h>

Collaboration diagram for olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY >:

Classes
struct	type

Public Types
using	parameters

Static Public Member Functions
static std::string	getName ()
template<typename CELL , typename PARAMETERS , typename OMEGA , typename RHO , typename PI , typename V = typename CELL::value_t>
static V	computeOmega (CELL &cell, PARAMETERS &parameters, OMEGA &omega0, RHO &rho, PI &pi) any_platform

Detailed Description

template<typename COLLISION, bool HERSCHELBULKLEY = false>
struct olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY >

Compute and update cell-wise OMEGA using Oswald-de-waele model.

Definition at line 77 of file powerLawBGKdynamics.h.

Member Typedef Documentation

parameters

template<typename COLLISION , bool HERSCHELBULKLEY = false>

using olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY >::parameters

Initial value:

 typename COLLISION::parameters::template include<
    descriptors::OMEGA, OMEGA_MIN, OMEGA_MAX, M, N, YIELD_STRESS, SHEAR_RATE_MIN
  >

Definition at line 78 of file powerLawBGKdynamics.h.

Member Function Documentation

computeOmega()

template<typename COLLISION , bool HERSCHELBULKLEY = false>

template<typename CELL , typename PARAMETERS , typename OMEGA , typename RHO , typename PI , typename V = typename CELL::value_t>

static V olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY >::computeOmega ( CELL & cell, PARAMETERS & parameters, OMEGA & omega0, RHO & rho, PI & pi )

inlinestatic

Definition at line 88 of file powerLawBGKdynamics.h.

  {
    using DESCRIPTOR = typename CELL::descriptor_t;
    V pre2 = V{0.5} * descriptors::invCs2<V,DESCRIPTOR>() * omega0 / rho; // strain rate tensor prefactor
    pre2 *= pre2;
    V gamma{};
    if constexpr (DESCRIPTOR::template provides<descriptors::SHEAR_RATE_MAGNITUDE>()) {
      gamma = cell.template getField<descriptors::SHEAR_RATE_MAGNITUDE>();
    }
    else {
      if constexpr (DESCRIPTOR::template provides<descriptors::FORCE>()) {
        // Cannot be done in just one line, it gives error - I don't know why. Davide Dapelo
        const auto force = cell.template getField<descriptors::FORCE>();
        gamma = util::sqrt(V{2}*pre2*lbm<DESCRIPTOR>::computePiNeqNormSqr(cell, force));
      }
      else {
        gamma = util::sqrt(V{2}*pre2*lbm<DESCRIPTOR>::computePiNeqNormSqr(cell));
      }
    }
    if constexpr(HERSCHELBULKLEY) {
      gamma = util::max(parameters.template get<SHEAR_RATE_MIN>(), gamma);
    }
    V m = parameters.template get<M>();
    V n = parameters.template get<N>();
    V nuNew = m * util::pow(gamma, n-V{1}); // Ostwald-de Waele relation
    if constexpr(HERSCHELBULKLEY) {
      // Second term necessary for Herschel-Bulkley relation
      nuNew += parameters.template get<YIELD_STRESS>() / gamma;
    }
    V newOmega = V{1} / (nuNew*descriptors::invCs2<V,DESCRIPTOR>() + V{0.5});
    V omegaMax = parameters.template get<OMEGA_MAX>();
    newOmega = util::min(newOmega, omegaMax);
    V omegaMin = parameters.template get<OMEGA_MIN>();
    newOmega = util::max(newOmega, omegaMin);
    return newOmega;
  }

References olb::lbm< DESCRIPTOR >::computePiNeqNormSqr(), olb::descriptors::invCs2(), olb::util::max(), olb::util::min(), olb::util::pow(), and olb::util::sqrt().

Here is the call graph for this function:

Here is the caller graph for this function:

getName()

template<typename COLLISION , bool HERSCHELBULKLEY = false>

static std::string olb::powerlaw::OmegaFromCell< COLLISION, HERSCHELBULKLEY >::getName ( )

inlinestatic

Definition at line 82 of file powerLawBGKdynamics.h.

  {
    return "powerlaw::OmegaFromCell<" + COLLISION::getName() + ">";
  }

---

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

• src/dynamics/powerLawBGKdynamics.h