Large-Eddy Simulations of Radiatively Driven Convection: Sensitivities to the Representation of Small Scales

Bjorn Stevens National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Bjorn Stevens in
Current site
Google Scholar
PubMed
Close
,
Chin-Hoh Moeng National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Chin-Hoh Moeng in
Current site
Google Scholar
PubMed
Close
, and
Peter P. Sullivan National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Peter P. Sullivan in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

Large-eddy simulations of a smoke cloud are examined with respect to their sensitivity to small scales as manifest in either the grid spacing or the subgrid-scale (SGS) model. Calculations based on a Smagorinsky SGS model are found to be more sensitive to the effective resolution of the simulation than are calculations based on the prognostic turbulent kinetic energy (TKE) SGS model. The difference between calculations based on the two SGS models is attributed to the advective transport, diffusive transport, and/or time-rate-of-change terms in the TKE equation. These terms are found to be leading order in the entrainment zone and allow the SGS TKE to behave in a way that tends to compensate for changes that result in larger or smaller resolved scale entrainment fluxes. This compensating behavior of the SGS TKE model is attributed to the fact that changes that reduce the resolved entrainment flux (viz., values of the eddy viscosity in the upper part of the PBL) simultaneously tend to increase the buoyant production of SGS TKE in the radiatively destabilized portion of the smoke cloud. Increased production of SGS TKE in this region then leads to increased amounts of transported, or fossil, SGS TKE in the entrainment zone itself, which in turn leads to compensating increases in the SGS entrainment fluxes. In the Smagorinsky model, the absence of a direct connection between SGS TKE in the entrainment and radiatively destabilized zones prevents this compensating mechanism from being active, and thus leads to calculations whose entrainment rate sensitivities as a whole reflect the sensitivities of the resolved-scale fluxes to values of upper PBL eddy viscosities.

Corresponding author address: Mr. Bjorn Stevens, Dept. of Atmospheric Sciences, University of California, Los Angeles, Los Angeles, CA 90095.

Abstract

Large-eddy simulations of a smoke cloud are examined with respect to their sensitivity to small scales as manifest in either the grid spacing or the subgrid-scale (SGS) model. Calculations based on a Smagorinsky SGS model are found to be more sensitive to the effective resolution of the simulation than are calculations based on the prognostic turbulent kinetic energy (TKE) SGS model. The difference between calculations based on the two SGS models is attributed to the advective transport, diffusive transport, and/or time-rate-of-change terms in the TKE equation. These terms are found to be leading order in the entrainment zone and allow the SGS TKE to behave in a way that tends to compensate for changes that result in larger or smaller resolved scale entrainment fluxes. This compensating behavior of the SGS TKE model is attributed to the fact that changes that reduce the resolved entrainment flux (viz., values of the eddy viscosity in the upper part of the PBL) simultaneously tend to increase the buoyant production of SGS TKE in the radiatively destabilized portion of the smoke cloud. Increased production of SGS TKE in this region then leads to increased amounts of transported, or fossil, SGS TKE in the entrainment zone itself, which in turn leads to compensating increases in the SGS entrainment fluxes. In the Smagorinsky model, the absence of a direct connection between SGS TKE in the entrainment and radiatively destabilized zones prevents this compensating mechanism from being active, and thus leads to calculations whose entrainment rate sensitivities as a whole reflect the sensitivities of the resolved-scale fluxes to values of upper PBL eddy viscosities.

Corresponding author address: Mr. Bjorn Stevens, Dept. of Atmospheric Sciences, University of California, Los Angeles, Los Angeles, CA 90095.

