On the correlation of air and pollutant exchange for street canyons in combined wind-buoyancy-driven flow

The ventilation and pollutant transport in a two-dimensional (2D) street canyon of building-height-to-street-width (aspect) ratio h/b = 1 under different unstable stratifications were examined. To characterize the combined wind-buoyancy-driven flow and pollutant transport at different Richardson number Ri, a computational fluid dynamics (CFD) model based on the Reynolds-averaged Navier–Stokes (RANS) equations with the Renormalization Group (RNG) k ? ε turbulence model was adopted. Unlike the isothermal condition, a secondary recirculation is initiated at the ground-level windward corner of the street canyon once the unstable stratification is switched on (Ri < 0). It traps the ground-level pollutant leading to elevated pollutant concentration there. As Ri further decreases, the enlarging secondary recirculation enables direct pollutant removal from its core to the shear layer that offsets the ground-level pollutant accumulation. The ventilation and pollutant removal performance under different unstable stratifications are compared by the air (ACH) and pollutant (PCH) exchange rates, and pollutant retention time (τ). Both the mean and turbulent components of ACH are found to increase with decreasing Ri, suggesting that unstable stratification promotes ventilation in street canyons. Moreover, the CFD results agree well with our theoretical model that ACH² varies linearly with Ri. Turbulent transport originally dominates the pollutant removal under isothermal condition. However, progressive domination of pollutant removal by mean wind can be observed with decreasing stability (decreasing Ri from 0 to ?10.6). The critical value is estimated to be Ri = ?8, below which mean wind is the major pollutant removal carrier. Reduction in τ is also observed with decreasing Ri. Hence, in unstable stratification, pollutant resides shorter time in the street canyon compared with its isothermal counterpart, and the ventilation and pollutant removal are more favorable.     

Three-dimensional modeling of air flow and pollutant dispersion in an urban street canyon with thermal effects

Effects of excess ground and building temperatures on airflow and dispersion of pollutants in an urban street canyon with an aspect ratio of 0.8 and a length-to-width ratio of 3 were investigated numerically. Three-dimensional governing equations of mass, momentum, energy, and species were modeled using the RNG k-epsilon turbulence model and Boussinesq approximation, which were solved using the finite volume method. Vehicle emissions were estimated from the measured traffic flow rates and modeled as banded line sources, with a street length and bandwidths equal to typical vehicle widths. Both measurements and simulations reveal that pollutant concentrations typically follow the traffic flow rate; they decline as the height increases and are higher on the leeward side than on the windward side. Three-dimensional simulations reveal that the vortex line, joining the centers of cross-sectional vortexes of the street canyon, meanders between street buildings and shifts toward the windward side when heating strength is increased. Thermal boundary layers are very thin. Entrainment of outside air increases, and pollutant concentration decreases with increasing heating condition. Also, traffic-produced turbulence enhances the turbulent kinetic energy and the mixing of temperature and admixtures in the canyon. Factors affecting the inaccuracy of the simulations are addressed.     

A methodology to urban air quality assessment during large time periods of winter using computational fluid dynamic models

The representativeness of point measurements in urban areas is limited due to the strong heterogeneity of the atmospheric flows in cities. To get information on air quality in the gaps between measurement points, and have a 3D field of pollutant concentration, Computational Fluid Dynamic (CFD) models can be used. However, unsteady simulations during time periods of the order of months, often required for regulatory purposes, are not possible for computational reasons. The main objective of this study is to develop a methodology to evaluate the air quality in a real urban area during large time periods by means of steady CFD simulations. One steady simulation for each inlet wind direction was performed and factors like the number of cars inside each street, the length of streets and the wind speed and direction were taken into account to compute the pollutant concentration. This approach is only valid in winter time when the pollutant concentrations are less affected by atmospheric chemistry. A model based on the steady-state Reynolds-Averaged Navier–Stokes equations (RANS) and standard k-? turbulence model was used to simulate a set of 16 different inlet wind directions over a real urban area (downtown Pamplona, Spain). The temporal series of NO_x and PM₁₀ and the spatial differences in pollutant concentration of NO₂ and BTEX obtained were in agreement with experimental data. Inside urban canopy, an important influence of urban boundary layer dynamics on the pollutant concentration patterns was observed. Large concentration differences between different zones of the same square were found. This showed that concentration levels measured by an automatic monitoring station depend on its location in the street or square, and a modelling methodology like this is useful to complement the experimental information. On the other hand, this methodology can also be applied to evaluate abatement strategies by redistributing traffic emissions.     

