首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 171 毫秒
1.
An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer   总被引：2，自引：27，他引：2  
Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.  相似文献   

2.
Wind-tunnel simulations of theatmospheric stable boundary layer (SBL) developedover a rough surface were conducted by using athermally stratified wind tunnel at the Research Institutefor Applied Mechanics (RIAM), Kyushu University. Thepresent experiment is a continuation of the workcarried out in a wind tunnel at Colorado StateUniversity (CSU), where the SBL flows were developed over asmooth surface. Stably stratified flows were createdby heating the wind-tunnel airflow to a temperature ofabout 40–50°and by cooling the test-section floor toa temperature of about 10°. To simulate therough surface, a chain roughness was placed over thetest-section floor. We have investigated the buoyancyeffect on the turbulent boundary layer developed overthis rough surface for a wide range of stability,particularly focusing on the turbulence structure andtransport process in the very stable boundary layer.The present experimental results broadly confirm theresults obtained in the CSU experiment with the smoothsurface, and emphasizes the following features: thevertical profiles of turbulence statistics exhibitdifferent behaviour in two distinct stability regimes with weak and strong stability,corresponding to the difference in the verticalprofiles of the local Richardson number. The tworegimes are separated by the critical Richardsonnumber. The magnitudes in turbulence intensities andturbulent fluxes for the weak stability regime aremuch greater than those of the CSU experiments becauseof the greater surface roughness. For the very stableboundary layer, the turbulent fluxes of momentum andheat tend to vanish and wave-like motions due to theKelvin–Helmholtz instability and the rolling up andbreaking of those waves can be observed. Furthermore,the appearance of internal gravity waves is suggestedfrom cross-spectrum analyses.  相似文献   

3.
The experiment IGLOS (Investigation of the Greenland Boundary Layer Over Summit) was conducted in June and July 2002 in the central plateau of the Greenland inland ice. The German research aircraft Polar2, equipped with the turbulence measurement system Meteopod, was used to investigate turbulence and radiation flux profiles near research station “Summit Camp”. Aircraft measurements are combined with measurements of radiation fluxes and turbulent quantities made from a 50 m tower at Summit Camp operated by Eidgenössische Technische Hochschule (ETH) Zürich. During all six flight missions, well-developed stable boundary layers were found. Even in high-wind conditions, the surface inversion thickness did not exceed roughly 100 m. The turbulent height of the stable boundary layer (SBL) was found to be much smaller than the surface inversion thickness. Above the surface layer, significant turbulent fluxes occurred only intermittently in intervals on the order of a few kilometres. Turbulent event fraction in the upper SBL shows the same dependence on gradient Richardson number as reported for near-surface measurements. Clear-air longwave radiation divergence was always found to contribute significantly to the SBL heat budget. In low-wind cases, radiative cooling even turned out to be dominant.  相似文献   

4.
We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.  相似文献   

5.
Large-eddy simulations (LES) of the Stable Atmospheric Boundary Layer (SBL) are difficult because the turbulence is not isotropic for strong stratification and the Kolmogorov theory might be no longer valid. This fact compells us to work on modifications to the subgrid turbulence schemes, although there is not any widely accepted theory on anisotropic turbulence. In this work, a LES model is used to see what range of stable stratification can still be simulated with a subgrid turbulence scheme using the Kolmogorov theory for the dissipation. Twenty simulations of increasing stability have been performed using a horizontal resolution of 5 m. The model is able to simulate weakly and moderately stable conditions and experiences runaway cooling for strong stability. The goodness of the successful simulations is inspected through comparison to observations from the experimental campaigns SABLES-98 and CASES-99. Other supplementary tests have been performed on the resolution and the surface boundary condition.  相似文献   

