Employing Heisenberg's turbulent spectral transfer theory to parameterize sub-filter scales in LES models

A turbulent subfilter viscosity for large eddy simulation (LES) models is proposed, based on Heisenberg's mechanism of energy transfer. Such viscosity is described in terms of a cutoff wave number, leading to relationships for the grid mesh spacing, for a convective boundary layer (CBL). The limiting wave number represents a sharp filter separating large and small scales of a turbulent flow and, henceforth, Heisenberg's model agrees with the physical foundation of LES models. The comparison between Heisenberg's turbulent viscosity and the classical ones, based on Smagorinsky's parameterization, shows that both procedures lead to similar subgrid exchange coefficients. With this result, the turbulence resolution length scale and the vertical mesh spacing are expressed only in terms of the longitudinal mesh spacing. Through the employment of spectral observational data in the CBL, the mesh spacings, the filter width and the subfilter eddy viscosity are described in terms of the CBL height. The present development shows that Heisenberg's theory naturally establishes a physical criterium that connects the subgrid terms to the large-scale dimensions of the system. The proposed constrain is tested employing a LES code and the results show that it leads to a good representation of the boundary layer variables, without an excessive refinement of the grid mesh.     

Numerical calculation of flow and stack-gas concentration fluctuation around a cubical building

A numerical simulation model was developed to predict the instantaneous concentration fluctuation of a plume and applied to stack-gas diffusion around a cubical building. The flow field, including an instantaneous velocity component, was predicted using the large eddy simulation (LES) method in the developed numerical model. Then, the instantaneous concentration fluctuation was predicted using the obtained unsteady flow field. Concentration was calculated using the finite difference method, in which the LES is expanded for concentration, and the puff method, in which small volumes of the tracer gas are divided and combined according to the calculation mesh sizes. In order to avoid numerical viscous effects, a puff method and finite difference method were applied separately in the regions near and far from the stack-gas release point, respectively. Then, the flow field around a cubical building and the diffusion of stack-gas, emitted from an elevated point source at an upstream position of the building, were calculated using the model mentioned above. Numerical calculation results were compared with those obtained in wind tunnel experiments in which concentration fluctuation was measured using high-response flame ionization detectors. Although there were some discrepancies in the flow field between the calculated results and those of wind tunnel experiments, e.g., the calculated windward length of a cavity region behind the building, the calculated mean velocity and turbulent intensity showed good agreement with those of the wind tunnel experiments. Furthermore, the calculated concentration fluctuation showed good agreement with that in the wind tunnel, not only regarding the features of fluctuating concentration signals, but also statistic quantities, viz., mean concentration, fluctuation intensity and high-concentration values.     

基于大涡模拟的旋风分离器内流场数值模拟研究

为了研究旋风分离器内流场的运动状况,借助Fluent软件,采用大涡模拟,基于曲线坐标系的SIMPLEC算法,对切向入口的Stairmand旋风分离器内流场进行了三维数值模拟研究。数值结果表明,大涡模拟适合于三维强旋流的流场模拟,分离器内部的流动空间可分为内、外2个流动区域,在不同的流动区域中,气体压力、速度场的分布有较大的差异。而且压力分布与速度分布不是绝对的轴对称分布。数值模拟结果对其设计优化研究有一定的应用价值。     

Dry deposition to a uniform canopy: Evaluation of a first-order-closure mathematical model

A mathematical model is developed for the transport of momentum and matter within a canopy consisting of identical elements protruding vertically from a smooth substrate. Turbulent flux is modelled using a mixing-length approach. The loss of momentum (or matter) to individual elements is related to the mean wind speed, and the element-element interaction via the turbulent wind field is represented by a sheltering factor. Careful consideration is given to the formulation of lower boundary conditions. The model assumptions are compared with those of other models.The model predictions are compared with measurements on a vertically- and horizontally-uniform artificial canopy in a wind-tunnel. The model reproduces well the observed relationship between the parameters of the logarithmic wind speed profile above the canopy and the observed deposition velocities of thorium-B (ThB) atoms and particles in the diameter range 0.08–32 μm, using a sheltering factor which is little dependent on wind speed and has the same magnitude for momentum, gas and particles. The predicted dependences of deposition velocity on friction velocity and, for particles, on diameter shed light on the performance of semi-empirical correlations proposed in the literature. For ThB atoms, the calculated deposition velocities are compared with those of other mathematical canopy models: a comparable degree of agreement is obtained here with fewer free parameters.The fraction of deposit on the substrate is underpredicted by an order of magnitude in some cases, pointing to the limitations of the modelling of conditions near the substrate in terms of quasi-shear flow.     

