Abstract
A theoretical and numerical small‐scale study of the evaporative cooling phenomenon that might appear in the stratocumulus‐topped boundary layers is presented. An ideal configuration of a cloud‐top mixing layer is considered as defined by two non‐turbulent horizontal layers, stably stratified and with buoyancy reversal within a certain range of mixture fractions due to the evaporative cooling. Linear stability analysis of the shear‐free configuration is employed to provide a new interpretation of the buoyancy reversal parameter, namely in terms of a time‐scale ratio between the stable and the unstable modes of the system. An incompressible high‐order numerical algorithm to perform direct numerical simulation of the configuration is described and two‐dimensional simulations of single‐mode perturbations are presented. These simulations confirm the role of the different parameters identified in the linear stability analysis and show that convoluted flow patterns can be generated by the evaporative cooling even for the low levels of buoyancy reversal found in stratocumulus clouds. They also show that there is no enhancement of turbulent entrainment of upper‐layer fluid in the shear‐free configuration, and turbulent mixing enhancement by the evaporative cooling is restricted to the lower layer. Copyright © 2009 Royal Meteorological Society