# On FTCS Approach for Box Model of Three-Dimension Advection-Diffusion Equation

Authors: Jeffry Kusuma; Agustinus Ribal; Andi Galsan Mahie
Journal: International Journal of Differential Equations
 (2018)

---

## Abstract

This paper describes a numerical solution for mathematical model of the transport equation in a simple rectangular box domain. The model of street tunnel pollution distribution using two-dimension advection and three-dimension diffusion is solved numerically. Because of the nature of the problem, the model is extended to become three-dimension advection and three-dimension diffusion to study the sea-sand mining pollution distribution. This model with various advection and diffusion parameters and the boundaries conditions is also solved numerically using a finite difference (FTCS) method.


---

## Body

## 1. Introduction

The solution of partial differential equation and their associated boundary and initial condition play an important role in modelling of phenomena in fields as diverse as physics, chemistry, geology, biology, engineering, and economics. The transport of pollutants occurs in a large variety of environmental, agricultural, and industrial processes. The phenomenon is usually modelled into partial differential equations with boundary and/or initial conditions. The models, however, in most cases have no analytical solution. Numerical solution becomes an alternative solution to models such as partial differential equation models in order to investigate, predict, and conclude the models.Numerical solution for an advection-diffusion equation or transport equation had become an interesting subject for many authors recently. Several improvements in finite difference approach had been noted. A stability limit for a finite difference scheme such as the forward time and space-centered numerical scheme applied the convection-diffusion equation is discussed in [1]. Three-dimension solution of advection-diffusion base on the two-level fully explicit and fully implicit finite difference approximation is discussed in [2]. The comparison of the two standard finite difference schemes such as FTCS and Crank-Nicolson methods is carried out by [3]. The method of time splitting to divide complicated time dependent partial differential equation into sets of simpler equations which could then be solved separately by numerical means over fraction time step had been done as in [4]. Numerical solution of a three-dimension advection-diffusion for models in a street tunnel is discussed in [5]. More recently, an application of the generalized finite difference method to solve the advection-diffusion equation using the explicit method is discussed in [6]. Also, practical use of finite difference methods in order to study the pollution distribution on Unhas lake had been carried out as in [7].The studying of a box model had been carried out numerically as in [5] using two-dimension advection and three-dimension diffusion. In this paper we extend the box model into three-dimension advection and three-dimension diffusion to determine the pollutants spread in the water or in the air. This is because of the nature of the pollutant particles. For certain dimension of particle, gravitational force due to the particle mass should be taken into account. We solve the three-dimension advection-diffusion equation by using the forward in time, center in space (FTCS) finite difference method. The domain of this study model is a rectangular box with length L, width W, and height H. Numerical results for several different pollutant source configurations are presented and discussed.The paper is organized as follows. In Section2, we introduce the basic equation and the problems. In Section 3, the finite difference schemes for the computation approach are introduced. The stability condition for FTCS scheme is discussed in Section 4. In Section 5, some numerical models, results, and discussions are presented. Finally, conclusions are found in Section 6.


## 2. Basic Equation and Problems

In this paper we consider the three-dimension advection-diffusion equation [2, 5, 6](1)
∂
C
∂
t
+
V
x
∂
C
∂
x
+
V
y
∂
C
∂
y
+
V
z
∂
C
∂
z
=
D
x
∂
2
C
∂
x
2
+
D
y
∂
2
C
∂
y
2
+
D
z
∂
2
C
∂
z
2
,
0
<
t
≤
T
,in the domains 0
≤
x
≤
L, 0
≤
y
≤
W, and 0
≤
z
≤
H, with initial condition(2)
C
x
,
y
,
z
,
0
=
f
x
,
y
,
z
,and the boundary conditions (3)
C
0
,
y
,
z
,
t
=
g
0
y
,
z
,
t
,
0
<
t
≤
T
,
(4)
C
L
,
y
,
z
,
t
=
g
L
y
,
z
,
t
,
0
<
t
≤
T
,
(5)
C
x
,
0
,
z
,
t
=
h
0
x
,
z
,
t
,
0
<
t
≤
T
,
(6)
C
x
,
W
,
z
,
t
=
h
W
x
,
z
,
t
,
0
<
t
≤
T
,
(7)
C
x
,
y
,
0
,
t
=
k
0
x
,
y
,
t
,
0
<
t
≤
T
,
(8)
C
x
,
y
,
H
,
t
=
k
H
x
,
y
,
t
,
0
<
t
≤
T
,where f
,
g
0
,
g
L
,
h
0
,
h
W
,
k
0, and k
H are known functions, while the function C is unknown. Here C
x
,
y
,
z
,
t denote materials concentration which is transported. The constants V
x
,
V
y
,
V
z represent speeds of advection with respect to x
-
a
x
i
s, y
-
a
x
i
s, and z
-
a
x
i
s, respectively. Also, the constants D
x
,
D
y
,
D
z, represent the speeds of diffusivities with respect to x
-
a
x
i
s, y
-
a
x
i
s, and z
-
a
x
i
s, respectively (see illustration in Figure 1).Figure 1
Domain of advection-diffusion.