CFD simulation of the atmospheric boundary layer: wall function problems

Accurate Computational Fluid Dynamics (CFD) simulations of atmospheric boundary layer (ABL) flow are essential for a wide variety of atmospheric studies including pollutant dispersion and deposition. The accuracy of such simulations can be seriously compromised when wall-function roughness modifications based on experimental data for sand-grain roughened pipes and channels are applied at the bottom of the computational domain. This type of roughness modification is currently present in many CFD codes including Fluent 6.2 and Ansys CFX 10.0, previously called CFX-5. The problems typically manifest themselves as unintended streamwise gradients in the vertical mean wind speed and turbulence profiles as they travel through the computational domain. These gradients can be held responsible—at least partly—for the discrepancies that are sometimes found between seemingly identical CFD simulations performed with different CFD codes and between CFD simulations and measurements. This paper discusses the problem by focusing on the simulation of a neutrally stratified, fully developed, horizontally homogeneous ABL over uniformly rough, flat terrain. The problem and its negative consequences are discussed and suggestions to improve the CFD simulations are made.     

Evaluation and investigation of the effects of ventilation layout,rate, and room temperature on pollution dispersion across a laboratory indoor environment

The presence of chemicals in laboratories and research centers exposes the staff working at such indoor environment to health risks. In this piece of research, a study was performed on the indoor environment of the Center for Environmental Engineering Research at Sahand University of Technology (Tabriz, Iran). For this purpose, the parameters affecting the dispersion of volatile organic compounds (VOCs), including ventilation rate, room temperature, pollution emission time, venting location, air flow regime within the indoor environment, and the number of vents, were simulated via CFD modeling. The CFD modeling was performed three-dimensionally in unsteady state. In case of turbulent flow within the indoor environment, k–ε turbulence model was used to obtain air velocity profile. Experimental data was used to validate the model. Results of the present research showed that when the venting location is on the ceiling, pollution concentration of 25 ppm can be achieved at some low temperature under a particular set of conditions. However, when the venting location was on the walls close to the pollution source, concentrations as low as 5 ppm and lower were observed within the laboratory indoor environment.

     

Application of computational fluid dynamics to closed-loop bioreactors: II. Simulation of biological phosphorus removal using computational fluid dynamics.

Based on the International Water Association's (London) Activated Sludge Model No. 2 (ASM2), biochemistry rate expressions for general heterotrophs and phosphorus-accumulating organisms (PAOs) were introduced to a previously developed, three-dimensional computational fluid dynamics (CFD) activated sludge model that characterized the mixing pattern within the outer channel of a full-scale, closed-loop bioreactor. Using acetate as the sole carbon and energy source, CFD simulations for general heterotrophs or PAOs individually agreed well with those of ASM2 for a chemostat with the same operating conditions. Competition between and selection of heterotrophs and PAOs was verified using conventional completely mixed and tanks-in-series models. Then, competition was studied in the CFD model. These results demonstrated that PAOs and heterotrophs can theoretically coexist in a single bioreactor when the oxygen input is appropriate to allow sufficient low-dissolved-oxygen zones to develop.     

Development of artificial neural network based metamodels for inactivation of anthrax spores in ventilated spaces using computational fluid dynamics

