Monte Carlo simulation of nitrogen oxides dispersion from a vehicular exhaust plume and its sensitivity studies

The pollutant dispersion behavior from the vehicular exhaust plume has a direct impact on human health, particularly to the drivers, bicyclists, motorcyclists, pedestrians, people working nearby and vehicle passengers. A two-dimensional pollutant dispersion numerical model was developed based on the joint-scalar probability density function (PDF) approach coupled with a k–ε turbulence model to simulate the initial dispersion process of nitrogen oxides, temperature and flow velocity distributions from a vehicular exhaust plume. A Monte Carlo algorithm was used to solve the PDF transport equations in order to obtain the dispersion distribution of nitrogen oxides concentration. The model was then validated by a series of sensitivity experimental studies in order to assess the effects of vehicular exhaust tailpipe velocities, wind speeds and chemistry on the initial dispersion of NO and NO₂ mass concentrations from the vehicular exhaust plume. The results show that the mass concentrations of nitrogen oxides decrease along the centerline of the vehicular exhaust plume in the downstream distance. The dispersion process can be enhanced when the vehicular exhaust tailpipe velocity is much larger than the wind speed. The oxidation reaction of NO plays an important role when the wind speed is large and the vehicular exhaust exit velocity is small, which leads to chemical reduction of NO, and the formation and accumulation of NO₂ in the exhaust plume. It is also found that the effect of vehicular exhaust-induced turbulence in the vicinity of the exhaust tailpipe exit is more dominant than the effect of wind turbulence, while the wind turbulence gradually shows a significant role for the dispersion of nitrogen oxides along with the development of exhaust plume. The range of dispersion of nitrogen oxides in the radial direction is increased along with the development of vehicular exhaust plume.     

Boundary Layer Emission Monitoring

A new methodology is described for determining the atmospheric emission rate of pollutants from large heterogeneous area sources, such as hazardous waste sites. The procedure hinges upon measuring average pollutant concentrations, at three or more different elevations, while traversing the plume downwind of the area source. A helium-filled tethersonde balloon is used to elevate the sampling lines to their appropriate height. During plume traversing the sampling rate is adjusted to be proportional to the sine of the angle between the wind vector and the direction of the traverse path. The average concentrations are corrected for any upwind, background concentration and then used to derive an average vertical concentration profile. This profile Is numerically integrated, with the wind velocity profile, over the pollutant boundary layer to yield the area source emission rate. The methodology was tested on several large industrial effluent lagoons and proved to be easy to use, robust, and precise.     

等效排气筒有效高度计算方法的比选及其应用于大气污染预测的可行性

目前,等效排气筒多用于大气污染物总量控制,其预测精度和范围的不明晰限制了其在污染物运移扩散领域的进一步应用。基于《大气污染物综合排放标准》,依据所预测的范围和浓度精度对8种典型等效计算方法进行了比选,并验证了将等效排气筒用于不同工况下污染物运移扩散预测的可行性。改进的有效高度等效算法 (源强加权算术平均法) 综合考虑了不同高度和源强参数特征,以2个排放同种污染物的相邻排气筒为例,所计算的高斯模式下等效后下风向污染物浓度场总体分布趋势与等效前叠加计算结果一致,且预测精度优于《大气污染物综合排放标准》中提出的均方根平均法和其他等效算法。对不同风速条件下 (1.5~4.5 m·s⁻¹) 等效前后下风向污染物浓度场分布计算比较,发现即使风速改变仍可保证较高的最大落地浓度预测精度 (−6.87%~−2.21%),特别是风速较大时其预测精度更高 (达到−2.21%) 。这验证了该方法的有效性和稳定性。本研究探讨的源强加权算术平均值算法,进一步提升了等效排气筒相关参数计算的合理性,并拓展了其在大气预测评价中的应用。     

The effects of roadside structures on the transport and dispersion of ultrafine particles from highways

