OpenLB 1.8.1
Loading...
Searching...
No Matches
olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 > Class Template Referencefinal

#include <advectionDiffusionBoundaries.h>

+ Inheritance diagram for olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >:
+ Collaboration diagram for olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >:

Public Types

using MomentaF = typename MOMENTA::template type<DESCRIPTOR>
 
using EquilibriumF = typename DYNAMICS::template exchange_momenta<MOMENTA>::EquilibriumF
 
using parameters = typename DYNAMICS::parameters
 
template<typename NEW_T >
using exchange_value_type
 
template<typename M >
using exchange_momenta = AdvectionDiffusionEdgesDynamics<T,DESCRIPTOR,DYNAMICS,M,plane,normal1,normal2>
 
- Public Types inherited from olb::dynamics::CustomCollision< T, DESCRIPTOR, MOMENTA >
using value_t = T
 
using descriptor_t = DESCRIPTOR
 
using MomentaF = typename MOMENTA::template type<DESCRIPTOR>
 
- Public Types inherited from olb::Dynamics< T, DESCRIPTOR >
using value_t = T
 
using descriptor_t = DESCRIPTOR
 

Public Member Functions

std::type_index id () override
 Expose unique type-identifier for RTTI.
 
AbstractParameters< T, DESCRIPTOR > & getParameters (BlockLattice< T, DESCRIPTOR > &block) override
 Parameters access for legacy post processors.
 
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V > collide (CELL &cell, PARAMETERS &parameters) any_platform
 
void computeEquilibrium (ConstCell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], T fEq[DESCRIPTOR::q]) const override
 Return iPop equilibrium for given first and second momenta.
 
std::string getName () const override
 Return human-readable name.
 
- Public Member Functions inherited from olb::dynamics::CustomCollision< T, DESCRIPTOR, MOMENTA >
void initialize (Cell< T, DESCRIPTOR > &cell) override
 Initialize dynamics-specific data for cell.
 
computeRho (ConstCell< T, DESCRIPTOR > &cell) const override
 Compute particle density.
 
void computeU (ConstCell< T, DESCRIPTOR > &cell, T u[DESCRIPTOR::d]) const override
 Compute fluid velocity.
 
void computeJ (ConstCell< T, DESCRIPTOR > &cell, T j[DESCRIPTOR::d]) const override
 Compute fluid momentum.
 
void computeStress (ConstCell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], T pi[util::TensorVal< DESCRIPTOR >::n]) const override
 Compute stress tensor.
 
void computeRhoU (ConstCell< T, DESCRIPTOR > &cell, T &rho, T u[DESCRIPTOR::d]) const override
 Compute fluid velocity and particle density.
 
void computeAllMomenta (ConstCell< T, DESCRIPTOR > &cell, T &rho, T u[DESCRIPTOR::d], T pi[util::TensorVal< DESCRIPTOR >::n]) const override
 Compute all momenta up to second order.
 
void defineRho (Cell< T, DESCRIPTOR > &cell, T rho) override
 Set particle density.
 
void defineU (Cell< T, DESCRIPTOR > &cell, const T u[DESCRIPTOR::d]) override
 Set fluid velocity.
 
void defineRhoU (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d]) override
 Define fluid velocity and particle density.
 
void defineAllMomenta (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], const T pi[util::TensorVal< DESCRIPTOR >::n]) override
 Define all momenta up to second order.
 
void inverseShiftRhoU (ConstCell< T, DESCRIPTOR > &cell, T &rho, T u[DESCRIPTOR::d]) const override
 Calculate population momenta s.t. the physical momenta are reproduced by the computeRhoU.
 
- Public Member Functions inherited from olb::Dynamics< T, DESCRIPTOR >
virtual ~Dynamics () any_platform
 
virtual CellStatistic< T > collide (Cell< T, DESCRIPTOR > &cell)
 Perform purely-local collision step on Cell interface (legacy, to be deprecated)
 
void iniEquilibrium (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d])
 Initialize to equilibrium distribution.
 
void iniRegularized (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], const T pi[util::TensorVal< DESCRIPTOR >::n])
 Initialize cell to equilibrium and non-equilibrum part.
 

Detailed Description

template<typename T, typename DESCRIPTOR, typename DYNAMICS, typename MOMENTA, int plane, int normal1, int normal2>
class olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >

Definition at line 125 of file advectionDiffusionBoundaries.h.

Member Typedef Documentation

◆ EquilibriumF

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
using olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::EquilibriumF = typename DYNAMICS::template exchange_momenta<MOMENTA>::EquilibriumF

Definition at line 128 of file advectionDiffusionBoundaries.h.

◆ exchange_momenta

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
template<typename M >
using olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::exchange_momenta = AdvectionDiffusionEdgesDynamics<T,DESCRIPTOR,DYNAMICS,M,plane,normal1,normal2>

