Journal of Applied Mathematics

Volume 2014 (2014), Article ID 839726, 11 pages

http://dx.doi.org/10.1155/2014/839726

Research Article

Kinetic-Monte-Carlo-Based Parallel Evolution Simulation Algorithm of Dust Particles

Xiaomei Hu, Zhifeng Xu, Hongxia Cai, and Junjun Hu

Shanghai Key Laboratory of Mechanical Automation and Robotics, School of Mechatronic Engineering and Automation, Shanghai University, Mailbox 232, No. 149 Yanchang Road, Shanghai 200072, China

Received 24 September 2013; Accepted 24 December 2013; Published 21 January 2014

Academic Editor: Bernard J. Geurts

Copyright © 2014 Xiaomei Hu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The evolution simulation of dust particles provides an important way to analyze the impact of dust on the environment. KMC-based parallel algorithm is proposed to simulate the evolution of dust particles. In the parallel evolution simulation algorithm of dust particles, data distribution way and communication optimizing strategy are raised to balance the load of every process and reduce the communication expense among processes. The experimental results show that the simulation of diffusion, sediment, and resuspension of dust particles in virtual campus is realized and the simulation time is shortened by parallel algorithm, which makes up for the shortage of serial computing and makes the simulation of large-scale virtual environment possible.

1. Introduction

Many ecological environmental problems have emerged during the process of urbanization. Sedimentation of lots of surface dust in cities caused by transportation is one of them. Dust has close relation with particulates in the atmosphere [1]. Ecological system can be hurt by dust covertly over a long period, while it is human body that dust can do most direct and greatest harm to. Dust (especially suspended particulate matter with aerodynamic diameters less than 10 μm) has been one of the most serious air pollutants in China for many years. The surface dust can be resuspension under some certain circumstances and the contaminant will make bad influence on our body. It is shown by some researches that Pb can retard children’s intellectual development and weaken their intelligence, while these acknowledged prisoners including Cu, Cd, Cr, Zn, As, and Hg are in a position to change human beings’ nervous and respiratory system. A lot of problems are caused by the pollution of dust, such as laze, respiratory disease, and lung cancer [2]. As a result, research on the surface dust is not only a crucial aspect of environment evaluation, but has also great significance on human health when the urbanization is accelerating.

The evolution process of surface dust includes sediment, diffusion, and resuspension. Nevertheless, it is quite difficult to study the evolution process of surface dust by experiment equipment currently because of the extreme complication of the evolution process. Therefore, simulation has become one of the most important means to research the dust evolution process. The current simulation research includes the diffusion of dust around buildings in cities and the relationship between dust deposition and wind power [3–7]. Dust diffusion is one of the hot topics of the dust evolution. Dust diffusion is a complex dynamic process. It represents a series of particles diffusion process in the wind, such as particles deposition, resuspension, and transit [8]. A number of investigations have been carried out for predicting turbulent diffusion around buildings using wind tunnel tests [9, 10]. Numerical methods for simulating flow and diffusion fields have developed rapidly. Particle deposition in vertical square ventilation duct flows by two different numerical models has been studied [11, 12]. The effects of particle diameter, dimensionless relaxation time, flow direction, and air speed in vertical upward and downward square duct flow on particle deposition velocities are discussed; lift and gravity have also been taken into account. Ali has investigated a time-dependent partial differential equation governing the transport of heavy dust into the atmosphere. Dust concentration is expressed in terms of a series of confluent hyper geometric functions [13]. Sun et al. investigated the microparticle deposition and distribution by employing the Eulerian approach with Reynolds stress turbulent model and a Lagrangian trajectory method [14]. With the rapid development in computer technology, computational fluid dynamics (CFD) method has matured to simulate the ventilation performance and contaminant dispersion and transfer in buildings. Saha et al. have used CFD to assess the effect of wind tunnel sizes on air velocity and concentration boundary layers and on ammonia emission estimation [15]. Tominaga et al. have used CFD to predict the air diffusion around a construction [16, 17]. Besides these, Schneider et al. proposed a semiempirical two-compartment, constant parameter model [8]. Roney and White compared the near-surface wind-tunnel fugitive dust concentration profiles arising from soil surfaces beds with those calculated numerically [18]. However, most of the studies paid attention to simulate the diffusion process of particles. The whole dust evolution process is seldom traced.

