首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Mesoscale models using a non-local K-scheme for parameterization of boundary-layer processes require an estimate of the planetary boundary layer (PBL) height z i at all times. In this paper, two-dimensional sea-breeze experiments are carried out to evaluate three different formulations for the advective contribution in the z i prognostic equation of Deardorff (1974).Poor representation of the thermal internal boundary layer in the sea breeze is obtained when z i is advected by the wind at level z i . However, significantly better results are produced if the mean PBL wind is used for the advecting velocity, or if z i is determined simply by checking for the first sufficiently stable layer above the ground.A Lagrangian particle model is used to demonstrate the effect of each formulation on plume dispersion by the sea breeze.  相似文献   

2.
By integrating the Fokker-Planck equation corresponding to a Lagrangian stochastic trajectory model, which is consitent with the selection criterion of Thomson (1987), an analytical solution is given for the joint probability density functionp(xi, ui, t) for the position (x i) and velocity (u i) at timet of a neutral particle released into linearly-sheared, homogeneous turbulence. The solution is compared with dispersion experiments conforming to the restrictions of the model and with a shortrange experiment performed in highly inhomogeneous turbulence within and above a model crop canopy. When the turbulence intensity, wind shear and covariance are strong, the present solution is better than simpler solutions (Taylor, 1921; Durbin, 1983) and as good as any numerical Lagrangian stochastic model yet reported.  相似文献   

3.
Convection in a quasi-steady, cloud-free, shear-free atmospheric boundary layer is investigated based on a large-eddy simulation model. The performed tests indicate that the characteristic (peak) values of statistical moments at the top of the mixed layer are proportional to the interfacial scales (from gradients of scalars in the interfacial layer). Based on this finding a parameterization is proposed for profiles of scalar variances. The parameterization employs two, semi-empirical similarity functions Fm(z/zi) andFi(z/zi), multiplied by a combination of the mixed-layer scales and the interfacial scales.  相似文献   

4.
Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The turbulent fields are obtained directly by solving the three-dimensional time-dependent Navier-Stokes equations (using the Boussinesq approximation to account for buoyancy effects), and no turbulence model is needed. Two levels of heating are considered, one quite vigorous, the other more moderate. Statistics for the vigorously heated case are found to agree reasonably well with laboratory, field, and large-eddy simulation results, when Deardorff's mixed-layer scaling is used. No indication of large-scale longitudinal roll cells is found in this convection-dominated flow, for which the inversion height to Obukhov length scale ratio –z i /L *=26. However, when heating is more moderate (so that –z i /L *=2), evidence of coherent rolls is present. About 10% of the total turbulent kinetic energy and turbulent heat flux, and 20% of the Reynolds shear stress, are estimated to be a direct consequence of the observed cells.  相似文献   

5.
A simple time-dependent one-dimensional model of the planetary boundary layer (PBL) is described and used to examine the degree to which model design decisions affect model output variables. The model's sensitivity to changes in the environmental conditions is also explored. Averages of the surface fluxes, near-ground wind speeds and other PBL properties from 48 h simulations are compared to control runs. The model-calculated surface fluxes are most sensitive, in decreasing order of importance, to the vertical grid spacing, the form of closure between the surface temperature and the atmosphere, the use of vertical diffusivity smoothing, the choice of maximum time step and choice of turbulence closure scheme. These fluxes are relatively insensitive to mixing-length scaling or choice of implicit time step weighting factor. Sensitivity to changes in soil type exceeds any of the design criteria tested. The modeled fluxes are moderately sensitive to small variations in the horizontal pressure gradient, to unsteadiness in the geostrophic wind and to variations in surface roughness. They are relatively insensitive to uncertainties in local vertical velocities and small (25%) variations applied separately to soil thermal diffusivity or heat capacity. The sensitivity of the average PBL depth (Z i ) to model and environmental changes are similar to those of surface fluxes except thatZ i is more sensitive to changes in mixing length, albedo and imposed vertical velocity then are the surface fluxes.  相似文献   

