Next Article in Journal
Exact Solution of Two-Dimensional Fractional Partial Differential Equations
Next Article in Special Issue
Comb Model with Non-Static Stochastic Resetting and Anomalous Diffusion
Previous Article in Journal
Fractional Kinetic Equations Associated with Incomplete I-Functions
Previous Article in Special Issue
Parametric Identification of Nonlinear Fractional Hammerstein Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Mathematical Model for the Description of Anomalous Migration of Soluble Substances in Vertical Filtration Flow

by
Vsevolod Bohaienko
* and
Volodymyr Bulavatsky
VM Glushkov Institute of Cybernetics of NAS of Ukraine, 03187 Kyiv, Ukraine
*
Author to whom correspondence should be addressed.
Fractal Fract. 2020, 4(2), 20; https://doi.org/10.3390/fractalfract4020020
Submission received: 11 April 2020 / Revised: 30 April 2020 / Accepted: 30 April 2020 / Published: 6 May 2020
(This article belongs to the Special Issue The Craft of Fractional Modelling in Science and Engineering III)

Abstract

:
Since the use of the fractional-differential mathematical model of anomalous geomigration process based on the MIM (mobile–immoble media) approach in engineering practice significantly complicates simulations, a corresponding simplified mathematical model is constructed. For this model, we state a two-dimensional initial-boundary value problem of convective diffusion of soluble substances under the conditions of vertical steady-state filtration of groundwater with free surface from a reservoir to a coastal drain. To simplify the domain of simulation, we use the technique of transition into the domain of complex flow potential through a conformal mapping. In the case of averaging filtration velocity over the domain of complex flow potential, an analytical solution of the considered problem is obtained. In the general case of a variable filtration velocity, an algorithm has been developed to obtain numerical solutions. The results of process simulation using the presented algorithm shows that the constructed mathematical model can be efficiently used to simplify and accelerate modeling process.

1. Introduction

We study the problem of mathematical modeling of anomalous dynamics of soluble substances’ convective diffusion under the conditions of vertical steady-state filtration of groundwater. Such problems arise while solving many problems related to the protection of water resources from pollution by industrial and domestic wastewater, as well as due to the need for desalination, leaching of soils during land reclamation, etc. It should be noted that an extensive amount of literature is devoted to these problems considering them within the framework of classical models of mass transfer (e.g., [1,2,3,4,5]). Empirical studies of these processes along with their investigation using differential models are presented in [6,7,8]. Soft computing techniques are also but rarely used for their analysis and prediction [9]. Regarding the issues of modeling anomalous migration of soluble substances in underground filtration flows, several problems have been posed and solved using the fractional-differential approach ([10,11,12,13,14]) that allows, in particular, taking into account the effects of memory and spatial correlations [15,16,17].
We consider the mathematical model of geomigration process based on the MIM (mobile–immobile media) approach [18,19,20]. If we take into account immobilization (entrapment of particles by soil’s skeleton or their penetration into the volume of a bound fluid), then the equation of convective diffusion takes the form [18,19,20]
C t + β C I t = L ( C ) ,
where C , C I are the volume concentrations of particles in mobile and immobile phases, β = θ I / θ , θ , θ I are the porosities in mobile and immobile zones, L ( C ) = d Δ C υ C , υ is the filtration velocity, d is the convective diffusion coefficient, ∇ is the Hamilton operator, and Δ is the Laplace operator with respect to the geometric variables.
According to Shumer et al. [20], the dynamics of particles’ outflow in the immobile phase can be described by the equation
C I t = C D t γ C
where C D t γ is the operator of Caputo–Gerasimov fractional derivative of the order γ ( 0 < γ < 1 ) with respect to the variable t [21,22]. Then, considering Equation (2), Equation (1) of the MIM model takes the form [18,19,20]
C t + β C D t γ C = L ( C )
with the total concentration calculated as C t o t = θ C + θ I C I . Considering that medium’s memory-related properties manifest themselves in the process of migration in mobile phase, Equation (3) takes the form
C D t α C + β C D t γ C = L ( C )
where C D t α is the operator of Caputo–Gerasimov fractional derivative of the order α ( 0 < α < 1 ) with respect to the variable t.
Let us note, that Equation (4) is formally similar to the equation of water exchange in mobile–immobile zones of swelling soils obtained in [23].
The use of the mathematical model of anomalous migration process that is based on the integro-differential Equation (4) in engineering practice with the aim of developing reliable methods for predicting process dynamics is associated with mathematical difficulties and complicates the simulation process. In some cases, one of the effective approaches to accelerate and simplify estimative engineering calculations when studying migration dynamics under the conditions of geofiltration consists in the simplification of the original mathematical model by an appropriate approximation of fractional derivatives and transition to a new, simplified mathematical model based on the classical convective diffusion equation. Thus, approximating the fractional derivatives in Equation (4) using the relation [24]
C D t α u ( t ) α u ( t ) + ( 1 α ) ( u ( t ) u ( 0 ) )
we obtain from Equation (4) the equation of a simplified mathematical model of the considered process of convective diffusion with particles’ immobilization in the form
a C t ( x , y , t ) = L ( C ) b ( C ( x , y , t ) C ( x , y , 0 ) )
where a = α + β γ , b = 1 α + β ( 1 γ ) .
Further, within the framework of the simplified mathematical model based on Equation (6), we state a two-dimensional initial-boundary value problem of convective diffusion of soluble substances under the conditions of vertical steady-state filtration of groundwater with a free surface from a reservoir into a coastal drain. In the case of averaged filtration velocity, we obtain an analytical solution of the considered problem. In the general case of variable filtration velocity, we present an algorithm for obtaining its numerical solution and the results of simulations.

