CC BY
12
-□ □-
Mathematical modeling and computer simulation methods have been used to investigate the extent of influence exerted by bio-clogging on the dynamics of excess head scattering in the soil massif. To this end, the classical equation of filtration consolidation has been modified for the case of variable porosity resulting from changes in the biomass. The numerical solution to the constructed mathematical model in the form of a nonlinear boundary problem was derived by a finite-element method. Numerical experiments were carried out and their analysis was performed. Specifically, this paper shows the charts of pressure differences in the soil array when neglecting bio-clogging and when estimating the effects exerted by bio-clogging at specific points in time. The numerical experiments demonstrated that in two years after the onset of the consolidation process in the neighborhood of the lower limit of the examined soil mass with a thickness of 10 meters, excess heads fall from the initial value of 10 m to 4 m. The greatest impact from the clogging of pores by microorganisms is revealed in the neighborhood of an upper limit. At a depth of 1 m, at t=180 days, the pressure difference reaches 2.4 m. This is about 200 % of the pressure distribution without taking into account the effects of bio-clogging. Over time, the effect of bacteria on the distribution of pressures in the neighborhood of the upper boundary decreases. However, this effect extends to the entire soil mass, up to the lower limit. Thus, at t=540 days, at the lower limit, the effect of bio-clogging leads to that excess heads are 1.8 m greater than for the case of pure water filtration (a relative increase of about 80 %).
Bio-clogging processes are intensified as a result of the development of microorganisms when organic chemicals enter the porous environment. Therefore, from a practical point of view, studying them is especially relevant for household waste storage facilities and the stability of their soil bases. It is advisable to undertake research by using the methods of mathematical modeling and computer simulation
Keywords: excess heads, bio-clogging, organic waste, finite element method, filtration
consolidation -□ □-
UDC 519.61.64:627.05
|DOI: 10.15587/1729-4061.2021.230238]
CONSTRUCTING A MATHEMATICAL MODEL AND STUDYING NUMERICALLY THE EFFECT OF BIO-CLOGGING ON SOIL FILTRATION CONSOLIDATION
Natalia Ivanchuk
PhD*
E-mail: [email protected] Petro Martyniuk
Doctor of Technical Sciences, Professor* Е-mail: [email protected] Olga Michuta PhD, Associate Professor* E-mail: [email protected] Yevhenii Malanchuk Doctor of Technical Sciences, Professor Department of Automation, Electrical and Computer-Integrated
Technologies** Е-mail: [email protected] Hanna Shlikhta PhD, Associate Professor Department of Information and Communication Technologies and Computer Science Teaching Methods Rivne State University of the Humanities Stepana Bandery str., 12, Rivne, Ukraine, 33028 Е-mail: [email protected] *Department of Computer Science and Applied Mathematics** **National University of Water and Environmental Engineering Soborna str., 11, Rivne, Ukraine, 33028
Received date 03.03.2021 How to Cite: Ivanchuk, N, Martyniuk, P., Michuta, O, Malanchuk, Y, Shlikhta, H. (2021). Constructing a mathematical model
Accepted date 19.04.2021 and studying numerically the effect of bio-clogging on soil filtration consolidation. Eastern-European Journal of Enterprise Tech-
Published date 26.04.2021 nologies, 2 (10 (110)), 27-34. doi: https://doi.org/10.15587/1729-4061.2021.230238
1. Introduction
The problem of soil compaction has remained acute up to now and is likely to become more relevant due to the intensification of construction work. Applying external loads on the soil causes the emergence of excess heads in the porous liquid of the porous medium. Therefore, the compaction dynamics can be monitored through the dynamics of scattering of these excess heads. The presence of excess heads, for example, in the porous liquid of fully saturated soil, means the incompleteness of its (soil) consolidation. This is described as part of the so-called soil consolidation filtration theory, which was first reflected in the mathematical models by Karl Terzagi and Victor Florin [1]. The main provisions of the theory developed in the 1930s have been used to this day. Mon-
itoring, as well as the parameters of filtration consolidation of soils, is a mandatory element in the State building codes and rules [2, 3]. All construction works are regulated by such official documents equated to the status of laws.
The dynamics of dissipation of excess heads over time, and, consequently, the dynamics of the compaction process, are determined by the filtration characteristics of a porous environment. The higher the filtration coefficient that is included in Darcy's law, the greater the rate of dissipation of excess heads. As a result, the soil enters a stabilized state faster. If some factors affect the reduction of the filtration coefficient (for example, the clogging of pores, an increase in the viscosity of the pore fluid, etc.), that then increases the time of dissipation of excess heads and the time of stabilization of the soil base of the structure.
Applied loads can include not only civilian or industrial buildings but also waste storage (both liquid and solid). If there are organic wastes in storage facilities, the development of microorganisms becomes one of the factors influencing the soil porosity, and, therefore, the filtration coefficient. The consequence of the development of microorganisms is the bio-clogging of pores in a porous environment. The presence of microorganisms in the soil cannot lead to additional excess heads but indirectly, due to a decrease in the filtration coefficient, leads to a slowdown in the dissipation of such pressures. Their (excess heads) presence leads to uneven soil subsidence in the bases of civil and industrial structures. Examples of field observations of uneven soil subsidence as a result of its filtration consolidation are given in paper [4] (research in Busan, South Korea, from 2007 to 2019). Paper [5] noted that the observed uneven soil subsidence due to consolidation was greater in industrial areas than in areas with civil buildings.
The authors of [6] investigated, among other factors, the impact of the development of microbes as one of the factors of denitrification in agricultural production. The cited work shows that bio-processes in porous environments, including the development of microorganisms, are important; studying them is an urgent task. This is evidenced by results reported in [7]. Paper [8] shows that filtration consolidation is a process that can lead to uneven soil subsidence in the foundations of civil and industrial structures.
Bio-clogging processes, as a result of the development of microorganisms, are intensified when organic chemicals enter the porous environment. Therefore, from a practical point of view, studying them is especially relevant for household waste storage facilities and the stability of their soil bases. And it is advisable to carry them out using the methods of mathematical modeling and computer simulation.
But, in fact, the compaction of the porous environment is also associated with a change in porosity due to the imposed loads. That is, if the porous environment falls into the zone of influence of organic waste, then, unlike the classic case, at least two factors affect the dynamics of changes in porosity - applied loads and microorganisms. Investigating this case requires the modification of the classical approach to the study of soil consolidation by mathematical modeling and computer simulation methods.
2. Literature review and problem statement
The studies of bio-clogging processes, the development of bio-films, and their impact on the parameters of porous environments are reported in detail in [9, 10]. However, the cited works do not build mathematical models of the specified processes; only experimental studies are described. The authors of [11, 12] derived some equations that are components of the mathematical model of the influence of microorganisms on filtration processes in soil environments. However, neither a general mathematical model of those biological processes nor solutions are given in the cited works. Also, now there are quite a lot of works where these effects are taken into consideration in mathematical and computer models in porous environments [13, 14].
However, among such models, there are no mathematical models of the impact of bio-clogging on the processes of dissipation of excess heads in the soils arising from the imposed external loads. Similarly, current studies of filtra-
tion consolidation processes by mathematical modeling and computer simulation methods [15-17] disregard the effects of bio-clogging. However, bio-clogging is associated with a decrease in the volume of pores, which, in turn, can affect the dynamics of consolidation processes that are also associated with changes in the volume of pore space. For a more detailed study, refutation or confirmation of the hypothesis suggested, more detailed studies of interrelated bio-clogging and consolidation processes are required.
Paper [18] investigated the processes of soil consolidation under the influence of inorganic chemicals. The peculiarity of the presence of inorganic chemical solutions in the pores of clay soils (for them the problem of consolidation is acute) is pronounced osmotic flows and the dependence of the filtration coefficient on their (chemicals) concentration. The latter feature is associated with the presence in clays of a large amount of bound water and filtration of electrolytes as a solution of an inorganic chemical. In addition, work [19] studied the influence of salinity-desalination processes on the dynamics of dissipation of excess heads. However, the cited works do not consider the peculiarities of soil consolidation in the case of propagation of organic substances. In the presence of organic substances in the soil, the development of microorganisms is intensified (which is fundamentally impossible in the presence of inorganic chemicals) and, as a result, bio-clogging of porous material occurs.
Bio-clogging leads to a change in the hydrophysical parameters of the porous environment. Specifically, the conductivity of the environment decreases, which is mathematically reflected by the reduction in the filtration coefficient. Examples of such studies confirming the above statement are given in works [20, 21]. In work [20], for example, it was established that the value of the soil filtration coefficient decreased by an order of magnitude by day 540 from the beginning of the experiment, while in [21] - already by day 12. A detailed overview of scientific papers on the effect of bio-clogging on the value of the filtration coefficient is given in [22]. Also, work [22] reported the corresponding mathematical model and studied the effect of bio-clogging on the value of pressure jumps in the geo barrier. The impact on the excess heads in the entire soil massif has not been investigated there. In addition, a classical equation is used as a consolidation equation. At the same time, work [23] proposes a methodology for deriving the consolidation equation taking into consideration the dynamic change in various impact factors. This procedure is also used in the study of cleaning processes in bio plateau filters. Article [24] focuses on the mathematical modeling of dynamic filtration consolidation processes in saturated geoporous environments based on the biparabolic equation of evolution and its fractional-differential analog. Work [25] reports a solution to the problem of compaction of the soil with a threshold gradient, which can be expressed by elementary functions. However, such an impact factor as bio-clogging is not taken into consideration.
Thus, bio-clogging, as one of the consequences of the development of microorganisms, leads to a change in the hydrological parameters of the soil in the bases of organic waste storage facilities. However, available mathematical models of soil compaction do not take into consideration such changes. That is why it is justified to conduct studies that would cover the issues on the assessment of the impact of microorganisms in soil pores on the dynamics of dissipation of excess heads. To this end, it is necessary to improve the mathematical model of filtration consolidation by taking
into consideration the effect of bio-clogging in the pressure equation; to find numerical finite-element solutions to the corresponding nonlinear boundary problem for a system of differential equations; to conduct a series of numerical experiments and analyze them.
3. The aim and objectives of the study
The purpose of this study is to quantitatively assess the impact of the presence of microorganisms in soil pores on the dynamics of dissipation of excess heads in the soil massif, which is in a state of filtration consolidation. This would make it possible to draw conclusions about the degree of influence of bio-clogging on the compaction process.
To accomplish the aim, the following tasks have been set:
- to build a mathematical model of filtration consolidation taking into consideration the effect of bio-clogging;
- to find numerical finite-element solutions to the corresponding nonlinear boundary problem for a system of differential equations;
- to conduct a series of numerical experiments and analyze the results.
4. The study materials and methods
Our derivation of the modified filtration consolidation equation for the case of the effect of bio-clogging is based on the results reported in work [23]. In it, the generalized equation of filtration consolidation, taking into consideration the influence of technogenic factors, takes the following form
ence, we would consider only bio-clogging, that is, a change in the biomass of microorganisms in a porous environment.
Numerical methods, namely a finite-element method (FEM), were used to find an approximate solution to the corresponding nonlinear boundary problem. FEM makes it possible to sample the problem for spatial variables. As a result, we obtain a Cauchy problem for a system of first-order nonlinear differential equations relative to unknown functions that depend only on time. To find a solution to the Cauchy problem, time-based sampling techniques are used. In the numerical experiments, piecewise quadratic functions were used as basic FEM functions while a completely implicit linearized difference scheme was applied as a sampling scheme in time.
Glucose solution was taken as an organic substance for our numerical experiments. Data from the scientific literature were employed in the dependences of the filtration coefficient on porosity and to convert the biomass to volume.
5. Results of investigating the problem of filtration consolidation, taking into consideration the effect of bio-clogging, using a finite-element method
5. 1. Construction of a mathematical model of the compaction of saturated soil, taking into consideration the effect of bio-clogging
Considering the accepted assumptions and neglecting the creeping soil skeleton, equation (1) produces the following:
dh de dB \ d (,, sdh ^
Jadt + dB~dt + e'dX[ (e)aXJ'
(2)
de de d& -A de ds-—+---+ >----+
dt d& dt ^o ^
ds- dt
dpp dt
yPP
dp™. Pm dt
= (1 + e)(V-(kh (t, S )Vh)-V-F0sm ),
(1)
where:
- e - soil porosity coefficient;
- 0 - the sum of the main stresses in the skeleton of a porous environment;
- S=(0, s2, ..., sn) - impact factors influencing a change in porosity and filtration coefficient kh;
- Fosm - osmotic flow;
- Pp, pm - water density and solid, water-solute components density;
- h - excess heads in the porous liquid.
The following assumptions have been accepted. First, a one-dimensional case with one spatial variable x is considered. Second, suppose that a chemical solution of organic matter enters the soil pores. Such chemical solutions may be the result of impregnation or surface effluents from waste storage facilities. The possibility of such a scenario in the case of disposal of solid municipal waste is indicated, for example, in review paper [26]. The possibility of such filtration flows is also mentioned in work [27], which investigates the BIOLEACH model for assessing the production of biogas in real time.
For organic substances, osmotic flows are not characteristic; Fosm= 0. Also, solid water-solute components are absent; dpm/dt= 0. The concentration of the chemical organic solution is small; dpp/dt=0. Besides, among the factors of influ-
where a is the compression coefficient of the soil in compression dependence (for a one-dimensional case, e=-a0+eo); B - the biomass of microorganisms in the soil (reduced to the unit of the volume of the porous environment).
The process of filtration consolidation for the case of the distribution of organic chemicals is investigated in the soil layer with a total thickness of l (Fig. 1).
0
n
Fig. 1. Examined soil layer
A mathematical model would include the equation above (2) of filtration consolidation. Also, its components should be the equation for the transfer of organic chemicals in the porous liquid of a porous medium and the equation of the dynamics of biomass of bacteria in a porous environment based on the Mono equation [22]. As a result, we obtain the following boundary problem:
I
x
dh de dB , ,d (,dh ^
Jadt + dB~dt + e'dXt (e)aXJ'
, ^, i = i^s, W e(0;T]. dt dx
x EÛ, t > 0, h(X't)|= h0 (t)' t *° «(xt)L = (-k(e)f)| = 0- i * 0-
h(x,0) = h0 (x), x eû,
dc d ( _ dc 1 dc Umax _
w— = —I ra.D— - u--L-EaL B
dt dx I dx I dx Y
x EÛ, t > 0, c (x't )l x=0=t > °
1 -
B
B
,k + c
7 c
(X, t)|=-nDdx
c(x,0) = c0 (x), x EÛ,
= 0, t > 0,
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
dB „ — B
dt Krna*
1 --
B
Ac + c
x EÛ, t > 0,
B( x,0) = B0 ( x ), X EÛ,
Moreover, the first two satisfy the boundary conditions of the first kind, as well as the functions h(x, t), c(x, t), respectively. Here, T>0.
Take {si(x); s2(x); s3(x)} belongs to H0. Multiply equation (3) and the initial condition (6) by s1(x), integrating them in the segment [0, l]; we obtain
j Ja ( x )dx + j dB ^i (x )dx +
tdB dt )dhdS( x )
dx dx
+j(l+e) k (e)^^^-dx = 0,
Vt e(0;T ], (13)
j h( x,0)s1 ( x ) dx = j h, ( x ) s1 ( x ) dx, Vt e(0;T]. (14)
Similarly, for the concentration of the organic chemicals and biomass, we obtain
in —s2 ( x )dx + inD——dx I dt Jn dx dx
+ i u—s2 ( x )dx +
2
dx
Vt e(0;T ], (15)
(11) (12)
J Y
B-s2 ( x ) dx = 0,
k + c
where:
- h - excess head;
- k - filtration coefficient that depends on the porosity coefficient e;
- u - filtration speed;
- n - soil porosity;
- qc - a flow of chemicals;
- Idec - a biomass decay coefficient;
- c - the concentration of the porous solution of organic chemical;
- D - an organic matter diffusion coefficient;
- |max - the maximum level of biomass growth;
- B - the concentration of bacteria biomass in a porous environment;
- Bmax - the maximum possible concentration of biomass of bacteria in a porous environment;
- kc - a parameter in the equation of the biomass dynamics at which the level of biomass growth equals half the maximum;
- Y- bacterial output.
5. 2. Numerical solution to a filtration consolidation problem taking into consideration the effect of bio-clogging using a finite-element method
Similarly to [28], let H0 be the space of vector-functions (si(x); s2(x); s3(x)}, whose each component in the interval (0;1) belongs to the Sobolev space W2 (Q), and the first two acquire zero values at the ends of the segment [0,l], where the functions h(x, t) and c(x, t) are given boundary conditions of the first kind, respectively. Let H be the space of functions (ui(x, t); u2(x, t); u3(x, t)}, whose each component is integrated with a square along with its first derivatives in the interval (0; 1)
J c ( x ,0) s2 ( x )dx = Jc0 (x)s2 (x)dx, Vt e(0;T], (16)
jf^s ( * )dx -
dt
-J
/ „ \
1 --
B
B„
V V max
k + c
Bs3 ( x )dx = 0,
Vt e(0;T ],
i i
J B ( x,0) s3 ( x ) dx = J B0 ( x ) s3 ( x )dx, Vi e (0;T ].
(17)
(18)
Definition 1. The function {h(x, t); c(x, t); B(x, t)} that belongs to H, which, for any {s1(x); s2(x); s3(x)}, satisfies the integrated ratios (13) to (18), is termed a generalized solution to boundary problem (3) to (12).
An approximate generalized solution to boundary problem (3) to (12) is to be found in the following form:
N N
h(X, t) = X h (t(X), C (X, t) = X C (t№2 (X),
B ( *, t ) = £ Bi (t № 3 ( X )
(19)
where h;(t), c;(t), B;(t) are the unknown coefficients that depend only on time. If we consider (19) as an approximate finite-element solution, then the functions tyn(x), 9i2(x), 9;3(x) are the polynomial basis functions with a finite carrier. Similarly to [22], in order to simplify fur-
ther considerations, and in accordance with the numerical experiments, we shall assume that for the approximate search for h(x, t), c(x, t), and B(x, t) the same grid of finite elements is used. In this case, tyn(x)=^i2(x)=^i3(x), and, in order to avoid double indexing, we denote the specified basis functions through q>t(x), i = 1, N.
Next, given the weak statement of (13) to (18) of problem (3) to (12), taking into consideration (19), we obtain
f+M"'§+L'(B)H {l >" °
m; ■ h (0)=F;,
M2 (B) -dC+L2 (H,B) C(t) = 0,
M2 C(0) = f2,
dB
M3-—+l (C, b).B (t )=0,
M3 ■ B(0) = F3, Where
H = (h(t))^ c = (ci(t))l B = (b(t))^ h(0)= (*(0))!=!, C=(c<(0))l, B(0)= (0))^,
^=«L- L=(rL- 1=(rL
F=№ ^ = h*
i _
m'w = irn m dx, k = 1,3, M13 = (m- , ,
i J TsY; > II 13 \ i ¡i -=i I
0
/ l = Jy«^ $ ;dx, 4(1) = J h¡(?iàx,
0 0
j = i «+•) » W dx,
„(13) -
1 = f.^^m m dx,
' f dB^'
l l m(2) = J$;dx, fP = J^^
0 0
j = dx + ^ è dx +
J v
1 -
B
B„
—è è. dx,
k + cjl
mf* j dx, ;=} s^dx,
/(3)=— 9 J
1 -
B
B„
+ c
(20) (21) (22)
(23)
(24)
(25)
The system of equations (20) to (25) is a Cauchy problem for the system of nonlinear differential equations of the first order. Finding a solution to it also requires the use of appropriate sampling schemes, which are described in detail in works [22, 29].
5. 3. Results of the numerical experiments to solve the problem of filtration consolidation taking into consideration the effect of bio-clogging
In the model problem, we considered a layer of soil with a thickness of l=10 m. The step on the variable x was 0.02 m. The step on the time t=1 day, a=5.12-10-S m2/H is the coefficient of soil compressibility, y =104 H/m3 is the specific weight of the pore solution. The initial distribution of heads is h0(x)=10 m, which corresponds to the application to the soil surface of the corresponding load - wastes that are a source of organic matter. At the upper limit, an unhindered outflow of pore fluid is provided, and there is no drainage at the lower limit. Silt Loam [30] was considered as the soil at: ¿0=0.108 m/day, n0=0.45, where the index "0" means the initial values. Also, e0=n0/(1-n0).
Glucose solution [22] was considered an organic solution. The initial distribution of glucose concentration in the porous water is c0(x)=1 mM. The glucose diffusion coefficient is D=0.518440-4 m2/day.
The following dependence from [31] is used for the dependence of filtration coefficient on porosity:
k_
19/6
where nB is the volume of biomass per unit of the volume of the porous medium. To convert the concentration of biomass to volume, we accept the assumption that the biovolume is 80 % water, and the remaining 20 % of the dry mass is 50 % carbon (bio-carbon) [32]. In other words,
B
0.8pw + 0.2 A
where pw is the density of water; pc - density of bio-carbon. In the numerical experiments, pc=100 kg/m3. The parameters of the biomass kinetics are as follows: |max=0.2 day-1, Bmax=100 kg/m3, |dec=0.01 day-1, kc=0.06 mM, F=50. Given the dependence e=n/(1-n), we obtain
de _ de dn dB ~dn dB'
1
1
(1 - n) 0.8p + 0.2 ^ Kn 2
where n=n0-nB.
Piecewise quadratic functions were used as the FEM basis functions. The results of our numerical experiments are shown in Fig. 2-4.
Fig. 2-4 shows, as an example, the results from the numerical experiments at t=180 days, 360 days, 540 days, and 720 days.
B
180 days
_|j360 days .540 days .720 days
x,m
Fig. 2. Distribution of excess heads under the influence of bio-clogging 'Ah.m1
----1--I- -I- r - r
1.5--
Fig. 3. Distribution of excess head difference considering and disregarding the effect of bio-clogging
f B, kg/m3
100- —^ — 1— - H - 4 -4-1- -1-4 - 4
1 1 ~ 1 ~ 1 T 1 ~r~r 1 1 "I" 4 1 1 ~ T 1
. _ _ 1 i- -1— 1 - -1 - 1 4 1 1 -4-1- 1 1 -1-4 1 - 4
50- ■-- ■ -1— - -i - 4 -4-1- -1-4 - 4
■ ~ T ~ n ~ T ~v~r 4 ~ T
. _ i -L- - J - 4 -4-1- -1-4 _ 4
10- ■ - 4 —1 - 4 -4-4 -1-4 - 4
0 2 4 6 8 10
Fig. 4. The distribution of bacteria biomass at a time of 720 days
6. Discussion of results of solving a filtration consolidation problem taking into consideration the effect of bio-clogging
Fig. 2 shows that in two years (720 days) after the onset of the consolidation process in the lower boundary of the soil massif, excess heads fall from the initial value of 10 m to 4 m. The effect of bio-clogging itself is best traced from Fig. 3. It shows the charts of excess head differences in case of taking into consideration the effects of bio-clogging and without taking them into consideration, respectively. Therefore, positive values of head differences in Fig. 3 mean that when bio-clogging is taken into consideration, the heads dissipate slower than in case of neglecting such a phenomenon. Slow-
ing down the head dissipation physically means unfinished consolidation of the soil mass. The greatest impact exerted by the clogging of pores by microorganisms is in the neighborhood of an upper limit. This is logical because nutrients enter the soil from the upper limit and it is at a shallow depth (Fig. 4) that intensive development of bacteria biomass occurs. Fig. 3 shows that at a depth of 1 m, at £=180 days, the head difference reaches 2.4 m. That is about 200 % of the head distribution without taking into account the effect of bio-clogging. The insignificant depth of penetration of bacteria biomass is explained by the nature of the spread of organic substances. Fig. 2 shows that the filtration rate of the pore fluid would be directed towards the upper limit (in the direction of the head anti-gradient). Therefore, the main role in the spread of organic matter in the soil would belong to the diffusion component. Given the low diffusion coefficient of the glucose solution, organic matter would penetrate the soil base only by 2 meters in two years. And that affects the dynamics of the development of microorganisms, the main role for which belongs to the presence of nutrients. Over time, the effect of bacteria on the distribution of heads in the neighborhood of the upper boundary decreases. However, this effect extends to the entire soil mass, up to the lower limit (charts in Fig. 3 at £=540 days, £=720 days). Thus, at £=540 days, at the lower limit, the effect of bio-clogging leads to the fact that excess heads are 1.8 m greater than for the case of pure water filtration (a relative increase of about 80 %).
From a physical point of view, bio-clogging, which means a decrease in the volume of pores, should lead to a slowdown in the dissipation of excess heads (after all, the filtration coefficient also decreases). Such conclusions are confirmed by the results from the numerical experiments (for example, Fig. 3). The distribution of the head difference in Fig. 3 actually confirms the physics of the process. However, to quantify the impact, one needs to find solutions to the refined mathematical model. An example of the numerical experiment shows that exposure to bio-clogging causes changes in heads, by 80-200 %, compared to pure water filtration. Despite the spread of microorganisms only to a depth of 2 meters, as shown by charts in Fig. 3, the effect of bio-clogging in 1.5-2 years extends to heads to a depth of 10 meters.
However, for the possibility of carrying out such calculations for actual conditions, the hydrological characteristics of soils and the parameters of microorganism growth dynamics should be known. This study could be advanced by increasing the dimensionality of the problem, improving the models of the development of microorganisms, and pro-
x,m
posing specific engineering solutions to reduce the effect of bio-clogging on the excess heads in the soil.
7. Conclusions
1. We have improved the mathematical model of filtration consolidation during the filtration of organic matter in the soil array, taking into consideration the effect of bio-clogging. At the same time, the equation of filtration consolidation has been modified taking into consideration the impact exerted on a change in porosity by the dynamics of biomass development. The Mono equation was used as an equation for the dynamics of bacteria biomass in a porous environment.
2. The numerical solutions to the corresponding nonlinear boundary problem were found by using a finite-element method. Computer simulation has shown that the presence of microorganisms in the soil's pores significantly affects the dynamics of dissipation of excess heads. Specifically, bio-clogging (based on the solution to the model problem) leads to a relative slowdown in the dissipation of heads by 80-200 % compared to the case of neglecting the microorganisms. Such influences should be taken into consideration in the studies of landslides that are exposed to organic chem-
icals because the presence of excess heads could lead to the intensification of shear processes. Estimating the extent of such impacts on landslide events is the essence of a separate task of research in the mechanics of porous environments.
3. The results of our experiments showed the adequacy of the improved mathematical model to the processes studied - a slowdown in the dissipation of pressures under the influence of bio-clogging. In practice, our research could be used in the preparation of design documentation and structural solutions for organic waste storage facilities. The development of microorganisms and the process of bio-clogging of pores of a porous environment with biomass leads to a slowdown in the dissipation of excess heads. After all, excess heads are a negative phenomenon in general, which leads to disrupting the stability of the soil foundations of buildings.
Acknowledgments
This research is funded within the framework of the study "Mathematical modeling and computer simulation of anthropogenic controlled processes in porous media with barriers under conditions of identification", State Registration No. 0120U102055.
References
1. Zaretskii, Yu. K. (1972). Theory of soil consolidation. Israel Program for Scientific Translation.
2. DBN V.2.4-3:2010 (2010). Derzhavni budivelni normy Ukrainy. Hidrotekhnichni, enerhetychni ta melioratyvni systemy i sporudy, pidzemni hirnychi vyrobky. Hidrolohichni sporudy. Osnovni polozhennia. Kyiv: Minrehionbud Ukrainy. Available at: https:// dbn.co.ua/load/normativy/dbn/1-1-0-802
3. DBN V.2.1-10-2009 (2009). Derzhavni budivelni normy Ukrainy. Obiekty budivnytstva ta promyslova produktsiya budivelnoho pryznachennia. Osnovy ta fundamenty budynkiv i sporud. Osnovy ta fundamenty sporud. Osnovni polozhennia proektuvannia. Kyiv: Minrehionbud Ukrainy. Available at: https://dbn.co.ua/load/normativy/dbn/dbn_v21_10_2009/1-1-0-319
4. Park, S., Hong, S. (2021). Nonlinear Modeling of Subsidence From a Decade of InSAR Time Series. Geophysical Research Letters, 48 (3). doi: http://doi.org/10.1029/2020gl090970
5. Widada, S., Zainuri, M., Yulianto, G., Satriadi, A., Wijaya, Y. (2020). Estimation of Land Subsidence Using Sentinel Image Analysis and Its Relation to Subsurface Lithology Based on Resistivity Data in the Coastal Area of Semarang City, Indonesia. Journal of Ecological Engineering, 21 (8), 47-56. doi: http://doi.org/10.12911/22998993/127394
6. 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: http://doi.org/10.1029/2018wr023202
7. Moshynsky, V., Riabova, O. (2013). Approaches to Aquatic Ecosystems Organic Energy Assessment and Modelling. Black Sea Energy Resource Development and Hydrogen Energy Problems. NATO Science for Peace and Security Series C: Environmental Security. Dordrecht: Springer, 125-135. doi: http://doi.org/10.1007/978-94-007-6152-0_12
8. Michuta, O. R., Vlasyuk, A. P., Martinyuk, P. N. (2013). Influence of chemical erosion on filtration consolidation of saline soils in nonisothermal conditions. Matematicheskoe modelirovanie, 25 (2), 3-18. Available at: http://www.mathnet.ru/php/archive.phtml ?wshow=paper&jrnid=mm&paperid=3327&option_lang=rus
9. Gui, R., Pan, Y., Ding, D., Liu, Y., Zhang, Z. (2018). Experimental Study on Bioclogging in Porous Media during the Radioactive Effluent Percolation. Advances in Civil Engineering, 2018, 1-6. doi: http://doi.org/10.1155/2018/9671371
10. Wang, Y., Huo, M., Li, Q., Fan, W., Yang, J., Cui, X. (2018). Comparison of clogging induced by organic and inorganic suspended particles in a porous medium: implications for choosing physical clogging indicators. Journal of Soils and Sediments, 18 (9), 2980-2994. doi: http://doi.org/10.1007/s11368-018-1967-6
11. Thullner, M., Regnier, P. (2019). Microbial Controls on the Biogeochemical Dynamics in the Subsurface. Reviews in Mineralogy and Geochemistry, 85 (1), 265-302. doi: http://doi.org/10.2138/rmg.2019.85.9
12. Glatstein, D. A., Montoro, M. A., Carro Pérez, M. E., Francisca, F. M. (2017). Hydraulic, Chemical and Biological Coupling on Heavy Metals Transport Through Landfills Liners. The Journal of Solid Waste Technology and Management, 43 (3), 261-269. doi: http://doi.org/10.5276/jswt.2017.261
13. Mohanadhas, B., Kumar, G. S. (2019). Numerical Experiments on Fate and Transport of Benzene with Biological Clogging in Vadoze Zone. Environmental Processes, 6 (4), 841-858. doi: http://doi.org/10.1007/s40710-019-00402-w
14. Lopez-Peña, L. A., Meulenbroek, B., Vermolen, F. (2019). A network model for the biofilm growth in porous media and its effects on permeability and porosity. Computing and Visualization in Science, 21 (1-6), 11-22. doi: http://doi.org/10.1007/s00791-019-00316-y
15. Bohaienko, V., Bulavatsky, V. (2020). Fractional-Fractal Modeling of Filtration-Consolidation Processes in Saline Saturated Soils. Fractal and Fractional, 4 (4), 59. doi: http://doi.org/10.3390/fractalfract4040059
16. Józefiak, K., Zbiciak, A., Brzeziñski, K., Maslakowski, M. (2021). A Novel Approach to the Analysis of the Soil Consolidation Problem by Using Non-Classical Rheological Schemes. Applied Sciences, 11 (5), 1980. doi: http://doi.org/10.3390/app11051980
17. Tian, Y., Wu, W., Jiang, G., El Naggar, M. H., Mei, G., Xu, M., Liang, R. (2020). One-dimensional consolidation of soil under multistage load based on continuous drainage boundary. International Journal for Numerical and Analytical Methods in Geomechanics, 44 (8), 1170-1183. doi: http://doi.org/10.1002/nag.3055
18. Vlasyuk, A. P., Martynyuk, P. M. (2010). Numerical solution of three-dimensional problems of filtration consolidation with regard for the influence of technogenic factors by the method of radial basis functions. Journal of Mathematical Sciences, 171 (5), 632-648. doi: http://doi.org/10.1007/s10958-010-0163-z
19. 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: http:// doi.org/10.1007/s10958-009-9518-8
20. Glatstein, D. A., Francisca, F. M. (2014). Hydraulic conductivity of compacted soils controlled by microbial activity. Environmental Technology, 35 (15), 1886-1892. doi: http://doi.org/10.1080/09593330.2014.885583
21. 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: http://doi.org/10.1016/j.jenvman.2018.05.039
22. Ulyanchuk-Martyniuk, O., Michuta, O., Ivanchuk, N. (2020). Biocolmation and the finite element modeling of its influence on changes in the head drop in a geobarrier. Eastern-European Journal of Enterprise Technologies, 4 (10 (106)), 18-26. doi: http:// doi.org/10.15587/1729-4061.2020.210044
23. Gerus, V. A., Martyniuk, P. M. (2015). Generalization of the soil consolidation equation taking into account the influence of physicochemical factors .Bulletin of VN Karazin Kharkiv National University. Series: Mathematical modeling. Information Technology. Automated control systems, 27, 41-52. Available at: http://nbuv.gov.ua/UJRN/VKhIMAM_2015_27_7
24. Bulavatsky, V. M., Bohaienko, V. O. (2020). Some Consolidation Dynamics Problems within the Framework of the Biparabolic Mathematical Model and its Fractional-Differential Analog. Cybernetics and Systems Analysis, 56 (5), 770-783. doi: http:// doi.org/10.1007/s10559-020-00298-7
25. Wang, H.-X., Xu, W., Zhang, Y.-Y., Sun, D.-A. (2021). Simplified solution to one-dimensional consolidation with threshold gradient. Computers and Geotechnics, 131, 103943. doi: http://doi.org/10.1016/jxompgeo.2020.103943
26. Araujo, F., Fantucci, H., Nunes, E., Santos, R. M. (2020). Geochemical Modeling Applied in Waste Disposal, and Its Relevance for Municipal Solid Waste Management. Minerals, 10 (10), 846. doi: http://doi.org/10.3390/min10100846
27. Rodrigo-Ilarri, J., Rodrigo-Clavero, M.-E., Cassiraga, E. (2020). BIOLEACH: A New Decision Support Model for the Real-Time Management of Municipal Solid Waste Bioreactor Landfills. International Journal of Environmental Research and Public Health, 17 (5), 1675. doi: http://doi.org/10.3390/ijerph17051675
28. Sergienko, I. V., Skopetskiy, V. V., Deyneka, V. S. (1991). Matematicheskoe modelirovanie i issledovanie protsessov v neodnorodnykh sredakh. Kyiv: Naukova dumka, 432.
29. Ulianchuk-Martyniuk, O. V. (2020). Numerical simulation of the effect of semipermeable 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: http://doi.org/10.32523/2306-6172-2020-8-1-91-104
30. Simunek, J., van Genuchten, M. T., Sejna, M. (2016). Recent Developments and Applications of the HYDRUS Computer Software Packages. Vadose Zone Journal, 15 (7). doi: http://doi.org/10.2136/vzj2016.04.0033
31. 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: http://doi.org/10.1111/j.1745-6584.1996.tb02088.x
32. Manfred, B., Jaap, B., Klaus, M., Rolf, M. (2006). Enumeration and Biovolume Determination od Microbial Cells. Microbiological Methods for Assessing Soil Quality. CABI Publishing, 93-113. doi: https://doi.org/10.1079/9780851990989.0093