6.
Forced convection in a quasi-steady atmospheric boundary layer is investigated based on a large-eddy simulation (LES) model. The performed simulations show that in the upper portion of the mixed layer the dimensionless (in terms of mixed layer scales) vertical gradients of temperature, humidity, and wind velocity depend on the dimensionless height z/z i and the Reech number Rn. The peak values of variances and covariances at the top of the mixed layer, scaled in terms of the interfacial scales, are functions of the interfacial Richardson number Ri. As a result expressions for the entrainment rates, in the case when the interfacial layer has a finite depth, and a condition for the presence of moistening or drying regimes in the mixed layer, are derived. Profiles of dimensionless scalar moments in the mixed layer are proposed to be expressed in terms of two empirical similarity functions F m and F i , dependent on dimensionless height z/z i , and the interfacial Richardson number Ri. The obtained similarity expressions adequately approximate the LES profiles of scalar statistics, and properly represent the impact of stability, shear, and entrainment. They are also consistent with the parameterization proposed for free convection in the first part of this paper.  相似文献   

7.
This paper examines the practical importance of stability, baroclinicity, and acceleration in the bulk ABL similarity formulations, in light of the random errors inherent in field measurements. This is done by propagating the measurement uncertainties through a theoretical model for the bulk ABL similarity functionsA 0 andB 0, under a range of assumed (but always unstable) conditions. It is shown that random measurement errors and acceleration effects may overwhelm most effects of baroclinicity and stability, once conditions are at least slightly unstable. Because of this, it is hard to discern a clear functional dependence ofA 0 andB 0 on degree of instability. Thus, for a given value ofh i/z0, whereh i is the inversion height andz 0 is the surface roughness, the geostrophic drag coefficient, which depends onA 0 andB 0, and weakly onh i/z0, may also be taken to be nearly independent of degree of instability.  相似文献   

8.
本文在GRAPES_TMM（Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model）——中国南海台风模式版（面向南海和东南亚）中发展和引进了KA95（Kim and Arakawa,1995）地形重力波拖曳参数化方案（GWDO）,并对2012年主要的9个登陆台风进行了试验对比研究,考察了不同标准Richardson数（Ric）的GWDO试验对台风路径和强度预报的影响。结果表明,在引入地形重力波拖曳参数化过程后,模式对台风登陆时路径和强度的预报能力均要有提高,对台风预报时长越长,GWDO的影响也更为显著。对双台风“SAOLA”和“DAMREY”试验结果表明,GWDO对台风外围距台风中心150 km的对流层中下层风速减弱较为明显,减弱了GRAPES区域模式对台风强度预报偏强的现象,对台风强度长时间预报改善更为明显。不同标准Ric对重力波拖曳力的计算较为敏感,当Ric取1.0时,动能迅速的在低层被频散,能量无法有效地上传;Ric取0.25时,大部分的能量在中高层被频散。总的来说,Ric取0.75时对台风路径和强度预报改进更为显著,其结果可为业务预报提供指导意义。  相似文献   

9.
It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u */w *, where u * and w * are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z i , where z i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.  相似文献   

10.
A Random Displacement Model (RDM) and a Langevin Equation Model (LEM) are used to simulate point releases in a complex flow around a building. The flow field is generated by a three-dimensional finite element model that uses the standardk- model to parameterize the turbulence. The RDM- and LEM-calculated concentration fields are compared, with particular emphasis on the structure in regions with high turbulence and/or recirculation. RDM and LEM results are similar qualitatively, but RDM tends to predict lower concentration levels. In part this is due to the higher early-time diffusion. However, the expected convergence at later times is prevented by the interaction of the diffusion with the strongly inhomogeneous mean flow.Notation a i coefficient in the Langevin equation - b ij coefficient in the Langevin equation - C 0 the universal constant associated with the Lagrangian structure function - H building height (22.5 m) - K eddy viscosity - K k eddy viscosity used in the definition of the off-diagonal Reynolds stresses - k turbulent kinetic energy - LEM Langevin Equation Model - p 1 local unit vector in thexy-plane, orthogonal tos - p 2 local unit vector, orthogonal to boths andp 1 - RDM Random Displacement Model - s local unit vector in the streamline direction - T local decorrelation time (Lagrangian time scale) - U magnitude of the local Eulerian mean wind velocity - u s total velocity in the streamline direction - u 1 velocity component in thexy-plane, orthogonal to the streamline direction - u 2 velocity component orthogonal to bothu s andu 1 - i mean Eulerian wind velocity - W i stochastic vector-valued Wiener process - x unit vector inx-direction - y unit vector iny-direction - z unit vector inz-direction - angle between thexy-plane and the mean wind streamline - angle between the projection in thexy-plane of the streamline and thex-axis - ij the Kronecker delta function - rate of turbulence dissipation - i/ga the part ofa i that contains mean wind and turbulence gradients - ij inverse of a Reynolds stress tensor component - ij shorthand for a quantity that defines a part of i/ga - i shorthand for a quantity that defines a part of i/ga - ij Reynolds stress tensor component  相似文献   

