Next Article in Journal
Forecasting Water Consumption in the Yangtze River Delta Based on Deformable Cumulative Multivariable Grey Model
Next Article in Special Issue
Robust Fisher-Regularized Twin Extreme Learning Machine with Capped L1-Norm for Classification
Previous Article in Journal
New Aspects on the Solvability of a Multidimensional Functional Integral Equation with Multivalued Feedback Control
Previous Article in Special Issue
Join Operation for Semantic Data Enrichment of Asynchronous Time Series Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex

by
Elena Corina Cipu
Department of Applied Mathematics, Centre for Research and Training in Innovative Techniques of Applied Mathematics in Engineering “Traian Lalescu” (CiTi), Faculty of Applied Sciences, University Politehnica of Bucharest, 313 Splaiul Independentei, 060042 Bucharest, Romania
Axioms 2023, 12(7), 654; https://doi.org/10.3390/axioms12070654
Submission received: 7 May 2023 / Revised: 24 June 2023 / Accepted: 27 June 2023 / Published: 30 June 2023

Abstract

:
The boundary element method (BEM) is used by applying Cauchy’s formula to the boundary of the water movement domain under a dam. By approximating the border with a polygon through linear interpolation, the relationships between the complex velocities on each edge of the polygon are analytically deduced. For the case of the flow domain described by a semi- circular closed contours, the numerical values of the velocity are computed and compared with those obtained only analytically. Conclusions on the analytical and numerical context are drawn.
MSC:
65N38; 30-04; 76-05; 31-00; 68W50

1. Introduction

Analysis of water infiltration under impermeable dams has important applications, and it deals with the study of how water seeps into the ground through permeable channels and how this affects the groundwater recharge and the stormwater management. Water infiltration is the most important way to replenish groundwater in the dam land. Groundwater is a vital source of freshwater for many ecosystems and human activities.
Permeable channels are unlined or lined with flexible materials that allow part of the runoff to infiltrate through their boundaries while conveying the rest of the runoff. Permeable channels can reduce the volume of runoff and lessen the effort needed to control it in the downstream basin. However, water infiltration also affects the hydraulic behavior of the flow in permeable channels, such as the water depth, velocity, and discharge. Therefore, analytic methods are needed to model and predict the infiltration capacity of permeable channels under different flow conditions and channel characteristics.
Cases of planar movement of groundwater are very frequently encountered in practice. They can be solved by exact methods, for simple cases, or through approximation methods, using simplified assumptions and numerical approximations (FEM, BEM). The movement of groundwater is a potential movement, and, if it is also planar, its study can be completed with the help of analytic functions of a complex variable (see [1,2]).
In a general case, for a potential flow, with the velocity potential φ the equation of equipotential lines is φ = c o n s t and the velocity is u = φ x , v = φ y . For a harmonic function φ , the complex potential could be considered through Cauchy-Riemann conditions. Denoting the analytic function that expresses the complex potential with f ( z ) = φ + i ψ , then φ x = ψ y , φ y = ψ x and it’s derivative
w ( z ) = d f d z = φ x + i ψ x = ψ y i φ y = u i v ,
the complex velocity is expressed, and ψ is the current function with ψ = c o n s t the equation of current lines.
Among the classic types of plane fluid flows that can be expressed by a complex potential we can enumerate: uniform flow; stagnation point flow; source and sink; vortex. The complex potential f ( z ) = U e i α z corresponds to uniform flow at speed U in a direction making an angle α with the x-axis. The stream function is ψ = V ( y cos α x sin α ) . Here we are interested in finding the velocity field v = ( u ( x , y ) , v ( x , y ) ) = V ( cos α , sin α ) . Furthermore, the complex potential f ( z ) = k 2 z 2 corresponds to stagnation point flow with strength k 0 .
For an arbitrary point ( a , b ) in the complex plane and c = a + i b , the complex potential f ( z ) = Q 2 π log ( z c ) represents a source of strength Q > 0 and a sink for Q < 0 .
A vortex of strength C at the origin is represented by the complex potential f ( z ) = i C 2 π log z . This is again a multi-valued function, but we consider the principal form. For C > 0 , rotation is anticlockwise, and for C < 0 rotation is clockwise.
Terzaghi’s theory is a classical method for analyzing seepage in earth dams. It involves the assumption of steady-state, two-dimensional flow through an isotropic, homogeneous soil medium. The theory uses Darcy’s law and provides a basic understanding of seepage patterns and flow velocities (see [3,4,5,6,7]). In some cases, simplified analytical solutions can be employed to estimate seepage quantities. We can assume a potential flow, meaning that the flow is irrotational and the streamlines are parallel to the dam surface.
The paper is organized as follows. In Section 2, analytical methods for water infiltration under an impermeable dam are discussed, and two cases are highlighted: semicircular and half plane water infiltration zone. The main results are then given in Section 3, in which Cauchy’s formula and the polygon decomposition of a curve are used for the analysis of water infiltration under an impermeable Dam. The solution of the system obtained leads to complex potential and complex velocity. At the end of the section, the numerical solution for the semicircular boundary is graphically represented. Section 4 end the paper.

