Основы метода наименьших квадратов

Содержание Введение Задание

1.Постановка задачи

1.1Метод наименьших квадратов

1.2Метод Брандона

1.3Реактор идеального вытеснения

1.4Синтез оптимальных тепловых систем

2. Расчетная часть

2.1Расчет значений К₀ и Е в уравнении Аррениуса с использованием метода наименьших квадратов

2.2Расчет статистической модели абсорбера методом Брандона

2.3Расчет абсорберов и реакторов идеального вытеснения

2.4Синтез системы теплообмена Выводы Список использованных источников Приложение абсорбер реактор вытеснение теплообмен

Введение

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

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

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

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

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

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

ХТС характеризуется структурой, т. е. составом элементов и их взаимными связями. Отличительная черта современных ХТС — большое количество внутренних связей, в том числе и более полного использования сырья, вторичных энергоресурсов и стремлением к созданию более гибких ХТС.

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

Введение

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

Задание Вариант № 4

Требуется синтезировать XTC работающую по следующей технологии:

Смесь, состоящая из компонентов A, Б и инертного компонента нагревается в системе теплообмена до 415 ⁰C, поступает в реактор, где протекает обратимая реакция:

A + 0,5 Б — C + q;

где q — тепловой эффект реакции q = 21 200кал/моль Реакция характеризуется константой скорости K=f (t) и константой равновесия Kp=f (t), для которых имеются экспериментальные данные.

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

После прохождения 4 реакторов смесь поступает в абсорбер для выделения компонента C и затем проходит 1 реактор и второй абсорбер.

Таким образом, операторная схема выглядит следующим образом:

Заданы температуры на входе в реакторы и абсорберы, и объемы реакторов и абсорберов. Заданы также плотности орошения в абсорберах, температура, расход и концентрации компонентов исходной смеси.

t₀ = 50 °C t₁ = 415 °C t₂ = 485 °C t₃ = 410 °C t₄ = 425 °C t₅ = 405°С

t_a1 = 180 °C t_a2 = 195 °C Расход смеси на входе в систему-132 000 м³/ч Объемы реакторов: V₁ = 79 м³, V₂ = 65 м³, V₃ = 40 м³, V₄ = 55 м³, V₅ = 35 м³.

Объемы абсорберов: V_a1 = 25 м³, V_a2 = 28 м³ .

Плотность орошения в первом абсорбере —16 м³/м²

Плотность орошения во втором абсорбере —16 м³/м²

Реакторы описываются моделями идеального вытеснения. Абсорберы описываются статистическими моделями по экспериментальным данным. Скорость реакции в реакторе описывается уравнением:

W=K (t)· a·b/(a+0,8·c)·(1-(c/K_p(t)·a·b ^0,5))²),

где а, b, с — концентрация компонентов (об.доли):

а = 0,085 об. доли

b = 0,085 об. доли с = 0,00 об. доли При построении системы теплообмена могут использоваться пар и вода со следующими характеристиками. Начальная температура воды = 20 ⁰C, конечная температура воды не более 90⁰C, температура пара равна 460 ⁰C, теплота конденсации греющего пара = 520 ккал/кг. Стоимость воды 0,7 ус.д.ед./кг. Стоимость пара =0,005 ус.д.ед./кг. Коэффициенты теплопередачи:

в теплообменниках 19 ккал/(м² · ч · °С) в холодильнике 20 ккал/(м² · ч · °С) в нагревателях 22 ккал/(м² · ч · °С) Теплоемкость в реакционной смеси 0,33 ккал/(м³· °С) Время работы установки 8800 час/год.

Для получения значений К₀ и Е в уравнении Аррениуса использовать данные приложения 1.

K_p(t) = 10^(49 051/T-4,6455)

Для получения статической модели абсорбера использовать данные приложения 2 и метод Брандона.

1. Постановка задачи

1.1 Метод наименьших квадратов

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

Для построения аппроксимирующей кривой используют уравнение

y = ?(x, a₀, a₁, a₂,…, a_m)(1.1)

Задачу определения параметров уравнения регрессии (1.1) решают обычно с использованием метода наименьших квадратов (MHK). Согласно методу наилучшими значениями параметров будут те, при которых сумма квадратов отклонений расчетных величин a₀,…, а_т от экспериментальных будет наименьшей, то есть:

Ф = (y_0i — ?(x, a₀, a₁, a₂,…, a_m))² = min

Условием минимума функции является равенство нулю ее частных производных по всем аргументам:

(1.2)

В теории метода система (1.2) называется нормальной системой: для неё число уравнений равно числу неизвестных.

Для рассматриваемой функции Ф условием минимума будет После преобразований нормальная система принимает вид:

(1.3)

Решение системы алгебраических уравнений (1.3) позволяет определить значения коэффициентов для конкретного вида зависимости ?(x, a₀, a₁, a₂,…, a_m) .

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

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

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

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

В общем случае современный технологический процесс представляется в виде многомерного объекта, блок-схема которого приведена на рис. 1.1. На объект действуют вектор входных параметров X, составляющие которого {x₁, x₂,…, x_l}, и вектор управления Z, составляющие которого{z₁, z₂,…z_k}.Выходные параметры {y₁, y₂,…y_p} составляют вектор выходных параметров Y. Общий вид статистической модели многомерного технологического объекта можно записать в виде системы алгебраических уравнений (1.4) или в векторной форме (1.5):

где X, Y — векторы входных и выходных параметров объекта.

B системе (1.4) параметры управления учтены как входные параметры

x_l+1, x_l+2,…, x_m (m = l+k)

Рис. 1 1. Блок-схема многомерного технологического процесса

B данной курсовой работе для построения модели многомерного технологического объекта используется метод Брандона.

Сущность метода заключается в следующем. Предполагается, что функция F₁{x₁, x₂,…, x_m} в системе (1.4) является произведением функций от входных параметров, то есть

y = y ?₁(x₁)?₂(x₂)???_m(x_m)(1.6)

или в более удобной форме:

y = y? ?_k(x_k)

где y — расчетное значение i — го выходного параметра; - средняя величина экспериментальных значений i — го выходного параметра; n — количество опытов в исходной выборке.

При использовании метода Брандона большое значение имеет порядок следования функций в уравнении (1.6). Чем больше влияния оказывает фактор на выходной параметр, тем меньшим должен быть его порядковый номер в указанном уравнении.

Оценить степень влияния к — го фактора на выходной параметр можно по величине частного коэффициента множественной корреляции:

(1.7)

где — величина частного коэффициента корреляции, учитывающая влияние k — го фактора на выходной параметр y при условии, что влияние всех прочих факторов исключено. D — определитель матрицы, построенной из парных коэффициентов корреляции. Матрица имеет вид:

D_lk — определитель матрицы с вычеркнутыми первой строкой и k — м столбцом.

D_llD_kk — определитель матрицы с вычеркнутыми первой и k — ой строками и первым и k — м столбцами соответственно. r_xy — парные коэффициенты корреляции определяемые по формуле:

Коэффициент корреляции по абсолютной величине не превышает единицы.

— 1? r_xy? 1.

Чем ближе абсолютное значение коэффициента ¦ r_xy ¦ к единице, тем сильнее линейная связь между величинами. Следует отметить, что коэффициент корреляции одинаково отмечает долю случайности и криволинейность связи между х и у. Зависимость x и у может быть близкой к функциональной, но существенно нелинейной; коэффициент корреляции при этом будет значительно меньше единицы.

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

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

B уравнении (1.6) каждая из функций ?₁(x₁)???_m(x_m) принимается либо линейной, либо нелинейной (степенной, показательной, экспоненциальной и т. д.). Прежде чем определять вид первой зависимости, следует представить исходные экспериментальные значения выходного параметра в каждом опыте y_эj в безразмерной форме y_э0j

где y — средняя величина выходного параметра.

Таким образом, исходными данными для поиска первой зависимости будут нормированные значения вектора выходных параметров y₀ и опытные значения первого влияющего фактора. Выбрав зависимость y₁ = ?₁(x₁) с помощью метода наименьших квадратов, определяют остаточный показатель У_э для каждого наблюдения:

Предполагая, что y_э1 не зависит от x₁, а зависит от x₂,…x_m, выбирают зависимость от второго фактора. Исходные данные для поиска — остаточный показатель у_э1 и опытные значения второго фактора. Получив расчетную зависимость y₂ = ?₂(x₂), находят остаточный показатель y_э2 для каждого наблюдения:

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

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

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

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

Математическая модель — это система уравнений, которая устанавливает связь входных и выходных параметров реактора.

aA + bB >cC

v = K· C_A^a· C_B^b

В нашем случае скорость реакции в реакторе описывается уравнением:

W=K (t)· a·b/(a+0,8·c)·(1-(c/K_p(t)·a·b ^0,5))²)

Для определения скорости реакции по каждому веществу для многоступенчатых химических реакций составляется стехиометрическая матрица размером m на n, где m — число стадий, n — число компонентов.

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

Скорость по i — му компоненту будет представлять собой сумму компонентов i — го столбца.

v_A = -2w_I

v_B = -w_I

v_C = 2w_I

v_D = w_II

Математическое описание:

К этой системе необходимо добавить уравнение теплового баланса:

где — коэффициенты адиабатического разогрева.

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

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

1.4 Синтез оптимальных тепловых систем

B наиболее традиционной постановке задача синтеза тепловых систем (TC) формулируется следующим образом: имеются т горячих и п холодных технологических потоков, которые называют основными технологическими потоками. Для каждого из этих потоков заданы начальные температуры t_нi^г и t_нj^х, конечные температуры t_кi^г, t_кj^х и значения водяных эквивалентов (произведение расхода на удельную теплоемкость) W_i^г, W_j^х. Здесь i=l, 2,…, m; j=l, 2,…, n. Индексы «г» и «x» относят соответствующую величину к горячему и холодному потокам.

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

Для решения задачи синтезируемую TC разделяют на две подсистемы: внутреннюю (рекуперативную), где и в теплообмене участвуют только основные технологические потоки;

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

При этом внешняя подсистема используется только тогда, когда во внутренней подсистеме не удается получить заданные конечные температуры.

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

З_пр = Е_н(З₁ + З₂) + З₃

З₁ — затраты на рекуперативные теплообменники, руб;

З₂ — затраты на вспомогательные теплообменники, руб;

З₃ — затраты на вспомогательные теплоносители, руб;

Е_н — нормативный коэффициент эффективности (E_н = 0,12).

Если во внутренней подсистеме используются К₁ теплообменных аппаратов, а во внешней — l₁, то:

При расчете стоимости i — го теплообменника любой подсистемы в данной работе используется зависимость где F — площадь поверхности теплообмена соответствующего i — го теплообменника, м²

а — стоимостной коэффициент, зависящий от типа теплообменника. Затраты на вспомогательные теплоносители определяются по формуле:

Где? — продолжительность годовой эксплуатации системы ч/год;

Ц_р — стоимость i — го вспомогательного теплоносителя,

G_pl — расходуго вспомогательного теплоносителя в l — м вспомогательном теплообменнике, кг/ч.

При синтезе TC используются формулы:

где Q — тепловая нагрузка теплообменника; k— соответствующий коэффициент теплопередачи.

Средняя разность температур для теплообменника:

где ?t_б и? t_м — разности температур на концах теплообменника.

Тепловая нагрузка теплообменника, или количество тепла, переданное в одном аппарате, определяется на основе концепции передачи максимально возможного тепла при минимально допустимой разности температур на концах теплообменника? t_min

Если (t_н^г — t_н^х — ?t_min)<0, то теплообмен невозможен.

Если t_к^г — t_н^х min, то t_к^г = t_н^х + ?t_min, иначе t_к^г = t_к^г

Если t_н^г — t_к^х< ?t_min, то t_к^х = t_н^г — ?t_min, иначе t_к^х = t_к^х.

Q_x = W_x (t_к^х — t_н^х)

Q_г = W_г (t_н^г — t_к^г)

Q = min (Q_х, Q_г)

t_к^х = t_н^х — Q/W_х

t_к^г = t_н^г — Q/W_г

Задача синтеза TC решается путем формирования множества возможных комбинаций исходных горячих и холодных потоков для проведения физически реализуемых операций теплообмена в теплообменном аппарате. Для этой цели строят таблицу пар взаимодействующих потоков исходя из условия Q > max. Из таблицы пар выбирается пара потоков, вступающих во взаимный теплообмен. Если в результате теплообмена данные потоки достигли заданных конечных температур, то они исключаются из рассмотрения. Иначе, начальным температурам этих потоков присваиваются значения конечных температур результирующих потоков, после чего таблица пар перестраивается, и выбирается новая пара потоков. Данная операция производится до тех пор, пока не останется потоков, способных вступать во взаимный теплообмен, или все потоки достигнут требуемых конечных температур.

При необходимости для достижения заданных конечных температур в теплообменных системах используются вспомогательные тепло — и хладагенты.

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

2. Расчетная часть

2.1 Расчет значений K_{0 и} E в уравнении Аррениуса с использованием метода наименьших квадратов Для расчета значений K₀ и E в уравнении Аррениуса использован метод наименьших квадратов. Вычисления проведены в Microsoft Excel:

Таблица 2.1 Зависимость константы скорости от температуры

Были получены значения К₀ =712 085 и Е = 82 914,21.

B соответствии с полученными данными уравнение Аррениуса имеет вид:

К = 712 085 · exp (82 914,21/RT)

2.2 Расчет статистической модели абсорбера методом Брандона Влияющие факторы:

Т_вх — температура на входе в абсорбер, ?C

с — плотность орошения, м³/м²

?-объем, м³

Выходные параметры:

Т_вых — температура на выходе из абсорбера, ?C

y — степень абсорбции, %.

На подготовительном этапе была создана таблица с экспериментальными данными (приложение 2).

Первый этап. Ранжирование влияющих факторов.

а). Построение матрицы коэффициентов парной корреляции.

Элементами матрицы являются коэффициенты парной корреляции, которые определяются с помощью надстройки Пакет анализа (на примере расчета Т_вых).

Таблица 2.2

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

Аналогично рассчитываем:

Частные коэффициенты множественной корреляции рассчитываются по формуле (1.7). Сведем вычисления в таблицу (2.3)

Таблица 2.3

В соответствии с убыванием величины коэффициентов корреляции определяем порядок расположения влияющих факторов в уравнении (1.6).

Таблица2.4 Экспериментальные и расчетные данные для построения статистической модели абсорбера для Т_вых

Второй этап. Выбор зависимостей выходных параметров от влияющих факторов.

Строим график зависимости T_вых / T_вых(ср) от с, добавим несколько линий тренда, найдем вид зависимости и искомые коэффициенты для с.

Рис. 2.2 График зависимости от с и линии тренда.

Выбираем зависимость с наибольшим значением R²: T_вых (V)= -0.0065· V+1.1602

Далее делим T _вых/T_вых(ср) на T (V). Построим график полученной функции от с. Далее все аналогично.

Вид зависимостей:

T_вых © =1,0612· e^-0,0035с

T_вых (T_вх) = 0,0004· T_вх+0.8367

Для y все операции аналогичны.

Таким образом, полученная экспериментально-статистическая модель имеет следующий Вид (уравнение (1.9)):

Т_вых = 333,2136· (-0.0065· V+1.1602)·(1,0612· e^-0,0035с)·(0,0004· T_вх+0.8367)

y = 91,56 864· (2,9903·е^{-0,0024Твх})·(0,7408·е^0,0121V)·(0,0089·с+0,8572)

Таблица 2.5 Экспериментальные и расчетные данные для построения статистической модели абсорбера для у

Третий этап. Проверка адекватности модели.

По формуле (1.9) определяем расчетные значения и сравниваем их с экспериментальными (см. таблицы 2.4 и 2.5). Для оценки точности аппроксимации найденных зависимостей вычислены корреляционные отношения з и средние относительные ошибки е:

Таблица 2.6

Таким образом, построенные статистические модели хорошо описывают экспериментальные данные.

2.3 Расчет абсорберов и реакторов идеального вытеснения

Ha основании полученных ранее расчетных зависимостей константы скорости и заданной константы равновесия от температуры проведены расчеты реакторов в MathCAD и Microsoft Excel. B результате были получены температуры, расходы и концентрации компонентов на выходе из реакторов.

Таблица 2.7 Характеристики реактора идеального вытеснения Для абсорберов также были рассчитаны температуры на выходе и с учетом степени абсорбции найденной по формуле (1.9), вычислены концентрации компонентов A, B и C.

Таблица 2.8 Характеристики абсорбера

2.4 Синтез системы теплообмена Таблица 2.9 Исходные данные для проектирования тепловой системы Значения водяного эквивалента W рассчитывают по формуле:

W = G C _р

где G — расход смеси, С _р — удельная теплоемкость.

Таблица 2.10 Таблица пар исходных потоков (I этап синтеза) В данной работе заданы две эвристики, следовательно, интервал от 0 до 1 разбивается на два равных. Выбираем наугад псевдослучайное число А, которое указывает интервал, а следовательно, и соответственное эвристическое правило, с помощью которого осуществляется данный этап синтеза. А=0,2881, а значит, на основании 1 эвристики выбраны 2-й холодный и 1-й горячий потоки, так как они имеют наибольшую начальную температуру.

Возможные комбинации теплообмена на втором этапе синтеза приведены в таблице 2.11

Таблица2.11 Таблица пар результирующих и исходных потоков (II этап синтеза)

А=0,9525, следовательно, 2 эвристика, согласно которой выбираются потоки: 3-й холодный и 1-й горячий, так как они обеспечивают максимальное количество передаваемого тепла.

Таблица 2.12 Таблица пар результирующих и исходных потоков (III этап синтеза)

А=0,1802, значит эвристика номер 1. Выбираем 3-й холодный и 2-й горячий потоки так как они имеют наибольшую начальную температуру.

Таблица 2.13

Таблица пар результирующих и исходных потоков (IV этап синтеза) А=0,5672, 2 эвристика. Выбираем 1-й холодный и 3-й горячий потоки, так как они обеспечивают максимальное количество передаваемого тепла.

B результате синтеза тепловой системы все 3 горячих потока достигли заданной конечной температуры, поэтому холодные потоки паром догревают до заданной конечной температуры.

Таким образом, тепловая схема имеет следующий вид:

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

G_п =Q _н / r ,

где Q _н — тепловая нагрузка, r — теплота конденсации насыщенного водяного пара.

G_п1=50,66 028· (698−507,5)/(520·4,1868)= 4,432 577кг/с = 15 957,28 кг/ч Поверхность нагревателя:

F₁= м²

Стоимость нагревателя:

Ц_н1=483· F₁^0,6=483·3691,28 161^0,6=66 719,6006 ус.ден. ед.

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

Для нагрева этого потока берем температуру пара равную 490⁰С, потому что тепло передается от более нагретого тела к менее нагретому. Но, в тоже время, стоимость пара не должна быть очень велика.

G_п2=49,224 138· (753−733)/(520 · 4,1868)= 0,45 217 кг/с = 1627,812 кг/ч Поверхность нагревателя:

F₂= м²

Стоимость нагревателя:

Ц_н2=483· F₂^0,6=483·2115,0752^0,6=47 768,6116 ус.ден. ед.

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

G_п3=45,133 704· (683−639,6)/(520 · 4,1868)= 0,899 672кг/с = 3238,821 кг/ч Поверхность нагревателя:

F₃= м²

Стоимость нагревателя:

Ц_н3=483· F ₃^0.6=483· 403,4 375^0.6=32 322,3385 ус.ден. ед.

Приведенные затраты синтезированной системы:

З_пр = 0,12(Ц ₁+Ц ₂+Ц ₃+Ц ₄+Ц _н1+Ц _н2+Ц _н3)+8800*(G _п1+G_п2+G _п3)*0,005

З_пр=0,12· (64 060,486+87 331,387+19 298,498+57 134,145+66 719,6006+47 768,6116+ 32 322,3385)+8800· 0,005·(15 957,28+1627,812+3238,821) = 961 208,2 ус.ден.ед.

Выводы

B курсовой работе проведена обработка экспериментальных данных методом наименьших квадратов, в результате которой были получены расчетные зависимости константы скорости от температуры К = 712 085 · exp (82 914,21/RT). C помощью метода Брандона осуществили построение модели абсорбера Т_вых = 333,2136· (-0.0065· V+1.1602)·(1,0612· e^-0,0035с)·(0,0004· T_вх+0.8367)

y = 91,56 864· (2,9903·е^{-0,0024Твх})·(0,7408·е^0,0121V)·(0,0089·с+0,8572)

Используя полученные данные, рассчитали входные и выходные расходы, температуры и концентрации в реакторах и абсорберах (таблица 2.6). Ha основании этих расчетов синтезировали теплообменную систему, которая состоит из внутренней и внешней систем.

Список использованных источников

Мухленов И. П. Химико-технологические системы. Синтез, оптимизация и управление.-Л.: Химия, 1986 — 444 с.

Саутин C.H., Иванова E.H., Петрик Ф. Ц., Князьков B.M. Составление экспериментально-статистических моделей по методу Брандона.

Методические указания/ СПбГТИСПб, 1988;24 с.

Мухленов И. П Общая химическая технология. B 2 т. — T. 2, M.; Высшая школа, 1984 г.

Кузичкин H.B., Саутин C.H., Холоднов B.A., Хартманн K. Синтез оптимальных тепловых систем. Методические указания/ СПбГТИСПб₅ 1985;19c.

Приложение 1

Зависимость константы скорости от температуры

Приложение 2

Экспериментальные данные по работе абсорберов