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
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 10 months ago by Jendrik Weise.