DOI: 10.14529/mmph190403
OPTIMAL CONTROL OF SOLUTIONS TO THE INITIAL-FINAL PROBLEM FOR THE MODEL OF LINEAR WAVES IN A PLASMA
A.A. Zamyshlyaeva, O.N. Tsyplenkova
South Ural State University, Chelyabinsk, Russian Federation E-mail: zamyshliaevaaa@susu.ru
The optimal control problem for a Sobolev type equation of higher order with a relatively polynomially bounded operator pencil is investigated in the paper. The results are applied to the study of the optimal control of solutions to the initial-final problem for the model of linear waves in plasma. The first results on the investigation of equation that describes the linear ion-acoustic waves in an unmagnetized plasma and on the study of some properties of these waves were obtained by Yu.D. Pletner. The initial-final conditions posed for the fourth-order Sobolev type equation are the generalization of the conditions in the Cauchy problem that is unsolvable at the arbitrary initial values. The research is based on the phase space method developed by G.A. Sviridiuk and the theory of relatively polynomially bounded operator pencil developed by A.A. Zamyshlyaeva. The article considers an equation that describes ion-acoustic waves in a plasma in an external magnetic field.
Keywords: Sobolev type equations of higher order with a relatively polynomially bounded operator pencil; model of linear waves in a plasma; optimal control problem; initial-final conditions.
Introduction
Let Q = following problem:
Let Q = (0,a) x (0,b) x (0, c) c R3. The article investigates the optimal control of solutions to the
82 x(s t)
(A - A)xm (s, t) = (A - A)xtt (s, t) + a-+ u(s, t), s eQ, t e (0, r), (1)
8S32
x(s, t) = 0, (s, t) e 8Q X (0, r).
Model (1), (2) describes ion-acoustic waves in plasma in an external magnetic field [1]. The parameters in equation (1) relate such physical quantities as the ionic hydro frequency, the Langmuir frequency, and the Debye radius. The function x(s, t) represents the generalized potential of an electric field, the function u(s, t) represents an external effect.
Problem (1), (2) is investigated in the framework of the theory of relatively polynomially bounded pencils of operators [2]. Consider a high-order abstract Sobolev type equation
Ax(n) = Bn_lX(n~l) +... + B0x + y + Cu, (3)
where the operators A, Bn_u...,B0 eL(X;Y\ CeL(U;Y), the functions M:[0,r)cS+^[/, y: [0, r) c R+ -» Y(t < oo), X, Y, U are Hilbert spaces.
Supplement equation (3) with initial-final conditions [3]
Pn (x№)(0) - x0 ) = 0,
. ^ , k = 0, n -1, (4)
Pfin ( x<k)(r) - xrk ) = 0,
where Pin(^ are some projectors in space X . Thus, the optimal control problem is to find a pair (x,u), where x is a solution of (1), (2), and u e Uad is the control for which the relation
J(xu) = min(x,u)eXXUad J(x,u) (5)
holds.
Here J(x, u) is some specially constructed penalty functional, and Uad is a closed convex set in the space of controls U.
Zamyshlyaeva A.A., Optimal Control of Solutions to the Initial-Final Problem
Tsyplenkova O.N. for the Model of Linear Waves in a Plasma
Many non-classical models of mathematical physics [4-9] are based on Sobolev type equations. For example, they occur in problems of hydromechanics, plasma physics, atmospheric physics, filtration theory, theory of electrical circuits, and others. In the work, to find the optimal control of solutions to linear Sobolev type equations of high-order, the ideas and methods obtained by G.A. Sviridyuk [10] and his students [11-13] in the study of first-order Sobolev type equations are used. Here the initial-final problem [14] is investigated. A distinctive feature of this problem is that one projection of the solution is specified at the initial moment of time, and the other at the final point. The initial-final problem for the first-order Sobolev type equations was considered by G.A. Sviridyuk and S.A. Zagrebina.
1. Polynomially A -bounded operator pencils and projectors. Strong solutions
By B denote the pencil formed by operators Bn_x,..., B0. The sets p '(B) = {// e(': (p1 A --pn~xBn_x -... - juB1 -B()f] eL(Y;X)} and crA(B) = C\ pA(B) are called an A - resolvent set and an A -spcctrum of pencil B. respectively. The operator function
RA (B) = (¿u"A - Hn_A -... - pB} -Boyl of a complex variable with domain p 1 (B) is called an A -
resolvent of pencil B .
Let the pencil be B polynomially A -bounded. Introduce an additional condition [15]
Ill't(B)dp = (), k = 0,n-2, (6)
where y is a contour that bounds the domain containing the relative spectrum of the pencil B . Then the operators
P = ^-\Rj(B)lun-lAdlu, Q = —!—[ pn~lARA(B)dp
2mh^ 2mh ^
are projectors in spaces X and Y . Denote X° = kerP, 7° = ker Q, X1 = imP and 71 = imQ . If the pencil of operators is polynomially A -bounded, oo is a pole of order e {OJ U A'' of the A -resolvent of pencil B, then the pencil of operators B is called (A, p) -bounded.
Let the pencil B be polynomially A -bounded, (6) hold and the following conditions be satisfied:
^ (3) = rf (3) U erf (3), <4 (3) *0,k = oTl; and there is a circuit rn c C,
, (7)
bounding the domain rocC such, that
r0 n CJA (3) = CJA (B), f 0 n df (B) = 0.
f ¿umRf(B)du = 0,m = 0,n-2 (8)
Operators
Pfin=^-\ju^(B)AdjueL(X)
r0
and Pin = P - Pjfin e L(X) are projectors in space X [15].
Definition 1. The vector-function x eHn(X)={x eL2(0,r;X): x(n) eL2(0,r;X)} is called a strong solution of linear non homogeneous Sobolev type equation
Ax™=Bn_1x<n-»+... + B0x + y, (9)
if it turns the equation to an identity almost everywhere on interval (0, r). A strong solution x = x(t) of (9) is called a strong solution to (4), (9), if (4) holds.
Theorem 1. If the pencil B is (A, p) -bounded, p e {()} U A''. conditions (6)-(8) are satisfied, then for arbitrary x°, xrk e X, k = 0, n -1 and y e Hp+n (Y) there exists a unique strong solution to (4) for equation (9).
Математика
2. Optimal control for the model of linear waves in plazma
Consider problem (4) for equation (3), where the functions x, y, u lie in X, Y and U, respectively.
Introduce the control space
p+n
H (U) = {MeI2(0,r;U) : uip+n) e L2(0,r;U),uÇq)(0) = 0,q = 0,p}, p e {0} U N. It is a Hilbert space with inner product
[V, w]= v(q\ w(q dt.
q=0
p+n
In the space H (I/) single out a closed convex subset Uad. A vector function u e Uad is called an optimal control of solutions to (3), (4), if (5) holds.
Let us prove the existence of optimal control U e Uad, minimizing penalty functional
n p+n
J(x,M) = //^J + dt. (10)
q=0 q=0 U
Here ju,v> 0, ¡u + v = 1, Nq&L(U), q = 0, 1, ..., p + n, are self-adjoint positively defined operators, and x(i) is the target state of the system.
Theorem 2. If the pencil B is polynomially A -bounded, conditions (6)-(8) are satisfied, then for arbitrary x0, xTk e X, k = 0, n -1 and y e Hp+n (Y) there exists a unique optimal control of solutions to problem (3), (4).
Reducing problem (1), (2) to equation (3) put
X = (x e Wl (Q) : x(s) = 0,5 e ÔQ}, Y = W (Q),
д 2
where W2 (Q) are Sobolev spaces. Define the operators A = A-A, B2= A- X, B0 = a—-,
Sx3
B3 =Bl=0. Operators A,B3,B2,Bl,B0 <eL(X;Y) for all /e{0}UN.
Denote by <pijk = J sin ^^^^sin sin 1 the eigenfunctions of the Dirichlet problem for the
J { a b c J
Laplace operator, where i, j, k e N, and the corresponding eigenvalues are denoted by
ЯЧ'к \
ni ]2+(5i Y+{жк 42
a
b
c
Since (q)ijk} с Cx (Q), then we get
/j4A - jU3B3 - fj2B2 - juBr - B0 = ^
i,j,k=1
(Л - j )/4 + (Л'- Äjk )/2 - «
жк c
< (Рук> (Рук,
where < •,• > is an inner product in
Lemma 1. [15] Let one of the following conditions be satisfied:
(i) ¿éo-(A) ;
(ii) (Aecr(
Then the pencil B is polynomially (A, 0) -bounded, and conditions (6) are fulfilled. The A -spectrum of the pencil B consists of solutions /.¿k ./ = 1.4 of equation
(j -Л)/4 + (j -Я ' )// -«ж ] =0.
(11)
Construct the projector
CO
P =
I, if (i) holds,
1 - Z {'^'Jk )v,jk,if (ii) holds.
To construct the projectorl'fin choose a domain rocC, containing, for example, a finite set cr0' (B) of points /.¿i- of the A - spectrum of the pencil B and such that oT0 f| cr()! (B) = 0. As it is easy
to see, the domain r0 can be chosen such that cT0 = /0 is a contour. Thus, condition (7) is satisfied. Consider the initial-final problem
Z <xt(-^)-xi,(Pyk >(Pyk = 0, Z <Xtt(-,0)-X2,<Pyk>(Pyk=0, Z _ <xm(-,°)-xt,<Pyk><Pyk=0,
Z <xM-xL<p,jk><p,jk = ^
Z <xt(-s)-xl,(Pyk >(P,jk = 0, Z < Xtt (•» r) - x2 , (Pyk > (Pyk = 0, Z _ <Xttt(-^)-Xl'<Pyk><Pyk=^
hk^A^m ]
for equation (1) with boundary conditions (2).
Theorem 3. For any a e R and X e R such that the conditions of Lemma 1 are fulfilled, and for
any re R+, x°, xk e X, k = 0,3, there exists a unique solution to the optimal control problem for equation (1) with conditions (2), (12) that minimizes functional (10).
References
1. Sveshnikov A.G., Al'shin A.B., Korpusov M.O., Pletner Yu.D. Lineynye i nelineynye uravneniya sobolevskogo tipa (Linear and Nonlinear the Sobolev Type Equations). Moscow, FIZMATLIT Publ., 2007, 736 p. (in Russ.).
2. Zamyshlyaeva A.A. The Higher-Order Sobolev-Type Models. Bulletin of the South Ural State University, Series "MathematicalModelling, Programming & Computer Software", 2014, Vol. 7, no. 2, pp. 5-28. (in Russ.). DOI: 10.14529/mmp140201
3. Sagadeeva M.A., Zagrebina S.A., Manakova N.A. Optimal Control of Solutions of a Multipoint Initial-Finish Problem for Non-Autonomous Evolutionary Sobolev Type Equation. Evolution equations and control theory, 2019, Vol. 8, no. 3, pp. 473-488. DOI: 10.3934/eect.2019023
4. Manakova N.A., Bogatyreva E.A. Mathematical Model of the Start Control of Electric Field Potential in Conducting Medium without Dispersion Considering Relaxation. 2016 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), Chelyabinsk, 2016, pp. 1-5. DOI: 10.1109/ICIEAM.2016.7911711
5. Zamyshlyaeva A.A., Tsyplenkova O.N. Optimal Control of Solutions to the Showalter-Sidorov Problem in a Model of Linear Waves in Plazma. Journal of Computational and Engineering Mathematics, 2018, Vol. 5, no. 4, pp. 46-57. DOI: 10.14529/jcem180404
Математика
6. Bogatyreva E.A., Manakova N.A. Numerical Simulation of the Process of Nonequilibrium Counterflow Capillary Imbiition. Computational mathematics and mathematical physics, 2016, Vol. 56, no. 1, pp. 132-139. DOI: 10.1134/S0965542516010085
7. Keller A.V., Ebel A.A. Parallelization of Numerical Algorithm for Optimum Dynamic Measurement Problem Solution. Proc. 2nd International Ural Conference on Measurements (UralCon), Chelyabinsk, 2017, pp. 372-377. DOI: 10.1109/URALC0N.2017.8120738
8. Zakirova G.A. Inverse Spectral Problems and Mathematical Models of Continuum Mechanics. Bulletin of the South Ural State University, Series "Mathematical Modelling, Programming & Computer Software ", 2019, Vol. 12, no. 2, pp. 5-24. DOI: 10.14529/mmp190201
9. Shafranov D.E., Kitaeva O.G. The Barenblatt-Zheltov-Kochina Model with the Showalter-Sidorov Condition and Additive "white noise" in Spaces of Differential Forms on Riemannian Manifolds without Boundary. Global and Stochastic Analysis, 2018, Vol. 5, no. 2, pp. 139-152.
10. Sviridyuk G.A., Fedorov V.E. Linear Sobolev Type Equations and Degenerate Semigroups of Operators, VSP, Utrecht-Boston-Köln-Tokyo, 2003, 216 p. DOI: 10.1515/9783110915501
11. Manakova N.A., Svirigyuk G.A. An Optimal Control of the Solutions of the Initial-Final Problem for Linear Sobolev Type Equations with Strongly Relatively p-Radial Operator. Semigroup of operators - theory and applications, 2015, pp. 213-224. DOI: 10.1007/978-3-319-12145-1_13
12. Sviridyuk G.A., Efremov A.A. Optimal Control of Sobolev Type Linear Equations with Relatively p-Sectorial Operators. Differential Equations, 1995, Vol. 31, no. 11, p. 1882-1890.
13. Zamyshlyaeva A.A., Tsyplenkova O.N., Bychkov E.V. Optimal Control of Solutions to the Showalter-Sidorov Problem for the Sobolev type Equation of Higher Order. Proc. 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), Chelyabinsk, 2016, pp. 1-4. DOI: 10.1109/ICIEAM.2016.7911725
14. Favini A., Zagrebina S.A., Sviridyuk G.A. Multipoint Initial-Final Value Problems for Dynamical Sobolev-type Equations in the Space of Noises. Electronic Journal of Differential Equations, 2018, Vol. 2018, p. 128.
15. Zamyshlyaeva A.A. Lineynye uravneniya sobolevskogo tipa vysokogo poryadka (Linear Sobolev Type Equations of Higher Order). Chelyabinsk, Publishing center of SUSU, 2012, 88 p. (in Russ).
Received October 5, 2019
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2019, vol. 11, no. 4, pp. 26-31
УДК 517.9 DOI: 10.14529/mmph190403
ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ РЕШЕНИЯМИ
НАЧАЛЬНО-КОНЕЧНОЙ ЗАДАЧИ
ДЛЯ МОДЕЛИ ЛИНЕЙНЫХ ВОЛН В ПЛАЗМЕ
А.А. Замышляева, О.Н. Цыпленкова
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: zamyshliaevaaa@susu.ru
В работе исследована задача оптимального управления для уравнения соболевского типа высокого порядка с относительно полиномиально ограниченным пучком операторов. Результаты применены к исследованию оптимального управления решениями начально-конечной задачи для модели линейных волн в плазме. Первые результаты по уравнению, которое описывает линейные ионно-звуковые волны в незамагниченной плазме, и изучению некоторых свойств этих волн были получены Ю.Д. Плетнером. Начально-конечные условия, поставленные для уравнения соболевского типа четвертого порядка, являются обобщением условий в задаче Коши, которые неразрешимы при произвольных начальных значениях. Работа основывается на методе фазового про-
странства, разработанного Г.А. Свиридюком, и теории относительно полиномиально ограниченных пучков операторов, разработанной А.А. Замышляевой. В статье рассмотрено уравнение, которое описывает ионно-звуковые волны в плазме во внешнем магнитном поле.
Ключевые слова: уравнения соболевского типа высокого порядка c относительно полиномиально ограниченным пучком операторов; модель линейных волн в плазме; оптимальное управление; начально-конечные условия.
Литература
1. Линейные и нелинейные уравнения соболевского типа / А.Г. Свешников, А.Б. Альшин, МО. Корпусов, Ю.Д. Плетнер. - М.: ФИЗМАТЛИТ, 2007. - 736 с.
2. Замышляева, А.А. Математические модели соболевского типа высокого порядка / А.А. Замышляева // Вестник ЮУрГУ. Серия: «Математическое моделирование и программирование». - 2014. - Т. 7, № 2. - С. 5-28.
3. Sagadeeva, M.A. Optimal Control of Solutions of a Multipoint Initial-Finish Problem for Non-Autonomous Evolutionary Sobolev Type Equation / M.A. Sagadeeva, S.A. Zagrebina, N.A. Manakova // Evolution Equations and Control Theory. - 2019. - Vol. 8, № 3. - P. 473-488.
4. Manakova, N.A. Mathematical Model of the Start Control of Electric Field Potential in Conducting Medium without Dispersion Considering Relaxation / N.A. Manakova, E.A. Bogatyreva // 2016 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). -Chelyabinsk, 2016. - P. 1-5.
5. Zamyshlyaeva, A.A. Optimal Control of Solutions to the Showalter-Sidorov Problem in a Model of Linear Waves in Plazma / A.A. Zamyshlyaeva, O.N. Tsyplenkova // Journal of Computational and Engineering Mathematics. - 2018. - Vol. 5, № 4. - P. 46-57.
6. Богатырева, Е.А. Численное моделирование процесса неравновесной противоточной капиллярной пропитки / Е.А. Богатырева, Н.А. Манакова // Журнал вычислительной математики и математической физики. - 2016. - Т. 56, № 1. - С. 125-132.
7. Keller, А^. Parallelization of Numerical Algorithm for Optimum Dynamic Measurement Problem Solution / A.V. Keller, A.A. Ebel // Proc. 2nd International Ural Conference on Measurements. -2017. - p. 372-377.
8. Zakirova, G.A. Inverse Spectral Problems and Mathematical Models of Continuum Mechanics / G.A. Zakirova // Вестник ЮУрГУ. Серия «Математическое моделирование и програмирование». - 2019. - Т. 12, № 2. - С. 5-24.
9. Shafranov, D.E. The Barenblatt-Zheltov-Kochina Model with the Showalter-Sidorov Condition and Additive «white noise» in Spaces of Differential Forms on Riemannian Manifolds without Boundary / D.E. Shafranov, O.G. Kitaeva // Global and Stochastic Analysis. - 2018. - Vol. 5, no. 2. -C.139-152.
10. Sviridyuk, G.A. Linear Sobolev Type Equations and Degenerate Semigroups of Operators / G.A. Sviridyuk, V.E. Fedorov. - Utrecht; Boston; Köln; Tokyo: VSP, 2003. - 216 p.
11. Manakova, NA. An Optimal Control of the Solutions of the Initial-Final Problem for Linear Sobolev Type Equations with Strongly Relatively ^-Radial Operator / NA. Manakova, GA. Sviridyuk // Semigroup of Operators - Theory and Applications. - 2015. - P. 213-224.
12. Свиридюк, Г.А. Оптимальное управление линейными уравнениями типа Соболева с относительно ^-секториальными операторами / Г.А. Свиридюк, А.А. Ефремов // Дифференциальные уравнения. - 1995. - Т. 31, № 11. - С. 1912-1919.
13. Zamyshlyaeva, A.A. Optimal Control of Solutions to the Showalter-Sidorov Problem for the Sobolev Type Equation of Higher Order / A.A. Zamyshlyaeva, O.N. Tsyplenkova, E.V. Bychkov // Proc. 2nd International Conference on Industrial Engineering, Applications and Manufacturing, ICIEAM 2016. - 2016. - P. 1-4.
14. Favini, A. Multipoint Initial-Final Value Problems for Dynamical Sobolev-type Equations in the Space of Noises / A. Favini, S.A. Zagrebina, G.A. Sviridyuk // Electronic Journal of Differential Equations. - 2018. - Vol. 2018. - P. 128.
15. Замышляева, А.А. Линейные уравнения соболевского типа высокого порядка / А.А. Замышляева. - Челябинск: Изд. центр ЮУрГУ, 2012. - 107 c.
Поступила в редакцию 5 октября 2019 г.