Influence of a Subsidiary Weir on the Stability of a Main Structure Built on a Finite Stratum

: Some dams globally have negatively affected downstream structures. Constructing sub-sidiary weirs may solve this problem. This novel study focuses on investigating the parameters of seepage beneath the original structure and the proposed subsidiary weir. Conformal mapping and finite element methods are used for the analysis. The proposed subsidiary weir consists of a sloping central apron, flat aprons on both the downstream and upstream ends, and upstream and downstream sheet piles of varying depths. The existing structure also has sheet piles of different depths at its upstream and downstream ends, with an impervious layer situated at a specific depth below both the structures. The study derives equations for the simulation of the upwards pressure on both the structures, seepage rate, and exit gradient along the downstream bed and the filter at an intermediate location. Our own developed software for the analysis and a commercial software for numerical methods named Finite Element Heat Transfer (FEHT)-version-1are used to calculate these parameters. The accuracy of the analytical and numerical methods is verified by comparing the results with experimental data, which demonstrate a good level of agreement. This study also simulates the impacts of various factors, such as sheet pile configurations, the depth of the stratum beneath the structure, the ratio of effective heads, and the length of the intermediate filter.


Introduction
The design, construction, maintenance, and rehabilitation of hydraulic structures present significant challenges due to the impact of seepage on their stability [1,2].Ongoing global research in this specialized field is evident in the studies conducted by [3][4][5][6][7][8][9][10][11][12][13].The construction of the Grand Ethiopian Renaissance Dam (GERD) on the Blue Nile River has caused changes in the downstream regime and sediment load [14], leading to degradation along the Nile River.With multiple barrages and dams located downstream of the GERD, the degradation of riverbeds during and after its construction has become a threat to the stability of existing structures.Egypt must implement various solutions, including the construction of new barrages downstream of the existing ones, to address this issue.However, this approach is costly, is not sustainable, and leaves many other barrages and dams at risk.Thus, an appropriate solution is urgently needed.The use of subsidiary weirs downstream of the existing barrage at a certain distance was proposed by [15,16] to ensure stability.This approach results in two hydraulic structures with intermediate filters.Previous studies using the conformal mapping technique, such as those conducted by [15][16][17][18][19][20], have investigated the seepage characteristics of such structures.However, these studies have limitations in terms of their generality.One study [6] analyzed the introduction of one or two intermediate filters in hydraulic structures and their impact on seepage flow lines, demonstrating reductions of up to 72% in uplift pressure on the downstream floor with the introduction of an intermediate filter.Another study [21] discussed the use of fully penetrating cutoff walls as vertical seepage barriers to prevent groundwater seepage through an aquifer.Another study [22] addressed the issue of steady-state seepage flow beneath diversion dams situated on permeable soil, involving a cutoff wall at the downstream end of a specific width.Another study [19] investigated an inclined middle apron without a cutoff.However, these studies are applicable only to a limited range of real-life problems.Other studies [16,[23][24][25][26] examined a more general real-life problem, although their consideration of an infinite depth of the stratum is physically unrealistic.
The current research aims to investigate a more general, reliable, and practical problem that exists downstream of the GERD in both Egypt and Sudan.The downstream structure, referred to as the proposed weir, is assumed to have a sloping middle apron and horizontal aprons on the downstream and upstream sides, along with upstream and downstream sheet piles.The existing structure, referred to as the upstream structure, also features sheet piles on its upstream and downstream ends.The problem includes an impervious layer at a certain depth below both structures, along with the floor thickness.
The objective is to explore both analytical and numerical approaches to solve the problem of seepage beneath two structures, where the analytical solution incorporates the use of conformal mapping transformation to address the presence of an intermediate filter.
Equations to calculate the distribution of the upwards pressure, gradient at the exit point, and seepage flow are derived, and a computer code is developed to determine the seepage parameters.According to the best knowledge of the authors, such issues related to the two structures founded on a finite stratum have not been addressed in the past [10,11,22,26].Therefore, these tasks constitute original and novel research work that will provide valuable results for the engineering community.The obtained results using both analytical and numerical methods are compared and verified using experimental data obtained from an electrical analogue modeling scheme.

Numerical Model
In Figure 1a, the arrangement of two structural elements, along with a filter at an intermediate location and the flow domain in the Z-plane, is depicted.The first structure, referred to as AUCDEF, serves as the upstream barrage and consists of two end sheet piles and a horizontal floor.Within this structure, L 1 represents the distance from point C to point D, S 10 indicates the depth at the front face of the sheet pile on the upstream end from point A to point U, and S 11 represents the depth at the back face of the sheet pile on the upstream end from point U to C. On the other hand, S 20 represents the front face depth of the sheet pile downstream from point D to point E, S 21 denotes the back face depth of the downstream sheet pile from point E to F, and L f signifies the length of the intermediate filter extending from point F to point G.
Figure 1a illustrates the configuration of the downstream structure, known as GJLM-NOQB, which consists of a sloping apron in the middle, two sheet piles, and two horizontal aprons.Within this structure, L 2 represents the length of the upstream apron from points L to M, while L 4 signifies the length of the downstream apron from points N to O.The middle apron, spanning from points M to N, is inclined at the "απ" angle with respect to the horizontal axis, and its projected length is denoted as L 3 .Additionally, S 30 indicates the depth at the front face of the sheet pile on the upstream end from points G to point J, while S 31 represents the depth at back face of the sheet pile on the upstream from point J to point Fluids 2023, 8, 319 3 of 20 L. Similarly, S 40 corresponds to the front face depth of the downstream sheet pile from point O to point Q, while S 41 signifies the downstream sheet pile depth at the back face from point Q to point B. The parameter d represents the difference in elevation between the upstream and downstream beds.Moreover, H 1 denotes the water level on the upstream side of the existing barrage, H 2 represents the water level within the filter zone, and H 3 signifies the water level on the downstream side of the subsidiary weir.Both structures are constructed upon a finite depth of homogeneous isotropic pervious soil, denoted as T, conforming to the actual field conditions.The complex potential w is mathematically expressed as: w = ϕ + iψ where ψ is the stream function; ϕ is the velocity potential function, which is equal to kH, in which H is the total head; and k is the permeability coefficient.Figure 1a illustrates the configuration of the downstream structure, known as GJLMNOQB, which consists of a sloping apron in the middle, two sheet piles, and two horizontal aprons.Within this structure, L2 represents the length of the upstream apron from points L to M, while L4 signifies the length of the downstream apron from points N  A comprehensive solution to the problem given as w = f(z) may be derived as given below: i.
The plane (z), where z = x + iy, is mapped onto the lower half of a supplementary semi-infinite plane t; ii.Plane (w) is mapped onto the lower half of the secondary semi-infinite plane t; iii.By combining the derived equations, the w-z relationship is obtained.

Boundary Conditions
According to the illustration in Figure 1a,c, the initial streamline, labeled as ψ = 0, aligns with the boundary at the subsurface of the upstream structural element A, U, C, D, E, F. The upstream segment A 1 A, the filter FG at the intermediate location, and the segment BB 1 at the downstream end correspond to lines of equal potential.Without losing the generality, the downstream water level can be chosen as a datum.Then, ϕ A1A = −kH 1 , ϕ FG = −kH 2 , and ϕ BB1 = 0.0, respectively.
Based on the sizes of both structural components, the filter length at the intermediate location, the depth of the stratum, and the ratio of (H 2 )/(H 1 ), three different situations of subsoil flow can occur: i.
The intermediate filter can be subdivided into two sections, namely from points F to P and from points P to G, where P represents a stagnation point.The first section, extending from point F to point P, serves as an outlet for a portion of the flow from seepage originating from the upstream end.The second section, spanning from points P to G, functions as an inlet face.Consequently, the stream function within this configuration fluctuates in between -q 1 and zero along the FP length, while it transitions in between −q 2 and −q 1 along the PG length (refer to Figure 1a,c).This net flow is decreasing as H 2 /H 1 is increasing to a certain limit at which the net flow is equal to zero.If H 2 /H 1 is increased beyond this point, the net flow will be negative and the entire filter at the intermediate location will work as the inlet face.ii.The streamline, denoted by ψ = −q 2 , originates from a point on the upstream bed, which will intersect the downstream structure GJLMNOQB at a certain point P. At this juncture, the streamline splits into two paths: one follows the PG trajectory, culminating at point G, while the other follows the PB trajectory, terminating at the B point.In such special situations, the entire length of the filter at the intermediate location acts as an outlet surface, while the stream function values change along its length, ranging from 0 to −q 2 .This information can be observed in Figure 2a,c.At the P point, as given in Figure 2, the head along the downstream structure has its peak value.From that point onwards, the hydraulic head gradually decreases in either direction until it reaches the value equal to H 2 at the G point and eventually decreases to 0 at the B point.iii.In this situation, a streamline originating from point A on the upstream bed, having a stream function value of zero, and the other streamline originating from the F point on the filter at the intermediate location, also with a stream function value of zero, intersect at a point P on the upstream structure.From there, they continue their movement towards the bed on the downstream end.The whole of the filter at the intermediate location is to function like an entrance face in this scenario, where the stream function (ψ) values range from 0 to −q 2 .Along the upstream floor, the hydraulic head is at its minimum at point P, as illustrated in Figure 3a,c.Moving away from point P, the value of the head rises in either direction until it attains the value of H 2 on the F point and a value equal to H1 on the A point, respectively.Additionally, the streamline corresponding to a stream function value of q 2 aligns with the contour of the subsurface of the structure downstream.
intersect at a point P on the upstream structure.From there, they continue their movement towards the bed on the downstream end.The whole of the filter at the intermediate location is to function like an entrance face in this scenario, where the stream function (ψ) values range from 0 to −q2.Along the upstream floor, the hydraulic head is at its minimum at point P, as illustrated in Figure 3a,c.Moving away from point P, the value of the head rises in either direction until it attains the value of H2 on the F point and a value equal to H1 on the A point, respectively.Additionally, the streamline corresponding to a stream function value of q2 aligns with the contour of the subsurface of the structure downstream.

Analytical Solution
As mentioned earlier, Figure 1a provides a representation of the structure's geometry, along with the values corresponding to the velocity potential (φ) and the ψ on the boundaries in plane z.The relationship between the quantities φ and ψ is illustrated in Figures 1c-3c in the w plane.To conveniently represent points in both plane z and plane w, complex variables are utilized, where z is equal to "x + iy" and w is equal to "φ + iψ".To solve the given equations, the points in the z plane and w plane are mapped onto points that are common in the secondary plane t.This mapping is demonstrated in Figures 1b-3b, employing the transformation described by Schwarz and Christoffel.By using this transformation, the relationships  =  1 () and  =  2 () are established.These

Analytical Solution
As mentioned earlier, Figure 1a provides a representation of the structure's geometry, along with the values corresponding to the velocity potential (φ) and the ψ on the boundaries in plane z.The relationship between the quantities φ and ψ is illustrated in Figures 1c-3c in the w plane.To conveniently represent points in both plane z and plane w, complex variables are utilized, where z is equal to "x + iy" and w is equal to "φ + iψ".To solve the given equations, the points in the z plane and w plane are mapped onto points that are common in the secondary plane t.This mapping is demonstrated in Figures 1b-3b, employing the transformation described by Schwarz and Christoffel.By using this transformation, the relationships z = f 1 (t).and w = f 2 (t) are established.These relationships allow for the determination of the values of w for all corresponding values of z in principle.

First Operation
The first operation, denoted as z = f 1 (t), involves mapping the polygon A 1 AUCDEFG JLMNOQBB 1 in the z plane onto the t plane using the Schwarz-Christoffel equation.This equation is represented as: Here, the points A and B 1 (Figure 1a,b) are located at 0 and 1, respectively, while We need to determine the values of these parameters.The inclination angle of the sloping floor MN is denoted by απ, and M 1 and N 1 are complex constants.At the A point, where z is equal to zero and t is equal to zero, N 1 is found to be zero.The value of M 1 is determined by integrating Equation ( 1) between the limits χ and σ, resulting in: where M 1r represents the magnitude of M 1 .
By integrating Equation ( 1) along the boundaries and dividing this by the length of the floor of the structural element L 1 upstream, the following equations are obtained: Fluids 2023, 8, 319 8 of 20 Here, I represents the integrand of Equation ( 1).At point B 1 (t = l), as t revolves around a semi-circle of a small radius (r = (t − 1)e −iθ ), resulting in a change in z denoted as (iT), we can substitute t = 1 + re −iθ and dt = re −iθ idθ into Equation ( 1).This substitution leads to the following equation: The residual for Equation ( 14) can be given as follows (by dividing it by the floor length (L 1 ), as was done previously for Equations ( 2)-( 13)): These equations (Equations ( 2)-( 13) and ( 15)) allow the determination of all the unknown variable because the number of equations matches the number of unknowns.
Since these equations are not explicitly expressed in terms of the transformation parameters, it is not possible to find a direct solution for these parameters.However, these can be solved simultaneously using an iterative method to determine the unknown variables by minimizing the sum of the absolute residuals E: Once the values of β, τ, δ, ε, ξ, γ, η, ν, χ, σ, λ, µ, and ζ are known, the variable M 1r can be evaluated using the following equation: The numerical integration process involved assuming parameters within the range of zero to 1.These parameters were then utilized to estimate the integration using Equations ( 1)- (13).Subsequently, the error sum, as described in Equation ( 16), was estimated.The iterative process was repeated until the minimum sum was achieved based on Equation (16).Due to the presence of stagnation points, the conventional numerical integration method was deemed unsuitable.As a result, the Chebyshev technique [25] was employed instead.The integration was carried out using the Chebyshev technique, while the random search method was utilized to assume different parameter values.A FORTRAN code was developed to perform the aforementioned tasks.

Second Operation: w Equal to f 2 (t)
The second operation involves the transformation of the w plane, as shown in Figure 1c, onto the t plane, depicted in Figure 1b.This conformal transformation is represented by the following equation: Here, M 2 is a complex constant and ρ is a state variable.Different points on plane w corresponding to plane t are illustrated in Figure 1b,c.
The integration of Equation ( 18) occurs along the contour at the upstream subsurface, A-U-C-D-E-F, from zero to t, and the following equation is obtained: In Equation ( 19), ϕ = −kH at point F, ϕ = −kH 2 , and t = ξ.Thus, we have: Similarly, by integrating Equation ( 18) over the G-J-L-M-N-O-Q-B part from "γ to t", the following equation is found: At point B, ϕ = −kH 3 = 0 and t = ζ.Therefore, we have: The constants M 2 and ρ. can be determined from (20) and (22), considering known values of ξ, γ, and ζ.

Potential
In both the first and second operations, the values of M 1r , β, τ, δ, ε, ξ, γ, η, ν, χ, σ, λ, µ, ζ, M 2 , and ρ. are determined.To calculate the potential at any point along the upstream structure (AUCDEF), one can substitute the corresponding value of t into Equation (19).Similarly, for any point along the downstream structure (GJLMNOQB), the potential can be determined by replacing the equivalent t value into Equation (21).Equation ( 1) can be utilized to determine the corresponding place of the point under consideration on plane z.

Exit Gradient
The exit gradient is an important parameter that describes the hydraulic gradient along the downstream bed or intermediate filter.According to [23], the exit gradient, denoted as I exit , is calculated at any point using the following equation: By employing Equations ( 1), (18), and ( 23), the exit gradient can be derived as: Fluids 2023, 8, 319 10 of 20 To determine the maximum values of the exit gradient, which may occur at points F, G, and B (as shown in Figure 1a), one can substitute t = ξ, γ, and ζ, respectively, into Equation ( 24).

Seepage Flow
Case 1: ξ < t < γ, and ζ The flow from seepage (q 1 ), infiltrating into the filter at the intermediate location through the portion FP, and the drained flow (q 2 ) from the intermediate filter through the portion PG can be determined using the following equations: The seepage flow (q 3 ), which represents the total flow infiltrating into the downstream bed through the portion BB 1 , can be obtained using the following equation: Case 2: ξ < γ < ρ The infiltration flow (q 2 ) through the filter at the intermediate location can be determined using the equation given below: The q 3 , representing the total flow infiltrating into the downstream bed through the portion BB 1 , can be determined with the help of the following equation: Case 3: ρ < ξ < γ The drained infiltration flow (q 2 ) from the filter at the intermediate location can be calculated using the following equation: The seepage flow (q 3 ) infiltrating into the downstream bed through portion BB 1 , minus the flow that seeps through the intermediate filter via portion FG, can be obtained from the following equation: The above-mentioned equations were utilized to calculate the pressures in the upward direction exerted on the downstream and upstream structures, the discharge from seepage, and the gradients at the exit point along the filter at the intermediate location and the bed at the downstream end.To carry out these calculations, numerical methods were employed and multiple integrals were evaluated using a digital computer.For the determination of the variables of Equation (1) corresponding to the actual structural dimensions, the random search method described in [24] was employed, in conjunction with the Gauss-Chebyshev formula [25].By obtaining the values of M 2 and ρ from Equations (20) and (22), respectively, the upward pressure values at critical points associated with the actual structural dimensions were derived using Equations ( 19) and (21).The obtained results were then used to illustrate the influence of variations in stratum depth (T), structural dimensions, and filter length along with the corresponding active head ration (H 2 /H 1 ) on the features of the seepage flow.

Finite Element Solution
The finite element technique is employed extensively to analyze subsurface flow in hydraulic structures.Previous studies [26,27] utilized this method to study the seepage underneath the foundation of the Diyala weir in Iraq, with a focus on identifying sheet pile defects resulting from corrosion.Another study [28] employed the finite element method to study the seepage flow under a dam featuring an inclined sheet pile, while [29] used the finite volume method to investigate the seepage under a concrete dam.In another study, ref. [6] used the finite element method to examine the seepage under hydraulic structures equipped with an intermediate filter.Additionally, ref. [30] used the finite element method to determine the optimal angle of inclination for the cutoff underneath the structural elements.Previously, researchers [5] have used slide rule programing to investigate the seepage below hydraulic structures.Other researchers [31] focused on developing a numerical method to investigate the performance of cutoff wall systems in relation to uplift pressure and the piping phenomenon.
In the present study, we utilized the finite element method to analyze a practical problem involving a barrage with a subsidiary weir, exploring all possible aspects of the problem.The finite element method is extensively described in [32], and for this study we employed the Finite Element Heat Transfer (FEHT) program.The results from FEHT are shown in Figure 4.It is worth mentioning here that the FEHT software has some limitations in drawing all possible graphs.It can be used only to estimate the potential heads, meaning it can draw the equipotential lines only.Additionally, the software can be used to estimate the streamlines but not those corresponding to the equipotential lines.Therefore, the streamlines are not shown in Figure 4. Figure 4 illustrates the equipotential lines for the analyzed sample data, with the following values: Fluids 2023, 8, x FOR PEER REVIEW 11 of 20 pile defects resulting from corrosion.Another study [28] employed the finite element method to study the seepage flow under a dam featuring an inclined sheet pile, while [29] used the finite volume method to investigate the seepage under a concrete dam.In another study, ref. [6] used the finite element method to examine the seepage under hydraulic structures equipped with an intermediate filter.Additionally, ref. [30] used the finite element method to determine the optimal angle of inclination for the cutoff underneath the structural elements.Previously, researchers [5] have used slide rule programing to investigate the seepage below hydraulic structures.Other researchers [31] focused on developing a numerical method to investigate the performance of cutoff wall systems in relation to uplift pressure and the piping phenomenon.
In the present study, we utilized the finite element method to analyze a practical problem involving a barrage with a subsidiary weir, exploring all possible aspects of the problem.The finite element method is extensively described in [32], and for this study we employed the Finite Element Heat Transfer (FEHT) program.The results from FEHT are shown in Figure 4.It is worth mentioning here that the FEHT software has some limitations in drawing all possible graphs.It can be used only to estimate the potential heads, meaning it can draw the equipotential lines only.Additionally, the software can be used to estimate the streamlines but not those corresponding to the equipotential lines.Therefore, the streamlines are not shown in Figure 4.

Experimental Verification
Various methods could be used for the verification of the results obtained from conformal mapping and finite element techniques.A frequency response analysis [33] could be one of the options, physical modeling could be the second one, and analogue experimental techniques could be the third one.The first one requires a large set of data and the second one is very expensive.Therefore, the experimental verification based on the analogue technique was conducted to validate the solution obtained through the implementation of conformal mapping and finite element techniques for the problem at

Experimental Verification
Various methods could be used for the verification of the results obtained from conformal mapping and finite element techniques.A frequency response analysis [33] could be one of the options, physical modeling could be the second one, and analogue experimental techniques could be the third one.The first one requires a large set of data and the second one is very expensive.Therefore, the experimental verification based on the analogue technique was conducted to validate the solution obtained through the implementation of conformal mapping and finite element techniques for the problem at hand.For this purpose, an electric analogue modeling system was prepared, as described by Vaidyanathan [34].A total of twelve experiments were carried out, utilizing the sample data provided in the previous section, in addition to varying values of the parameter "S 21 /L 1 " (0.10, 0.125, 0.15), along with H 2 /H 1 (0.1, 0.25, 0.5, 0.75).
During the experiments, the exerted uplift pressure beneath both the upstream and downstream structures was measured for each combination of S 21 /L 1 and H 2 /H 1 values.These measurements were then compared with the results obtained using the conformal mapping and finite element methods.

Validation of the Models
The proposed analytical model was verified for its applicability and accuracy by comparing it with both the electrical analogue method and the numerical method.The comparison of the analytical solution, finite element method, and electric analogue method was conducted for the twelve experiments mentioned earlier, and the results are presented in Figure 5a,b.The experimental setup is given in Figure 5c.by Vaidyanathan [34].A total of twelve experiments were carried out, utilizing the sample data provided in the previous section, in addition to varying values of the parameter "S21/L1" (0.10, 0.125, 0.15), along with H2/H1 (0.1, 0.25, 0.5, 0.75).
During the experiments, the exerted uplift pressure beneath both the upstream and downstream structures was measured for each combination of S21/L1 and H2/H1 values.These measurements were then compared with the results obtained using the conformal mapping and finite element methods.

Validation of the Models
The proposed analytical model was verified for its applicability and accuracy by comparing it with both the electrical analogue method and the numerical method.The comparison of the analytical solution, finite element method, and electric analogue method was conducted for the twelve experiments mentioned earlier, and the results are presented in Figure 5a,b.The experimental setup is given in Figure 5c.The agreement between the results obtained from the analytical analysis and using the finite element method was found to be excellent, with the maximum percentage deviation between the two results being less than 1%.This indicates a high level of agreement and confirms the accuracy of the analytical method.
Furthermore, the agreement between the data obtained from the analytical solution and the electric analogue experiment was found to be good, with the maximum percentage deviation between the two results being less than 5%.Although slightly higher than the deviation observed with the finite element method, this level of agreement still demonstrates the effectiveness and reliability of the proposed analytical technique.
Figure 5a,b visually represents a comparison between the analytical solution, finite element method, and electric analogue method for the twelve experiments, further supporting the excellent and good agreement observed between these approaches.
In Figure 5b, by going downstream between 0 and 23 m for H2/H1 = 0.1, the uplift pressure is increased, which is scenario number 2. The streamline, denoted by ψ = -q2, originates from a point on the upstream bed, which will intersect the downstream structure GJLMNOQB at a certain point P. At this juncture, the streamline splits into two paths: one follows the PG trajectory, culminating at point G, while the other follows the PB trajectory, terminating at the B point.In the present scenario, the filter at the intermediate location functions as the face at the outlet, while the stream function along its length varies between 0 and −q2 (refer to Figure 2a,c).The hydraulic head along the downstream structure reaches its maximum value at point P, as depicted in Figure 2, and subsequently decreases in either direction until it achieves a value equal to H2 at the G point and zero at the B point, respectively.It is noted that Figure 5 shows only the upstream and downstream structures; the impact of the intermediate filter is not presented here.

Effect of the Relative Stratum Depth (T/L1)
The same dimensions of the upstream and downstream structure as given in the verification section were used with S21/L1 = 0.125, H2/H1 = 0.25, and variable T/L1 to assess the effect of the relative stratum depth.The agreement between the results obtained from the analytical analysis and using the finite element method was found to be excellent, with the maximum percentage deviation between the two results being less than 1%.This indicates a high level of agreement and confirms the accuracy of the analytical method.
Furthermore, the agreement between the data obtained from the analytical solution and the electric analogue experiment was found to be good, with the maximum percentage deviation between the two results being less than 5%.Although slightly higher than the deviation observed with the finite element method, this level of agreement still demonstrates the effectiveness and reliability of the proposed analytical technique.
Figure 5a,b visually represents a comparison between the analytical solution, finite element method, and electric analogue method for the twelve experiments, further supporting the excellent and good agreement observed between these approaches.
In Figure 5b, by going downstream between 0 and 23 m for H 2 /H 1 = 0.1, the uplift pressure is increased, which is scenario number 2. The streamline, denoted by ψ = −q 2 , originates from a point on the upstream bed, which will intersect the downstream structure GJLMNOQB at a certain point P. At this juncture, the streamline splits into two paths: one follows the PG trajectory, culminating at point G, while the other follows the PB trajectory, terminating at the B point.In the present scenario, the filter at the intermediate location functions as the face at the outlet, while the stream function along its length varies between 0 and −q 2 (refer to Figure 2a,c).The hydraulic head along the downstream structure reaches its maximum value at point P, as depicted in Figure 2, and subsequently decreases in either direction until it achieves a value equal to H 2 at the G point and zero at the B point, respectively.It is noted that Figure 5 shows only the upstream and downstream structures; the impact of the intermediate filter is not presented here.

Effect of the Relative Stratum Depth (T/L 1 )
The same dimensions of the upstream and downstream structure as given in the verification section were used with S 21 /L 1 = 0.125, H 2 /H 1 = 0.25, and variable T/L 1 to assess the effect of the relative stratum depth.
The potential acting under the upstream structure was slightly decreased as the relative stratum depth (T/L 1 ) was increased (Figure 6a).The exit gradient at point F was also slightly decreased as the relative stratum depth (T/L 1 ) was increased.
The potential acting under the upstream structure was slightly decreased as the relative stratum depth (T/L1) was increased (Figure 6a).The exit gradient at point F was also slightly decreased as the relative stratum depth (T/L1) was increased.Figure 6b shows that the increase in relative stratum depth (T/L1) results in a larger increase in the upward pressure exerted on the downstream structure compared to that on the upstream structure.There was a 16% increase in potential at some points on the downstream structure due to an increase in stratum depth (T/L1) from 1.0 to 2.0.The exit gradient at the critical point B was appreciably increased as the relative stratum depth (T/L1) was increased.Figure 6b shows that the increase in relative stratum depth (T/L 1 ) results in a larger increase in the upward pressure exerted on the downstream structure compared to that on the upstream structure.There was a 16% increase in potential at some points on the downstream structure due to an increase in stratum depth (T/L 1 ) from 1.0 to 2.0.The exit gradient at the critical point B was appreciably increased as the relative stratum depth (T/L 1 ) was increased.
An increase in relative stratum depth (T/L 1 ) beyond the value of 3 has no effect on the uplift pressure on all points of the upstream and downstream structures.Additionally, the impact of increasing the relative stratum depth (T/L 1 ) beyond the value of 3 will diminish on the exit gradient at both the F and B points.Therefore, the structures can be dealt with as if founded on an infinite stratum for values of (T/L 1 ) equal to or greater than 3.

Effect of the Ratio of the Effective Heads H 2 /H 1
For a constant upstream head H1, increasing the ratio of the effective heads H 2 /H 1 resulted in a significant increase in potential values along the subsurface contour of the upstream structure (Figure 5a).The gradient at exit point F decreased as the relative head H 2 /H 1 increased.The gradient at exit point F continued to decrease until it reached zero, and for further increases in H 2 /H 1 , point F began to function as an inlet when H 2 /H 1 exceeded 0.5.
Figure 5b demonstrates that an increase in the ratio of effective heads led to a noticeable increase in the potential acting on the structure at the downstream end.The gradient at the exit point denoted by B increased as the head ratio H 2 /H 1 increased.

The Impact of the Relative Length of the Filter (L f /L 1 )
The relative filter length (L f /L 1 ) had a slight impact on the potential along the contour at the subsurface of the upstream structure when increased (Figure 7a).However, it significantly decreased the potential along the floor at the downstream end (Figure 7b).
uplift pressure on all points of the upstream and downstream structures.Additionally, the impact of increasing the relative stratum depth (T/L1) beyond the value of 3 will diminish on the exit gradient at both the F and B points.Therefore, the structures can be dealt with as if founded on an infinite stratum for values of (T/L1) equal to or greater than 3.

Effect of the Ratio of the Effective Heads H2/H1
For a constant upstream head H1, increasing the ratio of the effective heads H2/H1 resulted in a significant increase in potential values along the subsurface contour of the upstream structure (Figure 5a).The gradient at exit point F decreased as the relative head H2/H1 increased.The gradient at exit point F continued to decrease until it reached zero, and for further increases in H2/H1, point F began to function as an inlet when H2/H1 exceeded 0.5.
Figure 5b demonstrates that an increase in the ratio of effective heads led to a noticeable increase in the potential acting on the structure at the downstream end.The gradient at the exit point denoted by B increased as the head ratio H2/H1 increased.

The Impact of the Relative Length of the Filter (Lf/L1)
The relative filter length (Lf/L1) had a slight impact on the potential along the contour at the subsurface of the upstream structure when increased (Figure 7a).However, it significantly decreased the potential along the floor at the downstream end (Figure 7b).The exit gradient at both points F and B slightly increased with an increase in the filter length until Lf/L1 reached 0.05.Lengthening the intermediate filter further resulted in a reduction in the gradient at exit points F and B (Figure 7c).However, increasing the filter length beyond a certain point diminished its impact on the exit gradient, and at Lf/L1 = 3.0, the two structures acted as completely separate structures.

Effect of the Downstream Sheet Pile at the End of the Downstream Structure (S41)
The downstream sheet pile (S41) at the end of the downstream structure had no effect on the potential along the subsurface contour of the upstream structure, neither did the gradient at the exit point on the downstream area of the upstream structure (Figure 8a,b).The exit gradient at both points F and B slightly increased with an increase in the filter length until L f /L 1 reached 0.05.Lengthening the intermediate filter further resulted in a reduction in the gradient at exit points F and B (Figure 7c).However, increasing the filter length beyond a certain point diminished its impact on the exit gradient, and at L f /L 1 = 3.0, the two structures acted as completely separate structures.The downstream sheet pile (S41) at the end of the downstream structure had no effect on the potential along the subsurface contour of the upstream structure, neither did the gradient at the exit point on the downstream area of the upstream structure (Figure 8a,b).The exit gradient at both points F and B slightly increased with an increase in the filter length until Lf/L1 reached 0.05.Lengthening the intermediate filter further resulted in a reduction in the gradient at exit points F and B (Figure 7c).However, increasing the filter length beyond a certain point diminished its impact on the exit gradient, and at Lf/L1 = 3.0, the two structures acted as completely separate structures.

Effect of the Downstream Sheet Pile at the End of the Downstream Structure (S41)
The downstream sheet pile (S41) at the end of the downstream structure had no effect on the potential along the subsurface contour of the upstream structure, neither did the gradient at the exit point on the downstream area of the upstream structure (Figure 8a,b).However, it had a significant impact on the potential along the contour in the subsurface strata of the downstream structure and the exit gradient downstream of the structural Increasing the sheet pile depth (S41) increased the potential and decreased the exit gradient.The downstream sheet pile depth S41 was recommended to be the minimum essential for ensuring the stability of the subsidiary weir.

Summary and Conclusions
The research we conducted focused on addressing the seepage problem beneath a barrage with a subsidiary weir that incorporates an intermediate filter, constructed on finite pervious soil.To find a solution, conformal mapping and finite element methods were employed.The upstream structure, representing an existing barrage, featured a flat floor with two end sheet piles.The proposed subsidiary weir, on the other hand, consisted of upstream and downstream flat aprons, a middle-inclined apron, and two end sheet piles.
Through the use of conformal mapping, a set of equations was derived to simulate various factors.These factors included the potential along the subsidiary weir and existing barrage, the gradients at the exit points along the filter at the intermediate location and the bed at the downstream end, and the seepage flow or the flow draining from the filter at the intermediate location.A computer program was developed to solve these equations.
We also investigated the influence of different parameters on the seepage characteristics.Specifically, impact of the relative length of the filter (Lf/L1), the head ratio (H2/H1), the relative depth of the stratum (T/L1), and the presence of downstream sheet piles (S41) were examined.
Based on the conducted research, the following conclusions were made.The behavior of the present barrage and the secondary weir is influenced by the ratio Lf/L1.When this ratio is less than 3, it affects the behavior of both structures.However, if the ratio exceeds 3, the two structures can be treated as separate entities with independent behavior.
The depth of the stratum also plays a role in the seepage characteristics.If the ratio T/L1 is less than 3, the behavior of the structures is influenced by the stratum depth.However, if T/L1 is equal to or greater than 3, the structures can be considered as founded on an infinite stratum, and the impact of the stratum depth becomes negligible.
Regarding the downstream sheet pile S41, it does not affect the seepage characteristics of the upstream structure.Its influence is limited to the seepage characteristics of the downstream structure (subsidiary weir).Therefore, the design considerations for this sheet pile depend solely on ensuring the safety of the downstream structure.However, it had a significant impact on the potential along the contour in the subsurface strata of the downstream structure and the exit gradient downstream of the structural element.Increasing the sheet pile depth (S 41 ) increased the potential and decreased the exit gradient.The downstream sheet pile depth S 41 was recommended to be the minimum essential for ensuring the stability of the subsidiary weir.

Summary and Conclusions
The research we conducted focused on addressing the seepage problem beneath a barrage with a subsidiary weir that incorporates an intermediate filter, constructed on finite pervious soil.To find a solution, conformal mapping and finite element methods were employed.The upstream structure, representing an existing barrage, featured a flat floor with two end sheet piles.The proposed subsidiary weir, on the other hand, consisted of upstream and downstream flat aprons, a middle-inclined apron, and two end sheet piles.
Through the use of conformal mapping, a set of equations was derived to simulate various factors.These factors included the potential along the subsidiary weir and existing barrage, the gradients at the exit points along the filter at the intermediate location and the bed at the downstream end, and the seepage flow or the flow draining from the filter at the intermediate location.A computer program was developed to solve these equations.
We also investigated the influence of different parameters on the seepage characteristics.Specifically, the impact of the relative length of the filter (L f /L 1 ), the head ratio (H 2 /H 1 ), the relative depth of the stratum (T/L 1 ), and the presence of downstream sheet piles (S 41 ) were examined.
Based on the conducted research, the following conclusions were made.The behavior of the present barrage and the secondary weir is influenced by the ratio L f /L 1 .When this ratio is less than 3, it affects the behavior of both structures.However, if the ratio exceeds 3, the two structures can be treated as separate entities with independent behavior.
The depth of the stratum also plays a role in the seepage characteristics.If the ratio T/L 1 is less than 3, the behavior of the structures is influenced by the stratum depth.However, if T/L 1 is equal to or greater than 3, the structures can be considered as founded on an infinite stratum, and the impact of the stratum depth becomes negligible.
Regarding the downstream sheet pile S 41 , it does not affect the seepage characteristics of the upstream structure.Its influence is limited to the seepage characteristics of the downstream structure (subsidiary weir).Therefore, the design considerations for this sheet pile depend solely on ensuring the safety of the downstream structure.

Figure 4
illustrates the equipotential lines for the analyzed sample data, with the following values:

Figure 4 .
Figure 4.A sample of equipotential lines underneath the barrage with a subsidiary weir.

Figure 4 .
Figure 4.A sample of equipotential lines underneath the barrage with a subsidiary weir.

Fluids
Experimental set-up for electric analogue technique.

Figure 5 .
Figure 5.Comparison between experimental (ex) and theoretical (Ma) results for the potential values along the subsurface contour of the upstream and downstream structures for S21/L1 = 0.125: (a) upstream structure; (b) downstream structure; (c) experimental setup for the electric analogue technique.

Figure 5 .
Figure 5.Comparison between experimental (ex) and theoretical (Ma) results for the potential values along the subsurface contour of the upstream and downstream structures for S 21 /L 1 = 0.125: (a) upstream structure; (b) downstream structure; (c) experimental setup for the electric analogue technique.

Figure 6 .
Figure 6.Effect of the stratum depth T on the potential acting along the subsurface contour of the upstream structure and downstream structure for S21/L1 = 0.125, H1 = 4.0 m, and H2= 1.0 m.

Figure 6 .
Figure 6.Effect of the stratum depth T on the potential acting along the subsurface contour of the upstream structure and downstream structure for S 21 /L 1 = 0.125, H 1 = 4.0 m, and H 2 = 1.0 m.

Figure 7 .
Figure 7. Effect of the filter length L f /L 1 on the potential acting along the subsurface contour of the upstream structure and downstream structure and on the exit gradient for S 21 /L 1 = 0.125, H 1 = 4.0 m, and H 2 = 1.0 m.

3. 5 .
Effect of the Downstream Sheet Pile at the End of the Downstream Structure (S 41 )

Fluids 2023, 8 ,Figure 7 .
Figure 7. Effect of the filter length Lf/L1 on the potential acting along the subsurface contour of the upstream structure and downstream structure and on the exit gradient for S21/L1 = 0.125, H1 = 4.0 m, and H2 = 1.0 m.

Figure 8 .
Figure 8.Effect of the downstream sheet pile on the potential (H/H 2 ) along the subsurface contour of the upstream structure and downstream structure. Abbreviations