2. Analytical Methods for Water Infiltration under an Impermeable Dam

The mathematical physics governing water infiltration under impermeable dams can be described by Darcy’s Law, the Richards and Laplace’s equation (see [3,4,5,6,7]). We consider a planar potential flow with x the spatial variable and elevation y.
Darcy’s Law relates the flow rate of water through a porous medium to the hydraulic gradient and the hydraulic conductivity. It can be written as: then q = k A h or q = k ( h + γ y ) , where: q is the volumetric flow rate of water per unit area (debit); K is the hydraulic conductivity of the soil; A is the cross-sectional area; h is the hydraulic gradient in the direction of flow; h is the hydraulic head (water table elevation); and γ is the specific gravity of water.
The Richards Equation extends Darcy’s Law by also considering the unsaturated flow of water through porous media. It considers the change in water content and hydraulic head with respect to time and space. The one-dimensional form of the Richards Equation for water infiltration can be expressed as:
θ t = · k h + S or θ t = y [ k ( θ ) h y ] + S ,
where: θ is the volumetric water content of the soil and k ( θ ) is the unsaturated hydraulic conductivity of the soil, which is a function of θ ; S represents source/sink terms, such as precipitation or extraction. This equation describes the movement of water into and through the soil as a function of time, hydraulic conductivity, and the hydraulic head gradient.
The continuity equation ensures the conservation of mass during water flow. It’s expression, θ / t + · ( θ q ) = 0 , with ( · ) the divergence operator. It can be written, in the case of water infiltration, as · q = · k h = S , where S is the source/sink term representing any external water sources or sinks.
Laplace’s equation is derived from combining Darcy’s law and the continuity equation. For steady-state conditions without any external sources or sinks ( S = 0 ) , it simplifies to 2 h = 0 . This equation describes the potential distribution of hydraulic head throughout the soil domain.
Boundary conditions are essential to solving the problem. For impermeable dam surfaces, the hydraulic head at the dam surface is constant and equal to the elevation of the dam crest. Furthermore, the hydraulic head is known or prescribed at the ground surface, often determined by the initial conditions or external factors. Depending on the specific problem, appropriate boundary conditions need to be assigned at the lateral boundaries of the computational domain, such as no flow or a specified hydraulic head. Boundary conditions and initial conditions must be specified to solve the Richards Equation. The boundary conditions typically involve the hydraulic head or water content at the soil surface and at the interface with the impermeable dam. Solving the Richards Equation numerically, often using numerical methods such as finite difference or finite element techniques, allows for the prediction of water infiltration patterns and the assessment of factors such as seepage rates, water table variations, and potential risks to dam stability caused by water infiltration.
A linearized problem of the two-dimensional steady potential flow could be used for analyzing the water infiltration under an impermeable dam using complex analysis and BEM (see [8,9]), considering a closed domain for the water flow. To mathematically describe the potential complex of water infiltration under impermeable dams, we can use the concept of potential theory, in which case the governing equation for the potential function φ and for the stream function ψ , is Laplace’s equation:
2 φ = 2 φ / x 2 + 2 φ / y 2 = 0 ; 2 ψ = 2 ψ / x 2 + 2 ψ / y 2 = 0 .
Suppose that an impermeable dam separates two accumulations of water, according to Figure 1. In one of the accumulations (the one on the left) upstream, the water rises to the level h 1 , and in the other one (the one on the right) downstream, the water rises to the level h 2 . The curves A B and C D (dam sole) are waterproof lines. Using the notation q for the debit, the conditions imposed on them are
ψ A B = c , c = const ψ C D = c + q ,
usually c = 0 .
The curves B C and D A are power lines. As the potential calculus (both the complex and the real) is given through its derivatives, we can impose the boundary conditions by adding an arbitrary constant to the velocity potential (same on all the parts of the border). So we will put
φ D A = k h 1 , φ B C = k h 2 .
We have mentioned that we do not know the value of the constant q in advance, which is why this method is often not used. Finding this constant is particularly important because the debit characterizes the flow of water infiltrated under the dam. Indeed, noting with n the external normal vector at the boundary, the infiltration rate is
A D v · ( n ) d s = A D φ n d s = A D ψ s d s = ψ D ψ A = q .
In a problem of infiltration, the debit describes the seepage flow, meaning the amount of water that seeps or is lost through soil permeability or through cracks in the dam in this context ([10]). The seepage rate depends on the hydraulic pressure difference and the hydraulic characteristics of the medium through which the water flows, as stated before. If the dam is not impermeable, the infiltration area is considered as A = L H with L length of the dam and the complex potential take the form f ( z ) = C e i θ z + log z with C = k γ H 2 4 i and θ the angle of inclination of the water surface. In Figure 2, the infiltration speed and the seepage flow are depicted.
In Figure 3 and Figure 4 the potential and the stream function dependence of spatial variables is plotted, in case with inclination of water surface and without.

