On the Chernous’ko time-optimal problem for the equation of heat conductivity in a rod Текст научной статьи по специальности «Математика»

URAL MATHEMATICAL JOURNAL, Vol. 5, No. 1, 2019, pp. 13-23

DOI: 10.15826/umj.2019.1.002

ON THE CHERNOUS'KO TIME-OPTIMAL PROBLEM FOR THE EQUATION OF HEAT CONDUCTIVITY IN A ROD1

Abdulla A. Azamovt, Jasurbek A. Bakhramov^, Odiljon S. Akhmedov^tt

Institute of Mathematics, National University of Uzbekistan named after Mirzo Ulugbek, Durmon yuli st., 29 Tashkent, 100125, Uzbekistan

[email protected], [email protected], [email protected]

Abstract: The time-optimal problem for the controllable equation of heat conductivity in a rod is considered. By means of the Fourier expansion, the problem reduced to a countable system of one-dimensional control systems with a combined constraint joining control parameters in one relation. In order to improve the time of a suboptimal control constructed by F.L. Chernous'ko, a method of grouping coupled terms of the Fourier expansion of a control function is applied, and a synthesis of the improved suboptimal control is obtained in an explicit form.

Keywords: Heat equation, Time-optimal problem, Pontryagin maximum principle, Suboptimal control, Synthesis of control.

Introduction

It is known that a time-optimal problem occupied a very important place in the foundation and development of optimal control theory. Even for simple non-trivial cases, the problem required working-out new approaches and lead after all to Pontryagin's maximum principle [3, 10, 30]. Despite 70 years of development, the solution of concrete non-trivial examples of time-optimal control still needs considerable effort [2, 4, 19]. The problem becomes even more difficult when a control system is described by a partial differential equation [11, 24, 25, 34], particularly, for the heat conductivity equation [12, 22, 26, 29, 35, 36]. In [1], the correctness of parabolic equations for heat propagation is discussed and for that purpose, a parabolic equation with time delay is considered.

Here, the maximum principle can be formally written out as well, but it loses its effectiveness as compared with a finite-dimensional case or on cases when the time interval is fixed [2, 9, 18, 25, 32, 33]. Therefore, Chernous'ko suggested [13] another approach based on the Fourier expansion that allowed him to reduce the problem to an infinite system of one-dimensional problems whose control parameters are connected by a condition in the min-max form (see below (1.4)) generating a closed convex control set in a Hilbert space. Unfortunately, to deal with such a constraint is quite difficult (about other kinds of constraints see [17]). In order to overcome this complexity, the mentioned constraint was replaced [13] by an infinite system of separated conditions for scalar control parameters that can be interpreted as if one took Hilbert's brick inscribed into the control set. As a result, this approach made it possible to construct a suboptimal control and to give an explicit upper estimation for an optimal time. In [5], a co-Hilbert's brick inscribed into the control set was considered, and an improved suboptimal control function was constructed. In the present

1This work was supported by a grant from the Ministry of Innovative Development of the Republic of Uzbekistan (Project No. OT-$4-84).

paper, we suggest another way for constructing a suboptimal control function in the case of the heat conductivity equation in a rod.

1. Preliminaries

As it was noted above, Chernous'ko considered the time-optimal problem for an evolutional equation

du(t, x)

dt

with the initial and boundary conditions

= •)](t,x)+ v(t,x) (1.1)

u(0,x) = u0(x), Mu(t,s)= u* (t,s), (1.2)

where A is a uniformly elliptic differential operator, t > 0, x € D, D is a regular domain with Lyapunov boundary r, s € r, and M is a boundary operator [13].

The constraint on the control function in problem (1.1), (1.2) is bounded in the norm of the space L^,; i.e. |v(t,x)| < v0 for almost all t and every x € D, where v0 is a given positive number [31]. It is known that, for every control function v(t,x), problem (1.1), (1.2) has a unique solution u(t,x) [14, 21, 28].

If a solution u(t, x) of problem (1.1), (1.2) satisfies the condition u(T, x) = 0 at some T, T > 0, then the corresponding control function v(t, x) is called admissible, and the number T is called the transition time (from the initial state u0(-) into the equilibrium state u(t,x) = 0). Let V be the class of all admissible controls. Then the quantity T = T[v(-, ■)] will be a functional on V at every fixed u0(x) and u*(t, s).

