Skip to content

Reply To: Not clear which version of Shan Chen model is implemented

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB Bug Reports Not clear which version of Shan Chen model is implemented Reply To: Not clear which version of Shan Chen model is implemented

#7903
Jendrik Weise
Participant

Having investigated this further, point 2 is merely a result of “force” being actually an acceleration and simply confusingly named. Having now tested 1, correct results are indeed obtained by switching to dynamics along the lines of:

template <typename MOMENTA>
struct Identity {
template<typename DESCRIPTOR>
using type = typename MOMENTA::template type<DESCRIPTOR>;
};

template <typename T, typename DESCRIPTOR, typename MOMENTA=momenta::BulkTuple>
using IdentityForcedBGKdynamics = dynamics::Tuple<
T, DESCRIPTOR,
MOMENTA,
equilibria::SecondOrder,
collision::BGK,
forcing::Guo<Identity>
>;
and adding the forcing momentum adjustments within the processor. I believe altering the default behavior here and documenting this in more detail may indeed by advisable.

  • This reply was modified 7 months, 2 weeks ago by Jendrik Weise.