Linear, quadratic, and artificial neural network (ANN)-based metamodels were developed for predicting the extent of anthrax spore inactivation by chlorine dioxide in a ventilated three-dimensional space over time from computational fluid dynamics model (CFD) simulation data. Dimensionless groups were developed to define the design space of the problem scenario. The Hammersley sequence sampling (HSS) method was used to determine the sampling points for the numerical experiments within the design space. A CFD model, comprised of multiple submodels, was applied to conduct the numerical experiments. Large eddy simulation (LES) with the Smagorinsky subgridscale model was applied to compute the airflow. Anthrax spores were modeled as a dispersed solid phase using the Lagrangian treatment. The disinfectant transport was calculated by solving a mass transport equation. Kinetic decay constants were included for spontaneous decay of the disinfectant and for the reaction of the disinfectant with the surfaces of the three-dimensional space. To enhance the mixing of the disinfectant with the room air, a momentum source was included in the simulation. An inactivation rate equation accounted for the reaction between the spores and the disinfectant. The ANN-based metamodels were most successful in predicting the number of viable bioaerosols remaining in an arbitrary enclosed space. Sensitivity analysis showed that the mass fraction of the disinfectant, inactivation rate constant, and contact time had the most influence on the inactivation of the spores.     

Computational fluid dynamics modeling of laboratory flames and an industrial flare

A computational fluid dynamics (CFD) methodology for simulating the combustion process has been validated with experimental results. Three different types of experimental setups were used to validate the CFD model. These setups include an industrial-scale flare setups and two lab-scale flames. The CFD study also involved three different fuels: C₃H₆/CH₄/Air/N₂, C₂H₄/O₂/Ar, and CH₄/Air. In the first setup, flare efficiency data from the Texas Commission on Environmental Quality (TCEQ) 2010 field tests were used to validate the CFD model. In the second setup, a McKenna burner with flat flames was simulated. Temperature and mass fractions of important species were compared with the experimental data. Finally, results of an experimental study done at Sandia National Laboratories to generate a lifted jet flame were used for the purpose of validation. The reduced 50 species mechanism, LU 1.1, the realizable k-? turbulence model, and the EDC turbulence–chemistry interaction model were used for this work. Flare efficiency, axial profiles of temperature, and mass fractions of various intermediate species obtained in the simulation were compared with experimental data and a good agreement between the profiles was clearly observed. In particular, the simulation match with the TCEQ 2010 flare tests has been significantly improved (within 5% of the data) compared to the results reported by Singh et al. in 2012. Validation of the speciated flat flame data supports the view that flares can be a primary source of formaldehyde emission.

ImplicationsValidated computational fluid dynamics (CFD) models can be a useful tool to predict destruction and removal efficiency (DRE) and combustion efficiency (CE) under steam/air assist conditions in the face of many other flare operating variables such as fuel composition, exit jet velocity, and crosswind. Augmented with rigorous combustion chemistry, CFD is also a powerful tool to predict flare emissions such as formaldehyde. In fact, this study implicates flares emissions as a primary source of formaldehyde emissions. The rigorous CFD simulations, together with available controlled flare test data, can be fitted into simple response surface models for quick engineering use.     

City breathability and its link to pollutant concentration distribution within urban-like geometries

This paper is devoted to the study of pollutant concentration distribution within urban-like geometries. By applying efficiency concepts originally developed for indoor environments, the term ventilation is used as a measure of city “breathability”. It can be applied to analyse pollutant removal within a city in operational contexts. This implies the evaluation of the bulk flow balance over the city and of the mean age of air. The influence of building packing density on flow and pollutant removal is, therefore, evaluated using those quantities. Idealized cities of regular cubical buildings were created with packing density ranging from 6.25% to 69% to represent configurations from urban sprawl to compact cities. The relative simplicity of these arrangements allowed us to apply the Computational Fluid Dynamics (CFD) flow and dispersion simulations using the standard k–? turbulence model. Results show that city breathability within the urban canopy layer is strongly dependent from the building packing density. At the lower packing densities, the city responds to the wind as an agglomeration of obstacles, at larger densities (from about 44%) the city itself responds as a single obstacle. With the exception of the lowest packing density, airflow enters the array through lateral sides and leaves throughout the street top and flow out downstream. The air entering through lateral sides increases with increasing packing density.At the street top of the windward side of compact building configurations, a large upward flow is observed. This vertical transport reduces over short distance to turn into a downward flow further downstream of the building array. These findings suggest a practical way of identifying city breathability. Even though the application of these results to real scenarios require further analyses the paper illustrates a practical framework to be adopted in the assessment of the optimum neighbourhood building layout to minimize pollution levels.     

