olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 > Class Template Reference

Defines incoming directions on edge boundaries. More...

#include <rtlbmBoundaryDynamics.h>

Inheritance diagram for olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >:

Collaboration diagram for olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >:

Public Types
template<typename M >
using	exchange_momenta = RtlbmDiffuseConstEdgeBoundaryDynamics<T,DESCRIPTOR,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
	RtlbmDiffuseConstEdgeBoundaryDynamics (T omega_)
	Constructor.
T	computeEquilibrium (int iPop, T rho, const T u[DESCRIPTOR::d]) const override
	Compute equilibrium distribution function.
CellStatistic< T >	collide (Cell< T, DESCRIPTOR > &cell) override
	Collision step for edges.
T	getOmega () const
	place holder
void	setOmega (T omega_)
	place holder
Public Member Functions inherited from olb::legacy::BasicDynamics< T, DESCRIPTOR, MOMENTA >
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.
Public Member Functions inherited from olb::dynamics::CustomCollision< T, DESCRIPTOR, MOMENTA >
void	initialize (Cell< T, DESCRIPTOR > &cell) override
	Initialize dynamics-specific data for cell.
T	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 std::string	getName () const
	Return human-readable name.
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 MOMENTA, int plane, int normal1, int normal2>
class olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >

Defines incoming directions on edge boundaries.

Emposed dirichlet density is distributed equivalently to all incoming/unkown directions.

Parameters

plane	takes values 0, 1, 2; which denotes the alignment of edge either along x, y or z axis
normal1	takes values -1 or 1; normal orientation of a surface side 'right'
normal2	takes values -1 or 1; normal orientation of a surface side 'left'

Definition at line 134 of file rtlbmBoundaryDynamics.h.

Member Typedef Documentation

exchange_momenta

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

template<typename M >

using olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::exchange_momenta = RtlbmDiffuseConstEdgeBoundaryDynamics<T,DESCRIPTOR,M,plane,normal1,normal2>

Definition at line 137 of file rtlbmBoundaryDynamics.h.

Constructor & Destructor Documentation

RtlbmDiffuseConstEdgeBoundaryDynamics()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::RtlbmDiffuseConstEdgeBoundaryDynamics

(

T

omega_

)

Constructor.

Definition at line 229 of file rtlbmBoundaryDynamics.hh.

  : legacy::BasicDynamics<T,DESCRIPTOR,MOMENTA>()
{
  this->getName() = "RtlbmDiffuseConstEdgeBoundaryDynamics";
}

References olb::Dynamics< T, DESCRIPTOR >::getName().

Here is the call graph for this function:

Member Function Documentation

collide()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

CellStatistic< T > olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::collide ( Cell< T, DESCRIPTOR > & cell )

overridevirtual

Collision step for edges.

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

Definition at line 242 of file rtlbmBoundaryDynamics.hh.

{
  // For direction i \in I_in define
  // cell_i = w_i * dirichlet/sumWeights - w_i
  // For direction i \in I_out defube
  // cell_i = - w_i
  // This construction yields
  // sum_{i=0}^{q-1} cell_i == dirichlet - 1
 
  typedef DESCRIPTOR L;
 
  // shift all: cell_i = f_i - weight_i
  for ( int iPop = 0; iPop < L::q; ++iPop ) {
    cell[iPop] = - descriptors::t<T,L>(iPop);
  }
 
  constexpr auto missing_iPop = util::subIndexOutgoing3DonEdges<L,plane,normal1,normal2>();
  T sumWeights = 0;
  for ( int i : missing_iPop ) {
    sumWeights += descriptors::t<T,L>(i);
  }
 
  using MomentaF = typename MOMENTA::template type<DESCRIPTOR>;
  T dirichletTemperature = MomentaF().computeRho(cell);
  for ( int i : missing_iPop ) {
    cell[i] = descriptors::t<T,L>(i)*dirichletTemperature/sumWeights - descriptors::t<T,L>(i);
  }
  return {-1, -1};
}

computeEquilibrium()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

T olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::computeEquilibrium ( int iPop, T rho, const T u[DESCRIPTOR::d] ) const

overridevirtual

Compute equilibrium distribution function.

Reimplemented from olb::legacy::BasicDynamics< T, DESCRIPTOR, MOMENTA >.

Definition at line 236 of file rtlbmBoundaryDynamics.hh.

{
  return descriptors::t<T,DESCRIPTOR>(iPop)*rho - descriptors::t<T,DESCRIPTOR>(iPop);
}

getOmega()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

T olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::getOmega

(

)

const

place holder

Definition at line 273 of file rtlbmBoundaryDynamics.hh.

{
  return T(-1);
}

setOmega()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int plane, int normal1, int normal2>

void olb::RtlbmDiffuseConstEdgeBoundaryDynamics< T, DESCRIPTOR, MOMENTA, plane, normal1, normal2 >::setOmega

(

T

omega_

)

place holder

Definition at line 279 of file rtlbmBoundaryDynamics.hh.

{
}

---

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

• src/boundary/rtlbmBoundaryDynamics.h
• src/boundary/rtlbmBoundaryDynamics.hh