A comparison of mesh-adaptive LES with wind tunnel data for flow past buildings: Mean flows and velocity fluctuations

In this paper we address two important aspects of micro-scale urban airflow model evaluation: (a) the identification of key flow features as dictated by the physics of the problem and as captured by the simulations, and (b) the comparison of important model output parameters (mean flows and fluctuations) with experimental data. A series of mesh-adaptive large eddy simulations (LES) was carried out for the study of air flows within two intersecting street canyons with varying building configurations. The novelty of the approach lies in the combination of LES with mesh adaptivity, which allows a variable-filter length and the implementation of an anisotropic eddy-viscosity model. Both coarse and fine-mesh simulations were carried out, using single and parallel-processor systems respectively. The simulations showed clearly that the expected flow patterns such as the street canyon recirculation and the street-mouth vortices, as well as the exchange of air flow at the street intersections, can readily be captured by the mesh-adaptive LES.In addition, the detailed comparisons of mean flows and fluctuations of the resolved velocity field with the measured data showed that the simulation results agreed well with the patterns and trends of the wind tunnel measurements. In most cases the finer-mesh simulations improved considerably the accuracy of the mean flows, especially for the symmetrical configuration. The improvement in the predicted fluctuations was less obvious, with several detector locations underpredicting the measured values, although the overall comparison was also satisfactory. The typical errors for the mean flows for all three building configurations were less than 30%, whilst for the velocity fluctuations less that 40%. Both the simulated means flows and turbulence levels were generally more accurate in the streets parallel to the wind (streamwise direction) than in the streets normal to the wind.     

The impact of nudging coefficient for the initialization on the atmospheric flow field and the photochemical ozone concentration of Seoul,Korea

In order to incorporate correctly the large or local scale circulation in the model, a nudging term is introduced into the equation of motion. Nudging effects should be included properly in the model to reduce the uncertainties and improve the air flow field. To improve the meteorological components, the nudging coefficient should perform the adequate influence on complex area for the model initialization technique which related to data reliability and error suppression. Several numerical experiments have been undertaken in order to evaluate the effects on air quality modeling by comparing the performance of the meteorological result with variable nudging coefficient experiment. All experiments are calculated by the upper wind conditions (synoptic or asynoptic condition), respectively. Consequently, it is important to examine the model response to nudging effect of wind and mass information. The MM5–CMAQ model was used to assess the ozone differences in each case, during the episode day in Seoul, Korea and we revealed that there were large differences in the ozone concentration for each run.These results suggest that for the appropriate simulation of large or small-scale circulations, nudging considering the synoptic and asynoptic nudging coefficient does have a clear advantage over dynamic initialization, so appropriate limitation of these nudging coefficient values on its upper wind conditions is necessary before making an assessment. The statistical verifications showed that adequate nudging coefficient for both wind and temperature data throughout the model had a consistently positive impact on the atmospheric and air quality field. On the case dominated by large-scale circulation, a large nudging coefficient shows a minor improvement in the atmospheric and air quality field. However, when small-scale convection is present, the large nudging coefficient produces consistent improvement in the atmospheric and air quality field.     

A Simulation Model for Air Pollution over Connecticut

A different approach to mathematically modeling large-scale atmospheric processes is presented. Whereas past approaches have been to develop a model based on an accumulation of information from a specific geographical area, resulting in a model applicable to that area only, we have developed a general mathematical model applicable to any geographical area. The model’s applicability is controlled by specifying the input information describing the meteorological situation and pollution source configuration. A rectangular array of grid points is used to specify both the wind field, by using stream functions, and the average source strength of some pollutant for the area represented by the grid. The diffusion problem is divided into two areas: transport by the mean wind field, and dispersion based on travel time and distance as described by empirical equations. Trajectories of pollutants are traced backwards from the points of interest in the course of the calculations and the contributions of all sources that affect the points of interest are accumulated. The model requires an array of source strength information. An inventory of pollution sources in the State of Connecticut was compiled and maps of source strengths were prepared for five pollutants on a 5000-ft grid-square array. Maps of sulfur dioxide and carbon monoxide source strengths are presented with the resulting concentration distribution for “typical” meteorological conditions. The model permits the changing of meteorological or source values at predetermined intervals so that diurnal changes are incorporated in the calculations. The model has not been verified, but the values of pollution concentration are the right order of magnitude and the resulting patterns are as expected.     