If an admissible control v*(t,x) satisfies the condition T* = T[v*( ■, ■)] < T[v( ■, ■ )] for all v( ■, ■ ) € V, then v*( ■, ■ ) is called a time-optimal control, and the value T* is called optimal transition time.

The direct application of the Pontryagin maximum principle to problem (1.1), (1.2) is a very hard task, unlike optimization problems on a finite interval of time (see [4, 8, 15]). For example, in [25], only theorems on the existence of optimal control and the bang-bang principle are given, but no specific example of a solution was considered. In monograph [11], the time-optimal problem when a control parameter participates in boundary conditions was considered [11, Ch. 5, Sect. 1] and, instead of the necessary conditions, the method of the L-momentum of N.N. Krasovskii [19] was applied [11, Sect. 2]. In the recently published article [20], Butkovsky's approach was applied to the case of a fractional-order diffusion equation. It should be noted that the L-momentum method only allows one to simplify to some degree the time-optimal problem and rarely gives an explicit solution. Therefore, the approach suggested by Chernous'ko [13], where the method of expansions on the system of eigenfunctions of the operator A was used, seems to be more effective. That helped to reduce considering problem to the infinite system of one-dimensional control problems:

Vk = -Afc Vk + vk, /k(0) = Vk0, k = 0,1,2,.... (1.3)

(About solution of systems of this kind, see [16]).

In terms of system (1.3), the condition |v(t, s)| < v0 means that a counting system of the control parameters vk, k = 0,1,2,..., should satisfy the combined constraint

<x

maxlV ^k(x)vk < v0. (1.4)

xeD 'k=0

where ^k are eigenfunctions of the problem.

Condition (1.4) defines some closed convex set L in the Hilbert space l2, which is difficult to deal with. In this connection, it is natural to try to solve the problem of finding a suboptimal control. (It is essential to note that, if a time interval is fixed, then the method of penalty functions is enough effective for the construction of a suboptimal control. It would be interesting to apply this method for the time-interval problem as well [6].) For this purpose, in [13], constraint (1.4) was replaced by a more rigid system of constraints in the form

|vfc |< Ufc, k = 0,1, 2,..., (1.5)

where = max (x)|. Wherein, nonnegative numbers should be chosen satisfying the con-

xSD

dition Y,fc=0 afcUfc = V0.

Let T*fc be an optimal transition time in the problem

y/fc = -Afcyfc + vfc, yfc(0) = (1.6)

such that (T*k) = 0, k = 0,1,2.... In [13], it is shown that the numbers can be chosen so that all T*k coincide: T*k = T for some T. Let -0k(t) be the sequence of the corresponding optimal controls. Then T* < T and V* (t, x) = ^(x)v*k(t) may serve as the sought suboptimal control.

A new problem arises here: is it possible to use a more exact reduction of the constraint than (1.5)? As mentioned above in [5] it was used Hilbert's co-cube instead of (1.5). Here, we are going to follow another approach based on a special grouping of terms of (1.5). Effectiveness of this approach is tightly related to specific properties of eigenfunctions (■), so here it will be demonstrated for the operator A = d2/dx2 connected with the process of the heat conductivity in a rod.

2. A method of grouping terms of the Fourier expansion

Consider the following concretization of problem (1.1), (1.2):

du d2u , , . , ,.

~Qt = + V№> X)> rr) I < , t > 0, 0 < X < IT, ^21)

u(0,x) = u0 (x), u(t, 0) = 0, u(t, n) = 0.

The system of eigenfunctions (t) = sin kx, k = 1,2,... , of the operator d2/dx2 forms a complete orthogonal basis of the space L2[0, n] [21, 28].

Let u(t,x) = fc=1 sin kx and v(t,x) = ^vk sin kx be the Fourier expansions on the basis {sin kx} . Then the restriction (1.4) takes the form

max I > vk sin kx < v0. (2.2)

0<x<n I ^ k " 0 v 7

~ ~ k=i

Let us consider a more rigid restriction

max > Ivfc sin kx + v3k sin3kx < v0 (2.3)

