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
and the velocity is
. For a harmonic function
, the complex potential could be considered through Cauchy-Riemann conditions. Denoting the analytic function that expresses the complex potential with
, then
,
and it’s derivative
the complex velocity is expressed, and
is the current function with
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 corresponds to uniform flow at speed U in a direction making an angle with the x-axis. The stream function is . Here we are interested in finding the velocity field . Furthermore, the complex potential corresponds to stagnation point flow with strength .
For an arbitrary point in the complex plane and , the complex potential represents a source of strength and a sink for .
A vortex of strength C at the origin is represented by the complex potential . This is again a multi-valued function, but we consider the principal form. For , rotation is anticlockwise, and for 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 or , where: is the volumetric flow rate of water per unit area (debit); is the hydraulic conductivity of the soil; is the cross-sectional area; is the hydraulic gradient in the direction of flow; 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:
where:
is the volumetric water content of the soil and
is the unsaturated hydraulic conductivity of the soil, which is a function of
;
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, , with the divergence operator. It can be written, in the case of water infiltration, as , where 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 , it simplifies to . 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:
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
, and in the other one (the one on the right) downstream, the water rises to the level
. The curves
and
(dam sole) are waterproof lines. Using the notation
q for the debit, the conditions imposed on them are
usually
.
The curves
and
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
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
the external normal vector at the boundary, the infiltration rate is
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
with
L length of the dam and the complex potential take the form
with
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
and the power lines
Choosing the complex potential under the form
using the principal form of the logarithm, we have boundary conditions (
2) as well as the conditions
The infiltrated water flow rate is
In
Figure 5, the velocity field for
parameters is presented.
2.2. The Lower Half Plane Water Infiltration Zone
The infiltrated fluid flows in the half-plane and the affixes of the points in this half-plane are .
We assume that the segment
represents the waterproof sole of the dam, where
and the half-lines
and
represent power lines, on which
respectively
. From where
because the function
, defined on the inferior half-plane is vanishing at infinity and has the following values on the boundary:
we shall consider for the complex velocity the following expression
from where we obtain the complex potential
that on the real axis becomes
In
Figure 6, for water infiltration under a dam with segment
as an impermeable base, the velocity field for parameters
is presented.
In the vicinity of the points
and
the speed suddenly changes direction and becomes infinite. To visualize this, in
Figure 7 the streamlines of the velocity on the dependence of
and
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 under the boundary , we shall use the Cauchy’s formula, , , the domain determined by , to compute the value of 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
sections
; having the peaks
, called nodes or control points located on
. With linear interpolation for
from where
Analytically computing the integrals of (
17) results:
relations in which
was renamed with
z, that is the variable of functions
and
f. For each
,
, we introduce the coefficients
. These coefficients are calculated using (
18) except of
,
,
which have the expressions:
From (
18) we obtain the algebraic system
With notations
separating the real terms from the imaginary ones in (
20)
with
the Kronocker symbol, that us one for
.
The equations of the system (
22) are not independent, and below we will choose among them accordingly a number of
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
. Accordingly, the diagonal of the matrix of unknown coefficients will be dominant, and the matrix will be well conditioned. The system becomes
Using the Cauchy theorem
and based on (
16) approximation, one obtain
that is
and after integration
from where, separating the imaginary part and taking into account the boundary conditions described in
Table 1, results
equation that closes the system (
23)–(
26) in a matriceal form
with the vector of the unknowns expressed by
where
Using (
19), that for
and
becomes
and (
21)
1 leads to
Furthermore, for
further, the complex velocity in a point
with
is determined from
and using (
16)
with
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
We consider , with , , case in which , .
In
Figure 9, the velocity field for parameters
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
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]).