In the current code, the friction Reynolds number Retau is computed in the statistics module using velocity derivatives obtained through finite differences at the top and bottom of the domain. When IB walls are present, these velocity derivatives should be computed at the IB wall instead of a the boundary of the computational domain. Extra complications arise when the IB walls are non-stationary. We have to find a robust way to compute Retau in these cases, probably with different formulas depending on the value of body_type.
In the current code, the friction Reynolds number
Retauis computed in thestatisticsmodule using velocity derivatives obtained through finite differences at the top and bottom of the domain. When IB walls are present, these velocity derivatives should be computed at the IB wall instead of a the boundary of the computational domain. Extra complications arise when the IB walls are non-stationary. We have to find a robust way to computeRetauin these cases, probably with different formulas depending on the value ofbody_type.