首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 394 毫秒
1.
Large-eddy simulations (LES) of the continuously turbulent quasi-equilibrium stable boundary layer (SBL) are conducted with grid lengths in the range of 12.5 m to 2 m, in order to explore resolution sensitivity, and determine at what point grid convergence occurs. The structure of the mean potential temperature, winds, and turbulent fluxes varies significantly over this resolution range. The highest resolution simulations show a significant degree of convergence. The dimensionless momentum diffusivity asymptotes to a value of 0.06, corresponding to a limiting flux Richardson number of 0.15.Using the converged simulations, some scaling hypotheses underpinning first-order and second-order closure models are revisited. The effective Richardson number stability functions of the LES are compared with the forms often used in numerical weather prediction (NWP). The mixing implied by the LES is less than that used in NWP. The commonly used similarity profiles for heat and momentum fluxes, and the scalings for dissipation and pressure covariances are compared with the LES. This information could provide guidance for the next generation of SBL parametrization schemes.  相似文献   

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

3.
We describe and analyze the results of the third global energy and water cycle experiment atmospheric boundary layer Study intercomparison and evaluation study for single-column models. Each of the nineteen participating models was operated with its own physics package, including land-surface, radiation and turbulent mixing schemes, for a full diurnal cycle selected from the Cabauw observatory archive. By carefully prescribing the temporal evolution of the forcings on the vertical column, the models could be evaluated against observations. We focus on the gross features of the stable boundary layer (SBL), such as the onset of evening momentum decoupling, the 2-m minimum temperature, the evolution of the inertial oscillation and the morning transition. New process diagrams are introduced to interpret the variety of model results and the relative importance of processes in the SBL; the diagrams include the results of a number of sensitivity runs performed with one of the models. The models are characterized in terms of thermal coupling to the soil, longwave radiation and turbulent mixing. It is shown that differences in longwave radiation schemes among the models have only a small effect on the simulations; however, there are significant variations in downward radiation due to different boundary-layer profiles of temperature and humidity. The differences in modelled thermal coupling to the land surface are large and explain most of the variations in 2-m air temperature and longwave incoming radiation among models. Models with strong turbulent mixing overestimate the boundary-layer height, underestimate the wind speed at 200 m, and give a relatively large downward sensible heat flux. The result is that 2-m air temperature is relatively insensitive to turbulent mixing intensity. Evening transition times spread 1.5 h around the observed time of transition, with later transitions for models with coarse resolution. Time of onset in the morning transition spreads 2 h around the observed transition time. With this case, the morning transition appeared to be difficult to study, no relation could be found between the studied processes, and the variation in the time of the morning transition among the models.  相似文献   

4.
The dimensionless bottom-up and top-down gradient functions in the convective boundary layer (CBL) are evaluated utilizing long-term well-calibrated carbon dioxide mixing ratio and flux measurements from multiple levels of a 447-m tall tower over a forested area in northern Wisconsin, USA. The estimated bottom-up and top-down functions are qualitatively consistent with those from large-eddy simulation (LES) results and theoretical expectations. Newly fitted gradient functions are proposed based on observations for this forested site. The integrated bottom-up function over the lowest 4% of the CBL depth estimated from the tower data is about five times larger than that from LES results for a ‘with-canopy’ case, and is smaller than that from LES results for a ‘no-canopy’ case by a factor of 0.7. We discuss the uncertainty in the evaluated gradient functions due to stability, wind direction, and uncertainty in the entrainment flux and show that while all of these have a significant impact on the gradient functions, none can explain the differences between the modelled and observed functions. The effects of canopy features and atmospheric stability may need to be considered in the gradient function relations.  相似文献   

