# Numerical Simulation of an Air Pollution Model on Industrial Areas by Considering the Influence of Multiple Point Sources

Authors: Pravitra Oyjinda; Nopparat Pochai
Journal: International Journal of Differential Equations
 (2019)

---

## Abstract

A numerical simulation on a two-dimensional atmospheric diffusion equation of an air pollution measurement model is proposed. The considered area is separated into two parts that are an industrial zone and an urban zone. In this research, the air pollution measurement by releasing the pollutant from multiple point sources above an industrial zone to the other area is simulated. The governing partial differential equation of air pollutant concentration is approximated by using a finite difference technique. The approximate solutions of the air pollutant concentration on both areas are compared. The air pollutant concentration levels influenced by multiple point sources are also analyzed.


---

## Body

## 1. Introduction

A rapid growth of industrial sector can explain that air pollution affects the health of human being who lives around industrial areas. The air pollution has become a major problem of human life and environment. The purpose of this research was to study the air pollution assessment problem in two adjacent zones: industrial and urban zones by using the atmospheric diffusion model. In [1], the simulation of two-dimensional advection-diffusion model with a point source was presented. The numerical solutions were solved by using the finite difference techniques. In [2], the researchers used the mathematical model to simulate the dispersion of sulfur dioxide concentration with the wind and diffusion parameters regarding the reference atmospheric stability. In [3], the mass transport model consisted of the stream function, vorticity, and convection-diffusion equation. The smoke dispersion which released into the atmosphere from one and two-point sources was considered with obstacle domain. The approximated solutions were solved by using the finite element techniques. In [4], the researchers studied the smoke dispersion model in a two-dimensional space by considering two and three point sources with a two obstacles in the domain. In [5], the two-dimensional advection-diffusion equation with mesoscale wind, eddy diffusivity profiles, and removal mechanisms was introduced. Then, the primary pollutant released into the atmosphere from an area source, which was also considered. In [6], the researchers studied the two-dimensional advection-diffusion equation of primary and secondary pollutants. The area source with removal mechanisms and the point source considering on the boundary were proposed. The solutions of air pollution in [5, 6] were estimated by using the Crank-Nicolson implicit methods. In [7], the air-quality model in the three-dimensional with variations of the atmospheric stability classes and wind velocities from multiple sources was analyzed. The fractional step methods were used in order to predict the air pollutant concentration in [2, 4, 7]. In [8], the atmospheric diffusion model was used to describe the dispersion of air pollutant concentration near an industrial zone. The problem was considered by controlling the air pollution emission under a point source. From the numerical experiments, it was indicated that the air pollution control was necessary for air-quality management. In [9], the researchers studied the dispersion behavior of air pollution in the tunnel under a Bangkok sky train platform by using the simulation of a three-dimensional air-quality model. This model was considered varied cases on the wind inflow with obstacles. In [10], the three-dimensional advection-diffusion equation was considered to approximate the concentration of air pollutant in a heavy traffic area under the Bangkok sky train station. The numerical simulations were studied for three cases that were the average of source or sink emissions, the moving of source or sink emissions, and the mix of source and sink emissions. The explicit finite difference scheme was used to solve the air pollutant concentration in [8–10].The source that is smokestack of industrial factory or power plant discharges the air pollution into the system. The genesis of air pollution is the cause of problems. In this research, the simple finite difference methods are used for solving the atmospheric diffusion equation.


## 2. Governing Equation

### 2.1. The Atmospheric Diffusion Equation

