#include <rtlbmBoundaryDynamics.h>

Inheritance diagram for olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >:

Collaboration diagram for olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >:

Public Types
template<typename M >
using	exchange_momenta = RtlbmDirectedBoundaryDynamics<T,DESCRIPTOR,M,direction,orientation>

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
	RtlbmDirectedBoundaryDynamics (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 directed boundary walls.

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 direction, int orientation>
class olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >

Definition at line 179 of file rtlbmBoundaryDynamics.h.

Member Typedef Documentation

◆ exchange_momenta

template<typename T , typename DESCRIPTOR , typename MOMENTA , int direction, int orientation>

template<typename M >

using olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::exchange_momenta = RtlbmDirectedBoundaryDynamics<T,DESCRIPTOR,M,direction,orientation>

Definition at line 182 of file rtlbmBoundaryDynamics.h.

Constructor & Destructor Documentation

◆ RtlbmDirectedBoundaryDynamics()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int direction, int orientation>

olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::RtlbmDirectedBoundaryDynamics ( T omega_ )

Constructor.

Definition at line 345 of file rtlbmBoundaryDynamics.hh.

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

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 direction, int orientation>

CellStatistic< T > olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::collide ( Cell< T, DESCRIPTOR > & cell )

overridevirtual

Collision step for directed boundary walls.

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

Definition at line 358 of file rtlbmBoundaryDynamics.hh.

{
  typedef DESCRIPTOR L;
  using MomentaF = typename MOMENTA::template type<DESCRIPTOR>;
  T dirichletTemperature = MomentaF().computeRho(cell);
 
  for ( int iPop = 0; iPop < L::q; ++iPop ) {
    cell[iPop] = - descriptors::t<T,L>(iPop);
  }
 
  constexpr auto missingDiagonal = util::subIndexOutgoing<L,direction,orientation>();
  for ( int i : missingDiagonal ) {
    // compute norm of c_iPopMissing
    // is direction axis parallel
    if ( util::normSqr<int>({descriptors::c<L>(i,0), descriptors::c<L>(i,1), descriptors::c<L>(i,2)}) == 1 ) {
      if ( std::is_base_of<DESCRIPTOR,descriptors::D3Q7<> >::value ) {
        cell[i] = (1-descriptors::t<T,L>(0))*dirichletTemperature - descriptors::t<T,L>(i);
      }
      if ( DESCRIPTOR::template provides<descriptors::tag::RTLBM>() ) {
        cell[i] = dirichletTemperature - descriptors::t<T,L>(i);
      }
    }
  }
  if ( std::is_base_of<DESCRIPTOR,descriptors::D3Q7<> >::value ) {
    cell[0] = descriptors::t<T,L>(0)*dirichletTemperature - descriptors::t<T,L>(0);
  }
  return {-1, -1};
}

◆ computeEquilibrium()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int direction, int orientation>

T olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::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 352 of file rtlbmBoundaryDynamics.hh.

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

◆ getOmega()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int direction, int orientation>

T olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::getOmega ( ) const

place holder

Definition at line 388 of file rtlbmBoundaryDynamics.hh.

{
  return T(-1);
}

◆ setOmega()

template<typename T , typename DESCRIPTOR , typename MOMENTA , int direction, int orientation>

void olb::RtlbmDirectedBoundaryDynamics< T, DESCRIPTOR, MOMENTA, direction, orientation >::setOmega ( T omega_ )

place holder

Definition at line 394 of file rtlbmBoundaryDynamics.hh.

395{

396}

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

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

Public Types

Public Member Functions

Detailed Description

Member Typedef Documentation

◆ exchange_momenta

Constructor & Destructor Documentation

◆ RtlbmDirectedBoundaryDynamics()

Member Function Documentation

◆ collide()

◆ computeEquilibrium()

◆ getOmega()

◆ setOmega()