5.
In this paper, the model output machine learning (MOML) method is proposed for simulating weather consultation, which can improve the forecast results of numerical weather prediction (NWP). During weather consultation, the forecasters obtain the final results by combining the observations with the NWP results and giving opinions based on their experience. It is obvious that using a suitable post-processing algorithm for simulating weather consultation is an interesting and important topic. MOML is a post-processing method based on machine learning, which matches NWP forecasts against observations through a regression function. By adopting different feature engineering of datasets and training periods, the observational and model data can be processed into the corresponding training set and test set. The MOML regression function uses an existing machine learning algorithm with the processed dataset to revise the output of NWP models combined with the observations, so as to improve the results of weather forecasts. To test the new approach for grid temperature forecasts, the 2-m surface air temperature in the Beijing area from the ECMWF model is used. MOML with different feature engineering is compared against the ECMWF model and modified model output statistics (MOS) method. MOML shows a better numerical performance than the ECMWF model and MOS, especially for winter. The results of MOML with a linear algorithm, running training period, and dataset using spatial interpolation ideas, are better than others when the forecast time is within a few days. The results of MOML with the Random Forest algorithm, year-round training period, and dataset containing surrounding gridpoint information, are better when the forecast time is longer.  相似文献   

6.
The Role of Shear in the Morning Transition Boundary Layer   总被引:1,自引:1,他引:0  
We use large-eddy simulation (LES) to better define the early stages of the morning transition boundary layer. Previous LES studies relating to the morning transition boundary layer focus on the role of the entraining convective boundary layer (CBL). By using a combination of different domain sizes and grid lengths, the full evolution from the stable boundary layer (SBL) to the CBL is modelled here. In the early stages of the morning transition the boundary layer is shown to be a combination of a shallow mixed layer capped by a significant shear driven stable boundary layer (the so-called mixed CBL–SBL state). The mixed CBL–SBL state is the key to understanding the sensitivity to shear. Turbulent kinetic energy budgets also indicate that it is shear driven. The negative flux from the mixed CBL–SBL state extends much further above the minimum than is typically found for the CBL later in the day, and the depth of penetration scales as w m /N i , where w m is the combined friction and convective velocity scale and N i the static stability at the inversion top.  相似文献   

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

8.
Turbulence in a non-strongly stably stratified large-eddy simulation (LES) case is studied through probability density functions (PDFs) to obtain additional information than that provided by classical LES averages. The PDFs are computed for one hour within the steady-state regime at three different levels: near the surface, in the middle and at the top of the boundary layer, for the wind components and the temperature. The physical significance of these PDFs from LES is discussed and they are compared to those obtained from observations. The analysis of the eddy structures within the stably stratified boundary layer is made through the combined study of the fields, the spectra and the statistical moments obtained from the PDFs and joint PDFs. The homogeneity of the fields is inspected through a comparison of the ensemble to the temporal and the spatial PDFs, showing that the ergodic theorem is not fulfilled. To this end, the sensitivity of the PDF moments to the LES resolution is explored.  相似文献   

9.
Turbulence, Radiation and fog in Dutch Stable Boundary Layers   总被引:5,自引:1,他引:4  
The effect of longwave radiation on the structure the clear stable boundary layer (SBL) is examined. Special emphasis is given to radiative cooling near the surface and the top of the boundary layer and its impact on the heat flux profile. Further, the formation, growth and dissipation of fog in the SBL are studied both from observations and from a one-dimensional ensemble averaged turbulence closure model. The model is compared with detailed observations that were made for both a shallow (about 30 m) radiation fog and a deep (about 200 m) fog layer at the 200-m tower at Cabauw in the Netherlands. The model describes adequately the most important mechanisms occurring during the fog evolution. In this study special attention is given to the parameterization of the vegetation, which is important for a good representation of the (minimum) air temperature. The influence of turbulence transport, longwave radiative cooling and gravitational droplet settling on the fog evolution is described. The study demonstrates the need for more accurate measurements of turbulence quantities, especially the master length scale, in a variety of SBLs.  相似文献   

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