## 3. Finite Difference Schemes

The main idea behind the finite difference schemes for obtaining the solution of a given partial differential equation is to approximate the derivatives appearing in the equation by a set of values of the function at selected number of points. The most usual way to generate these approximations is through the use of Taylor series.The solution domain of the problem over a time0
≤
t
≤
T is covered by a mesh of uniformly spaced grid-lines parallel to the space and time coordinates axes, respectively.(9)
x
i
=
i
Δ
x
,
i
=
0,1
,
2
,
…
,
M
;
(10)
y
j
=
j
Δ
y
,
j
=
0,1
,
2
,
…
,
N
;
(11)
z
k
=
k
Δ
z
,
k
=
0,1
,
2
,
…
,
O
;
(12)
t
n
=
n
Δ
t
,
n
=
0,1
,
2
,
…
,
R
;Approximations C
i
,
j
,
k
n to C
(
i
Δ
x
,
j
Δ
y
,
k
Δ
z
,
n
Δ
t
) are calculated at the point of intersection of these lines according to the i
,
j
,
k
,
n grid points. The uniform spatial and temporal grid spacings are Δ
x
=
L
/
M, Δ
y
=
W
/
N, Δ
z
=
H
/
O, and Δ
t
=
T
/
R, where L is the length, W is the width, and H is the height of the domain of the interested rectangular box.Since the grid points are in three dimensions of the box form, several numerical approximations are needed for the grid points. These depend on the position of the grid points. The approximation for interior points is forward time center space (FTCS). The points obtained from intersection of two planes of the box are approximated using forward time center space forward or backward space while the edge points with intersection of three planes are approximated using forward time and combination of forward and backward approximation.Forward time centered space (FTCS) approximation of (1) for the interior points as in [2, 5, 7] is(13)
C
i
,
j
,
k
n
+
1
-
C
i
,
j
,
k
n
Δ
t
+
V
x
C
i
+
1
,
j
,
k
n
-
C
i
-
1
,
j
,
k
n
2
Δ
x
+
V
y
C
i
,
j
+
1
,
k
n
-
C
i
,
j
-
1
,
k
n
2
Δ
y
+
V
z
C
i
,
j
,
k
+
1
n
-
C
i
,
j
,
k
-
1
n
2
Δ
z
=
D
x
C
i
+
1
,
j
,
k
n
-
2
C
i
,
j
,
k
n
+
C
i
-
1
,
j
,
k
n
Δ
x
2
+
D
y
C
i
,
j
+
1
,
k
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
-
1
,
k
n
Δ
y
2
+
D
z
C
i
,
j
,
k
+
1
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
,
k
-
1
n
Δ
z
2
.Another finite difference approximation such as forward time backward space x center space y center space z is used for the boundary points on C
L
,
y
,
z
,
t
=
g
L
y
,
z
,
t
(14)
C
i
,
j
,
k
n
+
1
-
C
i
,
j
,
k
n
Δ
t
+
V
x
3
C
i
,
j
,
k
n
-
4
C
i
-
1
,
j
,
k
n
+
3
C
i
-
2
,
j
,
k
n
2
Δ
x
+
V
y
C
i
,
j
+
1
,
k
n
-
C
i
,
j
-
1
,
k
n
2
Δ
y
+
V
z
C
i
,
j
,
k
+
1
n
-
C
i
,
j
,
k
-
1
n
2
Δ
z
=
D
x
C
i
,
j
,
k
n
-
2
C
i
-
1
,
j
,
k
n
+
C
i
-
2
,
j
,
k
n
Δ
x
2
+
D
y
C
i
,
j
+
1
,
k
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
-
1
,
k
n
Δ
y
2
+
D
z
C
i
,
j
,
k
+
1
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
,
k
-
1
n
Δ
z
2
.Forward time backward space x backward space y center space z is used for the boundary points on C
L
,
W
,
z
,
t which is(15)
C
i
,
j
,
k
n
+
1
-
C
i
,
j
,
k
n
Δ
t
+
V
x
3
C
i
,
j
,
k
n
-
4
C
i
-
1
,
j
,
k
n
+
3
C
i
-
2
,
j
,
k
n
2
Δ
x
+
V
y
3
C
i
,
j
,
k
n
-
4
C
i
,
j
-
1
,
k
n
+
3
C
i
,
j
-
2
,
k
n
2
Δ
y
+
V
z
C
i
,
j
,
k
+
1
n
-
C
i
,
j
,
k
-
1
n
2
Δ
z
=
D
x
C
i
,
j
,
k
n
-
2
C
i
-
1
,
j
,
k
n
+
C
i
-
2
,
j
,
k
n
Δ
x
2
+
D
y
C
i
,
j
,
k
n
-
2
C
i
,
j
-
1
,
k
n
+
C
i
,
j
-
2
,
k
n
Δ
y
2
+
D
z
C
i
,
j
,
k
+
1
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
,
k
-
1
n
Δ
z
2
.Also forward time backward space x forward space y center space z is used for the boundary points on C
L
,
0
,
z
,
t which is(16)
C
i
,
j
,
k
n
+
1
-
C
i
,
j
,
k
n
Δ
t
+
V
x
3
C
i
,
j
,
k
n
-
4
C
i
-
1
,
j
,
k
n
+
3
C
i
-
2
,
j
,
k
n
2
Δ
x
+
V
y
C
i
,
j
,
k
n
-
4
C
i
,
j
+
1
,
k
n
+
3
C
i
,
j
+
2
,
k
n
2
Δ
y
+
V
z
C
i
,
j
,
k
+
1
n
-
C
i
,
j
,
k
-
1
n
2
Δ
z
=
D
x
C
i
,
j
,
k
n
-
2
C
i
-
1
,
j
,
k
n
+
C
i
-
2
,
j
,
k
n
Δ
x
2
+
D
y
C
i
,
j
,
k
n
-
2
C
i
,
j
+
1
,
k
n
+
C
i
,
j
+
2
,
k
n
Δ
y
2
+
D
z
C
i
,
j
,
k
+
1
n
-
2
C
i
,
j
,
k
n
+
C
i
,
j
,
k
-
1
n
Δ
z
2
.All approximations have error of first order in time interval and second order is spatial coordinate grid spacing.


