2 results
Structure of subfilter-scale fluxes in the atmospheric surface layer with application to large-eddy simulation modelling
- PETER P. SULLIVAN, THOMAS W. HORST, DONALD H. LENSCHOW, CHIN-HOH MOENG, JEFFREY C. WEIL
-
- Journal:
- Journal of Fluid Mechanics / Volume 482 / 10 May 2003
- Published online by Cambridge University Press:
- 13 May 2003, pp. 101-139
-
- Article
- Export citation
-
In the atmospheric surface layer, the wavelength of the peak in the vertical velocity spectrum $\Lambda_w$ decreases with increasing stable stratification and proximity to the surface and this dependence constrains our ability to perform high-Reynolds-number large-eddy simulation (LES). Near the ground, the LES filter cutoff $\Delta_f$ is comparable to or larger than $\Lambda_w$ and as a result the subfilter-scale (SFS) fluxes in LES are always significant and their contribution to the total flux grows with increasing stability.
We use the three-dimensional turbulence data collected during the Horizontal Array Turbulence Study (HATS) field program to construct SFS fluxes and variances that are modelled in LES codes. Detailed analysis of the measured SFS motions shows that the ratio $\Lambda_w/\Delta_f$ contains the essential information about stratification, vertical distance above the surface, and filter size, and this ratio allows us to connect measurements of SFS variables with LES applications. We find that the SFS fluxes and variances collapse reasonably well for atmospheric conditions and filter widths in the range $\Lambda_w/\Delta_f = [0.2,15]$. The SFS variances are anisotropic and the SFS energy is non-inertial, exhibiting a strong dependence on the stratification, large-scale shear, and proximity to the surface. SFS flux decomposition into modified-Leonard, cross-, and Reynolds terms illustrates that these terms are of comparable magnitude and scale content at large $\Lambda_w/\Delta_f$. As $\Lambda_w/\Delta_f \rightarrow 0$, the SFS flux approaches the-ensemble-average flux and is dominated by the Reynolds term. Backscatter of energy from the SFS motions to the resolved fields is small in the bulk of the surface layer, less than 20% for $\Lambda_w/\Delta_{f} < 2$.
A priori testing of typical SFS models using the HATS dataset shows that the turbulent kinetic energy and Smagorinsky model coefficients $C_k$ and $C_s$ depend on $\Lambda_w/\Delta_f$ and are smaller than theoretical estimates based on the assumption of a sharp spectral cutoff filter in the inertial range. $C_k$ and $C_s$ approach zero for small $\Lambda_w/\Delta_f$. Much higher correlations between measured and modelled SFS fluxes are obtained with a mixed SFS model that explicitly includes the modified-Leonard term. The eddy-viscosity model coefficients still retain a significant dependence on $\Lambda_w/\Delta_f$ with the mixed model. A dissipation model of the form $\epsilon = C_\epsilon E_s^{3/2}/\Delta_f$ is not universal across the range of $\Lambda_w/\Delta_f$ typical of atmospheric LES applications. The inclusion of a shear-stability-dependent length scale (Canuto & Cheng 1997) captures a large fraction of the variation in the eddy-viscosity and dissipation model coefficients.
Laboratory simulation of the rise of buoyant thermals created by open detonation
- ROGER S. THOMPSON, WILLIAM H. SNYDER, JEFFREY C. WEIL
-
- Journal:
- Journal of Fluid Mechanics / Volume 417 / 25 August 2000
- Published online by Cambridge University Press:
- 25 August 2000, pp. 127-156
-
- Article
- Export citation
-
Laboratory experiments were conducted in a water tank to investigate the rise through the atmosphere of thermals generated by the detonation of surplus military munitions. The fall of a dense volume of fluid through the water in the tank follows the same governing equations as the rise of a buoyant volume of gas in the atmosphere. By filling the tank with a layer of water and a layer of salt water, an elevated step change in density was obtained, simulating a temperature inversion or jump in the atmosphere. The growth of the linear dimensions of the thermal and its volume were determined and used in the development of a criterion for predicting when a thermal will fully penetrate the inversion. Replacing the second layer of fluid with water of gradually increasing salinity, an elevated constant density gradient was obtained. In these cases, the maximum penetration distance of the thermal was observed and the equilibrium position and vertical spread were determined experimentally.
These observations and the empirical relationships determined from them should prove useful in the development and evaluation of air pollution dispersion models for predicting the atmospheric transport and diffusion of material released during such detonations.