0<x<n ' I fceQ

instead of (2.2), thereby replacing the optimal control problem with a suboptimal control problem. System (1.6) takes the form

3/fe = -k2yfc + vfc, k € Z+. (2.4)

oo

Let Q be the set of all positive integers having the form 32pq, where p = 0,1,2,..., and q is relatively prime with 3. It is obvious that the set of all positive integers Z + is the union of the two disjoint sets Q and 3Q. Then (2.4) can be rewritten in the form

Vk = -k2Vk + vk, V3k = -9k2V3k + v3k, k € Q. (2.5)

After the substitutions

1k 1 l^k 2 + 1 1 2

yk = ' y3k = ' t = £2T' l'k = LlkW ' l'3k = tLkW '

all systems (2.5) will be reformulated to the following two-dimensional control system:

x1 = -x1 + w1, x2 = — 9x2 + w2. (2.6)

Now, following the Chernous'ko way, we replace (2.3) by the even more rigid restriction

max w sin kx + w| sin3kx| < 1, k € Q, (2.7)

0<x<n

that implies (2.3) if ^keg |k = v0. Thus, we have reduced the infinite dimensional control problem to the two-dimensional problem.

3. Solution of the auxiliary time-optimal problem on the plane

Let Pk denote the set of all pairs (w^, ) for which (2.7) holds. Setting

P = \ w = (w1, w2) € R2 : max |w1 sin t + w2 sin3t1 < 1 L I 0<t<n 1 1 J

we have Pk = | kP. As a result, the considered problem of constructing a suboptimal control reduces to the concrete problem of time-optimal control for the following two-dimensional system:

x1 = —x1 + w1, x2 = —9x2 + w2, (w1 ,w2) € P. (3.1)

Obviously, P is a convex and compact set with non-empty interior (i.e., a convex body). Since P is symmetric with respect to the origin, we may restrict ourselves to considering only the case w1 > 0. It is more convenient to set sin t = y. Then, by the formula

sin 3t = 3 sin t — 4 sin3 t,

we get