The diffusion model generally uses Gaussian plume idea, which is the well-known atmospheric diffusion equation. It represents the behavior of air pollution in industrial areas. The dispersion of pollutant concentration from multiple point sources is described by the following three-dimensional advection-diffusion equation:(1)
∂
c
∂
t
+
u
∂
c
∂
x
+
v
∂
c
∂
y
+
w
∂
c
∂
z
=
k
x
∂
2
c
∂
x
2
+
k
y
∂
2
c
∂
y
2
+
k
z
∂
2
c
∂
z
2
+
s
,where c
=
c
(
x
,
y
,
z
,
t
) is the concentration of air pollutant at (
x
,
y
,
z
) and time t (
k
g
/
m
3
), u
,
v, and w are the wind velocity component (
m
/
s
e
c
) in x
-
,
y
-
,
z-directions, respectively, k
x
,
k
y, and k
z are the diffusion coefficient (
m
2
/
s
e
c
) in x
-
,
y
-
,
z-directions respectively, and s is the sink rate of air pollutants (
s
e
c
-
1
).The assumptions of (1) are defined that the concentrations of air pollutant are emitted from continued point sources. The advection and diffusion in y-direction are laterally averaged. By the assumption, we can also eliminate all terms in y-direction. Therefore, the governing equation can be written as(2)
∂
c
∂
t
+
u
∂
c
∂
x
+
w
∂
c
∂
z
=
k
x
∂
2
c
∂
x
2
+
k
z
∂
2
c
∂
z
2
+
s
.The initial condition is assumed under the cold start assumption. That is,(3)
c
x
,
z
,
0
=
0
,for all x
>
0 and z
>
0. The boundary conditions are assumed that(4)
c
0
,
z
,
t
=
0
,
(5)
∂
c
∂
x
L
,
z
,
t
=
0
,
(6)
∂
c
∂
z
x
,
0
,
t
=
0
,
(7)
∂
c
∂
z
x
,
H
,
t
=
0
,for all t
>
0 where L is the length of the domain in x-direction and H is the height of the inversion layer. The concentration at the point sources is assumed to be the constant variables as(8)
c
x
p
,
0
,
t
=
c
s
p
,for p
=
1,2 where x
p is the position of the point source p in the x-direction and c
s
p is the concentration value at the point source of p.


### 2.2. The Nondimensional Form Equation

Now, we introduce the dimensionless form of equation (2). The nondimensional variables are denoted by letting C
=
c
/
c
m
a
x, X
=
x
/
l
x, Z
=
z
/
l
z, T
=
t
/
t
m
a
x, D
x
=
k
x
/
l
x
u
m
a
x, D
z
=
k
z
/
l
z
u
m
a
x, U
=
u
/
u
m
a
x, and W
=
β
w
m
a
x
/
u
m
a
x when β
=
w
/
w
m
a
x. We define c
m
a
x
=
max
⁡
c
x
,
z
,
t
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x, u
m
a
x
=
max
⁡
{
u
(
x
,
z
,
t
)
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x
}, w
m
a
x
=
m
a
x
{
w
(
x
,
z
,
t
)
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x
}, and t
m
a
x is a stationary time. Thus the nondimensional equation of air pollution is as follows:(9)
1
S
T
∂
C
∂
T
+
U
∂
C
∂
X
+
W
∂
C
∂
Z
=
D
x
∂
2
C
∂
X
2
+
D
z
∂
2
C
∂
Z
2
+
S
,where l
=
max
⁡
l
x
,
l
z and S
T
=
t
m
a
x
u
m
a
x
/
l when S
<
0 that means the air pollutant concentrations are absorbed from the atmosphere by the chemical reaction.


## 2.1. The Atmospheric Diffusion Equation

