**Chaojie Li, Liyan Liu , Wei Tan **

School of Chemical Engineering and Technology, Tianjin University, Tianjin 300350, China

Received:
September 2, 2018

Revised:
December 13, 2018

Accepted:
January 16, 2019

Download Citation:
||https://doi.org/10.4209/aaqr.2018.09.0328

Cite this article:

Li, C., Liu, L. and Tan, W. (2019). Evaluation of RSM for Simulating Dispersion of CO2 Cloud in Flat and Urban Terrains. *Aerosol Air Qual. Res.* 19: 390-398. https://doi.org/10.4209/aaqr.2018.09.0328

**HIGHLIGHTS**

- Simulations of CO
_{2}dispersion with stress-*ω*and SST*k*-*ω*models are conducted. - Simulation results are compared with experimental data for evaluation of RSM.
- Stress-
*ω*model exhibits better performance in flat terrain. - Performance of stress-ω model is worse than that of SST
*k*-*ω*model in urban terrain.

**ABSTRACT**

Accurate and reliable computational fluid dynamics (CFD) simulation of pollutant dispersion is essential for protecting human health, and the choice of turbulence model is an important parameter determining the accuracy of simulation results. This paper evaluates the ability of Reyonds stress model (RSM) to predict dispersion of carbon dioxide (CO_{2}) cloud, which is a typical type of heavy gas and similar to some particulate pollutants, in flat and urban terrains. The RSM simulation is conducted with stress-*ω* model, whereas SST *k*-*ω* two-equation model is selected as the benchmark. The simulation results are compared with the available wind tunnel measurements, and statistical performance indicators are used to obtain a comprehensive and quantitative evaluation of the performances of the two turbulence models. The results reveal that stress-ω model exhibits different capacities in flat terrain and urban terrain. Specifically, stress-ω model can present better results than SST *k*-*ω* model in flat terrain, and it performs better in the far-field region than in the near-field region. Although SST k-*ω* model can describe CO_{2} dispersion more accurately in urban terrain, the concentration distribution reproduced by stress-*ω* model is still within acceptable range.

Keywords:
Flat terrain; Urban terrain; Stress-ω; SST k-ω.

**I****NTRODUCTION**

Dispersion of pollutant represents a major threat to human health, and the precise prediction of pollutant concentration distribution is a prerequisitefor avoiding adverse air quality impacts (Zhang *et al**.*, 2011; Dodla *et al.*, 2017; Guttikunda and Jawahar, 2018). Compared with the simple Gaussian and other analytical dispersion models, computational fluid dynamics (CFD) models can provide results of the flow feature in complicated terrain simultaneously and have been widely employed to investigate pollutant dispersion (Gupta and Chan, 2016; Lateb *et al.*, 2016). Lin *et al.* (2016)studied the effects of street layout and wall heating on pollutant dispersion by CFD simulation. Huang *et al.* (2017) applied a CFD model toevaluate the impacts of the upstream building height and stack location on the dispersion of pollutant.

However, CFD models need to be established and evaluated correctly to accurately simulate pollutant dispersion. The selection of turbulence model is one of the prerequisites for CFD simulation, as it greatly affects the reproduction of flow field and the prediction of pollutant dispersion. Among the various turbulence models, the two-equation models are the main approach for turbulence modeling of the practical pollutant dispersion. However, these models are based on the isotropic Boussinesq hypothesis (Lane, 2017), which can make these models lack of the ability to predict accurately the complex flow structure along the ground, thus causing the results to be unreliable(Franke *et al.*, 2005; Longo *et al.*, 2017). To avoid the assumption of isotropy, the Reynolds stress model (RSM) can be used, which in theory should perform better than two-equation models. Yet, in some cases, the converged solution is far beyond satisfaction, and the approach cannot yield any notable gains in accuracy (Wang and McNamara, 2006; Tominaga and Stathopoulos, 2013; Peren *et al.*, 2015). For instance, Sklavounos and Rigas (2004) evaluated four turbulence models (i.e., standard *k*-*ε* model, standard *k*-*ω* model, SST *k*-*ω* model, and Speziale-Sarkar-Gatski (SSG) RSM) against Thorney Island field trials data, and found that SSG model prolonged calculation time with only constrained gains in accuracy. Bartzanas *et al.* (2007) tested four turbulence models (i.e., standard *k*-*ε* model, RNG *k*-*ε* model, realizable *k*-*ε* model, RSM) for the ventilated airflow, and found that RSM needed long convergence time and showed better performance only under low air velocities.