11.
云天地表总辐射和净辐射瞬时值的计算方法   总被引:1,自引:0,他引:1  
为减少计算机时,满足实时预报要求,全球数值预报模式中的辐射计算频率通常设定为三小时。这样处理会大大减少计算量,但也同时导致较大辐射日变化偏差,并影响模式对地面能量平衡,对流及降水的模拟。为改进这一缺陷,我们开发了一种辐射快速计算方案,可用于计算瞬时地面太阳总辐射和净辐射,使到达地面的太阳辐射计算可与模式积分同步进行,从而改善地面太阳辐射日变化模拟。本文介绍云天的计算方法。该方案所用的输入变量均为预报模式或卫星观测所能提供的量。结果表明:该方案既可用于数值预报模式也可利用观测资料独立计算地面太阳辐射。经与美国能源部大气辐射观测资料检验,该方案的精度很高,地面总辐射瞬时值的平均计算误差小于7%。  相似文献   

12.
Large-eddy simulation (LES) of a stable atmospheric boundary layer is performed using recently developed dynamic subgrid-scale (SGS) models. These models not only calculate the Smagorinsky coefficient and SGS Prandtl number dynamically based on the smallest resolved motions in the flow, they also allow for scale dependence of those coefficients. This dynamic calculation requires statistical averaging for numerical stability. Here, we evaluate three commonly used averaging schemes in stable atmospheric boundary-layer simulations: averaging over horizontal planes, over adjacent grid points, and following fluid particle trajectories. Particular attention is focused on assessing the effect of the different averaging methods on resolved flow statistics and SGS model coefficients. Our results indicate that averaging schemes that allow the coefficients to fluctuate locally give results that are in better agreement with boundary-layer similarity theory and previous LES studies. Even among models that are local, the averaging method is found to affect model coefficient probability density function distributions and turbulent spectra of the resolved velocity and temperature fields. Overall, averaging along fluid pathlines is found to produce the best combination of self consistent model coefficients, first- and second-order flow statistics and insensitivity to grid resolution.  相似文献   

13.
One-dimensional turbulence (ODT) is a single-column simulation in which vertical motions are represented by an unsteady advective process, rather than their customary representation by a diffusive process. No space or time averaging of mesh-resolved motions is invoked. Molecular-transport scales can be resolved in ODT simulations of laboratory-scale flows, but this resolution of these scales is prohibitively expensive in ODT simulations of the atmospheric boundary layer (ABL), except possibly in small subregions of a non-uniform mesh.Here, two methods for ODT simulation of the ABL on uniform meshes are described and applied to the GABLS (GEWEX Atmospheric Boundary Layer Study; GEWEX is the Global Energy and Water Cycle Experiment) stable boundary-layer intercomparison case. One method involves resolution of the roughness scale using a fixed eddy viscosity to represent subgrid motions. The other method, which is implemented at lower spatial resolution, involves a variable eddy viscosity determined by the local mesh-resolved flow, as in multi-dimensional large-eddy simulation (LES). When run at typical LES resolution, it reproduces some of the key high-resolution results, but its fidelity is lower in some important respects. It is concluded that a more elaborate empirically based representation of the subgrid physics, closely analogous to closures currently employed in LES of the ABL, might improve its performance substantially, yielding a cost-effective ABL simulation tool. Prospects for further application of ODT to the ABL, including possible use of ODT as a near-surface subgrid closure framework for general circulation modeling, are assessed.  相似文献   

14.
陈德辉 《气象》1994,20(9):9-13
作者介绍了近年来在数值天气预报领域的一个新动向-变网格一体化模式的进展。这种模式可以替代目前业务数值天气预报的有限区模式和全球谱模式,降低业务数值天气预报的计算成本。网格距离的可变性,以及采用半隐式-半拉格朗日式时间差分方案是该模式的主要特点。法国的变网格一体化模式自1992年已投入业务使用,其可行性已得到了预报检验的初步证实。  相似文献   