Understanding local-scale transport and dispersion of pollutants emitted from traffic sources is important for urban planning and air quality assessments. Predicting pollutant concentration patterns in complex environments depends on accurate representations of local features (e.g., noise barriers, trees, buildings) affecting near-field air flows. This study examined the effects of roadside barriers on the flow patterns and dispersion of pollutants from a high-traffic highway in Raleigh, North Carolina, USA. The effects of the structures were analyzed using the Quick Urban & Industrial Complex (QUIC) model, an empirically based diagnostic tool which simulates fine-scale wind field and dispersion patterns around obstacles. Model simulations were compared with the spatial distributions of ultrafine particles (UFP) from vehicular emissions measured using a passenger van equipped with a Differential Mobility Analyzer/Condensation Particle Counter. The field site allowed for an evaluation of pollutant concentrations in open terrain, with a noise barrier present near the road, and with a noise barrier and vegetation present near the road.Results indicated that air pollutant concentrations near the road were generally higher in open terrain situations with no barriers present; however, concentrations for this case decreased faster with distance than when roadside barriers were present. The presence of a noise barrier and vegetation resulted in the lowest downwind pollutant concentrations, indicating that the plume under this condition was relatively uniform and vertically well-mixed. Comparison of the QUIC model with the mobile UFP measurements indicated that QUIC reasonably represented pollutant transport and dispersion for each of the study configurations.     

A model for the maximum credible hourly impact on any ground receptor from point sources with thermal plume rise

A pollutant dispersion model is developed, allowing fast evaluation of the maximum credible 1-h average concentration on any given ground-level receptor, along with the corresponding critical meteorological conditions (wind speed and stability class) for stacks with buoyant plumes in urban or rural areas. Site-specific meteorological data are not required, as the computed concentrations are maximized against all credible combinations of wind speed, stability class, and mixing height. The analysis is based on the dispersion relations of Pasquill-Gifford and Briggs for rural and urban settings, respectively, the buoyancy induced dispersion correlation of Pasquill, the wind profile exponent values suggested by Irwin, the buoyant plume rise relations of Briggs, as well as the Benkley and Schulman's model for the minimum mixing heights. The model is particularly suited for air pollution management studies, as it allows fast screening of the maximum impact on any selected receptor and evaluation of the ways to have this impact reduced. It is also suited for regulatory purposes, as it can be used to define the minimum stack size requirements for a given source as a function of the exit gas volume and temperature, the pollutant emission rates and their hourly concentration standards, as well as the source location relative to sensitive receptors.     

Traffic atmospheric diffusion model

A digital computer model simulation of traffic flow using car-following theory and constrained random input has been developed to predict the concentration of gaseous pollutants (i.e. NO₂ and HC) in the atmosphere emitted from different types of vehicles on the roads, using a simple Gaussian dispersion point source plume.Measurements were made on different roads of Baghdad city. Atmospheric pollutant concentrations, traffic flow rate, traffic composition and climatic conditions were recorded.The prediction method was then tested against measurements. The predicted levels are in acceptable agreement with the measured values. The model was then used to examine various traffic and pollutant control strategies.     

道路路障对街道峡谷内污染物扩散的影响

采用数值模拟,研究不同风向角α (α=0°、45°、90°)及道路屏障位置(中间单路障和两侧双路障)对街道峡谷内机动车尾气污染物扩散的影响。数值模拟采用标准 k-ε 湍流模型且Sct选择0.3时,计算结果与风洞实验结果较好吻合。结果表明,2种路障布置方式可有效降低人行道内污染物浓度,特别是,当α=45°时,污染物浓度最多可降低46.23%。同时,风向角α对街道峡谷内污染物扩散影响较大。当 α=90°时,空气流通不良使得污染程度最为严重,且污染集中在背风侧近地面。单路障比双路障布置对污染物扩散影响更大,前者使污染物主要集中在街道中心背风侧,其他位置浓度明显降低;双路障时仅在一定范围内改善人行道内空气品质,但对街道整体污染物分布影响不大。     