2. Initial-Boundary Value Problem

We consider the problem of convective diffusion modeling with vertical steady-state filtration of groundwater from a reservoir to a buried semi-infinite drainage canal.
The filtration scheme corresponding to this problem has been studied quite well. For example, Lavryk et al. [5] provided a closed-form solution of the initial-boundary value problem obtained by conformal mapping of the filtration domain G z onto the domain G ω in the plane of complex flow potential ω = φ + i ψ (where φ is the potential of filtration velocity and ψ is the flow function). As a result of filtration problem’s solution, the characteristic flow function z = f ( ω ) was found and flow velocity field was determined. The corresponding relations are given in [5].
Within the framework of the original mathematical model based on Equation (4), the modeling of fractional-differential dynamics of convective diffusion in accordance with the considered filtration scheme can be performed finding in the domain G z × ( 0 , ) the solutions of the equation
C D t α C ( x , y , t ) + β C D t γ C ( x , y , t ) = d ( 2 C x 2 + 2 C y 2 ) υ x C x υ y C y ,
( 0 < α , γ < 1 )
with the following initial and boundary conditions:
C | y = 0 = C 1 , C n | γ 1 , γ 2 = 0 , C y | y = H = 0 ,
C ( x , y , 0 ) = C 0 ( x , y )
where C 1 is the given concentration of soluble substances on the filtration inflow y = 0 , C 0 ( x , y ) is the given initial distribution function, n is the outer normal to the corresponding curve, γ 1 , γ 2 are the depression curves, H is the depth of the drainage canal, and υ x , υ y are the projection of the filtration velocity vector on the axes O x and O y .
Since the filtration domain G z is a non-canonical domain of complex configuration, an efficient way to solve initial-boundary value problems of this type is to make transition to new variables ( φ , ψ ) —the points of geometrically simpler domain of complex flow potential, which is [5] for the considered problem the rectangle G ω with sides φ 0 , Q ( φ 0 = κ H where κ is the soil’s filtration coefficient and Q is the full filtration rate). Then, the problem in Equations (7)–(9) can be formulated in the domain of complex flow potential G ω in the form
( C D t α + β C D t γ ) C ( φ , ψ , t ) = υ 2 ( φ , ψ ) ( d Δ C ( φ , ψ , t ) C φ ) ,
( 0 < α , γ < 1 )
C | φ = 0 = C 1 , C ψ | ψ = 0 , ψ = Q = 0 , C φ | φ = φ 0 = 0 ,
C ( φ , ψ , 0 ) = C 0 ( φ , ψ )
where
υ 2 = υ x 2 + υ y 2 = ( φ x ) 2 + ( φ y ) 2 , Δ = 2 φ 2 + 2 ψ 2 .
Introducing the following variables and parameters:
φ = φ Q , ψ = ψ Q , t = ( υ 0 2 Q ) 1 α t , C = C C 1 , β = β ( υ 0 2 Q ) γ α 1 , φ 0 = φ 0 Q , d = d Q , υ = υ υ 0
( υ 0 is the characteristic parameter of velocity), we rewrite the problem in Equations (10)–(12) in the form (the prime sign is further omitted to simplify the formulas)
( C D t α + β C D t γ ) C ( φ , ψ , t ) = υ 2 ( φ , ψ ) ( d ( 2 C φ 2 + 2 C ψ 2 ) C φ ) , ( 0 < α , γ < 1 )
C | φ = 0 = 1 , C ψ | ψ = 0 , ψ = 1 = 0 , C φ | φ = φ 0 = 0 ,
C ( φ , ψ , 0 ) = C 0 ( φ , ψ ) .
Therefore, the modeling of migration process dynamics using this fractional differential model reduces to solving the initial-boundary value problem in Equations (15)–(17) in the domain G ω × ( 0 , + ) and the following transition from the domain G ω = { ( φ , ψ ) : 0 < φ < φ 0 , 0 < ψ < 1 } to the physical domain G z according to the solution of the corresponding filtration problem presented in [5]. Within the framework of the simplified mathematical model based on Equation (6), the modeling reduces to solving the initial-boundary value problem for the convective diffusion equation
C t = υ ˜ 2 ( φ , ψ ) ( d ( 2 C φ 2 + 2 C ψ 2 ) C φ ) ν ( C C 0 )
where
υ ˜ 2 = υ 2 α + β γ , ν = 1 α + β ( 1 γ ) α + β γ , C 0 = C ( φ , ψ , 0 ) ,
and the conditions in Equations (16)–(17) are met.