Definition at line 141 of file advectionDiffusionBoundaries.h.

◆ exchange_value_type

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
template<typename NEW_T >
using olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::exchange_value_type
Initial value:
AdvectionDiffusionEdgesDynamics<
NEW_T, DESCRIPTOR,
typename DYNAMICS::template exchange_value_type<NEW_T>,
MOMENTA,
plane, normal1, normal2
>
AdvectionDiffusionEdgesDynamics< NEW_T, DESCRIPTOR, typename DYNAMICS::template exchange_value_type< NEW_T >, MOMENTA, plane, normal1, normal2 > exchange_value_type

Definition at line 133 of file advectionDiffusionBoundaries.h.

◆ MomentaF

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
using olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 127 of file advectionDiffusionBoundaries.h.

◆ parameters

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
using olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::parameters = typename DYNAMICS::parameters

Definition at line 130 of file advectionDiffusionBoundaries.h.

Member Function Documentation

◆ collide()

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V > olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::collide ( CELL & cell,
PARAMETERS & parameters )
inline

Definition at line 152 of file advectionDiffusionBoundaries.h.

152 {
153 V temperature = MomentaF().computeRho(cell);
154 auto u = cell.template getField<descriptors::VELOCITY>();
155 // I need to get Missing information on the corners !!!!
157 // here I know all missing and non missing f_i
158
159 // The collision procedure for D2Q5 and D3Q7 lattice is the same ...
160 // Given the rule f_i_neq = -f_opposite(i)_neq
161 // I have the right number of equations for the number of unknowns using these lattices
162
163 for (unsigned iPop : unknownIndexes) {
164 cell[iPop] = equilibrium<DESCRIPTOR>::firstOrder(iPop, temperature, u)
167 temperature,
168 u));
169 }
170
171 // Once all the f_i are known, I can call the collision for the Regularized Model.
172 return typename DYNAMICS::template exchange_momenta<MOMENTA>::CollisionO().apply(cell, parameters);
173 };
typename MOMENTA::template type< DESCRIPTOR > MomentaF
constexpr int opposite(unsigned iPop) any_platform
Definition functions.h:95
constexpr auto subIndexOutgoing3DonEdges() any_platform
Definition util.h:292
static V firstOrder(int iPop, const RHO &rho, const U &u) any_platform
Computation of equilibrium distribution, first order in u.
Definition lbm.h:37

References olb::equilibrium< DESCRIPTOR >::firstOrder(), olb::descriptors::opposite(), and olb::util::subIndexOutgoing3DonEdges().

+ Here is the call graph for this function:

◆ computeEquilibrium()

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
void olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::computeEquilibrium ( ConstCell< T, DESCRIPTOR > & cell,
T rho,
const T u[DESCRIPTOR::d],
T fEq[DESCRIPTOR::q] ) const
inlineoverridevirtual

Return iPop equilibrium for given first and second momenta.

Implements olb::Dynamics< T, DESCRIPTOR >.

Definition at line 175 of file advectionDiffusionBoundaries.h.

175 {
176 for ( int iPop = 0; iPop < DESCRIPTOR::q; iPop++ ) {
177 fEq[iPop] = equilibrium<DESCRIPTOR>::firstOrder(iPop, rho, u);
178 }
179 };

References olb::equilibrium< DESCRIPTOR >::firstOrder().

+ Here is the call graph for this function:

◆ getName()

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
std::string olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::getName ( ) const
inlineoverridevirtual

Return human-readable name.

Reimplemented from olb::Dynamics< T, DESCRIPTOR >.

Definition at line 181 of file advectionDiffusionBoundaries.h.

181 {
182 return "AdvectionDiffusionEdgesDynamics";
183 };

◆ getParameters()

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
AbstractParameters< T, DESCRIPTOR > & olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::getParameters ( BlockLattice< T, DESCRIPTOR > & block)
inlineoverridevirtual

Parameters access for legacy post processors.

Implements olb::Dynamics< T, DESCRIPTOR >.

Definition at line 147 of file advectionDiffusionBoundaries.h.

147 {
148 return block.template getData<OperatorParameters<AdvectionDiffusionEdgesDynamics>>();
149 }

◆ id()

template<typename T , typename DESCRIPTOR , typename DYNAMICS , typename MOMENTA , int plane, int normal1, int normal2>
std::type_index olb::AdvectionDiffusionEdgesDynamics< T, DESCRIPTOR, DYNAMICS, MOMENTA, plane, normal1, normal2 >::id ( )
inlineoverridevirtual

Expose unique type-identifier for RTTI.

Implements olb::Dynamics< T, DESCRIPTOR >.

Definition at line 143 of file advectionDiffusionBoundaries.h.

143 {
144 return typeid(AdvectionDiffusionEdgesDynamics);
145 };

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