An analytical diffusion model for long distance transport of air pollutants

A steady-state two-dimensional diffusion model suitable for predicting ambient air pollutant concentrations averaged over a long time period (e.g., month, season or year) and resulting from the transport of pollutants for distances greater than about 100 km from the source is described. Analytical solutions are derived for the primary pollutant emitted from a point source and for secondary pollutant formed from it. Depletion effects, whether due to wet or dry deposition or chemical conversion to another species, are accounted for in these models as first order processes. Thus, solutions for multiple point sources may be superimposed.In this model the time-averaging of the random trajectories of pollutant-contaminated air parcels is represented by horizontal diffusion in a steady, two-dimensional flow field of the time-averaged wind. The resulting concentration isopleths for a point source show significant dispersion both upwind and cross wind of the source with respect to the mean wind field.The analytical theory for the dispersion of a primary pollutant is compared with the numerical predictions of a plume trajectory model for the case of steady emission from a point source. Good overall agreement between the two models is achieved whether or not depletion by wet and dry deposition is included.The theory for the dispersion of a secondary pollutant is compared with measurements of the annual average sulfate concentration in the U.S. Calculations are carried out using SO₂ emissions from electric power plants in the United States as a source inventory. Using optimum values of the dispersion parameters, the correlation coefficient of observed and calculated ambient concentrations is 0.87 for the eastern United States and 0.69 for the western region. The optimum dispersion parameters used are comparable to values quoted in the literature.The horizontal length scale characterizing the sulphate concentration distribution from a single source is about 500 km, being noticeably larger than that characterizing the primary (sulfur dioxide) distribution. Using optimum dispersion parameters, a point source of 33 kg s⁻¹ of sulfur dioxide would give rise to a maximum annual average sulfate concentration of 1 μ m⁻³.A calculation of annual average SO₂ concentrations in the United States is carried out using previously derived optimal values of the parameters from the sulfate calculation. The resulting isopleths are similar to measured values in the eastern U.S.     

Pollutant dispersion in urban street canopies

The effects of building configurations on pollutant dispersion around street canopies were studied numerically. The dispersion of pollutants emitted from ground sources was simulated by continuously discharging large number of particles into the computation domain. The mean wind velocities at each time-step were firstly computed by solving the time-dependent incompressible Navier–Stokes equations, while the fluctuated velocities were determined using a statistical procedure. The trajectories of the discharged particles were obtained from a Lagrangian particle model. Three categories of numerical simulation were conducted to study the effect of different canopy geometries on the pollutant dispersion. The computed wind field data were consistent with the wind field characteristics described in the previous wind tunnel studies. A counter-clockwise vortex was found resulting in high pollutant concentration at the windward side of the downstream building of the street canopy and low pollutant concentration at the leeward side of the upstream building. The increase in height of the urban roughness buildings would facilitate the pollutant dispersion in urban street canopy under certain building configurations. Two or more vortices stacked vertically in a street canopy were found when height of the upstream and downstream buildings of a street canopy was increased, preventing pollutants from escaping out of the canopy.     

Numerical modeling on air quality in an urban environment with changes of the aspect ratio and wind direction

Due to heavy traffic emissions within an urban environment, air quality during the last decade becomes worse year by year and hazard to public health. In the present work, numerical modeling of flow and dispersion of gaseous emissions from vehicle exhaust in a street canyon were investigated under changes of the aspect ratio and wind direction. The three-dimensional flow and dispersion of gaseous pollutants were modeled using a computational fluid dynamics (CFD) model which was numerically solved using Reynolds-averaged Navier–Stokes (RANS) equations. The diffusion flow field in the atmospheric boundary layer within the street canyon was studied for different aspect ratios (W/H?=?1/2, 3/4, and 1) and wind directions (θ?=?90°, 112.5°, 135°, and 157.5°). The numerical models were validated against wind tunnel results to optimize the turbulence model. The numerical results agreed well with the wind tunnel results. The simulation demonstrated that the minimum concentration at the human respiration height within the street canyon was on the windward side for aspect ratios W/H?=?1/2 and 1 and wind directions θ?=?112.5°, 135°, and 157.5°. The pollutant concentration level decreases as the wind direction and aspect ratio increase. The wind velocity and turbulence intensity increase as the aspect ratio and wind direction increase.     

