25#ifndef LATTICE_PHYS_SHEAR_RATE_MAG_3D_HH
26#define LATTICE_PHYS_SHEAR_RATE_MAG_3D_HH
44template<
typename T,
typename DESCRIPTOR>
49 this->
getName() =
"physShearRateMag";
50 const int maxC = this->
_sLattice.getLoadBalancer().size();
52 for (
int iC = 0; iC < maxC; iC++) {
61template<
typename T,
typename DESCRIPTOR>
67 this->
getName() =
"shearRateMag";
70template<
typename T,
typename DESCRIPTOR>
76 T omega = DESCRIPTOR::template provides<descriptors::OMEGA>()
77 ? cell.template getField<descriptors::OMEGA>()
78 : 1. / this->_converter.getLatticeRelaxationTime();
80 T pre2 =
util::pow(descriptors::invCs2<T,DESCRIPTOR>()/2.* omega/rho,2.);
82 if constexpr (DESCRIPTOR::template provides<descriptors::FORCE>()) {
83 const auto force = cell.template getField<descriptors::FORCE>();
89 output[0] = gamma / this->_converter.getConversionFactorTime();
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
functor returns pointwise phys shear rate magnitude on local lattice
bool operator()(T output[], const int input[])
has to be implemented for 'every' derived class
BlockLatticePhysShearRateMag3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter)
Platform-abstracted block lattice for external access and inter-block interaction.
Highest-level interface to Cell data.
void computeAllMomenta(T &rho, T u[descriptors::d< DESCRIPTOR >()], T pi[util::TensorVal< DESCRIPTOR >::n]) const
Compute all momenta on the celll, up to second order.
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
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
SuperLatticePhysShearRateMag3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter)
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 > sqrt(cpu::simd::Pack< T > value)
cpu::simd::Pack< T > pow(cpu::simd::Pack< T > base, cpu::simd::Pack< T > exp)
Top level namespace for all of OpenLB.
Collection of common computations for LBM.
Compute number of elements of a symmetric d-dimensional tensor.
Representation of a parallel 2D geometry – header file.