I am interested in the functions provided in superLatticeIntegralF2D, superLatticeLocalF2d for postprocessing, e.g., SuperLatticePhysWallShearStress2D, SuperLatticePhysDrag2D, etc.
I am also using DynOmegaD2Q9Descriptor for power-law fluid modeling. I am wondering whether those postprocessing functions (SuperLatticePhysWallShearStress2D and SuperLatticePhysDrag2D) can work with DynOmegaD2Q9Descriptor?
It seems that the postprocessing functions use a constant omega from the unit converter but not the varying omega from the external field.
BlockLatticePhysWallShearStress2D uses a constant relaxation frequency only in _physFactor for unit conversion. Its calculation can just be moved into the operator() for a cell dependant calculation.
AFAIK your descriptor DynOmegaD2Q9 stores omega in the external field and can therefore be accessed by a template variable.
Try changing the source of both Super- and BlockLatticePhysWallShearStress to provide a FIELD template argument. For a generic example see how its done in Super- and BlockLatticeField2D.
Then just replace omega in