Estimation of atmospheric trace metal emissions in Vilnius City, Lithuania, using vertical concentration gradient and road tunnel measurement data

A new approach for the estimation of trace metal emissions in Vilnius city was implemented, using vertical concentration profiles in the urban boundary layer and road tunnel measurement data. Heavy metal concentrations were examined in fine and coarse particle fractions using a virtual impactor (cut-off size diameter 2.5 μm). Negative vertical concentration gradients were obtained for all metals (Ba, Pb, V, Sb, Zn) and both fractions. It was estimated that the vertical concentration gradient was formed due to emissions from an area of about 12 km². Road tunnel measurements indicated that trace metal concentrations on fine particles were lower than those on coarse particles, which suggested that re-emitted road dust was highly enriched in trace metal due to historic emissions within the tunnel. Emission rates of different pollutants in the road tunnel were calculated using pollutant concentration differences at the tunnel entrance and exit and traffic flow data. Heavy metal emission rates from the area of Vilnius city were estimated using the vertical gradient of heavy metal concentrations and the coefficient of turbulent mixing, as derived from meteorological measurement data. The emission values calculated by the two different methods coincided reasonably well, which indicated that the main source of airborne trace metals in Vilnius city is traffic. The potential of the vertical concentration gradient method for the direct estimation of urban heavy metal emissions was demonstrated.     

On the prediction of concentration variations in a dispersing heavy-duty truck exhaust plume using k–ε turbulent closure

This work presents the computational fluid dynamic modeling of an exhaust plume dispersed from the exhaust pipe of a class-8 tractor truck powered by 330 hp Cummins M11 electronically controlled diesel engine. This effort utilizes an advanced CFD technique to accurately predict the variation of carbon dioxide concentration inside a turbulent plume using a k–ε eddy dissipation model. The simulation includes the “real-world” operation of a truck and its exhaust plume in a NASA, Langley aircraft testing wind tunnel, that had an effective volume of 226, 535 m³ (8,000,000 ft³). The predicted results show an excellent agreement with the experimentally measured values of CO₂ concentrations, dilution ratios, and the temperature variations inside the plume. A specific goal of this effort was to study the effect of recirculation region near the truck walls on dispersion of the plume. For this purpose, growth of the plume from the center of the exhaust pipe is also presented and discussed. This work also shows the benefits of CFD modeling in applications where dispersion correlations are not required a priori, instead the dispersion coefficients are calculated precisely by solving the turbulent kinetic energy and dissipation equations.     

Dispersion of Traffic and Space Heating Emissions In Urban Settings

A long-term dispersion model is presented for traffic and space heating emissions in urban areas, allowing fast assessment of the spatial-averaged and center-maximum pollutant concentrations.

The assumption of study areas with circular shape and normal emissions density profiles is made for the purpose of streamlining model inputs with the inventory data normally available. In addition, the rather typical assumptions of Gaussian dispersion, narrow plume, flat or gently rolling terrain, homogeneous wind field and nonreactive pollutants are made. Values of σz from Briggs correlation are used with an initial value of 30 to account for building effects.

Meterological data inputs are reduced to six parameters, inventory data inputs to two, while computations are simplified to a degree that use of a digital computer is not required.

The model is well suited to yield separate assessments for individual types of sources and control measures, as well as to reveal sensitivities from parameters such as city size, or emission density levels and distribution patterns. Its predictions are virtually identical to those of the CDM-2 UNAPMAR model for study areas with circular shape and normal emissions density profiles, and as results do not appear overly sensitive to shape and distribution patterns, the model is believed to be valid for most urban areas.     

