Main model. Setting a task. The flat established current from A’A Sthe 2 l-width lpit, enclosed by symmetrical S-length sprunts, isconsidered, through a waterproof layer of T power soil with a substline aquifer E’E containing underground or artesian waters, The pressure in which has a permanent value H₀ (Figure 1) On the roof of this layer is a waterproof area, D’D simulated horizontal segment length 2L. Due to thesymmetry of the motion picture, we will limit ourselves to the study of the right half of the ABCDEGR filtrationarea.6

Groundwater, dflowing ARG spool under the influence of pressure difference in the pit and the underlying well-permeable aquifer, rise behind it to a certain height of the RG and, overcoming the point of M zero speed on the spool, form a free surface of GE ,which receives infiltration waters with the intensity of the e (07

We will assume that the movement of groundwater is subject to the law of Darcy with a known filtration factor of z and occurs in a homogeneous and isotropic soil, which is considered incompressible.8

Under such conditions, as it is known, the basic fluid filtration equations can be recorded as

….(1)

Where φ is the potential for filtration speed; u and v - projections on the axis of filtration rate coordinates; h - pressure; p - pressure in the stream; γ is the specific weight of the liquid 9

From a mathematical point of view, the challenge is to find the complex potential of the flow of z (φ the potential of speed, ψ is the function of the current) as an analytical within the filtration area of the function of the complex coordinates z under the following boundary conditions:10

AB : y = 0, φ = −H; BC : x = 0, ψ = 0;

CD : y = −T, ψ = 0; AG : x = l, ψ = Q; ….(2)

DE : y = −T, φ = −H₀; GE : φ = −y − T, ψ = Q + ε(x − l),

The study is carried out in terms of the values of z and ω which are related to the actual values of the same through the means of equality.11

.

Building a solution to the edge problem. Let's turn to the area of complex velocity w (Figure 2) corresponding to the edge conditions (2). This area, which is a circular polygon with three incisions, tops N₁ and N₂ two of which correspond to the extremes of current function on impenetrable areas of AB and DE, belongs to the class of polygons in the polar grids. However, unlike the "15" limit, an additional corner point appears on the boundary of the traffic area, point B (see Figure 1); the total number of special points becomes nine, making the task much more difficult.12

To solve the regional problem, the Semibarin-Kochina method, which is based on the application of the analytical theory of linear differential equations of the Fuchs class. Entered

Auxiliary parametric variable ζ and z functions : conformally displaying the upper half-flatness on the z area when dots are matched 13

ζ_E = 0, ζ_G = 1, ζ_C = k⁻² (0 < k < 1), ζ_D⁼ ∞,

as well as derivatives

.

By determining the characteristic performance of the functions of the and the "regular special points" near the regular special points, we will find that they are linear combinations of the two branches of the next function of Riemann.14

….(3)

It can be seen that the points ζ = ζ_A and ζ= ζ_Bare the common points of function Y, representing the last symbol of Riemann, which corresponds to the linear differential equation of the Fuchs class with seven regular special points, very typical for the tasks of underground hydromechanics.15

…..(4)

Recall that along with ζ the ζ_F, ζ_N1 and ζ_N2 in the equation (4) the accessory parameters of_0,λ1 and 2remain unknown in setting the task and must be determined in the course of its solution.

Replacing variables

ζ = sn²(2Kτ,k)…(5)

translates the upper half-flatness ζ into a rectangle of the plane:

I'm where , ρK(k) is a full elliptical integral of the first kind with module ( k, 23, when dots are matched

τ_E = 0, τ_G = 1/2, τ_C = (1 + iρ)/2, τ_D = iρ/2,

and the Y equation integrals (4), which correspond to the symbol of Riemann (3) and are constructed according to the technique developed earlier, are transformed as follows:

Y_1,2(τ) = θ₀⁻³(τ)θ₁(τ ± iγ)θ₂(τ ± iβ)θ₂(τ * iβ)exp(±iπτ). (6)

Here, sn τ = sn (u.k) is Jacobi's elliptical function (sinus) at module k, θ₁((q),θ₀((q) - theta function with the parameter q q exp, which is uniquely associated with the module k , q, γ - some suitable permanent.16

Taking into account the ratios (3), (5) and (6), as well as the fact that the functionshas the same appearance.

...(7)

√ε = thπ(ρ/2 + β − α − γ), ..(8)

come to the addictions we're in.

….(9)

in which N is a large-scale constant simulation, sn- unknown residents of points A and B of the area of q.

(9) Permanent conformal displays of q, z and γ, which are bound by a ratio (8), are subject to the conditions

0 < α < r < β < m < a < b < ρ/2, 0 < γ < ρ/2,...(10)

Regulating the position on the boundary of the flow of points M and point of the spun R, as well as N₁ and N_2;; m and r are unknown residents of M and R points in the plane. τ.

You can verify that the functions (9) meet the conditions (2) reformulated in terms of functions, and thus are a parametric solution to the original edge problem.17

Recording ratios (9) for different parts of the boundary of the area, followed by integration across the entire loop of the auxiliary area, leading to a close of the movement area and thus serves as a control of the calculations.18

….(11)

….(12)

and coordinates of the EG depression curve points

….(13)

Account control is other expressions for the values d,L and filtration expense::

….(14)

In formulas (11)-(14) subtegrial functions are expressions of the right parts of equality (9) on the corresponding sections of the contour of theregion.19

Analysis of numerical results for the main filtration model. Views (9) --(14) contain seven unknown permanents: residency prototypes of points A,B,R in the plane τq , parameters of conformal display, γ, satisfying ratio (8) and inequities (10), as well as module k (0 < k < 1) and constant N modeling. N

The first means that the speed at the end of the spool turns into infinity, and the second is directly derived from the consideration of boundary conditions (2). After determining the unknown permanents are the exact values d and z by expressions (12) and finally, the formulas (13) are calculated coordinates of the points of the free surface EG.20

On rice. 1 is a picture of the current, calculated at q= 0. 6, To= 7,S,H= 3,H 3,L= 15,H, = 7L, 10 (basic values). and the rest are fixed by the baseline values, the results of calculations of the influence of the determining physical parameters are 0,T,S,H and L on the depth of D (negative values d means that the free surface rises above the absiss axis) and the expense of the q. Figure 3 are represented by the dependence of the d and filtration consumption of the 0,T,H,H,H,H.L.H. ε,T,S,H₀,L,H,l.21

Analysis of table and rice data. 3, and allows you to draw the following conclusions.

First of all, the same quality of dependencies of the values H _,L,d,S and l, associated with a decrease Q infiltration intensity and greater reservoir power, in which these filtration characteristics grow. 22

Increase infiltration intensity, impenetrable inclusion and pressure in the underlying layer and reduce the power of the layer, The length of the spool, the d pressure of the water in the pit ε,S,H and its width lead to a drop in the depth d,i.e. to the increase of the residents of the point G exit of the depression curve from under the spool. 23. Table 1.