2.1. Semicircular Water Infiltration Zone

We consider now the case in which the flow of the infiltrated fluid is in a semicircular domain.
We assume that the waterproof lines are
A B : z = R 2 e i θ , θ [ π , 0 ] , R 2 > R 1 , D C : z = R 1 e i θ , θ [ π , 0 ] ,
and the power lines
A D : y = 0 , x R 2 , R 1 , C B : y = 0 , x R 1 , R 2 .
Choosing the complex potential under the form
f ( z ) = k h 2 h 1 π i log z + k h 2 2 h 1 , C = 2 k h 2 h 1 > 0 ,
using the principal form of the logarithm, we have boundary conditions (2) as well as the conditions
ψ D C = k h 2 h 1 π log R 1 , ψ A B = k h 2 h 1 π log R 2 .
The infiltrated water flow rate is
q = ψ D C ψ A B = k h 2 h 1 π log R 1 R 2 .
The complex velocity is
u i v = i k H π z , H = h 1 h 2 .
In Figure 5, the velocity field for h 1 h 2 = 1 , R 2 = 2 , R 1 = 1 , k = 1 parameters is presented.

2.2. The Lower Half Plane Water Infiltration Zone

The infiltrated fluid flows in the half-plane y < 0 and the affixes of the points in this half-plane are z = r e i θ , 0 < r < , π < θ < 0 .
We assume that the segment y = 0 , x ( R , R ) represents the waterproof sole of the dam, where ψ = c and the half-lines y = 0 , x ( , R ) and y = 0 , x ( R , ) represent power lines, on which φ = k h 1 + c respectively φ = k h 2 + c . From where
v y 0 , x [ R , R ] = ψ x = 0
u y = 0 , x ( , R ) = u y = 0 , x ( R , ) = φ x = 0
because the function i a z 2 R 2 , a > 0 , defined on the inferior half-plane is vanishing at infinity and has the following values on the boundary:
i a z 2 R 2 = i a x 2 R 2 , y = 0 , x ( R , ) , a R 2 x 2 , y 0 , x ( R , R ) , i a x 2 R 2 , y = 0 , x ( , R ) ,
we shall consider for the complex velocity the following expression
w ( z ) = u ( x , y ) i v ( x , y ) = i a z 2 R 2
from where we obtain the complex potential
f ( z ) = i a log z + z 2 R 2 + K , K = k h 2 2 h 1
that on the real axis becomes
f ( z ) = i a log ( x + x 2 R 2 ) , y = 0 , x ( R , ) , i a log R arg x i R 2 x 2 , y 0 , x ( R , R ) , i a log x + x 2 R 2 + a π , y = 0 , x ( , R ) .
In Figure 6, for water infiltration under a dam with segment [ R , R ] as an impermeable base, the velocity field for parameters h 1 h 2 = 1 , R = 1 is presented.
In the vicinity of the points z = R and z = R the speed suddenly changes direction and becomes infinite. To visualize this, in Figure 7 the streamlines of the velocity on the dependence of r [ 0 , R 2 ] and θ [ π , 0 ] are plotted, which describe part of the half-plane. Also, for the same case, the potential function and stream function are depicted in in Figure 8.