A lot of methods can be used in the simulation of dynamic process, such as first principal method (FP), molecular dynamic method (MD), Monte Carlo method (MC), and finite element method (FEM) [19–22]. Among these methods, MC method has been well used in simulation of air pollution dispersion and river water pollution. Three decades ago, Monte Carlo (MC) method was used to study the dynamic process. The basic idea of the MC method is that the solution of the problem is equivalent to hypothetical statistical model parameters using random number, and the parameters are estimated by the statistical model of a sample [23]. Okin proposed MC method to simulate airborne dust diffusion model [24]. Yao et al. used MC method in simulation of emission height effects on building [25]. Tian presented the MC simulation of complex terrain effect on dust diffusion [26]. Dejun GU puts forward a model for the convective boundary layer of the line source diffusion along MC method. Xu et al. simulate gas dispersion based on Monte Carlo which could satisfy the requirement of the long-distance pipeline disasters emergency decision making [27]. The object researched by kinetic Monte Carlo (KMC) method is unbalanced or a relaxation process. Time evolution correct or not is the key factor in the simulation process; simulation time step must be the system real time step, so KMC method is an effective method of studying on the kinetic behavior. Sun et al. applied KMC analysis in environmental risk assessment of a chlorine release accident [28]. Peng and Yuan simulated gas-solid flow behavior in desulphurization tower based on KMC [29]. Because dust have a similar dynamic property with water and air pollutant, KMC method has widely been used in water flow simulation and atmospheric diffusion field. The kinetic Monte Carlo (KMC) related to the MC has the advantage of simulating in a long period. In addition, KMC is a stochastic process, which makes it fit for evolution simulation of the surface dust particles. In the current conditions, it is hard to simulate the KMC evolution of dust particles in a large-scale urban environment by one computer. Therefore, the parallel computing method can be used to shorten the simulation time [30–33].

Virtual reality (VR) has the property of intuition and interactive, which makes up the defeat of computer simulation and makes the results of simulation take on effect of 3D stereo display [34–38]. This paper will adopt visual reality technology (VRT) to visualize the results of dust evolution. A virtual Shanghai University campus is rendered by OpenGL, and a parallel simulation system of dust evolution in the virtual campus is created in order to study the impact of different parameters (such as nonrain period and wind) on dust evolution. By comparison with experimental data, the validity of the model is verified. The research aims at providing theory and quantitative reference for dust particles evolution.

2. Kinetic Monte Carlo Simulation of Surface Dust Evolution

2.1. The 3D Modeling of Virtual Campus

Shanghai University campus is simplified and the 3D model is shown in Figure 1.

Figure 1: 3D model of Shanghai University campus.

The campus is 120 meter in width and 470 meter in length. Since most of surface dust works in the height of 0 to 5 meters, the virtual campus is in the scale of 120 m * 470 m * 6 m. To verify the validity of the KMC simulation method, five positions are chosen as the sampling point and signed as 1, 2, 3, 4, and 5, shown in Figure 1. Location 1 is greatly influenced by the buildings, while the transportation has less impact on it. Location 2 is located on the campus road which is near the greensward. The flow of people has greater impact on it than the buildings. Location 3 and Location 4 are moderately influenced by the buildings and the flow of people. Location 5 is on a road where traffic is quite large and slightly impacted by buildings. By collecting and analyzing the weight of dust at the 5 positions, and comparing them with the KMC simulation result, the initial parameters of dust evolution simulation are defined to suit the real situation.

2.2. KMC Simulation Modeling of Surface Dust Particle in Virtual Campus

The surface dust modeling is dynamic, and the dust particles in the simulation system have three periods, including “emerging,” “moving,” and “vanishing.” With the passage of time, some of the existing particles are vanishing and the new ones are emerging. The survival particles move in the simulation process randomly, and three kinds of events occur in the process of movement: sediment, diffusion, and resuspension. Therefore, the KMC simulation modeling of dust particles consists of five events.

(1)Dust particles merge and join in the simulation system with attributes provided.(2)Dust particles diffuse in the air and update their attributes.(3)Certain dust particles sediment on the ground and update their attributes.(4)Some particles on the ground are resuspension and update their attributes.(5)Dust particles exceeding the life cycle are deleted from the simulation system.