The diffusion model generally uses Gaussian plume idea, which is the well-known atmospheric diffusion equation. It represents the behavior of air pollution in industrial areas. The dispersion of pollutant concentration from multiple point sources is described by the following three-dimensional advection-diffusion equation:(1)
∂
c
∂
t
+
u
∂
c
∂
x
+
v
∂
c
∂
y
+
w
∂
c
∂
z
=
k
x
∂
2
c
∂
x
2
+
k
y
∂
2
c
∂
y
2
+
k
z
∂
2
c
∂
z
2
+
s
,where c
=
c
(
x
,
y
,
z
,
t
) is the concentration of air pollutant at (
x
,
y
,
z
) and time t (
k
g
/
m
3
), u
,
v, and w are the wind velocity component (
m
/
s
e
c
) in x
-
,
y
-
,
z-directions, respectively, k
x
,
k
y, and k
z are the diffusion coefficient (
m
2
/
s
e
c
) in x
-
,
y
-
,
z-directions respectively, and s is the sink rate of air pollutants (
s
e
c
-
1
).The assumptions of (1) are defined that the concentrations of air pollutant are emitted from continued point sources. The advection and diffusion in y-direction are laterally averaged. By the assumption, we can also eliminate all terms in y-direction. Therefore, the governing equation can be written as(2)
∂
c
∂
t
+
u
∂
c
∂
x
+
w
∂
c
∂
z
=
k
x
∂
2
c
∂
x
2
+
k
z
∂
2
c
∂
z
2
+
s
.The initial condition is assumed under the cold start assumption. That is,(3)
c
x
,
z
,
0
=
0
,for all x
>
0 and z
>
0. The boundary conditions are assumed that(4)
c
0
,
z
,
t
=
0
,
(5)
∂
c
∂
x
L
,
z
,
t
=
0
,
(6)
∂
c
∂
z
x
,
0
,
t
=
0
,
(7)
∂
c
∂
z
x
,
H
,
t
=
0
,for all t
>
0 where L is the length of the domain in x-direction and H is the height of the inversion layer. The concentration at the point sources is assumed to be the constant variables as(8)
c
x
p
,
0
,
t
=
c
s
p
,for p
=
1,2 where x
p is the position of the point source p in the x-direction and c
s
p is the concentration value at the point source of p.


## 2.2. The Nondimensional Form Equation

Now, we introduce the dimensionless form of equation (2). The nondimensional variables are denoted by letting C
=
c
/
c
m
a
x, X
=
x
/
l
x, Z
=
z
/
l
z, T
=
t
/
t
m
a
x, D
x
=
k
x
/
l
x
u
m
a
x, D
z
=
k
z
/
l
z
u
m
a
x, U
=
u
/
u
m
a
x, and W
=
β
w
m
a
x
/
u
m
a
x when β
=
w
/
w
m
a
x. We define c
m
a
x
=
max
⁡
c
x
,
z
,
t
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x, u
m
a
x
=
max
⁡
{
u
(
x
,
z
,
t
)
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x
}, w
m
a
x
=
m
a
x
{
w
(
x
,
z
,
t
)
:
0
⩽
x
⩽
L
,
0
⩽
z
⩽
H
,
0
⩽
t
⩽
t
m
a
x
}, and t
m
a
x is a stationary time. Thus the nondimensional equation of air pollution is as follows:(9)
1
S
T
∂
C
∂
T
+
U
∂
C
∂
X
+
W
∂
C
∂
Z
=
D
x
∂
2
C
∂
X
2
+
D
z
∂
2
C
∂
Z
2
+
S
,where l
=
max
⁡
l
x
,
l
z and S
T
=
t
m
a
x
u
m
a
x
/
l when S
<
0 that means the air pollutant concentrations are absorbed from the atmosphere by the chemical reaction.


## 3. Numerical Methods

We use the finite difference methods for calculating the nondimensional form of the atmospheric diffusion equation. In (9), we get the concentration of C at each time T
n
+
1 from T
n when Δ
T is a time increment. The solution of concentration at (
X
,
Z
,
T
) is denoted by C
(
X
i
,
Z
j
,
T
n
)
=
C
i
,
j
n. The domain is divided by the grid spacing in X-direction and Z-direction which are Δ
X and Δ
Z, respectively, where X
i
=
i
Δ
X and Z
j
=
j
Δ
Z. The approximate solutions are obtained by using the following methods.


### 3.1. Forward Time Central Space Scheme

