olb::guoZhao::guoZhao_lbm< DESCRIPTOR > Struct Template Reference

#include <guoZhaoDynamics.h>

Collaboration diagram for olb::guoZhao::guoZhao_lbm< DESCRIPTOR >:

Static Public Member Functions
template<typename CELL , typename RHO , typename U , typename OMEGA , typename EPSILON , typename K , typename NU , typename FORCE , typename V = typename CELL::value_t>
static void	addExternalForce (CELL &cell, const RHO &rho, const U &u, const OMEGA &omega, const EPSILON &epsilon, const K &k, const NU &nu, const FORCE &force) any_platform
	Add a force term after BGK collision.
template<typename CELL , typename U , typename EPSILON , typename K , typename NU , typename FORCE , typename BODY_F , typename V = typename CELL::value_t>
static void	updateExternalForce (CELL &cell, const U &u, const EPSILON &epsilon, const K &k, const NU &nu, FORCE &force, BODY_F &bodyF) any_platform
	Updates the force terms with the reaction from the porous matrix, before applying it.

Detailed Description

template<typename DESCRIPTOR>
struct olb::guoZhao::guoZhao_lbm< DESCRIPTOR >

Definition at line 75 of file guoZhaoDynamics.h.

Member Function Documentation

addExternalForce()

template<typename DESCRIPTOR >

template<typename CELL , typename RHO , typename U , typename OMEGA , typename EPSILON , typename K , typename NU , typename FORCE , typename V = typename CELL::value_t>

static void olb::guoZhao::guoZhao_lbm< DESCRIPTOR >::addExternalForce ( CELL & cell, const RHO & rho, const U & u, const OMEGA & omega, const EPSILON & epsilon, const K & k, const NU & nu, const FORCE & force )

inlinestatic

Add a force term after BGK collision.

Definition at line 78 of file guoZhaoDynamics.h.

  {
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      V c_u = V();
      for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
        c_u += descriptors::c<DESCRIPTOR>(iPop,iD)*u[iD];
      }
      c_u *= descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()/epsilon;
      V forceTerm = V();
      for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
        forceTerm +=
          (   (epsilon*(V)descriptors::c<DESCRIPTOR>(iPop,iD)-u[iD]) * descriptors::invCs2<V,DESCRIPTOR>()/epsilon
              + c_u * descriptors::c<DESCRIPTOR>(iPop,iD)
          )
          * force[iD];
      }
      forceTerm *= descriptors::t<V,DESCRIPTOR>(iPop);
      forceTerm *= V(1) - omega/V(2);
      forceTerm *= rho;
      cell[iPop] += forceTerm;
    }
  }

Here is the caller graph for this function:

updateExternalForce()

template<typename DESCRIPTOR >

template<typename CELL , typename U , typename EPSILON , typename K , typename NU , typename FORCE , typename BODY_F , typename V = typename CELL::value_t>

static void olb::guoZhao::guoZhao_lbm< DESCRIPTOR >::updateExternalForce ( CELL & cell, const U & u, const EPSILON & epsilon, const K & k, const NU & nu, FORCE & force, BODY_F & bodyF )

inlinestatic

Updates the force terms with the reaction from the porous matrix, before applying it.

Definition at line 103 of file guoZhaoDynamics.h.

  {
    const V uMag = util::sqrt( util::normSqr<V,DESCRIPTOR::d>(u) );
    const V Fe = 0;//1.75/util::sqrt(150.*util::pow(epsilon,3));
 
    // Linear Darcy term, nonlinear Forchheimer term and body force
    for (int iDim=0; iDim <DESCRIPTOR::d; iDim++) {
      force[iDim] = -u[iDim]*epsilon*nu/k - epsilon*Fe/util::sqrt(k)*uMag*u[iDim] + bodyF[iDim]*epsilon;
    }
  }

References olb::util::sqrt().

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• src/dynamics/guoZhaoDynamics.h