2.2.1. The Initialization of Dust Particles

It is assumed that 10⁸ dust particles are released into the virtual campus environment. Every dust particle has six initial attributes: (1) initial position; (2) initial wind speed; (3) initial size; (4) initial color; (5) shape; (6) life cycle.

The initial positions of dust particles are uniformly distributed. These particles have the same size, shape, and color. Their density is 1800 kg/m³. The initial speed of each particle is the wind speed of its position which can be obtained by numerical simulation of wind field in the virtual campus. Common Computing Fluid Dynamics (CFD) software, such as Fluent, CFX, and Phonics Star-CD, can simulate the wind field.

Assuming air is incompressible viscous fluid, the type of flow is turbulent; density is regarded as constant. Control equation includes continuity and momentum equation. Related equations are as follows [39]:

and present the distance in each coordinate (), and present the velocity component in each coordinate (), is air pressure, is air density, and is time.

To solve the strongly swirling flow problem in the numerical simulation of wind field in the virtual campus, - model is used to produce certain distortion:

denotes inflow point velocity; denotes equivalent length.

Besides, in order to obtain the accurate simulation results of wind field, some improvements are made.

(1) There are two types of grid during Shanghai University virtual campus space grid division, including triangle and hybrid grid. Triangle grid has strong boundary adaptability and hybrid grid can save compute time. In order to simply model, campus local area and plane grid is used [40].

The simulation area contains two parts: building interval area and external wind field area, shown in Figure 2. The grid has been divided into unstructured grid and structured grid separately. In order to improve the accuracy of calculation, the building interval area is encrypted. Especially, the external wind field uses quadrilateral structure grid and the internal wind field uses triangle unstructured grid.

Figure 2: Grid division in virtual campus.

(2) There are three main boundaries in the simulation of wind field: flow inlet boundary, flow outlet boundary, and solid boundary. In order to simplify and optimize the boundary conditions, specific setting is listed in Table 1.

Table 1: Boundary conditions.

Boundary
Condition
Formula

Sky
No slip wall boundary

Inlet condition
Velocity inlet boundary
 m/s

Ground, building
Free-slip boundary
,         

Outlet condition
Pressure-type boundary
,         pa

(3) Before the simulation, it is necessary to make sure of its validity and judge its convergence. The experimental result shows that when iteration is around 1800 times, most kinds of iterative curves are close to our setting number under the residuals: . As shown in Figure 3, the wind field simulation result in this condition has good convergency and fidelity.

Figure 3: Residual curve simulated by Fluent.

According to the meteorological record, the average wind speed 3.2 m/s is set and wind direction is northeast. The wind speed in virtual campus is shown in Figure 4.

Figure 4: Profile of wind speed.

CFD postprocessing software, such as Ensight, Tecplot, and FieldView, can deal with the 3D grid data. Meanwhile, the spatial coordinates of dust particles are written on the text in a grid data form. With the help of CFD, wind speed at each grid node is available and the wind speed anywhere on campus can be drawn by interpolation.

2.2.2. Dust Particles’ KMC Movement

Once the particles emerge, they begin to move when they get the initial attributes. Their move attributes can be derived by the initial ones, considering sediment, diffusion, and resuspension incidents, respectively.

( 1) The Diffusion of Surface Dust Particles. Considering drag force, gravitational setting, Saffman lift force, and turbulent diffusions in the process of computation, motion equations of the particle can be written as [12]

is the drag force per unit particle mass

represents wind speed, is particle speed, is air density, is the density of particles, is particle diameter, is the molecular viscosity of the fluid, is drag coefficient, and is Reynolds number.

This event mainly deals with dust particles on airflow field, namely, particle moves on the MC lattice mentioned above. If particles move with no memory and have equal probability to each direction, the particle’s motion can be taken as Monte-Carlo motion [26].

The trajectory of particles is shown as follows:

is the coordinate direction of , , .

The motion of a particle is defined as

, , and are the mean values in this time step.

, , and are compounded by two parts: relevant part and stochastic part

, , and are the effectiveness of wind to the particles; is the probability to any of the direction, .

So the particles diffuse in the virtual campus environment according to the formulas above.

( 2) The Sediment of Surface Dust Particles. When the particles fulfill the following condition:

they can be taken as sediment on the surface of the earth in this model. At this time, .