Impact of building facades and ground heating on wind flow and pollutant transport in street canyons

This paper investigates the impacts of building facades and ground heating on the wind flow and pollutant transport in street canyons using the computational fluid dynamic (CFD) technique. Street canyons of H/W (H representing the building height and W the street width) varied from 0.1 to 2, which covered the basic flow regimes of skimming flow (H/W=1 or 2), wake interference flow (H/W=0.5), and isolated roughness flow (H/W=0.1), were examined in a series of sensitivity tests. Heating that occurred on different surfaces, including ground surface and building façades, posed considerable effects on the street canyon wind flow and pollutant transport compared with those under isothermal conditions. The CFD results showed that the mechanically induced wind flow and pollutant transport were complicated by the buoyancy under temperature stratification. Individual street canyons of different H/W and surface-heating scenarios exhibited their unique wind flow structure and pollutant transport behaviors. Two counter-rotating vortices were calculated in the street canyons of H/W=1, in which the zone of higher pollutant concentration under isothermal conditions was switched from the leeward side to the windward side. In the street canyon of H/W=2, the recirculating wind pattern was perturbed by surface heating that led to the development of either one primary vortex or three closely coupled vortices. Because of the complicated wind structure, the zones of higher pollutant concentration located either on the leeward or windward ground level were subjected to the surface-heating scenarios. Only two vortices were developed inside the street canyon of H/W=0.5. The large primary vortex, centered inside the street canyon, extended above the roof level of the street canyon. Meanwhile, a small secondary vortex was found at the ground-level windward corner whose size results as a function of surface-heating configurations. Finally, in the street canyon of H/W=0.1, an isolated clockwise-rotating vortex was developed beside the leeward building while the wind in the windward side blew in the prevailing wind direction. As a result, air pollutant emitted at the street centerline was unlikely to be carried into the leeward vortex. Instead, it was dispersed rapidly on the windward side before being removed from the street canyon.     

Comparative study of measured and modelled number concentrations of nanoparticles in an urban street canyon

This study presents a comparison between measured and modelled particle number concentrations (PNCs) in the 10–300 nm size range at different heights in a canyon. The PNCs were modelled using a simple modelling approach (modified Box model, including vertical variation), an Operational Street Pollution Model (OSPM) and Computational Fluid Dynamics (CFD) code FLUENT. All models disregarded any particle dynamics. CFD simulations have been carried out in a simplified geometry of the selected street canyon. Four different sizes of emission sources have been used in the CFD simulations to assess the effect of source size on mean PNC distributions in the street canyon. The measured PNCs were between a factor of two and three of those from the three models, suggesting that if the model inputs are chosen carefully, even a simplified approach can predict the PNCs as well as more complex models. CFD simulations showed that selection of the source size was critical to determine PNC distributions. A source size scaling the vehicle dimensions was found to better represent the measured PNC profiles in the lowest part of the canyon. The OSPM and Box model produced similar shapes of PNC profile across the entire height of the canyon, showing a well-mixed region up to first ≈2 m and then decreasing PNCs with increased height. The CFD profiles do correctly reproduce the increase from road level to a height of ≈2 m; however, they do not predict the measured PNC decrease higher in the canyon. The PNC differences were largest between idealised (CFD and Box) and operational (OSPM) models at upper sampling heights; these were attributed to weaker exchange of air between street and roof-above in the upper part of the canyon in the CFD calculations. Possible reasons for these discrepancies are given.     

