I have a branch of PyFR that does recirculation with body forcing to maintain a constant bulk velocity :
[soln-plugin-inletforcing]
inletname = 'periodic_0_l'
area = 9.01441
mdotstar = 1
'Inletname' is the inlet (or outlet) boundary name in the mesh file (which has to be a periodic boundary condition), 'area' is the total inlet area, and mdotstar is the desired mass flow rate normalized by area.
Best,
Tarik