25#ifndef LATTICE_PHYS_EFFECTIVE_DISSIPATION_3D_HH
26#define LATTICE_PHYS_EFFECTIVE_DISSIPATION_3D_HH
44template<
typename T,
typename DESCRIPTOR>
50 this->
getName() =
"physEffectiveDissipation";
51 int maxC = this->
_sLattice.getLoadBalancer().size();
53 for (
int iC = 0; iC < maxC; iC++) {
55 this->_converter, smagoConst, effectiveOmegaF));
59template<
typename T,
typename DESCRIPTOR>
64 _converter(converter), _smagoConst(smagoConst), _effectiveOmegaF(effectiveOmegaF)
66 this->
getName() =
"physEffectiveDissipation";
69template<
typename T,
typename DESCRIPTOR>
73 this->_blockLattice.get(input[0], input[1], input[2]).computeAllMomenta(rho,
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 dt = 1./ _converter.getConversionFactorTime();
84 auto cell = this->_blockLattice.get(input[0], input[1], input[2]);
85 T omegaEff = _effectiveOmegaF(cell);
86 T nuEff = ((1./omegaEff)-0.5)/descriptors::invCs2<T,DESCRIPTOR>();
88 output[0] = PiNeqNormSqr * nuEff
89 *
util::pow(omegaEff * descriptors::invCs2<T,DESCRIPTOR>() / rho, 2) / 2.
90 * _converter.getPhysViscosity() / _converter.getLatticeViscosity() / dt / dt;
represents all functors that operate on a DESCRIPTOR in general, e.g. getVelocity(),...
functor returns pointwise turbulent dissipation density on local lattices
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysEffectiveDissipation3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter, T smagoConst, std::function< T(Cell< T, DESCRIPTOR > &)> effectiveOmegaF)
Platform-abstracted block lattice for external access and inter-block interaction.
Highest-level interface to Cell data.
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
SuperLatticePhysEffectiveDissipation3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter, T smagoConst, std::function< T(Cell< T, DESCRIPTOR > &)> effectiveOmegaF)
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.