首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 218 毫秒
1.
Turbulent mixing induces variability in concentration that is important in many applications, such as reactive plumes, risk assessments or odour impact analyses (when the effects can have time scales on the order of a second). In urban canopies, the variability may be modified by the presence of buildings. Our purpose is to study concentration fluctuation variance in built-up areas using an Eulerian approach. We performed numerical simulations with the computational fluid dynamics model Mercure_Saturne, which is a three-dimensional model adapted to atmospheric flow and pollutant dispersion. We use a k − ϵ turbulence closure and predict the concentration variance with a transport equation model. The model performance is evaluated with the near-full scale experiment MUST (Mock Urban Setting Test), a field experiment conducted in Utah’s West Desert Test Center. The modelled root-mean-square of the concentration fluctuations is compared to measurements for 20 of the MUST trials. The model shows good agreement with the measurements, with the fraction of predictions within a factor of two of observations of 60.1%, with better results for horizontal lines of detectors than for the detectors on vertical masts (with fractions of predictions within a factor of two of observations of respectively 66.4% and 52.6%). The influence of different parameters on the fluctuation variance is also studied and we show the importance of taking into account the stability of the stratification when modelling the turbulent kinetic energy.  相似文献   

2.
A study of the neutrally-stratified flow within and over an array of three-dimensional buildings (cubes) was undertaken using simple Reynolds-averaged Navier—Stokes (RANS) flow models. These models consist of a general solution of the ensemble-averaged, steady-state, three-dimensional Navier—Stokes equations, where the k-ε turbulence model (k is turbulence kinetic energy and ε is viscous dissipation rate) has been used to close the system of equations. Two turbulence closure models were tested, namely, the standard and Kato—Launder k-ε models. The latter model is a modified k-ε model designed specifically to overcome the stagnation point anomaly in flows past a bluff body where the standard k-ε model overpredicts the production of turbulence kinetic energy near the stagnation point. Results of a detailed comparison between a wind-tunnel experiment and the RANS flow model predictions are presented. More specifically, vertical profiles of the predicted mean streamwise velocity, mean vertical velocity, and turbulence kinetic energy at a number of streamwise locations that extend from the impingement zone upstream of the array, through the array interior, to the exit region downstream of the array are presented and compared to those measured in the wind-tunnel experiment. Generally, the numerical predictions show good agreement for the mean flow velocities. The turbulence kinetic energy was underestimated by the two different closure models. After validation, the results of the high-resolution RANS flow model predictions were used to diagnose the dispersive stress, within and above the building array. The importance of dispersive stresses, which arise from point-to-point variations in the mean flow field, relative to the spatially-averaged Reynolds stresses are assessed for the building array.  相似文献   

3.
Further to our previous work—simulations of flow and dispersion in an oblique wind over the DAPPLE site (Xie and Castro, Atmos Environ 43:2174–2185, 2009)—large-eddy simulations of flows and dispersion over the same site in a wind perpendicular to Marylebone Road and the windward surfaces of most of the buildings were performed. The DAPPLE site is located at the intersection of Marylebone Road and Gloucester Place in central London. In order to investigate the effects of wind direction on flows and dispersion, the velocity and scalar fields in the perpendicular wind were compared with those in the oblique wind. Furthermore, realistic wind conditions measured on the BT Tower at 190 m above street level were processed and used to drive the numerical simulations of flows and dispersion at the DAPPLE site. This leads to significant predictive improvements of the dispersion compared with field measurements, which provides validation and confidence for coupling mesoscale meteorological models, e.g. the UK Met Office’s Unified Model and the NCAR’s Weather Research & Forecasting Model, with the street-scale large-eddy simulation of urban environments.  相似文献   

