You're right! Now in a two-phase flow, the operating density change works.
With the "mixture averaged" setting, the reverse flow disappears - while with any other it is almost 100%.
I've read that in Boussinesq's multiphase flow approach, "minimum phase average" should be selected. However, it gives me wrong results.