olb::legacy::OffDynamics< T, DESCRIPTOR > Class Template Reference

Dynamics for offLattice boundary conditions OffDynamics are basically NoLatticeDynamics with the additional functionality to store given velocities exactly at boundary links. More...

#include <dynamics.h>

Inheritance diagram for olb::legacy::OffDynamics< T, DESCRIPTOR >:

Collaboration diagram for olb::legacy::OffDynamics< T, DESCRIPTOR >:

Public Member Functions
	OffDynamics (const T _location[DESCRIPTOR::d])
	Constructor.
	OffDynamics (const T _location[DESCRIPTOR::d], T _distances[DESCRIPTOR::q])
	Constructor.
T	computeRho (ConstCell< T, DESCRIPTOR > &cell) const override
	Returns local stored rho which is updated if the bc is used as velocity!=0 condition.
void	computeU (ConstCell< T, DESCRIPTOR > &cell, T u[DESCRIPTOR::d]) const override
	Returns an average of the locally stored u.
void	setBoundaryIntersection (int iPop, T distance)
	Set Intersection of the link and the boundary.
bool	getBoundaryIntersection (int iPop, T intersection[DESCRIPTOR::d])
	Get Intersection of the link and the boundary.
void	defineRho (Cell< T, DESCRIPTOR > &cell, T rho) override
	Set particle density on the cell.
void	defineRho (int iPop, T rho)
	Set single velocity.
void	defineU (Cell< T, DESCRIPTOR > &cell, const T u[DESCRIPTOR::d]) override
	Set fluid velocity on the cell.
void	defineU (const T u[DESCRIPTOR::d])
	Set constant velocity.
void	defineU (int iPop, const T u[DESCRIPTOR::d])
	Set single velocity.
T	getVelocityCoefficient (int iPop)
	Get VelocitySummand for Bouzidi-Boundary Condition.
std::type_index	id () override
	Expose unique type-identifier for RTTI.
Public Member Functions inherited from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >
	NoLatticeDynamics (T rho=T(1))
	You may fix a fictitious density value on no dynamics node via this constructor.
T	computeEquilibrium (int iPop, T rho, const T u[DESCRIPTOR::d]) const override any_platform
	Yields 0;.
CellStatistic< T >	collide (Cell< T, DESCRIPTOR > &cell) override
	Collision step.
void	computeJ (ConstCell< T, DESCRIPTOR > &cell, T j[DESCRIPTOR::d]) const override
	Yields 0;.
void	computeStress (ConstCell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], T pi[util::TensorVal< DESCRIPTOR >::n]) const override
	Yields NaN.
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	defineRhoU (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d]) override
	Does nothing.
void	defineAllMomenta (Cell< T, DESCRIPTOR > &cell, T rho, const T u[DESCRIPTOR::d], const T pi[util::TensorVal< DESCRIPTOR >::n]) override
	Does nothing.
AbstractParameters< T, DESCRIPTOR > &	getParameters (BlockLattice< T, DESCRIPTOR > &block) override
	Parameters access for legacy post processors.
Public Member Functions inherited from olb::Dynamics< T, DESCRIPTOR >
virtual	~Dynamics () any_platform
virtual std::string	getName () const
	Return human-readable name.
virtual void	initialize (Cell< T, DESCRIPTOR > &cell)
	Initialize dynamics-specific data for cell.
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.
virtual void	inverseShiftRhoU (ConstCell< T, DESCRIPTOR > &cell, T &rho, T u[DESCRIPTOR::d]) const
	Calculate population momenta s.t. the physical momenta are reproduced by the computeRhoU.

Additional Inherited Members
Public Types inherited from olb::Dynamics< T, DESCRIPTOR >
using	value_t = T
using	descriptor_t = DESCRIPTOR

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::legacy::OffDynamics< T, DESCRIPTOR >

Dynamics for offLattice boundary conditions OffDynamics are basically NoLatticeDynamics with the additional functionality to store given velocities exactly at boundary links.

Definition at line 145 of file dynamics.h.

Constructor & Destructor Documentation

OffDynamics() [1/2]

template<typename T , typename DESCRIPTOR >

