25#ifndef LATTICE_DISSIPATION_3D_HH
26#define LATTICE_DISSIPATION_3D_HH
44template<
typename T,
typename DESCRIPTOR>
49 this->
getName() =
"dissipation";
50 int maxC = this->
_sLattice.getLoadBalancer().size();
52 for (
int iC = 0; iC < maxC; iC++) {
57template<
typename T,
typename DESCRIPTOR>
60 :
BlockLatticeF3D<T, DESCRIPTOR>(blockLattice, 1), _converter(converter)
62 this->
getName() =
"dissipation";
65template<
typename T,
typename DESCRIPTOR>
69 this->_blockLattice.get(input[0], input[1], input[2]).computeAllMomenta(rho,
73 T PiNeqNormSqr = pi[0] * pi[0] + 2. * pi[1] * pi[1] + pi[2] * pi[2];
75 PiNeqNormSqr += pi[2] * pi[2] + pi[3] * pi[3] + 2. * pi[4] * pi[4]
79 T nuLattice = _converter.getLatticeViscosity();
80 T omega = 1. / _converter.getLatticeRelaxationTime();
81 output[0] = PiNeqNormSqr * nuLattice
82 *
util::pow(omega * descriptors::invCs2<T,DESCRIPTOR>(), 2) / rho / 2.;
functor returns pointwise dissipation density on local lattices
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticeDissipation3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR in general, e.g. getVelocity(),...
Platform-abstracted block lattice for external access and inter-block interaction.
std::string & getName()
read and write access to name
std::vector< std::unique_ptr< BlockF3D< T > > > _blockF
Super functors may consist of several BlockF3D<W> derived functors.
SuperLatticeDissipation3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a SuperLattice in general, e.g. getVelocity(),...
SuperLattice< T, DESCRIPTOR > & _sLattice
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
Wrapper functions that simplify the use of MPI.
cpu::simd::Pack< T > pow(cpu::simd::Pack< T > base, cpu::simd::Pack< T > exp)
Top level namespace for all of OpenLB.
Compute number of elements of a symmetric d-dimensional tensor.
Representation of a parallel 2D geometry – header file.