Development of Artificial Neural Network Based Metamodels for Inactivation of Anthrax Spores in Ventilated Spaces Using Computational Fluid Dynamics

ABSTRACT

Linear, quadratic, and artificial neural network (ANN)-based metamodels were developed for predicting the extent of anthrax spore inactivation by chlorine dioxide in a ventilated three-dimensional space over time from computational fluid dynamics model (CFD) simulation data. Dimensionless groups were developed to define the design space of the problem scenario. The Hammersley sequence sampling (HSS) method was used to determine the sampling points for the numerical experiments within the design space. A CFD model, comprised of multiple submodels, was applied to conduct the numerical experiments. Large eddy simulation (LES) with the Smagorinsky subgrid-scale model was applied to compute the airflow. Anthrax spores were modeled as a dispersed solid phase using the Lagrangian treatment. The disinfectant transport was calculated by solving a mass transport equation. Kinetic decay constants were included for spontaneous decay of the disinfectant and for the reaction of the disinfectant with the surfaces of the three-dimensional space. To enhance the mixing of the disinfectant with the room air, a momentum source was included in the simulation. An inactivation rate equation accounted for the reaction between the spores and the disinfectant. The ANN-based metamodels were most successful in predicting the number of viable bioaerosols remaining in an arbitrary enclosed space. Sensitivity analysis showed that the mass fraction of the disinfectant, inactivation rate constant, and contact time had the most influence on the inactivation of the spores.

IMPLICATIONS This investigation presents a framework for the development of user-friendly models; metamodels for the prediction of the number of viable spores remaining in an indoor room during disinfection from accurate but time-consuming CFD studies. During any decontamination event, to know when to stop pumping in the disinfectant and to know what level of log reduction of the spores have been achieved before even starting decontamination would provide valuable guidance. The neural network based metamodels can be applied to obtain quick and relatively accurate answers. This would be necessary when immediate information is required during emergencies.     

Modeling reactive pollutant dispersion in an urban street canyon

Reactive pollutant dispersion in an urban street canyon with a street aspect ratio of one is numerically investigated using a computational fluid dynamics (CFD) model. The CFD model developed is a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–ε turbulence model and includes transport equations for NO, NO₂, and O₃ with simple photochemistry. An area emission source of NO and NO₂ is considered in the presence of background O₃ and street bottom heating (ΔT=5 °C) with an ambient wind perpendicular to the along-canyon direction. A primary vortex is formed in the street canyon and the line connecting the centers of cross-sectional vortices meanders over time and in the canyon space. The cross-canyon-averaged temperature and reactive pollutant concentrations oscillate with a period of about 15 min. The averaged temperature is found to be in phase with NO and NO₂ concentrations but out of phase with O₃ concentration. The photostationary state defect is small in the street canyon except for near the roof level and the upper downwind region of the canyon and its local minimum is observed near the center of the primary vortex. The budget analysis of NO (NO₂) concentration shows that the magnitude of the advection or turbulent diffusion term is much larger (larger) than that of the chemical reaction term and that the advection term is largely balanced by the turbulent diffusion term. On the other hand, the budget analysis of O₃ concentration shows that the magnitude of the chemical reaction term is comparable to that of the advection or turbulent diffusion term. The inhomogeneous temperature distribution itself affects O₃ concentration to some extent due to the temperature-dependent photolysis rate and reaction rate constant.     

Effects of Strength and Position of Pollutant Source on Pollutant Dispersion Within an Urban Street Canyon

