The atmospheric boundary layer (ABL) ranges from hundreds of meters to several kilometers depending on meteorological conditions, mainly wind, temperature, and humidity. Thus, structure of ABL is modified by the daily cycle of heating and cooling over Earth's surface producing three canonical types of boundary layers: convective or unstable, neutral, and stable boundary layers. Convective boundary layer is commonly observed during day when the surface is heated by the sun resulting in a positive buoyancy force, while stable boundary layer occurs during night when surface is cooled by radiation producing a negative buoyancy force, and neutral boundary layer is the case between the former two with little or no buoyancy.The structure of ABL has an important effect on anthropic activities such as mesoscale weather forecasting or pollutant dispersion in urban areas (Fernando et al., 2001). To better understand ABL and related urban processes, numerical simulation is a good complement to field measurements and wind tunnel experiments (Blocken, 2015). In the past much attention has been paid to the accurate CFD modeling of the ABL, both using Reynolds-averaged Navier-Stokes and large-eddy simulation (LES) approaches.LES (Sagaut, 2006), which is a high-fidelity approach for the unsteady simulation of turbulent flows, has been successfully applied to simulation of ABL (e.g.,