24#ifndef LATTICE_PHYS_DISSIPATION_2D_HH
25#define LATTICE_PHYS_DISSIPATION_2D_HH
44template<
typename T,
typename DESCRIPTOR>
49 this->
getName() =
"physDissipation";
50 int maxC = this->
_sLattice.getLoadBalancer().size();
52 for (
int iC = 0; iC < maxC; iC++) {
57template <
typename T,
typename DESCRIPTOR>
62 this->
getName() =
"physDissipation";
65template <
typename T,
typename DESCRIPTOR>
69 this->_blockLattice.get( input[0], input[1] ).computeAllMomenta(rho, uTemp, pi);
71 T PiNeqNormSqr = pi[0]*pi[0] + 2.*pi[1]*pi[1] + pi[2]*pi[2];
73 PiNeqNormSqr += pi[2]*pi[2] + pi[3]*pi[3] + 2.*pi[4]*pi[4] + pi[5]*pi[5];
76 T nuLattice = this->_converter.getLatticeViscosity();
77 T omega = 1. / this->_converter.getLatticeRelaxationTime();
78 T dt = this->_converter.getConversionFactorTime();
79 output[0] = PiNeqNormSqr*nuLattice*
util::pow(omega*descriptors::invCs2<T,DESCRIPTOR>()/rho,2)/2.*this->_converter.getPhysViscosity()/this->_converter.getLatticeViscosity()/dt/dt;
BlockLatticePhysDissipation2D returns pointwise physical dissipation density on local lattices.
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysDissipation2D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
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< BlockF2D< T > > > _blockF
Super functors may consist of several BlockF2D<W> derived functors.
SuperLattice< T, DESCRIPTOR > & _sLattice
SuperLatticePhysDissipation2D(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.
The description of a generic interface for all functor classes – header file.
This file contains indicator functions.
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.