However, careful inspection of the relevant studies, makes clear that solid conclusion cannot be drawn by simply evaluating the results of the limited number of cases. Therefore, additional research efforts are necessary to further evaluate the performance of RSM in certain cases. In light of this, the present research aims to test the ability of RSM to predict pollutant dispersion in two types of terrains. The first type is the flat terrain, where a relatively simple flow field is observed, while the second is the urban terrain, where the complicated flow structure is produced mainly due to geometric effects (Carpentieri *et al.*, 2015).

SST *k*-*ω* model is characterized by an advanced near-wall treatment, taking advantages of both the *k*-*ε* and *k*-*ω* models. It gets rid of the drawback of the standard *k*-*ε* model that unable to describe the near-wall zone, the drawback of the modified *k*-*ε* models that cannot capture information in the wake volume, and the drawback of the *k*-*ω* model that be sensitive to free stream (Li *et al**.*, 2006; Lateb *et al**.*, 2013; Yu and The, 2016). Thus, SST *k*-*ω* model is used as the reference point for the comparison of the results. Stress-*ω* model, one of the RSM type, uses the similar boundary condition for *ω* (Sharif and Guo, 2007). Thus, both models employ the same mesh densities and near-wall grid sizes, and the results are simply determined by the turbulence models themselves. Numerical simulation of CO_{2} cloud dispersion in flat and urban terrains were carried out by means of SST *k*-*ω* model and stress-*ω* model, and the results were tested against the available experimental data obtained from the literatures (Xing *et al.*, 2013; Tan *et al.*, 2018). Also, CO_{2} cloud is a typical type of heavy gas, and therefore the results can be conducive to the simulation of some particulate pollutants dispersion.

TURBULENCE MODEL

TURBULENCE MODEL

SST k-ω Model

SST k-ω Model

SST k-ω Model

SST *k*-*ω* turbulence model is a two-equation eddy-viscosity model based on Boussinesq eddy viscosity assumption, taking both the advantage of *k*-*ω* model that performs better in the near-wall region and the advantage of standard *k*-*ε* model that is accurate and reliable for the flow outside the boundary layer, with a limiter evaluating the proper transport behavior of turbulent shear stress and a blending function transiting between *k*-*ω*model and *k*-*ε* model on the basis of the distance to the surface and on the flow variables. It has been proved to be able to predict flow field and pollutant dispersion in the atmosphere nicely (Eca and Hoekstra, 2011; Yu and The, 2016). SST is a good benchmark for higher fidelity models such as RSM, to determine if they can accurately predict the flow field distribution and pollutant dispersion (O’Brien *et al.*, 2018). Equations of the model are presented in the following:* *

*k*-equation

*ω*-equation

The kinematic eddy viscosity is calculated from

where *S* is the strain rate magnitude.

The closure coefficients and auxiliary relations are given by

where *y* is the distance to the next surface.

The constants are summarized as *ϕ* = *ϕ*_{1}*F*_{1} + *ϕ*_{2}(1 – *F*_{1}), with *ϕ* = *σ _{k}*,

*σ*,

_{ω}*α*,

*β*,

*a*

_{1}= 0.31,

*α*

_{1}= 5/9,

*α*

_{2}= 0.44,

*β*

_{1}= 3/40,

*β*

_{2}= 0.0828,

*β*0.09,

^{*}=*σ*

_{k}_{1}= 0.85,

*σ*

_{k}_{2}= 1,

*σ*

_{ω}_{1}= 0.5,

*σ*

_{ω}_{2}= 0.856.

Stress-ω* Model*

Stress-ω

Stress-ω

*Model*

