CC BY
5
The influence of the biocolmation of the geobarrier for organic waste storage on the values of head drops has been investigated. A mathematical model of filtration of organic substances taking into consideration the biocolmation effect was formed. The mathematical model contains the equation of filtration under conditions of variable porosity. In addition, the mathematical model includes the equation of transfer of organic chemical substances in pore fluid and the equation of dynamics in bacteria biomass in a porous medium based on the Monod equation. The problem in the region with a thin inclusion was solved by the method of finite elements. The schematic algorithm of finding an approximate solution of the boundary problem, including the scheme of discretization over time, is presented. Numerical experiments were conducted with an analysis of their results. In particular, the tables have been given of the values of heads and their drops on inclusions when biocolmation is neglected and the values of heads and their drops on inclusions while taking into account biocolmation at specific moments of time. The numerical experiments showed that the existence of microorganisms in soil pores significantly influences the values of heads at the top and at the bottom of a geobarrier. In particular, relative changes in head drops, in comparison with the case of disregarding the influence of microorganisms, can reach 54.8 % towards an increase. Such differences, in turn, lead to a change in the estimation calculations of the propagation of waste storage contamination into groundwater. They can also cause negative changes in the stressed-strained state of a soil array in the vicinity of a geobarrier as a type of a thin inclusion and lead to the intensification of shear processes. At the same time, due to the nonlinear-ity of influences and complex interdependence among processes, it is not possible to predict such values and their differences without computer simulation and mathematical modeling
Keywords: biocolmation, organic waste, geobarrier, a finite element method, model of bacteria development
UDC 519.61/.64:627.05
pOI: 10.15587/1729-4061.2020.2100441
BIOCOLMATION AND THE FINITE ELEMENT MODELING OF ITS INFLUENCE ON CHANGES IN THE HEAD DROP IN A GEOBARRIER
O. U l y a n c h u k- M a r t y n i u k
Postgraduate Student* E-mail: o.v.ulianchukmartyniuk@nuwm.edu.ua
O. M i c h u t a PhD, Associate Professor* E-mail: o.r.michuta@nuwm.edu.ua N. Ivanchuk PhD*
E-mail: n.v.medvid@nuwm.edu.ua *Department of Computer Science and Applied Mathematics National University of Water and Environmental Engineering Soborna str., 11, Rivne, Ukraine, 33028
Received date 28.07.2020 Accepted date 13.08.2020 Published date 25.08.2020
Copyright © 2020, O. Ulyanchuk-Martyniuk, O. Michuta, N. Ivanchuk This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0)
1. Introduction
The development of industrial production and major cities actualizes the problem of today - waste storage and recycling. The design and construction of waste storage should perform the main function - to minimize the spreading of harmful substances from storage to the environment. One of the design solutions to this problem is the use of so-called geobarriers in soil storages [1, 2]. A geobarrier is a thin layer of natural or (and) synthetic material, which (layer) has pronounced non-penetration properties for harmful substances. A change in the physical-chemical protective properties of the geobarrier material due to the influence of external factors remains a relevant direction of scientific research, as it is related to the safety of life and health of people.
In the presence of organic wastes in the storage, one of the factors of influence on the geobarrier parameters is the development of microorganisms, and, as a consequence, biocolmation. Biocolmation may be the cause of changes in hy-drological parameters of the porous material of a geobarrier. This, in turn, leads to changes in estimation calculations of contamination propagation from waste storage to the
groundwater. At the same time, in virtue of non-linearity of the influences and complex interdependence of processes, it is impossible to predict the values of the amount of harmful substances that got to the environment through a geobarrier without computer simulation and mathematical modeling.
2. Literature review and problem statement
In paper [4], it was proposed to take into account a change in the hydrological parameters of geobarriers in the bases of waste storage in case of the influence of inorganic chemical substances and temperature. It was mathematically expressed in the modification of the conjugation conditions for humidity, contamination concentration and temperature in comparison with the classical case. In addition, in paper [5], the effect of chemical suffusion of the geobarrier material and its (material) possible structural changes under the influence of aggressive chemical solutions of inorganic origin were theoretically taken into consideration. In studies [3, 6], the method of finite elements was applied to the numerical solution of the corresponding boundary problems
with modified conjugation conditions. Using the example of numerical experiments, it was shown that in some cases it is not permissible to neglect nonlinear dependences of physical-chemical parameters of the material of thin inclusions, and consequently - modification of conjugation conditions.
However, these works do not explore the peculiarities of spreading organic substances from storages of domestic wastes, as a rule. One of these features is the development of microorganisms (which is fundamentally impossible in the presence of inorganic chemical waste) and, consequently, biocolmation of the porous material of a geobarrier. We aim to consider this factor of influence by constructing an appropriate mathematical model and its numerical research. While the main components of the mathematical model concerning the motion of the pore fluid, the contamination propagation, modification of conjugation conditions at a change in the properties of the material of a geobarrier are explored in research, for example, in papers [6] or [4], the development of microorganisms in the porous medium requires special attention.
The first step in this direction is an overview of available scientific research, which can be conditionally divided into two subgroups. The first one relates to the effect of microorganisms on the characteristics of porous media. The aim of such an overview is to show the existing scientific contribution at the level of field experiments, which once again emphasizes the relevance of the research and provides the grounds for mathematical and computer modeling. The second group refers to proper mathematical models of development and propagation of microorganisms in media, including porous ones.
Paper [7] provides a thorough overview of studies into the mechanisms of the biocolmation of porous media. In relation to mathematical modeling of biocolmation, its consequences, such as a decrease in filtration factor of a porous medium, were studied. The authors analyzed the analytical formulas for filtration factor, well-known in the scientific literature, however, mathematical models of the dynamic change of microorganisms in the porous medium were not considered and explored.
Paper [8] describes the field experiments on studying the influence of biocolmation on hydroconductivity of sandy soil. The author presented a thorough analysis of the known dependences of filtration factor on a change in porosity as a result of germs activity. Based on a comparison of his own experiments with the known ones, the author concluded that a decrease in hydroconductivity, in this case, is described by the models of colonies development, rather than the models based on biofilms. However, the author did not go further in his research - application of the obtained dependences in mathematical models.
The important problem of denitrification in agricultural production was studied in paper [9]. During the experiments, as a part of general dynamics, the development of microbes as one of the factors of denitrification (neutralization and decomposition of nitrates) was studied. This research shows that bio-processes in porous media, including the development of microorganisms, are important in different problems and their exploration is a relevant task of today. The results of research [10] also prove it. These studies can be applied to research into the problems of waste neutralization at their propagation from storage. However, it was not done in the above papers.
The effect of microbial development on the filtration factor of a mixture of sludge and bentonite (three samples
with 0 %, 5 %, 10 % of bentonite content) was studied in article [11]. When studying the filtration of nutrients solution on day 540 from the beginning of the research, the filtration factor of all soil samples decreased almost by an order of magnitude. If researchers subsequently replaced nutrients with antibiotics and antifungal drugs, on day 720 (180 days after the replacement), the filtration factor acquired the values characteristic to the beginning of the experiment. In the article, the results of experiments were compared with the known theoretical dependences of other researchers concerning the influence of biocolmation on filtration factor and good match both in qualitative and quantitative terms was obtained. The researchers see the benefits of their research in the fact that it shows the possibility of influencing and controlling in terms of hydro-conductivity and improving the efficiency of geobarriers of waste storage. However, the researchers did not apply the obtained results to the improvement of mathematical and computer models of the relevant processes.
Article [12] can be divided into two parts. In the first part, the authors perform a theoretical assessment of the impact of biomass existence in the porous medium on a change in permeability (in the sense of fluid filtration) of the mentioned medium. In particular, the authors used the approaches based on functional dependences by van Genuchten, on functional dependences by Brooks-Corey. And both cases gave the same result.
kb f
ka
\
1--
n0 7
19/6
where k0, kb are the factors of permeability of «clean» porous medium and the medium under the influence of germs activity; nf is the volume fraction of the biomass; n0 is the soil porosity. The authors termed the proposed approach a «microscopic approach». The other part of the article compares the obtained dependences with the models based on biofilm development. The studies showed the coincidence of results with relatively high accuracy. The mathematical model of the dynamics of microorganisms in a porous medium, proposed by the authors, includes the following three main equations. The first is the differential equation of nutrients propagation. The second is the differential equation for microbe propagation in the pore fluid. The third is a differential equation that describes the dynamics of a change in the number of microorganisms, absorbed by the skeleton of the porous medium. However, the porous medium considered by the authors is homogeneous and does not contain thin inclusions.
In paper [13], the field experiments revealed that the magnitude of the filtration factor of soil samples decreased by order of magnitude 12 days after the beginning of tests on the filtration of organic substances. Using the obtained values, the authors performed finite-element estimation calculations using a one-dimensional mathematical model of contamination propagation. The time of concentration breakthrough beyond the protective geobarrier of waste storage was studied. At the factor of clean water filtration, such breakthrough occurs after 25 years, and at the values of filtration factor that are by one order of magnitude lower - after 90 years. However, the considered model is linear, it was explored at various constant values of filtration factor, that is, its dynamic change in the process of numerical computer experiments is impossible. In addition, the medium was considered without the existence of fine inhomogeneous inclusions in it.
The authors of article [14] notice that the Monod model of microbe growth historically became the first model of dynamics of microorganisms:
d£ R D S
— = u,B = ux B--,
di P Pmat ks + S
where, according to the authors' designations and definitions, B is the biomass, |imax is the specific growth rate, S is the substrate concentration, and ks is the Monod constant of semi-saturation, that is, the substrate concentration when |i = 0.5^max. The specified point model has many classifications that include the existence of several types of microorganisms, their interaction, the existence of biochemical reactions, etc.
However, the model based on the classical Monod equation has a drawback - the biomass of microorganisms in the presence of a favorable medium can increase up to infinity. In paper [15], it was emphasized that one of the ways to eliminate this drawback is to use a logistic equation of development of microorganisms:
The specific volume of biophase is nbio. It is assumed that biomass X changes according to the Monod equation:
dX =m X H-^io--\i. X.
However, the proposed general models are not applicable to studying the processes in non-homogeneous media when there are geobarriers.
The author of article [17] notes that the biomass influences filtration processes in the soil due to a change in porosity, and hence it influences filtration factor. In the article, the author made a review of such dependences of filtration factor on porosity, known from scientific sources. However, the author did not apply the above-mentioned dependences in mathematical models for moisture transfer and filtration and did not evaluate the impact of such dependences on changes in heads and humidity in soil.
The model of biocolmation of a porous medium, developed and studied by the authors of article [18], includes: 1) the equation of mass transfer in a porous medium:
dc.
dt
' = №bacCh
1 _ bac
сГ"
bac
dCmob =_d_
dt dy
D uu mob
. L дУ
9c„
dy
where Cbac is the biomass of microorganisms, |imax is the growth factor, CaC is the maximum permissible biomass of microorganisms in a unit of volume of the medium (the author's designations).
The authors of article [16] proposed the so-called TBCmodel of propagation of organic substances in the porous medium. The essence of the model is that the elementary volume of a porous medium conditionally consists of three phases: liquid phase; solid phase; and biophase. In particular, a known equation of elastic filtration was used as the filtration equation:
where cmob is the concentration of a dissolved substance in the moving phase; DL is the coefficient of longitudinal dispersion; v is the average porous velocity; t is the time; y is the distance; s represents changes in the concentration of the moving phase of the types due to the chemical or biochemical reaction in a porous media;
2) the model of bacteria growth:
dX _ c1 c2 dt matksl + ci ks2 + c
X,
dt
w,
where s0 is the storage factor; w is the source; kf is the filtration factor; h is the piezometric drop. To change the concentration of the solution in pore water, the equation of the parabolic type was used:
Эс
dt
~=_V{uCmob ) + V{DVCmob ) + Г,
where cmob is the concentration of the dissolved substance in aqueous phase; D is the dispersed tensor, including molecular diffusion; u = v/n is the velocity of fluid motion; v is the filtration rate; n is the soil porosity; r is the function of sourced and drains of chemical substances. The exchange of dissolved substances between the water phase and the biophase is modeled with the help of [16]:
dc,
bio . dt
X
and
dCmob X
dt
=_( с _с
V bio m
n
where Cbio is the concentration of the dissolved substance in the biophase; X is the parameter of the changed course.
where vmax is the maximum growth rate; c1 is the concentration of basic substrate in biophases; c2 is the concentration of secondary substrate in biophases; ks1 is the constant of semi-saturation of the basic substrate; ks2 is the constant of semi-saturation of the secondary substrate; X is the general concentration of bacteria;
3) proper model of biocolmation which implies the dependence of a change in porosity on the bacteria concentration. That is,
(1 _£0 )PmQs
Ps (1 _ f) ,
£ = £0 _ V,
where v is the total volume of bacteria growth; £0 is the initial porosity before filtration; £ is the porosity after filtration; Q is the volume fraction of clogged or precipitated material (Qa=AX), A is the of sedimentation rate;pm is the density of soil-sawdust material; ps is the density of clogged or precipitated material; fa is the secondary porosity of precipitated material. The Kozani formula was used for the dependence of filtration factor on porosity.
Consequently, biocolmation as one of the consequences of microbial development leads to a change in hydrological parameters of the porous material of a geobarrier. However, in the known mathematical models of contamination propagation from waste storage facilities in the presence of a geobarrier, such changes were not taken into consideration. This is due to the difficulties of further exploration of such
a mathematical model. Firstly, it requires modification of conjugation conditions for unknown functions. Secondly, this will require the use of effective numerical algorithms for solving a non-linear mathematical model. However, without the development of mathematical models in this direction, it is impossible to take into account a large number of field experiments and, consequently, to specify forecasts.
stituents should also be the equations of transfer of organic chemical substances in the pore fluid of the porous medium and the equation of the dynamics of bacteria biomass in the porous medium based on the Monod equation. As a result, we obtain the following boundary problem:
„dh 9 ,dh^ „ „ . ...
P— = —I k(n)— , xeQ1 uQ,, t>0, (1)
dt dx ^ dx)
3. The aim and objectives of the study
The purpose of this study is to quantify the effect of the existence of microorganisms in soil pores on the value of head in a porous medium in general and, specifically, on the value of head drops during the transition through a thin geobarrier. This will make it possible to draw conclusions about the significant or insignificant changes in the predictive estimation of contamination propagation from waste storage facilities into groundwater.
To accomplish the aim, the following tasks have been set:
- to construct a mathematical model of the filtration of organic substances taking into consideration a biocolmation effect. Since the biocolmation effect causes a change in the hydrological parameters, in particular, a filtration factor, this phenomenon should be taken into consideration under conjugation conditions for heads;
- to find the numerical finite-element solutions to the corresponding non-linear boundary problem under the modified conjugation conditions, which describes the constructed mathematical model;
- to conduct a series of numerical experiments and perform their analysis.
4. Construction of a mathematical model of the filtration of organic substances taking into consideration the biocolmation effect
The article explored the process of the propagation of organic chemical substances in a soil layer of total thickness l with a thin inclusion of thickness d, which is placed at depth x=% (Fig. 1).
Ü,
Q,
x = 4
x = 4
h(xt)L = ho (t) t -
u (x, *L=(-k(n)t )
= 0, t - 0,
h(x,0) = h°(x), xeQi UÜ2,
dc д ( dc ) dc цта1
n— = —I nD— - u--LJmaL B
dt dx I Эх ) dx Y
x eQ1 uQ2, t > 0,
c (x, t)L = c0 (t) t - a
1 -
B
B
k, + c
qc(x,t L=- nD dx
= 0, t - 0,
c (x ,0) = c0 (x ), x eQi UÜ2, B
dB R
—=дт B dt
1 --
B
k + c
x eQ1 uQ2, t > 0, B(x,0) = B0 (x ), x eQ1 uQ2, ..±1 _ [h]
\x=%
\x=%
J
dx
k Y(
• dx
,nY DY
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9) (10) (11)
(12)
Fig. 1. Soil layer of thickness / with a thin inclusion of thickness d (d<< 1)
To formulate the mathematical model, we will partially use the findings by scientists from the reviewing chapter of the article. The mathematical model will include the filtration equation under conditions of variable porosity. Its con-
Here:
- Ц = (0;%), Q2 = l), 0 < l, ho (, ho (x), c0 (x), c0 (), B0 (x) are the known functions;
- P is the factor of the elastic capacity of soil;
- h is the head;
- k, kY are the factors of filtration of basic soil and inclusion soil, respectively;
- u is the filtration rate;
- n, n< are the porosity of soil and geobarrier material, respectively;
- qc is the flow of chemical substances;
- is the biomass decomposition factor;
- u1, q± are the values of the rates of filtration and of the flows of chemical substances at x= £,-0 and x= £,+0, respectively;
- [h] = h+- hT, [c] = c+ — С are the head drops and of concentration of organic chemical substances on thin inclusion;
- c is the concentration of the pore solution of the organic chemical substance;
4c,
0
x
- D is the factor of diffusion of an organic substance;
- ^max is the maximum rate of biomass growth;
- B is the concentration of bacteria biomass in the porous medium;
- Bmax is the maximum possible concentration of bacteria biomass in the porous medium;
- kc is the parameter in the equation of biomass dynamics, at which the rate of biomass growth is equal to half the maximum;
- Y is the bacterial yield.
Conjugation conditions (11), (12) differ from the classical ones and take into consideration the dependences of soil geobarrier parameters on porosity which, in turn, depend on the biomass concentration. Conditions (11), (12) are derived similarly to those from research [3].
5. Solving the problem on the filtration of organic substances taking into consideration the biocolmation effect in the area with inclusions by using a finite element method
Similarly to [19], let us assume that H0 is the space of vector-functions { (x); s2 (x); s3 (x)}, each component of which on each of the intervals (0;^), (£,;/) belong to the Sobolev space W2,1 (Q), in this case, the first two acquire zero values at the end of the segment [0,/], where the boundary conditions of the first kind are assigned for functions h(x, t) and c(x, t), respectively. Let us assume that H is the space of functions {v1 (x,t); v2 (x,t); v3 (x,t)), each component of which is integrated with the square along with its first variables dvjdt, dvjdx, i = 1,3 on each of intervals (0;^), (£,;/), Vt e(0;T], in this case, the first two satisfy the boundary conditions of the first kind like functions h (x, t), c (x, t), respectively. Here T>0.
Take {(x); s2 (x); s3 (x)}eH0. Multiply equation (1) and initial condition (4) by si(x), integrating them on the segment [0, /] and taking into account conjugation conditions (11), obtain:
JP^5! (x)dx +
+ jk („)¥ iM dx + iM- = 0, Vt e(0J ],
dx dx
d i:
dx
k Mn
(13)
J h (x,0)s1 (x )dx = J \ (x )s1 (x )dx, Vt e(0;T ]. (14)
00
Similarly, for the concentration of organic chemical substances and biomass, we obtain:
in —s, (x )dx + inD — d^2 (x) dx + J0 dt 2W 0 9x dx
+ Ju^ (x)dx + M +
dx
nD
Y
1 -
B
B
B——s2 (x)dx = 0, Vt e(0;T], (15) k + c
dt
-i
1-
B
B
k, + c
Bs3 (x)dx = 0, Vt e(0;T], (17)
J B(x,0)s3 (x )dx = J B0 (x )s3 (x )dx, Vt e(0;T ]. (18)
00
Attribute 1. Function {h(x,t);c(x,t);B(x,t)}eH, which for any {s1 (x); s2 (x); s3 (x)}eH0 satisfies integral ratios (13)-(18), is called the generalized solution to the boundary problem (1) to (12).
We seek for the approximated solution to the boundary problem (1) to (12) in the form of:
h (x, t ) = £ h (t ^ (x ), c (x, t ) = £ c (t )(|>!2 (x ),
i=1 i=1
B (x, t ) = £B (t )(i3 (x ),
(19)
where hi(t), ci(t), Bi(t), i = 1, N, are the unknown factors that depend only on time. If we consider (19) as the approximated finite-element solution, functions (x), ^i2 (x), ^i3 (x), i = 1, N, are the polynomial basis functions with a finite carrier. To simplify further calculations and in regard to numerical experiments, we will consider that the same grid of finite elements is used for approximated finding h(x, t), c (x, t) and B (x, t). In this case, (x) = ^i2 (x) = ^i3 (x) and to avoid double indexation, designate basic functions as ^ (x), i = 1, N.
Then, from a weak statement (13) to (18) of problem (1) to (12), taking into account (19), we will obtain:
dH
M1 ■ — + L1 (B)H (t ) = 0,
M 1 ■ H(0) = ]F1,
1
M2 (B) — + l2 (h, b) c (t ) = 0, M 2 ■ c(0) = f2,
dB
M3 ■ — + L3(C,B) B(t) = 0, M3 ■ B(0) = F3,
where
H = (( (t))) C = (ci (t))) B = ( (t))) H(0) =(( (0))) C(0) = (ci (0))), B(0) = ( (0))), Mk = (m(k))" Lk = (/(k))" Fk f))
(20) (21) (22)
(23)
(24)
(25)
jc(x,0)s2 (x)dx = jc0 (x)s2 (x)dx, Vt e(0;T], (16)
]k ft m=№ ^
X
c
-Jk) -
= J>s 4;dr, k -1,3,
< -J P4, 4 jdx, f^ -J h0$t dx,
' A [>]|>j]
j-J k (-)dx ^ dx +
dx
k Y(
m(2) - J «4i 4 j ^ fi2) - J C04i ^
0 q
r-Ld ^ dk dx +L ^ 4, dx+
J /-j-v rl -V J rl-r
J y
1 -
5
B„
0 dx
B ^ , [4, ][4j]
—4,4, àx + L ¿JL J, C dx
J n^D
A + C '
mf -J4i 4 j dx, f(3) -J Bo4i dx,
00
1 / / \ B
l(3) --J ij
V v
1-
B
k + c
4 j 4idx.
Step 2. To find H(p+1), perform discretization of the system of non-linear differential equations (20) in time with the use of completely implicit linearized differential scheme [6]:
Mi
H
( p+i) _ h( p)
Li (B
( p)\
H
(p+i) _
- 0, p -1,2,...
where H(p) - H (tp), B(p) - B (tp), tp - pT. After the formation of the SLAE:
G ■ H( p+1) - F, where
G -1M1 + L1(B(p)),
T
F -- M1 ■ H( p), T1
The system of equations (20) to (25) is the Cauchy problem for a system of non-linear differential first-order equations. Finding its solution also requires the use of appropriate discretization schemes.
the main boundary condition - the boundary condition of the first kind (2) - must be taken into consideration. It is considered by the implementation of the two following steps:
1) all the elements of the first SLAE equation, including a free term, are zeroed;
2) the first element of the first equation is equal to 1, and a free term is taken equal to ho (ip+1 ).
To find the integrals during the formation of G and F, we use numerical integration.
Step3. Similarly, we find C(p+1) from (22) and B(p+1) from (24).
Step 4. Increase p by 1. If (p+1) > M, complete calculation. Otherwise, proceed to step 2.
T
6. Schematic algorithm of the finite-element solution to problem (1) to (12) and time sampling
Step 0. The actions of this step are based on the ideas from [19]. Cover [0;k]u [k;l] with a finite-element grid. In this case, there are such two finite elements [xm-1; xm ]e[0; k] and
[ Xm+1; xm+2 ] e [l] that x m = xm+1 = k and xm , xm+1 £ [ l].
The basis functions 4m (x) and ^m+1 (x) are discontinuous at x=k, that is:
6 (x) = 1, A (x)| = 0,
4m+1 (x !x=k-0 = 0, Am+1 (x ^x=k+0 = 1
Step 1. Find H(0), C(0) and B(0). In practice, to perform this step of the system of linear algebraic equations (SLAE), (21), (23) are (25) not to be solved. Taking into consideration the properties of basic functions of the method of finite elements (4; (x; ) = 1, i = 1, N), we find the values of factors h(t), c(t), bi(t), i = 1, N at the initial moment of time t = 0 directly from initial conditions (4), (8), (10):
h (0)=h (x), c (0)=c0 (x), b (0)=b0 (x), i=1N.
Assign value 0 to index p. For simplification, time sampling will be performed on a uniform grid at pitch t. Assign the number of pitches M in time.
7. The results of numerical experiments on solving the problem of the filtration of organic substances taking into consideration the biocolmation effect
Parameters of soils for numerical experiments were taken from the freely available software Hydrus-1D. In particular, Silt Loam at k0 = 0.108 m/day, n0 = 0.45, where index «0» means initial values, was considered as the basic soil. The clay with the following parameters k0Y = 0.0048 m/day, n'Y = 0.36 was used as the thin inclusion soil:
The soil layer of thickness l = 10 m was considered for a model problem. The depth of the inclusion deposit was k=2 m, and its thickness d = 0.2 m. The pitch of variable x was 0.02 m. The time pitch t = 1 day.
The equation of filtration consolidation [20], where P = Yca/(1+e), e = n/(1 - n) is the soil porosity factor, a = = 5.12 ■ 10-5 m2/H is the soil compactness factor, yc = 104 H/m3 is the specific weight of pore solution, was considered as equation (1). The initial distribution of heads h0 (x ) = 10 m, which corresponds to the application to soil of the corresponding load - wastes that are the source of the organic substances. The free fluid outflow is ensured at the lower boundary, while there is no drain at the upper boundary.
Glucose solution was considered as an organic solution. The initial glucose distribution in pore water c0 (x ) = 1mM. In boundary condition (6) on the soil surface c0 (t) = 10 mM. The factor of glucose diffusion is D = 0.5184 ■ 10-4 m2/day.
Dependence [12] was used for the dependence of filtration factor on porosity:
19/6
where nB is the biomass volume in the unit of volume of the porous medium. To convert the biomass concentration to volume, use the assumption that 80 % of biovolume is water, the remaining 20 % is dry mass - 50 % carbon (bio-carbon) [21]. That is,
B
B 0.8pw + 0.2 ■ (pj 2)'
where pw is the density of water; pc is the density of bio-carbon. The complexity is that, for example, for bacteria colonies, value pc varies in the range from 50 to 650 kg/m3 [17]. In numerical experiments pc = 100 kg/m3 parameters of biomass kinetics are the following: |amax = 0.2 day-1, Bmax = = 100 kg/m3, iidec = 0.01day-1, kc = 0.06 mM, Y = 50.
Piecewise-quadratic functions were used on the basis functions of the MFE. The results of numerical experiments are given in Tables 1, 2.
Table 1
Values of heads and their drops on the inclusion when biocolmation is neglected
Moment of time h- h+ [h]
t=25 days 9.639 8.936 0.703
t=50 days 8.390 7.090 1.300
t = 100 days 5.689 4.550 1.139
t=200 days 2.431 1.925 0.506
t=300 days 1.032 0.817 0.215
Table 2
Values of heads and their drops on the inclusion when biocolmation is taken into consideration
Moment of time h- h+ [h]
t= 25 days 9.770 9.157 0.613 (-12.8 %)
t= 50 days 9.100 7.515 1.585 (+21.9 %)
t = 100 days 6.919 5.156 1.763 (+54.8 %)
t=200 days 3.214 2.445 0.769 (+52.0 %)
t=300 days 1.404 1.093 0.311 (+44.6 %)
Table 2 corresponds to the use of the modified conjugation conditions, and the results of Table 1 - to the classic ones [19]. Table 2 shows a relative increase (as a percentage with a «plus» sign) or a relative decrease (with a «minus» sign) in the values of head drops on a thin inclusion in comparison with the case when the biocolmation influence is neglected (the values in Table 1).
8. Discussion of results of solving the problem on the filtration of organic substances taking into consideration the biocolmation effect
As the comparison of data from Table 1 and Table 2 reveals, considering a dynamic change in the microorganisms in the porous medium affects the value of heads from the bottom and from the top of a protective geobarrier. And while the values of heads at the top (h-) are determined by the existence of microorganisms in the soil before a geobarrier, the
values at the bottom (h+) are determined by the development of microorganisms in the thin inclusion itself and the use of modified conjugation conditions (Table 2). These changes are clearly traced through the values of drops [h]: in the presence of microorganisms, excessive heads, as a consequence of the applied external loads from wastes, propagate more slowly in the vicinity of a geobarrier. In this case, due to the nonlineari-ty of the influences and complex interdependences of processes, it is impossible to predict such values and their differences without computer simulation and mathematical modeling. In practice, our research can be used to assess the propagation of organic contamination from waste storage facilities, because their dynamics depends on the head gradient (drop) in a geobarrier. In addition, since clays are the material for artificial geobarriers, the changes in heads in the neighborhood of a thin inclusion (especially, an increase in drops, which means an increase in the gradient of a drop) can lead to the formation of landslide hazardous zones, which is especially important for large waste storage facilities. Such clarification in the estimation calculations can be considered in the preparation of design documentation and engineering features of organic waste storage facilities. However, the hydrological characteristics of soils and the parameters of the dynamics of microorganisms' growth should be known to make such calculations possible. The problem of the relationships between biomass and biological volume of microorganisms is far from trivial. It imposes some constraints on the research results. This study will be developed in the following way:
1) the consideration of design features of waste storage facilities and the arrangement of a geobarrier, which is possible by increasing the dimensionality of the problem;
2) taking into account the models of the dynamics of microorganisms based on the differential equations in the partial second-order derivatives;
3) theoretical study into the accuracy of the derived finite element solutions for boundary problems with the modified conjugation conditions.
9. Conclusions
1. A mathematical model of the filtration of organic substances in the soil array in the presence of a thin geobarrier was constructed taking into consideration the biocolmation effect, which includes the filtration equation under conditions of variable porosity. Since the hydrological parameters of the materials of geobarriers depend greatly on the activity of microorganisms, the model uses the idea of modifying the conjugation conditions for heads, taking into consideration these effects. In addition, the components of the model are the equation of the transfer of organic chemical substances in the pore fluid of the porous medium and the equation of the dynamics of bacteria biomass in a porous medium based on the Monod equation.
2. Numerical solutions to the corresponding nonlinear boundary problem were derived by using the method of finite elements (MFE). The schematic algorithm of the finite-element solution to the obtained problem, including the scheme of time sampling, is presented. The application of the MFE to the problems of this type is substantiated by the necessity of finding generalized solutions, in which discontinuities of the first kind are possible. Based on the developed algorithms, a software suite for carrying out numerical experiments has been developed, which implies taking into account the modified conjugation conditions.
k
n
0
3. Computer simulation has revealed that the presence of microorganisms in the pores of soil significantly affects the values of heads at the top and at the bottom of a geobarrier, hence, the values of head drops. In particular, the biocolmation of a geobarrier (based on the solution to a model problem) can lead to a relative increase in a head drop by 54.8 % compared to the case of disregarding the activity of microorganisms. Such differences, in turn, lead to changes in the predictive calculations of contamination propagation from waste storage facilities into groundwater and can cause negative changes in the stressed-strained state of the soil array in the geobarrier neighborhood as a thin inclusion. At the same time, due to the nonlinearity of influences and complex interdependence of the processes, it is
impossible to predict such values and their differences without computer simulation and mathematical modeling. In addition, the hydrological characteristics of soils and the parameters of the dynamics of microorganism growth should be known.
Acknowledgment
This study is funded within the framework of the scientific-research work «Mathematical modeling and computer simulation of technogenic controlled processes in porous media with barriers under conditions of identification» (Ukraine State Registration No. SRN 0120U102055).
References
1. Kong, D.-J., Wu, H.-N., Chai, J.-C., Arulrajah, A. (2017). State-Of-The-Art Review of Geosynthetic Clay Liners. Sustainability, 9 (11), 2110. doi: https://doi.org/10.3390/su9112110
2. Scalia, J., Bareither, C. A., Shackelford, C. D. (2018). Advancing the use of geosynthetic clay liners as barriers. Geotechnical Engineering Journal of the SEAGS & AGSSEA, 49 (4), 100-114. Available at: https://www.researchgate.net/publication/ 329544978_Advancing_the_use_of_geosynthetic_clay_liners_as_barriers
3. Martyniuk, P. M., Michuta, O. R., Ulianchuk-Martyniuk, O. V., Kuzlo, M. T. (2018). Numerical investigation of pressure head jump values on a thin inclusion in one-dimensional non-linear soil mousture transport problem. International Journal of Apllied Mathematics, 31 (4), 649-660. doi: https://doi.org/10.12732/ijam.v31i4.10
4. Chui, Y., Martyniuk, P., Kuzlo, M., Ulianchuk-Martyniuk, O. (2019). The conditions of conjugation in the tasks of moisture transfer on a thin clay inclusion taking into account salt salutions and temperature. Journal of Theoretical and Applied Mechanics, 49 (1), 28-38. doi: https://doi.org/10.7546/jtam.49.19.01.03
5. Ulianchuk-Martyniuk, O., Michuta, O. (2020). Conjugation conditions in the problem of filtering chemical solutions in the case of structural changes to the material and chemical suffusion in the geobarrier. JP Journal of Heat and Mass Transfer, 19 (1), 141-154. doi: https://doi.org/10.17654/hm019010141
6. Ulianchuk-Martyniuk, O. V. (2020). Numerical simulation of the effect of semi-permeable properties of clay on the value of concentration jumps of contaminants in a thin geochemical barrier. Eurasian Journal of Mathematical and Computer Applications, 8 (1), 91-104. doi: https://doi.org/10.32523/2306-6172-2020-8-1-91-104
7. Baveye, P., Vandevivere, P., Hoyle, B. L., DeLeo, P. C., de Lozada, D. S. (1998). Environmental Impact and Mechanisms of the Biological Clogging of Saturated Soils and Aquifer Materials. Critical Reviews in Environmental Science and Technology, 28 (2), 123-191. doi: https://doi.org/10.1080/10643389891254197
8. Seki, K. (2013). Biological Clogging of Sand Columns. Open Journal of Soil Science, 03 (03), 148-152. doi: https://doi.org/ 10.4236/ojss.2013.33017
9. Knabe, D., Kludt, C., Jacques, D., Lichtner, P., Engelhardt, I. (2018). Development of a Fully Coupled Biogeochemical Reactive Transport Model to Simulate Microbial Oxidation of Organic Carbon and Pyrite Under Nitrate-Reducing Conditions. Water Resources Research, 54 (11), 9264-9286. doi: https://doi.org/10.1029/2018wr023202
10. Moshynsky, V., Riabova, O. (2013). Approaches to Aquatic Ecosystems Organic Energy Assessment and Modelling. NATO Science for Peace and Security Series C: Environmental Security, 125-135. doi: https://doi.org/10.1007/978-94-007-6152-0_12
11. Glatstein, D. A., Francisca, F. M. (2014). Hydraulic conductivity of compacted soils controlled by microbial activity. Environmental Technology, 35 (15), 1886-1892. doi: https://doi.org/10.1080/09593330.2014.885583
12. Clement, T. P., Hooker, B. S., Skeen, R. S. (1996). Macroscopic Models for Predicting Changes in Saturated Porous Media Properties Caused by Microbial Growth. Ground Water, 34 (5), 934-942. doi: https://doi.org/10.1111/j.1745-6584.1996.tb02088.x
13. Tang, Q., Gu, F., Zhang, Y., Zhang, Y., Mo, J. (2018). Impact of biological clogging on the barrier performance of landfill liners. Journal of Environmental Management, 222, 44-53. doi: https://doi.org/10.1016/j.jenvman.2018.05.039
14. Thullner, M., Regnier, P. (2019). Microbial Controls on the Biogeochemical Dynamics in the Subsurface. Reviews in Mineralogy and Geochemistry, 85 (1), 265-302. doi: https://doi.org/10.2138/rmg.2019.85.9
15. Bajracharya, B. M., Lu, C., Cirpka, O. A. (2014). Modeling substrate-bacteria-grazer interactions coupled to substrate transport in groundwater. Water Resources Research, 50 (5), 4149-4162. doi: https://doi.org/10.1002/2013wr015173
16. Thullner, M., Schroth, M. H., Zeyer, J., Kinzelbach, W. (2004). Modeling of a microbial growth experiment with bioclogging in a two-dimensional saturated porous media flow field. Journal of Contaminant Hydrology, 70 (1-2), 37-62. doi: https:// doi.org/10.1016/j.jconhyd.2003.08.008
17. Thullner, M. (2010). Comparison of bioclogging effects in saturated porous media within one- and two-dimensional flow systems. Ecological Engineering, 36 (2), 176-196. doi: https://doi.org/10.1016/j.ecoleng.2008.12.037
18. Eljamal, O., Jinno, K., Hosokawa, T. (2008). A mathematical model of biological clogging of soil-sawdust media. Journal of Environmental Hydrology, 16 (4), 1-12. Available at: https://www.researchgate.net/publication/281260984_A_mathematical_ model_of_biological_clogging_of_soil-sawdust_media
19. Sergienko, I. V., Skopetskiy, V. V., Deyneka, V. S. (1991). Matematicheskoe modelirovanie i issledovanie protsessov v neodnorodnyh sredah. Kyiv: Naukova dumka, 431.
20. Vlasyuk, A. P., Martynyuk, P. M., Fursovych, O. R. (2009). Numerical solution of a one-dimensional problem of filtration consolidation of saline soils in a nonisothermal regime. Journal of Mathematical Sciences, 160 (4), 525-535. doi: https://doi.org/ 10.1007/s10958-009-9518-8
21. Bölter, M., Bloem, J., Meiners, K., Möller, R. (2005). Enumeration and biovolume determination of microbial cells. Microbiological Methods for Assessing Soil Quality, 93-113. doi: https://doi.org/10.1079/9780851990989.0093
-□ □-
This paper reports a study on the application of aluminum sulfate solution, modified by the magnetic field and electrocoagulation, in the processes of drinking water preparation. The modification of the coagulant solution makes it possible to intensify water purification processes, to reduce reagent consumption by 25-30 %. It has been found that a dose of the modified aluminum sulfate solution of28-30 mg/dm3 improves the efficiency of removal of suspended substances and coloration by 35-40 %. The dosage of the conventional reagent solution was 40 mg/dm3 while reaching the same purification parameters.
Modifying a solution of aluminum sulfate with the magnetic field and electrocoagulation increases the hydraulic size of the coagulated suspension. A change in the hydraulic size in the suspension has been studied at different periods of the year. In winter, when treating water with the modified aluminum sulfate solution, there a decrease in the suspension content whose hydraulic size is 0.1 mm/s and less, from 89 % to 22 %. In this case, the content of suspended substances at settling decreases from 8.5-12.5 mg/dm3 to 5.6-8.3 mg/dm3. In spring, when using the modified coagulant solution, the content of suspension whose size is 0.1 mm/s and less decreased from 55 % to 15 %. In summer, there is an increase in the content of suspension whose size is 0.3-0.5 mm, from 58 % (a conventional reagent solution) to 66 % (the modified reagent solution). This indicates an intensification of the coagulation of impurities and the clarification of water.
The experimental data testify to an increase in the effectiveness of discoloration of natural low-turbid colored water to 63.3-63.9 % for the modified reagent solution at 45.5 % for a conventional reagent solution.
A change in the bacteriological parameters has been determined: the effectiveness of the decrease in a microbial number grows from 11.6-18.7 % to 18.6-25.1 %. In terms of a coli-index, the efficiency of purification grows from 16.6-23.1 % to 23.0-29.5 %
Keywords: drinking water quality, modification of reagent solution, coagulation, aluminum sulfate, hydraulic size, coloration, suspended substances, water clarification, magnetic field,
anode-dissolved iron -□ □-
UpC 628.16
pOI: 10.15587/1729-4061.2020.210096|
APPLYING A MODIFIED ALUMINUM SULFATE SOLUTION IN THE PROCESSES OF DRINKING WATER PREPARATION
S . D u s h ki n
PhD, Associate Professor Department of Applied Mechanics and Environmental Protection Technologies National University of Civil Defence of Ukraine Chernyshevska str., 94, Kharkiv, Ukraine, 61023 E-mail: d.akass@ukr.net T. Shevchenko PhD, Associate Professor Department of Water and Wastewater Engineering O. M. Beketov National University of Urban Economy in Kharkiv Marshala Bazhanova str., 17, Kharkiv, Ukraine, 61002 E-mail: Tamara.Shevchenko@kname.edu.ua
Received date 23.07.2020 Copyright © 2020, S. Dushkin, T. Shevchenko
Accepted date 14.08.2020 This is an open access article under the CC BY license
Published date 25.08.2020 (http://creativecommons.org/licenses/by/4.0)
1. Introduction
The purification of water from the surface sources of water supply, used to provide drinking needs, from coarsely-dispersed, colloidal, and other contaminants, mostly employs the physical-chemical technology. This technology includes
the processes of settling and filtration, which are important elements of water supply systems in the preparation of drinking water [1, 2].
One of the most common methods of water purification from coarsely-dispersed and colloidal contaminants is the method of water treatment by coagulants. However, under
