Применение пространственной аппроксимации распределенных сигналов при моделировании теплообменных аппаратов

Размещено на http://www.allbest.ru/

Применение пространственной аппроксимации распределенных сигналов при моделировании теплообменных аппаратов

М.А. Гусева

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

Ключевые слова: теплообменный аппарат, моделирование, система с распределенными параметрами, аппроксимация, численно-аналитическая модель.

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

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

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

, (1)

, (2)

, .

С начальными и граничными условиями:

, ,

для прямотока:

,

для противотока:

,

где Ti(x,t) - распределение температуры i-го потока по длине, град, ; vi - скорость i-го потока жидкости, м/c; вi - приведённый коэффициент конвективного теплообмена i-го потока с соседним потоком, 1/с; L - длина теплообменного аппарата, м.

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

Для потока, направление которого совпадает с осью абсцисс, передаточная и стандартизирующая функции имеют вид [1]:

, (3)

, (4)

где p - оператор преобразования по Лапласу, - единичная функция:

Tenv(x,p) - произвольное распределение теплоисточников по длине теплообменного аппарата (рис. 1). В уравнениях (1), (2) в роли Tenv(x,p) выступает температурное распределение соседнего потока.

Рис. 1. Пример распределения теплоисточников Tenv(x,p)

Распределенный сигнал можно аппроксимировать полиномиальной зависимостью N-го порядка (см. рис. 1):

; , (5)

где N - степень полинома, M - количество точек.

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

, где . (6)

Коэффициенты аппроксимации (5) могут быть найдены в общем случае по методу наименьших квадратов [4]:

, (7)

где - вектор неизвестных коэффициентов линеаризации, ;

- матрица размерности MЧ(N+1), ;

- вектор измеренных значений температур, .

Температурное распределение потока теплообменного аппарата может быть найдено по формуле [1]:

. (8)

Тогда значение температуры потока в определенной точке xm с учетом (3) - (5) рассчитывается по формуле:

Используя известное выражение для размещений , можно свернуть получившееся выражение, тогда температурное распределение примет вид:

(9)

На основе выражения (9) модель теплообменного аппарата, описываемого системой уравнений (1) - (2), можно представить в виде структурной схемы (рис. 2); , - векторы настраиваемых параметров теплообменника (скорость потока, длина теплообменника, коэффициент теплопередачи и т.д.).

Рис. 2. Структурная схема теплообменного аппарата

В блоке вычисления коэффициентов аппроксимации (БВКА) реализован алгоритм вычисления коэффициентов аппроксимации по выражениям (6) или (7).

Синтез и анализ численно-аналитических моделей теплообменного аппарата с различными способами аппроксимации распределенного сигнала. Полученная структурная схема (см. рис. 2) может быть легко реализована в пакетах численного моделирования динамических систем. Ранее, в работах [5, 6], на основе предложенного метода была разработана и реализована численно-аналитическая модель противоточного теплообменного аппарата с аппроксимацией распределенного сигнала линейной зависимостью вида

. (10)

Коэффициенты аппроксимации k1, k0 вычисляются по измеренным значениям температур теплообменного аппарата в фиксированных точках - на входе и выходе теплообменника.

В рамках проделанной работы были разработаны численно-аналитические модели противоточных теплообменников:

с аппроксимацией распределенного сигнала полиномиальной зависимостью третьего порядка:

. (11)

Коэффициенты аппроксимации k3, k2, k1, k0 вычисляются по выражению (6). Значения температур измерены в фиксированных точках - на входе и выходе теплообменника; на 1/3 и 2/3 его длины;

с аппроксимацией распределенного сигнала полиномиальной зависимостью второго порядка:

. (12)

Коэффициенты аппроксимации k2, k1, k0 вычисляются методом наименьших квадратов (7) по тем же четырем точкам (, , , ).

При расчете теплообменного аппарата были заданы следующие конструктивные и теплофизические характеристики:

- длина теплообменника L=6 м;

- скорость греющего потока воды v1=1,7 м/с;

- скорость нагреваемого потока воды v2=3 м/с;

- температуры теплоносителей на входе теплообменника: Tin 1=100 град; Tin 2=20 град.

Приведённые коэффициенты конвективного теплообмена между потоками: в1=0,2153 1/c; в2=0,3821 1/c.

В результате исследований трех реализованных моделей теплообменных аппаратов (10) - (12) были получены кривые температурного распределения на выходе нагреваемого потока после нанесения возмущения на входе (рис. 3) и переходный процесс температуры греющего потока при различных способах аппроксимации температурного распределения нагреваемого потока (рис. 4).

Рис. 3. Температурное распределение нагреваемого скачка на выходе после нанесения возмущения температурой на входе:

а - 0,6 сек после нанесения возмущения; б - 0,8 сек после нанесения возмущения;

в - 1,2 сек после нанесения возмущения; г - 1,4 сек после нанесения возмущения

Оценка точности предложенных способов аппроксимации распределенного сигнала была проведена по выражению (13) в разные моменты времени (значения температур были сняты в 11 точках - с шагом ): теплообмен распределенный сигнал

, . (13)

Полученные данные были сведены в таблицу.

Оценка погрешности аппроксимации распределенного сигнала

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

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Бутковский А.Г. Структурная теория распределённых систем. - М.: Наука, 1977.

2. Данилушкин И.А., Лежнев М.В. Структурное представление процесса теплообмена при встречном направлении взаимодействующих потоков // Вестник Самарского гос. техн. ун-та. Сер. Технические науки. Вып. 1(19)-2007. - Самара: СамГТУ, 2007. - С. 16-22.

3. Данилушкин И.А., Тимофеева О.Н. Структурное представление и передаточные функции прямоточного теплообменника// Вестник Самарского государственного технического университета. Сер. Технические науки. Вып. 1(21)-2008. - Самара: СамГТУ, 2008. - С. 9-13.

4. Линник Ю.В. Метод наименьших квадратов и основы теории обработки наблюдений. - М.: Гос. изд-во физико-математической литературы, 1958.

5. Данилушкин И.А., Гусева М.А. Построение численно-аналитических моделей теплообменных аппаратов // Международ. науч. конф. «Проблемы управления, передачи и обработки информации - АТМ-ТКИ-50»: Сб. трудов / Под ред. А.Г. Александрова, М.Ф. Степанова. - Саратов: Сарат. гос. техн. ун-т, 2009. - C. 168-170.

6. Данилушкин И.А., Гусева М.А. Численно-аналитическая модель нагрева потока жидкости // Высокие технологии, исследования, промышленность. Т. 1: IX Международная научно-практ. конференция «Исследование, разработка и применение высоких технологий в промышленности», 22-23 апреля 2010: Сб. трудов / Под. ред. А.П. Кудинова. - СПб.: Изд-во политехн. ун-та, 2010. - С. 297-298.

Размещено на Allbest.ru