Variable K-theory model for the dispersion of air pollutants in low-wind conditions in the surface-based inversion

A variable K-model has been proposed for the dispersion in low winds in the surface-based inversion by expressing the eddy diffusivities as a linear function of downwind distance from the source. The resulting partial differential equation with variable coefficients along with the physically relevant boundary conditions is solved analytically. For the accuracy of the so-obtained solution, an analysis of the convergence and error estimation has been carried out. It is shown that the series converges absolutely. An upper bound for the error based on the partial sum of the series is estimated and it is described that the error tends to zero as the number of terms in the expansion are sufficiently large.The solution has been used to simulate the field tracer data sets collected from Hanford and IIT diffusion experiments in stable and unstable conditions, respectively. It predicts 41% cases in stable and 35% cases in unstable conditions within a factor of two to observations.     

Impact of using prognostic and objective wind fields on the photochemical modeling of Athens, Greece

A three-dimensional Eulerian photochemical model is used to follow the dynamics of ozone, NO_x, and CO over the Athens area, for 25 May 1990, the day considered in the APSIS project. A unique aspect of this work lies in the study of the impacts of the wind field preparation methods on the concentrations predicted by the model. Three sets of wind fields are developed. The first one used is derived from a prognostic meteorological model. The second one is calculated from available wind observations using objective: methods. For these two cases, a previous day is simulated, using the same conditions, to develop preconditioned initial conditions for the following day. For the third simulation, again two days are simulated, this time using the observed winds for each of the two days modeled. The predictions using the prognostically derived and the objective analysis wind fields are significantly different, particularly for the primary pollutants. Comparing predictions to the observations did not favor any particular method of wind field preparation. In this case, when using the prognostically derived field, the simulations are very sensitive to boundary conditions. In contrast, when using the wind fields constructed by objective methods, the simulations became most sensitive to emissions and initial conditions. This comes directly from the different residence times in the domain, which are governed by the wind speed.     

A hybrid plume model for local-scale atmospheric dispersion

This paper presents a local-scale dispersion model, based on atmospheric boundary layer scaling theory. In the vicinity of the source, Gaussian equations are used in both the horizontal and vertical directions. After a specified transition distance, gradient transfer theory is applied in the vertical direction, while the horizontal dispersion is still assumed to be Gaussian. The dispersion parameters and eddy diffusivity are modelled in a form, which facilitates the use of a meteorological pre-processor. We present a novel model of the vertical eddy diffusivity (K_z), which is a continuous function of height in various atmospheric scaling regions. The model also includes a treatment of the dry deposition of gases and particulate matter. The accuracy of the numerical model was analysed by comparing the model predictions with two analytical solutions; the numerical deviations from these solutions were less than 2% for the computational regime. The model has been tested against the Kincaid experimental field data. The agreement of the predictions and the data is good on the average, although the internal variation of the predictions versus data scatter plot is substantial.     

An analytical model for crosswind integrated concentrations released from a continuous source in a finite atmospheric boundary layer

An analytical model for the crosswind integrated concentrations released from a continuous source in a finite atmospheric boundary layer is formulated by considering the wind speed as a power law profile of vertical height above the ground and eddy diffusivity as an explicit function of downwind distance from the source and vertical height. A closed form analytical solution of the resulting advection–diffusion equation for these profiles of wind speed and eddy diffusivity with the physically relevant boundary conditions is derived using the separation of variables technique that leads to a Sturm–Liouville eigen value problem. Various particular cases of the model are deduced.The model is evaluated with the observations obtained from Prairie Grass experiment in various stability classes varying from very unstable to neutral and stable conditions and Hanford diffusion experiment in stable conditions. The agreement is found to be good between the computed and observed concentrations in both the diffusion experiments. For Prairie Grass experiment, the model is predicting 78% cases with in a factor of two and gives a value of NMSE as 0.075. On the other hand, for Hanford observations in stable conditions, it predicts 70% cases with in factor of two. An extensive analysis of statistical measures with the downwind distances from the source reveals that the model is performing well close to the source.     

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.     

KAPPA-G,A Non-Gaussian Plume Dispersion Model: Description and Evaluation Against Tracer Measurements

