Skip to content

Reply To: Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor

Due to recent bot attacks we have chanced the sign-up process. If you want to participate in our forum, first register on this website and then send a message via our contact form.

Forums on OpenLB General Topics Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor Reply To: Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor

#4541
jonathan
Participant

An even easier way of achieving a static solution would be to just specify the field type inside the physFactor calculation.

_physFactor = -this->_blockLattice.get(input[0]+_overlap,input[1]+_overlap).template getFieldPointer<descriptors::OMEGA>()[0] * descriptors::invCs2<T,DESCRIPTOR>() / dt * this->_converter.getPhysDensity() * this->_converter.getPhysViscosity();

The method described in my first answer is more generic and necessary if you want to have easily swappable descriptors.

Best,
Jonathan