25#ifndef LATTICE_PHYS_VISCOSITY_3D_HH
26#define LATTICE_PHYS_VISCOSITY_3D_HH
44template<
typename T,
typename DESCRIPTOR>
49 this->
getName() =
"physViscosity";
50 const int maxC = this->
_sLattice.getLoadBalancer().size();
52 for (
int iC = 0; iC < maxC; iC++) {
56 this->_converter, logscale)
61template<
typename T,
typename DESCRIPTOR>
68 this->
getName() =
"physViscosity";
71template<
typename T,
typename DESCRIPTOR>
75 const auto omegaPL = cell.template getField<descriptors::OMEGA>();
77 output[0] = std::numeric_limits<double>::quiet_NaN();
80 output[0] = (1./omegaPL - 0.5) / descriptors::invCs2<T,DESCRIPTOR>()
81 * this->_converter.getPhysDeltaX() * this->_converter.getPhysDeltaX() / this->_converter.getPhysDeltaT();
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
functor returns pointwise phys viscosity on local lattices
BlockLatticePhysViscosity3D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter, bool logscale)
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
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
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
SuperLatticePhysViscosity3D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter, bool logscale=false)
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.
ADf< T, DIM > log10(const ADf< T, DIM > &a)
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.