( 3) The Resuspension of Surface Dust Particles. Many researches based on Reynolds stress have taken that the particles may be resuspended when the wind speed is increasing to critical friction velocity. The equation can be shown as follows:

is the critical friction velocity, is particle density, is the diameter of the particle, is the acceleration of gravity, is the density of wind, and is a random coefficient with values between 0.16 and 0.21.

When the particle is resuspended,

And the particle velocity is equal to the wind speed.

2.2.3. The Vanish of Surface Dust Particles

The particles have life cycle once they emerge in the virtual environment. They can be deemed as vanished when they move out of the boundary of the virtual campus or they move into the interior of the buildings.

3. KMC-Based Parallel Simulation of Dust Evolution

The premise of KMC-based simulation of dust evolution is that the surface dust location can be described by a point in the virtual campus. Thus, it is hard for a computer to complete the simulation task with the increase of virtual environment or dust particles. At this time, the virtual environment can be divided into several small subspaces and particles in the subspaces are assigned to multiprocessors to simulate concurrently.

3.1. The Way of Data Distribution

Since the evolution of surface dust particles is a random process, the virtual campus space should be divided into continuous space and every space contains the same number of dust particles. Each process simulates the evolution of dust particles in a subspace, which can effectively ensure load balancing of each process [41–43].

When virtual environment is divided into subspaces, the dust particles in regional boundary of a subspace may move to another subspace, which causes the communication between processes. When the particles are in the regional boundary, two adjacent regions need to communicate and confirm where particles are and their specific locations. In order to divide the virtual campus, block data distribution can be achieved by two ways: lattice and sheet divisions, as shown in Figure 5.

(a) Lattice division

(b) Sheet division

Figure 5: KMC division of virtual campus.

The division in Figure 5(a) makes each region be required to communicate with at least three adjacent regions which have the same boundary, while the division in Figure 5(b) makes each region be required to communicate with at most two adjacent regions which have the same boundary. So the division in Figure 5(b) will reduce the traffic and it is used in the division of the virtual campus.

3.2. Communication Optimization Strategy

Particles’ sediment, diffusion, and resuspension should be taken into account when the dust evolution based on KMC is used. To optimize the communication strategy and reduce the traffic between processors, the data storage space is divided into 2 parts: Local Store and Neighbor Copy. Local Store keeps the data of particles in the local subspace, while Neighbor Copy keeps the data of particles in the regional boundary of other neighboring subspace.

When the sediment and resuspension processes are simulated, the attributes of particles in Local Store need renewing and the particles which meet the condition sediment on the ground or re-suspend in the air. So the processors do not need to communicate with each other and reduce the communication frequency.

In the diffusion process, particles in the regional boundary update their attributes both in Local Store and Neighbor Copy. However, particles which are not in the regional boundary only update their attributes in Local Store, which will reduce the communication frequency and traffic between processors.

According to the description above, each process not only stores the data of particles in the local subspace, but also stores the data of particles in the neighborhood space. In Figure 6, the virtual environment is divided into three processors and dashed areas contain the Neighbor Copy space of process 2 because it needs to communicate with processes 1 and 3 in the diffusion process.

Figure 6: Dust data storage model in three processes.

Data structure of particle in the communication is described as follows:

where pos and vel are the particle’s location and velocity in the virtual campus at this time; size, color, shape, and lifecycle of each particle in the Neighbor Copy have the same value.

3.3. KMC-Based Dust Evolution Parallel Simulation Algorithm

Algorithm 1 is the frame of KMC-based dust evolution parallel simulation algorithm. Three main processes of dust evolution are shown in Algorithms 2, 3, and 4.

(1) Begin

(2) Initialization MPI, Define the number of processes and simulation time ;

(3) Master process reads the data structure of dust particles;

(4) Master process distributes the data to Local_Store of each sub-processors according to

the data distribution way of sheet division, and the simulation timer ;

(5) if , go to (6); otherwise, go to (10);

(6) Each sub-processor executes sediment of dust particles;

(7) Each sub-processor executes re-suspension of dust particles;

(8) Each sub-processor executes diffusion of dust particles;

(9) Each sub-processor updates the value of simulation timer , go to (5);

(10) Each sub-processor sends data structure of dust particles to the master process;

(11) The master processor collects the data sent by each sub-processor and writes to the file Output;