This paper presents a mixed methodology for the simulation of atmospheric disperson phenomena in which vertical diffusion is computed using an analytical solution of the K-theory equation, while horizontal diffusion is simulated by the Gaussian formula. This new formulation, while maintaining a simple analytical form for the concentration field, incorporates the effects of power-law vertical profiles of both wind speed and eddy diffusivity. The performance of this approach, which has been implemented into a full computer package (KAPPA-G), is evaluated by comparison with data from SF₆ tracer experiments.     

Turbulent diffusion modelling FQR wind flow and dispersion analysis

Within the framework of searching for relatively simple, reliable and universal eddy viscosity /diffusivity models, a new three dimensional general non-isotropic model is proposed applicable to any domain complexity and any atmospheric stability conditions. The model utilizes the transport equation for turbulent kinetic energy but introduces a new approach in effective length scale estimation based on the flow global characteristics and local atmospheric stability. The model is discussed in detail and predictions are given for flow field and boundary layer thickness. The results are compared with experimental data with satisfactory results.     

A Wind Tunnel Study to Design Large,Open-top Chambers for Whole-tree Pollutant Exposure Experiments

A wind tunnel study was conducted to determine the optimal design features of a large, open-top chamber, as needed for pollution exposure studies on mature trees. An optimally-designed, open-top chamber must minimize the incursion of ambient air through its opening and maintain a uniform treatment concentration throughout the chamber. The design features of interest are the diameter and height of the chamber and the deflection angle and opening size of any frustum that may be mounted on top of a model chamber.

Design specifications depend on the turbulence regime about the chamber, which is influenced by the nature of the surrounding vegetation. Consequently, our investigation was performed on scale-model, open-top chambers in a wind tunnel populated with a model coniferous forest. Turbulence measurements demonstrated the similarity between the turbulence regime of the model and a natural forest. A hydrocarbon tracer was injected into the wind tunnel flow to characterize chamber performance.

The main design features of open-top chambers are the velocity of air exiting through the top and the relationship between the length scale of the turbulence and the diameter of the chamber opening. As exit velocities increase, the proportion of eddies with sufficient force to penetrate into the chamber decrease. Therefore, for equal volumetric air flows, smaller opening sizes increase the exit velocities and reduce the number and extent of ambient air incursions. Almost total exclusion of ambient air is achieved as the exit velocity of the air exceeds the magnitude of one standard deviation of the vertical wind velocity measured at the chamber top. The incursion of ambient air is also reduced when the diameter of the chamber opening is smaller than the characteristic length scale of the turbulence, a measure of mean eddy size.

Frusta deflect air flow over the chamber. Three prototypes, with 30?, 45? and 60-degree angles were tested. A 30-degree frustum slightly improves the performance of the chamber and is more effective in preventing ambient air from entraining into the chamber opening than frusta with either a 45? or 60-degree angle. A flatter frustum allows for a smoother transition in the wind velocity streamline and is less apt to cause wake turbulence, as is the case with steeper frusta.

Knowledge of the turbulence characteristics of plant canopies are readily available in the literature and can aid scientists and engineers in designing the optimal chamber and frusta dimensions for their particular application. Therefore, the empirical approach to chamber design can be avoided, and substantial savings can be realized.     

Numerical simulation of flow and dispersion around an isolated cubical building: The effect of the atmospheric stratification

The aim of this work is to investigate atmospheric flow and dispersion of contaminants in the vicinity of single buildings under different stability conditions. The mathematical model used is based on the solution of equations of conservation of mass, linear momentum and energy with the use of a non-standard κ–? turbulence model. The modifications proposed in the κ–? model are the inclusion of the Kato and Launder correction in the production of turbulent kinetic energy and the use of a modified wall function. Results are presented of numerical simulations of dispersion around a cubical obstacle, under neutral, stable and unstable atmospheric conditions. Experimental data from wind tunnel and field trials obtained by previous authors are used to validate the numerical results. The numerical simulation results show a reasonable level of agreement with field and wind tunnel concentration data. The deviation between model results and field experimental data is of the same order as the deviation between field and wind tunnel data.     

Modelling the transport of suspended particulate matter by the Rhone River plume (France). Implications for pollutant dispersion