6.
Turbulent and mean meteorological data collected at five levels on a 20-m tower over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are analyzed to examine different regimes of the stable boundary layer (SBL). Eleven months of measurements during SHEBA cover a wide range of stability conditions, from the weakly unstable regime to very stable stratification. Scaling arguments and our analysis show that the SBL can be classified into four major regimes: (i) surface-layer scaling regime (weakly stable case), (ii) transition regime, (iii) turbulent Ekman layer, and (iv) intermittently turbulent Ekman layer (supercritical stable regime). These four regimes may be considered as the basic states of the traditional SBL. Sometimes these regimes, especially the last two, can be markedly perturbed by gravity waves, detached elevated turbulence (‘upside down SBL’), and inertial oscillations. Traditional Monin–Obukhov similarity theory works well in the weakly stable regime. In the transition regime, Businger–Dyer formulations work if scaling variables are re-defined in terms of local fluxes, although stability function estimates expressed in these terms include more scatter compared to the surface-layer scaling. As stability increases, the near-surface turbulence is affected by the turning effects of the Coriolis force (the turbulent Ekman layer). In this regime, the surface layer, where the turbulence is continuous, may be very shallow (< 5 m). Turbulent transfer near the critical Richardson number is characterized by small but still significant heat flux and negligible stress. The supercritical stable regime, where the Richardson number exceeds a critical value, is associated with collapsed turbulence and the strong influence of the earth’s rotation even near the surface. In the limit of very strong stability, the stress is no longer a primary scaling parameter.  相似文献   

7.
Turbulent statistics of neutrally stratified shear-driven flow within and above a sparse forest canopy are presented from a large-eddy simulation (LES) and compared with those from observations within and above a deciduous forest with similar height and foliage density. First- and second-order moments from the LES agree with observations quite well. Third-order moments from the LES have the same sign and similar vertical patterns as those from the observations, but the LES yields smaller magnitudes of such higher-order moments. Turbulent spectra and cospectra from the LES agree well with observations above the forest. However, at the highest frequencies, the LES spectra have steeper slopes than observations. Quadrant and conditional analyses of the LES resolved-scale flow fields also agree with observations. For example, both LES and observation find that sweeps are more important than ejections for the transport of momentum within the forest, while inward and outward interaction contributions are both small, except near the forest floor. The intermittency of the transport of momentum and scalar increases with depth into the forest. Finally, ramp structures in the time series of a passive scalar at multiple levels within and above the forest show similar features to those measured from field towers. Two-dimensional (height-time cross-section) contours of the passive scalar and wind vectors show sweeps and ejections, and the characteristics of the static pressure perturbation near the ground resemble those deduced from field tower-based measurements. In spite of the limited grid resolution (2 m × 2 m × 2 m) and domain size (192 m × 192 m × 60 m) used in this LES, we demonstrate that the LES is capable of resolving the most important characteristics of the turbulent flow within and above a forest canopy.  相似文献   

8.
为研究黄河源区边界层湍流特征及其对物质和能量输送的影响,本文首次采用大涡模拟的方法,对比分析了黄河源区两种不同下垫面上（鄂陵湖和湖边草地）对流边界层（CBL）中精细的湍流结构特征。使用资料为2012年夏季黄河源区鄂陵湖流域野外观测实验的GPS探空资料、涡动相关观测资料。分析表明,模拟的黄河源区草地和湖上CBL的平均结构与实测结果吻合较好,但草地和湖上CBL的湍流结构特征差异较明显。模拟结果显示,草地CBL内湍能收支、湍流特征量的时空分布和湍涡结构特征均与陆地上热力驱动CBL的研究结果一致;湖上CBL顶部存在明显的对流卷特征,且夹卷层的湍流强度比草地的强,而草地近地面湍强则更大。通过改变水平分辨率的模拟试验,发现两个不同下垫面上模拟结果对模式分辨率的敏感性不同,湖面CBL的模拟要选择较高的水平分辨率（50～100 m）,以提高近湖面和夹卷层对湍流动能和湍流通量模拟的精度,也充分模拟出各种尺度的波对湍流通量的累积贡献。考虑到计算时间等影响,模拟草地边界层精细的湍流结构时建议选择网格距为100～200 m。  相似文献   

