首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
A mesoscale planetary boundary layer (PBL) numerical model has been developed to study airflow over complex topography. Turbulence closures using the turbulent kinetic energy (TKE) and dissipation () equations are investigated in combination with the level 2.5 scheme of Mellor and Yamada (1982) to determine eddy diffusivities for momentum and heat. This modified E- closure is simpler than the level 3 one which requires more prognostic equations for moist turbulent transport.One-dimensional (1-D) model results show that the PBL mean flows under various stability conditions are not significantly sensitive to the modified Blackadar and Kolmogorov eddy mixing-length formulations used in this E- model, although the latter yields excessively large mixing lengths in the entrainment region of the upper PBL. Eddy mixing lengths in the Kolmogorov-type formulation can be better defined by introducing background dissipation. Using the same prognostic TKE equation, the 1-D model results are not significantly affected by different diagnostic formulations in the closures. The simulated results compare well with large-eddy simulations and those obtained using higher-order closure schemes including the level 3 one. The results are found to be insensitive to eddy Prandtl number, in contrast to the 2-D model results (see Part II).  相似文献   

2.
We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.  相似文献   

3.
Here we advance the physical background of the energy- and flux-budget turbulence closures based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth’s rotation. In accordance with modern experimental evidence, the closure implies the maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: “strong turbulence” at ${Ri \ll 1}$ typical of boundary-layer flows and characterized by the practically constant turbulent Prandtl number Pr T; and “weak turbulence” at Ri > 1 typical of the free atmosphere or deep ocean, where Pr T asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy viscosity and eddy conductivity, and a general non-gradient model based on prognostic equations for all the basic parameters of turbulence including turbulent fluxes.  相似文献   

4.
A nonlinear numerical model is developed for turbulent boundary-layer flowover a train of water waves of finite amplitude or slope. The airflow isassumed to be steady, two-dimensional, and neutrally-stratified. The wavesurface is assumed to be aerodynamically rough and flow conditions at thewave surface are prescribed. The numerical model used in this study adoptsthree turbulence closure schemes with different degrees of physicalcompleteness. Two of these are second-order schemes, whichare believed to describe turbulent flow more completely than thesimpler closures used in most previous studies. Although models with all turbulence closures agree qualitatively in the prediction of the amplitude of the surface normal stress perturbation, the lower- and higher-order closures differ significantly in predictions of phase, and hence the form drag and energy transfer rate between wind and waves. Our model results are in reasonable agreement with field and laboratory measurements, although predicted energy transfer rates are generally at the low end of the range of experimental values. Cases with airflow at various angles to the wave direction are also considered.  相似文献   

5.
Second-Order Modelling of Turbulence in Katabatic Flows   总被引：1，自引：1，他引：0  
A complete one-dimensional second-order closure model is used to simulate katabatic flows observed on glaciers and ice caps. The model is tested with two different closure assumptions for the viscous dissipation, one based on a prognostic equation for and the other on a diagnostic buoyant length scale. Both formulations give quite similar results. Model simulations are compared to observations made over sloping ice surfaces during periods dominated by katabatic flow. In general, good agreement is found for both mean wind and temperature profiles as well as eddy correlation measurements. It is also found that the turbulent transport terms play an important role in katabatic flows as opposed to the classical stable boundary layer where these terms are usually ignored. Even the turbulent transport of temperature variance, which leads to the well-known countergradient term in unstable boundary layers, is relatively important for modelling the observed temperature profiles. The effect of these terms on the flux-profile relationships, using observed and simulated profiles, is also discussed.  相似文献   

6.
The validity of a spectral cumulus parameterization (spectral scheme) for simulating a diurnal cycle of precipitation over the Maritime Continent (MC) was examined using a regional atmospheric model. The impacts of entrainment parameterization and each type of convective closure, i.e., non-equilibrium (or equilibrium) closure for deep convection, mid-level, and shallow convective closures, were also examined. When vertically variable entrainment and appropriate convective closures were employed, the model adequately simulated a diurnal cycle of precipitation over both land and ocean as compared to the observation. Analysis regarding the entrainment parameterization revealed that variable entrainment parameterization was needed not only for simulating better mean patterns of precipitation, but also for more realistic phases of diurnal cycles. The impacts of convective closures appeared in the differences in the precipitation amplitude. Analysis on diurnal cycles of convective properties and tendencies revealed that the cycles between boundary layer forcing and convective heating determined convection strength and were affected by each type of convective closure. It can be concluded that the spectral scheme with appropriate convective closures is able to simulate a realistic diurnal cycle over the MC.  相似文献   