The first method, we use the forward difference in transient term that is(10)
∂
C
∂
T
=
C
i
,
j
n
+
1
-
C
i
,
j
n
Δ
T
.Then, the centered difference for the advection and diffusion in X-direction and Z-direction is applied as follows:(11)
∂
C
∂
X
=
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
,
(12)
∂
C
∂
Z
=
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
,
(13)
∂
2
C
∂
X
2
=
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
,
(14)
∂
2
C
∂
Z
2
=
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
,respectively. We substitute (10)-(14) into (9). It will be(15)
1
S
T
C
i
,
j
n
+
1
-
C
i
,
j
n
Δ
T
+
U
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
+
W
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
=
D
x
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
+
D
z
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
+
S
.Thus, the forward time central space (FTCS) scheme of the nondimensional mathematical model is(16)
C
i
,
j
n
+
1
=
d
x
-
A
x
C
i
+
1
,
j
n
+
d
x
+
A
x
C
i
-
1
,
j
n
+
1
-
2
d
x
-
2
d
z
C
i
,
j
n
+
d
z
+
A
z
C
i
,
j
-
1
n
+
d
z
-
A
z
C
i
,
j
+
1
n
+
S
T
Δ
T
S
,where A
x
=
S
T
Δ
T
U
/
2
Δ
X, A
z
=
S
T
Δ
T
W
/
2
Δ
Z, d
x
=
S
T
(
Δ
T
)
D
x
/
Δ
X
2
,
d
z
=
S
T
(
Δ
T
)
D
z
/
Δ
Z
2. The stability of the forward time central space scheme can be investigated by using von Neumann stability analysis. We can obtain that the stability condition is 0
⩽
2
d
x
+
A
x
+
A
z
⩽
1.


### 3.2. Backward Time Central Space Scheme

The second method, we use the backward difference in transient term that is(17)
∂
C
∂
T
=
C
i
,
j
n
-
C
i
,
j
n
-
1
Δ
T
.Then, the centered difference for the advection and diffusion in X-direction and Z-direction is utilized as follows:(18)
∂
C
∂
X
=
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
,
(19)
∂
C
∂
Z
=
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
,
(20)
∂
2
C
∂
X
2
=
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
,
(21)
∂
2
C
∂
Z
2
=
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
,respectively. We substitute (17)-(21) into (9). It obtains that(22)
1
S
T
C
i
,
j
n
-
C
i
,
j
n
-
1
Δ
T
+
U
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
+
W
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
=
D
x
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
+
D
z
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
+
S
.Therefore, the backward time central space (BTCS) scheme of this research becomes(23)
A
x
-
d
x
C
i
+
1
,
j
n
+
1
+
A
x
+
d
x
C
i
-
1
,
j
n
+
1
+
1
+
2
d
x
+
2
d
z
C
i
,
j
n
+
1
-
A
z
+
d
z
C
i
,
j
-
1
n
+
1
+
A
z
-
d
z
C
i
,
j
+
1
n
+
1
=
C
i
,
j
n
+
S
T
Δ
T
S
.The stability of the implicit backward time central space scheme can be investigated by using von Neumann stability analysis. We can obtain that the method is an unconditionally stable method.


## 3.1. Forward Time Central Space Scheme

The first method, we use the forward difference in transient term that is(10)
∂
C
∂
T
=
C
i
,
j
n
+
1
-
C
i
,
j
n
Δ
T
.Then, the centered difference for the advection and diffusion in X-direction and Z-direction is applied as follows:(11)
∂
C
∂
X
=
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
,
(12)
∂
C
∂
Z
=
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
,
(13)
∂
2
C
∂
X
2
=
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
,
(14)
∂
2
C
∂
Z
2
=
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
,respectively. We substitute (10)-(14) into (9). It will be(15)
1
S
T
C
i
,
j
n
+
1
-
C
i
,
j
n
Δ
T
+
U
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
+
W
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
=
D
x
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
+
D
z
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
+
S
.Thus, the forward time central space (FTCS) scheme of the nondimensional mathematical model is(16)
C
i
,
j
n
+
1
=
d
x
-
A
x
C
i
+
1
,
j
n
+
d
x
+
A
x
C
i
-
1
,
j
n
+
1
-
2
d
x
-
2
d
z
C
i
,
j
n
+
d
z
+
A
z
C
i
,
j
-
1
n
+
d
z
-
A
z
C
i
,
j
+
1
n
+
S
T
Δ
T
S
,where A
x
=
S
T
Δ
T
U
/
2
Δ
X, A
z
=
S
T
Δ
T
W
/
2
Δ
Z, d
x
=
S
T
(
Δ
T
)
D
x
/
Δ
X
2
,
d
z
=
S
T
(
Δ
T
)
D
z
/
Δ
Z
2. The stability of the forward time central space scheme can be investigated by using von Neumann stability analysis. We can obtain that the stability condition is 0
⩽
2
d
x
+
A
x
+
A
z
⩽
1.


