CC BY
17
-□ □-
Побудовано математичну модель фть-траци з урахуванням кольматацшно-суфо-зшних процеыв в системi бюплато-фшьтра. Числовi розв'язки отриманог крайовог зада-чi знайдено методом сктченних елементiв з програмною реалiзацieю видтшдних алгорит-мiв в середовищ1 РгееРет+ + Проведено ряд числових ексnериментiв та здшснено аналiз впливу концентрации суспензи та кольмату-ючих частинок на процеси фшьтраци в бю-плато-фiльтрi
Ключовi слова: бюплато, проблема фшьтраци, кольматащя, суфозш, метод сктченних елементiв, РгееРет+ +
□-□
Построена математическая модель фильтрации с учетом кольматационно-суффози-онных процессов в системе биоплато-фильтра в двумерном случае. Числовые решения полученной краевой задачи найдены методом конечных элементов с программной реализацией соответствующих алгоритмов в среде РгееРет++. Проведен ряд численных экспериментов и осуществлен анализ влияния концентрации суспензии и кольматирующих частиц на процессы фильтрации в биоплато-фильтре на разных временных пластах
Ключевые слова: биоплато, проблема фильтрации, кольматация, суффозия, метод
конечных элементов, РгееРет+ + -□ □-
UDC 519.61/.64:627.05
|doi: 10.15587/1729-4061.2018.123455|
COMPUTER MODELING OF WATER CLEANING IN WETLAND TAKING INTO ACCOUNT OF SUFFOSION ANG COLMATATION
V. Moshynskyi
Doctor of Agricultural Sciences, Professor** Е-mail: rector@nuwm.edu.ua V. Fylypchuk Doctor of Technical Sciences, Рrofessor Department of Labor Protection and Life Safety** Е-mail: v.l.fylypchuk@nuwm.edu.ua N. Ivanchuk Senior Lecturer* E-mail: n.v.medvid@nuwm.edu.ua P. Martyniuk Doctor of Technical Sciences, Associate Professor* E-mail: p.m.martyniuk@nuwm.edu.ua *Department of Applied Mathematics** **National University of Water and Environmental Engineering Soborna str., 11, Rivne, Ukraine, 33028
1. Introduction
Bioplato that are facilities for the purification of household and industrial wastewater, as well as contaminated surface runoff, have been gaining popularity in recent years in different countries, specifically the United Kingdom, China, the Netherlands, Japan, Norway, Australia, etc. [1, 2]. The advantage of bioplato over a number of other water-treatment structures is the fact that they almost do not require electricity and chemical reagents, nor intensive operational maintenance [3, 4]. They also ensure the required quality of water purification from a wide range of pollutants of organic and mineral character [5].
The most common are the bioplato with an open water mirror (or similar structures, such as lagoons, ponds, channels, etc.) for water purification in countries with warm and moderate climate [6, 7]. The primary mechanism of water purification at a bioplato is the complex and interconnected life activities of heterotrophic and autotrophic organisms. These organisms develop in the thickness of filtration backfill, submerged in water, and on the surface of the root organs of plants in the form of immobilized bacterial medium. However, the open bioplato do not provide for the effective regulation of processes of mass exchange, gas saturation and aeration of water; the silt and sediment are not removed. This leads to the secondary water pollution and significantly reduces the effectiveness of purification from various, and particularly toxic, contaminants. In addition, in
areas with moderate climate the efficiency of open bioplato performance reduces in autumn-winter period by 30-40 % on average.
More effective are the closed bioplato of hydroponic type [8]. In them, water level is below the top level of the backfill, which has higher aquatic plants planted in it whose root system is constantly washed by water that moves vertically from top down or from bottom up. A special feature of such bioplato is the regulation of water quality through artificially created hydrobiocenosis, the characteristics of components of which form under the direct influence of plants or by using biopreparations.
However, known structures of bioplato undergo gradual colmatation of porous space in the filtering backfill and the bottom drainage with a biofilm and mineralized sediment, the accumulation of sludge in the ground part of the structure. The result is a decrease in the supply of oxygen to the root system of plants, which may reduce effectiveness of the work of structures, lead to sediment peptization and secondary pollution of water. Since the sediment in these bioplato is not removed, its accumulation gives rise to anaerobic biological processes occurring in the backfill, resulting in reduced sorption and detoxification of toxic impurities. Especially negative is feeding these bioplato with water with an elevated content of suspended solids particles, which greatly enhances negative processes that proceed in the thickness of a filtering backfill. Gradual colmatation of porous space in the filtering backfill leads to the process of suffusion, that is
©
extracting captured particles with purified water, which reduces efficiency of bioplato operation. To restore the work of bioplato, there should be a periodic pause in its operation to carry out complex and lengthy repair and restoration work related to washing and regeneration of the filtering backfill and drainage.
Therefore, it is an important task to explore suffu-sion-colmatation processes that occur in the thickness of a filtering backfill in a bioplato when filtering contaminated water, in terms of efficiency of its operation. Mathematical and computer modelling of such processes are associated with non-linear interdependence between parameters of the porous medium and parameters of the processes themselves. Specifically, suffusion and colmatation affect porosity of the backfill, which in turn is reflected on the magnitude of filtration coefficient. This impacts motion rate of the porous fluid, filtration-suffusion processes, and hence the performance of the bioplato. Therefore, ensuring the adequacy of research into these interconnected processes necessitates taking into consideration such non-linear mutual influences when constructing mathematical models of filtering in the bioplato.
2. Literature review and problem statement
In terms of the physics of a phenomenon, study of the processes of migration of insoluble particles in saturated porous media can be considered from two sides. The first is the problem of suffusion [9] (in the simplest case - carrying away movable particles from the skeleton of a porous medium [10]), related to problems of safety for hydro-technical, civil and industrial facilities [11]. Along with suffusion, the problems of colmatation are considered [12]. Processes of colmatation represent another side of the study [13]. These processes are physically opposite to the suffusion processes [14]. However, in terms of tools for mathematical and computer modeling, the specified processes are described by the same form of differential equations. The relations between parameters of these equations, assignment of the initial and boundary conditions for unknown functions, define mathematically the essence of the physical process. Specifically, when polluted water with a high concentration of suspended solid particles is fed to the soil massif, we observe predominance of colmatation processes over suffusion processes. If clean water is fed, suffusion processes will prevail. The importance of the process of migration of solid particles in porous media is emphasized in paper [15]. The observed processes are part of the comprehensive task on water purification from impurities by means of bioplato [7].
Mathematical model of suspension transport in porous media was constructed and applied to the processes of suffusion in paper [16]. Authors of work [17] applied for the specified processes a model of transfer (a system of first order differential equations), and a diffusion model when a degenerated equation of the parabolic type is obtained. Article [18] investigated iron nanoparticles transport in porous media and their application to clean up contaminated soil and groundwater. Results of laboratory study into colmata-tion and transfer of particles in a porous medium are given in paper [19].
Mathematical models of filtration-colmatation processes were constructed in some of the aforementioned scientific papers. Main attention is given to the processes of nanopar-ticles transport and the kinetics of colmatation processes.
However, the known studies failed to consider a dynamic change in the porosity of a porous medium as a result of its colmatation with particles and the non-linear impact of this process on the coefficient of filtration.
3. The aim and objectives of the study
The aim of present research is to quantify the influence of colmatation processes in filtering polluted water on the performance efficiency of the filter-bioplato by means of mathematical and computer modeling. This would make it possible to draw conclusions about the effectiveness or ineffectiveness of the work of bioplato of the applied design and to make a decision about the need (or lack of it) to introduce modifications to the engineering elements of the filter.
To accomplish the aim, the following tasks have been set:
- to account for the effect of concentration of colmatat-ing particles on a change in the porosity of a filtering backfill, which in turn would affect the modification of filtration equation;
- to take into consideration that the modified equation of filtration would affect via a convective component a change in the form of equation of mass transfer of suspension in the pores of the filter-bioplato;
- to build a mathematical model of filtering in the fil-ter-bioplato with respect to the nonlinearity of effect of suf-fusion-colmatation processes on the coefficient of filtration and porosity of the filtering backfill;
- find numerical solutions to the corresponding nonlinear boundary value problem, which describes a mathematical model to be constructed, applying the finite element method, and to implement software realization of the proposed algorithms using the software package FreeFem++.
4. Description of a model problem on the process of filtration in the filter-bioplato with respect to influence of suffusion and colmatation
Procedure of present research implies the construction of a mathematical model of the process of filtration in the filter-bioplato with respect to influence of the suffusion-col-matation processes. The effect of colmatating particles on filtration processes manifests itself in the nonlinear dependence of porosity coefficient on concentration. Specifically, for taking this dependence into consideration, one can apply the Kozeny-Karman's formulae [20]. Non-linearity of the appropriate boundary value problem, which describes the mathematical model constructed, requires the application of numerical methods, such as the method of finite elements. In order to automate the finite-element computation of model problems, we employed free specialized programming environment FreeFem++. The model problem for computer simulation was considered to be a bioplato with a length of 50 meters and a height of 2 meters with a gravel filling with fractions the size of 20 mm. Polluted water is fed through the upper drainage system made from perforated pipes and is discharged at the bottom of the filling through the perforated drainage deployed at the bottom of the bioplato. The bioplato was considered to be sufficiently wide to be able to apply the profile circuit of filtration and to reduce a spatial problem to the two-dimensional problem. Software implementation using the FreeFem++ makes it possible
to alter dimensions of the bioplato for model problems, to consider particularities of the backfill and porosity via the value of a filtration coefficient. We accepted the magnitude of 200 kg/m3 to be the maximum concentration of pollutants arriving at the bioplato.
As regards the equation for change s=s(X, t) - the concentration of colmatating (glued to the soil skeleton) particles, there are two fundamental approaches - linear and non-linear. According to [12], it is possible to use a linear equation of colmatation and suffusion
5. Construction of mathematical model of filtration with respect to colmatation-suffusion processes
Let a(X, t) be the porosity of soil, which is variable over time due to change in the concentration of colmatating particles; a0(X) is the porosity of the "skeleton" of the soil (at full absence of colmatating particles); s(X, t) is the mass concentration of colmatating particles (the mass of particles, which are associated with the skeleton of the soil and related to a unit volume of the porous environment); c is the concentration of suspension, which is filtered (the mass of suspended particles per a unit volume of a porous fluid); ps is density of a material of colmatating particles.
Let Vn be the volume of pores in separated volume V of the porous medium.
Then by definition we obtain
V V -V V V C(X, t ) = -п = п "s =1п.-Ijl,
V V V V
= °o(X ) -
Therefore
s(X, t) V ps
V
o(X, t ) = oo( X ) -
s( X, t )
ds
— = ac -p-s, dt F
or a non-linear equation of colmatation and suffusion
9s
— = a*-(s„ - s )-c -B-s,
dt
(4)
(5)
where a, a* are the coefficients of speed at which particles stick; p is a factor of speed of particle separation; smax is the maximum concentration of colmatating particles.
It was substantiated in [13] that it is more appropriate, rather than (5), to use a slightly modified equation. Specifically
ds 1 , s
(6)
where S=/0(a0-amitt), amin is the minimum porosity value (all passages are filled with colmatating particles); /0 is the average distance that a particle travels before being captured by a porous channel.
Paper [22] proposed, based on the study of filters, using the following equation for a change in the concentration of particles in the solid component:
9s
(1)
— = a-c -e-B(t) - s, dt FW
From [21], we obtain filtration equation
where s is a small parameter.
Let amin be the minimum porosity value. Then
v.(p, <c* (c,s, ^Hdc-J-p-1.
(2)
Here Vs is the volume of colmatating particles in separated volume V; pp=pp(c) is the density of porous fluid (suspension), which depends on the concentration of suspended solids particles c; kh=kh(c, s, a) is a filtration coefficient, which depends on the concentration of suspension and porosity; h is the head in a porous fluid.
Equation (2) is quasi-stationary relative to function h=h(X, t) and contains two unknown functions - c=c(X, t) and s=s(X, t), and, therefore, must be complemented with two more equations for the specified functions. To change the concentration of suspension c=c(X, t), we shall use a diffusion model. Specifically, according to [21], the modified equation of mass transfer takes the form
(7)
where as is the porosity of porous medium, which is composed of colmatating particles only. Then, from (7) and (1), we obtain
om¡„ =oo —
^ =Ps (1 )-°0-
(8)
Thus, the mathematical model of filtration with respect to the colmatation-suffusion processes is described by the following boundary value problem:
v.(p, «* («, „*»)=<,£.£-p. x .a
i - ± dpp
p p dc
эс
at"
1 - _l dpp p p dc
эс
at"
= v - ( A vc )-
u
v v
1 - dpp
Pp dc
- w / y
vc, dt
(3)
where Dc is the coefficient of dispersion of particles in the porous suspension; u=-kh(c, s, a)div h is the filtration rate of the porous suspension; w is the sedimentation rate of solids in suspension.
= v-(dcvc) -
9s
u
v
1 - ± dp p pdc
- w / y
ds
Vc-—, XеЦ dt
— = a-c-B-s, XеЦ dt P
u = -kh(c,s,o)vh, Xец
p
s
p
c(X ,0) = C0( X), s( X ,0) = s0(X ), X eQ = Qur,
h|ti = hi(X,t), X eyp (m,n)|y2 = 0,
c|Yi = Q(X,t), X ey1, (q,n)|T2 = 0.
Here, qc is the flow of suspended particles in the suspension; yt^y2=r is the bound of region Q, in which we examine the process; h^X, t), Ct(X, t) are the assigned functions.
6. Algorithm for finding an approximate solution to the obtained boundary value problem
In order to find an approximate solution to the stated boundary value problem, we employed the method of finite elements (FEM).
Weak statement of a boundary value problem implies the following. We shall multiply equation (2) by the test function
v1(X) eH0 = {v1(X): v1(X) e W21(Q), v1(X)| Yi = 0},
integrate over region Q, next we apply the Ostrogradsky-Gauss formula and obtain
JJ ( k,,ppVAVv1)d n-
,dPp dc
• —• v1
dc dt
dn-JJ
% A •
Ps dt
dn = 0.
Similarly, we shall multiply equation (3) and the initial condition for the concentration of suspension by the the test function
v 2(X) e H0 = { v 2(X) : v 2(X) e W21(Q), v 2(X) | y= o},
integrate over region Q, next we apply the Ostro-gradsky-Gauss formula and obtain
JJ DcVcVv2dQ-JJ
f f \ ^ c dp
u 1 v V
Ppdc
w / J
Vcv2dn+JJ— v2dn=
dt
=JJ
1-
V V
Pp dc
v2
dt
dn,
JJ c ( x, y,0 ) • v 2(X)dn = JJ c0 • v 2(X)d n.
For the concentration of particles in the solid component, we obtain
5« = (a • c'-1 -ps'-1 )-T + s'-1,
where t is the sampling time step, t=i-x, s'=s(X, t).
In order to find an approximated generalized solution to the obtained non-linear Cauchy problem, it is necessary to apply sampling over time. It is possible to use, for example, a predictor-corrector scheme [23], or fully implicit linearized difference scheme [24].
7. Results of the numerical experiment using the developed mathematical model of filtration with respect to the colmatation-suffusion processes
To perform a computer experiment using the developed mathematical model, we examined the filter-bioplato, which is described in the research procedure. Fig. 1 shows schematic of the region where the problem is being solved. At bound ri, there is a perforated pipe, which is used to feed polluted water to the surface of the bioplato. Although water is fed into the pipe under a certain pressure, when it is discharged out of the pipe it is freely filtered into the porous environment. Thus, we consider the pressure at which polluted water is fed to bound ri equal to the atmospheric pressure. Since pa=0, the boundary condition for heads on r1: h=2 m. At bound r2, there is a perforated pipe, out of which purified water is pumped out. The flows are accepted, respectively, at x=0 m, qmin=0.13 m/day, and at x=50 m, qmax=0.17 m/day. This enables performance efficiency of the bioplato at a volume of 160-200 m3 of water per day. Bounds r3, r4 are the impermeable boundaries.
Fig. 1. Schematic of the region where the problem is solved
At bound r2, the flow depends on the distance to the right edge:
I = gmin (x - X2 ) + ffmax (X - X1)
q XeF„ = + '
where q is the flow of fluid per unit length; x1=0 m; x2=50 m; at bound r1, c(X, t)=200 kg/m3. Other initial parameters are as follows:
o0 = 0.7, e0 = 2.3, ps = 2200 (kg/m3),
Dc = 9.5110-3 (m2/day), a = 0.1 (day-1),
P = 0.01 (day-1), w1 = 0 (m/day),
W2 = 3 (m/day,) C0 = 8 (kg/m3),
s0 = 0 (kg/m3).
Sampling time step t=30 days.
Fig. 2 shows difference in heads taking into account effect of the suspension concentration and the concentration of col-matating particles (filtration coefficient k=k0-(1+e0)/(1+e)x x((e)/(e0))3 [20]) and excluding it (filtration coefficient k=k0=100 m/day [25]), respectively.
Numerical experiments show that the maximum difference in heads considering and excluding effect of the suspension concentration and the concentration of colmatating particles amounts to 0.39 m in six months after launching the operation of a bioplato.
Fig. 2. Difference in heads in the bioplato backfill in 180 days
Fig. 3. Distribution diagram u2=—k(e)-3h/3y with respect to effect of the suspension concentration and the concentration of colmatating particles in 180 days
Fig. 4. Distribution diagram u2=—kdh/dy excluding effect of the suspension concentration and the concentration of colmatating particles in 180 days
8. Discussion of results of computer numerical experiments of filtration in the filter-bioplato
The research performed, although using a model problem as an example, showed that at the stage of designing it is important to investigate suffusion-colmatation processes that occur in the thickness of a filtering backfill in the bioplato during filtration of polluted water. These processes have a significant influence on the effectiveness of the filter-bio-plato. The tools of mathematical modeling and computer simulation, which were used, allow us to take into consideration non-linear mutual influences between parameters of a porous medium and parameters of the processes themselves. They make it possible to save significant resources and time required to carry out field experiments.
ith respect to that the maximum value of heads in the backfill of a bioplato is 2 m, one can see in Fig. 2 that the maximum relative difference reaches 19.5 %. Velocity diagrams of filtration along the 0y axis (component u2) demonstrate that over 180 days (Fig. 3, 4) maximum relative difference in the values of filtration velocity is about 25 %. In other words, the predicted performance efficiency of the bioplato will drop in six months by at least a quarter. Such an effect is explained by considering the nonlinear dependence of filtration coefficient on the concentration of colmatating par-
ticles. The increasing volume of particles that are captured by a backfill in the filter-bioplato leads to a decrease in the backfill porosity. This in turn results in the decreased filtration coefficient, hence a reduction in the rate and volume of fluid that is filtered in the bioplato. The built mathematical model accounts for the physical effects of dynamic change in porosity and the dependence of filtration coefficient on the concentration of pollutants, in contrast to known analogs. Taking these effects into consideration allows us to argue about improving the adequacy of the mathematical model to the examined physical processes. The negative impact of increasing concentration of colmatating particles on the performance efficiency of the filter could be quantified even excluding the mathematical model proposed. However, it is impossible to obtain specific forecast data without numerical study into constructed mathematical model and without conducting numerical experiments. Nonlinearity of the model built does not make it possible to argue about analytical solutions to the appropriate boundary value problem. Even more so for the would-be consideration of a three-dimensional problem. Mathematical modelling makes it possible to estimate the effect of impact of the examined factors, but it does not provide for the reduction of such a negative impact. Therefore, in order to maintain the operational performance of the bioplato, it is required to design engineering solutions that would reduce the impact of colmatation-suffusion processes on the effectiveness of the filter-bioplato.
9. Conclusions
1. We have constructed an improved mathematical model for the process of suspension filtration in the filter-biop-lato, which, in contrast to known analogues, accounts for a non-linear dependence of filtration coefficient on the concentration of colmatating particles, and a dynamic change in porosity in the process of filtration.
2. Numerical solutions to the corresponding nonlinear boundary value problem were derived applying the method of finite elements in the free software package FreeFem++. The use of FreeFem++ allowed us to partially automate programming implementation of algorithms for solving nonlinear boundary value problems for the systems of equations in partial derivatives. Specifically, one of the advantages of the employed programming environment is the automated coverage of regions with a grid of triangular finite elements, as well as the automatic reduction of the problem in weak statement to the system of linear algebraic equations.
3. Computer simulation revealed that, when taking into account the nonlinearity, the predicted performance efficiency of a bioplato reduces by 25 % in six months. This effect cannot be traced if the nonlinearity is neglected. Mathematical model allowed us to improve the adequacy of results in the examined physical processes; however, it does not make it possible to reduce the negative impact of colmatation. Therefore, such a task requires the development of engineering solutions to reduce the influence of colmatation-suffusion processes on the processes of filtration. This would require, in turn, accounting for these engineering solutions in the mathematical model, as well as further exploration of the efficiency of such solutions.
References
1. Healy M., Cawley A. M. Nutrient Processing Capacity of a Constructed Wetland in Western Ireland // Journal of Environment Quality. 2002. Vol. 31, Issue 5. P. 1739. doi: 10.2134/jeq2002.1739
2. Gleichman-Verheyc E. G., Putten W. H., Vander L. Alvalwaterzuvering met helofytenfilters, een haalbaarheidsstudie // Tijdschr. watervoorz. en. Afvalwater. 1992. Issue 3. P. 56-60.
3. Zavatskiy S. V., Kotelchuk L. S., Kotelchuk A. L. Bioengineering structures for low-waste sewage treatment // Construction, engineering systems and communications. 2012. Issue 1 (3). P. 57-63.
4. Purification of muddy waters in filtration regenerative wetland / Filipchuk V. L., Kuriluk M. S., Filipchuk L. V., Kuriluk O. M., Krylyuk V. M., Pochtar O. V. // Bulletin of the Engineering Academy of Ukraine. 2016. Issue 3. P. 150-155.
5. Water purification in filtration regenerative wetland / Filipchuk V. L., Bondar O. I., Kuriluk M. S., Ayaya A., Krivoshey P. P., Kuriluk O. M., Pochtar O. V. // Bulletin NUWEE. 2016. Issue 2 (74). P. 193-204.
6. Hosokova Y., Miyoshi E., Fukukawa K. Characteristics of the process of purification of coastal waters by reed beds // Rept. Partand Harbour. Res. 1991. Vol. 30, Issue 11. P. 206-257.
7. Blankenberg, A.-G. B., Braskerud B. C. «LIERDAMMEN» - a wetland testfield in Norway // Retention of nutrients, pesticides and sediments from a agriculture runoff: Diffuse Pollut. Conf. Dublin, 2003. P. 128-130.
8. Direnko A. A., Kotsar E. M. The use of higher aquatic plants in the practice of wastewater treatment and surface runoff // Plumbing, heating, air conditioning. 2006. Issue 4. P. 12-15.
9. Numerical simulation of backward erosion piping in heterogeneous fields / Liang Y., Yeh T.-C. J., Wang Y.-L., Liu M., Wang J., Hao Y. // Water Resources Research. 2017. Vol. 53, Issue 4. P. 3246-3261. doi: 10.1002/2017wr020425
10. Internal erosion in granular media: direct numerical simulations and energy interpretation / Sibille L., Lominé F., Poullain P., Sail Y., Marot D. // Hydrological Processes. 2014. Vol. 29, Issue 9. P. 2149-2163. doi: 10.1002/hyp.10351
11. Medvid N. V. Investigation of the influence of the conduit and the washout zone on the processes of filtration consolidation of the ground dam using the finite element method // Bulletin NUWEE. Series: Technical sciences. 2015. Issue 4 (72). P. 132-142.
12. Dobronravov O. O., Cremez V. S. Modeling of groundwater filtration taking into account the suffusion and colmatation // Problems of water supply, drainage and hydraulics. 2006. Issue 7. P. 141-146.
13. Kapranov U. I., Tropin N. M. The structure of the colmatation layer in the vicinity of the moving boundary // Applied Mechanics and Technical Physics. 2011. Vol. 52, Issue 6. P. 77-91.
14. Hydraulic Operating Conditions and Particle Concentration Effects on Physical Clogging of a Porous Medium / Alem A., Ahfir N.-D., Elkawafi A., Wang H. // Transport in Porous Media. 2014. Vol. 106, Issue 2. P. 303-321. doi: 10.1007/s11242-014-0402-8
15. Seetha N., Mohan Kumar M. S., Majid Hassanizadeh S. Modeling the co-transport of viruses and colloids in unsaturated porous media // Journal of Contaminant Hydrology. 2015. Vol. 181. P. 82-101. doi: 10.1016/j.jconhyd.2015.01.002
16. Chetti A., Benamar A., Hazzab A. Modeling of Particle Migration in Porous Media: Application to Soil Suffusion // Transport in Porous Media. 2016. Vol. 113, Issue 3. P. 591-606. doi: 10.1007/s11242-016-0714-y
17. Berres S., Bürger R., Wendland W. L. Mathematical Models for the Sedimentation of Suspensions // Lecture Notes in Applied and Computational Mechanics. 2006. Vol. 28. P. 7-44. doi: 10.1007/978-3-540-34961-7_1
18. Zhang M., He F., Zhao D., Hao X. Transport of stabilized iron nanoparticles in porous media: Effects of surface and solution chemistry and role of adsorption // Journal of Hazardous Materials. 2017. Vol. 322. P. 284-291. doi: 10.1016/j.jhazmat.2015.12.071
19. Chrysikopoulos C. V., Katzourakis V. E. Colloid particle size-dependent dispersivity // Water Resources Research. 2015. Vol. 51, Issue 6. P. 4668-4683. doi: 10.1002/2014wr016094
20. Francisca F. M., Glatstein D. A. Long term hydraulic conductivity of compacted soils permeated with landfill leachate // Applied Clay Science. 2010. Vol. 49, Issue 3. P. 187-193. doi: 10.1016/j.clay.2010.05.003
21. Gerus V. A., Kutya T. V., Martynyuk P. M. Generalization of the equation of filtration and heat and mass transfer in the case of suffusion processes // Mathematical and computer modeling. Series: Technical sciences. 2016. Issue 14. P. 48-63.
22. Safonyk A. P. Modelling the filtration processes of liquids from multicomponent contamination in the conditions of authentication of mass transfer coefficient // International Journal of Mathematical Models and Methods in Applied Sciences. 2015. Vol. 9. P. 189-192.
23. Sergienko I. V., Skopetskyi V. V., Deineka V. S. Mathematical modeling and research of processes in heterogeneous environments. Kyiv: Naukova Dumka, 1991. 432 p.
24. Vlasyuk A. P., Martynyuk P. M., Medvid N. V. Mathematical modeling of consolidation and filtration destruction of soils in the bases of hydraulic engineering and power structures. Rivne: NUWEE, 2017. 423 p.
25. Engineering geology and environmental protection: manual. Vinnitsa: VNTU, 2010. 262 p.
