CC BY
11
Vestnik ^AUNC. Fiz.-Mat. nauki. 2024. vol. 47. no. 2. P. 9-20. ISSN 2079-6641
MATHEMATICS
" https://doi.org/10.26117/2079-6641-2024-47-2-9-20 Research Article Full text in English MSC 35K05, 35K15
The Control Problem for a Heat Conduction Equation with Neumann Boundary Condition
F. N. Dekhkonov*
Namangan State University, 316, Uychi str., 160136, Namangan, Uzbekistan
Abstract. Previously, boundary control problems for a heat conduction equation with Dirichlet boundary condition were studied in a bounded domain. In this paper, we consider the boundary control problem for the heat conduction equation with Neumann boundary condition in a bounded one-dimensional domain. On the part of the border of the considered domain, the value of the solution with control parameter is given. Restrictions on the control are given in such a way that the average value of the solution in some part of the considered domain gets a given value. The studied initial boundary value problem is reduced to the Volterra integral equation of the first type using the method of separation of variables. It is known that the solution of Volterra's integral equation of the first kind cannot always be shown to exist. In our work, the existence of a solution to the Volterra integral equation of the first kind is shown using the method of Laplace transform. For this, the necessary estimates for the kernel of the integral equation were found. Finally, the admissibility of the control function is proved.
Key words: parabolic equation, integral equation, initial-boundary problem, admissible control, Laplace transform.
Received: 18.04.2024; Revised: 17.05.2024; Accepted: 08.06.2024; First online: 25.08.2024
For citation. Dekhkonov F. N. The control problem for a heat conduction equation with Neumann boundary condition. Vestmh KRAUNC. Fiz.-mat. nauki. 2024,47: 2,9-20. EDN: MNMAFB. https://doi.org/10.26117/2079-6641-2024-47-2-9-20.
Funding. The work was not carried out within the framework of funds
Competing interests. There are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. The author participated in the writing of the article and is fully responsible for submitting the final version of the article to the press.
* Correspondence: A E-mail: [email protected]
The content is published under the terms of the Creative Commons Attribution 4.0 International License © Dekhkonov F.N., 2024
© Institute of Cosmophysical Research and Radio Wave Propagation, 2024 (original layout, design, compilation)
Вестник КРАУНЦ. Физ.-мат. науки. 2024. Т. 47. №2. C.9-20. ISSN 2079-6641
МАТЕМАТИКА
" https://doi.org/10.26117/2079-6641-2024-47-2-9-20
Научная статья
Полный текст на английском языке
УДК 519.65
Задача управления для уравнения теплопроводности с граничным условием Неймана
Ф. Н. Дехконов*
Наманганский государственный университет, ул. Уйчи, 316, 160136, г. Наманган, Республика
Узбекистан
Аннотация. Ранее были исследованы задачи граничного управления для уравнения теплопроводности с граничным условием Дирихле в ограниченной области. В данной работе рассматривается задача граничного управления для уравнения теплопроводности с граничными условиями Неймана в ограниченной одномерной области. На части границы рассматриваемой области задано значение решения с управляющим параметром. Ограничения на управление задаются таким образом, чтобы среднее значение решения в некоторой части рассматриваемой области получало заданное значение. Исследуемая начально-краевая задача сводится к интегральному уравнению Вольтерра первого типа с использованием метода разделения переменных. Известно, что не всегда можно доказать существование решения интегрального уравнения Вольтерра первого рода. В нашей работе существование решения интегрального уравнения Вольтерра первого рода показано с помощью метода преобразования Лапласа. Для этого были найдены необходимые оценки ядра интегрального уравнения. Наконец, допустимость функции управления доказана.
Ключевые слова: параболическое уравнение, интегральное уравнение, начально-краевая задача, допустимое управление, преобразование Лапласа.
Получение: 18.04.2024; Исправление: 17.05.2024; Принятие: 08.06.2024; Публикация онлайн: 25.08.2024
Для цитирования. Dekhkonov F. N. The control problem for a heat conduction equation with Neumann boundary condition // Вестник КРАУНЦ. Физ.-мат. науки. 2024. Т. 47. № 2. C. 9-20. EDN: MNMAFB. https://doi.org/10.26117/2079-6641-2024-47-2-9-20. Финансирование. Работа не выполнялась в рамках фондов.
Конкурирующие интересы. Конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответственность. Автор участвовал в написании статьи и полностью несет ответственность
за предоставление окончательной версии статьи в печать.
* Корреспонденция: А E-mail: [email protected] Hgk ф-
Контент публикуется на условиях Creative Commons Attribution 4.0 International License © Dekhkonov F. N., 2024
© ИКИР ДВО РАН, 2024 (оригинал-макет, дизайн, составление)
1 Introduction
In this article, we consider the following heat conduction equation in the bounded domain O = {(x, t) : 0 < x < I, t > 0}:
du(x,t) = » fk(x) ), (x,t) 6 O, (1)
9t 9x \ 9x
with Neumann boundary conditions
ux(0,t) = -^(t), ux(l,t)= 0, t>0, (2)
and initial condition
u(x, 0) = 0, 0 < x < l, (3)
where |x(t) is control function.
Assume that the function k(x) G C2([0,l]) satisfies conditions
k(x) > ko > 0, k'(x) < 0, 0 < x < l. (4)
Definition 1. It is called that the function |x(t) G W](M+) is admissible control, if it fulfills the conditions p.(0) = 0 and |p.(t)| < 1 for all t > 0.
Control Problem. For the given function 0(t) Problem consists looking for the admissible control |x(t) such that the solution u(x, t) of the initial-boundary problem (1)-(3) exists and for all t > 0 satisfies the equation
u(x,t) dx = 0(t). (5)
Control problems for parabolic equations were first studied in [1, 2]. Control problems for the infinite-dimensional case were studied by Egorov [3], who generalized Pontryagin's maximum principle to a class of equations in Banach space, and the proof of a bang-bang principle was shown in the particular conditions.
The optimal time problem for second-order parabolic type equation in the bounded n—dimensional domain was studied in [4,5] and the optimal time estimate for achieving a given average temperature was found. The control problem for the heat equation associated with the Neumann boundary condition in a bounded three-dimensional domain is studied in [6]. In this work, an estimate of the optimal time was found when the average temperature is close to the critical value.
In [7,8], the control problems of the heat equation associated with the Dirichlet boundary condition in the two-dimensional domain are studied. In these articles, an estimate of the minimum time for achieving a given average temperature was found, and the existence of a control function is proved by the Laplace transform method. The boundary control problem related to the fast heating of the thin rod for the inhomogeneous heat conduction equation was studied in works [9] and the existence of the admissible control function was proved.
The minimal time problem for the heat conduction equation with the Neumann boundary condition in a one-dimensional domain is studied in [10]. The difference of this work from the previous works is that the required estimate for the minimum time is found with a non-negative definite weight function under the integral condition. In [11], the control problem for a second-order parabolic type equation with two control functions was studied and the existence of admissible control functions was proved by the Laplace transform method.
Boundary control problems for parabolic type equations are also studied in works [12-14].
A lot of information on the optimal control problems was given in detail in the monographs of Lions and Fursikov [15,16]. General numerical optimization and optimal boundary control have been studied in a great number of publications such as [17]. The practical approaches to optimal control of the heat equation are described in publications like [18].
In this work, the boundary control problem for the heat transfer equation is considered. The difference of this work from the previous works is that in this problem, the control problem for the heat conduction equation related to the Neumann boundary condition is studied. In Section 2, the boundary control problem studied in this work is reduced to the Volterra integral equation of the first kind by the Fourier method. In Section 3, the solution of Volterra's integral equation is proved using the Laplace transform method.
2 Main integral equation
Consider the following spectral problem
AkVkM= 0, 0 < x < I, (6)
with boundary condition
vk(0)= vkd)= 0, 0 < x < I. (7)
It is well-know that this problem is self-adjoint in L2(O) and there exists a sequence of eigenvalues {Ak} so that Ai < A2 < ... < Ak —> oo, k —> oo.
The corresponding eigenfuction vk form a complete orthonormal system {vk} in L2(O) and these function belong to C(O), where O = O U 9O (see [19,20]).
Definition 2. By the solution of the problem (1)-(3) we understand the function u(x, t) represented in the form
(l _ x)2
u(x,t) = ^^(t)_ w(x,t), (8)
where the function w(x, t) G CX,t(O) n C(O), wx G C(O) is the solution to the problem:
, 9w\ d / . , I — x\ r . (I — x)2 . ..
wt = 5X (k(x) dw ) + dx (k(x) -_-) ^(t) + ^ ^ '(t)> (9)
with boundary value conditions
Wx(0,t)= 0, Wx (l,t)= 0,
and initial value condition
w(x,0) = 0.
We set
Pk = (Akbk - ak) Ck, k = 1,2,..., where coefficients ak, bk and ck are as follows
(10) (11) (12)
=
dX fk(x) l-j-x ) vk(x)dx, bk =
(I - x)2 2l
Vk (x)dx,
(13)
and
Ck =
Vk(x)dx.
We understand the coefficients a0 and b0 as follows
(14)
ao =
A(k(x) LpW = - if-,
and
Thus, we have
bo =
(I - x)2 d I2 = 6.
( ^ l2 m k(0) w(x,t) = 12 M-(t) —2~
^(s) ds+
k=1
+ Z e-Ak(t-s) (^(s) ak + ^'(s) bk) ds ) vk(x),
(15)
where ak and bk defined by (13).
From (8) and (15), we get the solution of the mixed problem (1)-(3) (see [19]):
, > (l - x)2 , , l2 , k(0)
u(x, t) = 2j ; ^(t) -— ^(t) +
|x(s) ds-
- Z. { [ e-Ak(t-s) (^(s) ak + M-'(s) bk) ds^ vk(x) k=1 '
We know that the eigenvalues Ak of the boundary value problem (6)-(7) satisfies the following inequalities
Ak > 0, k = 0,1,... (16)
Indeed, since
then we get
dx(kWdVxr) + AkVk(x) = o, 0<x<l,
Ak =-
-MkW^ ) vk(x)dx =
k(x) |vk(x)|2dx > 0.
According to condition (5) and the solution of the problem (1)-(3), we may write
e(t) =
u(x,t)dx = |(t)
1 (l - x)2 , l^^k(0) lt -^dx - 12|(t) + ~2~
l(s) ds-
k=1
- W e-Ak(t-s) (|(sK + |'(s) bk) ds
Vk(x) dx =
= l(t)
l (l-x)^ l^^k(0)lt -^dx - 12|(t) + ~2~
| (s) ds-
k=1
-Y_ ak Ck e-Ak(t-s) |(s)ds -Y_ bk Ck e-Ak(t-s)|'(s)ds =
k=1
= l(t)
l (l-x)2 ^ l3 m ^(0)lt
-^dx -12|(t) + 2""
| (s) ds-
0
oo
- Y ak Ck k=1
e-Ak(t-s)|(s)ds - |(t)^ bkCk+
k=1
+ y Ak bk Ck
k=1
e
-Ak(t-s)
l(s)ds,
(17)
where ck defined by (14). Note that
(l - x)2 , l3 ~
2l
dx = + bk Ck,
k=1
(18)
where bk, Ck are defined by (13) and (14).
As a result, from (17) and (18), we obtain
0(t) =
k(0) I 2
+ (Akbk - flk) cke-Ak(t-sM ^(s)ds.
k=1
We set
k(0) I
B(t) = ^ ßk e-Akt, t > 0.
k=1
where |3k defined by (12).
Then we get the main integral equation
(19)
B(t - s) p.(s)ds = 0(t), t>0.
(20)
Lemma 1. For the coefficients {pk}k=i the following estimate is valid:
0 < Pk < C, k = 1,2,...,
where C is a positive constant.
Proof. First we calculate the following equality using (13)
Akbk =
(I - x)2 21
Ak Vk(x)dx = -
(l - x)2d^k(x) | dx =
21 dx
dx
(I - x)2 21
k(x) vk (x)
x=l
+
x=0
(l ^ x) k(x) vk(x)dx )
(I - x)
k(x) vk (x)dx = k(0)vk(0) +
dx mx) ) Vk(x)dx
= k(0)vk(0) + ak.
Then we have
Akbk — ak = k(0)vk(0), k = 1,2,.... We know that the following inequality is true (see [20])
(21)
Vk(0)
vk(x) dx > 0, k = 1,2,....
Thus, by (21) and (22), we have
ßk = (Akbk - ak) Ck = k(0) Vk(0)
vk(x)dx > 0.
I
It is clear that if k(x) e C1([0,l]), we may write the estimate (see [21,22])
max |vk(x)| < C.
0<x<l
From this we can obtain the following estimates
Pk < k(0) |vk(0) Ck| < C.
□
Lemma 2. Let 1/2 < a < 1. Then for the function B(t) defined by (19) the following estimate is valid:
C
0<B(t) < —, 0<t < 1,
where Ca is a constant depending only on a.
Proof. It is known from the general theory that if k(x) is a smooth function, the following estimate is valid (see [22]):
i
k2n2 „ _2 dx
Ak = + O(k ), p =
P
x/kW'
Let 1/2 < a < 1 and À > 0. Then the maximum value of the function h(t, A) = tae—At is reached at the point t = a and this value is equal to aae—a. As a result, for any 1/2 < a < 1, we get the estimate
a
B(t) < const - ^ Pktae-Akt < C ^ < CP
k=1 k=1 k
where
1
L Aa< +oc • k=1 Ak
□
3 Main result
In this section, we prove the existence of the control function.
22
Denote by W(M) the set of function 0 G W2 (—oo, +oo), which satisfies the condition
l|e||w|(R+) < m 0(t) = 0 for t < 0.
Theorem 1. There exists M > 0 such that for any function 0 e W(M) the solution |(t) of the equation (20) exists, belongs to C(R+) and satisfies condition
ll(t)| < 1.
We use the Laplace transform method to solve the integral equation (20). It is known
M-(p) =
e—pt B(t) dt.
Then we use Laplace transform obtain the following equation
0(p) =
e—pt dt
B(t - s)^(s)ds = B(p) b(p).
Thus, we get
B(p) = ~(p), where p = a + i£, a > 0, £ e R, B(p)
and
a+iTO __+00
B(t)= ' [ JBiPl ePtdp = 1
B( ) 2ni J B(p) P 2n
a—iw —00
B(a + ^ g(a+i£)td^ B(a + i£) .
Then we can write
B (p) =
B(t)e—pt dt =
k(0) I
2
k(0) I 1
e—pt dt + Y_ ßk
k=1
-(p+Ak)t dt =
+L
2 p k=rp +At'
where B(t) defined by (19) and
(23)
B (a + *) = k^ ' + L
2 a + f—■ a + Ak + i£, k=1
k(0) I a + £ ßk (a + Ak)
2 a2 + k=T (a + Ak)2 +
k(0) I i£,
" 2 a2 +
^ L
ßk
k=1
(a + Ak)2 +
= ReB(a + i£) + i ImB (a + i£),
where
nBr ^-n k(0) l_a ^ ßk (a + Ak)
ReB(a + i^) = —-— 2 , + )
2 a2 + k=T (a + Ak)2 +
ImB (a + i£) =
k(0)l
* L
2 a2 + "k=T (a + Ak)2 +
We know that
(a + Ak)2 + < [(a + Ak)2 + 1](1 + ^2),
e
ß
k
ß
k
ß
k
and we get the following inequalities
1
>
a2 + " 1 + 1 + a2'
and
1
>
1
1
(24)
(25)
(a + Ak)2 + " 1 + (a + Ak)2 + 1" Consequently, according to inequalities (24) and (25) we can obtain the following
estimates
|ReB (a + i£)| =
k(0) I a + j ßk (a + Ak) >
2 a2 + k=T (a + Ak)2 +
C
1a
1 /k(0) I a + ^ ßk (a + Ak) _ > 1 + £2 V 2 1 + a2 + k=7 (a + Ak)2 + 1) 1 +
and
|ImB (a + i£)| = |£,|
k(0) I 1
+ j
ßk
>
2 a2 + k"T (a + Ak)2 + |£| (k(0) I 1 , ^ ßk ^ C2a |£|
>
1 + £,2 V 2 1 + a2
+
^ (a + Ak)2 + 1) 1 + £,2'
where C1a, C2a as follows
r k(0) I a , ^ ßk (a + Ak)
C1a = -^-—-T +
2 1 + a2 ' (a + Ak)2 + 1'
and
C2a =
k(0) I _l_
2 1 + a2
+ j
ßk
k=1 (a + Ak)2 + T
From (26) and (27), we have the following estimate
|B(a + |2 = |ReB(a + |2 + |ImB(a + |2 > min(C2a, C2a)
and
|B(a + i£)| >
Ca
v/ttf
Then, when a —> 0 from (23), we obtain
+oo
1
1 +
, where Ca = min(C1a, C2a).
^(t) = 2n
0B™ e^d^.
B (i£)
(26)
(27)
(28)
(29)
Lemma 3. [9] Assume that 0(t) e W(M). Then for the image of the function 0(t) the following inequality
+oo
|9(i^)|V1 + d^ < C ||0||
W|(R+),
1
1
is valid.
Now we present the proof of the Theorem 1. Proof.
Now, we show that e W](R+). Indeed, according to (28) and (29), we obtain
+oo
+oo
|B(*)|2(1 + |^|2) d^ =
B (i£,)
(1 + №2) d^ <
< C
| B(i^)|2(1 +| 2)2d^ = C|9|W2(R).
Further,
|((t)- ((s)| =
('(t) dt
< Ml2Vt—
Hence, e Lip a, where a = 1/2. From (28), (29) and Lemma 3, we can write
where
□
+oo
((t)| < 2n
<
|B(i^)|^ 1
~-d^ < -——
|B(i£,)| < 2nC
C
+oo
0 J
|B(i£,)|V 1 + d^ <
2nCo
CM
I|0|w|(r+) < cm = 1.
2nCo
M =
2nCo C
4 Conclusion
In this paper, we have considered the boundary control problem for a parabolic-type equation in a one-dimensional bounded domain. By the method of separation of variables, the control problem was reduced to the Volterra integral equation of the first kind. Using the Laplace transform method, the existence of a solution to the integral equation was found and the admissibility of the control function was proved.
References
1. Fattorini H. O. Time-optimal control of solutions of operational differential equations, SIAM J. Control, 1964. vol. 2, pp. 49-65.
2. Friedman A. Optimal control for parabolic equations, J. Math. Anal. Appl., 1967. vol. 18, pp. 479491.
3. Egorov Yu.V. Optimal control in Banach spaces, Dokl. Akad. Nauk SSSR, 1963. vol. 150, pp. 241244 (In Russian).
4. Albeverio S., Alimov Sh. A. On a time-optimal control problem associated with the heat exchange process, Applied Mathematics and Optimization, 2008. vol.57, no. 1, pp. 58-68 DOI: 10.1007/s00245-007-9008-7.
5. Alimov Sh.A., Dekhkonov F.N. On the time-optimal control of the heat exchange process, Uzbek Mathematical Journal, 2019. no. 2, pp. 4-17.
2
—oo
—oo
—oo
6. Dekhkonov F. N. On the control problem associated with the heating process in the bounded domain, Vestnik KRAUNC. Fiz.-mat. nauki, 2022. vol.39, no. 2, pp. 20-31 DOI: 10.26117/20796641-2022-39-2-20-31.
7. Fayazova Z. K. Boundary control of the heat transfer process in the space, Russian Mathematics (Izvestiya VUZ. Matematika), 2019. vol.63, no. 12, pp. 71-79.
8. Dekhkonov F. N. On a time-optimal control of thermal processes in a boundary value problem, Lobachevskii journal of mathematics, 2022. vol.43, no. 1, pp. 192-198 D0I:10.1134/S1995080222040096.
9. Dekhkonov F. N., Kuchkorov E. I. On the time-optimal control problem associated with the heating process of a thin rod, Lobachevskii. J. Math., 2023. vol.44, no.3, pp. 1134-1144.
10. Dekhkonov F. N. On the time-optimal control problem for a heat equation, Bulletin of the Karaganda University. Mathematics Series, 2023. vol.111, no.3, pp. 28-38 DOI:10.31489/2023m3/28-38.
11. Dekhkonov F.N. Control problem concerned with the process of heating a thin plate, Vestnik KRAUNC. Fiz.-mat. nauki, 2023. vol.42, no. 1, pp. 69-79 DOI: 10.26117/2079-6641-2023-42-1-6979.
12. Chen N., Wang Y., Yang D. Time-varying bang-bang property of time optimal controls for heat equation and its applications, Syst. Control Lett, 2018. vol. 112, pp. 18-23.
13. Fattorini H. O. Time and norm optimal controls: a survey of recent results and open problems, Acta Math. Sci. Ser. B Engl. Ed., 2011. vol.31, pp. 2203-2218.
14. Fayazova Z. K. Boundary control for a Psevdo-Parabolic equation, Mathematical notes of NEFU, 2018. vol. 25, no. 2, pp. 40-45.
15. Lions J. L. Contrôle optimal de systèmes gouvernés par deséquations aux dérivées partielles. Dunod Gauthier-Villars: Paris, 1968.
16. Fursikov A.V. Optimal Control of Distributed Systems. Providence, Rhode Island: Math. Soc., 2000.
17. Altmüller A, Grüne L. Distributed and boundary model predictive control for the heat equation, Technical report, University of Bayreuth, Department of Mathematics, 2012.
18. Dubljevic S., Christofides P. D . Predictive control of parabolic PDEs with boundary control actuation: Chemical Engineering Science, 2006.
19. Tikhonov A.N., Samarsky A.A. Equations of Mathematical Physics. Moscow: Nauka, 1966 (In Russian).
20. Naimark M. A. Linear differential operators. Nauka: Moscow, 1962 (In Russian).
21. Yakubov V. Ya. Boundedness of normalized eigenfunctions for the Sturm-Liouville problem with minimal constraints on the smoothness of the coefficients, Differ. Equ., 1994. vol. 30, no. 8, pp. 13611364.
22. Vladykina V. E. Spectral characteristics of the Sturm-Liouville operator under minimal restrictions on smoothness of coefficients, Vestnik Moskov. Univ. Ser. 1. Mat. Mekh., 2019. vol.6, pp. 23-28 (In Russian).
Information about the author
Dekhkonov Farrukhjon Nuriddin oglifa - Ph.D. (Phys. & Math.), Namangan State University, Namangan, Uzbekistan, ORCID 0000-0003-4747-8557.
Информация об авторе
Дехконов Фаррухжон Нуриддин кандидат физико-математических наук, государственный университет, г. Наманган, СЖСГО 0000-0003-4747-8557.
угли& -Наманганский Узбекистан,