3. Closed-Form Solution of the Problem in the Case of Averaged Filtration Velocity

In the case when filtration velocity is averaged over the domain of complex flow potential as υ ˜ 2 ( φ , ψ ) = υ M 2 = c o n s t , we have instead of Equation (18) the following equation with constant coefficients
C t = υ M 2 ( d Δ C ( φ , ψ , t ) C φ ) ν ( C C 0 ) .
Let us obtain a closed-form solution of the problem in Equations (19), (16), and (17). Eliminating convective term from Equation (19) by the substitution
C ( φ , ψ , t ) = 1 exp ( φ 2 d ) w ( φ , ψ , t )
we obtain for the determination of w the following initial-boundary value problem:
w t ( φ , ψ , t ) = υ M 2 d Δ w ( φ , ψ , t ) ( ν + υ M 2 4 d ) w + g ( φ , ψ ) ,
w | φ = 0 = 1 , w ψ | ψ = 0 , ψ = 1 = 0 , ( w φ + 1 2 d w ) | φ = φ 0 = 0 ,
w ( φ , ψ , 0 ) = f ( φ , ψ )
where f ( φ , ψ ) = ( 1 C 0 ( φ , ψ ) ) exp ( φ 2 d ) , g ( φ , ψ ) = ν f ( φ , ψ ) .
Application of the Fourier finite integral sine transform [25] with respect to the variable φ in the form
w ¯ m ( ψ , t ) = 0 φ 0 w ( φ , ψ , t ) sin ( η m φ 0 φ ) d φ
where η m > 0 are the roots of the equation η c t g η = φ 0 2 d ( m = 1 , 2 , . . . ) to the problem in Equations (20)–(22) yields
d w ¯ m ( ψ , t ) d t = υ M 2 d 2 w ¯ m ( ψ , t ) ψ 2 κ m w ¯ m ( ψ , t ) + Φ m ( ψ ) ,
w ¯ m ( ψ , t ) ψ | ψ = 0 , ψ = 1 = 0 ,
w ¯ m ( ψ , 0 ) = F m ( ψ ) ( m = 1 , 2 , . . . )
where
F m ( ψ ) = 0 φ 0 f ( φ , ψ ) sin ( μ m φ 0 φ ) d φ , Φ m ( ψ ) = ν F m ( ψ ) ( m = 1 , 2 , . . . )
κ m = ν + υ M 2 ( 1 4 d + d ( η m φ 0 ) 2 ) ( m = 1 , 2 , . . . ) .
The following application of the Fourier finite cosine transform in the form
w ¯ ˜ m n ( t ) = 0 1 w ¯ m ( ψ , t ) cos ( n π ψ ) d ψ ( m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . )
to the obtained problem in Equations (24)–(26) yields the Cauchy problem
d w ¯ ˜ m n ( t ) d t + ρ m n w ¯ ˜ m n ( t ) = ω m n ,
w ¯ ˜ m n ( 0 ) = δ m n ( m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . )
where
ρ m n = κ m + d ( n π ν ) 2 , δ m n = 0 φ 0 0 1 f ( φ , ψ ) sin ( μ m φ 0 φ ) cos ( n π ψ ) d φ d ψ ,
ω m n = ν δ m n ( m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . ) .
The solution of the problem in Equations (30) and (31) has the form
w ¯ ˜ m n ( t ) = ( δ m n ω m n ρ m n ) e ρ m n t + ω m n ρ m n ( m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . ) .
Then, taking into consideration the formulas of the inversion [25] of Fourier transforms with respect to geometric variables, we obtain the following closed-form solution of the initial-boundary value problem:
C ( φ , ψ , t ) = 1 4 φ 0 exp ( φ 2 d ) m = 1 η m 2 + φ 0 2 4 d 2 η m 2 + φ 0 2 4 d 2 + φ 0 2 d [ 1 2 w ¯ ˜ m 0 ( t ) + + n = 1 w ¯ ˜ m n ( t ) cos ( n π ψ ) ] sin ( η m φ 0 φ )
where w ¯ ˜ m n ( t ) ( m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . ) is defined by Equation (34).
It is worth noting that if, e.g., the conditions
f C φ , ψ 3 ( G ¯ ω ) , f | φ = 0 , φ = φ 0 = f φ | φ = 0 , φ = φ 0 = f φ φ | φ = 0 , φ = φ 0 = 0 , f ψ | ψ = 0 , ψ = 1 = 0 ,
are met, from
| w ¯ ˜ m n ( t ) | C 1 ( | δ m n | + | ω m n ρ m n | ) ( C 1 = c o n s t > 0 , m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . ) ,
we get
| w ¯ ˜ m n ( t ) | C 2 θ n η m 3 ( C 2 = c o n s t > 0 , m = 1 , 2 , . . . , n = 0 , 1 , 2 , . . . , θ n = [ 1 , n = 0 ; n 3 , n = 1 , 2 , . . . ] ) .
Hence, we conclude that the series in Equation (35) converge absolutely and uniformly in the domain Ω ε : = G ¯ ω × ( ε , T ] ( ε > 0 ) .