9.
This case study introduces measurements of turbulent fluxes in a nocturnal boundary layer in North Germany with the new helicopter-borne turbulence measurement system HELIPOD, a detailed data analysis and examination in regard of systematic errors of the instrument, and some comparison with local similarity theory and experiments of the past, in order to confirm the occurrence of small vertical turbulent fluxes. The examined nocturnal boundary layer offered excellent conditions to analyse the quality of the measurement system. In this connection, a detailed look at a strong ground-based inversion disclosed small turbulent fluxes with a spectral maximum at ten metres wavelength or less, embedded in intermittent turbulence. For verification of these fluxes, the measurements were compared with well established results from past experiments. Local similarity theory was applied to calculate dimensionless variances of the turbulent quantities, which were found in good agreement with other observations. Since shear and stratification varied significantly on the horizontal flight legs due to global intermittency, a method was developed to determine vertical gradients on a horizontal flight pattern, by use of small fluctuations of the measurement height. With these locally determined gradients, gradient transport theory became applicable and the turbulent diffusivities for heat and momentum, the Richardson number, and the flux Richardson number were estimated within isolated strong turbulent outbursts. Within these outbursts the flux Richardson number was found between 0.1 and 0.2. The functional relationship between the gradient Richardson number and the turbulent Prandtl number agreed well with observations in past experiments and large eddy simulation. The impact of the stratification on the vertical turbulent exchange, as already described for the surface layer using Monin–Obukhov similarity, was analogously observed in the very stably stratified bulk flow when local scaling was applied.  相似文献   

10.
  相似文献   

11.
A method for performing nested grid calculations with a large-eddy simulation code is described. A common numerical method is used for all meshes, and the grid architecture consists of a single outer or coarse grid, and nested or fine grids, which overlap in some common region. Inter-grid communication matches the velocity, pressure and potential temperature fields in the overlap region. Resolved and sub-grid scale (SGS) turbulent fluxes and kinetic energy on the fine grid are averaged to the coarse grid using a conservation rule equivalent to Germano's identity used to develop dynamic SGS models.Simulations of a slightly convective, strong shear planetary boundary layer were carried out with varying surface-layer resolutions. Grid refinements in the (x, y, z) directions of up to (5, 5, 2) times were employed. Two-way interaction solutions on the coarse and fine meshes are successfully matched in the overlap region on an instantaneous basis, and the turbulent motions on the fine grid blend smoothly into the coarse grid across the grid interface. With surface-layer grid nesting, significant increases in resolved eddy fluxes and variances are found. The energy-scale content of the vertical velocity, and hence vertical turbulent fluxes, appear to be most influenced by increased grid resolution. Vertical velocity spectra show that the dominant scale shifts towards higher wavenumbers (smaller scales) and the magnitude of the peak energy is increased by more than a factor of 3 with finer resolution. Outside of the nested region the average heat and momentum fluxes and spectra are slightly influenced by the fine resolution in the surface layer. From these results we conclude that fine resolution is required to resolve the details of the turbulent motions in the surface layer. At the same time, however, increased resolution in the surface layer does not appreciably alter the ensemble statistics of the resolved and SGS motions outside of the nested region.  相似文献   

12.
Observations of surface-layer turbulence and turbulent fluxes were made over a desert in northwestern China as a part of HEIFE (HEIhe river Field Experiment). These show that the normalized variations of the vertical wind component and of the air temperature obey Monin-Obukhov similarity well, especially in free convective conditions. However, the variations of specific humidity do not obey Monin-Obukhov similarity.Mean bulk transfer coefficients of sensible heat and momentum flux are obtained as functions of stability over a wide stability range from the observed data of turbulent fluxes and mast profiles. However, the bulk transfer coefficient for water vapor could not be obtained because of the large scatter of the data. In free convective conditions, the sensible heat flux was found to be approximately proportional to the 1.4 power of temperature difference between the surface and 20m. The bulk transfer coefficient of sensible heat is also obtained as a function of the bulk Richardson number for practical convenience.  相似文献   