P = {w = (w1,w2) € R2 : max |(w1 + 3w2) y — 4w2 y3 | < ^. L 0<y<1 J

Just this transformation lay on the base of the separation Z + = Q U 3Q.

After elementary calculations, we find that the part of the boundary of the set P lying in the half-plane w1 > 0 is given by the formula

w1 =

w2 + 1 if —1 < w2 < 0.125,

3( Vw1 - w2) if 0.125 < w2 < 1,

while the other part is found by central symmetry (see Fig. 1).

Let us recall that, in the auxiliary problem (3.1), a unique optimal time-control function exists at each initial point (x0,x0) [7, 23, 27]. The existence follows from the property O € Int P. The

w2

1 KA B /'arctan2 .D 'E w1

. -1 • o "l F '

p

1 H

Figure 1. The straight ray AC is tangent at the point A to the curve AD, which is a part of the boundary of P.

uniqueness is a consequence of the following feature of P: the vector (1,1), which is orthogonal to the segment AC, is not an eigenvector of the matrix of system (3.1). Therefore, the optimal control problem (2.5) coincides with the extremal controls of Pontryagin's maximum principle [23, 30].

To calculate the latter, we prefer to use the "backward motion" principle. Let T(£¿,£0) be a transition time for the initial point (xq,®o) in the system (3.1). If we set t = T(£¿,£0) — t, then extremals of Pontryagin's maximum principle are defined by the system

dx1 ! dx2 2 -2 #i , #2 n/

—— = X - W , —— = 9x - W , —— = -01, —— = 902• (o o\

dT dT dT dT (3.2)

x(0) = y(0) = 0, 01(O)=cos s, 02(O) = sins, —n < s < n.

Since ^¿(t, s) = e-T cos s and 02(t, s) = e-9T sin s, an extremal control W(t, s) should be found by the Pontryagin's maximum principle, i.e., from the equation

ui1(t, s)e-T cos s + W2(t, s)e-9T sin s = max[w1 e-T cos s + w2e-9T sin s]. (3.3)

Equation (3.3) leads to the following construction of the extremal controls.

If 0(t, s) lies in the open angle AOB, then obviously W(t, s) = (0,1). Note that, if s = n/2, then 0i(t, s) = 0. Therefore, w(t, 7t/2) = (0,1). Similarly, if s = 0, then tp2(t,s) = 0; thus, w(t, 0) = (2 v/3/3, V/3/9).

Consider now the dynamics of 0(t,s). In the case 0 < s < n/2, the vector 0(t,s) lies in the quarter > 0, 02 > 0 and turns clockwise. Moreover, its direction tends to the axis of abscissas OE as t ^ (Similarly, if —n/2 < s < 0, then 0(t, s) lies in the quarter > 0, 02 < 0 and turns counterclockwise with the same limit direction.)

Thus, the extremal control has the following structure: if 0 < s < arctan2 (see Fig. 1), then 0(t, s) lies in the angle BOD for all t (t > 0), and, hence, W(t, s) is a point of the arc AD such that its projection to the direction 0(t, s) is maximal (the analytical expression for W(t, s) is given in Table 1).

Table 1. The analytical expression for w1(r, s) and w2(r, s).

w1{t, s) w2(t, s)

7T 7T --< s <--- 2 ~ 4 0 if 0 < t < t**, 3(2 - e"8rtans)M if r„ < r — 1 if 0 < t < r**, M if r„ < r

-7r/4 < s < 0 3(2 — e~8r tan s)M M

s = 0 2y/3/3 y/3/9

0 < s < arctan2 3(2 — e~8r tan s)M M

7t arctan2 < s < — 0 if 0<r<r*, 3(2 - e"8rtans)M if r* < r 1 if 0<r<r*, M if r* < r

s = tt/2 0 1

Table 2. The analytical expression for x1(r, s) and x2(r, s).

ir1(r, s) x2(r, s)

7T 7T --< S <--- 2 ~ 4 0 if 0<r<r„, n 3eT fl~'[)2io if - -4 tan s J p2q7 ** — m (e9r — 1) /9 if 0<r<r„, n e9r f qdp if r** < r 4 tan s ,/ m

—7r/4 < s < 0 n 3eT f 1~ p2 ip 4 tan s J p2q7 m 9 ™ /14- qdP 4 tan s ,/ m

s = 0 2V/3(1 - er)/3 >/3(1 — e9r)/81

0<s< arctan 2 n 3eT f 1~ p2 ip 4 tan s J p2q7 m n e9r f - / qdp 4 tan s J m

7r arctan 2<s<— 0 if 0<t<t*, n 3er i^- P'ip if - / -4 tan s 7 p2g7 * — m (1 — e9r) /9 if 0<r<r*, n e9r f qdp if r* < r 4 tan s ,/ m

s = tt/2 0 (1 - e9r) /9

Further, in the case arctan2 < s < n/2, we have wD(r,s) = (0,1) on the interval [0,r*), where r* = -1/8 ■ log(2cot s). At the time r = r*, the vector ^(r, s) becomes orthogonal to the right side tangent to the curve dP at the point (0,1) and it occurs "switching" of the extremal control from the value (0,1) to a continuous mode. Namely, wD(r, s) begins sliding along the arc AC (see Table 1) and tends to the point C as r ^

Similarly, if (—n/2 <s< —n/4), then wD(r, s) = (0, —1) at 0 < r < r**, where r** = — 1/8 ■ log(— cot s) and w(r, s) is a switching time. On the interval (r**, +rc>), wD(r, s) slides along the arc ED tending to the point D.

The entire synthesis of the extremal control is given in Table 1. Due to the central symmetry, the values of s are considered only on the range —n/2<s<n/2 and the following notation is used:

M = (3 — e-8T tan s)-3/2, m = (3e8T cot s — 1)-1/2, n = (3cot s — 1)-1/2,

p = (3 — e-8T tan s)-1/2, q = ((3 — p-2)cot s)1/8.

Now, extremal trajectories can be easily calculated by (3.2). The corresponding formulas are gathered in Table 2. They are illustrated in Fig. 2.

y

Figure 2. The extremal trajectories.

1.0 0.8 0.6 0.4 0.2

Figure 3. The graphs of the functions ww1(r, s) (the continuous line) and w2(t, s) (the dashed line) for arctan2 < s < n/2.

4. Construction of a suboptimal control in the initial problem

Let us now derive the solution of problem (2.1), (2.3), (2.5) basing one of reduced problems.

Extremals x(t, s), y(t, s) cover R2. Therefore, for every (x0, ) € R2 \ (0,0), there exists a pair (to,so) such that x1(ro,so) = £2(t0,s0) = Further, in the system (3.1), for every fixed (xo,x0) = 0, an optimal control is unique, which implies the uniqueness of the value to (while corresponding values of s0 may be not unique, but one can choose any of them).

Then T(£¿,£0) = t0 is the transition time and

v1(t) = v1 (to - t, so), v2(t) = v2(to - t, so)

1.0

0.5

0.5

T 1.0 1.5 2.0

v i

1 **

-0.5

-1.01-

Figure 4. The graphs of the functions w1(r, s) (the continuous line) and w2(t, s) (the dashed line) for —n/2 < s < arctan2.

is the suboptimal control for (2.1).

Let us now consider system (2.5). For an initial point (v0,v0k), the corresponding trajectory (yk(t),y3k(t)) satisfies the condition

Vk (Tk) = V3k (Tk) = 0,

where

The constructed synthesis implies that Tk is monotonically decreasing in |k, and it is easy to see that Tk ^ 0 as |k ^ and Tk ^ to as |k ^ 0. Therefore, for every k, there exists a unique value ik such that Tk(ik) is the same for all k. Moreover, can be chosen satisfying the condition = vo- One can easily see that

« ^ * ^ P

- ^k S -p

for some positive a and ft.

Finally, we consider the initial problem (2.1), (2.3), (2.5). Let u0(x) = ^u°k sinx be the Fourier expansion of the initial function u0(x). Taking (y°,y0k), k € Q, as an initial point for system (2.6), we find

w°k(t) = —v\{t), wlk[t) = —vltt), k e Q. Ik Ik

Thus, the following statement holds. Theorem 1. The function

<x

v(t, x) = ^ w^0(i) sin kx k=1

is a suboptimal control in problem (2.1) for the initial state u0(x).

5. Conclusion

The paper is devoted to the time-optimal problem for the process of heat conductivity in a rod when the control parameter is the intensity of external heat sources. A suboptimal control is

constructed by the combination of the Chernous'ko approach with the method of grouping terms of the Fourier expansion.

This method may be applied to the time-optimal control problem for other systems given in an evolutionary form.

The following question naturally arises: how effective is the method of grouping? First of all, let us bring general considerations. The set of all admissible controls in the initial problem (1.1)-(1.2) can be identified with the subset

Ulnitial = 1 u € I2 | sup | y] Uk sin kx| < Vo [.

0<x<n , „ — — k=0

As noted in Section 1, Chernous'ko restricted the set of controls using

Uch = {u € I2 | |uk |< Uk, k = 1, 2, 3,...},

where Uk is a sequence chosen from the condition ^ Uk < v0 and guaranteeing the equality u(t, x) = 0 for some T = TCh > 0.

The considerations in this paper are based on the set

Ugr = \ u € 12 | max |uk sin kx + u3kx sin 3kx| < Uk

I 0<x<n

taken as a region of admissible controls. One can easily see that

Ugr C UCh C UInitial.

These relations imply T2 < Ti < T0 for optimal and suboptimal times of transition respectively.

If one takes an initial point of the form (0,0,..., x^, 0,..., 0), i.e., in terms of the initial problem (1.1)-(1.2), <^(x) = (0,0, ...,xm sin mx, 0,..., 0), then, obviously, UInitial = UCh = Ugr and, thus, T2 = T1 = T0. But if an initial point is taken in the form (0,0,..., 0, x°k, 0,..., 0, x0k, 0,...), then UIn=Ugr while UInDUCh and, thus, T2 = T0 < T1. Table 3 contains values for specific cases.

Table 3

Initial point Ti T0 = T2

(1,0,1,0,0,0,...) 0.7 0.64

(1,0,2,0,0,0,...) 0.67 0.61

(1,0,-1,0,0,0,...) 0.72 0.65

(1,0,-2,0,0,0,...) 0.69 0.62

Obviously, xk = 0 at least for three values of the index k when T2 < T0 < Ti. The final note is that the method of grouping can be applied only if there some algebraic relations between the eigenfunctions of the operator A.

REFERENCES

1. Ahmed N. U. Optimal control of infinite dimensional systems governed by integro differential equations. In: Differential Equations: Dynamical Systems, and Control Science. Eds. K.D. Elworthy, W.N. Everitt, E.B. Lee. Ser. Lect. Notes Pure Appl. Math., 1994. Vol. 152. P. 383-402.

2. Alekseev V.M., Tikhomirov V.M., Fomin S.V. Optimal'noe upravlenie [Optimal Control]. Moscow: Nauka, 2005. 384 p. (in Russian)

3. Arutyunov A. V., Karamzin D. Y., Pereira F. M. The maximum principle for optimal control problems with state constraints by R.V. Gamkrelidze: Revisited. J. Optim. Theory Appl., 2011. Vol. 149, No. 3. P. 474-493. DOI: 10.1007/s10957-011-9807-5

4. Arutyunov A. V., Vinter R.B. A simple 'Finite Approximations' proof of the Pontryagin maximum principle under reduced differentiability hypotheses. Set-Valued Anal., 2004. Vol. 12, No. 1-2. P. 5-24. DOI: 10.1023/B:SVAN.0000023406.16145.a8

5. Azamov A. A., Ruzibayev M. R. The time-optimal problem for evolutionary partial differential equations. J. Appl. Math. Mech., 2013. Vol. 77, No. 2. P. 220-224. DOI: 10.1016/j.jappmathmech.2013.07.013

6. Barbu V. Analysis and Control of Nonlinear Infinite Dimensional Systems. Ser. Math. Sci. Eng., Vol. 190. Boston: Academic Press, MA, 1993. 475 p.

7. Blagodatskikh V.I. Vvedenie v optimal'noe upravlenie (linejnaya teoriya). [Introduction to Optimal Control (Linear Theory)]. Moscow: Vysshaya Shkola, 2001. 239 p. (in Russian)

8. Bongini F., Fornasier M., Rossi F., Solombrino F. Mean-field Pontryagin maximum principle. J. Optim. Theory Appl., 2017. Vol. 175, No. 1. P. 1-38. DOI: 10.1007/s10957-017-1149-5

9. Bonnet B., Rossi F. A Pontryagin Maximum Principle in Wasserstein Spaces for Constrained Optimal Control Problems. 2018. arXiv: 1810.13117v2 [math.OC]

10. Bryson A. E. Optimal Control — 1950 to 1985. IEEE Control Systems, 1996. Vol. 16, No. 3. P. 26-33. DOI: 10.1109/37.506395

11. Butkovskiy A. G. Metody upravleniya sistemami s raspredelennymi parametrami [Control Methods of the Systems with Distributed Parameters]. Moscow: Nauka, 1965. 474 p. (in Russian)

12. Carthel C., Glowinski R., Lions J.L. On exact and approximate boundary controllability for the heat equation: A numerical approach. J. Optim. Theory Appl., 1994. Vol. 82, No. 3. P. 429-484. DOI: 10.1007/BF02192213

13. Chernous'ko F.L. Bounded controls in distributed-parameter systems. J. Appl. Math. Mech., 1992. Vol. 56, No. 5. P. 707-723. DOI: 10.1016/0021-8928(92)90057-F

14. Evans L.C. Partial Differential Equations. 2-nd ed. Ser. Grad. Stud. Math., Vol. 19. Providence, Rhode Island: Amer. Math. Soc., 2010. 749 p.

15. Gong W., Hinze M., Zhou Z. A priori error analysis for finite element approximation of parabolic optimal control problems with pointwise control. SIAM J. Control Optim., 2014. Vol. 52, No. 1. P. 97-119. DOI: 10.1137/110840133

16. Ibragimov G., Risman M. H., Azamov A. A. Existence and uniqueness of the solution for an infinite system of differential equations. J. Karya Asli Lorekan Ahli Matematik, 2008. Vol. 1, No. 2. P. 9-14.

17. Ibragimov G. I. Optimal pursuit time for a differential game in the Hilbert space l2. Science Asia, 2013. Vol. 39S, No. 1. P. 25-30. DOI: 10.2306/scienceasia1513-1874.2013.39S.025

18. Ji G., Martin C. Optimal boundary control of the heat equation with target function at terminal time. Appl. Math. Comput., 2002. Vol. 127, No. 2-3. P. 335-345. DOI: 10.1016/S0096-3003(01)00011-X

19. Krasovskii N. N. Teoriya upravleniya dvizheniem [Theory of Control of Motion]. Moscow: Nauka, 1968. 476 p. (in Russian)

20. Kubyshkin V. A., Postnov S. S. Time-optimal boundary control for systems defined by a fractional order diffusion equation. Autom. Remote Control, 2018. Vol. 79, No. 5. P. 884-896. DOI: 10.1134/S0005117918050090

21. Ladyzhenskaya O. A. The boundary value problems of mathematical physics. Ser. Appl. Math. Sci., vol. 49. New York: Springer-Verlag, 1985. 322 p. DOI: 10.1007/978-1-4757-4317-3

22. Laykekhman D., Vexler B. Optimal a priori error estimates of parabolic optimal control problems with pointwise control. SIAM J. Numer. Anal., 2013. Vol. 51, No. 5. P. 2797-2821. DOI: 10.1137/120885772

23. Lee E.B., Markus L. Foundation of Optimal Control Theory. Malabar (FL): Krieger Pub Co, 1986. 586 p.

24. Lee M. J., Park J. Y. Pontryagin's maximum principle for optimal control of a non-well-posed parabolic differential equation involving a state constraint. ANZIAM J., 2004. Vol. 46, No. 2. P. 171-184. DOI: 10.1017/S1446181100013778

25. Lions J. L. Optimal Control of Systems Governed by Partial Differential Equations. Ser. Grundlehren Math. Wiss., Vol. 170. Berlin, Heidelberg: Springer-Verlag, 1971. 440 p.

26. Lii Q., Wang G. On the existence of time optimal controls with constraints of the rectangular type for heat equations. SIAM J. Control Optim., 2011. Vol. 49, No. 3. P. 1124-1149. DOI: 10.1137/10081277X

27. Magaril-Ilyaev G. G., Tikhomirov V. M. Convex Analysis: Theory and Applications. Ser. Transl. Math. Monogr. Providence, USA: Amer. Math. Soc., 2003. 183 p.

28. Mizohata S. The Theory of Partial Differential Equations. Cambridge: University Press, 1979. 490 p.

29. Pan L. P., Yong J. Optimal control for quasilinear retarded parabolic systems. ANZIAM J., 2001. Vol. 42, No. 4. P. 532-551. DOI: 10.1017/S1446181100012268

30. Pontryagin L. S., Boltyanskii V. G., Gamkrelidze R.V., Mishchenko E.F. The Mathematical Theory of Optimal Processes. NY, London: John Wiley & Sons, 1962. 320 p.

31. Raymond J. P., Zidani H. Pontryagin's Principle for state-constrained control problems governed by parabolic equations with unbounded controls. SIAM J. Control Optim., 1998. Vol. 36, No. 6. P. 18531879. DOI: 10.1137/S0363012996302470

32. Raymond J. P., Zidani H. Pontryagin's principle for time-optimal problems. J. Optim. Theory Appl., 1999. Vol. 101, No. 2. P. 375-402. DOI: 10.1023/A:1021793611520

33. Ross I. M. A Primer on Pontryagin's Principle in Optimal Control. Collegiate Publishers, 2009. 109 p.

34. Serag H. M. Distributed control for cooperative systems involving parabolic operators with an infinite number of variables. Pure Math. Appl., 2004. Vol. 15, No. 4. P. 439-451.

35. Tsachev T. An optimal control problem for the heat equation. Mathematica Balkanica, 1984. Vol. 3. P. 296-310. URL: http://www.math.bas.bg/infres/MathBalk/MB-03/MB-03-296-310.pdf

36. Zhang Y. On a kind of time optimal control problem of the heat equation. Adv. Difference Equ., 2018. Vol. 2018. Art. no. 117. P. 1-10. DOI: 10.1186/s13662-018-1577-z