The Cooling Box Problem: Convection with a quadratic equation of state

Abstract: We investigate the convective cooling of a fluid with a quadratic equation of state by performing three-dimensional direct numerical simulations of a flow with a fixed top-boundary temperature, which is lower than the initial fluid temperature. We consider fluid temperatures near the density maximum, where the nonlinearity is expected to be important. When the equation of state is nonlinear, the resultant vertical transport of heat is fundamentally different and significantly lower than the predictions derived for a linear equation of state. Further, three dimensionless groups parameterize the convective system: the Rayleigh number ($\hbox{Ra}0$), the Prandtl number (Pr), and the dimensionless bottom water temperature $(T_B)$. We further define an effective Rayleigh number ($\hbox{Ra}{eff} = Ra_0 \ T_B^2$), which is equivalent to the traditional Rayleigh number used with a linear equation of state. We present a predictive model for the vertical heat flux, the top boundary-layer thickness, and the turbulent kinetic energy of the system. We show that this model agrees well with the direct numerical simulations. This model could be used to understand how quickly freshwater lakes cool in high latitude environments.