首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 43 毫秒
1.
In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ -1. This mesoregion separates the inner region k ~ k s(where k s-1 and is the Kolmogorov length) and the outer region k k e(where k e -1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k -5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k -5/3-behavior in two regions (-1 k k s) and (k e k -1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k -5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k m-1.  相似文献   

2.
This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.  相似文献   

3.
A pair of parallel cold wires separated in either the vertical or lateral direction was used to obtain the three components x, y, z of the temperature derivative in the streamwise, lateral and vertical directions, respectively. The average absolute skewness values of x and z are nonzero and approximately equal, while the skewness of y is approximately zero. These results appear to be consistent with the presence of a large, three-dimensional organised structure in the surface layer. There is an apparent low-frequency contamination in the spectral density of y and z due mainly to small errors in estimating the sensitivity of the cold wires. The temperature derivatives were high-pass filtered, the filter being set to remove possible contributions from the large structure and to minimise low-frequency sensitivity contamination. The filtered rms ratios \~x/\~y and \~x/\~z were in the range 0.7 to 0.9, a result in qualitative agreement with that obtained in the laboratory boundary layer by Sreenivasan et al. (1977). The skewness of filtered x or z is negligible, consistent with local isotropy of small-scale temperature fluctuations and in support of the high wavenumber spectral isotropy discussed in Antonia and Chambers (1978).  相似文献   

4.
E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c 2 /e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c 2 = 1.92, = 1.3, e = 1). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and e to 1.1 and 1.6, respectively, giving, with c 2 = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.  相似文献   

5.
Summary The integral aerosol optical depths (k ) at the hour of 08:20 Local Standard Time (LST), are compared with those calculated previously at 11:20 and 14:20 LST, for clear days during summer in Athens over the period 1962–1988. The mean values at 08:20 LST were consistently lower than the values at 11:20 and 14:20 LST. The influence of the vertical wind profile on the values ofk was also investigated. A comparison was made of the wind profiles at 02:00 and 14:00 LST, for days in which the 11:20 and 14:20 LST values ofk were 0.200 andk 0.350, respectively. The corresponding bulk wind shear s was also found for the period 1980–1988. The most significant results occurred with the first category of days. The resultant wind velocities from the surface to the 900 hPa level, in each hour were higher by 2–4 m·s–1 with respect to the corresponding values for the second category. At 02:00 LST the bulk wind shear showed a considerable difference (1.8) between the two categories of days in the surface to 700 hPa layer at 02:00 LST. Finally, the associated weather conditions that appear to initiate a period of low values ofk (k 0.200) at 11:20 and 14:20 LST were examined for the period 1980–1988. Fifteen such cases were identified and it was found that they all occurred after the passage of weak cold fronts.With 6 Figures  相似文献   

6.
Summary The effect of the Alpine orography on prototype cold fronts approaching from the west is investigated by three-dimensional numerical model simulations. The numerical experiments cover a range of parameter constellations which govern the prefrontal environment of the front. Especially, the appearance and intensity of prefrontal northern Alpine foehn varies from case to case.The behaviour of a cold front north of the Alps depends much on the prefrontal condition it encounters. It is found that prefrontal foehn can either accelerate or retard the approaching front.An important feature is the pressure depression along the northern Alpine rim that results from the southerly foehn flow. In cases where this depression compensates the eastward directed pressure gradient associated with the largescale flow, the front tends to accelerate and the foehn breaks down as soon as the front passes. In contrast, the foehn prevents the front from a rapid eastward propagation if it is connected with a strong southerly wind component.No-foehn experiments are performed for comparison, where either the mountains are removed, or the static stability is set to neutral. Also shown are effects of different crossfrontal temperature contrasts.List of Symbols c F propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - z vertical grid spacing (cartesian system) - cross-frontal potential temperature difference - i potential temperature step at an inversion - E turbulent kinetic energy - f Coriolis parameter - FGP frontogenesis parameter (see section 2.2) - g gravity acceleration (g=9.81 m s–2) - vertical gradient of potential temperature - h terrain elevation (above MSL) - h i height of an inversion (h i =1000 m MSL) - H height of model lid (H=9000 m MSL) - K M exchange coefficient of momentum - K H exchange coefficient of heat and moisture - longitude - N Brunt-Väisäla-frequency - p pressure - Exner function (=T/) - latitude - q v specific humidity - R d gas constant of dry air (R d =287.06 J kg–1 K–1) - density of dry air - t time - T temperature - potential temperature - TFP thermal front parameter (see section 2.2) - u, v, w cartesian wind components - u g ,v g geostrophic wind components - horizontal wind vector - x, y, z cartesian coordinates Abbreviations GND (above) ground level - MSL (above) mean sea level - UTC universal time coordinated With 20 Figures  相似文献   

7.
Analytical solutions for the Ekman layer   总被引：1，自引：0，他引：1  
The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U *z (0.4,U *=friction velocity,z=elevation) to the geostrophic asymptoteK GU *d forzd is solved through an expansion in fd/U *1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU *d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U *d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U * (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U * and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U *d[1–exp(–z/d)] are developed in appendices.  相似文献   

8.
Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.  相似文献   

9.
Summary The standard equations for the theory of atmospheric tides are solved here by an integral representation on the continuous spectrum of free oscillations. The model profile of back-ground temperature is that of the U.S. Standard Atmosphere in the lower and middle atmosphere, and in the lower thermosphere, above which an isothermal top extends to arbitrarily great heights. The top is warm enough to bring both the Lamb and the Pekeris modes into the continuous spectrum.Computations are made for semidiurnal lunar tidal pressure at sea level at the equator, and the contributions are partitioned according to vertical as well as horizontal structure. Almost all the response is taken up by the Lamb and Pekeris modes of the slowest westward-propagating gravity wave. At sea level, the Lamb-mode response is direct and is relatively insensitive to details of the temperature profile. The Pekeris mode at sea level has an indirect response-in competition with the Lamb mode-and, as has been known since the time of its discovery, it is quite sensitive to the temperature profile, in particular to stratopause temperature. In the standard atmosphere the Lamb mode contributes about +0.078 mb to tidal surface pressure at the equator and the Pekeris mode about –0.048 mb.The aim of this investigation is to illustrate some consequences of representing the tide in terms of the structures of free oscillations. To simplify that task as much as possible, all modifying influences were omitted, such as background wind and ocean or earth tide. Perhaps the main defect of this paper's implementation of the free-oscillation spectrum is that, in contrast to the conventional expansion in the structures of forced oscillations, it does not include dissipation, either implicity or explicity, and thus does not satisfy causality. Dissipation could be added implicity by means of an impedance condition, for example, which would cause up-going energy flux to exceed downgoing flux at the base of the isothermal top layer. To achieve complete causality, however, the dissipation must be modeled explicity. Nevertheless, since the Lamb and Pekeris modes are strongly trapped in the lower and middle atmosphere, where dissipation is rather weak (except possibly in the surface boundary layer), more realistic modeling is not likely to change the broad features of the present results.Symbols a earth's mean radius; expansion coefficient in (5.3) - b recursion variable in (7.4); proximity to resonance in (9.2) - c sound speed in (2.2); specific heatc p in (2.2) - f Coriolis parameter 2sin in (2.2) - g standard surface gravity - h equivalent depth - i ; discretization index in (7.3) - j index for horizontal structure - k index for horizontal structure; upward unit vectork in (2.2) - m wave number in longitude - n spherical-harmonic degree; number of grid layers in a model layer - p tidal pressure perturbation; background pressurep 0 - q heating function (energy per mass per time) - r tidal state vector in (2.1) - s tidal entropy perturbation; background entropys 0 - t time - u tidal horizontal velocityu - w tidal vertical component of velocity - x excitation vector defined in (2.3); vertical coordinate lnp */p 0 [except in (3.8), where it is lnp /p 0] - y vertical-structure function in (7.1) - z geopotential height - A constant defined in (6.2) - C spherical-harmonic expansion coefficient in (3.6) - D vertical cross section defined in (5.6) and (5.9) - E eigenstate vector - F vertical-structure function for eigenstate pressure in (3.2) [re-defined with WKB scaling in (7.2)] - G vertical-structure function for eigenstate vertical velocity in (3.2) [re-defined with WKB scaling in (7.2)] - H pressure-scale height - I mode intensity defined in (8.1) - K quadratic form defined in (4.4) - L quadratic form defined in (4.4); horizontal-structure magnification factor defined in (5.11) - M vertical-structure magnification factor defined in (4.6) - P eigenstate pressure in (3.2); tidal pressure in (6.2) - R tidal state vector in (5.1) - S eigenstate entropy in (3.2); spherical surface area, in differential dS - T background molecular-scale (NOAA, 1976) absolute temperatureT 0 - U eigenstate horizontal velocityU in (3.2); coefficient in (7.3) - V horizontal-structure functionV for eigenstate horizontal velocity in (3.2); recursion variable in (7.3) - W eigenstate vertical velocity in (3.2) - X excitation vector in (5.1) - Y surface spherical harmonic in (3.7) - Z Hough function defined in (3.6) - +dH/dz - (1––)/2 - Kronecker delta; Dirac delta; correction operator in (7.6) - equilibrium tide elevation - (square-root of Hough-function eigenvalue) - ratio of specific gas constant to specific heat for air=2/7 - longitude - - - background density 0 - eigenstate frequency in (3.1) - proxy for heating functionq =c P/t - latitude - tide frequency - operator for the limitz - horizontal-structure function for eigenstate pressure in (3.2) - Hough function defined in (6.2) - earth's rotation speed - horizontal gradient operator - ()0 background variable - ()* surface value of background variable - () value at base of isothermal top layer - Õ state vector with zerow-component - , energy product defined in (2.4) - | | energy norm - ()* complex conjugate With 10 Figures  相似文献   

10.
Daily mean values of the Priestley-Taylor coefficient, ¯, are derived from a simple model of the daily growth of a convective boundary layer. For a particular control set of driving environmental variables, ¯ is related to the prescribed bulk surface resistance, rS by 1/¯ = 1/0 + mrS for parameters 0 and m. The dependence of the parameters 0 and m on weather is explored and a potential use of this linear relation to provide information about regional values of rS is indicated.  相似文献   

11.
Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO2 concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO2, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO2. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= cw for CO2, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux c w. Thus, when the net flux of CO2 is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., c w = – c w. Corrections for the mean convective flux are particularly significant for CO2 because cw and c w are often of similar magnitude. The correct measurement of the net CO2 flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.  相似文献   

12.
Drag and drag partition on rough surfaces   总被引：13，自引：0，他引：13  
An analytic treatment of drag and drag partition on rough surfaces is given. The aims are to provide simple predictive expressions for practical applications, and to rationalize existing laboratory and atmospheric data into a single framework. Using dimensional analysis and two physical hypotheses, theoretical predictions are developed for total stress (described by the square root of the canopy drag coefficient), stress partition (described by the ratio S/ of the stress s on the underlying ground surface to total stress ), zero-plane displacement and roughness length. The stress partition prediction is the simple equation S/= 1/(1+), where = CRCS the ratio of element and surface drag coefficients. This prediction agrees very well with data and is free of adjustable constants. Other predictions also agree well with a range of laboratory and atmospheric data.  相似文献   

13.
The carbon isotopic ratio of atmospheric carbon dioxide at Tsukuba,Japan   总被引：1，自引：0，他引：1  
To find out the secular and seasonal trends of the 13C value and CO2 concentration in the surface air and the determination of the 13C in the atmospheric CO2 collected at Tsukuba Science City was carried out during the period from July 1981 to October 1983. The monthly average of the 13C value of CO2 in the surface air collected at 1400 LMT ranged from -7.52 to \s-8.45 with an average of -7.96±0.25 and the CO2 concentration in the air varied from 334.5 l 1-1 to 359 l 1-1 with an average of 347.2±6.3 l 1-1. The 13C value is high in summer and low in winter and is negatively correlated with the CO2 concentration. In general, the relationship between the 13C and the CO2 concentration is explainable by a simple mixing model of two different constant carbon isotopic species but the relationship does not always follow the model. The correlation between the 13C value and the CO2 concentration is low during the plant growth season and high at other times. The observed negative deviation of the 13C value from the simple mixing model in the plant growth season is partly due to the isotopic fractionation process which takes place in the land biota.  相似文献   

14.
  相似文献   

15.
The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and  相似文献   

16.
Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /PE irrotational are zero-spin models, having d z=0, is proven  相似文献   

17.
Effect of finite sampling on atmospheric spectra   总被引：2，自引：0，他引：2  
The effect of a finite averaging time on variances is well known, but its effect on power spectra is less clearly understood. We present numerical solutions for the spectral distortion arising from sampling over a finite time interval T and show that the commonly used filter function (1 – sinc2f T), valid for variances, is a reasonable approximation for power spectra only when T 10 m , where f is the cyclic frequency, and m is the dominant time scale of the process. Our results exhibit an increasingly steeper low-frequency roll-off as T decreases relative to m , indicating that the measured spectrum is subject to a greater suppression of the lower frequencies (f > 1/T) than predicted by (1 – sinc2f T). This suppression is, in a sense, compensated by an overestimation of spectral estimates in the frequency range f 1/T.  相似文献   

18.
The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O2 of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O2 versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(Hdecomposition–HO 2 reaction), between these two reactions pathways. For (H)[22–0.5(HO 2 reaction)], alkoxy radical decomposition dominates over the reaction with O2 at room temperature and atmospheric pressure of air, while for (H)[25-0.5(HO 2 reaction)], the O2 reaction dominates over decomposition (where the units of H are in kcal mol–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.  相似文献   

19.
Analyses of concentration fluctuation (C) spectra from boundary-layer smoke plume experiments at six separate locations show that the spectra from these experiments generally exhibit an inertial subrange at high frequencies with a slope of -5/3 and indicate peak energy at a time period of about 50 to 100 s. These periods of peak energy are a factor of two to five less than those for the peak of the wind speed fluctuation (u or v) spectra. A general spectral formula fits normalized spectra from the U.S. and Australia, where the frequency, n, is made dimensionless by multiplying by the plume dispersion parameter, y , and dividing by the wind speed, u. Peak energy occurs at a dimensionless frequency of n y/u equal to about 0.15. The Kolmogorov constant in the inertial subrange is estimated from a set of averaged spectra. Cross-spectra indicate little relation between concentration and wind fluctuations. However, most of the correlation that exists is due to periods larger than about 10 or 20 s.  相似文献   

20.
Experimental evidence indicates that the diurnal behaviour of the fluxes of heat into the ground and into the atmosphere versus the net flux of radiation can be modelled by closed curves, the hourly values folowing one another in either a clockwise or counter clockfashion. A general formulation to express the different heat fluxes as a function of net radiation is proposed. This relationship between the different heat fluxes and can be expressed as a sum of three terms: the first indicates a direct proportionality to , the second gives the deviation from linearity and depends on /t, and the third gives the value of the flux when = 0. The formulae are then expressed as a function of time and the ratios between the heat fluxes and are evaluated. A comparison with the approximations generally used shows that the latter may be considered as particular cases of the more general equations proposed here.  相似文献   

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

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