4. Algorithm for Numerical Solution of the Initial-Boundary Value Problem in the Case of a Variable Filtration Velocity

In the case when υ ˜ 2 ( φ , ψ ) c o n s t , the approximated solution of the initial-boundary value problem
C t = υ ˜ 2 ( φ , ψ ) ( d ( 2 C φ 2 + 2 C ψ 2 ) C φ ) ν ( C C 0 ) , ( φ , ψ , t ) G ω × ( 0 , + )
C | φ = 0 = 1 , C ψ | ψ = 0 , ψ = 1 = 0 , C φ | φ = φ 0 = 0 ,
C ( φ , ψ , 0 ) = C 0 ( φ , ψ )
can be found on the base of finite-difference approach as follows.
Introducing the grid domain
ω h = { ( φ i , ψ k , t j ) : φ i = i h 1 ( i = 0 , m + 1 ¯ ) , ψ k = ( k 0 , 5 ) h 2 ( k = 0 , n + 1 ¯ ) , t j = j τ ( j = 0 , N + 1 ¯ ) }
where h 1 , h 2 , τ are the grid steps with respect to the geometric variables φ , ψ and the time t, correspondingly, we associate with Equation (39) the following analog of the locally one-dimensional [26] finite-difference scheme of A.A. Samarskii:
1 2 C ¯ t = υ ˜ 2 ( d C ¯ φ ¯ φ C ¯ φ 0 ) ν 2 ( C C 0 ) ,
1 2 C ^ t = υ ˜ 2 d C ^ ψ ¯ ψ ν 2 ( C C 0 )
where C ^ = C j + 1 , C ¯ = C j + 1 2 , C = C j , t j + 1 2 = t j + τ 2 , C ¯ t = C ¯ C τ / 2 , C ^ t = C ^ C ¯ τ / 2 , C φ 0 , C φ ¯ φ , C ψ ¯ ψ are central and second finite-difference derivatives with respect to φ and ψ , correspondingly.
Unwinding in Equation (43) the finite-difference operators taking introduced notations into consideration and collecting similar terms, we obtain on the half-integer time step t j + 1 2 the following equations’ system:
A i k C i + 1 , k j + 1 2 B i k C i k j + 1 2 + S i k C i 1 , k j + 1 2 = Φ i k j ( i = 1 , m ¯ ; k = 1 , n ¯ ; j = 0 , N ¯ )
where
A i k = υ ˜ i k 2 h 1 ( d h 1 0 , 5 ) , S i k = υ ˜ i k 2 h 1 ( d h 1 + 0 , 5 ) , B i k = 1 τ + A i k + S i k ,
Φ i k j = ( ν 2 1 τ ) C i k j ν 2 C 0 i k .
On the integer time step t j + 1 from Equation (44) we get
P i k C i , k + 1 j + 1 Q i k C i k j + 1 + P i k C i , k 1 j + 1 = Ω i k j ( i = 1 , m ¯ ; k = 1 , n ; ¯ j = 0 , N ¯ )
where
P i k = υ ˜ i k 2 d h 2 2 , Q i k = 1 τ + 2 P i k ,
Ω i k j = ν 2 C i k j 1 τ C i k j + 1 1 / 2 v 2 C 0 i k
Let us note that the difference equations of the systems in Equations (45) and (48) are three-diagonal and, using the difference analogs of the boundary conditions in Equations (40) and (41), are effectively solved by the Thomas algorithm [26,27]. The stability of the Thomas algorithm for Equations (45) and (48) follows from the fact of diagonal prevalence in the coefficient matrices of these linear algebraic equations systems.
The brief description of the main stages of the simulation algorithm for solving the considered initial-boundary value problem based on the finite-difference technique using the locally one-dimensional scheme is as follows.
At the first stage, based on the analytical dependencies from Lavryk et al. [5], the elements of the array of velocity field values that correspond to the considered filtration scheme are calculated. Having a filtration velocity field, at the second stage, the systems of algebraic Equations (45) and (48) are solved and, as a result, the arrays of concentration field values in the domain of complex flow potential are obtained. The final stage of the algorithm consists in the transition from the domain of complex flow potential G ω to the physical domain G z based on the characteristic flow function z = f ( ω ) .