11.
A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z 01)/ln(δ i /z 01), where z, δ i and z 01 are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.  相似文献   

12.
Our focus is the time evolution of the turbulent kinetic energy for decaying turbulence in the convective boundary layer. The theoretical model with buoyancy and inertial transfer terms has been extended by a source term due to mechanical energy and validated against large-eddy simulation data. The mechanical effects in a boundary layer of height z i at a convective surface-layer height z = 0.05z i are significant in the time evolution of the vertical component of the spectrum, i.e. they enhance the decay time scale by more than an order of magnitude. Our findings suggest that shear effects seem to feedback to eddies with smaller wavenumbers, preserving the original shape of the spectrum, and preventing the spectrum from shifting towards shorter wavelengths. This occurs in the case where thermal effects only are considered.  相似文献   

13.
Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z i. Neither the two-peak form at slight instability nor the dependence of peak position on z i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z s. The z s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (zsz)1/3eo2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where eo{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z i z s)1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z i z s)1/4 z 1/2; (4) for z/z s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.  相似文献   

14.
Aerodynamic roughness of vegetated surfaces   总被引：2，自引：0，他引：2  
Available experimental results indicate that as the density of roughness elements over a horizontally homogeneous surface is varied, the roughness length, z 0, varies in a manner that exhibits a maximum at intermediate density values. In an attempt to explain this behaviour, the available analytical solutions for the wind profile inside dense homogeneous canopies were reviewed. The review indicated that the variation of z 0 with density depends on the interrelationship between the leaf density, a, and the mixing length, l. In view of this finding, a numerical model was devised based on a simple rule for constructing mixing-length profiles in the canopy. The rule states that the actual value of l is the maximum possible under the two constraints: l l i and ¦dl/dz¦ k, where k is the von Karman constant and the intrinsic mixing length, l i, is a function of the local internal structure of the canopy. The model which ensures a smooth transition from dense to thin canopy, was used to reproduce the observed maximum of z 0. The model is also capable of handling vertically non-homogeneous canopies.  相似文献   

15.
It is proposed that the ratios of the standard deviations of the horizontal velocity components to the friction velocity in the surface layer under convective conditions depend only onz i /L wherez i is the height of the lowest inversion andL is the Monin-Obukhov length. This hypothesis is tested by using observations from several data sets over uniform surfaces and appears to fit the data well. Empirical curves are fitted to the observations which have the property that at largez i /-L, the standard deviations become proportional tow *, the convective scaling velocity.Fluctuations of vertical velocity obtained from the same experiments scale withz/L, wherez is the height above the surface, in good agreement with Monin-Obukhov theory.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

16.
The introduced mathematical model takes into account the role of the kinetic fractionation effect in a supersaturation environment at the ice surface as liquid and solid phases coexist in mixed cloud. Using the model, the temperature effect of stable isotopes in precipitation is simulated under different cooling conditions. The rate of change of δ18O against temperature in the process of wet adiabatic cooling is smaller than in the process of isobaric cooling under the same humidity. The increasing supersaturation ratio at the ice surface, Si, leads to the strengthening of the kinetic fractionation effect. The kinetic fractionation function makes the synthesis fractionation factor decreased and the change of δ18O with temperature flatted, compared with that in the equilibrium state. The simulated results show that the slope parameter b and the intercept d of the meteoric water line (MWL), δD = bδ18O+d, in wet adiabatic cooling are both greater than those in isobaric cooling. The global MWL lies between the two MWLs simulated under wet adiabatic and isobaric cooling processes, respectively. The magnitudes of b and d are directly proportional to Si. The greater the Si, the stronger the kinetic fractionation effect, and thus the greater the b and d, and vice versa. However, b and d have low sensitivity to the liquid-water contents in the cloud. Using the kinetic fractionation model, the variation of stable isotopes in precipitation at Urumqi is simulated. The simulated stable isotopic ratio vs temperature and the δD vs δ18O curves are very consistent with the actual regressions and MWL at Uruimqi, respectively.  相似文献   