## 3.2. Backward Time Central Space Scheme

The second method, we use the backward difference in transient term that is(17)
∂
C
∂
T
=
C
i
,
j
n
-
C
i
,
j
n
-
1
Δ
T
.Then, the centered difference for the advection and diffusion in X-direction and Z-direction is utilized as follows:(18)
∂
C
∂
X
=
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
,
(19)
∂
C
∂
Z
=
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
,
(20)
∂
2
C
∂
X
2
=
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
,
(21)
∂
2
C
∂
Z
2
=
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
,respectively. We substitute (17)-(21) into (9). It obtains that(22)
1
S
T
C
i
,
j
n
-
C
i
,
j
n
-
1
Δ
T
+
U
C
i
+
1
,
j
n
-
C
i
-
1
,
j
n
2
Δ
X
+
W
C
i
,
j
+
1
n
-
C
i
,
j
-
1
n
2
Δ
Z
=
D
x
C
i
+
1
,
j
n
-
2
C
i
,
j
n
+
C
i
-
1
,
j
n
Δ
X
2
+
D
z
C
i
,
j
+
1
n
-
2
C
i
,
j
n
+
C
i
,
j
-
1
n
Δ
Z
2
+
S
.Therefore, the backward time central space (BTCS) scheme of this research becomes(23)
A
x
-
d
x
C
i
+
1
,
j
n
+
1
+
A
x
+
d
x
C
i
-
1
,
j
n
+
1
+
1
+
2
d
x
+
2
d
z
C
i
,
j
n
+
1
-
A
z
+
d
z
C
i
,
j
-
1
n
+
1
+
A
z
-
d
z
C
i
,
j
+
1
n
+
1
=
C
i
,
j
n
+
S
T
Δ
T
S
.The stability of the implicit backward time central space scheme can be investigated by using von Neumann stability analysis. We can obtain that the method is an unconditionally stable method.


## 4. Numerical Experiment

The two-dimensional atmospheric diffusion equation (9) with a dimension 1,000
×
1,000 m
2 will be considered. The uniform wind velocities and constant diffusion coefficients are introduced. We choose that the wind velocities in x-direction and z-direction are 0.1 and 0.05
m
/
s
e
c, respectively. The diffusion coefficients in x-direction and z-direction are 4.5
×
1
0
-
1 and 4.5
×
1
0
-
5
m
2
/
s
e
c, respectively. The grid spacing is Δ
x
=
Δ
z
=
25 m. and the time interval is 20 s
e
c. In this research, we present two cases. The first case considers a point source when the concentration is 0.5
k
g
/
m
3. The second case considers two-point source when the concentration is 0.25 and 0.25
k
g
/
m
3. The air pollutants in (8) are released into our system. These examples are solved by using the forward time central space and the backward time central space schemes in (16) and (23), respectively, with the initial and boundary conditions (3) to (7).In Figure1, model of the problem is shown. The physical problem composed of two zones: an industrial zone and an urban zone with the stable wind along the x-axis and z-axis. The point sources are laid along the x-axis. We assume that the primary air pollutants are released from a factory smokestack by a single point source and coupled point sources on industrial zone. The emissions of air pollution are influenced on the urban zone by the rate of air pollutant absorption. In the numerical experiment, the considered domain of solutions is shown in Figure 2.Figure 1
Model of the problem.Figure 2
Domain of solutions.