RSM (Launder *et al.*, 1975) resolves directly a transport equation for each term of the Reynolds stress tensor separately using differential equations, taking the anisotropy of turbulence into account and making up for the deficiencies of the Boussinesq approximation. According to the modeling of the pressure-strain and dissipation-rate terms, the CFD-Fluent provides numerous types of RSM. Stress-*ω* model, a stress-transport model based on the omega equation and Launder-Reece-Rodi (LRR) model (Launder *et al.*, 1975), is considered in this paper. Stress-*ω* model is similar to *k*-*ω* model, mainly reflected in the transport equation for *ω* and the wall boundary condition for *ω *(Sharif and Guo, 2007)*.* The model equations are expressed as follows:

Stress-equation

*ω*-equation

The pressure-strain correlation is given as

where

The closure coefficients are defined as *C*_{1} = 9/5, *C*_{2} = 13/25, *α* = 13/25, *σ* = 1/2, *σ*^{*} = 1/2, *β*_{0} = 9/125,* β*^{*}_{0} = 9/100, *β *= *β*_{0}*f _{β}*,

*β*

^{*}=

*β*

^{*}_{0}

*f*, , ,.

_{β*}

STATISTICAL MODEL PERFORMANCE EVALUATION METHOD

STATISTICAL MODEL PERFORMANCE EVALUATION METHOD

Five performance measures (Hanna *et al.*, 2004), including the fractional bias (FB), the geometric mean bias (MG), the normalized mean square error (NMSE), the geometric variance (VG), and the fraction of predictions within a factor of two of observations (FAC2), defined as below, were used to obtain a comprehensive and quantitative evaluation of the performances of these two turbulence models for dispersion of CO_{2} cloud.

where *C _{p}* was model predictions of concentration,

*C*was observations of concentration, and

_{o}*C̅*was average over the data set.

According to the literature (Hanna *et al.*, 2004), FB = 0, MG = 1.0, NMSE = 0, VG = 1.0 and FAC2 = 1.0 represent the perfect model, and –0.3 < FB < 0.3, 0.7 < MG < 1.3, NMSE < 4, VG < 1.6 and 0.5 < FAC2 represent the acceptable range.

EVALUATION OF STRESS-*ω*** MODEL IN FLAT TERRAIN**

EVALUATION OF STRESS-

*ω*

Description of Wind Tunnel Experiment Considered

Description of Wind Tunnel Experiment Considered

Description of Wind Tunnel Experiment Considered

This experiment (Xing *et al.*, 2013) conducted by the Beijing Institute of Technology was simulated, and the simulation results were compared with experimental data to conduct the evaluation of stress-*ω* model for simulating dispersion of CO_{2} cloud in flat terrain. The experiment was implemented in a 6.4-m-wide, 5.2-m-high wind tunnel located in Zhangjiakou. The sonic anemometer was placed at the wind inlet at a height of 2 m to measure the wind speed. In the trials, the CO_{2} inlet had a radius of 0.01 m and CO_{2} of 99.9% by volume was oriented vertically upwards under atmospheric pressure, which could realize the simulation of pollutant dispersion from a point release at the ground level (Shen *et al**.*, 2017; Wingstedt *et al**.*, 2017). Four different cases with CO_{2} flow rates of 166.7 L min^{–1}, 200 L min^{–1}, 250 L min^{–1}, 300 L min^{–1}, corresponding to the flow speeds of 8.8 m s^{–1}, 10.6 m s^{–1}, 13.3 m s^{–1} and 15.9 m s^{–1}, were carried out. Ten sensors were arrayed in the central line through the jet nozzle along the downwind direction, with the distances of 0.5 m, 1 m, 1.5 m, 2 m, 2.5 m, 3 m, 4 m, 6 m, 8 m, 10 m, to acquire CO_{2} concentration data.

###

*Fluent CFD Model Set-up*

*Fluent CFD Model Set-up*