A model to simulate the transport of suspended particulate matter by the Rhone River plume has been developed. The model solves the 3D hydrodynamic equations, including baroclinic terms and a 1-equation turbulence model, and the suspended matter equations including advection/diffusion of particles, settling and deposition. Four particle classes are considered simultaneously according to observations in the Rhone. Computed currents, salinity and particle distributions are, in general, in good agreement with observations or previous calculations. The model also provides sedimentation rates and the distribution of different particle classes over the sea bed. It has been found that high sedimentation rates close to the river mouth are due to coarse particles that sink rapidly. Computed sedimentation rates are also similar to those derived from observations. The model has been applied to simulate the transport of radionuclides by the plume, since suspended matter is the main vector for them. The radionuclide transport model, previously described and validated, includes exchanges of radionuclides between water, suspended matter and bottom sediment described in terms of kinetic rates. A new feature is the explicit inclusion of the dependence of kinetic rates upon salinity. The model has been applied to 137Cs and 239,240Pu. Results are, in general, in good agreement with observations.     

Incinerator Conference Covers Emissions Control Aspects

The success of the application of computer modeling to decision-making will depend on the degree to which the scientifically valid “cause-and-effect” features of the air pollution system are represented. For this reason, dynamic simulation models are to be preferred to statistical and empirical models. A digital simulation model based on a stoichiometrically logical chemical mechanism and trajectory estimating routines was constructed, using Los Angeles source, meteorological and geographic input. The basic physical concept underlying the simulation model is the process of evolution of photochemical pollution in a parcel of air as it moves in a dynamic urban emission/meteorological environment along a given urban wind trajectory. Both the photochemical evolution and the trajectory are numerically integrated by a standard linear multistep predictor-corrector method. Concentrations of photochemical reactants and products (i.e., primary and secondary contaminants) are determined by this numerical integration, which also includes appropriate terms for relevant effects. In five preliminary validation runs, simulated NO₂, NO, and O₃ values were within 20% or 0.05 ppm of those observed at air monitoring stations located near the termini of the runs. The trajectories were plotted on the basis of hourly meteorological data for 22 stations. Six control strategy exercises were conducted to illustrate the application of the model to problem-solving situations.     

Validation of the lagrangian particle dispersion model FLEXPART against large-scale tracer experiment data

A comprehensive validation of FLEXPART, a recently developed Lagrangian particle dispersion model based on meteorological data from the European Centre for Medium-Range Weather Forecasts, is described in this paper. Measurement data from three large-scale tracer experiments, the Cross-Appalachian Tracer Experiment (CAPTEX), the Across North America Tracer Experiment (ANATEX) and the European Tracer Experiment (ETEX) are used for this purpose. The evaluation is based entirely on comparisons of model results and measurements paired in space and time. It is found that some of the statistical parameters often used for model validation are extremely sensitive to small measurement errors and should not be used in future studies. 40 cases of tracer dispersion are studied, allowing a validation of the model performance under a variety of different meteorological conditions. The model usually performs very well under undisturbed meteorological conditions, but it is less skilful in the presence of fronts. The two ETEX cases reveal the full range of the model’s skill, with the first one being among the best cases studied, and the second one being, by far, the worst. The model performance in terms of the statistical parameters used stays rather constant with time over the periods (up to 117 h) studied here. It is shown that the method used to estimate the concentrations at the receptor locations has a significant effect on the evaluation results. The vertical wind component sometimes has a large influence on the model results, but on the average only a slight improvement over simulations which neglect the vertical wind can be demonstrated. Subgrid variability of mixing heights is important and must be accounted for.     

Simulations of flow around a cubical building: comparison with towing-tank data and assessment of radiatively induced thermal effects

A three-dimensional (3-D) computational fluid dynamics (CFD) model, coupled with a meteorological radiation and surface physics package, is used to model the mean flow field and tracer dispersion in the vicinity of an idealized cubical building. We first compare the simulations with earlier numerical studies as well as towing-tank laboratory experiments, where radiation effects were not included. Our simulations capture most of the features revealed by the towing-tank data, including the variation of the flow reattachment point as a function of Froude number and the induction of a prominent lee wave in the low Froude number regime. The simulated tracer concentration also compares very favorably with the data.We then assess the thermal effects due to radiative heating on the ground and building including shading by the building, on the mean flow and tracer dispersion. Our simulations show that convergence within and beyond the cavity zone causes a substantial lofting of the air mass downstream from the building. This lofting results from the combination of thermal heating of the ground and building roof, and vortex circulation associated with the horseshoe eddy along the lateral sides of the building. The specific effect of shading on the flow field is isolated by comparing simulations for which the radiative heating and shading patterns are kept constant, but the environmental wind direction is altered. It is found that the shading exerts local cooling, which can be combined into the overall thermodynamic interaction, described above, to effectively alter the circulation downstream from the building.