13.
The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.  相似文献   

14.
This paper compares a number of one-dimensional closure models for the planetary boundary layer (PBL) that are currently in use in large-scale atmospheric models. Using the results of a large-eddy simulation (LES) model as the standard of comparison, the PBL models are evaluated over a range of stratifications from free convective to neutral and a range of surface shear stresses. Capping inversion strengths for the convective cases range from weakly to strongly capped. Six prototypical PBL models are evaluated in this study, which focuses on the accuracy of the boundary-layer fluxes of momentum, heat, and two passive scalars. One scalar mimics humidity and the other is a top-down scalar entrained into the boundary layer from above. A set of measures based on the layer-averaged differences of these fluxes from the LES solutions is developed. In addition to the methodological framework and suite of LES solutions, the main result of the evaluation is the recognition that all of the examined PBL parameterizations have difficulty reproducing the entrainment at the top of the PBL, as given by the LES, in most parameter regimes. Some of the PBL models are relatively accurate in their entrainment flux in a subset of parameter regimes. The sensitivity of the PBL models to vertical resolution is explored, and substantive differences are observed in the performance of the PBL models, relative to LES, at low resolution typical of large scale atmospheric models.  相似文献   

15.
A long-standing problem in large-eddy simulations (LES) of the planetary boundary layer (PBL) is that the mean wind and temperature profiles differ from the Monin-Obukhov similarity forms in the surface layer. This shortcoming of LES has been attributed to poor grid resolution and inadequate sub-grid-scale (SGS) modeling. We study this deficiency in PBL LES solutions calculated over a range of shear and buoyancy forcing conditions. The discrepancy from similarity forms becomes larger with increasing shear and smaller buoyancy forcing, and persists even with substantial horizontal grid refinement. With strong buoyancy forcing, however, the error is negligible.In order to achieve better agreement between LES and similarity forms in the surface layer, a two-part SGS eddy-viscosity model is proposed. The model preserves the usual SGS turbulent kinetic energy formulation for the SGS eddy viscosity, but it explicitly includes a contribution from the mean flow and a reduction of the contributions from the turbulent fluctuations near the surface. Solutions with the new model yield increased fluctuation amplitudes near the surface and better correspondence with similarity forms out to a distance of 0.1–0.2 times the PBL depth, i.e., a typical surface-layer depth. These results are also found to be independent of grid anisotropy. The new model is simple to implement and computationally inexpensive.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

16.
We use the mesoscale meteorological model Meso-NH, taking the drag force of trees into account under stable, unstable and neutral conditions in a real case study. Large-eddy simulations (LES) are carried out for real orography, using a regional forcing model and including the energy and water fluxes between the surface (mostly grass with some hedges of trees) and the atmosphere calculated using a state-of-the-art soil-vegetation-atmosphere-transfer model. The formulation of the drag approach consists of adding drag terms to the momentum equation and subgrid turbulent kinetic energy dissipation, as a function of the foliage density. Its implementation in Meso-NH is validated using Advanced Regional Prediction System simulation results and measurements from Shaw and Schumann (Boundary-Layer Meteorol, 61(1):47?C64, 1992). The simulation shows that the Meso-NH model successfully reproduces the flow within and above homogeneous covers. Then, real case studies are used in order to investigate the three different boundary layers in a LES configuration (resolution down to 2 m) over the ??Lannemezan 2005?? experimental campaign. Thus, we show that the model is able to reproduce realistic flows in these particular cases and confirm that the drag force approach is more efficient than the classical roughness approach in describing the flow in the presence of vegetation at these resolutions.  相似文献   