5. The Results Of Simulations

Numerical modeling of the dynamics of the considered migration process using the original MIM model and the proposed simplified mathematical model, as well as a comparative analysis of simulation results according to both of these models, were performed with respect to dimensionless variables determined by the relations in Equation (14) under the condition C 0 ( φ , ψ ) 0 .
Some obtained results for κ = 0 . 57 m/day, H = 0 . 45 m, Q = 6 m2/day, d = 0 . 1 m2/day, C 1 = φ 0 = 1 , τ = 0 . 0005 , and n = m = 30 are depicted in Figure 1, Figure 2 and Figure 3.
Figure 1 shows the curves of concentration field distribution along the streamline ψ = 0 . 5 at a fixed time t = 4 . 2 h that, according to Equation (14), corresponds to the values t = 0 . 03 for α = 1 , t = 0 . 212 for α = 0 . 84 , and t = 0 . 0256 for α = 0 . 92 . For the simulations, the results of which are shown in Figure 1, we used such values of fractional derivative exponents that yield accurate enough approximation of the initial model by the simplified one. Figure 2 shows the curves of the relative L 2 -error E R E ( α , γ ) = C C S C between the solutions according to the initial MIM model and the simplified model (here, C is the concentration that corresponds to the MIM model, C S is the concentration that corresponds to the simplified model, and · is the norm in L 2 ) depending on the value of the parameter γ for various fixed values of the parameter α . Similar curves for the E R E ( α , γ ) subject to the values of the parameter α for the fixed values of γ are presented in Figure 3.
Analysis of the numerical experiments’ results allows us to draw the following conclusions about the characteristics of the dynamics of soluble substances’ concentration field under the conditions of vertical steady-state groundwater filtration with a free surface:
  • The results of simulation using the simplified mathematical model are in good qualitative agreement with the corresponding results obtained using the original fractional-differential MIM model (Figure 1).
  • Simulation accuracy when the simplified model is used in comparison with the original MIM model is in many cases (e.g., when performing evaluative calculations in engineering practice) satisfactory. This can be illustrated by the following example: the inequality E R E ( α , γ ) < 0 . 3 in the conducted numerical experiments was met when ( α , γ ) [ 0 . 85 , 1 ] × [ 0 . 88 , 1 ] , and the inequality E R E ( α , γ ) < 0 . 2 was met when ( α , γ ) [ 0 . 94 , 1 ] × [ 0 . 97 , 1 ] (Figure 2 and Figure 3). Moreover, the relative L 2 -error between the solutions according to the original MIM model and the simplified model quickly decrease as the values of the derivatives’ orders α , γ approach to 1: α , γ 1 E R E ( α , γ ) E R E ( 1 , 1 ) = 0 . The dependency between the values of the relative L 2 -error and the obtained values of α , γ can be empirically represented (with a maximum absolute error of 0.07) as
    E R E ( α , γ ) = ( 1 γ γ 2 ) T ( 0 , 14 1 , 71 1 , 64 1 , 64 2 , 63 1 , 97 1 , 56 1 , 70 1 , 02 ) ( 1 α α 2 ) .
  • Computational complexity of the calculations using the original MIM model linearly depends on the number of time steps, while, for the simplified mathematical model, the time spent to obtain a solution on one step does not change during the simulation. Thus, the use of the simplified model linearly accelerates the modeling process.