3. BEM in Complex Analysis

Water infiltration can affect the stability of an impermeable dam. BEM (see [11,12,13,14,15]) can help evaluate the uplift pressure distribution on the dam foundation, which is caused by water infiltration. By modeling the dam and its surrounding soil as boundary elements, BEM can analyze the uplift pressures and assess their impact on the dam’s stability. Likewise, the seepage patterns and rates of water infiltration under an impermeable dam could be computed at various locations by discretizing the dam and the surrounding soil region into boundary elements and solving the Dirichlet problem expressed in Figure 1.
For each point z 0 under the boundary Γ = A B B C C D D A , we shall use the Cauchy’s formula, f ( z 0 ) = 1 2 π i Γ f ( z ) z z 0 d z , z 0 D , the domain determined by Γ , to compute the value of f ( z 0 ) through the values of the function on the boundary of the domain, Γ .
Within the boundary element method, the boundary Γ is approximated by a polygon consisting of N = 2 n + 2 m sections Γ j = z j , z j + 1 , j = 1 , 2 , , N ; having the peaks z j = x j + i y j , called nodes or control points located on Γ . With linear interpolation for f ( ζ )
f ( z ) Γ j = f z j z z j + 1 z j z j + 1 + f z j + 1 z j z z j z j + 1
from where
Γ f ( z ) z z 0 d z = j = 1 N f z j Γ j z z j + 1 z j z j + 1 d z z z 0 + f z j + 1 Γ j z j z z j z j + 1 d z z z 0 .
Analytically computing the integrals of (17) results:
f ( z ) = j = 1 N f ( z j ) g j ( z ) , 2 π i g j ( z ) = z z j 1 z j z j 1 log z z j z z j 1 + z z j + 1 z j z j + 1 log z z j + 1 z z j
relations in which z 0 D was renamed with z, that is the variable of functions ( g j ( z ) ) j { 1 , , N } and f. For each z : = z k = x k + i y k , j = 1 , 2 , , N , we introduce the coefficients G k j = g j ( z k ) . These coefficients are calculated using (18) except of G ( j 1 ) j , G j j , G ( j + 1 ) j which have the expressions:
2 π i G ( j 1 ) j = 2 π i lim z z j 1 g j ( z ) = z j 1 z j + 1 z j z j + 1 log z j 1 z j + 1 z j 1 z j , 2 π i G j j = 2 π i lim z z j g j ( z ) = log z j z j + 1 z j z j 1 2 π i G ( j + 1 ) j = 2 π i lim z z j + 1 g j ( z ) = z j + 1 z j 1 z j z j 1 log z j + 1 z j z j + 1 z j 1 .
From (18) we obtain the algebraic system
f ( z k ) = j = 1 N f z j G k j .
With notations
M k j = Re G k j , N k j = Im G k j , φ k = Re f ( z k ) , ψ k = Im f ( z k ) ,
separating the real terms from the imaginary ones in (20)
j = 1 N δ k j M k j φ j + j = 1 N N k j ψ j = 0 , k = 1 , 2 , , N j = 1 N N k j φ j + j = 1 N δ k j M k j ψ j = 0 , k = 1 , 2 , , N
with δ i j the Kronocker symbol, that us one for i = j .
The equations of the system (22) are not independent, and below we will choose among them accordingly a number of N = 2 n + 2 m independent equations. The choosing criterion for the equations is the following: we will look for the elements of the matrix with which the vector of the unknowns is multiplied to have diagonal elements of the form 1 M j j . Accordingly, the diagonal of the matrix of unknown coefficients will be dominant, and the matrix will be well conditioned. The system becomes
j = 1 n δ l j M l j φ j + j = n + 1 n + m N l j ψ j j = n + m + 1 2 n + m M l j φ j + j = 2 n + m + 1 2 n + 2 m N l j ψ j + ( j = 1 n N l j ) q = = k H j = 2 n + m + 1 2 n + 2 m M l j , l = 1 , , n
j = 1 n N l j φ j + j = n + 1 n + m δ l j M l j ψ j j = n + m + 1 2 n + m N l j φ j j = 2 n + m + 1 2 n + 2 m M l j ψ j ( j = 1 n M l j ) q = = k H j = 2 n + m + 1 2 n + 2 m N l j , l = n + 1 , , n + m
j = 1 n M l j φ j + j = n + 1 n + m N l j ψ j + j = n + m + 1 2 n + m δ l j M l j φ j + j = 2 n + m + 1 2 n + 2 m N l j ψ j + ( j = 1 n N l j ) q = = k H j = 2 n + m + 1 2 n + 2 m M l j , l = n + m + 1 , , 2 n + m
j = 1 n N k j φ j j = n + 1 n + m M k j ψ j j = n + m + 1 2 n + m N k j φ j j = 2 n + m + 1 2 n + 2 m δ k j M k j ψ j ( j = 1 n M l j ) q = = k H j = 2 n + m + 1 2 n + 2 m N l j , l = 2 n + m + 1 , , 2 n + 2 m .
Using the Cauchy theorem Γ f ( z ) d z = 0 and based on (16) approximation, one obtain
j = 1 2 n + 2 m Γ j f z j z z j + 1 z j z j + 1 + f z j + 1 z j z z j z j + 1 d z = 0
that is
j = 1 N Γ j f z j + 1 z j f z j z j + 1 z j z j + 1 d z + Γ j f z j f z j + 1 z j z j + 1 z d z = 0
and after integration
j = 1 N f z j z j z j + 1 + f z j + 1 z j z j + 1 = 0 , z N + 1 = z 1
from where, separating the imaginary part and taking into account the boundary conditions described in Table 1, results
j = 1 n y j + 1 y j 1 φ j + j = n + 1 n + m x j + 1 x j 1 ψ j + j = n + m + 1 2 n + m y j + 1 y j 1 φ j + j = 2 n + m + 1 2 n + 2 n x j + 1 x j 1 ψ j + q j = 1 n x j + 1 x j 1 = k H j = 2 n + m + 1 2 n + 2 n y j + 1 y j 1
equation that closes the system (23)–(26) in a matriceal form A X = B with the vector of the unknowns expressed by
X = [ ( φ j ) j = 1 , , n , ( ψ j ) j = n + 1 , , n + m , ( φ j ) j = n + m + 1 , , 2 n + m , ( ψ j ) j = 2 n + m + 1 , , 2 n + 2 m , q ] t ,
where
B l = k H j = 2 n + m + 1 2 n + 2 m M l j , l { 1 , , n } , B l = k H j = 2 n + m + 1 2 n + 2 m N l j , l { n + 1 , , n + m } , B l = k H j = 2 n + m + 1 2 n + 2 m M l j , l { n + m + 1 , , 2 n + m } , B l = k H j = 2 n + m + 1 2 n + 2 m N l j , l { 2 n + m + 1 , , 2 n + 2 m } .
Using (19), that for j = 1 and j = N becomes
2 π i G 11 = log z 1 z 2 z 1 z N , 2 π i G N N = log z N z 1 z N z N 1 2 π i G 1 N = z N z 2 z 1 z 2 log z N z 2 z N z 1 , 2 π i G N 1 = z 1 z N z 1 z N 1 log z 1 z N z 1 z N 1 ,
and (21)1 leads to
N 11 = 1 2 π log z 1 z 2 z 1 z N ; M 11 = 1 2 π arg z 1 z 2 z 1 z N N j j = 1 2 π log z j z j + 1 z j z j 1 ; M j j = 1 2 π arg z j z j + 1 z j i z j 1 , j { 2 , , N 1 } N N N = 1 2 π log z N z 1 z N z N 1 ; M N N = 1 2 π arg z N z 1 z N z N 1 .
Furthermore, for j { 2 , , N 1 }
G ( j 1 ) j = i 2 π z j 1 z j + 1 z j z j + 1 log z j 1 z j + 1 z j 1 z j , G ( j + 1 ) j = i 2 π z j + 1 z j 1 z j z j 1 log z j + 1 z j z j + 1 z j 1 , G k j = i 2 π z k z j 1 z j z j 1 log z k z j z k z j 1 + ( i ) 2 π z k z j + 1 z j z j + 1 log z k z j + 1 z k z j , k j 1 , j , j + 1 ,
further, the complex velocity in a point z 0 D with Γ = D is determined from
w ( z 0 ) = d f d z ( z 0 ) = 1 2 π i Γ f ( z ) ( z z 0 ) 2 d z
and using (16)
w ( z 0 ) = j = 1 2 n + 2 m h j ( z 0 ) f z j
with
h j ( z 0 ) = 1 2 π i z j z j + 1 f z j z z j + 1 z j z j + 1 + f z j + 1 z j z z j z j + 1 d z ( z z ) 2 = = 1 2 π i z j z j + 1 z z j + 1 z j z j + 1 d z ( z z 0 ) 2 + 1 2 π i z j 1 z j z j 1 z z j 1 z j d z ( z z 0 ) 2 = = 1 2 π i 1 z j z j 1 log z j z 0 z j 1 z 0 1 z j + 1 z j log z j + 1 z 0 z j z 0 .