17.
A three-dimensional numerical meteorological model is used to perform large-eddy simulations of the upslope flow circulation over a periodic ridge-valley terrain. The subgrid-scale quantities are modelled using a prognostic turbulence kinetic energy (TKE) scheme, with a grid that has a constant horizontal resolution of 50 m and is stretched along the vertical direction. To account for the grid anisotropy, a modified subgrid length scale is used. To allow for the response of the surface fluxes to the valley-flow circulation, the soil surface temperature is imposed and the surface heat and momentum fluxes are computed based on Monin–Obukhov similarity theory. The model is designed with a symmetrical geometry using periodic boundary conditions in both the x and y directions. Two cases are simulated to study the influence of along-valley geostrophic wind forcing with different intensities. The presence of the orography introduces numerous complexities both in the mean properties of the flow and in the turbulent features, even for the idealized symmetric geometry. Classical definitions for the height of the planetary boundary layer (PBL) are revisited and redefined to capture the complex structure of the boundary layer. Analysis of first- and second-moment statistics, along with TKE budget, highlights the different structure of the PBL at different regions of the domain.  相似文献   

18.
The surface layer of an atmospheric boundary layer (ABL) is most accessible to field measurements and hence its ensemble-mean structure has been well established. The Kansas field measurements were the first detailed study of this layer, providing numerous benchmark statistical profiles for a wide range of stability states. Large-eddy simulation (LES), in contrast, is most suitable for studying the mixed layer of the ABL where the energy-containing range of the vertical velocity field is well resolved. In the surface layer, typical large-eddy simulations barely resolve the energy-containing vertical-velocity fields and hence do not provide sufficient data for a detailed analysis.We carried out a nested-mesh simulation of a moderately convective ABL (-zi/L = 8) in which the lower 6% of the boundary layer had an effective grid resolution of 5123. We analyze the LES fields above the 6th vertical grid level (z = 23 m) where the vertical velocity field has a well formed inertial subrange, for a detailed comparison with the Kansas results. Various terms in the budgets of turbulent kinetic energy, temperature variance, Reynolds stress, temperature flux, and some higher order moments are compared. The agreement is generally quite good; however, we do observe certain discrepancies, particularly in the terms involving pressure fluctuations.  相似文献   

19.
Numerical weather prediction (NWP) model forecasts at horizontal grid lengths in the range of 100 m to 1 km are now possible. Within this range of grid lengths, the convective boundary layer (CBL) is partially resolved and thus in the so-called ‘grey zone’. For simulations in the grey zone, numerical dissipation sources from both the advection scheme and the subgrid model are likely to be significant. Until now, these effects have not been incorporated fully into our understanding of the grey zone. In order to quantify these effects, a dissipation length scale is defined based on the second moment of the turbulent kinetic energy (TKE) spectrum. An ensemble of simulations of a CBL are performed using a large-eddy model across the grey-zone resolutions and for a range of subgrid model, advection scheme and vertical grid configurations. The dissipation length scale distinguishes the effects of the different model configurations in the grey zone. In the middle of the boundary layer, the resolved TKE is strongly controlled by the numerical dissipation. This leads to a similarity law for the resolved TKE in the grey zone using the dissipation length scale. A new definition of the grey zone emerges where the inversion depth and dissipation length scale are the same size. This contrasts with the typical definition using the horizontal grid length. At the inversion, however, the variation of the dissipation length scale with grid length is less predictable, reflecting significant challenges for modelling entrainment in the grey zone. The dissipation length scale is thus a simple diagnostic to aid both NWP and large-eddy modellers in understanding the grey zone.  相似文献   

20.
A single-column model of the evolving stable boundary layer (SBL) is tested for self-similar properties of the flow and effects of ambient forcing. The turbulence closure of the model is diagnostic, based on the K-theory approach, with a semi-empirical form of the mixing length, and empirical stability functions of the Richardson number. The model results, expressed in terms of local similarity scales, are universal functions, satisfied in the entire SBL. Based on similarity expression, a realizability condition is derived for the minimum allowable turbulent heat flux in the SBL. Numerical experiments show that the development of “horse-shoe” shaped, fixed-elevation hodographs in the interior of the SBL around sunrise is controlled by effects imposed by surface thermal forcing.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号