25#ifndef LATTICE_PHYS_DISSIPATION_3D_HH
26#define LATTICE_PHYS_DISSIPATION_3D_HH
44template<
typename T,
typename DESCRIPTOR>
49 this->
getName() =
"physDissipation";
50 const int maxC = this->
_sLattice.getLoadBalancer().size();
52 for (
int iC = 0; iC < maxC; iC++) {
61template<
typename T,
typename DESCRIPTOR>
68 this->
getName() =
"physDissipation";
71template<
typename T,
typename DESCRIPTOR>
75 this->_blockLattice.get(input[0], input[1], input[2]).computeAllMomenta(rho, uTemp, pi);
77 T PiNeqNormSqr = pi[0] * pi[0] + 2. * pi[1] * pi[1] + pi[2] * pi[2];
79 PiNeqNormSqr += pi[2] * pi[2] + pi[3] * pi[3] + 2. * pi[4] * pi[4]
83 T nuLattice = this->_converter.getLatticeViscosity();
84 T omega = 1. / this->_converter.getLatticeRelaxationTime();
85 T dt = this->_converter.getConversionFactorTime();
86 output[0] = PiNeqNormSqr * nuLattice
87 *
util::pow(omega * descriptors::invCs2<T,DESCRIPTOR>() / rho, 2) / 2.
88 * this->_converter.getPhysViscosity() / nuLattice / dt / dt;
represents all functors that operate on a DESCRIPTOR in general, e.g. getVelocity(),...
functor returns pointwise dissipation density on local lattices
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysDissipation3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter)
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.
SuperLattice< T, DESCRIPTOR > & _sLattice
SuperLatticePhysDissipation3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
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.