functor to get the pointwise dual dissipation density on local lattices, if globIC is not on the local processor, the returned vector is empty More...

#include <dualFunctors3D.h>

Inheritance diagram for olb::opti::BlockLatticeDphysDissipationDf< T, DESCRIPTOR >:

Collaboration diagram for olb::opti::BlockLatticeDphysDissipationDf< T, DESCRIPTOR >:

Public Member Functions
	BlockLatticeDphysDissipationDf (BlockLattice< T, DESCRIPTOR > &blockLattice, int overlap, const UnitConverter< T, DESCRIPTOR > &converter)

bool	operator() (T output[], const int input[])

Protected Attributes
const int	_overlap

const UnitConverter< T, DESCRIPTOR > &	_converter

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::opti::BlockLatticeDphysDissipationDf< T, DESCRIPTOR >

functor to get the pointwise dual dissipation density on local lattices, if globIC is not on the local processor, the returned vector is empty

Definition at line 39 of file dualFunctors3D.h.

Constructor & Destructor Documentation

◆ BlockLatticeDphysDissipationDf()

template<typename T , typename DESCRIPTOR >

olb::opti::BlockLatticeDphysDissipationDf< T, DESCRIPTOR >::BlockLatticeDphysDissipationDf	(	BlockLattice< T, DESCRIPTOR > &	blockLattice,
		int	overlap,
		const UnitConverter< T, DESCRIPTOR > &	converter )

Definition at line 40 of file dualFunctors3D.hh.

  : BlockLatticeF<T,DESCRIPTOR>(blockLattice,DESCRIPTOR::q),
    _overlap(overlap),
    _converter(converter)
{
  this->getName() = "dPhysDissipationDf";
}

References olb::getName().

Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename DESCRIPTOR >

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

Definition at line 52 of file dualFunctors3D.hh.

{
  const T omega = _converter.getLatticeRelaxationFrequency();
  const T dt = _converter.getConversionFactorTime();
 
  T rho, u[DESCRIPTOR::d], pi[util::TensorVal<DESCRIPTOR >::n];
  auto cell = this->_blockLattice.get(latticeR);
  cell.computeAllMomenta(rho, u, pi);
 
  for (int jPop=0; jPop < DESCRIPTOR::q; ++jPop) {
    dDdf[jPop] = T();
    T dpidf[util::TensorVal<DESCRIPTOR >::n];
    dualLbMomentaHelpers<DESCRIPTOR>::dPiDf(cell, dpidf, jPop);
    for (int iAlpha=0; iAlpha < DESCRIPTOR::d; ++iAlpha) {
      for (int iBeta=0; iBeta < DESCRIPTOR::d; ++iBeta) {
        const int iPi = util::serialSymmetricTensorIndex<DESCRIPTOR::d>(iAlpha, iBeta);
        dDdf[jPop] += pi[iPi] * dpidf[iPi] - pi[iPi] * pi[iPi] / rho;
      }
    }
    dDdf[jPop] *= T(2)*util::pow(omega*descriptors::invCs2<T,DESCRIPTOR>()/rho,2)/2.*_converter.getPhysViscosity()/dt/dt;
  }
  return true;
}