(12) End

Algorithm 1: Frame of KMC-based dust evolution parallel simulation algorithm.

(1) Begin

(2) Each processor updates the position of each particle in the Local_Store;

(3) If , go to (4); Otherwise, go to (5);

(4) Particle is sediment on the ground, ;

(5) If all particles in Local_Store are finished scanning, go to (6); Otherwise, go to (3);

(6) End

Algorithm 2: The sediment simulation algorithm of dust particles.

(1) Begin

(2) Each processor updates the position of each particle in the Local_Store;

(3) If , go to (5); Otherwise go to (5);

(4) If the wind speed in the location of particle is greater than the critical friction velocity,

Particle is re-suspended and ;

(5) If all particles in Local_Store are finished scanning, go to (6); Otherwise, go to (3);

(6) End

Algorithm 3: The re-suspension simulation algorithm of dust particles.

(1) Begin

(2) Each processor sends the update data of particles in the border to other processors, then

receives the update data of particles from other processors and copies them in Neighbor_Copy;

(3) Calculate the diffusion probability and diffusion direction of the particle j in Local_Store;

(4) Execute the diffusion of particles, update , mark particles which enter other Neighbor_Copy;

(5) If all particles in Local_Store are finished scanning, go to (6); Otherwise, go to (3);

(6) Communicate with other processor to update the data of particles in Neighbor_Copy;

(7) End

Algorithm 4: The diffusion simulation algorithm of dust particles.

4. Visualization on Surface Dust Evolution

The coordinates of particles obtained from the concurrent calculation are recorded in the text, and OpenGL makes the dust evolution visible. To see the dust evolution clearly, the boundaries of virtual campus are drawn as transparent.

In the visualization on surface dust evolution, the dust particles are drawn pro rata because of the large amount, and the coefficient scale is 10⁻³:

Suppose that the total particle amount in output at moment is . Then the visible amount of dust particles is at the simulation of moment.

The visualization algorithm of dust particles’ KMC evolution is shown in Algorithm 5. And Figure 7 shows the visualization result. It is clear to see the evolution of dust particles in the virtual campus.

(1) Begin

(2) Create a transparent border of virtual campus;

(3) Render a yellow ground;

(4) Render the buildings in the virtual campus;

(5) Set the attributes of dust particles, including their size, color and shape;

(6) Read the coordinates of dust particles from file Output, calculate NParts;

(7) If the particles are in the area of sampling points ,

calculate the accumulation amount of dust particles ;

(8) Render the particles according to the value of NParts;

(9) Output the accumulation amount of dust particles in five sampling points;

(10) End

Algorithm 5: The visualization algorithm of dust particles’ KMC evolution.

Figure 7: The visualization result of dust evolution in virtual campus.

5. Results and Analysis

In order to evaluate the effectiveness of KMC-based parallel simulation algorithm of dust evolution in virtual campus environment, the experiment is designed as follows.

Dust in five collection areas of campus is collected each nonrain day during four months. The weight of dust is gained by a delicate electronic balance and recorded. At the same time, the weather condition like wind scale, rainy day, and nonrain period is marked. According to the records, the northeast wind is the most frequent wind during the experimental period. So the following analysis is based on the condition of northeast wind. The record of dust accumulation in the northeast wind is shown in Table 2.

Table 2: The record of dust accumulation in different locations.

Nonrain period (day)
Location 1
Location 2
Location 3
Location 4
Location 5

1
0.406
1.095
2.095
0.745
0.171

1
0.037
1.378
0.782
0.215
0.663

1
0.098
1.790
3.964
1.846
0.285

1
0.347
2.524
1.208
2.558
0.705

1
0.166
1.251
0.920
1.974
0.520

1
0.234
2.372
1.272
2.705
0.984

1
0.336
1.128
2.449
1.612
2.056

1
0.282
0.877
0.920
0.614
0.551

1
0.112
0.654
0.346
1.115
1.554

1
0.193
0.511
0.628
1.016
2.063

1
0.205
0.960
0.922
0.838
0.924

2
0.056
1.226
2.198
0.985
0.221

2
0.296
1.865
1.088
1.084
0.547

2
0.294
1.384
0.887
2.015
0.552

2
0.259
0.597
0.611
1.224
1.866

2
0.180
0.599
0.528
1.054
1.586

