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

#include <forcing.h>

Collaboration diagram for olb::forcing::LinearVelocity::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)

Detailed Description

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

Definition at line 539 of file forcing.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 541 of file forcing.h.

MomentaF

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

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

Definition at line 540 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::LinearVelocity::combined_collision< DESCRIPTOR, MOMENTA, EQUILIBRIUM, COLLISION >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 544 of file forcing.h.

                                                               {
      V rho, u[DESCRIPTOR::d], pi[util::TensorVal<DESCRIPTOR >::n];
      MomentaF().computeAllMomenta(cell, rho, u, pi);
      auto force = cell.template getFieldPointer<descriptors::FORCE>();
      constexpr int nDim = DESCRIPTOR::d;
      V forceSave[nDim];
      // adds a+Bu to force, where
      //   d=2: a1=v[0], a2=v[1], B11=v[2], B12=v[3], B21=v[4], B22=v[5]
      //   d=2: a1=v[0], a2=v[1], a3=v[2], B11=v[3], B12=v[4], B13=v[5], B21=v[6], B22=v[7], B23=v[8], B31=v[9], B32=v[10], B33=v[11]
      auto v = cell.template getFieldPointer<descriptors::V12>();
      for (int iDim=0; iDim<nDim; ++iDim) {
        forceSave[iDim] = force[iDim];
        force[iDim] += v[iDim];
        for (int jDim=0; jDim<nDim; ++jDim) {
          force[iDim] += v[jDim + iDim*nDim + nDim]*u[jDim];
        }
      }
      for (int iVel=0; iVel<nDim; ++iVel) {
        u[iVel] += force[iVel] / V{2.};
      }
 
      auto statistics = CollisionO().apply(cell, parameters);
      V newOmega = parameters.template get<descriptors::OMEGA>();
      lbm<DESCRIPTOR>::addExternalForce(cell, rho, u, newOmega, force);
      // Writing back to froce fector
      for (int iVel=0; iVel<nDim; ++iVel) {
        force[iVel] = forceSave[iVel];
      }
      return statistics;
    };

References olb::lbm< DESCRIPTOR >::addExternalForce().

Here is the call graph for this function:

---

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

• src/dynamics/forcing.h