Flow and dispersion in an urban cubical cavity

Flow and dispersion in an urban cubical cavity are numerically investigated using a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–? turbulence closure model. The urban cubical cavity is surrounded by flank walls that are parallel to the streamwise direction, called end-walls, as well as upstream and downstream walls. A primary vortex and secondary vortices including end-wall vortices are formed in the cavity. Because of the end-wall drag effect, the averaged mean-flow kinetic energy in the cavity is smaller than that in an urban street canyon that is open in the along-canyon direction. A trajectory analysis shows that the end-wall vortices cause fluid particles to move in the spanwise direction, indicating that flow in the cavity is essentially three-dimensional. The iso-surfaces of the Okubo–Weiss criterion capture cavity vortices well. The pollutant concentration is high near the bottom of the upstream side in the case of continuous pollutant emission, whereas it is high near the center of the primary vortex in the case of instantaneous pollutant emission. To get some insight into the degree of pollutant escape from the cavity according to various meteorological factors, extensive numerical experiments with different ambient wind speeds and directions, inflow turbulence intensities, and cavity-bottom heating intensities are performed. For each experiment, we calculate the time constant, which is defined as the time taken for the pollutant concentration to decrease to e^?1 of its initial value. The time constant decreases substantially with increasing ambient wind speed, and tends to decrease with increasing inflow turbulence intensity and cavity-bottom heating intensity. The time constant increases as the ambient wind direction becomes oblique. High ambient wind speed is found to be the most crucial factor for ventilating the cavity, thus improving air quality in an urban cubical cavity.     

Size and Concentration Measurements of Particles Produced in Commercial Chromium Plating Processes

Abstract

Optical measurements of particle size and concentration were made at the chromium plating tank and exhaust system at a commercial hexavalent chromium plating facility. Particles were examined at three locations in the exhaust system: 1) directly at the hexavalent chromium plating bath surface, 2) at the exit of a cyclone separator located in the exhaust system approximately three to four meters downstream of the bath, and 3) in the exhaust stack, downstream of the induced draft fan and all abatement devices. Particle diameters at the bath surface ranged from 0.3 to 25 μm. Downstream of the cyclone exit and mesh pad filters, particle top sizes were approximately 5 and 0.7 mm, respectively. On a mass basis, the collection efficiency of all abatement devices was 99.997%. Assuming that droplets in the flow consist primarily of water and chromium, correcting the total particle mass flow against water content gives a chromium emission rate of 64,000 μg/hr, which compares favorably with a value of 77,000 μg/hr measured with EPA methods. This initial agreement, which should be validated through additional measurements over a broad range of flow conditions, raises the possibility of continuous monitoring for chromium metal emissions using particle size/mass as a surrogate.     

Measurement and Three-Dimensional Modeling of Airflow and Pollutant Dispersion in an Undersea Traffic Tunnel

Abstract

Airflow and pollutant dispersion in a cross-harbor traffic tunnel were experimentally and numerically studied. Concentrations of the gaseous pollutants CO, NO_x, and total hydrocarbons (THC) at three axial locations in the tunnel, together with traffic flow rate, traffic speed, and types of vehicle were measured. Three-dimensional (3D) turbulent flow and dispersion of air pollutants in the tunnel were modeled and solved numerically using the finite volume method. Traffic emissions were modeled accordingly as banded line sources along the tunnel floor. The results reveal that cross-sectional concentrations are nonuniformly distributed and that concentrations rise with downstream distance. The piston effect of vehicles alone can provide 9–23% dilution of air pollutants in the tunnel, compounded to a 23–74% dilution effect according to the ventilation condition.     

Turbulent Plume in a Turbulent Cross Flow: Comparison of Wind Tunnel Tests with Field Observations