7.
Second-order closure models for the canopy sublayer (CSL) employ aset of closure schemes developed for `free-air' flow equations andthen add extra terms to account for canopy related processes. Muchof the current research thrust in CSL closure has focused on thesecanopy modifications. Instead of offering new closure formulationshere, we propose a new mixing length model that accounts for basicenergetic modes within the CSL. Detailed flume experiments withcylindrical rods in dense arrays to represent a rigid canopy areconducted to test the closure model. We show that when this lengthscale model is combined with standard second-order closureschemes, first and second moments, triple velocity correlations,the mean turbulent kinetic energy dissipation rate, and the wakeproduction are all well reproduced within the CSL provided thedrag coefficient (CD) is well parameterized. The maintheoretical novelty here is the analytical linkage betweengradient-diffusion closure schemes for the triple velocitycorrelation and non-local momentum transfer via cumulant expansionmethods. We showed that second-order closure models reproducereasonably well the relative importance of ejections and sweeps onmomentum transfer despite their local closure approximations.Hence, it is demonstrated that for simple canopy morphology (e.g.,cylindrical rods) with well-defined length scales, standard closureschemes can reproduce key flow statistics without much revision.When all these results are taken together, it appears that thepredictive skills of second-order closure models are not limitedby closure formulations; rather, they are limited by our abilityto independently connect the drag coefficient and the effectivemixing length to the canopy roughness density. With rapidadvancements in laser altimetry, the canopy roughness densitydistribution will become available for many terrestrialecosystems. Quantifying the sheltering effect, the homogeneity andisotropy of the drag coefficient, and more importantly, thecanonical mixing length, for such variable roughness density isstill lacking.  相似文献   

8.
A TKE-dissipation model for the atmospheric boundary layer   总被引：1，自引：0，他引：1  
The dissipation, , of turbulent kinetic energy (TKE) is a key parameter in atmospheric boundary-layer (ABL) models. Besides being a sink for momentum, it is often used together with the TKE to define an internal turbulence time scale for closure relations. A prognostic formulation for the dissipation of TKE is formulated, based on isotropic tensor modeling methods. The formulation is coupled to a level 2.5 second-order closure model and evaluated against measurements taken in horizontally homogeneous conditions, as well as against a tailored length-scale formulation. A formulation suitable for convective as well as neutral and stable ABLs is suggested.On leave from Department of Meteorology, Uppsala University, P.O. Box 516, S-751 20 Uppsala, Sweden.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

9.
A Second-Order Closure for Neutrally Stratified Vegetative Canopy Flows   总被引：1，自引：1，他引：0  
An existing second-order closure model is modified to include the effects on mean and turbulent motions of form and viscous drag in vegetative canopies. The additional physical mechanisms represented by the closure are viscous and pressure drag on canopy elements, their role in momentum absorption, in the creation of fine scale turbulent eddies and in enhancing the total viscous dissipation in the canopy airspace. Viscous dissipation is split into a standard 'isotropic contribution associated with the spectral eddy cascade and a foliage contribution associated with work against pressure and viscous drag on the foliage. Changes in the turbulent time scale that result from these mechanisms are included in the standard parameterisations of third moments and of the eddy cascade contribution to dissipation. The model is tested against a wind- tunnel 'wheat canopy, a corn canopy and a eucalypt canopy, a height range from 50 mm to 12.6 m. Model results show that the parameterisations of foliage interaction used in the closure are sufficiently robust to reproduce second-moment profiles within and above vegetative canopies to a high degree of accuracy without resorting to 'tuning of the model constants. The model also shows the natural emergence of two length scales, one associated with the familiar eddy cascade isotropic contribution to total dissipation and the other associated with the length scales of the canopy elements.  相似文献   

10.
  相似文献   

11.
The Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary-layer (PBL) scheme is a second-order turbulence closure model that is an improved version of the Mellor–Yamada scheme based on large-eddy simulation data. It simulates PBL structure and evolution well, particularly over the ocean surface. However, when used with various underlying surfaces in China, the scheme overestimates the turbulent momentum flux and the sensible heat flux. Based on observations of surface fluxes in China, we attempt to improve the MYNN model by modifying the parameters and representation of the turbulence scale. Closure constants and empirical expressions in the diagnostic equation are chosen first, and an additional component of the turbulent heat flux is considered in the potential temperature prognostic equation to improve the surface heat-flux modelling. The modified MYNN scheme is incorporated into a three-dimensional mesoscale model and is evaluated using various underlying surface observations. Amelioration of the surface turbulent fluxes is confirmed at five observational sites in China over different land-use types.  相似文献   

12.
An analytical model that predicts how much of the temperature–humidity covariance within the marine atmospheric surface layer (ASL) originates just above the ASL and just near the surface is proposed and tested using observations from the Risø Air Sea Experiment (RASEX). The model is based on a simplified budget for the two-scalar covariance that retains three basic terms: production, dissipation, and vertical transport. Standard second-order closure formulations are employed for the triple moments and the dissipation terms, and the canonical mixing length for the closure model is assumed linear with height (z) from the surface. Despite the poor performance of the gradient–diffusion closure in reproducing the measured triple moment, the overall covariance model was shown to be sufficiently robust to these assumptions. One of the main findings from the analytical treatment is the origin of the asymmetry in how the top and bottom boundary conditions affect the two-scalar covariance in the ASL. The analytical model reveals that ‘bottom-up’ boundary-condition variations scale with \(z^{-\sqrt{a}}\) , while ‘top-down’ variations scale with \(z^{\sqrt{a}}\) , where a is a constant that can be derived from similarity and closure constants. The genesis of this asymmetry stems from the flux-transport term but is modulated by the dissipation, and persists even in the absence of any inhomogeneity in the local production function. It is shown that the local production function acts to adjust the relative proportions of these two boundary conditions with weights that vary with the Obukhov length. The findings here do not provide ‘finality’ to the discussions on the covariance between humidity and temperature or the role of entrainment in modulating the turbulence within the ASL. Rather, they are intended to guide new hypotheses about interpretations of existing field data and identify needs for future field and numerical experiments.  相似文献   

13.
Significant improvements are occurring in the representation of physical processes in atmospheric convection models. They should go along with parallel improvements in the parameterization of subgrid scale turbulent processes. This problem appears to be particularly delicate in the presence of clouds, due to the local release of latent heat.Two important points are the choice of adequate turbulent thermodynamic variables and of the method for truncating the statistical moment equations. These topics are discussed here within the framework of the three-dimensional convection model under development at the Laboratoire de Météorologie Dynamique. Assuming the need for at least a simplified second-order closure, two improvements are tested on a numerical simulation of the Porto Rico experiment conducted by the National Center for Atmospheric Research (U.S.A.) in 1972. They concern the use of a rate equation for sub-grid scale turbulent kinetic energy and of specific variables which are approximately conserved in the condensation process.  相似文献   

14.
An Investigation of Higher-Order Closure Models for a Forested Canopy   总被引：11，自引：10，他引：1  
Simultaneous triaxial sonic anemometer velocity measurements vertically arrayed at six levels within and above a uniform pine forest were used to examine two parameterization schemes for the triple-velocity correlation tensor employed in higher-order closure models. These parameterizations are the gradient-diffusion approximation typically used in second-order closure models, and the full budget for the triple-velocity correlation tensor typically employed in third-order closure models. Both second- and third-order closure models failed to reproduce the measured profiles of the triple-velocity correlation within and above the canopy. However, the Reynolds stress tensor profiles (including velocity variances) deviated greatly from the measurements only within the lower levels of the canopy. It is shown that the Reynolds stresses are most sensitive to the parameterization of the triple-velocity correlation in these lower canopy regions where local turbulent production is negligible and turbulence is mainly sustained by the flux transport term. The failure of the third-order closure model to reproduce the measured third moments in the upper layers of the canopy-top contradicts conclusions from a previous study over shorter vegetation but agrees with another study for a deciduous forest. Whether the third-order closure model failure is due to the zero-fourth-cumulant closure approximation is therefore considered. Comparisons between measured and predicted quadruple velocity correlations suggest that the zero-fourth-cumulant approximation is valid close to the canopy-atmosphere in agreement with recent experiments.  相似文献   

15.
A high-resolution, nonhydrostatic, three-dimensional diagnostic PBL model over small-scale concave terrain was established in this paper. A two-dimensional prognostic model was developed based on the diagnostic model. The hydrostatic approximation was abandoned and the simple energy (E-ε) closure scheme was used in both mod-els. Using the two models, characteristics of PBL structure and its evolution were fully studied. The main characteris-tic of the PBL is the circulation, and it fairly affects the distribution of the pollutant in the pit.  相似文献   

16.
An analytical one-dimensional second-order closure model is developed to describe the within canopy velocity variances, turbulent intensities, dissipation rates, Lagrangian time scale and Lagrangian far field diffusivities for vegetation canopies of arbitrary structure and density. The model incorporates and extends the model of momentum transfer developed by Massman (1997) and the model of within canopy velocity variances developed by Weil (unpublished) from the second-order closure model of Wilson and Shaw (1977). Model predictions of within and above canopy velocity variances, turbulent intensities, dissipation rates and the Lagrangian time scale are in reasonable agreement with previously measured or estimated values for these parameters. The present model suggests that the Lagrangian time scale and the far field diffusivity could be strongly dependent upon foliage structure and density through the foliage effects on the velocity variances. A simple formulation for the Lagrangian time scale at canopy height is derived from model results. Taken as a whole, the present model may provide a relatively simple way to incorporate turbulence parameters into models of soil/canopy/atmosphere mass transfer.  相似文献   

17.
Marine stratocumulus observations show a large variability in cloud droplet number concentration (CDNC) related to variability in aerosol concentration. Changes in CDNC modify the cloud reflectivity, but also affect cloud water content, cloud lifetime, and cloudiness, through changes in precipitation. In mesoscale models and general circulation models (GCMs), precipitation mechanisms are parameterized. Here we examine how the precipitation parameterization can affect the simulated cloud. Simulations are carried out with the one-dimensional version of the hydrostatic primitive equation model MAR (Modéle Atmosphérique Régional) developed at the Université catholique de Louvain. It includes a E- turbulence closure, a wide-band formulation of the radiative transfer, and a parameterized microphysics including prognostic equations for water vapour, cloud droplets and rain drops concentrations. In a first step, the model is used to simulate a horizontally homogeneous stratocumulus deck observed during the Atlantic Stratocumulus Transition Experiment (ASTEX) on the night of 12–13 June 1992. The observations show that the model is able to realistically reproduce the vertical structure of the cloud-topped boundary layer. In a second step, several precipitation parameterizations commonly used in mesoscale models and GCMs are tested. It is found that most parameterizations tend to overestimate the precipitation, which results in an underestimation of the vertically integrated liquid water content. Afterwards, using those parameterizations that are sensitive to CDNC, several simulations are performed to estimate the effect of CDNC variations on the simulated cloud. Based upon the simulation results, we argue that currently used parameterizations do not enable assessment of such a sensitivity.  相似文献   

18.
In this work, three turbulence closure models, Mellor andYamada level 2.5, E - l and E - implemented in a circulation model, are compared in neutral condition over complex terrain. They are firstly applied to a one-dimensional case on flat terrain and then to a schematic two-dimensional valley. The simulation results, in terms of wind field and turbulent kinetic energy, are tested against measurements from a wind-tunnel experiment. The empirical constants defining the characteristic length scales of the closures are modified based on turbulence parameters estimated in the experiment. The formulation of the diffusion coefficients is analysed to explain the differences among the various closures in the simulation results. Regarding the mean flow, both on flat and complex terrain, all the closures yield satisfactory results. Concerning the turbulent kinetic energy, the best results are obtained by E - l and E - closures.  相似文献   

19.
We perform large-eddy simulation (LES) of a moderately convective atmospheric boundary layer (ABL) using a prognostic subfilter-scale (SFS) model obtained by truncating the full conservation equations for the SFS stresses and fluxes. The truncated conservation equations contain production mechanisms that are absent in eddy-diffusivity closures and, thus, have the potential to better parametrize the SFS stresses and fluxes. To study the performance of the conservation-equation-based SFS closure, we compare LES results from the surface layer with observations from the Horizontal Array Turbulence Study (HATS) experiment. For comparison, we also show LES results obtained using an eddy-diffusivity closure. Following past studies, we plot various statistics versus the non-dimensional parameter, Λ w /Δ, where Λ w is the wavelength corresponding to the peak in the vertical velocity spectrum and Δ is the filter width. The LES runs are designed using different domain sizes, filter widths and surface fluxes, in order to replicate partly the conditions in the HATS experiment. Our results show that statistics from the different LES runs collapse reasonably and exhibit clear trends when plotted against Λ w /Δ. The trends exhibited by the production terms in the modelled SFS conservation equations are qualitatively similar to those seen in the HATS data with the exception of SFS buoyant production, which is underpredicted. The dominant production terms in the modelled SFS stress and flux budgets obtained from LES are found to approach asymptotically constant values at low Λ w /Δ. For the SFS stress budgets, we show that several of these asymptotes are in good agreement with their corresponding theoretical values in the limit Λ w /Δ → 0. The modelled SFS conservation equations yield trends in the mean values and fluctuations of the SFS stresses and fluxes that agree better with the HATS data than do those obtained using an eddy-diffusivity closure. They, however, underpredict considerably the level of SFS anisotropy near the wall when compared to observations, which could be a consequence of the shortcomings in the model used for the pressure destruction terms. Finally, we address the computational cost incurred due to the use of additional prognostic equations.  相似文献   

20.
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.  相似文献   

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

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