The obtained data indicate that in some cases the simplified mathematical model of the geomigration process can be effectively used for carrying out evaluative calculations in engineering practice while developing constructive solutions in the field of designing environmentally hazardous engineering objects.

6. Conclusions

The use of a fractional-differential mathematical model of an anomalous geomigration process under the conditions of filtration taking into account the MIM approach in engineering practice presents certain severe difficulties associated with the complexity of simulation process. One of the effective approaches to simplify and accelerate the process of evaluative engineering calculations in this case is an approach based on the simplification of the original mathematical model and the transition to a simpler model based on the classical convective diffusion equation. We constructed such a simplified model and within its framework formulated the two-dimensional initial-boundary value problem of convective diffusion of soluble substances under the conditions of vertical steady-state groundwater filtration with a free surface from a reservoir into a coastal drain. In the case of averaged filtration velocity, an analytical solution of this problem is obtained. In the general case of a variable filtration velocity, an algorithm has been developed for obtaining its numerical solution. The results of simulations based on this algorithm indicate that the results of concentration fields modeling within the framework of the simplified mathematical model are in qualitative agreement with the corresponding results obtained using the original MIM model. The constructed model can be effectively used in carrying out evaluative computations in engineering practice for simplification and acceleration of simulation process.

Author Contributions

Methodology, V.B. (Volodymyr Bulavatsky); investigation, V.B. (Vsevolod Bohaienko) and V.B. (Volodymyr Bulavatsky); software, V.B. (Vsevolod Bohaienko); visualization, V.B. (Vsevolod Bohaienko); and writing—Original draft, V.B. (Volodymyr Bulavatsky) and V.B. (Vsevolod Bohaienko). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Polubarinova-Kochina, P.I. Theory of Groundwater Movement; Nauka: Moscow, Russia, 1977. (In Russian) [Google Scholar]
  2. Mystetsky, G.E. Hydroconstruction. Automation of Calculations of Mass Transfer in Soils; Budivelnyk: Kyiv, Ukraine, 1985. (In Russian) [Google Scholar]
  3. Liashko, I.I.; Demchenko, L.I.; Mystetsky, G.E. Numerical Solution of the Problems of Heat and Mass Transfer in Porous Media; Naukova Dumka: Kyiv, Ukraine, 1991. (In Russian) [Google Scholar]
  4. Bulavatsky, V.M.; Kryvonos, I.G.; Skopetsky, V.V. Non-Classical Mathematical Models of Heat and Mass Transfer Processes; Naukova Dumka: Kyiv, Ukraine, 2005. (In Ukrainian) [Google Scholar]
  5. Lavryk, V.I.; Filchakova, V.P.; Iashyn, A.A. Conformal Mappings of Physical Topological Models; Naukova Dumka: Kyiv, Ukraine, 1990. (In Russian) [Google Scholar]
  6. Berkowitz, B.; Scher, H. Theory of anomalous chemical transport in random fracture networks. Phys. Rev. E 1998, 57, 5858–5869. [Google Scholar] [CrossRef] [Green Version]
  7. Adams, E.E.; Gelhar, L.W. Field study of dispersion in a heterogeneous aquifer: 2. Spatial moments analysis. Water Resour. Res. 1992, 28, 3293–3397. [Google Scholar] [CrossRef]
  8. Boggs, J.M.; Adams, E.E. Field study of dispersion in a heterogeneous aquifer: 4. Investigation of adsorption and sampling bias. Water Resour. Res. 1993, 28, 3325–3336. [Google Scholar] [CrossRef]
  9. Chen, Z.; Zhang, R.R.; Han, S.S. An enhanced environmential multimedia modelling system (FEMMS): Part I—Development and model verification. Environ. Eng. Manag. J. (EEMJ) 2017, 16, 317–328. [Google Scholar] [CrossRef]
  10. Bohaenko, V.O.; Bulavatsky, V.M. Mathematical Modeling of Solutes Migration under the Conditions of Groundwater Filtration by the Model with the k-Caputo Fractional Derivative. Fractal Fract. 2018, 2, 28. [Google Scholar] [CrossRef] [Green Version]
  11. Bulavatsky, V.M. Mathematical modeling of dynamics of the process of filtration convection diffusion under the condition of time nonlocality. J. Autom. Inf. Sci. 2012, 44, 13–22. [Google Scholar] [CrossRef]
  12. Bulavatsky, V.M. Numerical modeling of the dynamics of a convection diffusion process locally non-equilibrium in time. Cybern. Syst. Anal. 2012, 48, 861–869. [Google Scholar] [CrossRef]
  13. Bulavatsky, V.M.; Bogaenko, V.A. Mathematical modeling of the fractional differential dynamics of the relaxation process of convective diffusion under conditions of planned filtration. Cybern. Syst. Anal. 2015, 51, 886–895. [Google Scholar] [CrossRef]
  14. Bulavatsky, V.M.; Bogaenko, V.A. Mathematical modeling of the dynamics of nonequilibrium in time convective-diffusion processes in domains with free boundaries. Cybern. Syst. Anal. 2016, 52, 427–440. [Google Scholar] [CrossRef]
  15. Hristov, J. Double Integral-Balance Method to the Fractional Subdiffusion Equation: Approximate Solutions, Optimization Problems to be Resolved and Numerical Simulations. J. Vib. Control 2015, 23, 2795–2818. [Google Scholar] [CrossRef]
  16. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity; Imperial College Press: London, UK, 2010. [Google Scholar]
  17. Sandev, T.; Tomovsky, Z. Fractional Equations and Models. Theory and Applications; Springer: Cham, Switzerland, 2019. [Google Scholar]
  18. Deans, H.A. A mathematical model for dispersion in the direction of flow in porous media. Soc. Petrol. Eng. J. 1993, 3, 49–52. [Google Scholar] [CrossRef]
  19. van Genuchten, M.T.; Wierenga, P.J. Mass transfer studies in sorbing porous media, 1: Analytical solutions. Soil Sci. Soc. Am. J. 1976, 40, 473–480. [Google Scholar] [CrossRef] [Green Version]
  20. Shumer, R.; Benson, D.A.; Meershaert, M.M.; Baeumer, B. Fractal mobile-immobile solute transport. Water Resour. Res. 2003, 39, 1296–1309. [Google Scholar] [CrossRef]
  21. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  22. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  23. Su, N. Distributed-order infiltration, absorption and water exchange in mobile and immobile zones of swelling soils. J. Hydrol. 2012, 468–469, 1–10. [Google Scholar] [CrossRef]
  24. Ahmadian, A.; Salahshour, S.; Ali-Akbari, M.; Ismail, F.; Baleanu, D. A novel approach to approximate fractional derivative with uncertain conditions. Chaos Solitons Fractals 2017, 104, 68–76. [Google Scholar] [CrossRef]
  25. Sneddon, I. The Use of Integral Transform; Mc. Graw-Hill Book Comp.: New York, NY, USA, 1973. [Google Scholar]
  26. Samarskii, A.A. The Theory of Difference Schemes; CRC Press: New York, NY, USA, 2001. [Google Scholar]
  27. Samarskii, A.A.; Vabishchevich, P.N. Computational Heat Transfer; Wiley: New York, NY, USA, 1995; Volume 2. [Google Scholar]
