olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR > Class Template Reference

functor to scale particle distributions to a time step More...

#include <latticeTimeStepScale3D.h>

Inheritance diagram for olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR >:

Collaboration diagram for olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR >:

Public Member Functions
	BlockLatticeTimeStepScale3D (BlockLattice< T, DESCRIPTOR > &blockLattice, T oldTau, const UnitConverter< T, DESCRIPTOR > &converter)
bool	operator() (T output[], const int input[]) override
	has to be implemented for 'every' derived class
Public Member Functions inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > &	getBlock ()
	Copy Constructor.
Public Member Functions inherited from olb::BlockF3D< T >
	~BlockF3D () override
	virtual destructor for defined behaviour
virtual BlockStructureD< 3 > &	getBlockStructure () const
BlockF3D< T > &	operator- (BlockF3D< T > &rhs)
BlockF3D< T > &	operator+ (BlockF3D< T > &rhs)
BlockF3D< T > &	operator* (BlockF3D< T > &rhs)
BlockF3D< T > &	operator/ (BlockF3D< T > &rhs)
Public Member Functions inherited from olb::GenericF< T, int >
virtual	~GenericF ()=default
int	getSourceDim () const
	read only access to member variable _m
int	getTargetDim () const
	read only access to member variable _n
std::string &	getName ()
	read and write access to name
std::string const &	getName () const
	read only access to name
bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above
bool	operator() (T output[], int input0)
bool	operator() (T output[], int input0, int input1)
bool	operator() (T output[], int input0, int input1, int input2)
bool	operator() (T output[], int input0, int input1, int input2, int input3)

Additional Inherited Members
Public Types inherited from olb::GenericF< T, int >
using	targetType
using	sourceType
Public Attributes inherited from olb::GenericF< T, int >
std::shared_ptr< GenericF< T, int > >	_ptrCalcC
	memory management, frees resouces (calcClass)
Protected Member Functions inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
	BlockLatticeF3D (BlockLattice< T, DESCRIPTOR > &blockLattice, int targetDim)
Protected Member Functions inherited from olb::BlockF3D< T >
	BlockF3D (BlockStructureD< 3 > &blockStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< T, int >
	GenericF (int targetDim, int sourceDim)
Protected Attributes inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > &	_blockLattice
Protected Attributes inherited from olb::BlockF3D< T >
BlockStructureD< 3 > &	_blockStructure

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR >

functor to scale particle distributions to a time step

Definition at line 59 of file latticeTimeStepScale3D.h.

Constructor & Destructor Documentation

BlockLatticeTimeStepScale3D()

template<typename T , typename DESCRIPTOR >

olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR >::BlockLatticeTimeStepScale3D	(	BlockLattice< T, DESCRIPTOR > &	blockLattice,
		T	oldTau,
		const UnitConverter< T, DESCRIPTOR > &	converter )

Definition at line 64 of file latticeTimeStepScale3D.hh.

  : BlockLatticeF3D<T,DESCRIPTOR>(blockLattice, DESCRIPTOR::q), _tau_old(oldTau), _converter(converter)
{
  this->getName() = "latticeTimeStepScale";
}

References olb::GenericF< T, int >::getName().

Here is the call graph for this function:

Member Function Documentation

operator()()

template<typename T , typename DESCRIPTOR >

bool olb::BlockLatticeTimeStepScale3D< T, DESCRIPTOR >::operator() ( T output[], const int input[] )

overridevirtual

has to be implemented for 'every' derived class

Implements olb::GenericF< T, int >.

Definition at line 72 of file latticeTimeStepScale3D.hh.

{
 
  Cell<T,DESCRIPTOR> cell = this->_blockLattice.get(input[0], input[1], input[2]);
  Dynamics<T,DESCRIPTOR>* dynamics = this->_blockLattice.getDynamics(input[0], input[1], input[2]);
  T rho_old, rho_new, u_old[DESCRIPTOR::d], u_new[DESCRIPTOR::d];;
  cell.computeRhoU(rho_old,u_old);
  T fNeq[DESCRIPTOR::q];
 
  T uSqr_old = util::normSqr<T,DESCRIPTOR::d>(u_old);
 
  T tau_new = _converter.getLatticeRelaxationTime();
 
  T physDeltaT_new = _converter.getPhysDeltaT();
  T physDeltaT_old = physDeltaT_new * (_tau_old-0.5) / (tau_new-0.5);
  T neqScaling = ((physDeltaT_new*tau_new)/(physDeltaT_old*_tau_old));
 
  for (int iDim=0; iDim < DESCRIPTOR::d; iDim++) {
    u_new[iDim] = u_old[iDim] * physDeltaT_new / physDeltaT_old ;
  }
 
  rho_new = (rho_old-1.0) * (physDeltaT_new / physDeltaT_old) * (physDeltaT_new / physDeltaT_old) + 1.0;
 
  T uSqr_new = util::normSqr<T,DESCRIPTOR::d>(u_new);
 
 
  for (int iPop = 0; iPop < DESCRIPTOR::q; ++iPop) {
    fNeq[iPop] = cell[iPop] - dynamics -> computeEquilibrium(iPop,rho_old,u_old,uSqr_old);
//    output[iPop] = cell[iPop];
    output[iPop] = (dynamics -> computeEquilibrium(iPop,rho_old,u_old,uSqr_old)+descriptors::t<T,DESCRIPTOR>(iPop)// * ((dynamics -> computeEquilibrium(iPop,rho_old,u_new,uSqr_new)+descriptors::t<T,DESCRIPTOR>(iPop))/(dynamics -> computeEquilibrium(iPop,rho_old,u_old,uSqr_old)+descriptors::t<T,DESCRIPTOR>(iPop)))
                    + fNeq[iPop]*neqScaling)
                   * ((dynamics -> computeEquilibrium(iPop,rho_new,u_new,uSqr_new)+descriptors::t<T,DESCRIPTOR>(iPop))/(dynamics -> computeEquilibrium(iPop,rho_old,u_old,uSqr_old)+descriptors::t<T,DESCRIPTOR>(iPop))) - descriptors::t<T,DESCRIPTOR>(iPop);
  }
  return true;
}

References olb::ConstCell< T, DESCRIPTOR >::computeRhoU().

Here is the call graph for this function:

---

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

• src/functors/lattice/latticeTimeStepScale3D.h
• src/functors/lattice/latticeTimeStepScale3D.hh