Plume rise downwind of a large stationary gas turbine was measured in the field and the conditions were then scaled in the laboratory. For the laboratory, the plume exit conditions, wind velocity and temperature profiles, and wind direction were matched. It was found that for high temperature exhaust, the buoyancy is best matched by calculating a dimensionless density difference. With properly calculated buoyancy length scales, the plume trajectories were compared and were found to agree quite well. The probability distributions of the entrainment constant and the average values of the entrapment constant with downwind distance were compared. The field data showed about 15% greater plume rise. The median entrainment constant was about 10% greater for the lab test and the shape of the probability distribution matched very closely.     

Investigation of respirable suspended particulate trend and relevant environmental factors in Hong Kong downtown areas

As the addressing of high demand of good air quality in urban area, a study on air pollutant dispersion and distribution resulting from vehicular exhaust emission is strongly required. In particular, vehicular exhaust emission has become a major air pollution source in metropolitan city like Hong Kong, which is characterized with the heavy, dense traffic flow and the limited land resources. Respirable suspended particulate (RSP) is one of main pollutants resulted from vehicular exhaust emission in urban area. Hence, in this study, we focus on analyzing the variation of RSP levels including diurnal, monthly and annual patterns at selected roadsides in Hong Kong during the period of 1998--2005. Furthermore, the relationships between RSP level and the relevant meteorological factors such as temperature, rainfall and wind conditions in Hong Kong territory have been discussed as well.     

Panel Discussion

A stack design procedure is developed which accounts for the effect of plume interception by downwind buildings, and which provides information on effluent concentrations in a form useful to planning authorities. The information presented in this paper is directed to engineers carrying out stack designs for locations where downwind buildings are of comparable height to the stack. A wind tunnel investigation using tracer gas techniques indicates that, for a plume at building height, downwash on the upwind face of a building causes the high concentrations observed near the roof to be transported to ground level. The effect of a plume on elevated points is determined by the concept of the minimum descent height of the maximum allowable ambient concentration isopleth. This minimum descent height, computed using Gaussian plume dispersion theory, defines a building height below which pollutant concentrations will always lie within safe limits. A case study is presented for the use of the design procedure for a small thermal power plant in an urban area.     

Analysis of Motorcycle Exhaust Regular Testing Data—A Case Study of Taipei City

Abstract

In Taiwan, a continuous increase in the number of motorcycles has made exhaust pollution one of the major emission sources of air pollutants. The regular testing program carried out by the Republic of China Environmental Protection Agency was designed to reduce air pollutant emissions by enhancing maintenance and repair. During the execution period, abundant testing results were accumulated to discuss pollutant emissions from motorcycles. Exhaust testing data of motorcycles in Taipei City from 1996 to 2005 were chosen as the basic data to survey changes in motorcycle exhaust. Effects of motorcycle age and mileage on exhaust pollution were studied. The introduction of advanced emission standards enhances the elimination of high-emitting motorcycles. The testing data indicate that the testing rate rose from approximately 50 to 70% and the failure rate changed from approximately 15 to 10%. The operation cycles of two-stroke motorcycles make them high-emitting vehicles. Concentrations of carbon monoxide and hydrocarbons are higher in two-stroke motorcycle exhaust than that in four-stroke motorcycles. In contrast, the concentration of carbon dioxide produced from complete oxidation processes is lower in exhaust from two-stroke motorcycles. Therefore, failure rates of two-stroke motorcycles are higher than those of four-stroke motorcycles and were also observed to deactivate more easily. On the basis of analytical results of testing data, we found that failure rates show a gradually increasing trend for motorcycles older than 3 yr or used for mileages greater than 10,000 km, and failure rates are highly correlated to the age/mileage of motorcycles. We reason that the accumulation of age or mileage means accumulating usage time of engines and emission control systems. Concentrations of pollutant emissions would increase because of engine wear and emission control system deactivation. After discussing changes of failure rates and pollutant emissions, some suggestions are proposed to improve the testing rate and effectiveness of regular testing.