A three‐dimensional numerical investigation of the idealized planetary boundary layer

Abstract

The nonlinear equations of motion are integrated numerically in time for a region of x‐y‐z space of volume 3h × h × h, where h turns out to be a height slightly above the level where the wind first attains the geostrophic flow direction. Only the ideal case is treated of a horizontal lower boundary, neutral stability, horizontal homogeneity of all dependent mean variables except the mean pressure, and statistically steady state. The resulting flow patterns are turbulent and the eddies transport required amounts of momentum vertically. Topics which are investigated include the relative directions of stress, wind shear and wind; differences in Ekman wind spirals for the neutral numerical case and a stable atmospheric case; profiles of dimensionless turbulence statistics; effect of allowing the mean density to be either constant or to decrease with height; effect of the wind direction or latitude upon the turbulence intensities; and characteristic structure of the eddies in the planetary boundary layer.