I
nte
rna
t
io
na
l J
o
urna
l o
f
E
lect
rica
l a
nd
Co
m
p
ute
r
E
ng
in
ee
ring
(
I
J
E
CE
)
Vo
l.
7
,
No
.
3
,
J
u
n
e
201
7
,
p
p
.
1522
~
1
5
2
9
I
SS
N:
2
0
8
8
-
8708
,
DOI
: 1
0
.
1
1
5
9
1
/
i
j
ec
e
.
v
7
i
3
.
p
p
1
5
2
2
-
1
529
1522
J
o
ur
na
l ho
m
ep
a
g
e
:
h
ttp
:
//ia
e
s
jo
u
r
n
a
l.c
o
m/o
n
lin
e/in
d
ex
.
p
h
p
/I
JE
C
E
P
a
ra
m
e
ter
O
p
ti
mis
a
tion
o
f
Aerob
ic
G
ra
nula
r
Sludg
e
a
t
H
ig
h
Te
m
p
erature
Usi
ng
Respo
nse Sur
f
a
ce Methodo
lo
g
y
Sy
a
hira
I
bra
hi
m
1
,
No
rha
liza
Abdu
l Wa
ha
b
2
,
Azna
h No
r
Anua
r
3
,
M
us
t
a
f
a
B
o
b
4
1,
2
De
p
a
rtm
e
n
t
o
f
Co
n
tro
l
&
M
e
c
h
a
tro
n
ics
En
g
in
e
e
rin
g
,
F
a
c
u
lt
y
o
f
El
e
c
tri
c
a
l
En
g
in
e
e
rin
g
,
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
lay
sia
,
8
1
3
1
0
S
k
u
d
a
i
Jo
h
o
r
Ba
h
ru
,
M
a
lay
s
ia
3
De
p
a
rtm
e
n
t
o
f
En
v
iro
n
m
e
n
tal
En
g
in
e
e
rin
g
,
F
a
c
u
lt
y
o
f
Civ
il
En
g
in
e
e
rin
g
,
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
la
y
sia
4
De
p
a
rtm
e
n
t
o
f
Civ
il
En
g
in
e
e
rin
g
,
Co
ll
e
g
e
o
f
En
g
in
e
e
rin
g
,
T
a
ib
a
h
Un
iv
e
rsity
,
M
a
d
in
a
h
Cit
y
,
Kin
g
d
o
m
o
f
S
a
u
d
i
A
ra
b
ia
Art
icle
I
nfo
AB
ST
RAC
T
A
r
ticle
his
to
r
y:
R
ec
eiv
ed
J
an
26
,
2
0
1
7
R
ev
i
s
ed
Ap
r
1
,
2
0
1
7
A
cc
ep
ted
A
p
r
2
6
,
2
0
1
7
T
h
is
p
a
p
e
r
p
ro
p
o
se
s
a
n
im
p
ro
v
e
d
o
p
ti
m
isa
ti
o
n
o
f
se
q
u
e
n
c
i
n
g
b
a
t
c
h
re
a
c
to
rs
(
S
BR)
f
o
r
a
e
ro
b
ic
g
ra
n
u
lar
slu
d
g
e
(AG
S
)
a
t
h
ig
h
tem
p
e
r
a
tu
re
-
lo
w
h
u
m
id
it
y
f
o
r
d
o
m
e
stic
wa
ste
wa
ter
trea
t
m
e
n
t
u
sin
g
re
sp
o
n
se
su
rf
a
c
e
m
e
th
o
d
o
l
o
g
y
(RS
M
).
T
h
e
m
a
in
a
d
v
a
n
tag
e
s
o
f
RS
M
a
re
le
ss
n
u
m
b
e
r
o
f
e
x
p
e
rime
n
t
re
q
u
ired
a
n
d
su
it
a
b
le
f
o
r
c
o
m
p
le
x
p
ro
c
e
ss
.
T
h
e
slu
d
g
e
f
ro
m
a
c
o
n
v
e
n
ti
o
n
a
l
a
c
ti
v
a
ted
slu
d
g
e
w
a
ste
wa
ter
tre
a
tme
n
t
p
lan
t
a
n
d
t
h
re
e
se
q
u
e
n
c
in
g
b
a
tch
re
a
c
to
rs
(S
BRs
)
w
e
re
f
e
d
w
it
h
sy
n
th
e
ti
c
w
a
st
e
w
a
ter.
T
h
e
e
x
p
e
rime
n
t
we
r
e
c
a
rried
o
u
t
a
t
d
if
fe
re
n
t
h
ig
h
t
e
m
p
e
ra
tu
re
s
(3
0
,
4
0
a
n
d
5
0
°
C)
a
n
d
th
e
f
o
r
m
a
ti
o
n
o
f
AG
S
f
o
r
si
m
u
lt
a
n
e
o
u
s
o
rg
a
n
ics
a
n
d
n
u
tri
e
n
ts
re
m
o
v
a
l
we
re
e
x
a
m
in
e
d
in
6
0
d
a
y
s.
RS
M
is
u
se
d
to
m
o
d
e
l
a
n
d
t
o
o
p
ti
m
ize
th
e
b
io
l
o
g
ica
l
p
a
ra
m
e
ters
f
o
r
c
h
e
m
ica
l
o
x
y
g
e
n
d
e
m
a
n
d
(COD
)
a
n
d
t
o
tal
p
h
o
sp
h
o
r
u
s
re
m
o
v
a
l
in
S
BR s
y
ste
m
.
T
h
e
si
m
u
latio
n
re
su
lt
s sh
o
w
e
d
th
a
t
a
t
te
m
p
e
ra
tu
re
o
f
4
5
.
3
3
°
C
g
iv
e
th
e
o
p
ti
m
u
m
c
o
n
d
it
io
n
f
o
r
th
e
to
tal
re
m
o
v
a
l
o
f
COD
a
n
d
p
h
o
s
p
h
o
ru
s,
w
h
ich
c
o
rre
sp
o
n
d
t
o
p
e
rf
o
rm
a
n
c
e
in
d
e
x
R
2
o
f
0
.
9
5
5
a
n
d
0
.
9
1
,
re
sp
e
c
ti
v
e
l
y
.
K
ey
w
o
r
d
:
A
er
o
b
ic
g
r
an
u
lar
s
l
u
d
g
e
Hig
h
t
e
m
p
er
at
u
r
e
Op
ti
m
izatio
n
R
esp
o
n
s
e
s
u
r
f
ac
e
m
eth
o
d
o
lo
g
y
Seq
u
en
ci
n
g
b
atch
r
ea
cto
r
Co
p
y
rig
h
t
©
2
0
1
7
In
stit
u
te o
f
A
d
v
a
n
c
e
d
E
n
g
i
n
e
e
rin
g
a
n
d
S
c
ien
c
e
.
Al
l
rig
h
ts
re
se
rv
e
d
.
C
o
r
r
e
s
p
o
nd
ing
A
uth
o
r
:
No
r
h
aliza
A
b
d
u
l W
ah
ab
,
Dep
ar
t
m
en
t o
f
C
o
n
tr
o
l &
Me
c
h
atr
o
n
ics E
n
g
in
ee
r
i
n
g
,
Un
i
v
er
s
iti T
ek
n
o
lo
g
i M
ala
y
s
ia
,
8
1
3
1
0
Sk
u
d
ai
J
o
h
o
r
B
ah
r
u
,
Ma
la
y
s
ia
.
E
m
ail:
aliza
@
f
k
e.
u
t
m
.
m
y
1.
I
NT
RO
D
UCT
I
O
N
Seq
u
en
ci
n
g
b
atch
r
ea
cto
r
(
SB
R
)
is
a
m
o
d
if
ied
ac
ti
v
ated
s
lu
d
g
e
p
r
o
ce
s
s
t
h
at
h
as
b
ee
n
ap
p
lied
f
o
r
n
u
m
er
o
u
s
w
as
te
w
ater
tr
ea
t
m
en
t
p
lan
ts
s
u
c
h
as
s
o
y
s
au
ce
[
1
]
,
d
air
y
[
2
]
,
in
d
u
s
tr
ial
est
ate
[
3
]
,
tex
tile
[
4
]
s
e
w
ag
e,
s
y
n
t
h
etic
w
aste
w
ater
[5
]
,
p
alm
o
il
m
ill
e
f
f
l
u
en
t
(
P
OM
E
)
[
6
]
an
d
s
o
o
n
.
T
h
is
s
y
s
te
m
is
o
p
er
ated
in
a
s
in
g
le
r
ea
cto
r
a
n
d
it
h
a
s
f
i
v
e
b
asic
o
p
er
atin
g
m
o
d
es
in
c
lu
d
i
n
g
f
ill,
r
ea
ct,
s
ettle,
d
r
a
w
a
n
d
id
le
[
7
]
.
C
o
m
p
ar
ed
w
it
h
t
h
e
co
s
tl
y
a
n
d
lar
g
e
-
s
ize
o
f
co
n
v
en
tio
n
al
ac
ti
v
ated
s
l
u
d
g
e
p
r
o
ce
s
s
,
SB
R
is
a
f
u
ll
b
i
o
l
o
g
ical
tr
ea
t
m
e
n
t,
less
p
o
w
er
co
n
s
u
m
p
tio
n
,
r
eq
u
ir
ed
n
o
c
h
e
m
ical
ad
d
itio
n
a
n
d
its
o
p
er
ates
i
n
t
i
m
e
r
ath
er
th
a
n
i
n
s
p
ac
e
[
8
]
.
I
n
d
u
s
tr
ial
w
aste
w
ater
co
n
tai
n
s
v
ar
io
u
s
ar
o
m
atic
h
y
d
r
o
ca
r
b
o
n
s
an
d
i
n
o
r
g
a
n
ic
s
u
b
s
ta
n
ce
s
s
u
c
h
as
c
h
e
m
ical
o
x
y
g
en
d
e
m
an
d
(
C
OD)
,
p
h
o
s
p
h
o
r
u
s
(
P
)
,
n
itro
g
en
(
N)
an
d
a
m
m
o
n
ia
n
itro
g
e
n
(
NH
3
-
N)
p
lu
s
it
v
ar
y
i
n
g
i
n
p
H,
te
m
p
er
atu
r
e
a
n
d
co
lo
r
.
E
x
ce
s
s
iv
e
o
f
th
e
s
e
n
u
tr
ien
t
s
i
n
aq
u
at
i
c
w
a
ter
lead
s
to
s
e
v
er
e
p
r
o
b
lem
o
f
eu
tr
o
p
h
ica
tio
n
o
n
th
e
w
ater
s
u
r
f
ac
e
f
o
r
lak
es
an
d
o
ce
an
s
[
9
]
.
T
h
e
d
ev
elo
p
m
en
t
o
f
ae
r
o
b
ic
g
r
an
u
lar
s
lu
d
g
e
(
A
G
S)
in
SB
R
(
A
G
S
-
SB
R
)
h
as
b
ee
n
r
ep
o
r
ted
s
in
ce
1
9
9
0
s
.
A
GS
is
m
ad
e
u
p
o
f
a
d
en
s
e
cl
u
s
ter
o
f
s
y
m
b
io
tic
o
r
g
a
n
is
m
s
w
i
th
g
o
o
d
b
io
l
o
g
ical
ac
tiv
it
y
p
er
f
o
r
m
a
n
ce
an
d
ex
ce
l
len
t
m
a
s
s
tr
an
s
f
er
e
f
f
icie
n
c
y
.
Us
u
all
y
,
A
G
S
ca
n
b
e
cu
lti
v
ated
u
s
i
n
g
s
in
g
le
r
ea
ct
o
r
o
f
SB
R
s
y
s
te
m
s
w
it
h
air
li
f
t o
r
b
u
b
b
le
co
l
u
m
n
r
ea
cto
r
s
b
ec
au
s
e
n
i
tr
if
ica
tio
n
,
d
en
itri
f
icatio
n
a
n
d
b
io
lo
g
ical
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
co
u
l
d
ca
r
r
ied
o
u
t
w
ith
in
t
h
e
g
r
an
u
le
s
[
9
]
.
T
h
er
e
is
a
s
tu
d
y
w
i
th
i
n
lab
s
ca
le
t
h
at
i
s
d
e
m
o
n
s
tr
ated
to
id
en
tify
t
h
e
m
o
s
t
v
i
tal
f
ea
t
u
r
e
th
at
in
f
l
u
e
n
ci
n
g
t
h
e
d
ev
elo
p
m
e
n
t
o
f
AGS
s
u
c
h
as
o
r
g
an
ic
lo
ad
in
g
r
ate,
h
y
d
r
o
d
y
n
a
m
ic
s
h
ea
r
Evaluation Warning : The document was created with Spire.PDF for Python.
I
J
E
C
E
I
SS
N:
2
0
8
8
-
8708
P
a
r
a
mete
r
Op
timis
a
tio
n
o
f A
ero
b
ic
Gra
n
u
la
r
S
lu
d
g
e
a
t H
ig
h
Temp
era
tu
r
e
Usi
n
g
.
.
.
.
(
S
ya
h
ir
a
I
b
r
a
h
im)
1523
f
o
r
ce
,
s
ettli
n
g
ti
m
e,
a
n
d
s
u
b
s
t
r
ate
co
m
p
o
s
it
io
n
[
1
0
]
.
A
G
S
f
o
r
m
at
io
n
is
a
co
m
p
le
x
ec
o
lo
g
ical
p
r
o
ce
s
s
an
d
it
m
a
y
in
f
l
u
en
ce
d
b
y
o
t
h
er
p
ar
a
m
eter
s
t
h
at
n
ee
d
to
b
e
f
u
r
t
h
er
ex
a
m
i
n
ed
.
Or
d
in
ar
il
y
t
h
e
s
tu
d
ies
o
f
ae
r
o
b
ic
g
r
an
u
lar
s
l
u
d
g
e
u
s
i
n
g
SB
R
s
y
s
te
m
(
A
G
S
-
SB
R
)
w
er
e
p
er
f
o
r
m
ed
at
2
0
–
2
5
°C
[
5
]
o
r
lo
w
er
[
1
0
]
.
T
h
er
ef
o
r
e,
it is
n
o
t
s
u
itab
le
f
o
r
th
e
co
u
n
tr
y
w
i
th
h
o
t
w
ea
th
er
o
r
cli
m
ate
lik
e
s
A
r
ab
Sau
d
i.
I
t
is
a
n
ee
d
o
f
f
u
r
th
er
in
v
es
tig
a
tio
n
o
n
th
e
e
f
f
ec
t
o
f
ae
r
o
b
ic
g
r
an
u
l
atio
n
p
r
o
ce
s
s
at
d
i
f
f
er
en
t
h
i
g
h
te
m
p
er
atu
r
e
(
3
0
°C
,
4
0
°C
an
d
5
0
°C
)
u
s
i
n
g
SB
R
s
y
s
te
m
.
Mo
d
elin
g
o
f
t
h
e
n
u
tr
ien
t
r
e
m
o
v
al
i
n
SB
R
s
y
s
te
m
is
v
ital
b
u
t
to
d
e
v
elo
p
th
e
m
o
d
el
u
s
i
n
g
cla
s
s
ica
l
m
at
h
e
m
a
tical
ap
p
r
o
ac
h
is
co
m
p
licated
.
I
t
is
b
ec
a
u
s
e,
p
r
o
d
u
ctio
n
o
f
A
G
S
i
n
v
o
l
v
ed
a
v
e
r
y
co
m
p
le
x
p
r
o
ce
s
s
w
it
h
m
an
y
in
f
l
u
en
tial
f
ac
to
r
s
.
Usu
al
l
y
,
co
n
v
en
t
io
n
al
o
p
ti
m
iz
atio
n
m
et
h
o
d
w
as
ca
r
r
ied
o
u
t
b
y
v
ar
y
i
n
g
a
s
i
n
g
le
f
ac
to
r
w
h
ile
k
ee
p
i
n
g
o
th
er
f
ac
to
r
s
f
ix
ed
at
a
s
p
ec
if
ic
v
al
u
e
.
I
n
th
is
ca
s
e,
r
esp
o
n
s
e
s
u
r
f
ac
e
m
e
t
h
o
d
(
R
SM)
is
a
m
o
n
g
th
e
b
est
tec
h
n
iq
u
e
to
d
esig
n
th
e
m
o
d
el
an
d
to
f
in
d
t
h
e
o
p
tim
u
m
v
al
u
e
o
f
r
esp
o
n
s
e
s
.
T
h
is
is
d
u
e
to
th
e
ca
p
ab
ilit
y
o
f
R
SM
in
e
s
ti
m
ati
n
g
th
e
e
f
f
ec
ts
o
f
s
ev
er
al
f
ac
to
r
s
an
d
ex
a
m
i
n
i
n
g
o
p
ti
m
u
m
co
n
d
itio
n
s
f
o
r
d
esire
d
r
esp
o
n
s
es
[
1
1
]
.
Sev
er
al
s
tu
d
ie
s
also
e
m
p
lo
y
ed
R
SM
i
n
SB
R
f
o
r
th
e
te
x
tile d
y
e
w
a
s
te
w
ater
tr
ea
t
m
e
n
t
[
4
]
,
p
alm
o
il
m
ill
[
6
]
,
s
y
n
t
h
etic
w
a
s
te
w
a
ter
[
1
2
]
an
d
h
as
s
h
o
w
n
g
o
o
d
r
esu
lt
s
o
f
m
o
d
elin
g
.
Fro
m
th
e
i
d
en
tifie
d
m
o
d
el,
it
ca
n
b
e
s
ee
n
t
h
at
t
h
er
e
w
er
e
s
o
m
e
li
m
itatio
n
in
m
ea
s
u
r
i
n
g
t
h
e
d
y
n
a
m
ic
o
f
t
h
e
b
io
lo
g
ical
p
r
o
ce
s
s
[
1
3
]
.
B
esid
es,
it
w
i
ll h
elp
t
h
e
d
ev
elo
p
m
e
n
t o
f
a
s
i
m
p
le
m
o
d
el
in
o
r
d
er
to
co
n
tr
o
l th
e
SB
R
s
y
s
te
m
i
n
a
n
ef
f
icie
n
t
w
a
y
.
T
h
e
f
o
r
m
at
io
n
o
f
A
G
S
in
SB
R
p
r
o
ce
s
s
at
h
ig
h
te
m
p
er
atu
r
e
h
as
b
ee
n
s
t
u
d
ied
in
[
5
]
.
T
h
en
,
[
1
3
]
h
as
f
u
r
t
h
er
in
v
es
tig
a
ted
th
e
A
GS
p
r
o
ce
s
s
u
s
i
n
g
Ne
u
r
al
Net
w
o
r
k
w
it
h
Gr
av
itat
i
o
n
al
s
ea
r
c
h
al
g
o
r
ith
m
(
GS
A
)
an
d
I
n
er
tia
w
eig
h
t
p
ar
ticle
s
w
ar
m
o
p
tim
izatio
n
(
I
W
-
P
SO)
to
u
n
d
er
s
tan
d
t
h
e
d
y
n
a
m
ic
b
eh
a
v
io
u
r
o
f
th
e
A
G
S.
T
h
is
p
ap
er
aim
s
to
m
o
d
el
a
n
d
to
d
eter
m
i
n
e
th
e
o
p
ti
m
u
m
v
al
u
es
o
f
th
e
b
io
lo
g
ical
p
ar
a
m
eter
s
f
o
r
A
GS
p
r
o
d
u
ctio
n
o
n
SB
R
s
y
s
te
m
a
t
h
ig
h
te
m
p
er
atu
r
e
u
s
in
g
R
SM
ap
p
r
o
ac
h
.
2.
RE
S
E
ARCH
M
E
T
H
O
D
2
.
1
.
E
x
peri
m
e
nta
l Set
up
a
nd
An
a
ly
t
ica
l
M
e
t
ho
ds
E
x
p
er
i
m
e
n
ts
w
er
e
ca
r
r
ied
o
u
t
in
t
h
r
ee
d
o
u
b
le
-
w
alled
c
y
li
n
d
r
ical
co
lu
m
n
b
io
r
ea
cto
r
w
it
h
w
o
r
k
i
n
g
v
o
lu
m
e
o
f
3
L
.
D
u
r
in
g
t
h
e
s
tar
t
-
u
p
p
er
io
d
,
1
.
5
L
o
f
ac
ti
v
ated
s
lu
d
g
e
f
r
o
m
a
m
u
n
icip
al
s
e
w
a
g
e
tr
ea
t
m
en
t
p
la
n
t
w
a
s
ad
d
ed
in
to
th
e
b
io
r
ea
cto
r
as
in
o
cu
lu
m
s
.
P
r
o
g
r
a
m
m
ab
le
l
o
g
ic
co
n
tr
o
ller
(
P
L
C
)
w
er
e
ap
p
lied
to
co
n
tr
o
l
t
h
e
SB
R
o
p
er
ati
on
(
f
ee
d
in
g
p
u
m
p
,
d
is
ch
ar
g
e
p
u
m
p
an
d
ae
r
ato
r
p
u
m
p
w
it
h
th
e
s
etti
n
g
t
i
m
e
f
o
r
ea
ch
p
h
ase)
.
T
h
e
w
o
r
k
i
n
g
te
m
p
er
at
u
r
e
f
o
r
th
e
b
io
r
ea
cto
r
w
er
e
co
n
tr
o
lled
at
th
r
ee
d
if
f
er
en
t
te
m
p
er
atu
r
e
(
3
0
,
4
0
,
an
d
5
0
±
1
°C
)
u
s
i
n
g
w
a
ter
b
ath
s
leev
e
s
an
d
th
e
th
er
m
o
s
ta
t
w
ith
o
u
t
co
n
tr
o
llin
g
t
h
e
o
x
y
g
en
a
n
d
p
H
lev
el.
C
OD
an
d
p
h
o
s
p
h
ate
an
a
l
y
s
is
w
a
s
p
er
f
o
r
m
ed
as
d
escr
ib
ed
in
th
e
S
tan
d
ar
d
Me
th
o
d
s
f
o
r
th
e
E
x
a
m
in
atio
n
o
f
W
ater
an
d
W
aste
w
a
ter
(
A
P
HA
,
2
0
1
2
)
.
T
h
e
d
ata
h
as
b
ee
n
co
llected
ev
er
y
t
h
r
ee
d
a
y
s
i
n
6
0
d
a
y
s
.
T
h
e
ap
p
r
o
p
r
iate
ex
p
er
i
m
e
n
tal
an
d
a
n
al
y
tical
m
et
h
o
d
s
ca
n
b
e
r
ef
er
at
th
i
s
p
ap
er
[
5
]
T
h
e
to
tal
C
OD
an
d
P
h
o
s
p
h
o
r
u
s
(
P
)
r
e
m
o
v
al
e
f
f
icien
c
y
w
as c
alc
u
l
ated
u
s
i
n
g
E
q
u
at
io
n
(
1
)
:
C
OD
re
m
ov
a
l
(
%
)
=
CO
D
0
−
CO
D
CO
D
0
(
1
)
w
h
er
e
C
OD
₀
a
n
d
C
OD
ar
e
th
e
in
itia
l
co
n
ce
n
tr
atio
n
o
f
m
u
n
i
cip
al
s
e
w
ag
e
w
aste
w
ater
an
d
t
h
e
co
n
ce
n
tr
atio
n
o
f
p
er
m
ea
te
p
r
o
d
u
ce
d
.
2
.
2
.
Respo
ns
e
Su
rf
a
ce
M
et
ho
do
l
o
g
y
(
RS
M
)
T
h
e
te
m
p
er
atu
r
e
an
d
ti
m
e
w
er
e
s
elec
ted
as
i
n
d
ep
en
d
en
t
f
ac
t
o
r
s
an
d
th
e
to
tal
C
OD
r
e
m
o
v
a
l
an
d
to
tal
P
h
o
s
p
h
o
r
u
s
r
e
m
o
v
a
l
w
a
s
c
h
o
s
en
a
s
t
h
e
r
e
s
p
o
n
s
e
f
o
r
t
h
e
o
p
ti
m
izatio
n
m
o
d
el.
T
h
e
t
w
o
i
n
f
l
u
en
ts
ar
e
d
esi
g
n
ated
at
h
i
g
h
(
+1
)
,
m
id
d
le
(
0
)
an
d
lo
w
(
-
1
)
lev
els
ac
co
r
d
in
g
to
B
o
x
-
B
eh
n
k
en
Desi
g
n
(
B
B
D)
.
T
h
e
r
an
g
es
a
n
d
le
v
el
s
u
s
ed
i
n
th
e
e
x
p
er
i
m
e
n
t
w
er
e
d
eter
m
i
n
ed
b
y
p
r
eli
m
i
n
ar
y
e
x
p
er
i
m
en
ts
b
ased
o
n
li
ter
atu
r
e
(
s
ee
T
ab
le
1
)
,
w
h
er
e
X
1
an
d
X
2
d
en
o
ted
b
y
te
m
p
er
a
tu
r
e
an
d
ti
m
e
r
esp
ec
tiv
el
y
.
T
ab
le
1
.
L
ev
el
o
f
v
ar
iab
les o
f
B
o
x
-
B
eh
n
k
e
n
Des
ig
n
d
esi
g
n
V
a
r
i
a
b
l
e
s
R
a
n
g
e
a
n
d
l
e
v
e
l
s
-
1
0
+1
(X
1
)
T
e
mp
e
r
a
t
u
r
e
,
°
C
30
40
50
(X
2
)
T
i
me
,
d
a
y
1
30
60
Fo
r
s
tatis
t
ical
ca
lc
u
latio
n
,
th
e
p
u
r
p
o
s
e
o
f
co
n
v
er
ted
in
d
ep
en
d
en
t
v
ar
iab
les
i
n
to
d
i
m
en
s
io
n
les
s
co
d
if
ied
v
al
u
es
is
to
let
t
h
e
co
m
p
ar
i
s
o
n
o
f
f
ac
to
r
s
o
f
d
i
f
f
er
en
t
n
at
u
r
es
w
i
th
d
if
f
er
e
n
t
u
n
its
.
B
esid
es,
it
h
elp
s
to
r
ed
u
ce
th
e
er
r
o
r
in
th
e
p
o
l
y
n
o
m
ial
f
it
[
1
4
]
.
T
h
e
v
ar
iab
les
X
i
w
er
e
co
d
ed
as
x
i
ac
co
r
d
in
g
to
th
e
f
o
llo
w
i
n
g
E
q
u
atio
n
(
2
)
:
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2
0
8
8
-
8708
I
J
E
C
E
Vo
l.
7
,
No
.
3
,
J
u
n
e
2
0
1
7
:
1
5
2
2
–
1
5
2
9
1524
x
i
=
X
i
−
X
0
δ
X
(
2
)
w
h
er
e
X
i
is
t
h
e
u
n
co
d
ed
v
alu
e
o
f
th
e
i
t
h
in
d
ep
en
d
e
n
t
v
ar
ia
b
le,
X
0
th
e
v
al
u
e
o
f
X
i
at
t
h
e
ce
n
tr
e
p
o
in
t
o
f
t
h
e
in
v
e
s
ti
g
ated
ar
ea
an
d
δX
is
th
e
s
tep
ch
an
g
e.
T
h
e
ex
p
er
im
e
n
ts
w
er
e
p
er
f
o
r
m
ed
r
an
d
o
m
l
y
t
o
av
o
id
s
y
s
te
m
atic
er
r
o
r
s
[
1
5
]
.
A
s
ec
o
n
d
-
o
r
d
er
p
o
l
y
n
o
m
ial
m
o
d
el
b
y
t
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
r
eg
r
ess
io
n
p
r
o
ce
d
u
r
e
th
at
id
en
ti
f
ied
all
o
f
th
e
p
o
s
s
ib
le
i
n
ter
ac
tio
n
s
b
etw
ee
n
th
e
s
elec
te
d
f
ac
to
r
s
w
a
s
as E
q
u
atio
n
(
3
)
:
Y
=
b
0
+
∑
b
i
X
i
k
i
=
1
+
∑
b
ii
X
i
2
k
i
=
1
+
∑
∑
b
ij
X
i
X
j
j
i
<
i
(
i
=
1
,
2
,
…
,
j
)
(
3
)
w
h
er
e
Y:
p
r
ed
icted
r
esp
o
n
s
e
,
b
0
:
co
n
s
tan
t
co
ef
f
icie
n
t,
b
i
:
lin
ea
r
ef
f
ec
t
co
ef
f
icien
t,
b
ij
:
in
ter
ac
tio
n
ef
f
ec
t
co
ef
f
icie
n
t,
b
ii
: q
u
ad
r
atic
ef
f
ec
t c
o
ef
f
icie
n
t a
n
d
k
:
n
u
m
b
er
o
f
f
ac
to
r
s
[
1
1
]
.
Mo
d
el
ter
m
s
w
er
e
d
ev
elo
p
a
n
d
ev
al
u
ated
b
y
an
al
y
s
is
o
f
v
ar
ian
ce
(
A
N
OV
A
)
f
o
r
ea
ch
r
esp
o
n
s
e.
T
h
e
an
al
y
s
is
w
a
s
e
v
al
u
ated
u
s
i
n
g
s
o
f
t
w
ar
e
De
s
ig
n
E
x
p
er
t©
v
er
s
i
o
n
6
.
0
.
4
State
E
ase
i
n
c.
,
US
A
.
Sev
er
al
r
esp
o
n
s
e
s
w
er
e
ex
a
m
i
n
ed
in
t
h
is
s
t
u
d
y
a
n
d
th
e
d
if
f
er
en
t
d
eg
r
ee
p
o
ly
n
o
m
ial
m
o
d
els
w
er
e
u
s
ed
f
o
r
d
ata
f
itti
n
g
.
I
n
o
r
d
er
to
q
u
an
ti
f
y
t
h
e
c
u
r
v
atu
r
e
e
f
f
ec
ts
,
th
e
d
ata
f
r
o
m
t
h
e
e
x
p
er
im
e
n
tal
r
es
u
lts
w
er
e
f
itte
d
t
o
h
i
g
h
er
d
e
g
r
ee
p
o
ly
n
o
m
ia
l
eq
u
atio
n
s
s
u
ch
a
s
t
w
o
-
f
ac
to
r
in
ter
ac
tio
n
(
2
FI)
,
q
u
ad
r
atic
m
o
d
el,
q
u
ar
tic
m
o
d
el
etc.
[
1
2
]
,
[
1
6
]
.
T
h
en
,
th
e
h
i
g
h
er
p
o
l
y
n
o
m
ial
m
o
d
el
ter
m
s
i
n
v
o
l
v
ed
in
t
h
e
eq
u
atio
n
s
ar
e
th
o
s
e
t
h
at
ar
e
s
elec
ted
af
ter
t
h
e
eli
m
i
n
atio
n
o
f
i
n
s
i
g
n
if
ica
n
t
v
a
r
iab
les an
d
th
eir
i
n
ter
ac
tio
n
[
6
]
,
[
1
2
]
,
[
1
6
]
.
3.
RE
SU
L
T
S
A
ND
AN
AL
Y
SI
S
3
.
1
.
ANO
VA
Ana
ly
s
is
E
x
p
er
i
m
e
n
ts
w
er
e
ca
r
r
ied
o
u
t
u
s
i
n
g
B
B
D
m
et
h
o
d
to
g
e
n
er
ate
a
co
r
r
elatio
n
b
et
w
ee
n
t
w
o
d
if
f
er
e
n
t
p
r
o
ce
s
s
v
ar
iab
les
(
te
m
p
er
atu
r
e
(
X
1
)
an
d
ti
m
e
(
X
2
)
)
a
n
d
t
w
o
r
esp
o
n
s
es
(
%
o
f
C
OD
r
e
m
o
v
al
an
d
%
o
f
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al)
.
T
h
er
e
ar
e
1
3
r
u
n
s
o
f
ex
p
er
i
m
en
t
s
w
er
e
co
n
d
u
cted
an
d
th
e
r
esp
o
n
s
es
o
b
t
ain
ed
w
er
e
s
p
ec
if
ied
i
n
T
ab
le
2
.
T
h
e
r
esp
o
n
s
e
o
f
e
x
p
er
i
m
e
n
tal
r
e
s
u
l
ts
s
h
o
w
s
t
h
at
t
h
e
p
er
ce
n
ta
g
e
o
f
C
OD
r
e
m
o
v
al
w
as
i
n
th
e
r
an
g
e
o
f
4
1
%
to
9
8
.
1
7
%,
w
h
ile
th
e
p
er
ce
n
tag
e
o
f
P
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
w
as
in
b
et
wee
n
1
5
%
an
d
7
2
%.
Fro
m
T
ab
le
2
,
th
e
p
r
ed
icted
r
e
s
p
o
n
s
es
f
o
llo
w
ed
th
e
e
x
p
er
i
m
en
tal
r
esp
o
n
s
e
s
w
ell.
T
ab
le
2
.
E
x
p
er
im
e
n
tal
an
d
p
r
ed
icted
v
alu
es o
f
to
tal
C
OD
a
n
d
p
h
o
s
p
h
o
r
u
s
in
SB
R
s
y
s
te
m
R
u
n
V
a
r
i
a
b
l
e
s
R
e
sp
o
n
se
s (%)
T
e
mp
e
r
a
t
u
r
e
T
i
me
T
o
t
a
l
C
O
D
r
e
mo
v
a
l
T
o
t
a
l
P
h
o
sp
h
o
r
u
s
r
e
mo
v
a
l
X
1
(
°
C
)
X
2
(
d
a
y
)
Ex
p
t
.
P
r
e
d
.
Ex
p
t
.
P
r
e
d
.
1
50
1
4
1
.
0
0
0
3
7
.
0
8
1
5
.
0
0
0
1
9
.
9
4
2
50
15
5
6
.
0
0
0
5
7
.
3
4
4
1
.
6
1
0
3
3
.
4
5
3
50
45
9
1
.
0
0
0
8
9
.
2
3
6
5
.
0
0
0
6
2
.
4
1
1
4
50
60
9
8
.
1
7
1
9
9
.
2
8
7
2
.
4
6
0
7
6
.
8
8
5
40
3
4
7
.
5
0
0
5
3
.
6
5
3
0
.
0
0
0
3
6
.
0
4
6
40
18
7
1
.
7
9
5
7
0
.
0
1
3
8
.
7
6
2
4
5
.
8
7
7
40
39
8
9
.
7
4
4
8
6
.
3
2
6
3
.
7
9
0
5
9
.
6
4
8
40
57
9
5
.
1
2
2
9
4
.
1
8
6
6
.
2
5
0
7
1
.
4
4
9
30
1
6
2
.
5
0
0
6
0
.
5
1
5
3
.
5
0
0
4
9
.
5
2
10
30
21
7
6
.
1
5
4
7
5
.
9
1
5
5
.
7
1
4
5
6
.
4
4
11
30
45
8
3
.
4
1
5
8
5
.
1
7
7
0
.
0
0
0
6
4
.
7
4
12
30
60
8
5
.
3
6
6
8
5
.
8
5
6
8
.
5
0
0
6
9
.
9
3
13
50
30
7
2
.
0
0
0
7
5
.
2
5
5
3
.
6
4
0
4
7
.
9
3
T
h
e
r
esu
lt
s
w
er
e
a
n
al
y
ze
d
u
s
i
n
g
An
al
y
s
is
o
f
Var
ian
ce
(
A
N
OV
A
)
an
d
s
u
m
m
ar
ized
i
n
T
ab
le
3
.
T
h
e
an
al
y
s
is
s
h
o
w
s
t
h
e
q
u
ad
r
atic
r
eg
r
ess
io
n
m
o
d
el
f
o
r
to
tal
r
e
m
o
v
al
o
f
C
OD,
w
h
ile
to
tal
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
w
a
s
t
w
o
-
f
ac
to
r
in
ter
ac
tio
n
r
e
g
r
es
s
i
o
n
m
o
d
el
(
E
q
u
atio
n
s
4
an
d
5
)
.
T
h
is
t
w
o
r
e
g
r
ess
io
n
m
o
d
el
s
i
n
d
icate
d
th
e
m
o
d
el
ar
e
h
ig
h
l
y
s
i
g
n
i
f
ica
n
t
d
u
e
to
th
e
P
v
al
u
e
(
P
r
o
b
>F)
w
as
v
e
r
y
lo
w
(
<0
.
0
0
0
1
)
.
I
t
is
b
ec
au
s
e,
in
o
r
d
er
f
o
r
th
e
m
o
d
el
to
b
e
s
i
g
n
if
ican
t
at
9
5
%
co
n
f
id
en
ce
lev
el,
th
e
“
P
r
o
b
>F”
m
u
s
t
les
s
t
h
an
0
.
0
5
[
6
]
.
T
h
e
s
m
aller
th
e
m
ag
n
it
u
d
e
o
f
t
h
e
P
v
al
u
e,
t
h
e
m
o
r
e
s
ig
n
i
f
ica
n
t
w
a
s
t
h
e
co
r
r
esp
o
n
d
in
g
co
ef
f
icie
n
t.
T
h
er
e
f
o
r
e,
it
ca
n
b
e
s
ee
n
th
at
th
e
s
ig
n
i
f
ica
n
t
ter
m
s
f
o
r
C
OD
r
e
m
o
v
a
l
ar
e
X
2
,
X
2
2
an
d
X
1
X
2
,
w
h
er
ea
s
X
1
an
d
X
1
2
ar
e
tr
iv
ial
f
o
r
C
OD
r
e
m
o
v
al.
W
h
ile,
all
ter
m
s
u
s
e
d
f
o
r
th
e
to
tal
p
h
o
s
p
h
o
r
u
s
(
X
1
,X
2
an
d
X
1
X
2
)
ar
e
s
i
g
n
i
f
ica
n
t.
T
h
is
P
v
a
lu
e
s
w
er
e
u
s
ed
as
a
to
o
l
to
co
n
f
ir
m
th
at
th
e
q
u
ad
r
atic
o
r
t
w
o
-
f
ac
to
r
i
n
ter
ac
tio
n
m
o
d
el
w
a
s
s
tati
s
tic
all
y
s
i
g
n
i
f
ican
t
an
d
s
u
f
f
icie
n
t
to
e
x
p
lain
th
e
ac
tu
al
r
elatio
n
s
h
ip
b
et
w
ee
n
t
h
e
r
esp
o
n
s
e
a
n
d
th
e
s
i
g
n
if
ican
t
v
ar
ia
b
les
[
1
7
]
.
Ho
w
e
v
er
,
m
o
d
el
F
-
v
al
u
e
w
as
ca
lcu
la
ted
as
a
r
atio
o
f
m
ea
n
s
q
u
ar
e
r
eg
r
ess
io
n
a
n
d
m
ea
n
s
q
u
ar
e
r
es
id
u
al.
T
h
e
F
-
v
alu
e
an
d
th
e
co
r
r
esp
o
n
d
in
g
P
v
al
u
e
s
w
er
e
g
i
v
en
i
n
T
ab
le
3
.
Evaluation Warning : The document was created with Spire.PDF for Python.
I
J
E
C
E
I
SS
N:
2
0
8
8
-
8708
P
a
r
a
mete
r
Op
timis
a
tio
n
o
f A
ero
b
ic
Gra
n
u
la
r
S
lu
d
g
e
a
t H
ig
h
Temp
era
tu
r
e
Usi
n
g
.
.
.
.
(
S
ya
h
ir
a
I
b
r
a
h
im)
1525
T
h
e
s
q
u
ar
e
o
f
co
r
r
elatio
n
co
ef
f
icie
n
t
f
o
r
ea
ch
r
esp
o
n
s
e
w
a
s
co
m
p
u
te
as
a
co
ef
f
icie
n
t
o
f
d
eter
m
in
at
io
n
,
R
2
a
n
d
t
h
e
v
al
u
e
i
s
al
w
a
y
s
b
et
w
ee
n
0
an
d
1
[
1
2
]
,
T
h
e
v
al
u
e
o
f
R
2
>
0
.
9
0
in
d
icate
s
t
h
at
t
h
e
m
o
d
el
is
e
x
ce
lle
n
t.
A
cc
o
r
d
in
g
to
S.
Sat
h
ian
et
al.
[
4
]
,
a
g
o
o
d
s
tatis
tical
m
o
d
el
s
h
o
u
ld
h
av
e
R
2
v
a
lu
e
clo
s
es
to
1
.
0
.
I
n
th
i
s
ca
s
e,
t
h
e
p
r
ed
icted
R
2
v
al
u
e
f
o
r
to
tal
r
e
m
o
v
al
o
f
C
OD
a
n
d
p
h
o
s
p
h
o
r
u
s
r
esp
o
n
s
es
w
as
i
n
r
ea
s
o
n
ab
le
ag
r
ee
m
en
t
s
i
n
ce
th
e
f
it o
f
th
e
m
o
d
el
w
as
f
o
u
n
d
t
o
b
e
0
.
9
7
7
an
d
0
.
9
1
r
esp
ec
tiv
el
y
.
I
t c
o
n
cl
u
d
es t
h
at
th
e
p
r
ed
ictio
n
m
o
d
el
o
b
tain
e
d
u
s
i
n
g
B
B
D
tech
n
iq
u
e
o
f
R
SM
ap
p
r
o
ac
h
ca
n
r
elativ
el
y
r
ef
lect
t
h
e
ef
f
ec
t
s
o
f
te
m
p
er
atu
r
e
a
n
d
ti
m
e
o
n
t
h
e
p
er
f
o
r
m
a
n
ce
o
f
SB
R
s
y
s
te
m
o
n
t
h
e
e
x
a
m
in
a
tio
n
o
f
s
y
n
th
e
tic
w
as
te
w
ater
tr
ea
t
m
e
n
t b
ased
o
n
n
u
tr
ien
t r
e
m
o
v
al.
T
ab
le
3
.
A
NOV
A
an
al
y
s
is
f
o
r
to
tal
C
OD
an
d
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
S
o
u
r
c
e
T
o
t
a
l
C
O
D
r
e
mo
v
a
l
T
o
t
a
l
P
r
e
mo
v
a
l
.
F
V
a
l
u
e
P
r
o
b
>
F
F
V
a
l
u
e
P
r
o
b
>
F
M
o
d
e
l
5
9
.
4
7
<
0
.
0
0
0
1
3
0
.
5
2
<
0
.
0
0
0
1
X
1
-
t
e
mp
e
r
a
t
u
r
e
4
.
1
5
0
.
0
8
1
1
7
.
8
4
0
.
0
2
0
7
X
2
-
t
i
me
2
4
4
.
4
6
<
0
.
0
0
0
1
6
9
.
8
5
<
0
.
0
0
0
1
X
1
2
1
.
1
6
0
.
3
1
6
8
X
2
2
8
.
9
7
0
.
0
2
0
1
X
1
X
2
3
1
.
0
4
0
.
0
0
0
8
1
1
.
1
4
0
.
0
0
8
7
S
u
m o
f
s
q
u
a
r
e
s
3
9
4
8
.
4
6
3
3
2
8
.
1
4
M
e
a
n
sq
u
a
r
e
s
7
8
9
.
6
9
1
1
0
9
.
3
8
S
t
d
.
D
e
v
.
3
.
6
4
6
.
0
3
M
e
a
n
7
4
.
6
0
5
3
.
4
0
R
2
0
.
9
7
7
0
0
.
9
1
0
5
I
n
o
r
d
er
to
in
v
esti
g
ate
t
h
e
e
f
f
ec
t
o
f
h
i
g
h
te
m
p
er
at
u
r
e
o
n
t
h
e
f
o
r
m
atio
n
o
f
AGS
i
n
SB
R
s
y
s
te
m
a
n
d
o
v
er
all
p
er
f
o
r
m
a
n
ce
o
f
b
io
lo
g
ical
tr
ea
t
m
e
n
t,
to
tal
C
OD
a
n
d
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
e
f
f
icie
n
c
y
w
as
ca
lcu
lated
a
s
a
r
esp
o
n
s
e.
F
ig
u
r
e
1
(
a)
an
d
(
b
)
d
e
m
o
n
s
tr
ate
s
a
v
er
y
g
o
o
d
co
n
v
er
s
io
n
b
et
w
ee
n
t
h
e
e
x
p
er
i
m
en
tal
a
n
d
p
r
ed
icted
v
alu
e
s
o
f
t
w
o
r
esp
o
n
s
es
(
C
O
D
an
d
p
h
o
s
p
h
o
r
u
s
)
,
w
it
h
r
e
m
o
v
al
ef
f
icie
n
c
y
0
.
9
7
7
an
d
0
.
9
1
r
esp
ec
tiv
el
y
.
A
ct
u
al
v
alu
e
s
ar
e
th
e
m
ea
s
u
r
ed
r
esp
o
n
s
e
d
ata
f
o
r
a
p
ar
ticu
lar
r
u
n
a
n
d
th
e
p
r
ed
icted
v
al
u
es
w
er
e
ev
alu
a
ted
f
r
o
m
t
h
e
m
o
d
el
an
d
g
en
er
ated
b
y
u
s
i
n
g
th
e
ap
p
r
o
x
i
m
ati
n
g
f
u
n
ctio
n
.
I
t
s
h
o
w
s
t
h
at
t
h
e
p
r
ed
icted
v
al
u
es
o
f
t
h
e
r
esp
o
n
s
e
f
r
o
m
t
h
e
m
o
d
el
a
g
r
ee
d
w
it
h
t
h
e
o
b
s
er
v
ed
v
al
u
es
a
n
d
co
u
ld
b
e
p
r
o
p
er
ly
ap
p
lied
to
n
av
i
g
a
te
th
e
d
esi
g
n
s
p
ac
e
[
1
8
]
.
B
y
ap
p
l
y
in
g
m
u
l
tip
le
r
eg
r
ess
io
n
an
al
y
s
is
o
n
t
h
e
ex
p
er
im
e
n
tal
d
ata,
b
o
th
r
esp
o
n
s
es
w
er
e
f
itted
a
n
d
p
r
o
d
u
ce
q
u
ad
r
atic
an
d
tw
o
-
f
a
cto
r
in
ter
ac
tio
n
m
o
d
els.
T
h
e
e
m
p
ir
ical
m
o
d
els
in
ter
m
s
o
f
c
o
d
ed
f
ac
to
r
s
f
o
r
th
e
q
u
ad
r
a
tic
r
eg
r
ess
io
n
eq
u
a
tio
n
o
f
C
OD
r
e
m
o
v
al
w
er
e
s
tat
ed
as
E
q
u
atio
n
(
4
)
,
w
h
ile
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
r
esp
o
n
s
e
w
as p
r
esen
t a
s
t
w
o
-
f
ac
to
r
in
ter
ac
tio
n
m
o
d
el
(
E
q
u
at
io
n
(
5
)
)
:
To
ta
l
C
OD
r
e
m
ova
l
(
%
)
=
80
.
65
−
2
.
50
X
1
+
21
.
88
X
2
−
2
.
37
X
1
2
−
7
.
59
X
2
2
+
9
.
22
X
1
X
2
(
4
)
To
ta
l
phos
phor
us
r
e
m
ova
l
(
%
)
=
54
.
07
−
5
.
65
X
1
+
19
.
34
X
2
+
9
.
13
X
1
X
2
(
5
)
w
h
er
e
X
1
an
d
X
2
w
er
e
th
e
co
d
ed
v
alu
es
o
f
th
e
tes
t
v
ar
iab
l
es
te
m
p
er
atu
r
e
(
°C
)
an
d
ti
m
e
(
d
ay
)
,
r
esp
ec
tiv
el
y
.
T
h
e
p
o
s
itiv
e
s
i
g
n
a
n
d
n
eg
ati
v
e
s
ig
n
r
ep
r
esen
t
s
th
e
s
y
n
er
g
is
t
ic
ef
f
ec
t
an
d
an
tag
o
n
is
tic
e
f
f
e
cts
b
et
w
ee
n
m
u
t
u
al
in
ter
ac
tio
n
a
n
d
in
d
i
v
id
u
al
p
ar
a
m
eter
s
r
esp
ec
ti
v
el
y
[
1
7
]
.
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2
0
8
8
-
8708
I
J
E
C
E
Vo
l.
7
,
No
.
3
,
J
u
n
e
2
0
1
7
:
1
5
2
2
–
1
5
2
9
1526
Fig
u
r
e
1
.
P
r
ed
icte
d
v
s
.
ac
tu
al
v
alu
e
s
p
lo
t f
o
r
C
OD
r
e
m
o
v
al
(
a)
an
d
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
(
b
)
.
3
.
2
.
E
f
f
ec
t
o
f
v
a
ria
bles
o
n t
re
a
t
m
ent
o
f
s
y
nthet
ic
wa
s
t
ewa
t
er
i
n SB
R
s
y
s
t
e
m
T
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
a
n
d
c
o
n
to
u
r
p
lo
ts
w
er
e
g
e
n
er
ated
f
o
r
d
if
f
er
e
n
t
i
n
ter
ac
tio
n
s
o
f
a
n
y
t
w
o
in
d
ep
en
d
en
t
v
ar
iab
les,
w
h
ile
h
o
ld
in
g
th
e
v
al
u
e
o
f
t
h
e
o
th
er
v
ar
iab
les
as
a
co
n
s
ta
n
t
[
4
]
.
T
h
e
o
p
ti
m
izatio
n
o
f
p
r
o
ce
s
s
v
ar
iab
les
w
as
ai
m
ed
at
f
i
n
d
in
g
t
h
e
lev
el
s
o
f
i
n
d
e
p
en
d
en
t
v
ar
iab
le
s
,
w
h
ic
h
w
o
u
ld
g
i
v
e
m
ax
i
m
u
m
p
er
ce
n
tag
e
o
f
to
tal
r
e
m
o
v
al
f
o
r
C
OD
an
d
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
a
l.
T
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
p
lo
ts
an
d
th
e
co
r
r
esp
o
n
d
in
g
co
n
to
u
r
c
u
r
v
es
d
escr
ib
ed
b
y
th
e
r
eg
r
e
s
s
io
n
e
q
u
atio
n
ar
e
s
h
o
w
n
i
n
Fi
g
u
r
e
2
an
d
3
,
(
a)
an
d
(
b
)
.
Fig
u
r
e
2
(
a)
s
h
o
w
s
t
h
e
i
n
t
er
ac
tio
n
b
et
w
ee
n
ti
m
e
a
n
d
te
m
p
er
atu
r
e
o
n
p
er
ce
n
ta
g
e
o
f
to
tal
C
OD.
T
h
e
o
b
s
er
v
atio
n
r
e
v
ea
ls
th
at
t
h
e
p
er
ce
n
tag
e
o
f
to
tal
C
OD
r
em
o
v
al
in
cr
ea
s
ed
w
h
e
n
th
e
ti
m
e
w
a
s
i
n
cr
ea
s
e
a
t
te
m
p
er
atu
r
e
in
cr
ea
s
e.
T
h
e
m
a
x
i
m
u
m
to
tal
C
OD
r
e
m
o
v
a
l
(
8
6
.
3
3
%)
w
as
o
b
s
er
v
ed
at
te
m
p
er
atu
r
e
(
4
0
°C
)
a
n
d
o
p
er
atio
n
al
ti
m
e
o
f
3
9
d
a
y
s
.
Af
ter
t
h
at,
u
n
d
er
co
n
s
ta
n
t
o
p
er
atio
n
al
ti
m
e,
a
s
l
ig
h
t
i
n
cr
ea
s
e
in
to
tal
C
OD
re
m
o
v
al
w
as
o
b
s
er
v
ed
f
o
r
h
ig
h
er
te
m
p
er
atu
r
e
(
>4
0
°C
)
.
T
h
is
h
as
b
ee
n
s
u
p
p
o
r
ted
b
y
M.
H.
A
b
Hali
m
et
al.
[
5
]
,
w
h
er
e
h
i
g
h
te
m
p
er
at
u
r
e
w
as
l
ea
d
f
o
r
th
e
h
i
g
h
as
h
co
n
te
n
t
i
n
b
io
g
r
an
u
les
a
n
d
ac
co
r
d
in
g
l
y
in
cr
ea
s
e
th
e
s
ize
.
T
h
e
r
esu
lts
ca
n
b
e
m
o
r
e
clea
r
l
y
r
ev
ea
led
b
y
r
e
f
er
r
in
g
to
T
ab
le
3
,
w
h
er
e
th
e
F
-
v
al
u
e
o
f
ti
m
e
(
2
4
4
.
4
6
)
w
a
s
h
ig
h
er
t
h
an
t
h
e
te
m
p
er
atu
r
e
(
4
.
1
5
)
,
s
u
g
g
e
s
ti
n
g
t
h
at
t
h
e
o
p
er
atio
n
al
ti
m
e
h
ad
m
o
r
e
i
m
p
ac
t
o
n
r
e
m
o
v
in
g
C
OD.
B
esid
es
th
a
t,
th
e
s
i
g
n
i
f
ican
t
a
s
P
-
v
al
u
e
le
s
s
th
a
n
0
.
0
5
w
er
e
f
o
u
n
d
i
n
li
n
ea
r
ter
m
o
f
ti
m
e
(
X
2
)
,
s
q
u
ar
e
ter
m
(X
2
2
)
o
f
ti
m
e
a
n
d
q
u
ad
r
atic
ter
m
(
X
1
X
2
)
o
f
te
m
p
er
at
u
r
e
an
d
ti
m
e.
T
h
ese
r
esu
lts
d
e
m
o
n
s
tr
ate
th
at
te
m
p
er
at
u
r
e
an
d
ti
m
e
w
er
e
r
eq
u
ir
ed
b
y
C
O
D
r
e
m
o
v
al
to
m
a
x
i
m
ize
A
G
S
f
o
r
m
atio
n
in
SB
R
s
y
s
te
m
.
T
h
er
e
is
a
s
i
m
ilar
tr
en
d
o
f
3
D
g
r
ap
h
f
o
r
t
h
e
to
tal
C
OD
a
n
d
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
(
Fi
g
u
r
e
2
(
a)
an
d
Fig
u
r
e
3
(
a)
)
,
w
h
er
e
t
h
e
p
er
ce
n
tag
e
o
f
to
tal
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
in
cr
ea
s
ed
w
h
e
n
th
e
ti
m
e
w
as
in
cr
ea
s
e
at
te
m
p
er
atu
r
e
i
n
cr
ea
s
e.
T
h
e
n
,
t
h
e
o
p
er
atio
n
al
t
i
m
e
w
as
h
i
g
h
l
y
s
ig
n
i
f
ica
n
t
o
n
t
h
e
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
a
l
co
m
p
ar
e
to
ti
m
e.
Ho
w
e
v
er
,
th
e
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
w
as
i
n
cr
ea
s
e
s
g
r
ad
u
all
y
w
i
th
t
i
m
e
w
h
e
n
t
h
e
te
m
p
er
at
u
r
e
w
a
s
in
cr
ea
s
e.
T
h
e
s
a
m
e
tr
e
n
d
w
er
e
ar
ch
iv
ed
f
r
o
m
th
e
p
r
ev
io
u
s
s
tu
d
y
r
ep
o
r
ted
b
y
M.
H.
A
b
H
ali
m
e
t
al.
[
5
]
,
w
h
en
th
e
a
v
er
ag
e
r
e
m
o
v
al
r
ate
f
o
r
5
0
°C
is
t
h
e
h
i
g
h
e
s
t,
f
o
llo
w
ed
b
y
3
0
°C
a
n
d
4
0
°C
.
I
n
ad
d
itio
n
,
t
h
e
p
er
ce
n
ta
g
e
o
f
C
OD
r
e
m
o
v
a
l
w
a
s
s
ee
n
to
b
e
h
ig
h
er
th
a
n
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al.
T
h
ese
w
as
d
u
e
to
ac
etate
(
C
OD)
w
a
s
u
til
ized
to
3
4
.
5
%
at
3
0
°C
,
5
1
.
2
2
%
at
4
0
°C
an
d
3
8
.
7
9
%
at
5
0
°C
,
af
ter
t
h
at,
th
e
r
e
s
t
o
f
ac
eta
te
w
a
s
u
t
ilized
i
n
th
e
p
er
io
d
o
f
ae
r
o
b
ic
cy
cle
w
h
il
e
p
h
o
s
p
h
ate
w
as
r
e
m
o
v
ed
.
As
ill
u
s
tr
ated
in
A
NO
V
A
a
n
al
y
s
i
s
,
t
h
e
m
u
tu
a
l
in
ter
ac
tio
n
b
et
w
ee
n
b
o
th
p
ar
am
eter
w
as
h
ig
h
l
y
s
i
g
n
i
f
ica
n
t,
it
w
a
s
o
b
s
er
v
ed
b
y
a
v
er
y
s
m
all
“
P
r
o
b
>F”
v
alu
e
(
0
.
0
0
8
7
)
.
T
h
u
s
,
it
m
a
y
co
n
clu
d
e
th
at
t
h
e
i
n
ter
ac
tio
n
o
f
te
m
p
er
atu
r
e
an
d
ti
m
e
w
er
e
r
eq
u
ir
ed
in
r
e
m
o
v
in
g
t
h
e
to
tal
C
O
D
a
n
d
p
h
o
s
p
h
o
r
u
s
i
n
SB
R
s
y
s
te
m
.
B
esid
es
t
h
at,
it
s
h
o
w
s
th
at
te
m
p
er
at
u
r
e
is
n
o
t
a
cr
itical
f
ac
to
r
i
n
th
e
f
o
r
m
a
tio
n
o
f
A
GS
i
n
SB
R
,
h
o
w
e
v
er
it
w
as
ab
le
to
in
d
u
c
e
th
e
m
o
r
p
h
o
lo
g
y
a
n
d
s
ettlin
g
ab
ilit
y
o
f
g
r
an
u
les
[
1
9
]
.
A
lth
o
u
g
h
,
th
is
ex
p
er
i
m
en
tal
r
esu
lts
w
as
d
i
f
f
er
w
it
h
th
e
s
tate
m
e
n
t
m
ad
e
b
y
So
n
g
Z
h
i
w
ei
et
al.
[
1
9
]
,
w
h
ic
h
d
ec
lar
ed
th
at
th
e
b
io
ac
tiv
it
y
o
f
m
icr
o
b
es
w
il
l
b
ec
o
m
e
w
o
r
s
e
an
d
co
u
ld
also
in
f
l
u
e
n
ce
th
e
b
io
lo
g
ical
r
ea
ctio
n
ef
f
icie
n
tl
y
w
h
e
n
t
h
e
t
e
m
p
er
atu
r
e
ex
ce
ed
3
0
°C
.
D
ESI
G
N
-EXPER
T
Pl
o
t
T
o
t
a
l
P
re
mo
va
l
A
c
t
u
a
l
P
r
e
d
i
c
t
e
d
P
re
d
i
ct
e
d
v
s
.
A
ct
u
a
l
1
5
.0
0
3
0
.4
7
4
5
.9
4
6
1
.4
1
7
6
.8
8
1
5
.0
0
3
0
.4
7
4
5
.9
4
6
1
.4
1
7
6
.8
8
D
ESI
G
N
-EXPER
T
Pl
o
t
T
o
t
a
l
C
O
D
re
mo
va
l
A
c
t
u
a
l
P
r
e
d
i
c
t
e
d
P
re
d
i
ct
e
d
v
s
.
A
ct
u
a
l
3
7
.0
8
5
2
.6
3
6
8
.1
8
8
3
.7
3
9
9
.2
8
3
7
.0
8
5
2
.6
3
6
8
.1
8
8
3
.7
3
9
9
.2
8
(
a)
(
b
)
Evaluation Warning : The document was created with Spire.PDF for Python.
I
J
E
C
E
I
SS
N:
2
0
8
8
-
8708
P
a
r
a
mete
r
Op
timis
a
tio
n
o
f A
ero
b
ic
Gra
n
u
la
r
S
lu
d
g
e
a
t H
ig
h
Temp
era
tu
r
e
Usi
n
g
.
.
.
.
(
S
ya
h
ir
a
I
b
r
a
h
im)
1527
(
a)
(
b
)
Fig
u
r
e
2
.
Gr
ap
h
d
ep
ictin
g
t
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
p
lo
t (
a)
an
d
co
n
to
u
r
p
lo
t (
b
)
s
h
o
w
i
n
g
t
h
e
e
f
f
ec
t o
f
te
m
p
er
atu
r
e
(X
1
)
an
d
ti
m
e
(
X
2
)
an
d
th
eir
m
u
tu
a
l in
ter
ac
tio
n
o
f
to
tal
C
OD
r
e
m
o
v
al
i
n
SB
R
s
y
s
te
m
(
a)
(
b
)
Fig
u
r
e
3
.
Gr
ap
h
d
ep
ictin
g
t
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
p
lo
t (
a)
an
d
co
n
to
u
r
p
lo
t (
b
)
s
h
o
w
i
n
g
t
h
e
e
f
f
ec
t o
f
te
m
p
er
atu
r
e
(X
1
)
an
d
ti
m
e
(
X
2
)
an
d
th
eir
m
u
tu
a
l in
ter
ac
tio
n
o
f
to
tal
p
h
o
s
p
h
o
r
u
s
r
e
m
o
v
al
i
n
SB
R
s
y
s
te
m
T
h
e
ex
p
er
i
m
en
tal
co
n
d
itio
n
w
er
e
p
r
o
p
o
s
ed
to
f
in
d
t
h
e
o
p
ti
m
u
m
p
o
i
n
t
i
n
o
r
d
er
to
o
b
tain
i
n
g
th
e
h
ig
h
e
s
t
co
n
v
er
s
io
n
o
n
m
a
x
i
m
i
ze
th
e
to
tal
r
em
o
v
al
co
n
d
itio
n
.
T
h
e
r
esu
lt
o
f
th
e
p
r
o
p
o
s
ed
co
n
d
itio
n
is
s
h
o
w
n
at
T
ab
le
4
.
T
h
e
o
p
tim
ize
co
n
d
i
tio
n
v
al
u
es
s
u
g
g
e
s
ted
b
y
R
S
M
ap
p
r
o
ac
h
w
er
e
4
5
.
3
3
°C
a
t
6
0
d
ay
s
,
w
h
ich
co
r
r
esp
o
n
d
ed
to
th
e
ac
tu
al
e
x
p
er
i
m
e
n
tal
v
a
lu
e
s
o
f
9
7
.
8
4
%
an
d
7
5
.
2
6
%
f
o
r
to
tal
r
e
m
o
v
al
o
f
C
OD
a
n
d
p
h
o
s
p
h
o
r
u
s
r
esp
ec
ti
v
el
y
.
T
h
e
ex
p
er
i
m
e
n
tal
v
a
lu
e
s
w
er
e
r
ela
tiv
el
y
c
lo
s
e
to
th
e
p
r
ed
icted
v
alu
es
b
y
t
h
e
m
o
d
el
an
d
w
er
e
t
h
er
ef
o
r
e
co
n
f
ir
m
i
n
g
th
e
ad
eq
u
ac
y
a
n
d
v
alid
it
y
o
f
t
h
e
p
r
ed
icted
m
o
d
el.
T
ab
le
4
.
Op
tim
ized
v
al
u
e
o
f
t
h
e
p
r
o
ce
s
s
S
o
l
u
t
i
o
n
T
e
mp
e
r
a
t
u
r
e
(
°
C
)
T
i
me
(
d
a
y
)
T
o
t
a
l
C
O
D
r
e
mo
v
a
l
(
%)
T
o
t
a
l
P
r
e
mo
v
a
l
(
%)
D
e
si
r
a
b
i
l
i
t
y
1
4
5
.
3
3
6
0
.
0
0
9
7
.
8
4
2
2
7
5
.
2
5
7
8
0
.
9
9
8
D
E
S
I
G
N
-
E
X
P
E
R
T
P
l
o
t
T
o
t
a
l
C
O
D
r
e
m
o
v
a
l
X
=
A
:
T
e
m
p
e
r
a
t
u
r
e
Y
=
B
:
T
i
m
e
3
7
.
0
8
5
2
.
6
3
6
8
.
1
8
8
3
.
7
3
9
9
.
2
8
T
o
t
a
l
C
O
D
r
e
m
o
v
a
l
3
0
.
0
0
3
5
.
0
0
4
0
.
0
0
4
5
.
0
0
5
0
.
0
0
1
.
0
0
1
5
.
7
5
3
0
.
5
0
4
5
.
2
5
6
0
.
0
0
A
:
T
e
m
p
e
r
a
t
u
r
e
B
:
T
i
m
e
D
ESI
G
N
-EXPER
T
Pl
o
t
T
o
t
a
l
C
O
D
re
mo
va
l
D
e
si
g
n
Po
i
n
t
s
T
o
t
a
l
C
O
D
re
mo
va
l
=
8
6
.
3
3
R
u
n
#
7
X:
A:
T
e
mp
e
ra
t
u
re
=
4
0
.
0
0
Y:
B:
T
i
me
=
3
9
.
0
0
T
o
t
a
l
C
O
D
re
m
o
v
a
l
A
:
T
e
m
p
e
r
a
t
u
r
e
B
:
T
i
m
e
3
0
.
0
0
3
5
.
0
0
4
0
.
0
0
4
5
.
0
0
5
0
.
0
0
1
.
0
0
1
5
.
7
5
3
0
.
5
0
4
5
.
2
5
6
0
.
0
0
4
7
.
4
4
5
6
5
7
.
8
1
2
9
6
8
.
1
8
0
3
7
8
.
5
4
7
6
8
8
.
9
1
4
9
8
6
.
2
9
1
8
9
4
.
4
8
6
9
D
ESI
G
N
-EXPER
T
Pl
o
t
T
o
t
a
l
C
O
D
re
mo
va
l
D
e
si
g
n
Po
i
n
t
s
T
o
t
a
l
C
O
D
re
mo
va
l
=
8
6
.
3
3
R
u
n
#
7
X:
A:
T
e
mp
e
ra
t
u
re
=
4
0
.
0
0
Y:
B:
T
i
me
=
3
9
.
0
0
T
o
t
a
l
C
O
D
re
m
o
v
a
l
A
:
T
e
m
p
e
r
a
t
u
r
e
B
:
T
i
m
e
3
0
.
0
0
3
5
.
0
0
4
0
.
0
0
4
5
.
0
0
5
0
.
0
0
1
.
0
0
1
5
.
7
5
3
0
.
5
0
4
5
.
2
5
6
0
.
0
0
4
7
.
4
4
5
6
5
7
.
8
1
2
9
6
8
.
1
8
0
3
7
8
.
5
4
7
6
8
8
.
9
1
4
9
8
6
.
2
9
1
8
9
4
.
4
8
6
9
D
ESI
G
N
-EXPER
T
Pl
o
t
T
o
t
a
l
P
re
mo
va
l
D
e
si
g
n
Po
i
n
t
s
T
o
t
a
l
P
re
mo
va
l
=
7
1
.
4
4
R
u
n
#
8
X:
A:
T
e
mp
e
ra
t
u
re
=
4
0
.
0
0
Y:
B:
T
i
me
=
5
7
.
0
0
T
o
t
a
l
P
re
m
o
v
a
l
A
:
T
e
m
p
e
r
a
t
u
r
e
B
:
T
i
m
e
3
0
.0
0
3
5
.0
0
4
0
.0
0
4
5
.0
0
5
0
.0
0
1
.0
0
1
5
.7
5
3
0
.5
0
4
5
.2
5
6
0
.0
0
2
9
.
4
3
3
8
.
9
2
4
8
.
4
1
5
7
.
9
0
6
7
.
3
9
7
1
.
4
4
6
2
.
4
9
D
E
S
I
G
N
-
E
X
P
E
R
T
P
l
o
t
T
o
t
a
l
P
r
e
m
o
v
a
l
X
=
A
:
T
e
m
p
e
r
a
t
u
r
e
Y
=
B
:
T
i
m
e
1
9
.
9
4
3
4
.
1
8
4
8
.
4
1
6
2
.
6
5
7
6
.
8
8
T
o
t
a
l
P
r
e
m
o
v
a
l
3
0
.
0
0
3
5
.
0
0
4
0
.
0
0
4
5
.
0
0
5
0
.
0
0
1
.
0
0
1
5
.
7
5
3
0
.
5
0
4
5
.
2
5
6
0
.
0
0
A
:
T
e
m
p
e
r
a
t
u
r
e
B
:
T
i
m
e
Evaluation Warning : The document was created with Spire.PDF for Python.
I
SS
N
:
2
0
8
8
-
8708
I
J
E
C
E
Vo
l.
7
,
No
.
3
,
J
u
n
e
2
0
1
7
:
1
5
2
2
–
1
5
2
9
1528
4.
CO
NCLU
SI
O
N
T
h
e
r
esp
o
n
s
e
s
u
r
f
ac
e
m
et
h
o
d
o
lo
g
y
f
o
r
p
ar
a
m
eter
o
p
ti
m
i
s
at
io
n
o
f
ae
r
o
b
ic
g
r
an
u
lar
s
l
u
d
g
e
SB
R
h
a
s
b
ee
n
p
r
esen
ted
.
T
h
e
s
i
m
u
lat
io
n
r
es
u
lt
s
s
h
o
w
ed
t
h
at
R
SM
ca
n
p
r
ed
ict
t
h
e
ae
r
o
b
ic
g
r
a
n
u
lar
s
l
u
d
g
e
p
ar
a
m
eter
s
w
ell
at
d
i
f
f
er
en
t
te
m
p
er
at
u
r
e
co
n
d
itio
n
(
3
0
,
4
0
an
d
5
0
0
C
)
.
T
h
e
r
esu
lt
also
s
h
o
w
ed
t
h
at
th
e
h
i
g
h
te
m
p
er
atu
r
e
in
f
lu
e
n
ce
d
th
e
m
o
r
p
h
o
lo
g
y
a
n
d
s
ettli
n
g
ab
ilit
y
o
f
g
r
an
u
les
,
alth
o
u
g
h
less
s
i
g
n
i
f
ica
n
t
i
n
b
io
lo
g
ical
r
e
m
o
v
al.
B
esid
es
th
at,
w
it
h
o
n
l
y
a
f
e
w
n
u
m
b
er
o
f
e
x
p
er
i
m
e
n
ts
,
R
SM
ca
n
d
ev
elo
p
ed
an
ex
ce
llen
t
C
OD
r
e
m
o
v
al
ef
f
icien
c
y
.
T
h
i
s
s
tu
d
y
d
e
m
o
n
s
tr
ates
t
h
at
R
SM
ap
p
r
o
ac
h
is
s
u
itab
le
f
o
r
p
r
ed
ictin
g
t
h
e
h
i
g
h
co
n
d
itio
n
o
f
A
G
S
in
SB
R
s
y
s
te
m
b
ased
o
n
t
h
e
p
er
ce
n
tag
e
o
f
to
tal
r
e
m
o
v
al
o
f
C
OD
an
d
p
h
o
s
p
h
o
r
u
s
.
T
h
er
ef
o
r
e,
t
h
e
m
o
d
el
d
ev
elo
p
ed
in
th
is
w
o
r
k
ca
n
b
e
u
s
e
f
u
l f
o
r
p
r
ed
ic
tin
g
t
h
e
f
o
r
m
a
tio
n
o
f
A
GS i
n
SB
R
s
y
s
te
m
at
h
o
t c
li
m
ate.
ACK
NO
WL
E
D
G
E
M
E
NT
S
T
h
e
au
th
o
r
s
w
o
u
ld
li
k
e
to
t
h
an
k
t
h
e
R
esear
ch
Un
iv
er
s
it
y
Gr
an
t
(
GUP
)
v
o
te
1
3
H7
0
,
Un
i
v
er
s
iti
T
ek
n
o
lo
g
i M
ala
y
s
ia
f
o
r
t
h
e
f
i
n
an
cial
s
u
p
p
o
r
t.
RE
F
E
R
E
NC
E
S
[1
]
H.
Ha
sn
id
a
a
n
d
N.
A
z
n
a
h
,
“
D
e
v
e
lo
p
m
e
n
t
a
n
d
Util
iza
ti
o
n
o
f
A
e
ro
b
ic
G
ra
n
u
les
f
o
r
S
o
y
S
a
u
c
e
W
a
ste
wa
ter
T
re
a
t
m
e
n
t :
Op
ti
m
iza
ti
o
n
b
y
Re
sp
o
n
se
S
u
rf
a
c
e
M
e
th
o
d
o
l
o
g
y
,
”
J
.
T
e
k
n
o
l.
,
v
o
l
/i
ss
u
e
:
69
(
5
)
,
p
p
.
3
1
–
3
7
,
2
0
1
4
.
[2
]
M
.
Am
in
i,
“
P
h
o
sp
h
o
ru
s
Re
m
o
v
a
l
f
ro
m
Da
ir
y
W
a
ste
wa
ter
in
Ba
tch
S
y
ste
m
s
u
n
d
e
r
S
im
u
lt
a
n
e
o
u
s
A
e
ro
b
ic
/
A
n
a
e
ro
b
ic
Co
n
d
it
i
o
n
s :
A
p
p
li
c
a
ti
o
n
o
f
Re
sp
o
n
se
S
u
rf
a
c
e
M
e
th
o
d
o
lo
g
y
,
”
In
t.
J
.
En
g
.
,
v
o
l
/i
ss
u
e
:
28
(
6
)
,
p
p
.
8
5
5
–
8
6
3
,
2
0
1
5
.
[3
]
A
.
A
sa
d
i
a
n
d
A
.
A
.
L
.
Zi
a
n
ti
z
a
d
e
h
,
“
S
tatisti
c
a
l
A
n
a
ly
sis a
n
d
Op
ti
m
iza
ti
o
n
o
f
a
n
A
e
ro
b
ic S
BR T
re
a
ti
n
g
a
n
In
d
u
strial
Estate
Was
te
w
a
t
e
r
Us
in
g
R
e
sp
o
n
se
S
u
rf
a
c
e
M
e
th
o
d
o
l
o
g
y
(
RS
M
)
,
”
Ira
n
.
J
.
En
e
rg
y
En
v
iro
n
tme
n
t
,
v
o
l
/i
ss
u
e
:
2
(
4
)
,
p
p
.
3
5
6
–
3
6
5
,
2
0
1
1
.
[4
]
S
.
S
a
th
ian
,
e
t
a
l.
,
“
P
e
rf
o
rm
a
n
c
e
o
f
S
BR
f
o
r
th
e
trea
tm
e
n
t
o
f
tex
ti
le
d
y
e
wa
ste
w
a
ter :
Op
ti
m
iza
t
io
n
a
n
d
k
i
n
e
ti
c
stu
d
ies
,
”
Al
e
x
a
n
d
ri
a
E
n
g
.
J
.
,
v
o
l.
5
3
,
p
p
.
4
1
7
–
4
2
6
,
2
0
1
4
.
[5
]
A
.
H.
M
.
Ha
k
i
m
,
e
t
a
l.
,
“
A
e
ro
b
ic
slu
d
g
e
g
ra
n
u
latio
n
a
t
h
ig
h
te
m
p
e
ra
tu
re
s
f
o
r
d
o
m
e
stic
w
a
ste
wa
ter
trea
t
m
e
n
t,
”
Bi
o
re
so
u
r.
T
e
c
h
n
o
l.
,
v
o
l.
1
8
5
,
p
p
.
4
4
5
–
4
4
9
,
2
0
1
5
.
[6
]
K.
G
o
b
i,
e
t
a
l
.
,
“
De
v
e
lo
p
m
e
n
t
a
n
d
u
ti
li
z
a
ti
o
n
o
f
a
e
ro
b
ic
g
ra
n
u
les
f
o
r
th
e
p
a
lm
o
il
m
il
l
(
P
O
M
)
w
a
ste
wa
ter
trea
t
m
e
n
t,
”
Ch
e
m.
En
g
.
J
.
,
v
o
l.
1
7
4
,
p
p
.
2
1
3
–
2
2
0
,
2
0
1
1
.
[7
]
P
.
G
.
P
a
ti
l,
e
t
a
l.
,
“
Ae
ro
b
ic
S
e
q
u
e
n
c
in
g
Ba
tch
Re
a
c
to
r
f
o
r
w
a
ste
wa
ter
tre
a
t
m
e
n
t:
A
r
e
v
ie
w
,
”
In
t.
J
.
En
g
.
Res
.
T
e
c
h
n
o
l
.
,
v
o
l
/
issu
e
:
2
(
10
)
,
p
p
.
5
3
4
–
5
5
0
,
2
0
1
3
.
[8
]
W
.
L
.
L
i
e
w
,
e
t
a
l.
,
“
Co
n
v
e
n
ti
o
n
a
l
m
e
th
o
d
s
a
n
d
e
m
e
r
g
in
g
w
a
st
e
wa
ter
p
o
li
sh
i
n
g
tec
h
n
o
lo
g
ies
f
o
r
p
a
lm
o
il
m
il
l
e
ff
lu
e
n
t
trea
t
m
e
n
t :
A
re
v
i
e
w
,
”
J
.
En
v
iro
n
.
M
a
n
a
g
e
.
,
v
o
l.
1
4
9
,
p
p
.
2
2
2
–
2
3
5
,
2
0
1
5
.
[9
]
Y.
L
i,
e
t
a
l.
,
“
A
e
ro
b
ic
g
ra
n
u
lar
slu
d
g
e
f
o
r
si
m
u
lt
a
n
e
o
u
s
a
c
c
u
m
u
latio
n
o
f
m
in
e
ra
l
p
h
o
sp
h
o
r
u
s
a
n
d
re
m
o
v
a
l
o
f
n
it
ro
g
e
n
v
ia n
it
rit
e
i
n
w
a
ste
w
a
ter,
”
Bi
o
re
so
u
r.
T
e
c
h
n
o
l.
,
v
o
l.
1
5
4
,
p
p
.
1
7
8
–
1
8
4
,
2
0
1
4
.
[1
0
]
S
.
L
o
c
h
m
a
tt
e
r
a
n
d
C.
Ho
ll
ig
e
r,
“
Op
ti
m
iza
ti
o
n
o
f
o
p
e
ra
ti
o
n
c
o
n
d
it
io
n
s
f
o
r
t
h
e
sta
rtu
p
o
f
a
e
ro
b
ic
g
ra
n
u
lar
slu
d
g
e
re
a
c
to
rs b
io
l
o
g
ica
ll
y
re
m
o
v
in
g
c
a
rb
o
n
,
n
it
r
o
g
e
n
,
a
n
d
p
h
o
s
p
h
o
ro
u
s
,
”
W
a
ter
Res
.
,
v
o
l.
5
9
,
p
p
.
5
8
–
7
0
,
2
0
1
4
.
[1
1
]
L
.
Ch
e
n
a
n
d
C.
Ca
o
,
“
Ch
a
ra
c
te
risti
c
s
a
n
d
sim
u
latio
n
o
f
so
lu
b
l
e
m
icro
b
ial
p
ro
d
u
c
ts
in
m
e
m
b
r
a
n
e
b
io
re
a
c
to
rs
c
o
u
p
le
d
w
it
h
m
o
v
in
g
c
a
rriers
(
M
BR
-
M
C
),
”
De
sa
li
n
.
W
a
ter
T
re
a
t
.
,
v
o
l.
4
0
,
p
p
.
4
5
–
5
5
,
2
0
1
2
.
[1
2
]
M
.
A
.
M
o
h
a
m
m
a
d
,
e
t
a
l.
,
“
S
tat
isti
c
a
l
a
n
a
l
y
sis
a
n
d
o
p
ti
m
iza
ti
o
n
o
f
si
m
u
lt
a
n
e
o
u
s
b
io
l
o
g
ica
l
n
u
t
rien
ts
re
m
o
v
a
l
p
ro
c
e
ss
in
a
n
in
term
it
ten
tl
y
a
e
ra
te
d
S
BR,”
K
o
re
a
n
J
.
Ch
e
m.
E
n
g
.
,
v
o
l
/i
ss
u
e
:
31
(
1
)
,
p
p
.
8
8
–
9
7
,
2
0
1
4
.
[1
3
]
Y.
Zak
a
riah
,
e
t
a
l.
,
“
M
o
d
e
li
n
g
o
f
S
BR A
e
ro
b
ic
G
ra
n
u
lar S
lu
d
g
e
Us
in
g
Ne
u
ra
l
Ne
t
w
o
rk
w
it
h
G
S
A
a
n
d
IW
-
P
S
O,” i
n
Co
n
tro
l
Co
n
fer
e
n
c
e
(
AS
CC),
2
0
1
5
1
0
th
Asia
n
,
p
p
.
1
–
6
,
2
0
1
5
.
[1
4
]
R.
Da
v
a
rn
e
jad
a
n
d
A
.
S
a
h
ra
e
i,
“
In
d
u
str
ial
w
a
ste
wa
ter
trea
t
m
e
n
t
u
sin
g
a
n
e
lec
tro
c
h
e
m
ica
l
tec
h
n
iq
u
e
:
a
n
o
p
ti
m
ize
d
p
ro
c
e
ss
,
”
De
sa
li
n
.
W
a
ter
T
re
a
t.
,
p
p
.
1
–
1
3
,
2
0
1
5
.
[1
5
]
S
.
N.
S
h
im
,
e
t
a
l.
,
“
Ev
a
lu
a
ti
o
n
o
f
m
e
c
h
a
n
ica
l
m
e
m
b
ra
n
e
c
le
a
n
in
g
with
m
o
v
in
g
b
e
a
d
s in
M
BR
u
sin
g
Bo
x
–
Be
h
n
k
e
n
re
sp
o
n
se
su
rf
a
c
e
m
e
th
o
d
o
l
o
g
y
,
”
Da
sa
li
n
a
ti
o
n
W
a
ter
T
re
a
t.
,
v
o
l.
5
6
,
p
p
.
2
7
9
7
–
2
8
0
6
,
2
0
1
5
.
[1
6
]
Z.
S
ir
u
s,
e
t
a
l.
,
“
Hig
h
f
re
q
u
e
n
c
y
u
lt
ra
so
u
n
d
-
in
d
u
c
e
d
se
q
u
e
n
c
e
b
a
tch
re
a
c
to
r
a
s
a
p
ra
c
ti
c
a
l
so
lu
ti
o
n
f
o
r
h
ig
h
ra
t
e
w
a
ste
w
a
ter t
re
a
t
m
e
n
t,
”
J
.
En
v
iro
n
.
Ch
e
m.
E
n
g
.
,
v
o
l.
3
,
p
p
.
2
1
7
–
2
2
6
,
2
0
1
5
.
[1
7
]
C.
M
.
N
.
Ha
z
iq
a
h
,
e
t
a
l.
,
“
M
o
d
e
ll
in
g
a
n
d
o
p
ti
m
iza
ti
o
n
o
f
Ca
n
d
id
a
ru
g
o
sa
n
a
n
o
b
io
c
o
n
j
u
g
a
tes
c
a
t
a
l
y
se
d
s
y
sth
e
sis
o
f
m
e
th
y
l
o
lea
te b
y
re
sp
o
n
se
su
rf
a
c
e
m
e
th
o
d
o
lo
g
y
,
”
Bi
o
tec
h
n
o
l.
B
io
te
c
h
n
o
l
.
E
q
u
i
p
.
,
p
p
.
1
–
1
3
,
2
0
1
5
.
[1
8
]
M
.
N
.
R
o
y
h
a
il
a
,
e
t
a
l.
,
“
Re
sp
o
n
s
e
su
rf
a
c
e
m
e
th
o
d
o
l
o
g
y
a
p
p
ro
a
c
h
f
o
r
o
p
ti
m
izin
g
p
ro
d
u
c
ti
o
n
o
f
g
e
r
a
n
y
l
p
ro
p
io
n
a
te
c
a
tal
y
s
e
d
b
y
c
a
rb
o
n
n
a
n
o
tu
b
e
s n
a
n
o
b
io
c
o
n
j
u
g
a
tes
,
”
Bi
o
tec
h
n
o
l.
B
i
o
tec
h
n
o
l.
E
q
u
ip
.
,
p
p
.
1
–
7
,
2
0
1
5
.
[1
9
]
S
.
Z
h
iw
e
i,
e
t
a
l.
,
“
In
f
lu
e
n
c
e
o
f
tem
p
e
ra
tu
re
o
n
t
h
e
c
h
a
ra
c
teristics
o
f
a
e
ro
b
ic
g
ra
n
u
lati
o
n
in
se
q
u
e
n
c
in
g
b
a
tch
a
irl
if
t
re
a
c
to
rs,”
J
.
En
v
iro
n
.
S
c
i
.
,
v
o
l.
2
1
,
p
p
.
2
7
3
–
2
7
8
,
2
0
0
9
.
Evaluation Warning : The document was created with Spire.PDF for Python.
I
J
E
C
E
I
SS
N:
2
0
8
8
-
8708
P
a
r
a
mete
r
Op
timis
a
tio
n
o
f A
ero
b
ic
Gra
n
u
la
r
S
lu
d
g
e
a
t H
ig
h
Temp
era
tu
r
e
Usi
n
g
.
.
.
.
(
S
ya
h
ir
a
I
b
r
a
h
im)
1529
B
I
O
G
RAP
H
I
E
S
O
F
AUTH
O
RS
S
y
a
h
ira
Ib
ra
h
im
re
c
e
i
v
e
d
h
e
r
B.
En
g
.
Ho
n
s
(El
e
c
tri
c
-
Co
n
tr
o
l
a
n
d
In
stru
m
e
n
tatio
n
)
a
n
d
M
S
c
.
o
f
En
g
in
e
e
rin
g
(El
e
c
tri
c
a
l)
f
ro
m
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
lay
sia
in
2
0
1
1
a
n
d
2
0
1
5
re
sp
e
c
ti
v
e
ly
.
S
h
e
is
c
u
rre
n
tl
y
w
o
rk
in
g
to
wa
rd
h
e
r
P
h
D
i
n
p
ro
c
e
ss
c
o
n
tr
o
l
a
t
F
a
c
u
lt
y
o
f
El
e
c
rica
l
En
g
in
e
e
rin
g
,
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
la
y
sia
.
He
r
c
u
rre
n
t
re
se
a
rc
h
in
tere
sts
in
c
lu
d
e
m
o
d
e
li
n
g
o
f
n
e
a
r
-
in
f
ra
re
d
sp
e
c
tro
sc
o
p
y
a
n
d
m
e
m
b
ra
n
e
f
il
tratio
n
u
si
n
g
a
rti
f
icia
l
in
telli
g
e
n
t
sy
s
tem
.
Ir.
Dr.
No
rh
a
li
z
a
A
b
d
u
l
W
a
h
a
b
is
c
u
rre
n
tl
y
a
n
A
ss
o
c
iate
P
ro
f
e
ss
o
r
a
t
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
la
y
sia
(UT
M
).
S
h
e
is
c
u
rr
e
n
tl
y
th
e
He
a
d
De
p
a
rtm
e
n
t
o
f
Co
n
tro
l
a
n
d
M
e
c
h
a
tro
n
ic
s
En
g
in
e
e
rin
g
De
p
a
rtme
n
t
a
t
th
e
F
a
c
u
lt
y
o
f
El
e
c
tri
c
a
l
En
g
in
e
e
rin
g
,
UT
M
.
S
h
e
c
o
m
p
lete
d
h
e
r
P
h
D
in
El
e
c
tri
c
a
l
E
n
g
in
e
e
rin
g
m
a
jo
rin
g
in
Co
n
tro
l
i
n
Ju
ly
2
0
0
9
.
S
h
e
is
a
c
ti
v
e
ly
in
v
o
lv
e
d
i
n
re
se
a
rc
h
in
g
a
n
d
tea
c
h
in
g
in
th
e
f
ield
o
f
in
d
u
str
ial
p
r
o
c
e
ss
c
o
n
tr
o
l.
He
r
e
x
p
e
rti
se
is
in
m
o
d
e
ll
i
n
g
a
n
d
c
o
n
tr
o
l
o
f
in
d
u
strial
p
r
o
c
e
ss
p
lan
t.
Re
c
e
n
tl
y
sh
e
h
a
s
w
o
rk
e
d
p
rim
a
ril
y
o
n
d
if
f
e
r
e
n
t
t
y
p
e
s
o
f
d
o
m
e
stic
a
n
d
in
d
u
strial
w
a
ste
wa
ter
trea
t
m
e
n
t
tec
h
n
o
lo
g
y
to
wa
rd
s
o
p
ti
m
iza
ti
o
n
a
n
d
e
n
e
rg
y
sa
v
in
g
s
y
st
e
m
.
Dr.
A
z
n
a
h
No
r
A
n
u
a
r
is
S
e
n
i
o
r
L
e
c
tu
re
r
a
t
De
p
a
rtm
e
n
t
o
f
En
v
iro
n
m
e
n
tal
En
g
in
e
e
rin
g
,
Un
iv
e
rsiti
T
e
k
n
o
lo
g
i
M
a
la
y
sia
(
UT
M
).
S
h
e
a
lso
se
rv
e
s
a
s
Re
s
e
a
rc
h
F
e
ll
o
w
f
o
r
In
stit
u
te
o
f
En
v
iro
n
m
e
n
tal
a
n
d
W
a
ter
Re
so
u
rc
e
M
a
n
a
g
e
m
e
n
t
(IP
A
S
A
),
U
T
M
.
S
h
e
o
b
tai
n
e
d
h
e
r
Ba
c
h
e
lo
r
o
f
En
g
in
e
e
rin
g
,
M
a
s
ter
a
n
d
P
h
D
d
e
g
re
e
s
in
En
v
iro
n
m
e
n
tal
En
g
in
e
e
rin
g
f
ro
m
U
T
M
.
He
r
re
se
a
r
c
h
in
tere
sts
in
c
lu
d
e
w
a
ter
a
n
d
w
a
s
tew
a
ter
trea
t
m
e
n
t
p
a
rti
c
u
larly
in
G
ra
n
u
lar
A
c
ti
v
a
ted
S
lu
d
g
e
T
e
c
h
n
o
lo
g
y
a
n
d
su
sta
in
a
b
le en
v
ir
o
n
m
e
n
tal
m
a
n
a
g
e
m
e
n
t.
M
u
sta
f
a
Bo
b
h
a
s
c
o
m
p
lete
d
h
is
P
h
D
a
t
t
h
e
a
g
e
o
f
3
5
y
e
a
rs
f
ro
m
T
h
e
Oh
io
S
tate
U
n
iv
e
rsity
a
n
d
h
e
wa
s
a
p
o
std
o
c
t
o
ra
l
f
e
ll
o
w
w
it
h
th
e
Na
ti
o
n
a
l
Re
se
a
rc
h
Co
u
n
c
il
(NRC)
in
USA
.
H
e
is
c
u
rre
n
tl
y
a
n
a
ss
istan
t
p
ro
f
e
ss
o
r
a
t
T
a
ib
a
h
Un
iv
e
rsit
y
in
S
a
u
d
i
A
ra
b
ia.
He
h
a
s
p
u
b
l
ish
e
d
m
o
re
th
a
n
1
0
p
a
p
e
rs i
n
re
p
u
t
e
d
jo
u
rn
a
ls
a
n
d
h
a
s b
e
e
n
se
rv
in
g
a
s a re
v
ie
w
e
r
f
o
r
a
n
u
m
b
e
r
o
f
re
p
u
ted
j
o
u
r
n
a
ls.
Evaluation Warning : The document was created with Spire.PDF for Python.