Figure 1. Concentrationfield along the streamline ψ = 0 . 5 in the fixed moment of time t = 4 . 2 h in the case of MIM (curves 1,4), simplified (curves 2,5) and classical (curves 3,6) mathematical models (1,2, α = γ = 0 . 92 , β = 0 . 1 ; 3, α = γ = 1 , β = 0 . 1 ; 4,5, α = 0 . 84 , γ = 0 . 96 , β = 0 . 8 ; and 6, α = γ = 1 , β = 0 . 8 ).
Figure 1. Concentrationfield along the streamline ψ = 0 . 5 in the fixed moment of time t = 4 . 2 h in the case of MIM (curves 1,4), simplified (curves 2,5) and classical (curves 3,6) mathematical models (1,2, α = γ = 0 . 92 , β = 0 . 1 ; 3, α = γ = 1 , β = 0 . 1 ; 4,5, α = 0 . 84 , γ = 0 . 96 , β = 0 . 8 ; and 6, α = γ = 1 , β = 0 . 8 ).
Fractalfract 04 00020 g001
Figure 2. Curvesof the relative L 2 -error E R E ( α , γ ) subject to the value of the parameter γ for: (1) α = 0 . 8 ; (2) α = 0 . 85 ; (3) α = 0 . 9 ; (4) α = 0 . 95 ; (5) α = 1 . 0 ; and (1–5) β = 0 . 8 , t = 0 . 01 .
Figure 2. Curvesof the relative L 2 -error E R E ( α , γ ) subject to the value of the parameter γ for: (1) α = 0 . 8 ; (2) α = 0 . 85 ; (3) α = 0 . 9 ; (4) α = 0 . 95 ; (5) α = 1 . 0 ; and (1–5) β = 0 . 8 , t = 0 . 01 .
Fractalfract 04 00020 g002
Figure 3. Curvesof the relative L 2 -error E R E ( α , γ ) subject to the value of the parameter α for: (1) γ = 0 . 8 ; (2) γ = 0 . 85 ; (3) γ = 0 . 9 ; (4) γ = 0 . 95 ; (5) γ = 1 . 0 ; and (1–5) β = 0 . 8 , t = 0 . 01 .
Figure 3. Curvesof the relative L 2 -error E R E ( α , γ ) subject to the value of the parameter α for: (1) γ = 0 . 8 ; (2) γ = 0 . 85 ; (3) γ = 0 . 9 ; (4) γ = 0 . 95 ; (5) γ = 1 . 0 ; and (1–5) β = 0 . 8 , t = 0 . 01 .
Fractalfract 04 00020 g003

Share and Cite

MDPI and ACS Style

Bohaienko, V.; Bulavatsky, V. Simplified Mathematical Model for the Description of Anomalous Migration of Soluble Substances in Vertical Filtration Flow. Fractal Fract. 2020, 4, 20. https://doi.org/10.3390/fractalfract4020020

AMA Style

Bohaienko V, Bulavatsky V. Simplified Mathematical Model for the Description of Anomalous Migration of Soluble Substances in Vertical Filtration Flow. Fractal and Fractional. 2020; 4(2):20. https://doi.org/10.3390/fractalfract4020020

Chicago/Turabian Style

Bohaienko, Vsevolod, and Volodymyr Bulavatsky. 2020. "Simplified Mathematical Model for the Description of Anomalous Migration of Soluble Substances in Vertical Filtration Flow" Fractal and Fractional 4, no. 2: 20. https://doi.org/10.3390/fractalfract4020020

APA Style

Bohaienko, V., & Bulavatsky, V. (2020). Simplified Mathematical Model for the Description of Anomalous Migration of Soluble Substances in Vertical Filtration Flow. Fractal and Fractional, 4(2), 20. https://doi.org/10.3390/fractalfract4020020

Article Metrics

Back to TopTop