## 5. Discussion

The air pollutant emission from multiple point sources above an industrial zone to the urban area is presented. The finite difference techniques introduced two methods for calculating the air pollutant concentrations. Figures3 and 4 compare the air pollutant concentrations between two cases: a single point source and coupled point sources, respectively. From the both figures, it is apparent that the results of the forward time central space scheme are close to the results of the backward time central space scheme, when there is no sink of pollutant absorption (
s
=
0
). Figures 5 and 6 illustrate that the sink of pollutant absorption (
s
=
1
0
-
4
) is added to the base of urban zone. The air pollutant concentration near human living goes down and the two methods also give the close result. In Figures 7 and 8, the computed approximate solutions which are calculated by using the forward time central space and the backward time central space schemes are compared. We can see that the results of added sink case and without sink case are quite similar. These graphs also indicate that the forward time central space scheme gives the computed solutions close to the backward time central space scheme.Figure 3
The air pollutant concentration levels after 2 hours passed which are computed by the forward time central space scheme (s
=
0): (a) one-point source and (b) two-point source.
(a)
(b)Figure 4
The air pollutant concentration levels after 2 hours passed which are computed by the backward time central space scheme (s
=
0): (a) one-point source and (b) two-point source.
(a)
(b)Figure 5
The air pollutant concentration levels after 2 hours passed which are computed by the forward time central space scheme (s
=
1
0
-
4): (a) one-point source and (b) two-point source.
(a)
(b)Figure 6
The air pollutant concentration levels after 2 hours passed which are computed by the backward time central space scheme (s
=
1
0
-
4): (a) one-point source and (b) two-point source.
(a)
(b)Figure 7
The air pollutant concentration between the forward time central space and the backward time central space schemes (s
=
0) at z
=
0 m. and x
=
600 m.Figure 8
The air pollutant concentration between the forward time central space and the backward time central space schemes (s
=
1
0
-
4) at z
=
0 m. and x
=
400 m.Figures9 and 10 demonstrate that the air pollutant concentration at the height z
=
25 m. and z
=
50 m. are solved by using the forward time central space scheme. The added sink case is less concentration than the without sink case. Therefore, the sink can lower the overall pollutant levels. Figure 11 establishes the various concentrations when we take more sink rate into our system. The comparison of computing time shows that the forward time central space is faster than the backward time central space scheme in Table 1.Table 1
Computing time comparison of forward time central space and backward time central space schemes.
Simulation Time
FTCS (sec.)
BTCS (sec.)
30 minutes
1.49
22.48
1 hour
1.68
42.66
2 hours
2.05
84.18Figure 9
The air pollutant concentration between 2 cases: added sink and without sink (computed by the forward time central space scheme) atz
=
25 m. and x
=
600 m.Figure 10
The air pollutant concentration between 2 cases: added sink and without sink (computed by the forward time central space scheme) atz
=
50 m. and x
=
600 m.Figure 11
The air pollutant concentration with the variant values of sink rate (computed by the forward time central space scheme) atz
=
0 m. and x
=
600 m.


## 6. Conclusion

The simple air pollution measurement models which are released air pollutants by a single point source and coupled point sources are proposed. The traditional finite difference methods such as forward time central space and backward time central space schemes can be used to approximate the air pollutant levels for each points and times. The results of this study show that the air pollutant concentrations of forward time central space are close to the air pollutant concentrations of backward time central space. In the case of a coupled point sources problem, the overall concentration levels of air pollution are less than a single point source problem. Therefore, the influence of multiple point sources and the variable rate of sink are also considered. It obtains that the higher sink rate does decrease pollutant levels around human living. The both finite difference methods are used to compute the numerical solutions of air pollution by MATLAB. The forward time central space has advantages that the method gives less computing time than the backward time central space computing time. On the other hand, the forward time central space also has disadvantages that are the limitation of the grid spacing due to the stability condition.



---

*Source: 2319831-2019-02-11.xml*