## 4. Stability

To ensure obtained solutions have a nonpropagate error, the approximation or schemes should meet certain conditions. The approximation schemes on the boundary usually have a nonpropagate error since boundary exact conditions are supplied on the boundary. However, for the interior points the error might propagate. For the FTCS approximations, (13) might be rewritten as(17)
C
i
,
j
,
k
n
+
1
=
D
x
Δ
t
Δ
x
2
+
V
x
Δ
t
2
Δ
x
C
i
-
1
,
j
,
k
n
+
D
y
Δ
t
Δ
y
2
+
V
y
Δ
t
2
Δ
y
C
i
,
j
-
1
,
k
n
+
D
z
Δ
t
Δ
z
2
+
V
z
Δ
t
2
Δ
z
C
i
,
j
,
k
-
1
n
+
D
x
Δ
t
Δ
x
2
-
V
x
Δ
t
2
Δ
x
C
i
+
1
,
j
,
k
n
+
D
y
Δ
t
Δ
y
2
-
V
y
Δ
t
2
Δ
y
C
i
,
j
+
1
,
k
n
+
D
z
Δ
t
Δ
z
2
-
V
z
Δ
t
2
Δ
z
C
i
,
j
,
k
-
1
n
+
1
-
2
D
x
Δ
t
Δ
x
2
+
D
y
Δ
t
Δ
y
2
+
D
z
Δ
t
Δ
z
2
C
i
,
j
,
k
nor simply (18)
C
i
,
j
,
k
n
+
1
=
s
x
+
c
x
2
C
i
-
1
,
j
,
k
n
+
s
y
+
c
y
2
C
i
,
j
-
1
,
k
n
+
s
z
+
c
z
2
C
i
,
j
,
k
-
1
n
+
s
x
-
c
x
2
C
i
+
1
,
j
,
k
n
+
s
y
-
c
y
2
C
i
,
j
+
1
,
k
n
+
s
z
-
c
z
2
C
i
,
j
,
k
+
1
n
+
1
-
2
s
x
+
s
y
+
s
z
C
i
,
j
,
k
nin which (19)
s
x
=
D
x
Δ
t
Δ
x
2
,
s
y
=
D
y
Δ
t
Δ
y
2
,
s
z
=
D
z
Δ
t
Δ
z
2
,
(20)
c
x
=
V
x
Δ
t
Δ
x
,
c
y
=
V
y
Δ
t
Δ
y
,
c
z
=
V
z
Δ
t
Δ
z
.The stability of this three-dimensional difference scheme may be investigated using the von Neumann method. As in [2, 5] application of this method of stability analysis shows that (18) will be stable if it satisfies both equations(21)
s
x
+
s
y
+
s
z
≤
1
2and(22)
c
x
2
s
x
+
c
y
2
s
y
+
c
z
2
s
z
≤
3
.For one-dimension version of the forward time centered space (FTCS) type formula for the case that s
x
=
s
y
=
s
z
=
s and c
x
=
c
x
=
c
x
=
c, it should also satisfy(23)
c
2
≤
s
≤
1
6
.This clearly is much more strict than one-dimension advection-diffusion equation of the forward time centered space (FTCS) approximation which is s
≤
1
/
2
.


