Implementation of turbulence damping in the OpenFOAM multiphase flow solver interFoam

In the presented work Egorov’s approach (adding a source term to the ω -equation in the k - ω model, which mimics the damping of turbulence close to a solid wall) was implemented in on the subclass of shear stress transport models. Hence, turbulence damping is available for all shear stress transport type models, including hybrid models that are based on the ω -equation. It is shown that turbulence damping improves the prediction of the axial velocity profile not only for Reynolds-averaged Navier–Stokes simulation but also for detached eddy simulation and delayed detached eddy simulation models. Furthermore, it leads to a more realistic estimation of the pressure drop and, hence, to a more correct prediction of the liquid level. In this paper, simulation results for four different turbulence models are presented and validated by comparison with experimental data. Furthermore, the influence of the magnitude of the damping factor on the pressure drop in the channel is investigated for a variety of different gas-to-liquid flow rate ratios. These investigations show that higher gas-to-liquid flow rate ratios require higher damping factors to correctly predict the pressure drop. In the end, advice is formulated on how an appropriate damping factor can be determined for a specific test case.