15.
数值天气预报业务模式现状与展望   总被引:36,自引:5,他引:31  
陈德辉  薛纪善 《气象学报》2004,62(5):623-633
20 0 4年是数值天气预报理论提出 10 0周年 ,同时也是数值天气预报业务化应用 5 0周年。经过百年的发展历程 ,数值天气预报学科有了飞跃的发展。特别是最近 10多年来 ,大气科学以及地球科学的研究进展 ,高速度、大容量的巨型计算机及网络系统的快速发展 ,更加快了数值天气预报的发展步伐。文中从模式动力框架、物理过程参数化、模式程序软件等方面对数值天气预报业务模式现状进行了简要综述 ,对存在的问题进行了探讨 ,并对数值天气预报业务模式的未来发展作了展望。当前数值天气预报业务模式发展的特点有 :(1) 2 0世纪 90年代中期以来 ,各国的全球和区域模式水平和垂直分辨率都有明显提高 ,且模式物理过程也同步进行改进 ;数值预报业务模式已进入了大规模并行计算的阶段 ;(2 )主要发达国家和中国都正在致力于研发各自的新业务数值预报模式———非静力 (多尺度 )一体化模式或非静力中尺度模式 ,部分国家的新一代天气 气候一体化数值模式已业务运行 ;(3)业务数值预报模式正在朝着不断完善的方向发展。随着模式分辨率的提高 ,云物理过程、陆面过程和湍流过程、考虑坡度 -坡向因子的辐射过程等在模式中的参数化方案 ,以及模式垂直坐标的选择越来越受重视 ,这些物理过程的描述成为业务数值模式改进的重点和  相似文献   

16.
Large-eddy simulation (LES) is a well-established numerical technique, resolving the most energetic turbulent fluctuations in the planetary boundary layer. By averaging these fluctuations, high-quality profiles of mean quantities and turbulence statistics can be obtained in experiments with well-defined initial and boundary conditions. Hence, LES data can be beneficial for assessment and optimisation of turbulence closure schemes. A database of 80 LES runs (DATABASE64) for neutral and stably stratified planetary boundary layers (PBLs) is applied in this study to optimize first-order turbulence closure (FOC). Approximations for the mixing length scale and stability correction functions have been made to minimise a relative root-mean-square error over the entire database. New stability functions have correct asymptotes describing regimes of strong and weak mixing found in theoretical approaches, atmospheric observations and LES. The correct asymptotes exclude the need for a critical Richardson number in the FOC formulation. Further, we analysed the FOC quality as functions of the integral PBL stability and the vertical model resolution. We show that the FOC is never perfect because the turbulence in the upper half of the PBL is not generated by the local vertical gradients. Accordingly, the parameterised and LES-based fluxes decorrelate in the upper PBL. With this imperfection in mind, we show that there is no systematic quality deterioration of the FOC in the strongly stable PBL provided that the vertical model resolution is better than 10 levels within the PBL. In agreement with previous studies, we found that the quality improves slowly with the vertical resolution refinement, though it is generally wise not to overstretch the mesh in the lowest 500 m of the atmosphere where the observed, simulated and theoretically predicted stably stratified PBL is mostly located. The submission to a special issue of the “Boundary-Layer Meteorology” devoted to the NATO advanced research workshop “Atmospheric Boundary Layers: Modelling and Applications for Environmental Security”.  相似文献   

17.
In the stable atmospheric boundary layer (SBL), the interaction of processes in numerical weather prediction (NWP) models seems to test their representation more stringently than their separate 'validation.Some SBL schemes derived from micrometorological research seem to allow a 'decoupling behaviour when implemented in NWP. That is, turbulence dies out from the ground upwards. Such 'decoupling of the surface from atmospheric fluxes can permit dramatic and possibly unrealistic falls in surface temperature. This study traces the mechanism of model decoupling, asks whether this behaviour is in any sense correct and considers the implications.It is shown that decoupling can occur in idealized single-column models, originating from an unstable boundary-mode. This behaviour can depend critically on parameters such as surface roughness and soil thermal diffusivity as well as turbulence. But the turbulence dependence arises through the response of the boundary layer as a whole, and not just the surface-layer scheme. Such decoupling arises from the 'physics, rather than the finite-difference schemes, and appears to occur sometimes in the real atmosphere.  相似文献   

