Turbulent flow can be numerically resolved with different levels of accuracy. Many numerical approaches for solving turbulence have been proposed, such as the Reynolds-Averaged Navier–Stokes (RANS), the Large Eddy Simulation (LES), and Direct Numerical Simulation (DNS) approaches. Among these numerical methods, the RANS approach, specifically the Eddy Viscosity Model (EVM), is widely used for calculating turbulent flows thanks to its relatively high accuracy in predicting the mean flow features and its more limited computational demands. However, this approach suffers from several weaknesses, e.g., compromised accuracy and uncertainties due to assumptions in the model construction and insufficient incorporation of the fluid physics. In the LES approach, the whole eddy range is separated into two parts, namely, the large-scale eddy and subgrid-scale (SGS) eddy. The former can be directly resolved, while the latter is computed using the SGS model. As the computing power rapidly increases, this approach is extensively used to study turbulence physics and to resolve low-to-medium Reynolds number flows.