Abstract
A 2D higher-order turbulence closure model for research on the structure of the thermal internal boundary layer (TIBL) has been developed in this paper. The mean quantities (temperature and wind), as well as their turbulent moments and their distribution under the TIBL, were computed. Results of numerical simulation show that under the initial condition of onshore flow and surface temperature on land being higher, than on water. 1) the profile of the TIBL on shore can be identified by the distributions of the mean wind and temperature, and during the integration hours there is an unstable stratified region over land that extends upward and inland gradually; 2) the shape of the profiles of the TIBL is roughly in concordance with observed profiles, but there are some differences, obviously, between the results computed by the formula of h ∼ x1/2 and the results of the numerical experiment; and 3)
