首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 187 毫秒
1.
A numerical study of stably stratified flow over a three-dimensional hill is presented. Large-eddy simulation is used here to examine in detail the laboratory experimental flows described in the landmark work of Hunt and Snyder about stratified flow over a hill. The flow is linearly stratified and U/Nh is varied from 0.2 to 1.0. Here N and U are the buoyancy frequency and freestream velocity respectively, and h is the height of the hill. The Reynolds number based on the hill height is varied from 365 to 2968. The characteristic flow patterns at various values of U/Nh have been obtained and they are in good agreement with earlier theoretical and experimental results. It is shown that the flow field cannot be predicted by Drazin's theory when recirculation exists at the leeside of the hill even at UNh 1. The wake structure agrees well with a two-dimensional wake assumption when U/Nh 1 but lee waves start to influence the wake structure as U/Nh increases. The dividing-streamline heights obtained in the simulation are in accordance with experimental results and Sheppard's formula. The energy loss along the dividing streamline due to friction/turbulence approximately offsets the energy gained from pressure field. When lee waves are present, linear theory always underestimates the amplitude and overestimates the wavelength of three-dimensional lee waves. The simulated variations of drag coefficients with the parameterK (=ND/ U) are qualitatively consistent with experimental data and linear theory. Here D is the depth of the tank.  相似文献   

2.
A novel and readily applicable Structure-Activity Relationship (SAR) for predicting the barrier height Eb to decomposition by C-C scission of (substituted) alkoxy radicals is presented. Alkoxy radicals are pivotal intermediates in the atmospheric oxidation of (biogenic) volatile organic compounds, and their fate is therefore of crucial importance to the understanding of atmospheric VOC degradation mechanisms. The SAR is based on available theoretical energy barriers and validated against barriers derived from experimental data. The SAR is expressed solely in terms of the number(s) Ni of alkyl-, hydroxy- and/or oxo-substituents on the - and -carbons of the breaking bond: Eb(kcal/mol) =17.5 – 2.1 × N(alk) – 3.1 ×N(alk) – 8.0 × N,(OH) – 8.0 × N(O=) – 12 × N(O=). For barriers below 7 kcal/mol, an additional, second-order term accounts for the curvature. The SAR reproduces the available experimental and theoretical data within 0.5 to 1 kcal/mol. The SAR generally allows conclusive predictions as to the fate of alkoxy radicals; several examples concerning oxy radicals from prominent atmospheric VOC are presented. Specific limitations of the SAR are also discussed. Using the predicted barrier height Eb, the high-pressure rate coefficient for alkoxy decomposition k diss (298 K) can be obtained from k diss (298 K) = L ×1.8 × 1013 exp(–Eb/RT) s–1, with L the reaction path degeneracy.  相似文献   

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

4.
Wind and temperature profiles in the stable boundary layer were analyzed in the context of MoninObukhov similarity. The measurements were made on a 60-m tower in Kansas during October 1999 (CASES-99). Fluxprofile relationships, obtained from these measurements in their integral forms, were established for wind speed and temperature. Use of the integral forms eliminates the uncertainty and accuracy issues resulting from gradient computations. The corresponding stability functions, which were nearly the same for momentum and virtual sensible heat, were found to exhibit different features under weakly stable conditions compared to those under strongly stable conditions. The gradient stability functions were found to be linear, namely m = 1+ 5.8 and h = 1 + 5.4 up to a limit of the MoninObukhov stability parameter = 0.8; this is consistent with earlier findings. However, for stronger stabilities beyond a transition range, both functions were observed gradually to approach a constant, with a value of approximately 7. To link these two distinct regimes, a general but pliable functional form with only two parameters is proposed for the stability functions, covering the entire stability range from neutral to very stable conditions.  相似文献   

5.
A one-dimensional numerical model based on the equations of mean motion and turbulent kinetic energy (TKE), with Delage's (1974) mixing-length parameterization has been used to simulate the mean and turbulent structure of the evolving stably stratified nocturnal boundary layer (NBL). The model also includes a predictive equation for the surface temperature and longwave radiational cooling effects.In the absence of advective and gravity wave effects, it is found that the model-simulated structure, after a few hours of evolution, could be ordered fairly well by a similarity scaling (u *0, *0, L 0, and h) based on surface fluxes and the NBL height. Simple expressions are suggested to describe the normalized profiles of momentum and heat fluxes, TKE, eddy-viscosity and energy dissipation. A good ordering of the same variables is also achieved by a local scaling (u *0, * and L) based on the height-dependent local fluxes. The normalized TKE, eddy viscosity and energy dissipation are unique functions of z/L and approach constant values as z/L , where L is the local Monin-Obukhov length. These constants are close to the values predicted for the surface layer as z/L , thus suggesting that the Monin-Obukhov similarity theory can be extended to the whole NBL, by using the local (height-dependent) scales in place of surface-layer scales. The observed NBL structure has been shown to follow local similarity (Nieuwstadt, 1984).  相似文献   

