20MATHEMATICAL MODELS AND OPTIMAL CONTROL ALGORITHMS FOR CHANNELS OF IRRIGATION SYSTEMS, TAKING INTO ACCOUNT THE DISCRETENESS OF WATER SUPPLY
A. J. Seytov
PHD, Chirchik State Pedagogical Institute of the Tashkent Region
B. R. Khanimkulov
Chirchik State Pedagogical Institute of the Tashkent Region
M. R. Sherbaev
Researcher Research Institute of Irrigation and Water Problems
G. U. Muzaffarova
Intern Researcher, Research Institute of Irrigation and Water Problems
A. A. Kudaybergenov
Doctoral student National University of Uzbekistan
ABSTRACT
The object of research was selected section of the southern Golodnostep main canal, located in the North-East of the Republic of Uzbekistan. In recent years, scientific research has been carried out in many countries around the world, aimed at developing mathematical models and algorithms for solving problems of optimal management of water management systems, with the use of modern information systems. Methods of mathematical modeling and algorithmization of optimal control problems for systems with distributed parameters, what is considered the main channel. Optimal changes in water consumption over time and along the length of the main channel section are obtained, opening its gates allows you to increase the amount of water flow along the length of the site. During t=34,7 min. the water flow rate at the end of the main channel section reaches a value of Q=100 m3/sec and stabilization which important for canals and water gates
Keywords: optimal control, information systems, numerical methods, channels, water distribution, discrete water supply.
INTRODUCTION
In our republic, measures are being taken to solve the problem of water delivery to consumers, which must be solved. The urgency of this problem for canals of irrigation systems is to save water resources by optimally managing the process of water delivery to consumers by providing water to a specific consumer at a set time. This problem is paid much attention to in the national economy of the republic.
Consider a section of the channel (Fig. 1, a, b) with five water intakes. The problem of water distribution taking into account the discreteness of water supply will be considered as ensuring the supply of water flow for each water intake at time point T, i.e. stepwise change in water flow, with minimal fluctuation of the water level in the channel [1, 2].
Í Qi H* ? t
1 i h H1 ¿ t
J L H" J t
J L Qi h1-" t
1 i Qs Q -i í
1 ""
Fig. 1. Discrete water supply in the canal sections
METHODOLOGY
The direct wave model. Consider a model statement of the problems of water distribution in a given section of the channel, taking into account the discreteness of water supply, taking into account only the delay in the distribution of water flow along the length of the channel [2,3].
The discreteness of water flow is formulated using a discrete unit function
qt ■ l ■ (t - T).
As a mathematical model of the channel section, the one-dimensional differential equation in the form
^^M^ = q(xA (1)
where Q(x,t) - change in water flow in the channel section, v- flow rate. Initial condition:
Q(x,0) = Q(x), (2)
where Q(x) - initial distribution of water flow in the channel section. Boundary condition:
Q(0,t) = Q(t), (3)
where Q(t) - change in water flow at the beginning of the channel section. Variable definition area
x > 0, t > 0, v > 0. (4)
Water flows at the points of water intake of the channel section q(x,t) under the conditions of discreteness of water distribution have the form
q(x,t) = -Zq. S(x - a )l(t -T). (5)
i=1
Equation (5) takes into account the main properties of the irrigation canal, such as the delay in water flow along the length of the canal. The change in water flow in the initial section of the channel section leads to its change in other sections of the channel section after a certain time, this is the delay. The farther the target is considered from the initial alignment, the greater the delay in water flow. The fundamental solution to equation (1) - the Green's function has the form [4,5,6]
G(x, 4,t) = l(x -4 )S(vt - (x -4)), (6)
where: S(x) - Dirac delta function
The analytical solution of equation (1) in the presence of its fundamental solution is determined as follows
11
Q( x,t) = J J G(x,^,t - t )w(Z,t )d&T, (7)
t 0 0
where co(x, t) - standardizing function for the boundary value problem (1) - (4), which has the form [7]
co(x, t) = q(x, t) + Q (x, t )g(x, t )&(t) - vö(x)Q1 (t). (8)
RESULTS
Consider the properties of the solution of equation (1). The change in water flow at the beginning of the channel section under zero initial conditions propagates with speed along the channel length. The wave front in this case does not change in length.
Kinematic wave model. The equation of the kinematic wave in the case of
lateral taps is written as follows [8]
d®(x,t) | dQ(x,t) = q(xt)
dt dx q ' J (9)
Q(x,t) = co(x,t )c( x,t )Jr( x,t)i .
Here Q(x,t) - change in water flow in the channel section, v - flow rate. Initial conditions
Q(x,0) = Q0(x), co(x,0) = a0(x), (10)
where: Q0(x) - initial distribution of water flow in the channel section. Boundary condition:
Q(0 , t) = Q1(t), (11)
where: Q0(x) - change in water flow at the beginning of the channel section. Variable definition area
x > 0, t > 0, v> 0. (12)
The water flow at the points of intake of the channel q(x, t) in the conditions of discreteness of the distribution of water has the form [1,9]
N
q(x, t) = -£q. S(x - a )l(t - T). (13)
Equation (9) takes into account the main properties of the irrigation canal, such as the delay and transformation of water flow along the length of the canal. In this case, the water flow rate changed in the initial section of the channel section leads to a change in water flow at other sections of the channel section after a certain time, this is called a delay, which gradually changes in time.
Q(x,t)
Fig. 2. The propagation of a square wave along the length of the channel
The convection-diffuse model is based on neglect of the inertial terms of the equations and has the form [10, 11]
dQ JdQ dK)dQ _= afxt) dt [dK dh )dx 2b\Q\ dx2 a( ' 7 ;
(14)
where: K _ flow module.
The flow module K(x,z) characterizes the magnitude of the friction forces and is determined by the following well-known formula
K = CjR,
where: R - hydraulic radius of the channel; x - wetted channel perimeter; C - Shezy coefficient.
To determine the Shezy coefficient, there is a whole series of empirical formulas. As one of them, Pavlovsky's formula can be adopted [1,12,13]
C =1Ry, y = 2,5jn _ 0,13 _ 0,75y[R [fn _ 0,l)
(15)
where n - channel roughness coefficient. Initial condition:
Q(x,0) = Qo(x), co(x,0) = a*o(x), (16)
where Q0(x) - initial distribution of water flow in the channel section. Boundary condition:
Q(0 , t) = Q1 (t),
where Q(t) - change in water flow at the beginning of the channel section. Variable definition area
x > 0, t > 0, v> 0.
(17)
(18)
Water consumption at the points of water intake of the channel section in the conditions of discreteness of water distribution
N
q(x,t) = -£q S(x - a )l(t - T).
i=1
The characteristic form of equations (1) has the form [7,14]
(19)
f + ^ ± ^ f " B(V + ^
dz , / -L — + (v ± c)—
dt 1 ydx
(20)
K
2
gm-(v + c)q.
Here (p =
i +-
1 (dm
B
dx )u _
h = const
v
,2
The initial conditions are specified as
z(x,0) = z0(x), Q( x,0) = Q(x), (21)
where Q (x), z0(x)- known features.
The boundary conditions at points xx = 0 and x2 = l are written as follows
Q(0,t)=u/t), Q(l,t)=u(t). (22)
Water flows at the points of water intake of the channel section, the right side of equation (1) under the conditions of discreteness of water distribution, has the form [2,15]
N
q(x,t) = -£ q S(x - a )l(t - T).
i=l
(23)
In these models, there is no analytical solution of equations (9), (14) and (20) under the indicated boundary conditions, since the hydraulic parameters of the water flow is a nonlinear function that depends on the shape of the cross section of the channel section.
From the expression of the lateral water intakes (23) it is seen that consumers are provided with a discrete water supply in time in the form of a step function. With stepwise functions, to solve the problem of optimal control of water distribution, it is necessary to formulate criteria for the quality of water distribution in the channels of irrigation systems under conditions of discrete water supply to consumers and a system of restrictions.
Below, we analyze methods that use the necessary optimality conditions for the optimal distribution of water in irrigation systems under optimal water supply conditions.
The gradient minimization method is used to approximate the solution of the problem [16.17] I (ujt))^min; ujtjeH, (24)
where H - Hilbert space, based on the construction of a minimizing sequence, which can be written in the form
UmJt) = Um(t)-aj'(Um(t), t0 < t < T,m = 0,1,2,...,. (25) where u0(t) - some given initial value of the control action, am - positive value, I' -functional gradient, calculated by known value um(t). If I' (um(t)) ^0, then a can choose so, so that I(um+1(t)) < I(um(t)).
From the condition of twice differentiability of the optimality criterion, we have
I(Um+1 (t)) - I(Um (t)) = am (-11' (Um (t))f + o( a m ))< 0 (26)
For all sufficiently small am> 0, if I'(um(t)) ^ 0, then process (26) stops when the given accuracy of the problem is satisfied and, if necessary, additional studies are conducted of the behavior of the functionum(t) in the vicinity of the control to determine whether it belongs toum(t) the optimal control or not. Since for optimal water distribution problems in irrigation systems under optimal water supply conditions, all optimality criteria are convex functions, therefore, the resulting control is close to optimal control.
There are various ways to select a value am. We list some of them [18.19].
1) am is selected from the condition
fm(am) = inf fm(am) = fn, , fm(am) = I(um(t)-amt(um(t))) (27)
am>0
this version of the gradient method is commonly called the steepest descent method.
An exact determination of the quantity am from (27) is not always possible, therefore, in practice, instead of (28), it uses the condition
fm*< fm(am) = fm,^5m, 5m > 0, £ ^ , (28)
t =1
or
fm,< fm(am) = (1 -K)fJ0) + f, 0 <X<Xm < 1 , (29)
Quantities Xm, dm here they characterize the error in fulfilling the condition (29).
2) consider am=a>0, then check the condition of monotony I(um+J<I(um), and if necessary, crush the value a, achieving the fulfillment of the monotonicity condition.
3) sometimes a determines from conditions
0 <s0 <am < 2/(L + 2s) , (30) where L - Lipschitz constant s, sm - positive numbers, which are the parameters of the method.
4) possible choice a from the condition
I(um ) - I(um - a J (um )) > sa\r (um )f, s < 0, (31)
to determine this a, they usually ask am =aand then crush it until the inequality written out is fulfilled.
5) a priori assignment of am from the conditions
am > 0, m = 0,l,2,...., X «m = 0, S «m <œ (32)
m ' / j m
m=0 m=0
For example, you can take am = c(k +1) b,c = const > 0,1 <b < 1. Such a choiceam
is simple to implement, but does not guarantee the conditions for the fulfillment of monotony and converges slowly.
6) in cases where the value is known in advance h = inf I(um) > -&>, can take
am = I(um(t)) - i:/\i(um(t))}[, (33)
In practice, for iterative methods, iterations continue until a condition is satisfied, i.e. criterion for ending an account. Here it is possible to use such criteria for ending the account as||um(t)-umjt)\<e, \i(um(t))-i(umjt))\ <s, u, ^\i'(um(t)^<y,
(34)
where s,5 u y - given numbers. Sometimes the number of iterations is specified.
The gradient projection method is used to approximate solution of problem (1), (2) and (3) in the case of imposing restrictions on the control functions, i.e. [20]
um(t) GÜ, (35)
where U - convex closed set from Hilbert space U.
The minimizing sequence of solving the optimization problem based on the gradient projection method is constructed according to the rule
um+l(t) = Pu («m(t)-aml'(um(t)), to < t < T,m = 0,1,2,.... (36)
where p - operator of design in the field of management. The projection um(t) onto the setU is defined as follows
IK (t) - COm (t)\\ = fum (t) - »m (t)\\ (37)
Expression (37) means that if the calculated control value belongs to the set U, then the calculated control value is accepted. If the calculated control value leaves does not belong to the set, then close values of the set U are taken. This method is convenient when there is a formula for projecting control onto many constraints U [21].
If
U = {Um(t); Umm < Umt < UmaJ , (38)
then the design operator is defined as follows
Um,n nPU Um(t) -amI(Um(t)) < Um,n, Um+1(t) =\Um(t)-amI(Um(t)) nPU Um,n < Um(t) - amI(Um(t)) < Umax, (39)
u npu u (t) -a I (u (t)) > u .
max r m\ / m \ m\ // max
If the given accuracy of the problem is satisfied and, if necessary, additional studies are carried out on the behavior of the function in the vicinity of the control um(t) to determine whether it will belong to um(t) the optimal control.
For the direct wave model, we apply the gradient projection method.
при um(x,t) -am (X)
0 < Un
min
um+1(x,t) = {
um(x,t) -am(¿0 при umin < um(x,t) -am(Л)т0 < umax, (40)
Umax nPU um(x,t) am W > umax •
vimrn при Sv0(t)-amA(0,t) m<vlmin,
to?1 (t) =W(t)-am x(0,t)m при umm < Sv0(t)-am 0,t)m < v^ (41)
Vlmax nPU SV7 (t)- amW(0, tУ > V1max■
Here X(x,t) - the conjugate variable determined by the solution of the conjugate boundary value problem for the kinematic wave model.
DISCUSSION
It should be noted that the boundary value problem for the main variables for the simplified model is solved in direct time, and the boundary value problem for conjugate variables in the opposite. These boundary value problems are solved by numerical methods, for example, finite-difference or finite elements [6,7]. With a good choice of the initial condition and initial controls, the proposed algorithm quickly converges. The initial controls are selected from the solution of the optimal water distribution problem under optimal water supply between consumers obtained using simplified models of unsteady water movement in the canal sections [1], and
the initial conditions are determined from the solution of the steady water movement in the canal sections.
The optimal distribution of water in irrigation systems under optimal water supply conditions is formulated as the problem of optimal control of systems with distributed parameters [4, 11].
The optimal control of unsteady water movement in the channels of irrigation systems was tested on the example of the section of the South Golodnostepsky main canal (SGMC), which is located in the Syrdarya and Jizzakh regions in the north-east of the Republic of Uzbekistan
The hydraulic and morphometric parameters of the SGMC section are as
"5
follows: water consumption in the channel section Q0=101 m/c; depth of water flow in this section of the channel H0 =5,05 m; gravitational constant g=9,8 m/c ; Shezy coefficient y=1/6; channel length l=17,7 km; channel bottom slope i =0,00006; the width of the water flow at the bottom of the living section of the channel B0 =39,4 m; water flow rate v0 =0,85 m /c; efficiency E=0,9.
Q (m \'C)
140 130 " 120 " no ■ 100 ■ 90 80
4000 8000 12000 15000 18000 18000 20000 1 (c)
Fig. 3. Change in water flow over time and along the length of the main canal section.
It can be seen from the figure that, after the gates are opened, the increased flow rate at the beginning of the channel section allows increasing the water consumption along the length of the specified section of the main channel. During t =
"5
20 823 c (34.7 min.), the water flow at the end of the section increases to Q = 100 m /c and stabilizes.
CONCLUSION
The results of numerical experiments show that the water flow at the end of the channel section is stabilized, which is necessary for the water intake from the channel located there.
Comparison of the results of numerical experiments and field studies conducted in this section of the main canal shows that the water level parameters in them differ slightly, their error is no more than 3-5%, this is of great importance in the national economy.
References
1. Rakhimov, S., Seytov, A., Nazarov, B., Buvabekov, B., Optimal control of unstable water movement in channels of irrigation systems under conditions of discontinuity of water delivery to consumers. IOP Conf. Series: Materials Science and Engineering 883 (2020) 012065, Dagestan, 2020, IOP Publishing DOI: 10.1088/1757-899X/883/1/012065 (№5, Scopus, IF=4,652)
2. А. Kabulov, I. Normatov, A. Seytov and A. Kudaybergenov, "Optimal Management of Water Resources in Large Main Canals with Cascade Pumping Stations," 2020 IEEE International IOT, Electronics and Mechatronics Conference (IEMTRONICS), Vancouver, BC, Canada, 2020, pp. 1-4, DOI: 10.1109/IEMTRONICS51293.2020.9216402 (№ 6, Scopus, IF= 9.936).
3. Shavkat Rakhimov, Aybek Seytov, Nasiba Rakhimova, Bahrom Xonimqulov. Mathematical models of optimal distribution of water in main channels. 2020 IEEE 14th International Conference on Application of Information and Communication Technologies (AICT), INSPEC Accession Number: 20413548, IEEE Access, Tashkent, Uzbekistan, DOI:10.1109/AICT50176.2020.9368798 (AICT) pp. 1-4,(№ 6, Scopus, IF=3,557)
4. А.V. Kabulov, A.J. Seytov, A.A. Kudaybergenov, Classification of mathematical models of unsteady water movement in the main canals of irrigation systems, International Journal of Advanced Research in Science, Engineering and Technology Vol. 7, Issue 4 , April 2020, ISSN: 2350-0328, India, pp. 13392- 13401.(№ 6, Web of science, IF=3,98)
5. Sh.Kh.Rakhimov, A.J. Seytov, A.A. Kudaybergenov, Optimal control of unsteady water movement in the main canals. International Journal of Advanced Research in Science, Engineering and Technology Vol. 7, Issue 4 , April 2020, India, ISSN: 2350-0328, pp. 13380-13391. (№ 6, Web of science, IF=3,98).
6. A.J. Seytov, A.R. Kutlimuradov, R.N. Turaev,N.K. Muradov,A.A. Kudaybergenov, Mathematical model of optimal control of the supply canal to the first pumping station of the cascade of the Karshi main canal. International Journal of Advanced Research in Science, Engineering and Technology Vol. 8, Issue 3 , March 2021. India. ISSN: 2350-0328. pp. 16790- 16797. (№5, web of science IF=6,646)
7. А. V. Kabulov, A. J. Seytov & A. A. Kudaybergenov. Mathematical models of the optimal distribution of water in the channels of irrigation systems. International Journal of Mechanical and Production Engineering Research and Development (IJMPERD) ISSN(P): 2249-6890; ISSN(E): 2249-8001 Vol. 10, Issue 3, Jun 2020, pp. 14193-14202 (№5 Scopus IF = 9.6246)
8. Sh. Kh. Rakhimov, A. J. Seytov, D. K. Jumamuratov & N. K. Rakhimova. Optimal control of water distribution in a typical element of a cascade of structures of a machine canal pump station, hydraulic structure and pump station. India. International Journal of Mechanical and Production Engineering Research and Development (IJMPERD) ISSN (P): 2249-6890; ISSN (E): 2249-8001 Vol. 10, Issue 3, Jun 2020, pp. 11103-11120. (№5 Scopus IF = 9.6246)
9. A Zh Seitov, BR Khanimkulov. Mathematical models and criteria for water distribution quality in large main irrigation canals. Academic research in educational sciences. Uzbekistan. Ares.uz. Vol. 1. №2, 2020. ISSN 2181-1385. Pp.405-415. (№5, web of science IF=5.723)
10. А. Ж. Сейтов А. Р. Кутлимурадов Р. Н. Тураев Э. М. Махкамов Б. Р. Хонимкулов. Оптимальные управления водных ресурсов крупных магистральных каналов с каскадом насосных станций ирригационных систем. academic research in educational sciences volume 2 | ISSUE 2 | 2021 ISSN: 21811385 Scientific Journal Impact Factor (SJIF) 2021: (№5, web of science IF=5.723)
11. Seytov Aybek Jumayevich, Solaeva Mehribon Norimonovna, TadjibayevIkram Uralbaevich. The product of a function and its place in physics. Solid State Technology. Vol. 63 No. 4 (2020). (№5 scopus IF=0.3)
12. Aybek Jumabayevich Seytov, Mamatqobil Nurmamatovich Esonturdiyev, Obid Sherqul Ogli Qarshiboyev, Gulhayo Baxodirovna Quzmanova. academic research in educational sciences volume 1 | ISSUE 3 | 2020 ISSN: 2181-1385 Scientific Journal Impact Factor (SJIF) 2020: 4.804. pp. 784- 788.
13. Ш. Х. Рахимов, А. Ж. Сейтов, М. Р. Шербаев, Д. Жумамурадов, Ф. Ж. Дусиёров. Структура базы данных и программные модули для моделирования управления водными ресурсами каскада насосных станций каршинского магистрального канала. Мелиорация 2019 3(89) стр. 85-91. (№5, web of science IF=0.144)
14. А.В.Кабулов, А.Ж.Сейтов, А.А.Кудайбергенов, Критерий управления задач оперативного управления водными ресурсами объектов водохозяйственных систем. ILIM ham JAMIYET. science and society Scientific-methodical journal
Series: Natural-technical sciences. Social and economic sciences. Philological scienes №2 2020. Pp.6-7.
15. Ш.Х. Рахимов, А.Ж. Сейтов, А.А. Кудайбергенов. Оптимальное управление распределением воды в магистральных каналах ирригационных систем. ILIM ham JAMIYET. SCIENCE and SOCIETY Scientific-methodical journal Series: Natural-technical sciences. Social and economic sciences. Philological scienes. pp. 810.
16. Mekhriban Salaeva, Kakhramon Eshkaraev, Aybek Seytov. Solving mathematical problems in unusual ways with excellent limits. European Scientific Conference. Пенза, 17 мая 2020 года рр. 254-257
17. Рахимов Ш. Х., Сейтов А. Ж., Кудайбергенов А. А. Критерии управления задач оперативного управления водными ресурсами объектов водохозяйственных систем. Abstracts of IX International Scientific and Practical Conference Kharkiv, Ukraine 2-4 August 2020. Стр. 125-131.
18. Сейтов А. Ж., Кудайбергенов А. А., Хонимкулов Б. Р. Моделирования двумерного неустановившегося движения воды на открытых руслах на основе проекционного метода. сборник докладов Республиканской научно-технической конференции «Инновационные идеи в разработке информационно-коммуникационных технологий и программных обеспечений» 15-16 мая 2020 года. САМАРКАНД. Стр. 60-63.