2
0.181
0.497
0.314
0.215
1.470

3
0.264
1.426
1.149
0.446
0.623

3
0.226
1.495
0.808
0.791
0.753

3
0.184
1.064
0.839
0.439
3.740

4
0.045
1.804
2.656
1.575
0.219

4
0.055
1.452
0.484
0.895
1.098

4
0.214
3.577
1.203
1.388
0.315

4
0.214
0.833
0.804
1.034
1.401

4
0.195
0.804
0.510
0.964
0.854

4
0.073
1.506
1.565
0.921
1.286

4
0.156
0.691
1.476
1.366
1.481

5
0.082
1.533
1.194
0.271
0.148

5
0.323
0.602
1.958
0.978
0.694

5
0.194
0.399
1.430
1.138
2.630

6
0.049
2.416
2.467
1.561
0.453

6
0.116
0.385
0.348
0.859
1.806

7
0.080
0.804
1.070
0.999
0.575

7
0.150
0.482
0.602
0.471
1.404

Figure 8 compares the experimental results and simulation results based on KMC serial and parallel algorithm of dust evolution by the effect of different nonrain periods in five collection areas.

(a) Relationship of nonrain period and dust fall at Location 1

(b) Relationship of nonrain period and dust fall at Location 2

(c) Relationship of nonrain period and dust fall at Location 3

(d) Relationship of nonrain period and dust fall at Location 4

(e) Relationship of nonrain period and dust fall at Location 5

Figure 8: Comparison of experimental and simulation results of dust fall accumulation.

Figures 8(a) to 8(e) show that the dust fall accumulation became heavier and heavier with the increase of the nonrain period, which proves the effectiveness of simulation algorithm. From Figure 8, serial and parallel simulations had the same accumulation amount of dust particles in five points, which shows the accuracy of parallel simulation algorithm.

To evaluate the validity of parallel simulation algorithm, its acceleration ratio and efficiency are calculated and results are shown in Table 3.

Table 3: Comparison study of results of parallel calculation in different processors.

The number of processes

/hour

1
1
26.31

2
1.49
17.66

4
1.91
13.77

8
2.87
9.17

16
6.08
4.33

The parallel acceleration ratio is defined as

where is the time used by serial algorithm and is the time used by parallel algorithm in processes.

From Table 2, the value of the acceleration ratio is small because the algorithm is related to the text operation, while the acceleration ratio increases with the number of the processes and the computation time reduces evidently. This indicates that the parallel algorithm on dust evolution can promote the efficiency although KMC evolution algorithm needs lots of boundary exams on the particles, which makes the simulation of large-scale virtual environment possible.

6. Conclusion

It is efficient to use the parallel algorithm to simulate the KMC evolution of surface dust particles in large-scale virtual environment. A parallel simulation algorithm of particles’ KMC evolution is proposed. It is useful to balance the load of every process and reduce the communication expense among processes with the help of data distribution way of sheet division and communication optimizing strategy. The experiment results show that simulation operation time is shortened enormously, the acceleration ratio is easy to get, and the parallel efficiency is promoted due to the reasonable process numbers in the parallel simulation algorithm which also compensates the disability of single computer. With the 3D visible simulation result, researchers can have a good understanding of the segmentation, diffusion, and resuspension of dust particles and analyze their movement discipline to lay a theoretical foundation for the dust prevention.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported in part by NSFC (Project no. 41101454), the Grand Science & Technology Program Shanghai China (no. 13111101300), and Industrial Innovation Grand Projects (no. 07CH-008).

References

J. Chang, M. Liu, L.-J. Hou, S.-Y. Xu, X. Lin, and S. Ballo, “Concept, pollution character and environmental effect of urban surface dust,” Chinese Journal of Applied Ecology, vol. 18, no. 5, pp. 1153–1158, 2007.

N. Li, Effect of Haze on Respiratory Healthy in Guangzhou, Wuhan University Of Technology, 2009.

L.-L. Da, T.-W. Yang, Y.-Y. Li, and X.-T. Lu, “Accelerating volume rendering of 3D datasets based on PC hardware,” Journal of System Simulation, vol. 17, no. 10, pp. 2422–2425, 2005.

D.-P. Xi and W.-P. Jiang, “Research and application of three-dimensional visibility based on digital map,” Earth Science, vol. 27, no. 3, pp. 278–284, 2002.