Semicircular Water Infiltration Zone

For the semicircular boundary, we shall compare the results with those obtained through analytical methods.
As expressed in Table 2, for the nodes in Cartesian coordinates, we have
z j A B , x j = R 2 cos ( n + 1 j ) π n , y j = R 2 sin ( n + 1 j ) π n z j B C , x j = R 2 R 2 R 1 ( j n 1 ) m , y j = 0 z j C D , x j = R 1 cos ( j n m 1 ) π n , y j = R 1 sin ( j n m 1 ) π n z j D A , x j = R 1 R 2 R 1 ( j 2 n m 1 ) m , y j = 0 .
We consider z 0 D , with x 0 = r cos ( θ 0 ) , y 0 = r sin ( θ 0 ) , case in which R 1 < r = | z 0 | < R 2 , θ 0 [ π , 0 ] .
In Figure 9, the velocity field for parameters h 1 h 2 = 1 , R 2 = 2 , R 1 = 1 , k = 1 in case of a semicircular water infiltration zone, is presented.

4. Conclusions

The infiltration process of check dams is complex and is affected by many factors such as rainfall, soil characteristics, land cover, slope, and evapotranspiration. Analytical methods based on potential theory and complex analysis can help to simulate and predict the process of water infiltration from dams, allowing for more precise and effective control than empirical models (see [16]). This can provide a theoretical basis for the optimal use of soil water and the design of permeable channels for stormwater management (see [17]).
Different types of complex potentials expressed in literature were used to simulate the water infiltration, and the quantities derived from knowing the potential were determined using the original MatLab R2022b codes. A detailed determination of complex potential for particular cases, such as semicircular crowns and lower half-plane water infiltration zones were considered. The velocity field for both cases was precised.
Based on the analytical calculus made through the boundary element method in Section 3, it is possible to obtain the numerical values for the complex potential and the stream function as well as for the flow rate (debit constant q) using the (23)–(26) system in the matriceal form. The solution could be used if the boundary Γ of the flow domain is known, such as in the semicircular case. Using the two functions, the complex velocity is analytically obtained through Cauchy’s formula.The system could be used for any smooth, closed curve describing the boundary of the flow domain when the parameterization of the curve is known.
If only the correspondences of the points ( x j , y j ) are known through experimental measurements then curve fitting could be used to improve the accuracy of the results more than obtained with the linearization used in (16). The accuracy of the original solution of the Dirichlet problem proposed in Figure 1 obtained in Section 3 can be improved by replacing the linear interpolation, at least with Spline functions whose analytical form can be used. Other studies, such as ([18]), are made using artificial neural networks instead of BEM ([19]).