The rectangular domain with a cross-section of 6.4 m × 5.2 m was developed in this work, as shown in Fig. 1. The front of the domain was air inlet, and the rear, named outlet, remained open to the atmosphere. The wind direction was parallel to the bottom from the front to the rear, and a fully developed flow can be established in the computational domain. Structural grid was adopted to discrete the computational domain. Three quantities of cells (i.e., 1133136, 1890864, and 3828672) were utilized to perform the grid sensitivity analysis. SST *k*-*ω* model was chosen here. CO_{2} concentration in the central ground at 200 s after release with the flow speed of 8.8 m s^{–1} was compared, as shown in Fig. 2. It was found that the change in the results was limited and the grid with 1890864 cells was chosen. The minimum horizontal grid size was about 0.0025 m at CO_{2}inlet. The minimum grid size in the vertical direction was 0.02 m near the ground and expanding away from it with an expansion factor of 1.1 to 0.08 m.

Fig**. 1.** Computational grid for flat terrain.**Fig. 2.** Grid independence analysis for flat terrain.

Difference between the wind condition in the experiment and in the simulation can also lead to the deviation between the experimental and simulated cloud concentration (Garcia-Sanchez *et al.*, 2014; Weerasuriya *et al.*, 2018; Zhang *et al.*, 2018). The inlet boundary condition in the CFD simulation was obtained considering the local atmospheric stability class of “*F*” at the field site. The wind velocity profile was defined by the power law (Tominaga *et al.*, 2008; Wu *et al**.*, 2014) to describe the boundary layer.

where, the reference velocity *U*_{0} and the reference height *Z*_{0} were respectively 0.37 m s^{–}^{1} and 2 m, and the value of the power-law exponent *α* was set as 0.3 according to the literature (Xing *et al.*, 2013). The power law profile was included into the FLUENT code through user-defined functions (UDFs). Air temperature and CO_{2} temperature

were obtained based on experimental data. The initial and boundary conditions are summarized in Table 1.

Comparison of CO_{2} Concentration with Various Turbulence Models

Comparison of CO

Comparison of CO

_{2}Concentration with Various Turbulence ModelsThe two cases with CO_{2} flow speeds of 8.8 m s^{–1} (166.7 L min^{–1}) and 13.3 m s^{–1} (250 L min^{–1}) were modeled, respectively. The simulation strategy involved a steady-state simulation for flow field and a transient simulation for CO_{2} cloud dispersion, with the time step of 0.05 s. CO_{2}concentration data in the central ground line parallel to the wind velocity at 200 s, when the flow field has been relatively steady, were measured. Figs. 3 and 4 show the comparisons between the experimental data and the simulation results using various turbulence models.

Fig**. ****3****. **CO_{2} concentration in the central ground line parallel to wind flow with CO_{2} flow speed of 8.8 m s^{–}^{1}.**Fig****. ****4****.** CO_{2} concentration in the central ground line parallel to wind flow with CO_{2} flow speed of 13.3 m s^{–}^{1}.

The comparisons between experimental data and simulation results show that both of the two turbulence models had a weak accuracy in the near-field region, where the source flow showed strong interaction with wind flow. With the vertical momentum of the source, vortices can be formed in both the upstream and downstream of CO_{2} inlet, and the flow in the near-field region was dominated by the vortices (Wingstedt *et al**.*, 2017). This created highly anisotropic turbulent structure. For the case with CO_{2} flow speed of 8.8 m s^{–1}, stress-*ω* model provided the better results in both near- and far-field regions in comparison to SST *k*-*ω* model, although SST *k*-*ω* model captured the maximum value and the position corresponding to the maximum value. Especially in the near-field region, SST *k*-*ω* model significantly underestimated CO_{2} concentration and the advantage of stress-*ω* model was more pronounced. For the case with CO_{2} flow speed of 13.3 m s^{–1}, the simulation results from stress-*ω* model were in good agreement with the experimental data with respect to both the magnitude and trend. It successfully captured the position corresponding to the maximum value and the variation of CO_{2} concentration in the near-field region that decreased firstly then increased with distance. Whereas, the results predicted by SST *k*-*ω* model exhibited different law that CO_{2} concentration increased with distance. Also, the position of maximum point appeared a certain degree of deviation.

As shown in Table 2, all five performance measures for stress-*ω* model were within acceptable range. While, some measures for SST *k*-*ω* model were outside, as highlighted in gray. And the main reason is that the strong interaction between wind flow and source flow creates highly anisotropic turbulent structure in the near-field region, and SST *k*-*ω* model lacks the capacity to predict accurately the complex flow structure, causing the larger deviation between simulation result and experimental data. The results of the validation metrics support that stress-*ω* model can produce much better results compared with SST *k*-*ω* model in flat terrain. Therefore, the enhancement of accuracy of simulation results can be considered to be worth the increased computational time.

