Large-eddy simulation of flow past a circular cylinder for Reynolds numbers 400 to 3900

Hongyi Jiang,Liang Cheng

DOI: https://doi.org/10.1063/5.0041168

IF: 4.6

2021-03-01

Physics of Fluids

Abstract:The benchmarking case of flow past a circular cylinder at the Reynolds number (<i>Re</i>) of 3900 is computed with two open-source codes, OpenFOAM and Nektar++, which are based on the conventional finite volume method (FVM) and the high-order spectral/<i>hp</i> element method, respectively. By using the Nektar++ model, mesh convergence for the case <i>Re </i>=3900 is demonstrated (perhaps for the first time) through a systematic mesh dependence study, which includes separate examinations of the spanwise domain length (<i>L_z</i>/<i>D</i>), spanwise resolution, and the resolution in the plane perpendicular to the spanwise direction. The computational efficiencies for the Nektar++ and OpenFOAM approaches are then compared. This benchmarking study adds value to the broad Nektar++ and OpenFOAM communities and to the numerical modeling of bluff-body flows in general. Based on the Nektar++ approach, the computations are then generalized to a range of <i>Re </i>=400–3900. It is found that <i>L_z</i>/<i>D </i>=3 is adequate for <i>Re </i>=2500–3900, while an increased <i>L_z</i>/<i>D </i>=6 is recommended for <i>Re </i>=400–2000. Based on the present high-fidelity numerical data, the physical mechanisms for the variations in the wake recirculation length and the hydrodynamic forces and pressure on the cylinder with <i>Re</i> are explored. In particular, the physics behind the inverse correlation between the root mean square lift coefficient (<i>C_L′</i>) and the wake recirculation length, which includes a significant decrease in <i>C_L′</i> over <i>Re </i>=270–1500, is highlighted.

mechanics,physics, fluids & plasmas