C. Ye, J.-S. Wang, and X. I. Li, “Field measurement and numerical simulation for pollutant dispersion from vehicular exhaust in street canyon,” Environmental Chemistry, vol. 25, pp. 364–366, 2006.

A. Zhang, L. Zhang, and J. Zhou, “Numerical simulation of wind environment around two adjacent buildings,” Chinese Journal of Computational Mechanics, vol. 20, no. 5, pp. 553–558, 2003.

T. L. Chan, G. Dong, C. W. Leung, C. S. Cheung, and W. T. Hung, “Validation of a two-dimensional pollutant dispersion model in an isolated street canyon,” Atmospheric Environment, vol. 36, no. 5, pp. 861–872, 2002.

T. Schneider, J. Kildes, and N. O. Breum, “A two compartment model for determining the contribution of sources, surface deposition and resuspension to air and surface dust concentration levels in occupied rooms,” Building and Environment, vol. 34, no. 5, pp. 583–595, 1999.

R. Yao, Q. Qiao, and X. Yu, “Wind tunnel simulation of flow and dispersion around complex buildings,” Radiation Protection Bulletin, vol. 22, pp. 1–6, 2002.

D.-P. Guo, Q.-D. Qiao, and R.-T. Yao, “Examining the k-ε(RNG) model and LES of flow feature and turbulence dispersion around a building by means of wind tunnel tests,” Journal of Experiments in Fluid Mechanics, vol. 25, no. 5, pp. 55–63, 2011.

J. Zhang and A. Li, “Study on particle deposition in vertical square ventilation duct flows by different models,” Energy Conversion and Management, vol. 49, no. 5, pp. 1008–1018, 2008.

R. Gao and A. Li, “Modeling deposition of particles in vertical square ventilation duct flows,” Building and Environment, vol. 46, no. 1, pp. 245–252, 2011.

I. Ali, S. L. Kalla, and H. G. Khajah, “A time dependent model for the transport of heavy pollutants from ground-level aerial sources,” Applied Mathematics and Computation, vol. 105, no. 1, pp. 91–99, 1999.

K. Sun, L. Lu, and H. Jiang, “A numerical study of bend-induced particle deposition in and behind duct bends,” Building and Environment, vol. 52, pp. 77–87, 2012.

C. K. Saha, W. Wu, G. Zhang, and B. Bjerg, “Assessing effect of wind tunnel sizes on air velocity and concentration boundary layers and on ammonia emission estimation using computational fluid dynamics (CFD),” Computers and Electronics in Agriculture, vol. 78, no. 1, pp. 49–60, 2011.

Y. Tominaga, S. Murakami, and A. Mochida, “CFD prediction of gaseous diffusion around a cubic model using a dynamic mixed SGS model based on composite grid technique,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 67-68, pp. 827–841, 1997.

W. Nie, W.-M. Cheng, G. Zhou, and Y. Yao, “The numerical simulation on the regularity of dust dispersion in whole-rock mechanized excavation face with different air draft amount,” Procedia Engineering, vol. 26, pp. 961–971, 2011.

J. A. Roney and B. R. White, “Comparison of a two-dimensional numerical dust transport model with experimental dust emissions from soil surfaces in a wind tunnel,” Atmospheric Environment, vol. 44, no. 4, pp. 512–522, 2010.

P. Song, J. Lu, Q. Hu, M. Zhao, and B. Yang, “Application and development of computer simulation in thin film deposition,” Materials Review, vol. 17, pp. 154–157, 2003.

S. Razmyan and F. Hosseinzadeh Lotfi, “An application of Monte-Carlo-based sensitivity analysis on the overlap in discriminant analysis,” Journal of Applied Mathematics, vol. 2012, Article ID 315868, 14 pages, 2012.

Y. X. Jie, H. N. Yuan, and H. D. Zhou, “Bending moment calculations for piles based on the finite element method,” Journal of Applied Mathematics, vol. 2013, Article ID 784583, 19 pages, 2013.

L. Zhang and Z. Chen, “A stabilized mixed finite element method for single-phase compressible flow,” Journal of Applied Mathematics, vol. 2011, Article ID 129724, 16 pages, 2011.