olb::legacy::OffDynamics< T, DESCRIPTOR >::OffDynamics

(

const T

_location[DESCRIPTOR::d]

)

Constructor.

Definition at line 287 of file dynamics.h.

{
  this->getName() = "OffDynamics";
  typedef DESCRIPTOR L;
  for (int iD = 0; iD < L::d; iD++) {
    location[iD] = _location[iD];
  }
  for (int iPop = 0; iPop < L::q; iPop++) {
    distances[iPop] = -1;
    velocityCoefficient[iPop] = 0;
    for (int iD = 0; iD < L::d; iD++) {
      boundaryIntersection[iPop][iD] = _location[iD];
      _u[iPop][iD] = T();
    }
  }
  _rho=T(1);
}

OffDynamics() [2/2]

template<typename T , typename DESCRIPTOR >

olb::legacy::OffDynamics< T, DESCRIPTOR >::OffDynamics	(	const T	_location[DESCRIPTOR::d],
		T	_distances[DESCRIPTOR::q] )

Constructor.

Definition at line 306 of file dynamics.h.

{
  this->getName() = "OffDynamics";
  typedef DESCRIPTOR L;
  for (int iD = 0; iD < L::d; iD++) {
    location[iD] = _location[iD];
  }
  for (int iPop = 0; iPop < L::q; iPop++) {
    distances[iPop] = _distances[iPop];
    velocityCoefficient[iPop] = 0;
    for (int iD = 0; iD < L::d; iD++) {
      boundaryIntersection[iPop][iD] = _location[iD] - _distances[iPop]*descriptors::c<L>(iPop,iD);
      _u[iPop][iD] = T();
    }
  }
  _rho=T(1);
}

Member Function Documentation

computeRho()

template<typename T , typename DESCRIPTOR >

T olb::legacy::OffDynamics< T, DESCRIPTOR >::computeRho ( ConstCell< T, DESCRIPTOR > & cell ) const

overridevirtual

Returns local stored rho which is updated if the bc is used as velocity!=0 condition.

Reimplemented from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >.

Definition at line 325 of file dynamics.h.

{
  /*typedef DESCRIPTOR L;
  T rhoTmp = T();
  T counter = T();
  int counter2 = int();
  for (int iPop = 0; iPop < L::q; iPop++) {
    if (distances[iPop] != -1) {
      rhoTmp += (cell[iPop] + descriptors::t<T,L>(iPop))*descriptors::t<T,L>(iPop);
      counter += descriptors::t<T,L>(iPop);
      counter2++;
    }
  }
  //if (rhoTmp/counter + 1<0.1999) std::cout << rhoTmp/counter2 + 1 <<std::endl;
  //if (rhoTmp/counter + 1>1.001) std::cout << rhoTmp/counter2 + 1 <<std::endl;
  return rhoTmp/counter/counter;*/
  return _rho;
}

computeU()

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::computeU ( ConstCell< T, DESCRIPTOR > & cell, T u[DESCRIPTOR::d] ) const

overridevirtual

Returns an average of the locally stored u.

Reimplemented from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >.

Definition at line 345 of file dynamics.h.

{
  typedef DESCRIPTOR L;
  for (int iD = 0; iD < L::d; iD++) {
    u[iD] = T();
  }
  int counter = 0;
  for (int iPop = 0; iPop < L::q; iPop++) {
    if ( !util::nearZero(distances[iPop]+1) ) {
      for (int iD = 0; iD < L::d; iD++) {
        u[iD] += _u[iPop][iD];
      }
      counter++;
    }
  }
  if (counter!=0) {
    for (int iD = 0; iD < L::d; iD++) {
      u[iD] /= counter;
    }
  }
  return;
}

References olb::util::nearZero().

Here is the call graph for this function:

defineRho() [1/2]

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::defineRho ( Cell< T, DESCRIPTOR > & cell, T rho )

overridevirtual

Set particle density on the cell.

Reimplemented from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >.

Definition at line 394 of file dynamics.h.

{
  _rho=rho;
}

defineRho() [2/2]

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::defineRho	(	int	iPop,
		T	rho )

Set single velocity.