Funding

The work has been founded through the University Politehnica of Bucharest, “PubArt” program.

Data Availability Statement

Not applicable.

Acknowledgments

The work has been founded through the University Politehnica of Bucharest, “PubArt” program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fetter, C.W. Applied Hydrogeology, 4th ed.; Prentice-Hall, Inc.: Upper Saddle River, NJ, USA, 2001; ISBN 0-13-088239-9. [Google Scholar]
  2. Mateescu, C. Hidraulica, Hydraulics; Didactica si Pedagogica: Bucharest, Romania, 1961. [Google Scholar]
  3. Bear, J.; Cheng, A.H.D. Seepage Drainage and Flow Nets. Rev. Geophys. Space Phys. 1980, 18, 729. [Google Scholar] [CrossRef]
  4. Barry, D.A.; Chandler, J.H.; Groenevelt, J.P.E. Analytic Elements in the Solution of Laplace’s Equation: A Review. Adv. Water Resour. 2004, 27, 931–947. [Google Scholar] [CrossRef]
  5. Rubin, Y.; Langevin, C.D.; Wolfsberg, A.W. MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model—User Guide to the LMT6 Package. In Techniques and Methods, 6th ed.; U.S. Geological Survey: Reston, VA, USA, 2003. [Google Scholar]
  6. Salas, J.D.; Hite, R.A.; Bausch, R.E. Two-Dimensional Flow in Heterogeneous Aquifers: A Review of Analytical Methods. Adv. Water Resour. 1991, 14. [Google Scholar] [CrossRef]
  7. van Genuchten, M.T. A Closed-Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  8. Carabineanu, A.; Dinu, A.; Oprea, I. The Application of the Boundary Element Method to the mAgnetohydrodynamic Duct Flow. Z. Für Angew. Math. Und Phys. 1995, 46, 971–981. [Google Scholar] [CrossRef]
  9. Carabineanu, A. The study of the potential flow past a submerged hydrofoil by the complex boundary element method. Eng. Anal. Bound. Elem. 2014, 39, 23–35. [Google Scholar] [CrossRef]
  10. Saffman, P.G. Vortex Dynamics; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  11. Brebbia, C.A.; Dominguez, J. Boundary Elements: An Introductory Course; WIT Press: Billerica, MA, USA, 1992; ISBN 978-1-85312-349-8. [Google Scholar]
  12. Beskos, A.; Kolmogorov, K. Boundary Element Methods in Engineering Science; Springer Science and Business Media: Berlin/Heidelberg, Germany; Imperial College Press: London, UK, 2008; ISBN 978-1-84816-579-3. [Google Scholar]
  13. Huang, Y.S.; Ling, P.J. The Application of Boundary Element Method to the Seepage Analysis of Water Conservation Dams. J. Hydraul. Res. 1992, 30. [Google Scholar] [CrossRef]
  14. Leung, C.M.; Tse, P.W. Application of the Boundary Element Method to Flow through Anisotropic Porous Media. Eng. Anal. Bound. Elem. 1995, 15. [Google Scholar] [CrossRef]
  15. Willis, J.R.; Brebbia, C.A.; Santos, G.D. Boundary Element Techniques: Theory and Applications in Engineering; Springer Science and Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  16. Scrădeanu, D.; Gheorghe, A. Hidrogeologie Generală, General Hydrogeology; Bucharest University Press: Bucharest, Romania, 2007; ISBN 978-973-737-367-0. [Google Scholar]
  17. Song, Z.; Zhao, F.; Cui, Q.H.; Wang, J. Stability Analysis of Tailings Dam under Muddy Water Infiltration. In Proceedings of the 2015 International Conference on Architectural, Civil and Hydraulics Engineering (ICACHE 2015), Guangzhou, China, 28–29 November 2015; pp. 279–283. [Google Scholar]
  18. Zhang, H.; Song, Z.; Peng, P.; Sun, Y.; Ding, Z.; Zhang, X. Research on seepage field of concrete dam foundation based on artificial neural network. Alex. Eng. J. 2021, 60, 1–14. [Google Scholar] [CrossRef]
  19. Carabineanu, A.; Dinu, A. The study of the incompressible flow past a smooth obstacle in a channel by the boundary element method. Rev. Roum. Sci. Techn. Mec. Appl. 1993, 38, 601–616. [Google Scholar]
