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

functor returns pointwise phys force for PSM dynamics More...

#include <latticePSMPhysForce3D.h>

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

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

Public Member Functions
	BlockLatticePSMPhysForce3D (BlockLattice< T, DESCRIPTOR > &blockLattice, 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::BlockLatticePhysF3D< T, DESCRIPTOR >
	BlockLatticePhysF3D (BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter, int targetDim)
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::BlockLatticePhysF3D< T, DESCRIPTOR >
const UnitConverter< T, DESCRIPTOR > &	_converter
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::BlockLatticePSMPhysForce3D< T, DESCRIPTOR >

functor returns pointwise phys force for PSM dynamics

Definition at line 59 of file latticePSMPhysForce3D.h.

Constructor & Destructor Documentation

BlockLatticePSMPhysForce3D()

template<typename T , typename DESCRIPTOR >

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

Definition at line 60 of file latticePSMPhysForce3D.hh.

  : BlockLatticePhysF3D<T,DESCRIPTOR>(blockLattice, converter, 3)
{
  this->getName() = "physPSMForce";
}

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::BlockLatticePSMPhysForce3D< 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 69 of file latticePSMPhysForce3D.hh.

{
  for (int i = 0; i < this->getTargetDim(); ++i) {
    output[i] = T();
  }
 
  T epsilon = 1. - this->_blockLattice.get(input).template getField<descriptors::POROSITY>();
 
  //if ((epsilon > 1e-5 && epsilon < 1 - 1e-5)) {
  if ((epsilon > 1e-5)) {
    T rho, u[DESCRIPTOR::d], u_s[DESCRIPTOR::d];
 
    for (int i = 0; i < DESCRIPTOR::d; i++) {
      u_s[i] = this->_blockLattice.get(input).template getFieldComponent<descriptors::VELOCITY_SOLID>(i);
    }
    T paramA = this->_converter.getLatticeRelaxationTime() - 0.5;
    // speed up paramB
    T paramB = (epsilon * paramA) / ((1. - epsilon) + paramA);
 
    T omega_s;
    T omega = 1. / this->_converter.getLatticeRelaxationTime();
 
    this->_blockLattice.get(input).computeRhoU(rho, u);
 
    const T uSqr_s = util::normSqr<T,DESCRIPTOR::d>(u_s);
    T uSqr = util::normSqr<T,DESCRIPTOR::d>(u);
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      //switch (mode) {
      //case M2:
        omega_s = (lbm< DESCRIPTOR>::equilibrium(iPop, rho, u_s, uSqr_s)
                   - this->_blockLattice.get(input[0], input[1], input[2])[iPop])
                  + (1 - omega)
                  * (this->_blockLattice.get(input[0], input[1], input[2])[iPop]
                     - lbm< DESCRIPTOR>::equilibrium(iPop, rho, u, uSqr));
      //  break;
      //case M3:
      //  omega_s =
      //    (this->_blockLattice.get(input[0], input[1], input[2])[descriptors::opposite<DESCRIPTOR>(iPop)]
      //     - lbm< DESCRIPTOR>::equilibrium(
      //       descriptors::opposite<DESCRIPTOR>(iPop), rho, u_s, uSqr_s))
      //    - (this->_blockLattice.get(input[0], input[1], input[2])[iPop]
      //       - lbm< DESCRIPTOR>::equilibrium(iPop, rho, u_s, uSqr_s));
      //}
 
      for (int i = 0; i < this->getTargetDim(); ++i) {
        output[i] -= descriptors::c<DESCRIPTOR>(iPop,i) * omega_s;
      }
    }
 
    for (int i = 0; i < this->getTargetDim(); ++i) {
      output[i] = this->_converter.getPhysForce(output[i] * paramB);
    }
  }
  return true;
}

---

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

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