18.
边界层参数化方案在“灰色区域”尺度下的适用性评估   总被引:2,自引:0,他引:2  
随着数值预报模式分辨率的提高,当模式网格距与含能湍涡的长度尺度相当时,模式动力过程可解析一部分湍流运动,而剩余的湍流运动仍需参数化,此时便产生了湍流参数化的“灰色区域”问题。对传统的PBL(Planetary Boundary Layer)方案在“灰色区域”下的适用性评估,是改进PBL方案以使其能够适应分辨率变化的前提和基础。本研究基于干对流边界层的大涡模拟试验,比较了WRF(Weather Research and Forecast Model)模式中四种常用的边界层参数化方案[YSU(Yonsei University)、MYJ(Mellor-Yamada-Janjic)、MYNN2.5(Mellor-Yamada-Nakanishi-Niino Level 2.5)、MYNN3)]在“灰色区域”尺度下的表现。研究表明,混合层内总热通量对所使用的参数化方案和水平分辨率均不敏感。不同参数化方案中次网格与网格通量的比例表现出对水平网格距不同的依赖性。局地PBL方案(MYJ、MYNN2.5)在混合层内的平均位温随网格距减小而增大,次网格通量随网格距减小而减小,较参考湍流场对次网格通量有所低估。YSU方案的非局地项几乎不随水平格距改变而变化,对次网格通量的表征并未表现出较强的分辨率依赖性,且过强的非局地次网格输送使混合层内温度层结呈弱稳定,抑制了可分辨湍流输送,不易于激发次级环流。MYNN3方案的非局地次网格通量(负梯度输送项)随网格距减小而减小,使其对次网格通量的表征具有较好的分辨率依赖性。PBL方案在“灰色区域”尺度下的适用性与具体分辨率有关。以分辨率500 m为例,四种PBL方案中不存在一种最佳方案,能对边界层的热力结构和湍流统计特征均有准确的描述。  相似文献   

19.
We present the main results from the second model intercomparison within the GEWEX (Global Energy and Water cycle EXperiment) Atmospheric Boundary Layer Study (GABLS). The target is to examine the diurnal cycle over land in today??s numerical weather prediction and climate models for operational and research purposes. The set-up of the case is based on observations taken during the Cooperative Atmosphere-Surface Exchange Study-1999 (CASES-99), which was held in Kansas, USA in the early autumn with a strong diurnal cycle with no clouds present. The models are forced with a constant geostrophic wind, prescribed surface temperature and large-scale divergence. Results from 30 different model simulations and one large-eddy simulation (LES) are analyzed and compared with observations. Even though the surface temperature is prescribed, the models give variable near-surface air temperatures. This, in turn, gives rise to differences in low-level stability affecting the turbulence and the turbulent heat fluxes. The increase in modelled upward sensible heat flux during the morning transition is typically too weak and the growth of the convective boundary layer before noon is too slow. This is related to weak modelled near-surface winds during the morning hours. The agreement between the models, the LES and observations is the best during the late afternoon. From this intercomparison study, we find that modelling the diurnal cycle is still a big challenge. For the convective part of the diurnal cycle, some of the first-order schemes perform somewhat better while the turbulent kinetic energy (TKE) schemes tend to be slightly better during nighttime conditions. Finer vertical resolution tends to improve results to some extent, but is certainly not the solution to all the deficiencies identified.  相似文献   

20.
We have analyzed measurements of vertical velocity w statistics with the NOAA high resolution Doppler lidar (HRDL) from about 390 m above the surface to the top of the convective boundary layer (CBL) over a relatively flat and uniform agricultural surface during the Lidars-in-Flat-Terrain (LIFT) experiment in 1996. The temporal resolution of the zenith-pointing lidar was about 1 s, and the range-gate resolution about 30 m. Vertical cross-sections of w were used to calculate second- to fourth-moment statistics of w as a function of height throughout most of the CBL. We compare the results with large-eddy simulations (LES) of the CBL and with in situ aircraft measurements. A major cause of the observed case-to-case variability in the vertical profiles of the higher moments is differences in stability. For example, for the most convective cases, the skewness from both LES and observations changes more with height than for cases with more shear, with the observations changing more with stability than the LES. We also found a decrease in skewness, particularly in the upper part of the CBL, with an increase in LES grid resolution.  相似文献   

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

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