Definition at line 400 of file dynamics.h.

{
  _rho=rho;
}

defineU() [1/3]

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::defineU ( Cell< T, DESCRIPTOR > & cell, const T u[DESCRIPTOR::d] )

overridevirtual

Set fluid velocity on the cell.

Reimplemented from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >.

Definition at line 406 of file dynamics.h.

{
  defineU(u);
}

defineU() [2/3]

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::defineU

(

const T

u[DESCRIPTOR::d]

)

Set constant velocity.

Definition at line 414 of file dynamics.h.

{
  typedef DESCRIPTOR L;
  for (int iPop = 0; iPop < L::q; iPop++) {
    if ( !util::nearZero(distances[iPop]+1) ) {
      defineU(iPop, u);
    }
  }
}

References olb::util::nearZero().

Here is the call graph for this function:

defineU() [3/3]

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::defineU	(	int	iPop,
		const T	u[DESCRIPTOR::d] )

Set single velocity.

Bouzidi velocity boundary condition formulas for the Coefficients:

2* invCs2*weight*(c,u) for dist < 1/2 1/dist*invCs2*weight*(c,u) for dist >= 1/2

Definition at line 430 of file dynamics.h.

{
  OLB_PRECONDITION(distances[iPop] != -1)
  typedef DESCRIPTOR L;
  velocityCoefficient[iPop] = 0;
  // scalar product of c(iPop) and u
  for (int sum = 0; sum < L::d; sum++) { // +/- problem because of first stream than postprocess
    velocityCoefficient[iPop] -= descriptors::c<L>(iPop,sum)*u[sum];
  }
  // compute summand for boundary condition
  velocityCoefficient[iPop] *= 2*descriptors::invCs2<T,L>() * descriptors::t<T,L>(iPop);
 
  for (int iD = 0; iD < L::d; iD++) {
    _u[iPop][iD] = u[iD];
  }
}

References OLB_PRECONDITION.

getBoundaryIntersection()

template<typename T , typename DESCRIPTOR >

bool olb::legacy::OffDynamics< T, DESCRIPTOR >::getBoundaryIntersection	(	int	iPop,
		T	intersection[DESCRIPTOR::d] )

Get Intersection of the link and the boundary.

Definition at line 381 of file dynamics.h.

{
  typedef DESCRIPTOR L;
  if ( !util::nearZero(distances[iPop]+1) ) {
    for (int iD = 0; iD < L::d; iD++) {
      intersection[iD] = boundaryIntersection[iPop][iD];
    }
    return true;
  }
  return false;
}

References olb::util::nearZero().

Here is the call graph for this function:

getVelocityCoefficient()

template<typename T , typename DESCRIPTOR >

T olb::legacy::OffDynamics< T, DESCRIPTOR >::getVelocityCoefficient

(

int

iPop

)

Get VelocitySummand for Bouzidi-Boundary Condition.

Definition at line 449 of file dynamics.h.

{
  return velocityCoefficient[iPop];
}

Here is the caller graph for this function:

id()

template<typename T , typename DESCRIPTOR >

std::type_index olb::legacy::OffDynamics< T, DESCRIPTOR >::id ( )

inlineoverridevirtual

Expose unique type-identifier for RTTI.

Reimplemented from olb::legacy::NoLatticeDynamics< T, DESCRIPTOR >.

Definition at line 172 of file dynamics.h.

                              {
    return typeid(OffDynamics<T,DESCRIPTOR>);
  }

setBoundaryIntersection()

template<typename T , typename DESCRIPTOR >

void olb::legacy::OffDynamics< T, DESCRIPTOR >::setBoundaryIntersection	(	int	iPop,
		T	distance )

Set Intersection of the link and the boundary.

direction points from the fluid node into the solid domain distance is the distance from the fluid node to the solid wall

Definition at line 369 of file dynamics.h.

{
  typedef DESCRIPTOR L;
  distances[iPop] = distance;
  for (int iD = 0; iD < L::d; iD++) {
    boundaryIntersection[iPop][iD] = location[iD] - distance*descriptors::c<L>(iPop,iD);
  }
}

---

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

• src/dynamics/legacy/dynamics.h