A two-dimensional numerical model for simulating airflow and pollutant dispersion inside an urban street canyon was first developed using the FLUENT code, and then it was validated against a wind tunnel experiment. Then the effects of strength and position of pollutant sources on pollutant dispersion within an urban street canyon were investigated numerically. The numerical results showed that the dimensionless pollutant concentrations within the urban street canyon were independent from the source strength. The results also revealed that the pollutant distributions inside the urban street canyon with a two-lane road were influenced significantly by the positions of the two sources: 1) the closer the two sources were to the street center of the canyon, the lower the pollutant concentrations on the leeward wall and at the human respiration level in the leeward footpath became; 2) the pollutant concentrations on the windward wall and at the human respiration level in the windward footpath were not sensitive to the locations of the two sources as long as the source on the windward lane was situated outside the small recirculation zone at the bottom corner of the canyon windward wall; 3) the pollutant concentrations on the lower parts of the windward and leeward walls as well as in the two footpaths increased greatly when the two sources were moved from outside into the small recirculation zones.     

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.     

Analytical solutions to compartmental indoor air quality models with application to environmental tobacco smoke concentrations measured in a house

This paper derives the analytical solutions to multi-compartment indoor air quality models for predicting indoor air pollutant concentrations in the home and evaluates the solutions using experimental measurements in the rooms of a single-story residence. The model uses Laplace transform methods to solve the mass balance equations for two interconnected compartments, obtaining analytical solutions that can be applied without a computer. Environmental tobacco smoke (ETS) sources such as the cigarette typically emit pollutants for relatively short times (7-11 min) and are represented mathematically by a "rectangular" source emission time function, or approximated by a short-duration source called an "impulse" time function. Other time-varying indoor sources also can be represented by Laplace transforms. The two-compartment model is more complicated than the single-compartment model and has more parameters, including the cigarette or combustion source emission rate as a function of time, room volumes, compartmental air change rates, and interzonal air flow factors expressed as dimensionless ratios. This paper provides analytical solutions for the impulse, step (Heaviside), and rectangular source emission time functions. It evaluates the indoor model in an unoccupied two-bedroom home using cigars and cigarettes as sources with continuous measurements of carbon monoxide (CO), respirable suspended particles (RSP), and particulate polycyclic aromatic hydrocarbons (PPAH). Fine particle mass concentrations (RSP or PM3.5) are measured using real-time monitors. In our experiments, simultaneous measurements of concentrations at three heights in a bedroom confirm an important assumption of the model-spatial uniformity of mixing. The parameter values of the two-compartment model were obtained using a "grid search" optimization method, and the predicted solutions agreed well with the measured concentration time series in the rooms of the home. The door and window positions in each room had considerable effect on the pollutant concentrations observed in the home. Because of the small volumes and low air change rates of most homes, indoor pollutant concentrations from smoking activity in a home can be very high and can persist at measurable levels indoors for many hours.     

Critical review of the building downwash algorithms in AERMOD

The only documentation on the building downwash algorithm in AERMOD (American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model), referred to as PRIME (Plume Rise Model Enhancements), is found in the 2000 A&WMA journal article by Schulman, Strimaitis and Scire. Recent field and wind tunnel studies have shown that AERMOD can overpredict concentrations by factors of 2 to 8 for certain building configurations. While a wind tunnel equivalent building dimension study (EBD) can be conducted to approximately correct the overprediction bias, past field and wind tunnel studies indicate that there are notable flaws in the PRIME building downwash theory. A detailed review of the theory supported by CFD (Computational Fluid Dynamics) and wind tunnel simulations of flow over simple rectangular buildings revealed the following serious theoretical flaws: enhanced turbulence in the building wake starting at the wrong longitudinal location; constant enhanced turbulence extending up to the wake height; constant initial enhanced turbulence in the building wake (does not vary with roughness or stability); discontinuities in the streamline calculations; and no method to account for streamlined or porous structures.

Implications: This paper documents theoretical and other problems in PRIME along with CFD simulations and wind tunnel observations that support these findings. Although AERMOD/PRIME may provide accurate and unbiased estimates (within a factor of 2) for some building configurations, a major review and update is needed so that accurate estimates can be obtained for other building configurations where significant overpredictions or underpredictions are common due to downwash effects. This will ensure that regulatory evaluations subject to dispersion modeling requirements can be based on an accurate model. Thus, it is imperative that the downwash theory in PRIME is corrected to improve model performance and ensure that the model better represents reality.     