W. Zhu, G. Hu, X. Hu, L. Hongbo, and W. Zhang, “Visual simulation of GaInP thin film growth,” Simulation Modelling Practice and Theory, vol. 18, no. 1, pp. 87–99, 2010.

G. Okin, “The Role of Spatial Variability in Wind Erosion and Dust Emission,” Geophysical Research Abstracts 12583, 2003.

R. Yao, H. Hao, and E. Hu, “Comparison of two kinds of atmospheric dispersion model chains in rodos,” Radiation Protection, vol. 23, pp. 146–155, 2003.

R. Tian, “Monte-Carlo model simulates the influence of complex terrain on diffusion,” Scientia Atmospherica Sinica, vol. 18, pp. 37–42, 1994.

K. Xu, H. G. He, and Y. C. Zhu, “Study on dispersion simulation of long-distant pipeline leaked gas based on Monte-Carlo,” Journal of Safety Science and Technology, vol. 8, pp. 18–23, 2012.

Y. Sun, Y. Qian, and Y. Zhang, “Application of Monte Carlo analysis in environmental risk assessment of a chlorine release accident,” Acta Scientiae Circumstantiae, vol. 31, no. 11, pp. 2570–2577, 2011.

Z.-B. Peng and Z.-L. Yuan, “Numerical simulation of gas-solid flow behaviours in desulfurization tower based on Monte Carlo,” Proceedings of the Chinese Society of Electrical Engineering, vol. 28, no. 14, pp. 6–14, 2008.

S. Tanaka, T. Nishide, and K. Sakurai, “Efficient implementation for QUAD stream cipher with GPUs,” Computer Science and Information Systems, vol. 10, no. 2, pp. 897–911, 2013.

Y. Shang, G. Lu, and L. Shang, “Parallel processing on block-based Gauss-Jordan algorithm for desktop grid,” Computer Science and Information Systems, vol. 8, no. 3, pp. 739–759, 2011.

F. H. Pereira and S. I. Nabeta, “A parallel wavelet-based algebraic multigrid black-box solver and preconditioner,” Journal of Applied Mathematics, vol. 2012, Article ID 894074, 15 pages, 2012.

C. Han, T. Feng, G. He, and T. Guo, “Parallel variable distribution algorithm for constrained optimization with nonmonotone technique,” Journal of Applied Mathematics, vol. 2013, Article ID 295147, 7 pages, 2013.

Z. Wenhua, F. Xiong, H. Guihua, S. Yupeng, and X. Hu, Virtual Reality Technology and Application, Intellectual Property Press, Beijing, China, 2007.

Y.-M. Chen, J.-S. Bao, Y. Jin, C.-C. Xu, and Y.-C. Yang, “Key technology study and application of engineering analysis data's immersive visualization,” Journal of System Simulation, vol. 16, no. 10, pp. 2309–2312, 2004.

A. Attenberger and K. Buchenrieder, “Modeling and visualization of classification-based control schemes for upper limb prostheses,” Computer Science and Information Systems, vol. 10, no. 1, pp. 349–367, 2013.

B. Jin-Song, J. Ye, M. Deng-Zhe, and Y. Jun-Qi, “Immersive scientific visualization with realist geometry,” Journal of System Simulation, vol. 15, pp. 653–655, 2003.

S. Sang, J. Zhao, H. Wu, S. Chen, and Q. An, “Modeling and simulation of a spherical mobile robot,” Computer Science and Information Systems, vol. 7, no. 1, pp. 51–62, 2010.

P. F. Li, M. Y. Xu, and F. F. Wang, FLUENT GAMBIT ICEM CFD Tecplot, Beijing Institute of Technology Press, Beijing, China, 2005.

J. F. Xie, “A comparative study of various numerical simulation approaches to wind environment within urban green space,” Agriculture Network Information, vol. 26, no. 07, pp. 18–21, 2011.

K. Lu and X. Lin, “Implementing load balance in MPI parallel program,” Microcomputer Information, vol. 05X, pp. 226–227, 2007.

L. Lu, “Research on parallel program design strategy,” Electronic Computers, vol. 141, no. 6, pp. 2–8, 1999.

B. Zhou, J. Shen, and Q. Peng, “Communication scheme of parallel clustering algorithm for PCs cluster,” Computer Engineering, vol. 30, no. 7, pp. 20–21, 2004.