Save
  • Bretherton, C. S., and Coauthors, 1999: An intercomparison of radiatively driven entrainment and turbulence in a smoke cloud, as simulated by different numerical models. Quart. J. Roy. Meteor. Soc.,125, 391–423.

  • Clark, R. A., J. H. Ferziger, and W. C. Reynolds, 1979: Evaluation of subgrid-scale models using an accurately simulated turbulent flow. J. Fluid Mech.,91, 1–16.

  • Deardorff, J. W., 1980: Stratocumulus-capped mixed layers derived from a three-dimensional model. Bound.-Layer Meteor.,18, 495–527.

  • ——, G. E. Willis, and B. H. Stockton, 1980: Laboratory studies of the entrainment zone of a convectively mixed layer. J. Fluid Mech.,100, 41–64.

  • Ghosal, S., 1996: An analysis of numerical errors in large-eddy simulations of turbulence. J. Comput. Phys.,125, 187–206.

  • Kraus, E. B., 1963: The diurnal precipitation change over the sea. J. Atmos. Sci.,20, 551–556.

  • Lenderink, G., M. C. vanZanten, and P. G. Duynkerke, 1999: Can an e–l turbulence closure simulate entrainment in convective boundary layers? J. Atmos. Sci.,56, 3331–3337.

  • Lewellen, D., and W. Lewellen, 1998: Large-eddy boundary layer entrainment. J. Atmos. Sci.,55, 2645–2665.

  • Lilly, D. K., 1967: The representation of small-scale turbulence in numerical simulation experiments. IBM Scientific Computing Symp. on Environmental Sciences, Yorktown Heights, NY, IBM DP Division, 195–210.

  • ——, 1968: Models of cloud topped mixed layers under a strong inversion. Quart. J. Roy. Meteor. Soc.,94, 292–309.

  • Lock, A. P., and M. K. MacVean, 1999: A parametrization of entrainment driven by surface heating and cloud-top cooling. Quart. J. Roy. Meteor. Soc.,125, 271–300.

  • Mason, P., 1989: Large eddy simulation of the convective atmospheric boundary layer. J. Atmos. Sci.,46, 1492–1516.

  • ——, and A. R. Brown, 1999: On subgrid models and filter operations in large eddy simulations. J. Atmos. Sci.,56, 2101–2114.

  • Moeng, C.-H., 1984: A large-eddy simulation for the study of planetary boundary layer turbulence. J. Atmos. Sci.,41, 2052–2062.

  • ——, and U. Schumann, 1991: Composite structure of plumes in stratus-topped boundary layers. J. Atmos. Sci.,48, 2280–2291.

  • ——, P. P. Sullivan, and B. Stevens, 1999: Including radiative effects in an entrainment rate formula for buoyancy-driven PBLs. J. Atmos. Sci.,56, 1031–1049.

  • ——, and J. C. Wyngaard, 1989: Spectral analysis of large-eddy simulations of the convective boundary layer. J. Atmos. Sci.,45, 3573–3587.

  • Nieuwstadt, F. T. M., P. J. Mason, C. H. Moeng, and U. Schumann, 1991: Large-eddy simulation of the convective boundary layer:A comparison of four computer codes. Selected Papers From the 8th Symposium on Turbulent Shear Flows, Springer-Verlag, 343–367.

  • Patankar, S. V., 1980: Numerical Heat Transfer and Fluid Flow. Series in Computational Methods in Mechanics and Thermal Sciences. McGraw Hill, 197 pp.

  • Pielke, R. A., and Coauthors, 1992: A comprehensive meteorological modeling system—RAMS. Meteor. Atmos. Phys.,49, 69–91.

  • Rouse, H., and J. Dodu, 1955: Turbulent diffusion across a density discontinuity. Houille Blanche,10, 522–532.

  • Sayler, B. J., and R. E. Breidenthal, 1998: Laboratory simulations of radiatively induced entrainment in stratiform clouds. J. Geophys. Res.,103, 8827–8837.

  • Schmidt, H., and U. Schumann, 1989: Coherent structure of the convective boundary layer derived from large-eddy simulations. J. Fluid Mech.,200, 511–562.

  • Schumann, U., 1991: Subgrid length-scales for large-eddy simulation of stratified turbulence. Theoret. Comput. Fluid Dyn.,2, 279–290.

  • Sommeria, G., 1976: Three-dimensional simulation of turbulent processes in an undisturbed trade wind boundary layer. J. Atmos. Sci.,33, 216–241.

  • Stevens, B., W. R. Cotton, G. Feingold, and C.-H. Moeng, 1998: Large-eddy simulations of strongly precipitating, shallow, stratocumulus-topped boundary layers. J. Atmos. Sci.,55, 3616–3638.

  • ——, C.-H. Moeng, and P. P. Sullivan, 1999: Entrainment and subgrid lengthscales in large-eddy simulations of atmospheric boundary layer flows. Developments in Geophysical Turbulence, R. Kerr and Y. Kimura, Eds., Kluwer, in press.

  • Stevens, D., and C. S. Bretherton, 1999: Effects of resolution on the simulation of stratocumulus entrainment. Quart. J. Roy. Meteor. Soc.,125, 425–439.

  • Sullivan, P. P., J. C. McWilliams, and C.-H. Moeng, 1996: A grid nesting method for large-eddy simulation of planetary boundary-layer flows. Bound.-Layer Meteor.,80, 167–202.

  • ——, C.-H. Moeng, B. Stevens, D. H. Lenschow, and S. D. Mayor, 1998: Entrainment and structure of the inversion layer in the convective planetary boundary layer. J. Atmos. Sci.,55, 3042–3064.

  • Tremback, C., J. Powell, W. R. Cotton, and R. A. Pielke, 1987: The forward-in-time upstream advection scheme: Extension to higher orders. Mon. Wea. Rev.,115, 3540–555.

  • Turner, J. S., 1986: Turbulent entrainment: The development of the entrainment assumption, and its application to geophysical flows. J. Fluid Mech.,173, 431–471.

  • vanZanten, M. C., P. G. Duynkerke, and J. W. M. Cuijpers, 1999: Entrainment parameterization in convective boundary layers derived from large eddy simulations. J. Atmos. Sci.,56, 813–828.

  • Zalesak, S. T., 1979: Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys.,31, 335–362.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 837 223 36
PDF Downloads 504 137 11