Large eddy simulation (LES) is widely used to investigate the aerodynamics and convective heat transfer (CHT) at the surfaces of sharp-edged bluff bodies for a wide range of Reynolds (Re) numbers. Due to the heavy computational costs associated with implicit filtering in LES at high Reynolds number flows (Re β₯ 105), wall-modeled (WM) rather than wall-resolved (WR) LES is often adopted. However, the performance of LES-WM for such applications has not yet been systematically investigated. Therefore, this study evaluates the performance of LES-WM and LES-WR for the flow and thermal field at the facades of a low-rise building immersed in an atmospheric boundary layer. Four grids are constructed for LES-WM, each employing different resolution at the building surfaces reaching maximum non-dimensional wall distance y+ = 43, 57, 70, and 95. In addition, the performance of two wall functions, namely the Werner and Wengle and the enhanced wall function is investigated. The results show that the use of LES-WM can result in significant deviations in the predicted near-facade flow pattern and the surface convective heat transfer coefficient (CHTC). Grid resolution significantly impacts the CHTC results and deviations go up to 88% (at the base of the windward facade). Considerable deviations among the employed wall functions are apparent only on the finest grid. In this case, the implementation of the enhanced wall function indicates better performance compared to the non-blended law of the wall (combined with the Werner and Wengle) for CHTC in the regions of the leeward facade where the flow remains attached to the wall. The deviation of the enhanced wall function for surface-averaged CHTC is found to be 10.8% against the wall-resolved LES results, while for the non-blended law of the wall this is 19.2%.</p