4.
Air flow inside an array of cubes is simulated. Cubes (edge length 0.15 m) are arranged in a regular array, separated by 0.15 m in the streamwise and spanwise directions. Numerical simulations are performed based on Reynolds-averaged Navier–Stokes equations (RANS), solved in a computational fluid dynamics model (CFD), with standard k–ε turbulent closure (two prognostic equations are solved for the turbulent kinetic energy k and its dissipation ε, respectively). Simulations are validated against wind-tunnel data using a technique based on hit-rate calculations, and calculated statistical parameters. The results show that the horizontal velocity is very well modelled, and despite some discrepancies, the model that fulfils the hit-rate test criteria gives useful results that are used to investigate three-dimensional (3-D) flow structures. The 3-D analysis of the flow shows interesting patterns: the centre of the canyon vortex is at 3/4 of the canyon height, and stronger downward than upward motions are present within the canyon. Such behaviour is explained by the presence of a compensation flow through the side of the canyon, which enters the canyon from the upper part and exits from the lower part. This complex 3-D structure affects the tracer dispersion, and is responsible for pollutant transport and diffusion.  相似文献   

5.
Simulating turbulent flows in a city of many thousands of buildings using general high-resolution microscopic simulations requires a grid number that is beyond present computer resources. We thus regard a city as porous media and divide the whole hybrid domain into a porous city region and a clear fluid region, which are represented by a macroscopic k–e{\varepsilon} model. Some microscopic information is neglected by the volume-averaging technique in the porous city to reduce the calculation load. A single domain approach is used to account for the interface conditions. We investigated the turbulent airflow through aligned cube arrays (with 7, 14 or 21 rows). The building height H, the street width W, and the building width B are the same (0.15 m), and the fraction of the volume occupied by fluid (i.e. the porosity) is 0.75; the approaching flow is parallel to the main streets. There are both microscopic and macroscopic simulations, with microscopic simulations being well validated by experimental data. We analysed microscopic wind conditions and the ventilation capacity in such cube arrays, and then calculated macroscopic time-averaged properties to provide a comparison for macroscopic simulations. We found that the macroscopic k–e{\varepsilon} turbulence model predicted the macroscopic flow reduction through porous cube clusters relatively well, but under-predicted the macroscopic turbulent kinetic energy (TKE) near the windward edge of the porous region. For a sufficiently long porous cube array, macroscopic flow quantities maintain constant conditions in a fully developed region.  相似文献   

6.
  相似文献   

7.
In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.  相似文献   

8.
To investigate the suitability of computational fluid dynamics (CFD) with regard to windbreak aerodynamics, simulations are performed with a state-of-the-art numerical scheme (Fluent) and compared against experimental data for two- and three-dimensional disturbances, namely the case of a long straight porous shelter fence and the case of a shelter fence erected in a square about an enclosed plot. A thorough sensitivity study quantifies the impact of numerical choices on the simulation (e.g. grid-point density, domain size, turbulence closure), and leads to guidelines that should ensure objective simulation of windbreak flows. On a fine grid Fluent’s “realizable k–ε closure” gives results that are in qualitative accord with the observed mean winds.  相似文献   

9.
The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Re τ  = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0°, 30°, 45°). The results agree well with available experimental data measured in a water channel for a flow angle of 0°. Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays.  相似文献   

10.
The turbulence closure in atmospheric boundary-layer modelling utilizing Reynolds Averaged Navier–Stokes (RANS) equations at mesoscale as well as at local scale is lacking today a common approach. The standard kɛ model, although it has been successful for local scale problems especially in neutral conditions, is deficient for mesoscale flows without modifications. The k–ɛ model is re-examined and a new general approach in developing two-equation turbulence models is proposed with the aim of improving their reliability and consequently their range of applicability. This exercise has led to the replacement of the ɛ-transport equation by the transport equation for the turbulence inverse length scale (wavenumber). The present version of the model is restricted to neutrally stratified flows but applicable to both local scale and mesoscale flows. The model capabilities are demonstrated by application to a series of one-dimensional planetary boundary-layer problems and a two-dimensional flow over a square obstacle. For those applications, the present model gave considerably better results than the standard k–ɛ model.  相似文献   

11.
A study of the oxidation mechanism of N-methyl pyrrolidinone (C5H9NO, NMP) initiated by hydroxyl radicals was made at EUPHORE at atmospheric pressure (1000 ± 10) mbar of air and ambient temperature (T = 300 ± 5 K). The main products were N-methyl succinimide (NMS) (52 ± 4)% and N-formyl pyrrolidinone (FP) (23 ± 9)%. The relative rate technique was used to determine the rate constants of OH with NMP, NMS and FP, the measured values were (in units of cm3 molecule − 1 s− 1): kNMP = (2.2 ± 0.4) × 10− 11, kNMS = (1.4 ± 0.3) × 10− 12 and kFP = (6 ± 1) × 10− 12. The results are presented and discussed in terms of the atmospheric impact.  相似文献   

12.
The SF6 gas tracer observations for puffs released near the ground during the Joint Urban 2003 (JU2003) urban dispersion experiment in Oklahoma City have been analysed. The JU2003 observations, at distances of about 100–1,100 m from the source, show that, at small times, when the puff is still within the built-up downtown domain, the standard deviation of the concentration time series, σt, is influenced by the initial puff spread due to buildings near the source and by hold-up in the wakes of large buildings at the sampler locations. This effect is parameterised by assuming an initial σto of about 42 s, leading to a comprehensive similarity formula: σt = 42 + 0.1t. The second term, 0.1t, is consistent with an earlier similarity relation, σt = 0.1t, derived from puff observations in many experiments over rural terrain. The along-wind dispersion coefficient, σx, is assumed to equal σt u, in which u is the puff speed calculated as the distance from the source to the sampler, x, divided by the time after the release that the maximum concentration is observed at the sampler. σx can be expressed as σx = σxo + 0.14x, with the initial σxo of 45 m. This initial σxo agrees with the suggestion of an initial plume spread of about 40 m, made by McElroy and Pooler from analysis of the 1960s’ St. Louis urban dispersion experiment. The puff speeds, u, are initially only about 20% of the observed wind speed, averaged over about 80 street-level and rooftop anemometers in the city, but approach the mean observed wind speed as the puffs grow vertically. The scatter in the σt data is about ± a factor of two or three at any given travel time. The maximum σt is about 250 s, and the maximum duration of the puff over the sampler, Dt, sometimes called the retention time, is about 1,100 s or 18 min for these puffs and distances.  相似文献   

13.
应用非地转湿Q*,对2001年5月31日至6月2日发生在云南的持续性暴雨天气过程进行诊断分析。结果表明,非地转湿Q矢量辐合区是暴雨发生的有利区域;非地转湿Q矢量的垂直分布反映了次级环流的方向和强弱,云南的暴雨落区位于次级环流的上升支附近。  相似文献   

14.
The reaction kinetics of S(IV) autoxidation catalyzed by Mn(II) in the pH range 3–5 typical for atmospheric liquid water, was investigated. For reactions with pH maintained constant during the reaction course, the predictions obtained by a simple integral approach cover kinetic results only for concentrations of HSO 3 − up to 0.2 mM at pH 4.5. Thus, a generalized simple kinetic model, which can be used for predicting the reaction kinetics in wider concentration, pH and temperature ranges, was derived. This model is based on the assumption that the reaction rate is proportional to the concentration of a transient manganese-sulfito complex formed in the initial step of a radical chain mechanism. In the proposed power law rate equation
the concentration of complex is calculated from the stability constant K and concentrations of reactants at a specific reaction time. This rate equation adequately predicts the reaction kinetics in the pH range 3–5, in the concentration ranges 0.1 ≤ [HSO 3 − ] ≤ 0.4 mM and 2 ≤ [Mn(II)] ≤ 14.6 μM. For the temperature range 15–35 °C, the estimated value for activation energy is 92.0 ± 0.1 kJ mol−1 and the Gibbs free energy of formation of the manganese-sulfito complex is −20.4 ± 0.3 kJ mol−1. Furthermore, the kinetics for catalytic reactions with pH maintained constant during the reaction course as well as with initial pH adjusted only at the start of the reaction, is described satisfactorily by the present model.  相似文献   

15.
应用含湍流频散效应的近地面层的运动方程求解了不同层结下的风速廓线，着重分析了湍流频散效应对近地面层平稳运动的影响。分析指出：湍流的频散效应对经典的幂律廓线一对数修正，该修正在不稳定层结时比稳定层时明显；利用相似理论也得到了该常数。  相似文献   

16.
The atmospheric boundary layer (ABL) model of Weng and Taylor with E−ℓ turbulence closure is applied to simulate the one-dimensional stably stratified ABL. The model has been run for nine hours from specified initial wind, potential temperature and turbulent kinetic energy profiles, and with a specified cooling rate applied at the surface. Different runs are conducted for different cooling rates, geostrophic winds and surface roughnesses. The results are discussed and compared with other models, large-eddy simulations and published field data.  相似文献   

