Reply To: Different Diffusion Coefficients in one Simulation

[Adrian]

A basic introduction into the dynamics system is available in section 2.3 of the user guide.

What you will want to do is to wrap your collision operator in the modifier:

using LocalOmegaBGK = dynamics::Tuple<
  T, DESCRIPTOR,
  momenta::BulkTuple,
  equilibria::SecondOrder,
  collision::ParameterFromCell<descriptors::OMEGA, collision::BGK>
>;

Don’t forget to assign a value to the descriptors::OMEGA field using e.g.

sLattice.defineField<descriptors::OMEGA>(omegaF); // some functor providing the values