Evaluation of Third-Order Weighted Essentially Non-Oscillatory Scheme Within Implicit Large Eddy Simulation Framework Using OpenFOAM

The current study investigates the performance of implicit Large Eddy Simulation (iLES) incorporating an unstructured third-order Weighted Essentially Non-Oscillatory (WENO) reconstruction method, alongside conventional Large Eddy Simulation (LES) using the Wall-Adapting Local Eddy Viscosity (WALE) model, for wall-bounded flows. Specifically, iLES is applied to the flow around a NACA0012 airfoil at a Reynolds number which involves key flow phenomena such as laminar separation, transition to turbulence, and flow reattachment. Simulations are conducted using the open-source computational fluid dynamics package OpenFOAM, with a second-order implicit Euler scheme for time integration and the Pressure-Implicit Splitting Operator (PISO) algorithm for pressure–velocity coupling. The results are compared against direct numerical simulation (DNS) for the same flow conditions. Key metrics, including the pressure coefficient and reattached turbulent velocity profiles, show excellent agreement between the iLES and DNS reference results. However, both iLES and LES predict a thinner separation bubble in the transitional flow region then DNS. Notably, the iLES approach achieved a 35% reduction in mesh resolution relative to wall-resolving LES, and a 70% reduction relative to DNS, while maintaining satisfactory accuracy. The study also captures detailed instantaneous flow evolution on the airfoil’s upper surface, with evidence suggesting that three-dimensional disturbances arise from interactions between separating boundary layers near the trailing edge.