17.
Summary k-day extreme precipitation depths (k=1,2,3, … 30) for the climatological network of Belgium (165 stations) are analysed to detect a possible evolution in the occurrence of extreme rainfall events during the 1951–1995 reference period. The calendar year and the hydrological summer and winter are considered separately. Spearman’s rank correlation coefficient shows a strong spatial correlation between extreme k-day precipitation events, depending on the time of the year (lower during summer than during winter) and increasing with k. In some cases the distances of de-correlation exceed 200 km which is comparable to the size of the country. Due to this correlation, tests for trends have been carried out on the leading principal components (PC) derived from the covariance matrix. Various PC selection rules have been applied to identify the number of components to analyse. The number of components needed to reproduce a given proportion of the total variance varies, with larger values for summer than for winter and a decrease with growing k. The Fisher test is used as a global test. It combines the individual Mann-Kendall trend tests carried out on the selected PC scores. Significant trends have been found in extreme winter k-day precipitation for all the values of k and none in extreme summer precipitation. The results for the annual k-day precipitation depths are between those for the two seasons: no trend for small k because summer events dominate and a significant trend for k larger than 7 due to the winter events. Analysis of a few stations with long-term series shows no significant trend for the period 1910–1995, these series also reproduce almost the same trends as those found for the shorter 1951–1995 period. Received April 23, 1999 Revised December 6, 1999  相似文献   

18.
A comprehensive model for the prediction of concentration fluctuations in plumes dispersing in the complex and highly disturbed wind flows in an urban environment is formulated. The mean flow and turbulence fields in the urban area are obtained using a Reynolds-averaged Navier-Stokes (RANS) flow model, while the standard k-ϵ turbulence model (k is the turbulence kinetic energy and ϵ is the viscous dissipation rate) is used to close the model. The RANS model provides a specification of the velocity statistics of the highly disturbed wind flow in the urban area, required for the solution of the transport equations for the mean concentration and concentration variance (both of which are formulated in the Eulerian framework). A physically-based formulation for the scalar dissipation time scale t d , required for the closure of the transport equation for , is presented. This formulation relates t d to an inner time scale corresponding to “internal” concentration fluctuation associated with relative dispersion, rather than an outer time scale associated with the entire portion of the fluctuation spectrum. The two lowest-order moments of concentration ( and ) are used to determine the parameters of a pre-chosen functional form for the concentration probability density function (clipped-gamma distribution). Results of detailed comparisons between a water-channel experiment of flow and dispersion in an idealized obstacle array and the model predictions for mean flow, turbulence kinetic energy, mean concentration, concentration variance, and concentration probability density function are presented.  相似文献   

19.
An analysis was performed of the turbulent data obtained from Yucheng experimental station in the Shandong Province in 1984. [t is shown that at variant wind speed, the spectra of streamwise velocity remain similar and the intensity of wind fluctuations is proportional to wind speed in the downwind area of shelter belt. Therefore, we may decide the similarity of wind fluctuations by a speed scale and a length scale which is not correlated with stability, σu /V0 = F(X / H). The -5/3 power range of temperature spectra extends to lower frequency. The variation of ratio σ0 /T. with stability becomes σ0 / T . = C(X / H)( - Z / L)-1/3 . There is not such an extension of -5 / 3 power range in the humidity spectra.  相似文献   

20.
The note presents a rational approach to modelling the source/sink due to vegetation or buoyancy effects that appear in the turbulent kinetic energy, E, equation and a supplementary equation for a length-scale determining variable, φ, when two-equation closure is applied to canopy and atmospheric boundary-layer flows. The approach implements only standard model coefficients C φ1 and C φ2 in the production and destruction terms of the φ equation, respectively. Numerical tests illustrate the practical applicability of the method, where, for example, simulations with the E–ω model (where is the specific dissipation and is the dissipation rate of E) properly reproduce both the surface-layer wind profile estimated from the Monin-Obukhov similarity theory and the mixing-height evolution observed above forested terrain in Southern Finland.  相似文献   

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

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