EVALUATION OF STRESS-*ω*** MODEL IN URBAN TERRAIN**

EVALUATION OF STRESS-

*ω*

Description of Wind Tunnel Experiment Considered

Description of Wind Tunnel Experiment Considered

Description of Wind Tunnel Experiment Considered

In this work, urban terrain was referred to as the street canyon, which was the basic element of urban fabric. In the street canyon, cloud concentration field was mainly governed by the wind field induced by the geometric effects, consisting of wakes, stagnating zones, recirculation, and preferential paths (Park *et al.*, 2015; Mei* et al.*, 2016; Tominaga and Stathopoulos, 2016; Nosek *et al.*, 2017). The experiment (Tan *et al.*, 2018) focused on CO_{2} cloud dispersion in the street canyon was simulated for the evaluation of stress-*ω* model in cloud dispersion in urban terrain. This experiment was carried out in a 0.75-m-wide, 0.5-m-high wind tunnel. Spires and roughness elements were installed in front of the building arrays to develop the boundary layer. The wind direction was perpendicular to the target street canyon. The incident velocity was measured in the wind tunnel with anemometer at the middle line upstream of the test section to calculate mean velocity and turbulence intensity profile. The reference mean velocity at the building height was 0.66 m s^{–1} and the exponent α in the power law profile was 0.308. The street canyons were modeled using box-shaped blocks and the CO_{2} inlet was a circle with a radius of 0.0065 m located in the target canyon. CO_{2} of 99.9% by volume was jetted vertically upwards at ground level with the flow rate of 4.48 L min^{–1} under atmospheric pressure, corresponding to the flow speed of 0.5628 m s^{–1}. These two box-shaped blocks used in the experiment had dimensions of 0.2 × 0.05 × 0.05 m^{3} (L × W × H) and 0.05 × 0.05 × 0.05 m^{3} (L × W × H) respectively, forming five streets connected through a streamwise street. Four sensors based on Infra-Red and detector technology were respectively located near the CO_{2} inlet, on the leeward side, on the windward side and above the windward building to acquire CO_{2} concentration data. Considering CO_{2} inlet as the origin with *y* measured as positive in the downwind direction and *z* positive upwards, the four measuring points’ coordinates were respectively Point 1 (0, –0.025, 0.025), Point 2 (0.125, –0.025, 0.025), Point 3 (0.125, 0.025, 0.025), and Point 4 (0, 0.05, 0.05).

Fluent CFD Model Set-up

Fluent CFD Model Set-up

Fluent CFD Model Set-up

The size of the domain was determined following the wind-tunnel geometry, with width 0.75 m, and height 0.5 m, as shown in Fig. 5. A fully structured grid was built. SST *k*-*ω* model was also chosen to perform grid-sensitivity analysis among three types of grids, with the number of1380774, 3183462 and 6135432. Fig. 6 shows the influence

of grid size on CO_{2} concentration in the centerline of the street at z = 0.025 m and at 15 s after release. The grid with 3,183,462 cells was selectedkeeping in view both the computing time and the solution accuracy. The grid size in the canyon was 0.005 m and expanding away from it to 0.01 m. Further refinement close to walls and ground was performed.

**Fig****. ****5****. **Computational grid for urban terrain.**Fig. 6.** Grid independence analysis for urban terrain.

The inlet boundary conditions used in the simulation were based on the measured incident vertical profiles of mean wind velocity and turbulence intensity in the wind tunnel experiment. The vertical profiles of mean wind velocity, turbulent kinetic energy and its specific dissipation rate were used to describe the variation of wind as a function of height.

where *v̅*(*z*) was the mean velocity at z.

where *I* was the turbulence intensity.

where *L* was turbulence integral scale.

According to the experiment, CO_{2} was released at an average speed of 0.5628 m s^{–1}. The initial and boundary conditions used in the simulation are shown in Table 3.

Comparison of CO_{2} Concentration with Various Turbulence Models