Figure 1. Water infiltration under an impermeable Dam.
Figure 1. Water infiltration under an impermeable Dam.
Axioms 12 00654 g001
Figure 2. Infiltration speed and seepage flow: L = 2, H = 1, k = 0.1, θ = 0 .
Figure 2. Infiltration speed and seepage flow: L = 2, H = 1, k = 0.1, θ = 0 .
Axioms 12 00654 g002
Figure 3. Potential function for different inclination angles (L = 2, H = 1).
Figure 3. Potential function for different inclination angles (L = 2, H = 1).
Axioms 12 00654 g003
Figure 4. Stream function for different inclination angles (L = 2, H = 1, k = 0.1).
Figure 4. Stream function for different inclination angles (L = 2, H = 1, k = 0.1).
Axioms 12 00654 g004
Figure 5. Velocity field for semicircular case.
Figure 5. Velocity field for semicircular case.
Axioms 12 00654 g005
Figure 6. Velocity field for half-plane case.
Figure 6. Velocity field for half-plane case.
Axioms 12 00654 g006
Figure 7. Velocity streamlines for half-plane case.
Figure 7. Velocity streamlines for half-plane case.
Axioms 12 00654 g007
Figure 8. Potential function and stream function for half-plane case.
Figure 8. Potential function and stream function for half-plane case.
Axioms 12 00654 g008
Figure 9. Velocity field for a semicircular water infiltration zone, numerical result.
Figure 9. Velocity field for a semicircular water infiltration zone, numerical result.
Axioms 12 00654 g009
Table 1. Boundary conditions specifying.
Table 1. Boundary conditions specifying.
Boundary Γ Nodes z j UnknownsConditions
A B : waterproof bed j = 1 , , n φ j ψ j = q
B C : feeding surface j = n + 1 , , n + m ψ j φ j = 0
C D : bottom dam j = n + m + 1 , , 2 n + m φ j ψ j = 0
D A : feeding surface j = 2 n + m + 1 , , 2 n + 2 m ψ j φ j = k H
Table 2. Boundary parametrization for semicircular case.
Table 2. Boundary parametrization for semicircular case.
Boundary Γ Nodes z j Nodes Values
A B j = 1 , , n z j = R 2 exp ( i ( n + 1 j ) π / n )
B C j = n + 1 , , n + m z j = R 2 R 2 R 1 ( j n 1 ) m
C D j = n + m + 1 , , 2 n + m z j = R 1 exp ( i ( j n m 1 ) π / n )
D A j = 2 n + m + 1 , , 2 n + 2 m z j = R 1 R 2 R 1 ( j 2 n m 1 ) m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cipu, E.C. Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex. Axioms 2023, 12, 654. https://doi.org/10.3390/axioms12070654

AMA Style

Cipu EC. Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex. Axioms. 2023; 12(7):654. https://doi.org/10.3390/axioms12070654

Chicago/Turabian Style

Cipu, Elena Corina. 2023. "Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex" Axioms 12, no. 7: 654. https://doi.org/10.3390/axioms12070654

APA Style

Cipu, E. C. (2023). Analysis of Water Infiltration under Impermeable Dams by Analytical and Boundary Element Methods in Complex. Axioms, 12(7), 654. https://doi.org/10.3390/axioms12070654

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop