Метод моделирования водогрейных котельных малой мощности и их тепловой нагрузки

Обложка


Цитировать

Полный текст

Аннотация

Представлен новый метод моделирования автономных водогрейных котельных малой мощности и их тепловой нагрузки. Предложенный подход основан на представлении основных элементов тепловой схемы (котел, теплообменник, потребители отопления и ГВС) в виде группы взаимодействующих тепловых масс и использовании системы дифференциальных уравнений для описания динамики теплообмена. Для определения коэффициентов теплопередачи между тепловыми массами применяется анализ установившегося режима с учетом граничных условий на внешние потоки энергии. Компьютерная реализация модели выполнена в среде MATLAB Simulink, параметры идентифицированы по реальным данным котельной станции с двумя котлами по 1 МВт и КПД 90 %. В вычислительном эксперименте проведены статичные и динамические тесты: прогрев котла при номинальной нагрузке и его остывание при отключенном теплоносителе. Результаты моделирования продемонстрировали высокую степень согласования с эксплуатационными данными (средняя абсолютная ошибка температур составила менее 0,5 °C), что подтверждает адекватность и надежность предложенного метода. Универсальность подхода позволяет масштабировать модель на различное оборудование и режимы работы без существенного роста вычислительной нагрузки. Предложенный метод может служить эффективным инструментом инженерного анализа, проектирования и последующей оптимизации систем теплоснабжения, обеспечивая сокращение затрат времени и ресурсов на опытную проверку.

Полный текст

Введение

Котельные широко используются для отопления и подачи горячей воды потребителям. В качестве топлива используются природный газ, уголь, мазут и другие энергоносители. В данном исследовании рассматривается моделирование работы автономной водогрейной котельной с нагрузкой в виде потребителей тепла и горячего водоснабжения. Автономные котельные обычно включают в себя один или несколько котлов, теплообменники, насосы, трубопроводную арматуру и систему управления, которая выполняет функции регулирования температуры воды, подаваемой потребителю, управления котлами и аварийной защиты оборудования. Современные водогрейные котельные все чаще автоматизируются для повышения их надежности и энергоэффективности [1, 2]. При проектировании и оптимизации работы котельных применяются компьютерные и математические модели, которые используются для решения задач повышения их экономичности и эффективности [3, 4], а также снижения воздействия на окружающую среду [5].

Описание объекта и постановка задачи

Основной функцией котельной является нагрев воды и ее распределение по двум основным контурам: системе отопления и системе горячего водоснабжения (ГВС). Ее упрощенная схема представлена на рис. 1. Нагретая в группе котлов вода поступает к потребителям. В процессе прохождения по системе вода отдает тепло радиаторам в жилых помещениях или теплообменникам для подогрева воды в системе горячего водоснабжения. Остывшая вода, потерявшая часть своей тепловой энергии, возвращается по обратным трубопроводам в котельную для повторного нагрева.

Основной принцип функционирования котельной - это поддержание постоянной циркуляции теплоносителя по замкнутому контуру. Насосы обеспечивают движение воды по системе, создавая необходимое давление для того, чтобы теплоноситель достигал каждого потребителя. Для управления температурой используются трехходовые клапаны, которые позволяют смешивать потоки теплоносителя разной температуры для поддержания заданной температуры в системе.

 

Рис. 1. Схема котельной: Н1, Н2, Н3 - группы насосов; К1, К2 - трехходовые клапаны

 

В данном исследовании разрабатывается математическая модель котельной, основанная на системе дифференциальных уравнений и предназначенная для анализа ее работы и последующей оптимизации. Такая модель позволяет с необходимой для задач автоматизации точностью описывать процессы, происходящие в различных узлах системы, и рассчитывать их поведение при изменении параметров.

Исследование тепловой схемы (см. рис. 1) как единой системы позволяет решить целый ряд вопросов:

  • оценка правильности выбора оборудования котельной;
  • определение динамических характеристик процессов тепломассопереноса и расчет оптимальных регуляторов систем автоматического управления;
  • исследование различных режимов работы оборудования котельной и их влияния на динамические характеристики системы;
  • сокращение затрат времени и ресурсов за счет выполнения расчетов вместо проведения экспериментальных испытаний.

Обзор методов моделирования

Общий подход к моделированию тепловых процессов предполагает представление процессов как системы трех уравнений: уравнения закона сохранения вещества, уравнения закона сохранения количества движения, уравнения закона сохранения энергии [6].

В случае моделирования котельных законы сохранения вещества и сохранения количества движения определяют гидравлическую модель потоков теплоносителя, закон сохранения энергии лежит в основе тепловых моделей. Динамическая модель гидравлической системы в данной работе не рассматривается, поскольку переходные процессы в ней затухают очень быстро и могут быть приняты безынерционными по сравнению с тепловыми. К тепловым объектам относятся котел, теплообменник, потребитель отопления и потребитель горячего водоснабжения (ГВС).

Методы моделирования тепловых объектов можно разделить на три группы: аналитические, экспериментальные и численные [7-9]. Аналитические методы основаны на построении математических моделей объекта с использованием фундаментальных законов физики [10]. К данным методам относятся модели с распределенными параметрами [11] и модели с сосредоточенными параметрами [12]. При этом для сосредоточенных моделей может применяться метод термоэлектрической аналогии [13].

Экспериментальные методы базируются на основе реальных измерений, полученных при эксплуатации объектов, которые используются для создания моделей. Статические модели могут быть построены при помощи регрессии и корреляции с установлением эмпирических зависимостей между входными и выходными величинами; динамические модели - при помощи процедур идентификации, когда параметры модели определяются посредством специально организованных экспериментов.

Численные методы применяются для решения нелинейных дифференциальных уравнений, описывающих тепловые и гидравлические процессы, с учетом сложной геометрии и временных характеристик объекта [14].

Для достижения высокой точности моделирования больше всего подходят численные методы. Высокая точность численных моделей чаще всего достигается за счет ресурсоемких вычислений, снижение вычислительных затрат достигается за счет внесения допущений в модель с учетом ограничений по желаемой точности. Аналитический метод, предполагающий получение замкнутых математических уравнений, не может быть применен для полного описания котельной станции ввиду ее высокой структурной и функциональной сложности, а экспериментальные методы дают высокую точность только в пределах наблюдаемой опытами области и сильно зависят от исходных данных.

В целом исследования по моделированию тепловых систем можно разделить на аналитические и эмпирические. Аналитические имеют высокую сложность и требуют серьезных исследований для применимости в каждом новом случае, эмпирические же более просты и широко применимы, но зависят от наличия и качества данных об объекте моделирования.

В современных исследованиях активно обсуждаются задачи моделирования процессов теплопередачи, производится поиск эффективных методов. Статья [15] посвящена разработке низкоразмерной динамической модели компактного пластинчатого теплообменника. Авторы используют теорию подобия Нуссельта для описания процессов конвекции и теплопроводности в каналах теплообменника, а также метод конечных объемов для пространственной дискретизации. Применяя метод, основанный на аппроксимированном стационарном решении исходных уравнений, авторы получили модель, которая требует меньшей вычислительной мощности и при этом точно представляет динамические процессы и нелинейности. В [16] для моделирования динамики теплообменных систем используется метод термоэлектрической аналогии. Подход, предложенный в [16], эффективен для описания динамики отдельных теплообменников или их компактных объединенных систем, где можно принять однородность температурного поля и линейность динамических процессов. Статьи [15, 16] посвящены моделированию конкретных тепловых процессов и не учитывают взаимодействие моделируемого объекта с окружающей средой.

Работа [17] посвящена построению агрегированных математических моделей прямоточного парового котла с помощью аппарата дифференциально-алгебраических уравнений в частных производных. В работе построены математические модели отдельных подсистем прямоточного парового котла, реализованные с использованием как сосредоточенных, так и распределенных параметров, а также алгебраические уравнения, отражающие гидравлические и энергетические балансы. Автор проводит анализ разрешимости дифференциально-алгебраических уравнений в частных производных, формулирует и доказывает теоремы существования и единственности решений. Предложенный подход обладает высокой математической сложностью, применение подхода требует тонкой настройки параметров и значительных вычислительных ресурсов. Подход может использоваться для описания агрегатов котельной, однако имеет все недостатки аналитических методов.

Представленный в статье [18] метод матричной формализации основывается на составлении системы уравнений баланса массы и энергии, решаемой методами математического программирования. Модель [18] описывает почти все подсистемы ТЭС, за исключением турбоустановки, и системы регенеративного подогрева питательной воды, что позволяет получить единое агрегированное математическое описание станции. Метод матричной формализации позволяет объединить расчеты различных подсистем станции в одну систему уравнений. В сравнении с распределенными моделями агрегированное матричное представление снижает вычислительную сложность и упрощает алгоритм расчета. Но объединение множества подсистем в единую модель приводит к утрате подробной динамики отдельных устройств, что снижает точность описания локальных процессов и общей модели. Для корректной работы модели требуется значительное количество экспериментально определяемых параметров, что затрудняет построение модели. Метод матричной формализации относится к эмпирическим и подходит для упрощенного описания всей системы при наличии достаточного количества реальных данных с исследуемого объекта.

Математическая модель

Для построения модели котельной в данной работе будет использоваться разработанный метод, который относится к группе численных методов и основывается на законах сохранения энергии и сохранения масс. В основе данного метода лежит представление тепловых объектов в виде группы взаимодействующих тепловых масс.

Понятийно подход можно описать следующим образом. Тепловая масса - это вещество или материал, которые входят в состав теплового объекта и являются необходимым элементом для моделирования этого объекта. Количество элементов тепловой массы должно быть достаточным для моделирования всех ключевых процессов теплопередачи, которые описывают назначение теплового объекта, а также его взаимодействие с окружающей средой. Взаимодействие элементов тепловых масс происходит через передачу энергии от элемента с большей температурой к элементу с меньшей температурой; один элемент может обмениваться энергией с несколькими соседними элементами. Также тепловая масса может обмениваться энергией с окружающей средой, которая может как вносить, так и забирать энергию.

Каждая тепловая масса имеет свою массу и теплоемкость, а также главный параметр моделирования - температуру. Из уравнения теплового баланса [19] температура i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaaaa@36D3@ -той тепловой массы определяется как

d dt T i (t)= 1 m i c i j=1 N q ij (t)+ Q i (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGKbaabaGaamizaiaadshaaaGaamivamaaBaaaleaacaWGPbaabeaa kiaacIcacaWG0bGaaiykaiabg2da9maalaaabaGaaGymaaqaaiaad2 gadaWgaaWcbaGaamyAaaqabaGccaWGJbWaaSbaaSqaaiaadMgaaeqa aaaakmaabmaabaWaaabCaeaacaWGXbWaaSbaaSqaaiaadMgacaWGQb aabeaakiaacIcacaWG0bGaaiykaiabgUcaRiaadgfadaWgaaWcbaGa amyAaaqabaGccaGGOaGaamiDaiaacMcaaSqaaiaadQgacqGH9aqpca aIXaaabaGaamOtaaqdcqGHris5aaGccaGLOaGaayzkaaaaaa@54E3@ , i= 1,N ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabg2 da9maanaaabaGaaGymaiaacYcacaWGobaaaaaa@3A28@ , (1)

где m i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaacaWGPbaabeaaaaa@37F1@  - тепловая масса; c i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4yamaaBa aaleaacaWGPbaabeaaaaa@37E7@  - ее удельная теплоемкость, N MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOtaaaa@36B8@  - количество тепловых масс объекта; q ij (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyCamaaBa aaleaacaWGPbGaamOAaaqabaGccaGGOaGaamiDaiaacMcaaaa@3B40@  - тепловой поток, приходящий от j MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOAaaaa@36D4@ -й тепловой массы; Q i (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGPbaabeaakiaacIcacaWG0bGaaiykaaaa@3A31@ - количество тепла, привнесенное или убывшее от внешней среды.

Уравнение (1) можно записать в другой форме, используя алгебраизированный оператор дифференцирования p= d dt MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiCaiabg2 da9maalaaabaGaamizaaqaaiaadsgacaWG0baaaaaa@3ABB@  [20], а также выразив T i (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaWGPbaabeaakiaacIcacaWG0bGaaiykaaaa@3A34@ . Тогда уравнение примет вид

T i (t)= 1 m i c i 1 p j=1 N q ij (t)+ Q i (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaWGPbaabeaakiaacIcacaWG0bGaaiykaiabg2da9maalaaa baGaaGymaaqaaiaad2gadaWgaaWcbaGaamyAaaqabaGccaWGJbWaaS baaSqaaiaadMgaaeqaaaaakmaalaaabaGaaGymaaqaaiaadchaaaWa aeWaaeaadaaeWbqaaiaadghadaWgaaWcbaGaamyAaiaadQgaaeqaaO GaaiikaiaadshacaGGPaGaey4kaSIaamyuamaaBaaaleaacaWGPbaa beaakiaacIcacaWG0bGaaiykaaWcbaGaamOAaiabg2da9iaaigdaae aacaWGobaaniabggHiLdaakiaawIcacaGLPaaaaaa@53C8@ , i= 1,N ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabg2 da9maanaaabaGaaGymaiaacYcacaWGobaaaaaa@3A28@ . (2)

Тепловые потоки между тепловыми массами можно описать как

q ij (t)= β ij (t) T i (t) T j (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyCamaaBa aaleaacaWGPbGaamOAaaqabaGccaGGOaGaamiDaiaacMcacqGH9aqp cqGHsislcqaHYoGydaWgaaWcbaGaamyAaiaadQgaaeqaaOGaaiikai aadshacaGGPaWaaeWaaeaacaWGubWaaSbaaSqaaiaadMgaaeqaaOGa aiikaiaadshacaGGPaGaeyOeI0IaamivamaaBaaaleaacaWGQbaabe aakiaacIcacaWG0bGaaiykaaGaayjkaiaawMcaaaaa@4E4E@ , i,j= 1,N ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiaacY cacaWGQbGaeyypa0Zaa0aaaeaacaaIXaGaaiilaiaad6eaaaaaaa@3BC7@ , (3)

где β ij (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqOSdi2aaS baaSqaaiaadMgacaWGQbaabeaakiaacIcacaWG0bGaaiykaaaa@3BEB@  - приведенный коэффициент теплообмена между элементами. Для случая, когда коэффициенты теплообмена не зависят от времени, можно определить матрицу коэффициентов

B= β ij N×N MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaCOqaiabg2 da9maadmaabaGaeqOSdi2aaSbaaSqaaiaadMgacaWGQbaabeaaaOGa ay5waiaaw2faamaaBaaaleaacaWGobGaey41aqRaamOtaaqabaaaaa@4145@ . (4)

Систему уравнений (2), (3) можно использовать для построения динамической модели процесса теплообмена. Для этого необходимо определить значения коэффициентов матрицы B. Эту задачу можно решить исходя из информации о состоянии объекта в установившемся режиме, который можно определить через типовые значения температур тепловых масс T i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaDa aaleaacaWGPbaabaGaey4fIOcaaaaa@38C8@ . В матричном виде можно записать

T= T i N×1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaCivaiabg2 da9maadmaabaGaamivamaaDaaaleaacaWGPbaabaGaey4fIOcaaaGc caGLBbGaayzxaaWaaSbaaSqaaiaad6eacqGHxdaTcaaIXaaabeaaaa a@4078@ , (5)

где T i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaDa aaleaacaWGPbaabaGaey4fIOcaaaaa@38C8@  - типичная температура i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaaaa@36D3@ -той тепловой массы. Аналогично (5) можно представить поток Q i (t) MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGPbaabeaakiaacIcacaWG0bGaaiykaaaa@3A31@ для каждой тепловой массы в статике:

Q= Q i * N×1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaCyuaiabg2 da9maadmaabaGaamyuamaaDaaaleaacaWGPbaabaGaaiOkaaaaaOGa ay5waiaaw2faamaaBaaaleaacaWGobGaey41aqRaaGymaaqabaaaaa@4031@ , (6)

где Q i * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaDa aaleaacaWGPbaabaGaaiOkaaaaaaa@3884@  - статичный приток или отток тепла i MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaaaa@36D3@ -той тепловой массы.

В установившемся режиме dT dt 0 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGKbGaamivaaqaaiaadsgacaWG0baaaiabggMi6kaaicdaaaa@3C1D@ , тогда из (1)

j=1 N q ij + Q i * =0 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaeWaaeaada aeWbqaaiaadghadaWgaaWcbaGaamyAaiaadQgaaeqaaOGaey4kaSIa amyuamaaDaaaleaacaWGPbaabaGaaiOkaaaaaeaacaWGQbGaeyypa0 JaaGymaaqaaiaad6eaa0GaeyyeIuoaaOGaayjkaiaawMcaaiabg2da 9iaaicdaaaa@457C@ , i= 1,N ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabg2 da9maanaaabaGaaGymaiaacYcacaWGobaaaaaa@3A28@ ; (7)

отсюда с учетом (3)

j=1 N β ij T j * T i * = Q i * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaabCaeaacq aHYoGydaWgaaWcbaGaamyAaiaadQgaaeqaaOWaaeWaaeaacaWGubWa a0baaSqaaiaadQgaaeaacaGGQaaaaOGaeyOeI0IaamivamaaDaaale aacaWGPbaabaGaaiOkaaaaaOGaayjkaiaawMcaaaWcbaGaamOAaiab g2da9iaaigdaaeaacaWGobaaniabggHiLdGccqGH9aqpcqGHsislca WGrbWaa0baaSqaaiaadMgaaeaacaGGQaaaaaaa@4BC9@ , i= 1,N ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabg2 da9maanaaabaGaaGymaiaacYcacaWGobaaaaaa@3A28@ . (8)

С учетом (4) - (8) можно определить задачу расчета приведенных коэффициентов теплообмена в общем виде, при заданных типовых значениях температур и тепловых потоков для всех взаимодействующих тепловых масс:

BTdiag(T)B1=Q MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaCOqaiabgw SixlaahsfacqGHsislcaqGKbGaaeyAaiaabggacaqGNbGaaiikaiaa hsfacaGGPaGaeyyXICTaaCOqaiabgwSixlaahgdacqGH9aqpcqGHsi slcaWHrbaaaa@4981@ , (9)

где 1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaCymaaaa@369F@  - единичная матрица размером N×1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOtaiabgE na0kaaigdaaaa@398A@ . В общем виде данное уравнение решить нельзя, т. к. при четырех и более элементах группы система будет недоопределена. В таком случае решение может быть получено с помощью постановки задачи минимизации функционала диссипации, где система (9), дополненная рядом дополнительных условий (симметрия матрицы B и т. п.), будет выступать в виде ограничений. Но во множестве практически полезных случаев можно рассматривать группу тепловых масс с некоторыми ограничениями.

Первое ограничение накладывается на количество связей в системе и определяется как

NE MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOtaiabgw MiZkaadweaaaa@3948@ , (10)

где N MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOtaaaa@36B8@  - количество тепловых масс; E MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyraaaa@36AF@  - количество связей в группе. Второе ограничение относится к случаю, когда в одну тепловую массу приходят два или более потоков. Пусть D k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiramaaBa aaleaacaWGRbaabeaaaaa@37CA@  - множество тепловых масс, которые привносят энергию k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@36D5@ -й тепловой массе, тогда

q ik T i * T k * = C k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGXbWaaSbaaSqaaiaadMgacaWGRbaabeaaaOqaaiaadsfadaqhaaWc baGaamyAaaqaaiaacQcaaaGccqGHsislcaWGubWaa0baaSqaaiaadU gaaeaacaGGQaaaaaaakiabg2da9iaadoeadaWgaaWcbaGaam4Aaaqa baaaaa@4230@ , i D k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabgI GiolaadseadaWgaaWcbaGaam4Aaaqabaaaaa@3A3C@ , (11)

где C k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGRbaabeaaaaa@37C9@  - постоянная величина, C k >0 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGRbaabeaakiabg6da+iaaicdaaaa@3995@ , откуда следует

β ik = C k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqOSdi2aaS baaSqaaiaadMgacaWGRbaabeaakiabg2da9iaadoeadaWgaaWcbaGa am4Aaaqabaaaaa@3C84@ , i D k MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyAaiabgI GiolaadseadaWgaaWcbaGaam4Aaaqabaaaaa@3A3C@ . (12)

Граничные условия описываются следующим образом: все N тепловых масс, описывающих объект, можно разделить на внешние и внутренние. Множество S MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36BD@  внешних тепловых масс включает в себя два подмножества: S IN MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4uamaaBa aaleaacaWGjbGaamOtaaqabaaaaa@388A@  - входные и S OUT MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4uamaaBa aaleaacaWGpbGaamyvaiaadsfaaeqaaaaa@3970@  - выходные тепловые массы. Для внешних входных тепловых масс энергия, привнесенная извне, будет положительной:

Q s >0 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGZbaabeaakiabg6da+iaaicdaaaa@39AB@ , s S IN MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4CaiabgI GiolaadofadaWgaaWcbaGaamysaiaad6eaaeqaaaaa@3B06@ , (13)

тогда сумма всей энергии, привнесенной извне,

Q = s S IN Q s MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacqGHris5aeqaaOGaeyypa0ZaaabuaeaacaWGrbWaaSbaaSqa aiaadohaaeqaaaqaaiaadohacqGHiiIZcaWGtbWaaSbaaWqaaiaadM eacaWGobaabeaaaSqab0GaeyyeIuoaaaa@42DA@ . (14)

Внешним выходным тепловым массам будет присущ отток энергии во внешнюю среду:

Q s <0 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGZbaabeaakiabgYda8iaaicdaaaa@39A7@ , s S OUT MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4CaiabgI GiolaadofadaWgaaWcbaGaam4taiaadwfacaWGubaabeaaaaa@3BEC@ , (15)

с учетом (14) распределение энергии во внешнюю среду определяется

Q s = γ s Q MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGZbaabeaakiabg2da9iabeo7aNnaaBaaaleaacaWGZbaa beaakiaadgfadaWgaaWcbaGaeyyeIuoabeaaaaa@3E6A@ , s S OUT MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4CaiabgI GiolaadofadaWgaaWcbaGaam4taiaadwfacaWGubaabeaaaaa@3BEC@ , (16)

где γ s MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeq4SdC2aaS baaSqaaiaadohaaeqaaaaa@38B0@  - доля суммарной привнесенной энергии, уходящая вовне от s-ной внешней выходной тепловой массы,

s S OUT γ s =1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaabuaeaacq aHZoWzdaWgaaWcbaGaam4CaaqabaaabaGaam4CaiabgIGiolaadofa daWgaaadbaGaam4taiaadwfacaWGubaabeaaaSqab0GaeyyeIuoaki abg2da9iaaigdaaaa@42A6@ . (17)

Таким образом, вся энергия, привнесенная в систему через внешние входные тепловые массы, убывает через внешние выходные тепловые массы. Этот же баланс соблюдается и во внутренних тепловых массах. Для v-той внутренней тепловой массы, взаимодействующей со множеством тепловых масс R v IN MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOuamaaDa aaleaacaWG2baabaGaamysaiaad6eaaaaaaa@3985@ , привносящих энергию, и отдающей энергию множеству тепловых масс R v OUT MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOuamaaDa aaleaacaWG2baabaGaam4taiaadwfacaWGubaaaaaa@3A6B@ , справедливо

i R v IN β iv T i * T v * = i R v OUT β vi T i * T v * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaabuaeaacq aHYoGydaWgaaWcbaGaamyAaiaadAhaaeqaaOWaaeWaaeaacaWGubWa a0baaSqaaiaadMgaaeaacaGGQaaaaOGaeyOeI0IaamivamaaDaaale aacaWG2baabaGaaiOkaaaaaOGaayjkaiaawMcaaaWcbaGaamyAaiab gIGiolaadkfadaqhaaadbaGaamODaaqaaiaadMeacaWGobaaaaWcbe qdcqGHris5aOGaeyypa0JaeyOeI0YaaabuaeaacqaHYoGydaWgaaWc baGaamODaiaadMgaaeqaaOWaaeWaaeaacaWGubWaa0baaSqaaiaadM gaaeaacaGGQaaaaOGaeyOeI0IaamivamaaDaaaleaacaWG2baabaGa aiOkaaaaaOGaayjkaiaawMcaaaWcbaGaamyAaiabgIGiolaadkfada qhaaadbaGaamODaaqaaiaad+eacaWGvbGaamivaaaaaSqab0Gaeyye Iuoaaaa@6080@ . (18)

Энергия может поступать в группу тепловых масс и убывать из нее двумя способами. Первый способ предполагает обмен энергией вследствие теплопередачи от более горячей к более холодной тепловой массе. При втором способе приток или отток энергии может происходить из-за замещения вещества, принятого за одну из тепловых масс модели. Когда в объекте моделирования присутствует проходящий через него теплоноситель, в качестве тепловой массы принимается его постоянный объем, находящийся непосредственно во взаимодействии с остальными тепловыми массами объекта моделирования. В этом случае необходимо учитывать замену теплоносителя в этом объеме (модель идеального вытеснения) и его температуру на входе. Из этого следует, что типичное значение температуры тепловой массы теплоносителя на выходе, используемое для определения коэффициентов β ij MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqOSdi2aaS baaSqaaiaadMgacaWGQbaabeaaaaa@398F@ , можно рассчитать по формуле

T W * = F W * c W T IN * + Q W * F W * c W MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaDa aaleaacaWGxbaabaGaaiOkaaaakiabg2da9maalaaabaGaamOramaa DaaaleaacaWGxbaabaGaaiOkaaaakiaadogadaWgaaWcbaGaam4vaa qabaGccaWGubWaa0baaSqaaiaadMeacaWGobaabaGaaiOkaaaakiab gUcaRiaadgfadaqhaaWcbaGaam4vaaqaaiaacQcaaaaakeaacaWGgb Waa0baaSqaaiaadEfaaeaacaGGQaaaaOGaam4yamaaBaaaleaacaWG xbaabeaaaaaaaa@496F@ , (19)

где T W * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaDa aaleaacaWGxbaabaGaaiOkaaaaaaa@3875@  - типичная температура тепловой массы-теплоносителя; F W * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOramaaDa aaleaacaWGxbaabaGaaiOkaaaaaaa@3867@  - расход теплоносителя в статике; T IN * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaDa aaleaacaWGjbGaamOtaaqaaiaacQcaaaaaaa@393A@  - температура замещающего теплоносителя в статике; c W MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4yamaaBa aaleaacaWGxbaabeaaaaa@37D5@  - удельная теплоемкость; Q W * MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaDa aaleaacaWGxbaabaGaaiOkaaaaaaa@3872@  - энергия, приходящая или уходящая на тепловую массу-теплоноситель, т. к. тепловая масса, принятая за теплоноситель с расходом, всегда будет являться внешней тепловой массой. Отнесение внешней тепловой массы-теплоносителя к входной или выходной будет зависеть от конкретного объекта моделирования и роли теплоносителя в нем.

Решая систему из уравнений (8), составленных для каждой тепловой массы системы в статике, с учетом граничных условий (13), (14), (17) и ограничений (10), (12), а также баланса для внешних (17) и внутренних (18) тепловых масс можно однозначно определить все коэффициенты β ij MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqOSdi2aaS baaSqaaiaadMgacaWGQbaabeaaaaa@398F@ .

Предложенный подход к построению математической модели позволяет описывать различные элементы котельной. При расчете коэффициентов в статике направления тепловых потоков известны и взаимодействие тепловых масс для каждого элемента может быть описано направленным графом. В динамике взаимодействия между тепловыми массами происходят в направлениях, которые определяются температурами тепловых масс.

Моделирование тепловых объектов котельной

Водогрейный котел. Основное назначение котла - это нагрев теплоносителя. Передача тепла к теплоносителю осуществляется от металла тепловоспринимающих поверхностей котла, которые, в свою очередь, нагреваются от продуктов сгорания топлива. В процессе работы часть тепла уходит на нагрев кожуха котла, который, в свою очередь, отдает тепло окружающей среде. Тепловые массы модели котла и тепловые потоки между ними изображены на рис. 2. Для расчета приведенных коэффициентов теплообмена между тепловыми массами можно использовать оценочные и типичные значения температур для каждой тепловой массы. При этом первая тепловая масса является внешней входной, т. к. на нее поступает энергия от горелки, а вторая и третья - внешними выходными, причем распределение выходных потоков зависит от КПД котла. Поток тепла к поверхностям нагрева определяется комбинацией воздействий: излучения пламени и теплопередачи от продуктов сгорания. При этом температура продуктов сгорания на порядок выше, чем температура металла поверхностей нагрева. Поэтому при расчете принято, что тепловой поток Q 1 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaaIXaaabeaaaaa@37A3@  не зависит от разности температур и определяется мощностью горелки, зависящей от расхода топлива.

 

Рис. 2. Представление котла: 1 - металл тепловоспринимающих поверхностей; 2 - теплоноситель; 3 - кожух; T 1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaIXaaabeaaaaa@39C8@ - T 3 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaIZaaabeaaaaa@39CA@  - температуры соответствующих тепловых масс; T 4 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaI0aaabeaaaaa@39CB@ - температура окружающей среды

 

Тепловой поток Q 2 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaaIYaaabeaaaaa@37A4@  определяется расходом теплоносителя через котел, соответствующим номинальному режиму работы котла; температура теплоносителя на выходе котла в номинальном режиме оценивается по формуле (19). Окружающая среда фактически выступает тепловой массой с бесконечно большой теплоемкостью, и ее температура T 4 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaI0aaabeaaaaa@37A9@  не изменяется, сколько бы тепла на нее ни поступило.

Теплообменный аппарат. Функция теплообменника - обмен тепла между греющим и нагреваемым теплоносителями. Теплоносители не взаимодействуют напрямую друг с другом, они передают тепло стенке, разделяющей потоки. Эксперименты показали, что повысить точность модели теплообменника позволяет представление стенки, разделяющей потоки, в виде двух взаимодействующих тепловых масс, обменивающихся энергией друг с другом и с внешней поверхностью, которая отдает тепло вовне. Тепловые массы модели теплообменника и тепловые потоки между ними изображены на рис. 3. Привнесение энергии в систему происходит через первую тепловую массу (объем греющего теплоносителя), а отток - через четвертую (объем нагреваемого теплоносителя) и пятую (внешнюю поверхность теплообменника, взаимодействующую с внешней средой). В модели теплообменника наблюдается полная аналогия с моделью котла: распределение выходных потоков зависит от КПД теплообменника, окружающая среда выступает тепловой массой с бесконечно большой теплоемкостью.

 

Рис. 3. Представление теплообменника:

1 - греющий теплоноситель; 2 - металл греющего теплоносителя; 3 - металл нагреваемого теплоносителя; 4 - нагреваемый теплоноситель; 5 - внешняя поверхность; T 1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaIXaaabeaaaaa@39C8@ - T 5 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaI1aaabeaaaaa@39CC@  - температуры соответствующих тепловых масс; T 6 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamivamaaBa aaleaacaaI2aaabeaaaaa@39CD@  - температура окружающей среды

 

Представление потребителей отопления и ГВС в виде графа имеет линейный вид: тепло, привносимое в систему теплоносителем полностью (по цепочке), переходит в последующие тепловые массы потребителей. Для потребителя отопления важен нагрев воздуха в помещениях с помощью радиаторов, отток же тепла происходит через стены помещений во внешнюю среду. Таким образом, модель потребителя отопления представляется цепочкой тепловых масс, взаимодействующих только с соседними тепловыми массами: теплоноситель, металл радиаторов, воздух в помещении, стены помещений, окружающая среда.

Для моделирования потребителя ГВС достаточно меньшего количества тепловых масс. Основной отбор тепла осуществляется за счет отбора некоторого объема теплоносителя и замещения его холодной водой.

Вычислительный эксперимент

Для оценки предложенного подхода была реализована компьютерная модель в пакете MATLAB Simulink. Параметры Simulink-модели идентифицировались по реальным данным эксплуатации котельной, состоящей из двух котлов мощностью в 1 МВт каждый, с КПД 90 %, подключенных параллельно; каждый котел имеет собственную систему управления температурой теплоносителя на выходе. Она представлена в виде двухступенчатого регулятора, управляющего мощностью горелки (0 %, 70 %, 100 %). Выбор количества работающих котлов реализован в виде отдельной подсистемы управления. Пластинчатый противоточный теплообменник, используемый в котельной, имеет мощность 600 кВт и КПД 90 %. Моделируемая нагрузка - это 5700 м ² MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2Caerbov2D09 MBdbqedmvETj2BSbqeduuDJXwAKbYu51MyVXgarqqr1ngBPrgifHhD YfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9 Ff0dmeaabaqaciaacaWaaeqabaqabeGaeaaakeaaruWqHXwAIjxAaG qbaKqzGfaeaaaaaaaaa8qacaWFYcaaaa@3AC4@ кирпичных жилых помещений и 3300 м ² MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2Caerbov2D09 MBdbqedmvETj2BSbqeduuDJXwAKbYu51MyVXgarqqr1ngBPrgifHhD YfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9 Ff0dmeaabaqaciaacaWaaeqabaqabeGaeaaakeaaruWqHXwAIjxAaG qbaKqzGfaeaaaaaaaaa8qacaWFYcaaaa@3AC4@ панельных жилых помещений, находящихся в умеренно-континентальном регионе. Необходимые затраты на отопление [21] при температуре теплоносителя 95 °C и температуре наружного воздуха -35 °C составляют 70 и 100 Вт/м ² MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2Caerbov2D09 MBdbqedmvETj2BSbqeduuDJXwAKbYu51MyVXgarqqr1ngBPrgifHhD YfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9 Ff0dmeaabaqaciaacaWaaeqabaqabeGaeaaakeaaruWqHXwAIjxAaG qbaKqzGfaeaaaaaaaaa8qacaWFYcaaaa@3AC4@ для кирпичных и панельных помещений соответственно. На данных площадях находятся около 180 квартир.

Гидравлическая система для данной котельной станции представлена в виде функций, реализующих моментальные перераспределения потоков при изменении положений трехходовых клапанов, управление которыми происходит с помощью релейных регуляторов и по отклонению от заданных значений температур. Для отопления температура воды на потребителя определяется погодным графиком, для контура ГВС температура на выходе из котельной фиксирована и составляет 55 °C. Расход теплоносителя в гидравлической системе определяется насосами и составляет 50 т/ч для потребителя отопления и 20 т/ч для внутреннего контура ГВС. Расход воды на потребителя ГВС принят постоянным и составляет 10 т/ч. Общий объем теплоносителя для контура потребителя отопления составляет 30 м ³ MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2Caerbov2D09 MBdbqedmvETj2BSbqeduuDJXwAKbYu51MyVXgarqqr1ngBPrgifHhD YfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9 Ff0dmeaabaqaciaacaWaaeqabaqabeGaeaaakeaaruWqHXwAIjxAaG qbaKqzGfaeaaaaaaaaa8qacaWFZcaaaa@3AC5@ , для контура ГВС - 10 м ³ MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2Caerbov2D09 MBdbqedmvETj2BSbqeduuDJXwAKbYu51MyVXgarqqr1ngBPrgifHhD YfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9 Ff0dmeaabaqaciaacaWaaeqabaqabeGaeaaakeaaruWqHXwAIjxAaG qbaKqzGfaeaaaaaaaaa8qacaWFZcaaaa@3AC5@ .

В качестве примера результатов моделирования котла рассматривается его нагрев и остывание при отсутствии расхода теплоносителя. На рис. 4 приведены температуры тепловых масс модели котла (см. рис. 2), горелка была включена на полную мощность в течение 370 с, после чего выключена.

При включенной горелке быстрее всего нагревался металл поверхностей нагрева котла, теплоноситель и кожух нагревались медленнее, после отключения горелки температура металла достаточно быстро снижается из-за постоянного теплообмена с теплоемким теплоносителем. После 800 с температура металла становится меньше температуры теплоносителя. Видно, что теплоноситель постепенно остывает, взаимодействуя с остывающим металлом поверхностей нагрева, которые, в свою очередь, остывают от оттока тепла к кожуху котла. Кожух котла взаимодействует с внешней средой - воздухом внутри котельной с температурой 25 °C.

Остывание котла, изображенное на рис. 4, длится около 24 ч, что соответствует реальным данным по температуре теплоносителя после остановки котла мощностью 1 МВт.

 

Рис. 4. Нагрев и остывание котла

 

Результаты моделирования котельной вместе с ее тепловой нагрузкой представлены на рис. 5, 6. Исходные данные представляют собой четыре фрагмента разной продолжительности. Для проведения эксперимента фрагменты были соединены в один с ограничением на скорость изменения температур в 1 °C/ч. На рис. 6 для наглядности представлен период, когда работал один котел и проводилась беспрерывная запись параметров котельной. При моделировании учитывалось изменение температуры наружного воздуха, которая влияла на потребление тепла на отопление.

Графики (см. рис. 5, 6) свидетельствуют о корректности построения моделей ввиду соблюдения тепловой нагрузки и динамики реальной системы, а также верной реализации системы управления. Точность моделирования определялась как средняя абсолютная ошибка (MAE) и максимальная абсолютная ошибка (MaxAE) между модельными и реальными данными (см. таблицу). Полученное значение средней абсолютной ошибки соизмеримо с погрешностью датчиков, использующихся в реальной системе.

 

Оценка точности модели

Параметр

MAE

MaxAE

16 суток

6 суток

16 суток

6 суток

T подачи к потребителю контура отопления, °C

0,4743

0,3890

5,9820

4,2910

T возврата от потребителя контура отопления, °C

0,2564

0,1929

2,3605

1,4354

 

Рис. 5. Моделирование за 16 суток

 

Рис. 6. Моделирование за 6 суток

 

На графиках (см. рис. 5, 6) наблюдаются кратковременные расхождения (выбросы) между поведением реальной системы и модели, которые могут быть объяснены в первую очередь недостаточной точностью моделирования потребления в контуре ГВС: увеличение потребления в утренние и вечерние часы. Разработка более точной модели потребления ГВС имеет смысл только для решения задачи оптимизации локальной автоматической системы поддержания температуры подачи теплоносителя в контуре ГВС и в данной работе не рассматривалась.

Заключение

Разработанный метод моделирования водогрейных котельных малой мощности на основе тепловых масс позволяет получить адекватное описание тепловых процессов с учетом взаимодействия элементов между собой и с окружающей средой. Предложенная математическая постановка обеспечивает универсальность и масштабируемость модели, а также возможность ее применения к различным элементам тепловой схемы, включая котлы, теплообменники и потребителей тепла. Применение численного подхода на базе законов сохранения массы и энергии позволило избежать чрезмерной математической сложности, характерной для аналитических моделей, сохранив при этом необходимую точность описания и обеспечив возможность реализации динамической модели системы.

Результаты моделирования, полученные в пакете MATLAB Simulink, показали высокую степень соответствия с реальными эксплуатационными данными котельной, что подтверждается малыми значениями среднеквадратичной абсолютной ошибки температур теплоносителя. Проведенный вычислительный эксперимент подтвердил возможность применения модели для задач анализа, оптимизации и прогнозирования работы котельных установок. Таким образом, предложенный метод может быть рекомендован в качестве инженерного инструмента для проектирования и совершенствования автономных систем теплоснабжения без необходимости привлечения ресурсоемких вычислительных средств или сложных экспериментальных исследований.

×

Об авторах

Л. П. Строкин

Самарский государственный технический университет

Автор, ответственный за переписку.
Email: strokin.samara@mail.ru

магистр

Россия, 443100, г. Самара, ул. Молодогвардейская, 244

И. А. Данилушкин

Самарский государственный технический университет

Email: idanilushkin@mail.ru

кандидат технических наук, доцент, доцент кафедры автоматики и управления в технических системах

Россия, 443100, Самара, ул. Молодогвардейская, 244

Список литературы

  1. Милевская Т.А. Анализ технологий, повышающих энергоэффективность котельных // Научно-исследовательский центр «Technical Innovations». 2021. № 5. С. 215–219. EDN: SRMSBM
  2. Глухова А.В., Трускова Е.В. Анализ основных направлений развития системы теплоснабжения в Кузбассе // Россия молодая: Сб. материалов XIV Всерос. науч.-практ. конф. с межд. участием, Кемерово, 19–21 апр. 2022 г. Кемерово: КузГТУ им. Т.Ф. Горбачева, 2022. С. 21206.1–21206.5. EDN: RBTTVO
  3. Babadi K.A. et al. Numerical simulations and an experimental study for optimal design of a 1500 kW water-tube condensing boiler // Energy. 2024. Vol. 307. P. 132675. doi: 10.1016/j.energy.2024.132675. EDN: PKCKFB
  4. Якупова К.Р. Исследование эффективности теплообмена при изменении форм поперечного сечения конвективных труб газового котла КВГА-0,8 // Инновационные науч. исследования. 2021. № 5–1(7). С. 54–62. doi: 10.5281/zenodo.5007206. EDN:XRISEX
  5. Arowoiya V.A., Moehler R.C., Fang Y. Digital twin technology for thermal comfort and energy efficiency in buildings: A state-of-the-art and future directions // Energy and Built Environment. 2024. Vol. 5, № 5. Pр. 641–656. doi: 10.1016/j.enbenv.2023.05.004. EDN: GMTKEU
  6. Пащенко Ф.Ф., Пикина Г.А. Основы моделирования энергетических объектов. М.: Физматлит, 2011. 464 с. EDN: UGPDPP
  7. Аветисян А.Р., Пащенко А.Ф., Пащенко Ф.Ф. и др. Теплогидравлические модели оборудования электрических станций. М.: Физматлит, 2013. 448 с. EDN: UGLECL
  8. Hosain M.L. Fluid flow and heat transfer simulations for complex industrial applications: from Reynolds averaged Navier-Stokes towards smoothed particle hydrodynamics: дис. ... канд. техн. наук. Mälardalen Univ., 2018.
  9. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Машиностроение, 1984.
  10. Кудинов А.А. Математическое моделирование теплового и термонапряженного состояний элементов конструкций теплогенерирующих установок: дис. ... д-ра техн. наук. Казань, 1999. 474 с. EDN: QDALTF
  11. Раппопорт Э.Я. Структурное моделирование объектов и систем управления с распределенными параметрами: учеб. пособие. М.: Высшая школа, 2003. 302 с.
  12. Ramallo-González A.P., Eames M.E., Coley D.A. Lumped parameter models for building thermal modelling: An analytic approach to simplifying complex multi-layered constructions // Energy and Buildings. 2013. Vol. 60. Pр. 174–184. doi: 10.1016/j.enbuild.2013.01.014
  13. Фокин В.М., Бойков Г.П., Видин Ю.В. Основы технической теплофизики. М.: Машиностроение-1, 2004. EDN: QJMPQB
  14. Кузнецов В.А. Математическое моделирование процессов тепло- и массообмена для перспективных технологий энергетического использования угольного топлива: дис. ... канд. техн. наук. Красноярск, 2019. 125 с. EDN: ZKECKC
  15. Michel A., Kugi A. Accurate low-order dynamic model of a compact plate heat exchanger // International Journal of Heat and Mass Transfer. 2013. Vol. 61. Pр. 323–331. doi: 10.1016/j.ijheatmasstransfer.2013.01.072
  16. He K.-L. et al. A transient heat current model for dynamic performance analysis and optimal control of heat transfer system // International Journal of Heat and Mass Transfer. 2019. Vol. 145. P. 118767. doi: 10.1016/j.ijheatmasstransfer.2019.118767
  17. Нгуен Х.Д. Моделирование теплоэнергетических установок на основе теории дифференциально-алгебраических уравнений в частных производных: автореф. дис. ... канд. техн. наук: 05.13.18. Иркутск, 2014. 22 с. EDN: ZPNXBH
  18. Барочкин А.Е., Жуков В.П., Барочкин Е.В., Ледуховский Г.В. Матричное представление модели тепловой схемы электрической станции // Вестник ИГЭУ. 2018. № 6. С. 66–72. doi: 10.17588/2072-2672.2018.6.066-072. EDN: VPQIQN
  19. Лыков А.В. Теория теплопроводности. М., 1967.
  20. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. 3-е изд. М.: Наука, 1975. 768 с. EDN: WGDAWV
  21. СНиП 23-02-2003. Тепловая защита зданий / Госстрой России. Изд. офиц. М.: ФГУП ЦПП, 2004. 26 с.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Схема котельной: Н1, Н2, Н3 - группы насосов; К1, К2 - трехходовые клапаны

Скачать (82KB)
3. Рис. 2. Представление котла: 1  металл тепловоспринимающих поверхностей; 2 - теплоноситель; 3 - кожух; -  - температуры соответствующих тепловых масс; - температура окружающей среды

Скачать (23KB)
4. Рис. 3. Представление теплообменника: 1 - греющий теплоноситель; 2 - металл греющего теплоносителя; 3 - металл нагреваемого теплоносителя; 4 - нагреваемый теплоноситель; 5 - внешняя поверхность; -  - температуры соответствующих тепловых масс;  - температура окружающей среды

Скачать (35KB)
5. Рис. 4. Нагрев и остывание котла

Скачать (131KB)
6. Рис. 5. Моделирование за 16 суток

Скачать (151KB)
7. Рис. 6. Моделирование за 6 суток

Скачать (137KB)

© Строкин Л.П., Данилушкин И.А., 2025

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.