In the present paper, a hybrid RANS/LES model with the wall-modelled LES capability (called WM-HRL model) is developed to perform the high-fidelity CFD simulation investigation for complex flow phenomena around a bluff body with coherent structure in subcritical Reynolds number region. The proposed method can achieve a fast and seamless transition from RANS to LES through a filter parameter r_k which is only related to the space resolution capability of the local grid system for various turbulent scales. Furthermore, the boundary positions of the transition region from RANS to LES, as well as the resolving capabilities for the turbulent kinetic energy in the three regions, i.e. RANS, LES and transition region, can be preset by two guide index parameters nr_k1-Q and nr_k2-Q. Through a series of numerical simulations of the flow around a circular cylinder at Reynolds number Re = 3900, the combination conditions are obtained for such two guide index parameters nr_k1-Q and nr_k2-Q that have the capability of high-fidelity resolving and capturing temporally- and spatially-developing coherent structures for such complex three-dimensional flows around such a circular cylinder. The results demonstrate that the new WM-HRL model is capable of accurately resolving and capturing the fine spectral structures of the small-scale Kelvin-Helmholtz instability in the shear layer for flow around such a circular cylinder. Furthermore, under a consistent grid system, through different combinations of these two guide index parameters r_k1 and r_k2, the fine structuresof the recirculation zones with two different lengths and the U-shaped and V-shaped distribution of the average stream-wise velocity in the cylinder near the wake can also be obtained.