17.
Meteorological data obtained from sensors mounted on a 300-m antenna mast, including wind components and temperature fluctuations together with high resolution profile measurements, were used to calculate quantities which describe the most important processes during the development and decay of ground-based inversions. The vertical distribution of air temperature during the developing phase of radiation inversions has been measured continuously and compared with an existing simple prognostic model. Free and forced convection conditions have been investigated for the inversion lifting phase by measuring the heat fluxes at the surface (H s) and at the inversion base (H i), the heating rate of the mixed layer (/z) m , the friction velocity (u *) and the entrainment rate (dh/dt). Comparisons between the directly measured entrainment rates and those calculated by models show good agreement. For nearly free convection conditions, the entrainment parameter (k = -H i/Hs) has been determined to be about 0.25.  相似文献   

18.
High-resolution water vapour measurements made by the Atmospheric Radiation Measurement (ARM) Raman lidar operated at the Southern Great Plains Climate Research Facility site near Lamont, Oklahoma, U.S.A. are presented. Using a 2-h measurement period for the convective boundary layer (CBL) on 13 September 2005, with temporal and spatial resolutions of 10 s and 75 m, respectively, spectral and autocovariance analyses of water vapour mixing ratio time series are performed. It is demonstrated that the major part of the inertial subrange was detected and that the integral scale was significantly larger than the time resolution. Consequently, the major part of the turbulent fluctuations was resolved. Different methods to retrieve noise error profiles yield consistent results and compare well with noise profiles estimated using Poisson statistics of the Raman lidar signals. Integral scale, mixing-ratio variance, skewness, and kurtosis profiles were determined including error bars with respect to statistical and sampling errors. The integral scale ranges between 70 and 130 s at the top of the CBL. Within the CBL, up to the third order, noise errors are significantly smaller than sampling errors and the absolute values of turbulent variables, respectively. The mixing-ratio variance profile rises monotonically from ≈0.07 to ≈3.7 g2 kg−2 in the entrainment zone. The skewness is nearly zero up to 0.6 z/z i , becomes −1 around 0.7–0.8 z/z i , crosses zero at about 0.95 z/z i , and reaches about 1.7 at 1.1 z/z i (here, z is the height and z i is the CBL depth). The noise errors are too large to derive fourth-order moments with sufficient accuracy. Consequently, to the best of our knowledge, the ARM Raman lidar is the first water vapour Raman lidar with demonstrated capability to retrieve profiles of turbulent variables up to the third order during daytime throughout the atmospheric CBL.  相似文献   

19.
Summary Sublayer-Stanton numbers, Bi, of heat and matter for the interfacial sublayer over aerodynamically smooth surfaces determined for forced convective conditions by elementary and numerical integration are reviewed and evaluated. The results are based on Roth’s modified Heisenberg model for the spectral energy transfer in the equilibrium range under locally isotropic conditions and the approaches of Reichardt, Elser, Deissler, van Driest, Rannie, Sheppard, and Spalding for the normalized eddy diffusivity Km/v. The results substantiate that with the exception of Sheppard’s Km/v-approach all formulations are appropriate to provide sublayer-Stanton numbers with a sufficient degree of accuracy. From a theoretical point of view the Km/v-relationships of Roth, Reichardt, van Driest, and Spalding are to be preferred, when a turbulent Prandtl number Prt = 1 is presumed. Since within the framework of mesoscale meteorological modelling numerical integration techniques would consume too much CPU-time because of the large number of near-wall grid points, a parameterization formula for the sublayer-Stanton number is proposed and evaluated. Compared to the Bi −1-results obtained by numerical integration, this kind of parameterization leads to a relative error of less than 5 percent for roughness Reynolds numbers, ηr, ranging from 30 to 600. Received January 2001 Revised November 5, 2001  相似文献   

20.
We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z m (defined by the maximum of the potential temperature gradient) and its standard deviation σ m , the mixed-layer height z i (defined by the minimum of the potential temperature flux), entrainment velocity W e, and the entrainment flux H i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.  相似文献   

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

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