6.
From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - Ci constants i=u, v, w, (see Equation (5) - kme,kmT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - kmi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - ksb, kc characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - m nondimensional mean shear - H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - D nondimensional divergence of vertical transports of turbulent linetic energy  相似文献   

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

8.
Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T * ). Normal atmospheric monsoon circulations are defined by long-term average winds (u n,v n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u L,v L); 7–20 day filtered winds (u M,v M); and 2–6 day filtered winds (u S,v S). Three-month running mean anomaly kinetic energy (signified asK L , K M , andK S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK L , K M , andK S perturbations propagate slowly eastward with same phase speed (0.3 ms–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K L , K M , K S ) and interannualu (OLR) modes. Namely, (K L , K M , K S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K L , K M , K S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T * 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT * variations on interannual time scales. HighestT * occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT * variations are also controlled by the prevailing monsoonal zonal windsu n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.  相似文献   

9.
The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (Fa) and atmospheric moisture diffusivity (Kq). In a large ensemble of spinup experiments, the values of Fa and Kq are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of Fa and Kq, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.  相似文献   

10.
During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - h , v horizontal and vertical turbulent viscosity - K h ,K v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h s deviation from the average water surface elevation - L *,H * length and depth scale - U *,W * horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u * friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C 1,C 2,C 3 dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T bx ,T by bottom stress components - c bottom drag coefficient  相似文献   

11.
A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO3 (0.01), N2O5 (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.  相似文献   

12.
Parameterizing turbulent diffusion through the joint probability density   总被引：3，自引：3，他引：0  
The convective mass flux parameterization often used in meteorological modeling expresses the vertical flux of a transported scalar as proportional to the product of the difference in mean values of the scalar in updrafts and downdrafts and their characteristic velocity. The proportionality factor is a constant to be specified. We show that this proportionality factor also appears in the relaxed eddy accumulation technique of Businger and Oncley. That associates the surface-layer flux of a scalar with the product of the standard deviation of vertical velocity and the mean concentration difference between updrafts and downdrafts.We show that this constant (b) is determined uniquely by the joint probability density (jpd) of vertical velocity and the scalar. Using large-eddy simulation, we generate this jpd for a conservative scalar diffusing through a convective boundary layer. It has quite different forms in top-down and bottom-up diffusion geometries. The bottom-up jpd is fairly well represented by a jointly Gaussian form and implies b ~ 0.6, in good agreement with the surface-layer value reported by Businger and Oncley. The top-down jpd is strikingly non-Gaussian and gives b ~ 0.47. Updrafts carry the bulk of the scalar flux - 70% in the bottom-up case, 60% in the top-down case.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

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

14.
It is shown that the ratio of standard deviation of lateral velocity to the friction velocity, /u *, and therefore wind direction fluctuations, are sensitive to mesoscale terrain properties. Under neutral conditions, /u * is almost 40% larger in rolling terrain than over a horizontal surface. In the lee of a low mountain, the fluctuations may be 2.5 times as strong as over horizontal terrain. In contrast, vertical velocity fluctuations are little influenced by mesoscale terrain features.Now with Air Weather Service, Offutt AFB, Omaha, Nebraska.  相似文献   

15.
Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - 0 zonal average - vertical mean of the zonal average - 0s zonal average at the surface - 0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m2/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - * saturated humidity mixing ratio - opacity of the atmosphere - 1,2 factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - D subsurface temperature (K) - a specific volume - 0xs ,0ys eastward and northward components of surface frictional stress - * vertical velocity at the top of the boundary layer (N/m2/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P a pressure of the model atmosphere (N/m2) - P s pressure at the surface (N/m2) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m2/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C p specific heat of air at constant pressure (J/Kg/K) - k R/C p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H s H 0s (1) +H 0s (2) +H 0s (3) +H 0s (4) +H 0s (5) (J/m2/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H a H 0a (1) +H 0a (2) +H 0a (3) +H 0a (4) (J/m2/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H 0s (1) ,H 0a (1) heat flux due to short wave radiation - H 0s (2) ,H 0a (2) heat flux due to long wave radiation - H 0s (3) ,H 0a (3) heat flux due to small scale convection - H 0s (4) heat flux due to evaporation - H 0a (4) heat flux due to condensation - H 0s (5) heat flux due to subsurface conduction and convection - e * saturation vapor pressure - R solar constant (W/m2) - r a albedo of the atmosphere - r s albedo of the surface - b 2 empirical constant (J/m2/sec) - c 2 empirical constant (J/m2/sec) - e 2 nondimensional empirical constant - f 2 empirical constant (J/m2/sec) - factor proportional to the conductive capacity of the surface medium - a s constant used in Sellers model - b s positive constant of proportionality used in the Sellers model (kg m2/J/sec2) - K HT coefficient for eddy diffusivity of heat (m2/sec) - K HE exchange coefficient for water vapor (m2/sec) - h depth of the water column (m) - z height (m) - V 0ws meridional component of surface current (m/sec) - n cloud amount - G 0,n long wave radiation form the atmosphere for cloud amount n (W/m2) - B 0 long wave radiation from the surface (W/m2) - S 0,n short wave radiation from the atmosphere for cloud amount n (W/m2) - A n albedo factor for a cloud amount n - R f1 large scale rainfall (mm/day) - R f2 small scale rainfall (mm/day) With 22 Figures  相似文献   

16.
This paper considers the ground area which affects the properties of fluid parcels observed at a given spot in the Planetary Boundary Layer (PBL). We examine two source-area functions; the footprint, giving the source area for a measurement of vertical flux: and the distribution of contact distance, the distance since a particle observed aloft last made contact with the surface. We explain why the distribution of contact distance extends vastly farther upwind than the footprint, and suggest for the extent of the footprint the inequalities: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOabaeqabaGaamyvam% aalaaabaGaamiAaaqaaiabeo8aZnaaBaaaleaacaWGxbaabeaakiaa% cIcacaWGObGaaiykaaaacqGH8aapcaWG4bGaeyipaWJaamyvaKazaa% iadaGabaqaamaaDaaajqwaacqaaiaadIgacaGGVaGabmOEayaacaGa% aiilaiaabccacaGGVbGaaiiDaiaacIgacaGGLbGaaiOCaiaacEhaca% GGPbGaai4CaiaacwgaaeaacaWGubWaaSbaaKazcaiabaGaamitaaqa% baqcKfaGaiaacIcacaWGObGaaiykaiaabYcacaqGGaGaaeiAaiaabc% cacaGGHbGaaiOyaiaac+gacaGG2bGaaiyzaiaabccacaGGZbGaaiyD% aiaackhacaGGMbGaaiyyaiaacogacaGGLbGaeyOeI0IaaiiBaiaacg% gacaGG5bGaaiyzaiaackhaaaaajqgaacGaay5EaaaakeaaaeaacaGG% 8bGaamyEaiaacYhacqGH8aapcqaHdpWCdaWgaaWcbaGaamODaaqaba% GccaGGOaGaamiAaiaacMcadaWcaaqaaiaadIhaaeaacaWGvbaaaaaa% aa!7877!\[\begin{array}{l} U\frac{h}{{\sigma _W (h)}} < x < U\left\{ {_{h/\dot z,{\rm{ }}otherwise}^{T_L (h){\rm{, h }}above{\rm{ }}surface - layer} } \right. \\ \\ |y| < \sigma _v (h)\frac{x}{U} \\ \end{array}\] where U is the mean streamwise (x) velocity, h is the observation height, L is the Lagrangian timescale, v and w are the standard deviations of the cross-stream horizontal (y) and vertical (z) velocity fluctuations, and is the Lagrangian Similarity prediction for the rate of rise of the centre of gravity of a puff released at ground.Simple analytical solutions for the contact-time and the footprint are derived, by treating the PBL as consisting of two sub-layers. The contact-time solutions agree very well with the predictions of a Lagrangian stochastic model, which we adopt in the absence of measurements as our best estimate of reality, but the footprint solution offers no improvement over the above inequality.  相似文献   

17.
Summary A simple parameterization for the estimation of turbulent kinetic energy (TKE) and momentum flux profiles under near-neutral stratification based on sodar measurements of the vertical velocity variance has been tested using data from the LINEX-2000 experiment. Measurements included operation of a phased-array Doppler sodar DSDPA.90 and of a sonic anemometer USA-1 mounted at a meteorological tower at a height of 90m. Good agreement has been found between the TKE and momentum flux values derived from the sonic and sodar data (with correlation coefficients r>0.90 and a slope of the regression lines of about 1.01.1) suggesting the possible use of sodar measurements of w 2 to derive turbulence parameter profiles above the tower range.  相似文献   

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

19.
The present study explores the extent to which the logarithmic region of the adiabatic atmospheric boundary layer can be modeled using a three-dimensional large eddy simulation. A value of the von Kármán constant (LES) is obtained by determining the slope of a logarithmic portion of the velocity profile. Its numerical value is found to be dependent on the value of the Smagorinsky-Model Reynolds number, ReSM: the value of LES increases with ReSM. Results indicate that LES approaches a value of 0.35 as ReSM reaches about 7.75 × 105 for the largest domain. The sensitivity of LES to the profile region over which it is evaluated has been tested. Results show that LES is not sensitive to the depth of this evaluation region when we employ five grids above the sub-grid buffer layer where sub-grid-scale effects dominate. The maximum LES is obtained when the lower boundary of the evaluation region is just above the top of the sub-grid-scale buffer layer. This result is consistent with modelled mean speed and resolved-scale shear stress profiles.  相似文献   

20.
Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.  相似文献   

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

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