olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION > Struct Template Reference

#include <forcing.h>

Collaboration diagram for olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >:

Public Types
using	MomentaF = typename MOMENTA::template type<DESCRIPTOR>
using	CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Public Member Functions
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM, typename COLLISION>
struct olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >

Definition at line 401 of file forcing.h.

Member Typedef Documentation

CollisionO

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::CollisionO = typename COLLISION::template type<DESCRIPTOR,MOMENTA,EQUILIBRIUM>

Definition at line 403 of file forcing.h.

MomentaF

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

using olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 402 of file forcing.h.

Member Function Documentation

apply()

template<typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM , typename COLLISION >

template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>

CellStatistic< V > olb::forcing::Wagner::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 409 of file forcing.h.

                                                                            {
      V rho, u[DESCRIPTOR::d];
      MomentaF().computeRhoU(cell, rho, u);
      CollisionO().apply(cell, parameters);
      const V omega = cell.template getField<descriptors::OMEGA>();
      const auto force = cell.template getField<descriptors::FORCE>();
      const V d2_rho = cell.template getField<descriptors::SCALAR>();
      Vector<V,DESCRIPTOR::d * DESCRIPTOR::d> psi_w{};
 
      // Forcing scheme
      for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
        V c_u{};
        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>();
        V forceTerm{};
        for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
          forceTerm +=
            (   (descriptors::c<DESCRIPTOR>(iPop,iD) - u[iD]) * descriptors::invCs2<V,DESCRIPTOR>()
                + c_u * descriptors::c<DESCRIPTOR>(iPop,iD)
            ) * force[iD];
        }
 
        V delta_ij = 0;
        for (int n = 0; n < DESCRIPTOR::d; ++n) {
          for (int m = 0; m < DESCRIPTOR::d; ++m) {
            int ind = n*DESCRIPTOR::d + m;
 
            delta_ij = (n == m) ? 1 : 0;
 
            psi_w[ind] = (V(1.)-omega/V(4.))*force[n]*force[m] + omega/V(12.)*( d2_rho )/rho*delta_ij;
 
            // Computing the source term
            forceTerm += V(9.) * ( descriptors::c<DESCRIPTOR>(iPop)[n]
                * descriptors::c<DESCRIPTOR>(iPop)[m] - V(1./3.)*delta_ij ) / V(2.)
                * psi_w[ind];
          }
        }
 
        forceTerm *= descriptors::t<V,DESCRIPTOR>(iPop);
        forceTerm *= rho;
        cell[iPop] += forceTerm;
      }
 
      return {rho, util::normSqr<V,DESCRIPTOR::d>(u)};
    };

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

Here is the call graph for this function:

---

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

• src/dynamics/forcing.h