## 5. Numerical Results and Discussions

Problem 1.
In order to verify the accuracy of the procedure that has been built, we consider the three-dimension advection-diffusion equation (1), with initial condition (2) and boundary conditions (3)-(8). By taking V
y
=
V
z
=
0, L
=
1
,
W
=
1
,
H
=
1, this procedure simply reduced to three-dimension advection-diffusion equation for transport of pollutants in street tunnel as discussed in [5](24)
∂
C
∂
t
+
V
x
∂
C
∂
x
=
D
x
∂
2
C
∂
x
2
+
D
y
∂
2
C
∂
y
2
+
D
z
∂
2
C
∂
z
2
,
0
<
t
<
T
.Furthermore, taking appropriate function for f
,
g
0
,
g
L
,
h
0
,
h
W
,
k
0
,
k
H, the initial condition and boundary conditions yield(25)
C
x
,
y
,
z
,
0
=
0
0
≤
x
≤
1
;
0
≤
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
0
0
≤
y
<
0.5
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
1
0.5
≤
y
≤
1
;
0
≤
z
≤
1
C
1
,
y
,
z
,
t
=
0
0
≤
y
≤
1
;
0
≤
z
≤
1
∂
C
∂
y
x
,
0
,
z
,
t
=
0
t
>
0
∂
C
∂
y
x
,
1
,
z
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
0
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
1
,
t
=
0
t
>
0This partial differential equation and its initial and boundary conditions come from pollution distribution on a street tunnel, where the wind flows steadily only in the x direction. There is no disperse flux of the pollutant through the solid side-walls nor through the base and the roof. By using nondimensionalised parameters Δ
x
=
Δ
y
=
Δ
z
=
0.1, Δ
t
=
0.01, D
x
=
D
y
=
D
z
=
0.1, V
x
=
0.02, and for the time T
=
20. These given values will certainly satisfy both conditions s
x
+
s
y
+
s
z
=
3
/
10
≤
1
/
2 and c
x
2
/
s
x
+
c
y
2
/
s
y
+
c
z
2
/
s
z
=
3
/
25000
≤
3
. Thus numerical results are stable and can be found as in the Figure 2.Figure 2
Pollutant distribution in a street tunnel with advection only inx direction.There is no significant difference from the results obtained from [5], as from Figure 2 left one can see the mesh plot on z
=
0. This is reasonable, with the boundary conditions satisfied, the concentration decreases away from the source and is less than one-half of the source value over more than three-quarters of the tunnel. Also, the solutions are independent of height since all contour plots on z
=
0
,
z
=
0.2
,
…
,
z
=
1 are the same as Figure 2 right. The difference may be noted from the previous work as in [5] lies only on x
=
0.1.Problem 2.
Here, we consider the three-dimension advection-diffusion equation (1), with initial condition (2) and boundary conditions (3)-(8). By taking V
z
=
0, L
=
1
,
W
=
1
,
H
=
1, this simply reduces to two-dimension advection three-dimension diffusion equation for transport of pollutants in street tunnel problem as discussed in [5](26)
∂
C
∂
t
+
V
x
∂
C
∂
x
+
V
y
∂
C
∂
y
=
D
x
∂
2
C
∂
x
2
+
D
y
∂
2
C
∂
y
2
+
D
z
∂
2
C
∂
z
2
,
0
<
t
<
T
.Furthermore, taking appropriate function for f
,
g
0
,
g
L
,
h
0
,
h
W
,
k
0
,
k
H, the initial condition and boundary conditions yield(27)
C
x
,
y
,
z
,
0
=
0
0
≤
x
≤
1
;
0
≤
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
0
0
≤
y
<
0.5
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
1
0.5
≤
y
≤
1
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
0
0
≤
y
<
0.3
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
1
0.3
≤
y
≤
0.6
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
0
0.6
<
y
≤
1
;
0
≤
z
≤
1
∂
C
∂
y
x
,
0
,
z
,
t
=
0
t
>
0
∂
C
∂
y
x
,
1
,
z
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
0
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
1
,
t
=
0
t
>
0This model of partial differential equation and its initial and boundary conditions comes from pollution distribution on a street tunnel, where the wind flows steadily only in the x and y directions. There is no disperse flux of the pollutant through the solid side-walls nor through the base and the roof. By using nondimensionalised parameters Δ
x
=
Δ
y
=
Δ
z
=
0.1, Δ
t
=
0.005, D
x
=
D
y
=
D
z
=
0.2, V
x
=
0.6, V
y
=
0.4 and for the time T
=
20. These values satisfy both conditions s
x
+
s
y
+
s
z
=
3
/
10
≤
1
/
2 and c
x
2
/
s
x
+
c
y
2
/
s
y
+
c
z
2
/
s
z
=
13
/
1000
≤
3
. The numerical results are found as in Figure 3.Figure 3
Pollutant distribution in a street tunnel with advection inx and y directions.Figure3 shows a mesh plot and a contour plot on z
=
0.2. The other mesh plot and contour plot on the other z show no difference in z
=
0.2. This is reasonable, with the boundary conditions satisfied, the concentration decreases away from the source faster in the x direction rather than in the y direction.Problem 3.
We consider a three-dimensional advection-diffusion equation for transport of pollutants in the water such as a sea-sand mining or in the air such as limestone rock mining to produce cement. Such activity usually spread pollutants to the environment nearby. By assuming that the domain of the physical problem had been nondimensionalized, and taking the box model, it will satisfy three-dimension advection and three-dimension diffusion equation. There is no disperse flux of the pollutant through the base and the roof. The boundary condition on top of the box will be sea surface and the bottom of the box will be sea floor and satisfy∂
C
/
∂
z
x
,
y
,
H
,
t
=
0 and ∂
C
/
∂
z
x
,
y
,
0
,
t
=
0
. The advection constant V
z
≠
0 due to the gravitational force acting on the particle. The constants V
x and V
y are the sea current speed in x and y directions. Another initial condition and boundary conditions yield(28)
C
x
,
y
,
z
,
0
=
0
0
≤
x
≤
1
;
0
≤
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
0
0
≤
y
<
0.4
;
0.6
<
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
1
0.4
≤
y
≤
0.6
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
0
0
≤
x
<
0.4
;
0.6
<
x
≤
1
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
1
0.4
≤
x
≤
0.6
;
0
≤
z
≤
1
C
1
,
y
,
z
,
t
=
0
0
≤
y
≤
1
;
0
≤
z
≤
1
∂
C
∂
y
x
,
0
,
z
,
t
=
0
t
>
0
∂
C
∂
y
x
,
1
,
z
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
0
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
1
,
t
=
0
t
>
0Using nondimensionalised parameters Δ
x
=
Δ
y
=
Δ
z
=
0.1, Δ
t
=
0.005; D
x
=
D
y
=
D
z
=
0.1, V
x
=
1.0, V
y
=
0.3, V
z
=
-
0.3 and for the time T
=
2
. Here, the values satisfy also both stability conditions s
x
+
s
y
+
s
z
=
3
/
20
≤
1
/
2 and c
x
2
/
s
x
+
c
y
2
/
s
y
+
c
z
2
/
s
z
=
59
/
1000
≤
3
. The numerical results for surface plot and contour plot on z
=
0.5 and several values of time t are found as in Figure 4. Surface plot and contour plot for other values of z are similar as in Figure 4.Figure 4
Pollutant distribution in a box model under the sea with advection inx, y, and z directions.Problem 4.
We consider a three-dimensional advection-diffusion equation for transport of pollutants in the water such as a sea-sand mining again. Such activity usually lasts for certain period of time and spreads pollutants to the environment nearby. Of our interest is to find out the pollutant distribution of this activity. By assuming that the domain of the physical problem had been nondimensionalized and taking the box model, it will satisfy three-dimension advection and three-dimension diffusion equation. There is no disperse flux of the pollutant through the base and the roof. The boundary condition on top of the box will be sea surface and the bottom of the box will be sea floor and satisfy∂
C
/
∂
z
x
,
y
,
H
,
t
=
0 and ∂
C
/
∂
z
x
,
y
,
0
,
t
=
0
. Other initial condition and boundary conditions yield(29)
C
x
,
y
,
z
,
0
=
0
0
≤
x
≤
1
;
0
≤
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
0
0
≤
y
<
0.4
;
0.6
<
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
1
0.4
≤
y
≤
0.6
;
0
≤
z
≤
1
;
0
≤
t
≤
0.5
C
x
,
0
,
z
,
t
=
0
0
≤
x
<
0.4
;
0.6
<
x
≤
1
;
0
≤
z
≤
1
C
x
,
0
,
z
,
t
=
1
0.4
≤
x
≤
0.6
;
0
≤
z
≤
1
;
0
≤
t
≤
0.5
C
1
,
y
,
z
,
t
=
0
0
≤
y
≤
1
;
0
≤
z
≤
1
C
0
,
y
,
z
,
t
=
0
0.4
≤
y
≤
0.6
;
0
≤
z
≤
1
;
t
>
0.5
C
x
,
0
,
z
,
t
=
0
0.4
≤
x
≤
0.6
;
0
≤
z
≤
1
;
t
>
0.5
∂
C
∂
y
x
,
0
,
z
,
t
=
0
t
>
0
∂
C
∂
y
x
,
1
,
z
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
0
,
t
=
0
t
>
0
∂
C
∂
z
x
,
y
,
1
,
t
=
0
t
>
0Using nondimensionalised parameters Δ
x
=
Δ
y
=
Δ
z
=
0.1, Δ
t
=
0.005; D
x
=
D
y
=
D
z
=
0.2, V
x
=
1.0, V
y
=
0.3, V
z
=
-
0.3 and for the time T
=
2
. These values also satisfy both stability conditions s
x
+
s
y
+
s
z
=
3
/
10
≤
1
/
2 and c
x
2
/
s
x
+
c
y
2
/
s
y
+
c
z
2
/
s
z
=
295
/
10000
≤
3
. There are two pollutant sources on the x-axis and y-axis and they last for certain period of time 0
≤
t
≤
0.5
. The numerical results for surface plot and contour plot on z
=
0.4for several different times are found as in Figure 5. Unlike the previous Problem 3, surface and contour plot for other values of z show that pollutant materials move to x
,
y
,
z directions as expected.Figure 5
Pollutant distribution in a box model under the sea with advection inx, y, and z directions.


## 6. Conclusions

We have employed a standard finite difference scheme to study the pollution distribution for two-dimension advection and three-dimension diffusion equation and extend into three-dimension advection and three-dimension diffusion equation. The stability conditions ensure that the solution of all the interior points can be obtained. The scheme works well as one can see numerical results obtained in solutions of the problems above.



---

*Source: 7597861-2018-11-01.xml*