Hi,rnrncheck out our example “”cylinder3d”” and the implementation of the “”drag”” functor. We have two ways for the computation of momentum exchange implemented. The first one is the standard one (google “”momentum exchange LBM””). The second one is the average method which was proposed by Junk some years ago.rnrnBest regardsrnMathias