Turbulence effects on chemical reactions in smog chamber flows

Turbulent mixing effects on the reaction rate of a non-premixed flow are presented for a moderately slow second-order irreversible chemical reaction. The turbulent mixing process leads to inhomogeneities in the concentration of the reactants. Chemical reactions are normally highly non-linear and large errors can result from using average concentrations in the computation of mean reaction rates. A brief review of the literature on this problem and its application areas is made with particular emphasis placed on near isothermal flows where the fluctuation in the reaction rate constant can be neglected. The reaction between mixing air jets containing dilute nitric oxide, NO, and ozone, O₃, is studied in a large Turbulent Smog Chamber (TSC) under conditions of high Reynolds number, three-dimensional flow. The measured reactant concentrations, obtained with good time and space resolution, are used to investigate the effects of the reactant species concentration fluctuations on the mean reactant species concentration field through their contribution to the mean reaction rate. Two flow geometries have been studied: two opposed jets issuing into a large chamber and a plume-like source issuing into a turbulent background flow. A three-dimensional finite-difference computation has been carried out for the flow in the chamber using the k-ε-g model of turbulence and mixing. The reacting species concentration field is calculated by considering the transport of a ‘perturbation variable’ equal to the mean of the difference between the species concentration and its corresponding fast chemistry value. A closure for the mean chemical reaction rate based on this quantity is presented and its experimental validation discussed.     

Removal of lead and zinc ions from aqueous solution using Amasya zeolites from Turkey

A 3-D Eulerian-Lagrangian approach to moving vehicles is presented that takes into account the traffic-induced flow rate and turbulence. The method is applied to pollutant dispersion in an individual street canyon and a system of two street canyons forming a perpendicular intersection. The approach is based on computational fluid dynamics (CFD) calculations using a Eulerian approach for continuous phase and a Lagrangian approach for moving vehicles. The wind speed was assigned values of 4, 7 and 12 m/s. One-way and two-way traffic with different traffic rates per lane is considered. In the case of the intersection, a longitudinal wind direction was assumed. Predictions show differences in the pollutant dispersion in the case of one-way and two-way traffic.     

Numerical simulation of an array of fences in Saemangeum reclaimed land

This paper discusses about the quantitative effect of windbreak fences on wind velocity in the reclaimed land at Saemangeum in South Korea. Windbreak fences were constructed in the reclaimed land to reduce the wind velocity to prevent the generation and diffusion of dust. However, up to this time, no in-depth studies were conducted to quantitatively measure the effect of the windbreak fences on wind velocity thus an optimum windbreak structure is not yet determined. Using CFD simulations, the effects of fence porosity, fence height, and the distance between the adjacent fences were investigated. A wind tunnel experiment was initially conducted and data gathered were used to develop the CFD models. From the experiments and CFD simulations, the overall percentage difference of the measured velocities was 7.20% which is generally acceptable to establishing the reliability of the CFD models. The reduction effect on wind velocity was measured in between the adjacent fences up to a height of 0.6 m from the ground surface. In terms of porosity ( = 0, 0.2, 0.4, 0.6), 0.2 was found to be the optimum value. Conversely, the effect of fence height (0.6, 0.8 and 1.0 m) showed no significant difference; therefore, 0.6 m height is recommended. In addition, the reduction effect of distance between the adjacent fences (2, 4 and 6 m) on wind velocity having a 0.2 porosity has decreased to about 75% regardless of the distance. In the case of the reclaimed land in Saemangeum, a decrease of 75% can prevent the generation and diffusion of dusts. However, the source of dusts is very large. Therefore, constructing an array of windbreak with 6 m distance between them is deemed necessary.