Comparison of CO

Comparison of CO

_{2}Concentration with Various Turbulence ModelsThe time step was 1 × 10^{–3} s in the transient simulation to assure Courant number was less than 1. The calculation was carried out up to *t* =15 s, by which time CO_{2} concentration had reached a relative stable value. CO_{2} concentration data at the measuring points were calculated. As shown in Fig. 7, CO_{2} concentration in urban terrain was larger than in flat terrain due to the blocking effect of spreading downstream. Simulation results reproduced by stress-*ω* model agreed worse with experimental data at all the four measuring points. It underestimated CO_{2} concentration near the CO_{2} inlet and on the windward side, and overestimated the concentration on the leeward side and above the windward building. However, SST *k*-*ω*model underestimated the concentration above the windward building, overestimated the concentration on both the leeward side and windward building, and had the closest coincidence with experimental data near the CO_{2} inlet.

**Fig****. ****7****. **CO_{2} concentration at the measuring points.

As shown in Table 4, the simulation results from the two turbulence models were all within acceptable range. And the results also reveal that the performance of stress-*ω* model was worse than that of SST *k*-*ω* model. This may be due to two reasons. The first is that stress-*ω* model requires more empirical data which can vary from one flow to another, and the advantage itself for complex turbulent flow in urban terrain may be not significant. The second is that an advanced near-wall treatment may play a more important role than the consideration of turbulence anisotropy. In other words, the advantage of SST *k*-*ω* that automatically switches between the low- and high-Re number *ω* formulation is more significant for the prediction accuracy of flow structure in urban terrain. Coupled with the extra computational effort for the higher number of differential equations, stress-*ω* model still needs further investigation and optimization for accounting for various effects of complex turbulent flow in urban terrain.

C**ONCLUSIONS**

C

This paper performs a comparison of CO_{2} concentration data obtained by SST *k*-*ω* model (as the benchmark) and stress-*ω *model with the available wind tunnel test data, to evaluate the applicability and accuracy of RSM for simulating cloud dispersion in flat terrain and urban terrain. The results reflect that stress-*ω* model shows better agreement with experimental data compared to SST *k*-*ω* model in flat terrain. And it performs better in the far field than in the near-field region. The performance measures for stress-*ω* model are all within acceptable range, while some measures for SST *k*-*ω* model are outside. For urban terrain, the performance of stress-*ω *turbulence model is worse than that of SST *k*-*ω* model, although the simulation results are still within acceptable range. This leads to the conclusion that the capacity of RSM is concerned with its application scenario. Specially, stress-*ω* model exhibits better performance for simulation of the relatively simple flow field in flat terrain and SST *k*-*ω* model is still ascendant for the complicated flow structure in urban terrain. The conclusions are of primary importance when deciding which turbulence model to use and drafting guidelines on the use of CFD for the simulation of pollutant dispersion in flat and urban terrains.

It is recommended that the comparison of two-equation model and RSM be made in different terrain conditions. Also, additional experimental data examining the flow structure are necessary to provide a quantitative explanation for the limitation of turbulence models.** **

A**CKNOWLEDGMENTS**

A

This work was supported by The National Key Research and Development Program of China (grant number 2016YFC0801900); Tianjin Science and Technology Planning Project (grant number 15ZCZDSF00550).

NOMENCLATURE

NOMENCLATURE

*k* Turbulent kinetic energy*t* Time*ε** *Turbulent dissipation rate*ω* Specific turbulent dissipation rate*δ** _{ij}* Kronecker delta function

*δ*Boundary layer thickness

*v*Turbulence dynamic viscosity

*v*Turbulence kinematic viscosity

_{t}*ρ*

*Density*

*τ*

*Reynolds stress tensor*

_{ij}*S*The strain rate magnitude

*P*The production limiter to prevent turbulence buildup in stagnation regions

_{k}*U*Velocity vector

*i*,

*j*Grid node in

*x*and

*y*direction

Ω

*The mean-strain-rate tensor*

_{ij}*S*The mean-rotation-tensor

_{ki}*Aerosol Air Qual. Res. *19
:390
-398
. https://doi.org/10.4209/aaqr.2018.09.0328