CC BY
2
УДК 519.6
БОТ: 10.18698/2309-3684-2021-1-5566
Численное моделирование нестационарной задачи оптимального размещения источников тепла минимальной мощности в однородной среде
© Б.Х. Хайиткулов
Национальный Университет Узбекистана, г. Ташкент, 100174, Узбекистан
Данная работа посвящена численному решению нестационарной задачи оптимального размещения источников тепла минимальной мощности. Постановка задачи требует одновременного выполнения двух условий. Первое условие — обеспечить нахождение температуры в пределе минимальных и максимальных температур за счет оптимального размещения источников тепла с минимальной мощностью в прямоугольнике. Второе условие заключается в том, чтобы суммарная мощность источников тепла, используемых для обогрева, была минимальной. Эта задача изучалась в стационарных условиях в работах других учёных. Однако в нестационарном случае задача не рассматривалось. Поскольку найти непрерывное решение краевой задачи сложно, то ищем численное решение задачи. Трудно найти интегральный оператор с непрерывным ядром (функция Грина). Найдено численное значение функции Грина в виде матрицы. Предложен новый алгоритм численного решения ностационарной задачи оптимального управления размещением источников тепла с минимальной мощностью в процессах, описываемых дифференциальными уравнениями с частными производными параболического типа. Предложена новая методика численного решения. Построена математическая и численная модель процессов, описываемых уравнением теплопроводности с постоянными коэффициентами, заданными для первой краевой задачи. Краевая задача изучается для двумерного случая. Для численного решения задачи использовалась неявная конечно-разностная схема. По этой схеме была создана система разностных уравнений. Сформированная система разностных уравнений приведена к задаче линейного программирования. Задача линейного программирования решается с помощью М-метода. При каждом значении времени решается задача линейного программирования. Предложен новый подход к численному решению задач. Приведена общая блок-схема алгоритма решения нестационарной задачи оптимального управления размещением источников тепла с минимальной мощностью. Разработан алгоритм и программное обеспечение для численного решения задачи. Приведено краткое описание программного обеспечения. На конкретных примерах показано, что численное решение краевой задачи находится в заданных пределах, сумма оптимально размещенных источников тепла с минимальной мощностью дает минимум функционалу. Визуализированы результаты вычислительного эксперимента.
Ключевые слова: уравнение теплопроводности, нестационарные задачи, оптимальное размещение, источник тепла, неявные схемы, М-метод
Введение. Одним из самых распространенных объектов в различных сферах деятельности человека является система источников тепла, тепловой баланс в отапливаемых помещениях.
Математическое моделирование таких систем ставит задачу об оптимальном размещении источников тепла в отапливаемых помещениях, что связано с ресурсосберегающими инженерными технологиями. Задача оптимального размещения источников тепла отапливаемых площадей всегда была актуальной при проектных работах в строительстве, теплицах и других технических и технологических сферах. Задачи этого типа были изучены для стационарного случая. Однако задача оптимального размещения источников тепла с минимальной мощностью для нестационарного случая не изучена. Трудно найти точное решение этих задач. Поэтому в статье показано численное решение задачи.
В работах [1-4] изучалась постановка задачи оптимального управления краевой задачей параболического типа, существование решения и решение приближенными методами. В работе [5] рассматривается задача оптимизации плотности источников тепла в стационарных процессах, описываемых эллиптическими уравнениями, задаваемыми третьим граничным условием. В работе
[6] разработаны новые методы решения задачи оптимального размещения источников в неоднородных средах, описываемых уравнениями эллиптического и параболического типа. Показан метод выбора оптимального варианта размещения источника из допустимого множества на основе метода ветвей и границ. В работе
[7] изучено влияние эффекта Коанда на скорость теплопередачи в одиночном цилиндре с учетом расстояния до стенки. Найдено оптимальное расстояние цилиндра от стенки. В работе [8] влияние управляющих параметров на тепловое явление в случае смешанного конвективного теплообмена в закрытом помещении с дискретными источниками тепла исследовано численно методом коллокации кубических сплайнов. В работе [9] рассмотрены задачи оптимального нагрева помещения на основе принципа максимума Понтрягина. В работе [10] рассмотрена задача энергоэффективного теплоснабжения здания в системе центрального отопления. В работах [11-13] предложен метод численного решения нестационарной задачи оптимального размещения источников тепла с минимальной мощностью в процессах, описываемых уравнениями параболического типа. Разработан алгоритм и комплекс программ для численного решения нестационарных задач оптимального управления расположением источников тепла и визуализации полученных результатов. В работе [14] рассмотрена краевая задача параболического типа. Распределение тепла в рассматриваемом теле контролируется функцией, которая находится на границе тела. В работе [15] изучена задача оптимального управления процессами, описываемыми уравнением теплопроводности. Управляющий параметр задан в граничном условии и достиг минимума
функционала, задаваемого интегральным квадратичным выражением. Показан метод нахождения допустимого управления, дающего минимум функционалу.
В настоящей работе рассмотрена задача управления теплопроводности на основе оптимизации линейного целевого функционала с учетом ограничений, которая решается на основе аппроксимации и сведения к задаче линейного программирования. В работе предлагаются метод и алгоритм решения нестационарной задачи поддержания температуры внутри области в заданных пределах, путем оптимального размещения источников тепла в прямоугольнике. Для проведения вычислительных экспериментов было разработано программное обеспечение.
Математическая постановка задачи. В области Б = {а < х < Ь, с < у < й, 0 < г < Т} требуется найти функцию /(х, у, г) > 0 такую, что для любого г е [0, Т] линейный функционал
Ь й
} = | \/(х, у, г)йуйх ^ Ш1П, (1)
достигал минимума и удовлетворялись следующие условия:
ды ¥ = %
{ д 2ы д 2ыЛ
кдх2 ду2 ,
+ /(х, у, г), х е (а, Ь), у е (с, й), г е (0, Т],
ы(х, у, 0) = ы (х, у), а < х < Ь, с < у < й, (2)
ы(а, у, г) = м(у, г), ы(Ь, у, г) = м(у, г), с < у < й, 0 < г < Т, ы(х,с,г) = м(х,г), ы(х,й,г) = м(х, г), а < х < Ь, 0 < г < Т,
т(х, у, г) < ы(х, у, г) < М(х, у, г), (х, у, г) е Б, (3)
где ы = ы (х, у, г) — температура в точке (х, у) прямоугольника в момент г; % — коэффициент температуропроводности;
ы0(х, у), м (у, г), Му, г), Мз(х, г), х, г), т(х, у, г), М(х, у, г)
— заданные непрерывные функции, удовлетворяющие условия сопряжения
м(с, г) = Мз(а, г), м(й, г) = м(а, г), М (с, г) = м (Ь, г), м(й, г) = м (Ь, г).
Функции т( х, у, г), М (х, у, г) имеют смысл функций минимального и максимального профиля температуры в области Б соответственно. Источник тепла описывается квадратично интегрируемой функцией /(х,у,г) в пространстве Ь2(П).
Пусть
ды
Ьы =--у
дг 4
д ы д ы
2, Л
дх2 ду2)
Оператор Ь, определенный в Ь2 (Б), имеет обратный Ь1. Здесь
Ь"1 интегральный оператор с непрерывным ядром (функция Грина). Используя его, можно записать задачу (1)-(3) в следующем виде:
/(•, •, •) е Ь2 (Б), /(х, у, г) > 0,
т(х, у, г) < (ь)(х, у, г) < м(х, у, г).
(4)
Конечно-разностная аппроксимация задачи. Поскольку трудно найти непрерывное решение задачи (1)-(4), ищем численное решение задачи. В этом случае, используя неявную схему, заменяем задачу (2) конечно-разностным уравнением.
Введем в Б равномерную по трем переменным разностную сетку
Къ = ёь Х®Т = {(хг = у, = А, гк = ктХ
I = 0,N1, } = 0,N2, к = 0,N3}
с шагами
Ъ = (ь - а)/N, Ъ = (а - с)/N, г = т / N.
Неявная разностная схема для задачи (2) имеет вид:
ы;1 - ы;
с ык, +1 - 2ык+1 + ык+1 ыкг! - 2ык+1 + ык+1 >
= х
1+1}
1 -1}.+- }+1
*у -1
к2
+ Г. +1 ¿11 '
I = 1,2,...,N -1, 1 = 1,2,...,N -1, к = 0,1,...,N -1,
ы0 = ы0(х,у}), I = 0,1,.,N1, } = 0,1,.,N2,
0} ПУ-Т р1к+1Ь "Щ} ык0+1 ,г.+1), ы™ =^1,г.+1), I = 0Ж, к = 0,N3 -1
Введем обозначения
Х7 =
Г к12 к2 )
X = -4, 7 = -4. к12' й?
(5)
ы0к11 = (у 1,гк+1), ы£1 = ^2(у,,гк+1), } = 0,к = 0,N3 -1,
Рассмотрим расширенную матрицу системы:
г
<
A =
ХУ У 0 ... О X 0 ...... О
У ХУ У О ... О X О ... О
О ... О X О ... О У ХУ У ^ О ...... О X О ... О У ХУ
Получим
G = Л"1.
Аппроксимируем задачу (1)-(5) в виде задачи линейного программирования. Разделим область Б по х, у, г соответственно на N,
N, N равных частей:
»3 1»! 1»2
д=ииич
к=1 ¿=1 у=1
где
Бк = {(х,у,г), хм < х < х,, у,_1 < у < у,, гк_ < г < гк}
, Ч _1 - ь - Ч)
I=1,N1, ]=1,N2, к=1,N3.
В пространстве Ь2(П) функции
/, = / (х,, у,, гк), (^ ^г) е Б
(I = 1, N _ 1, , = 1, N _ 1, к = 1, N)
определяются как кусочно-постоянные функции. Отсюда получим
N5 N _1 щ _1
/ (х, у, г) «2: е Е /!.
к=1 г=1 ,=1
Пусть
*„=С,т!;=т(х1,уМ, м;=М(х„уМ, =/Ц,
N = (N1 _ 1)(N2 _ 1), р = д, д = (I _ 1)(#2 _ 1) + ], р = 1, N д = 1, N
1 = 1,N1 _ 1, ] = 1,N2 _1, к = 1,N3.
Подставим выражение / (х, у, г) в (1) и заменим неравенство (4) на сеточные функции.
После этого получим следующую задачу линейного программирования:
N -1 N2-1
Л{f} = Z ZCmesDf)fj ^ min, k = 1,2,...,N3,
i=i j=i
N _ _______
9=1
/;>0, </ = 1,2,к = 1,2,...,Ы3.
Задача (6) решается M-методом [16]. Численное решение задачи (2) находится с помощью
N
< = 2> J?, 0" = 0> ■-1) ■- W -1) +1, J = ■-1) mod (N2 -1) +1),
9=1
где обозначение -5- есть символ целочисленного деления, a mod есть
символ остатка от деления. Найденная ц является функцией,
дающей минимум функционалу (1).
Численные эксперименты. Для приближенного решения задачи (1)-(6) разработано программное обеспечение на языке C#. Для визуализации результатов использовались 3D-графики MathCAD.
Была расмотрена следующая задача. Требуется найти оптимальное расположение источников тепла с минимальной мощностью в прямоугольнике. Задача решалась при следующих входных данных: х, y е [0,1], коэффициент температуропроводности
Х = 1 м2 /с, начальное и граничные условия определяются функциями:
U (х, y) = 2 + х 2+ y2 м/с, д (y, t) = 2 + y2 + t2 м/с, Д(у,t) = 3 + y2 +t2 м/с, д(х,t) = 2 + х2 +t2 м/с, д (х, t) = 3 + х2 +t2 м/с, минимальная и максимальная температуры задаются функциями
т(х, y, t) = 1 + х 2+ y2 +t2 К, M (х, y, t) = 4 + х2 + y2 +t2 K, окончание времени T = 1. Расчетная сетка с числом источников (N -1) х (N -1) X N = 9 X 9 х10.
Минимальное значение функционала при численном решении равно Jmin = 92.36 K • м/с.
60
На рис. 1 представлены результаты численного решения задачи (6). Представлены результаты с минимальным (границы с синим цветом, ниже), максимальным (границы с красным цветом, выше) и приближенным (зеленым цветом, посередине) значением температуры. На рис. 2 показано оптимальное расположение источников тепла с минимальной мощностью в виде гистограммы.
т,и,М |_7
М (х, у, г)
6
ы (х, у, г)
~и\\\\Л\\////7/7У/
10 ^т----т(х,у,1) ^-----
Рис. 1. График решения задачи (6) при г = Т
/ (х, у, г)
600 400 200
0
Рис. 2. Оптимальное расположение источников тепла / (х, у, г)
4
т, ы, М
М (х, у, г)
М (х, у, г)
т, ы, М 6 5
8 6 х. - , ч
10.. .8-- " т (х,у,г)
10 8
~6 .8 ш тх .
0 "----
б)
Рис. 3. График решения задачи (6): а — при г = 0.1; б — при г = 0.5
6 8 ю
Из рисунков 1 и 3 видно, что решение задачи лежит в заданном пределе, т.е. решение удовлетворяет неравенству (3). Видно, что значение ы (х, у, г) практически равно минимальной температуре. Это
означает, что функционал {/} достигает минимума.
а) б)
Рис. 4. Оптимальное расположение источников тепла: а — при г = 0.1; б — при г = 0.5
На рис. 2 и 4 мощность оптимально размещенных источников тепла отображается в виде гистограммы. Все источники мощности выделены зеленым (рис. 2) и разными цветами (рис. 4). В этом случае источники большой мощности расположены в основном на границе внутренних узловых точек. Черные линии на рисунке образованы наложением источников тепла в каждый фиксированный момент времени.
Выводы. Предложены методика и алгоритм решения нестационарной задачи поддержания температуры внутри области в
заданных пределах путем оптимального размещения источников тепла в прямоугольнике. Задача решена на основе численного моделирования процесса распространения тепла и последовательного решения задач линейного программирования. Результаты вычислительного эксперимента подтверждают достижение функционалом минимума и решение основной задачи. Реализация комплекса программ позволит анализировать оптимальность систем отопления на строительных площадках, визуализировать размещение источников тепла разной мощности, сократить время расчетов, широко внедрить энергосберегающие технологии, резко снизить энергозатраты и повысить энергоэффективность. Это способствует повышению точности и эффективности проектных работ.
Работа выполнена при финансовой поддержке Узбекского фонда фундаментальных исследований (проект ОТ-Ф4-33).
ЛИТЕРАТУРА
[1] Бутковский А.Г., Малый С.А., Андреев Ю.Н. Управление нагревом металла. Москва, Металлургия, 1981, 272 с.
[2] Лионс Ж.Л. Оптимальное управление системами, описываемыми уравнениями с частными производными. Москва, Мир, 1972, 416 с.
[3] Егоров А.И. Оптимальное управление тепловыми и диффузионными процессами. Москва, Наука, 1978, 464 с.
[4] Федоренко Р.П. Приближённое решение задач оптимального управления. Москва, Наука, 1978, 488 с.
[5] Осипов О.В., Брусенцев А.Г. Оптимальное расположение источников тепла внутри областей сложной геометрической формы. Математическое моделирование, 2019, т. 31, № 4, с. 3-16.
[6] Ахметзянов А.В, Кулибанов В.Н. Задача оптимального выбора координат доразбуривания добывающих скважин на нефтяных месторождениях. Автоматика и телемеханика, 2002, № 11, с. 3-12.
[7] Kapjor A., Durcansky P., Vantuch M. Effect of heat source placement on natural convection from cylindrical surfaces. Energies, 2020, vol. 13, iss. 17, p. 1-13.
[8] Hsu T.H., W Ang S.G. Mixed convection in a rectangular enclosure with discrete heat sources. Numerical Heat Transfer, Part A: Applications, 2010, vol. 38, iss. 6, p. 627-652.
[9] Мирская С.Ю, Сидельников В.И. Экономичный обогрев помещения как задача оптимального управления. Технико-технологические проблемы сервиса, 2014, № 4(30), с. 75-78.
[10] Сабденов К.О, Байтасов Т.М. Оптимальное (энергоэффективное) теплоснабжение здания в системе центрального отопления. Известия Томского политехнического университета. Инжиниринг георесурсов, 2015, т. 326, № 8. с. 53-60.
[11] Khaitkulov B.Kh. Homogeneous different schemes of the problem for optimum selection of the location of heat sources in a rectangular body. Solid State Technology, 2020, vol. 63, iss. 4, pp. 583-592.
[12] Хайиткулов Б.Х. Консервативные разностные схемы по оптимальному
выбору местоположения источников тепла в стержне. Математическое моделирование и численные методы, 2020, № 3, с. 85-98.
[13] Хайиткулов Б.Х. Конечно-разностный метод решения нестационарных задач управления конвекцией-диффузией. Вестник Томского государственного университета. Управление, вычислительная техника и информатика, 2021, № 57, с. 45-52.
[14] Тухтасинов М.Т., Абдуолимова Г.М., Хайиткулов Б.Х. Граничное управление распространением тепла в ограниченном теле. Бюллетень Института математики, 2019, № 5, с. 1-10.
[15] Егоров А.И., Знаменская Л.Н. Об управлении процессом теплопроводностис квадратичным функционалом качества. Журнал вычислительной математики и математической физики, 2017, т. 57, № 12, с. 2053-2064.
[16] Dantzig G.B. Linear programming and extensions. Princeton, Princeton University Press, 2016, 656 p.
Статья поступила в редакцию 31.07.2023
Ссылку на эту статью просим оформлять следующим образом:
Хайиткулов Б. Х. Численное моделирование нестационарной задачи оптимального размещения источников тепла минимальной мощности в однородной среде. Математическое моделирование и численные методы, 2024, № 1, с. 55-66. Хайиткулов Бобомурод Хайитович — канд. физ.-мат. наук, доцент кафедры «Алгоритмы и технологии программирования» Национального Университета Узбекистана. e-mail: [email protected].
Numerical modeling of the non-stationary problem of optimal placement of heat sources of minimum power in a homogeneous medium
© B.Kh. Khayitkulov National University of Uzbekistan, Tashkent, 100174, Uzbekistan
This work is devoted to the numerical solution of the non-stationary problem of optimal placement of heat sources of minimum power. The statement of the problem requires the simultaneous fulfillment of two conditions. The first condition is to ensure that the temperature is within the limits of minimum and maximum temperatures due to the optimal placement of heat sources with a minimum power in the rectangle. The second condition is that the total power of the heat sources used for heating is minimal. This problem was studied under stationary conditions in the works of other scientists. However, the problem was not considered in the non-stationary case. Since it is difficult to find a continuous solution to the boundary value problem, we are looking for a numerical solution to the problem. It is difficult to find an integral operator with a continuous kernel (Green's function). We find the numerical value of the Green's function in the form of a matrix. A new algorithm for the numerical solution of a non-stationary optimal control problem for the placement of heat sources with a minimum power in processes described by parabolic partial differential equations is proposed. A new technique for numerical solution is proposed. A mathematical and numerical model of the processes described by the heat conduction equation with constant coefficients given
for the first boundary value problem is constructed. The boundary value problem is studied for the two-dimensional case. An implicit finite difference scheme was used to solve the problem numerically. According to this scheme, a system of difference equations was created. The formed system of difference equations is reduced to a linear programming problem. The problem of linear programming is solved using the M-method. For each time value, a linear programming problem is solved. A new approach to the numerical solution of problems is proposed. A general block diagram of the algorithm for solving the non-stationary problem of optimal control of the placement of heat sources with a minimum power is given. An algorithm and software for the numerical solution of the problem have been developed. A brief description of the software is given. On specific examples, it is shown that the numerical solution of the boundary value problem is within the specified limits, the sum of optimally placed heat sources with a minimum power gives a minimum to the functional. The results of the computational experiment are visualized.
Keywords: heat equation, non-stationary problems, optimal placement, heat sources, implicit schemas, big M method
REFERENCES
[1] Butkovskii A.G., Malyi S.A., Andreev Yu.N. Upravlenie nagrevom metalla [Metal heating control]. Moscow, Metallurgiya Publ., 1981, 272 p.
[2] Lions J.-L. Controle optimal de systemes gouvernes par des equations aux derivees partielles. Paris, Dunod, 1968, 416 p.
[3] Egorov A.I. Optimal'noe upravlenie teplovymi i diffuzionnymi protsessami [Optimal control of thermal and diffusion processes]. Moscow, Nauka Publ., 1978, 464 p.
[4] Fedorenko R.P. Priblizhennoe reshenie zadach optimal'nogo upravleniya [Approximate solution of optimal control problems]. Moscow, Nauka Publ., 1978, 488 p.
[5] Osipov O.V., Brusentsev A.G. Optimal location of heat sources inside areas with complex geometric shapes. Mathematical Models and Computer Simulations, 2019, vol. 11, no. 6, pp. 905-913.
[6] Akhmetzyanov A.V., Kulibanov V.N. Optimal choice of coordinates for oil well drilling. Automation and Remote Control, 2002, vol. 63, no. 11, pp. 1699-1706.
[7] Kapjor A., Durcansky P., Vantuch M. Effect of heat source placement on natural convection from cylindrical surfaces. Energies, 2020, vol. 13, iss. 17, pp. 1-13.
[8] Hsu T.H., W Ang S.G. Mixed convection in a rectangular enclosure with discrete heat sources. Numerical Heat Transfer, Part A: Applications, 2010, vol. 38, iss. 6, pp. 627-652.
[9] Mirskay S.Yu., Sidelnikov V.I. Efficient heating of the room as the optimal control problem. Tehniko-tehnologicheskie problemy servisa [Technical and Technological Problems of Service], 2014, no. 4(30), pp. 75-78.
[10] Sabdenov K.O., Baytasov T.M. Optimal (energy efficient) heat supply to buildings in central heating system. Izvestija Tomskogo politehnicheskogo universiteta. Inzhiniring georesursov [News of Tomsk Polytechnic University. Geo-Resource Engineering], 2015, vol. 326, no. 8, pp. 53-60.
[11] Khaitkulov B.Kh. Homogeneous different schemes of the problem for optimum selection of the location of heat sources in a rectangular body. Solid State Technology, 2020, vol. 63, no. 17, pp. 583-592.
[12] Khayitkulov B.Kh. Conservative difference schemes for the optimal selection of the location of heat sources in the rod. Matematicheskoe modelirovanie i
chislennye metody [Mathematical Modeling and Computational Methods], 2020, no. 3, pp. 85-98.
[13] Khayitkulov B.Kh. Finite-difference method for solving non-stationary problems of convection-diffusion control. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika, 2021, no. 57, pp. 45-52.
[14] Tukhtasinov M.T., Abduolimova G.M., Khayitkulov B.Kh. Boundary control of heat propagation in a bounded body. Bjulleten' Instituta matematiki [Bulletin of the Institute of Mathematics], 2019, no. 1, pp. 1-10.
[15] Egorov A.I., Znamenskaya L.N. Control of a heat conduction process with a quadratic cost functional. Computational Mathematics and Mathematical Physics, 2017, vol. 57, no. 12, pp. 2005-2016.
[16] Dantzig G.B. Linear programming and extensions. Princeton, Princeton University Press, 2016, 656 p.
Khayitkulov B.Kh., Cand. Sc. (Phys. — Math.), Associate Professor, Department of Algorithms and Programming Technologies of the National University of Uzbekistan. email: [email protected]
