99532

Разработка и исследование математической модели статики процесса ректификации многокомпонентной смеси (объект моделирования – двухколонная установка ректификации дивинила из многокомпонентной смеси)

Курсовая

Экономическая теория и математическое моделирование

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

Русский

2016-09-21

1.21 MB

0 чел.

Разработка и исследование математической модели статики процесса ректификации многокомпонентной смеси (объект моделирования – двухколонная установка ректификации дивинила из многокомпонентной смеси)

Введение

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

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

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

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

Бутан CH3CH2CH2CH3 при нормальных условиях — газ, а под давлением может находиться в жидком состоянии и при обыкновенной температуре.

Бутилен CH2=CHCH2CH3 при нормальных условиях также газообразен.

Циклопентадиен образуется при пиролизе циклопентена:

Введение второй двойной связи в молекулу и её положение мало сказывается на температуре кипения tкип= 42,5оС.

Дивинил (бутадиен)  CH2=CHCH=CH2 – бесцветный газ с характерным неприятным запахом, сжижающийся при -5оС. Является одним из важнейших мономеров  для производства синтетических каучуков таких как: СКБ (бутадиеновые каучуки), СКС (Бутадиен-стирольный каучук), СКН (Бутадиен-нитрильный каучук), а так же: латексов, пластмасс и других органических соединений.

1. Математическое описание объекта, допущения.

Объект моделирования - двухколонная установка ректификации дивинила из многокомпонентной смеси.

     При составлении математического описания процесса принимаются следующие допущения:

- пар, покидающий отдельные тарелки, находится в  равновесии с жидкостью, т.е. рассматриваются теоретические тарелки;

- температурный профиль по высоте колонны принимаем линейным;

- разделительная способность куба колонны эквивалентна  теоретической тарелке;

- весь пар, поступающий с последней тарелки, полностью конденсируется, и часть жидкости поступает на орошение в качестве флегмы;

-  унос жидкости паром с тарелок отсутствует;

- жидкость на тарелках колонны, а также в кубе и дефлегматоре идеально перемешиваются;

- расходы жидкости и пара по высоте отдельной секции колонны постоянны;

- сырье и рецикл поступают в колонну 1 в жидкой фазе;

- коэффициенты относительной летучести постоянны и не зависят от температуры.

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

Для колонны 1 общий материальный баланс в статическом режиме описывается уравнениями:

                                     F1=D1+W1-FR              (1)

                                     V1N1=R1+D1                  (2)

Расходы пара и жидкости в отдельных секциях колонны:

V1i=R1+D1                     0  i N1                   (3)

L1i=F1+R1+FR              1  i NFR                (4)

L1i=F1+R1                     NFR< i NF1            (5)

L1i=R1                           NF1< i N1+1           (6)

L10=W1=F1-D1+FR                                          (7)

Уравнения покомпонентного материального баланса в стационарном состоянии для каждой из тарелок колонны, от куба до дефлегматора:

  Номер тарелки

                                                                       0

                                                      1

     NFR

        NF1    (8)

                                                i

                                N1

                                                             N1+1

   где  - расход пара и жидкости, уходящие с i -той  тарелки колонны 1.

Для колонны 2 общий материальный баланс в статическом режиме описывается уравнениями:

W1=D2+W2+FR             (9)

V2N2=R2+D2+FR           (10)

Расходы пара и жидкости в отдельных секциях колонны:

V2i=R2+D2+FR             0  i N2            (11)

L2i=W1+R2                    1  i NF2           (12)

L2i=R2                           NF2< i N2+1    (13)

L20=W2=W1-D2-FR                                  (14)

Уравнения покомпонентного материального баланса в стационарном состоянии для каждой из тарелок колонны, от куба до дефлегматора:

  Номер тарелки

                                                                        0

                                                      1

                                                i

          NF2 (15)                                                                      

                                N2

                                         N2+1

   где  - расход пара и жидкости, уходящие с i -той  тарелки колонны 2.

При известных потоках   пара  и жидкости  (рассчитанных по формулам (3-7))  уравнения покомпонентного баланса (8) представляют собой систему (N1+2)k уравнений относительно

При этом общее число неизвестных концентраций  на каждой тарелке равно (N1+2)k, а концентраций   - (N1+1)k. Таким образом, общее число неизвестных (2N1+3)k превышает число уравнений. Используя уравнение Генри

,                    (16)

где Yj, Xj - равновесные концентрации j-го компонента в газовой и жидкой фазах;

Kj - константа фазового равновесия;

j - коэффициент относительной летучести j-го компонента,;

T, P - температура и общее давление смеси.

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

Одним из способов решения поставленной задачи является метод независимого определения концентраций (метод Льюиса-Матиссона), согласно которому система уравнений (8) приводится к линейному виду относительно определяемых концентраций  заменой в (8)

,                     (17)

где - константа фазового равновесия на i - ой тарелке j - го компонента в колонне 1.

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

 

                      (18)

                         

В качестве начального распределения принимаем концентрации, равные составу исходной смеси = (, )

Ввиду произвольности начального профиля концентраций получающиеся составы  для каждой ступени после решения системы (18) в сумме не равны единице (или 100%). Поэтому последние корректируются  (с помощью нормирования) и используются в качестве нового приближения для последующей итерации. Нормирование производится по формуле:

      (i=, j=)       (19)

Критерием окончания решения является выполнение стехиометрических соотношений  (i=)  с заданной степенью точности для всех ступеней разделения

          (),                (20)

  где - точность расчета состава.

Матрица коэффициентов системы (18) является трехдиагональной:

                  (21)

Аналогично изложенному выше для колонны 2 систему (15) преобразуем к линейному виду:

 

                       (22)

где K2i,j - константа фазового равновесия на i - ой тарелке j - го компонента во второй колонне, рассчитывается аналогично K1i,j .

В качестве начального профиля концентраций берется состав кубового продукта колонны 1 после первого просчета колонны 1.

Нормирование концентраций  проводится аналогично колонне 1.

Последовательность расчета комплекса колонн следующая:

  1.   рассчитываем профиль концентраций колонны 1,
  2.   рассчитываем профиль концентраций колонны 2,
  3.   проверяем условие окончания вычислений:

               (23)

XFR,1 - концентрация ведущего компонента (дивинила) в потоке рецикла;

X2N2+1,1 - концентрация дивинила в дефлегматоре колонны 2.

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

Проверка правильности расчета осуществляется по общему покомпонентному материальному балансу системы в соответствии с уравнением:

.              (24)


2
Алгоритм решения Схема алгоритма расчёта стационарного режима ректификационных колонн


3. Программирование решения

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

Программа написана в среде MathCAD 2003.

Программа позволяет выполнить:

  •  Ввод исходных данных
  •  Просмотр исходных данных
  •  Расчёт по текущим исходным данным
  •  Вывод результатов расчёта и графиков на экран

Листинг программы приведён в приложении.

4. Результаты исследования моделируемого объекта:

Рис. 1. График концентраций компонентов в жидкой фазе по тарелкам, для колонны 1.

Рис. 2. График концентраций компонентов в жидкой фазе по тарелкам, для колонны 2.

Рис. 3. График концентраций компонентов в паровой фазе по контактным ступеням, для колонны 1.

Рис. 4. График концентраций компонентов в паровой фазе по контактным ступеням, для колонны 2.

5. Анализ результатов исследования моделируемого объекта

Исходная смесь состоит из низкокипящих компонентов имеющих tкип:

- дивинил –                    - 4,5 оС,

- бутилен –                     - 6,3 оС,

- бутан –                         - 0,5 оС,

и высококипящих компонентов:

- циклопентадиен –        42,5 оС;

Анализируя полученные результаты получаем:

  1.  Для первой колонны концентрации компонентов в жидкой фазе изменяются:
    •  циклопентадиен от 2.3% в кубе  до 1.6510-6% в дефлегматоре;
    •  дивинил от 65.1% в кубе до 90.4% в дефлегматоре;
    •  бутилен от 32.6% в кубе  до 32.6%  на 39-ой тарелке, 1% в дефлегматоре;
    •  бутан от 510-3% в кубе  до 10% в дефлегматоре;

Т.е. в процессе протекания ректификации в первой колонне концентрация высококипящего компонента (циклопентадиен) уменьшается при  движении снизу вверх колонн, а концентрация низкокипящих компонентов (дивинил, бутан) увеличивается, что соответствует физическому смыслу. Уменьшение концентрации бутилена в дефлегматоре связано с тем, что сырьё с концентрацией бутилена 56%, подаётся на 39 тарелку первой колонны рис. 1.

  1.  Для второй колонны концентрации компонентов в жидкой фазе изменяются:
  •  циклопентадиен от 99.9% в кубе  до 2.2% в дефлегматоре;
  •  дивинил от 3.610-8% в кубе  до 60% в дефлегматоре;
  •  бутилен от 4.510-3% в кубе  до 32% в дефлегматоре;
  •  бутан от 0% в кубе  до 5.3210-3% в дефлегматоре;

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


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

    1. Для первой колонны концентрации компонентов в жидкой фазе изменяются:

  •  циклопентадиен от 2.29% в кубе  до 1.710-7% в дефлегматоре;
    •  дивинил от 95.386% в кубе до 0.569% в дефлегматоре;
    •  бутилен от 0.389% в кубе  до 81.009% в дефлегматоре;
    •  бутан от 1.94% в кубе  до 18.417% в дефлегматоре;

Уменьшение концентрации низкокипящего компонента дивинила связано с подачей на 39 тарелку первой колонны 44 т/ч дивинила с концентрацией 26%.

    2. Для второй колонны концентрации компонентов в жидкой фазе изменяются:

  •  циклопентадиен от 99.9% в кубе  до 2.21% в дефлегматоре;
  •  дивинил от 6.4710-7% в кубе  до 95.46% в дефлегматоре;
  •  бутилен от 2.69510-13% в кубе  до 0.39% в дефлегматоре;
  •  бутан от 1.4910-14% в кубе  до 1.938% в дефлегматоре;

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

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

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

  1.  Кафаров В. В.  Глебов М. Б. Математическое моделирование основных процессов химических производств. М:  Высшая школа, 1991. С. 261-264, 386 .
  2.  Методические указания    134 по дисциплине «Моделирование объектов и систем управления»  В. С. Кудряшов, А. Н. Гаврилов.

 

 

PAGE  16


 

А также другие работы, которые могут Вас заинтересовать

48702. Гломестный канал (УМК) наземного фазового моноимпульсного радиолокатора (РЛ) дальнего обнаружения объектов с ЭПР 1.16 MB
  Расчёт параметров сигнала. Параметры РЛ: дальность действия от соответствующей длительности импульса зондирующего сигнала до где длительность прямой видимости при высоте цели сектор обзора по углу места при разрешающей способности по углу. При расстоянии до цели погрешность измерения угла места не должна превышать заданного значения при коэффициенте шума приемника равном 3 и потерях энергии сигнала по высокой частоте и при обработке. 2 Определить параметры антенны; зондирующего сигнала; трактов формирования...
48703. Расчет изменения частоты вращения вала 1.32 MB
  Опорами ротора служат подшипники скольжения 8 с жидкой принудительной смазкой (под давлением) от маслоустановки агрегатов. Остаточное осевое усилие ротора воспринимают два упорных подшипника 9. Рабочее колесо литое, одностороннего входа. Направляющий аппарат – литой.
48704. Влияние формы контура области питания скважины. Возможность использования формулы радиального притока в случае нерадиального движения жидкости к скважине 891.5 KB
  Представим себе, что в однородный горизонтальный пласт весьма больших (теоретически неограниченных) размеров и постоянной мощности проведены гидродинамически совершенные равнодебитные нагнетательная и эксплуатационная скважины одинакового радиуса R.
48705. Анализ технических возможностей способов сварки плавлением барабана изготовленного из стали 10 216 KB
  Введение Сварка широко применяется в основных отраслях производства так как резко сокращает сроки выполнения работ и трудоемкость производственных процессов. Сварка позволяет уменьшить затраты на единицу продукции сократить длительность производственного цикла улучшить качество изделий. Для данного изделия возможны следующие способы сварки плавлением: ручная дуговая сварка; сварка в защитных газах плавящимся электродом; –плазменная сварка; лазерная сварка; электроннолучевая сварка; газовая сварка. ручная дуговая сварка покрытым...
48706. Сохранение и укрепление здоровья обслуживаемого населения 903.5 KB
  Номер – тип: int Тип стрипа – тип: int 8луночные и 12луночные Режим измерения – тип: chr измерение оптической плотности Единицы измерения – тип: chr бел Выходные значения анализатора Объект представляет собой строку байт. Тип – тип: byte Преобразованные значения Объект представляет собой данные строкового формата. Тип – тип: chr Документы предметной области также представлены в виде классов. Номер – тип: int Дата – тип: dte ФИО – тип: chr Возраст – тип: int Учреждение – тип: chr Отделение – тип: chr Исследовать – тип: chr Диагноз...
48707. Генеалогическое древо. Информационная система 3.95 MB
  Прежде чем приступать к разработке информационной системы, необходимо представить себе схему настоящего генеалогического древа. Генеалогическое древо - схематичное представление родственных связей, родословной росписи в виде условно-символического «дерева»
48708. Стратегическое планирование на ОАО «МТС» 1.13 MB
  О предприятии ОАО МТС. В курсовой работе в качестве исследуемого предприятия была выбрана компания ОАО МТС в качестве продукта – мобильный телефон. ОАО Мобильные ТелеСистемы МТС российская телекоммуникационная компания оператор сотовой связи в форматах GSM и UMTS оказывающая услуги в России странах СНГ и Индии под торговой маркой МТС. Компания МТС образована как закрытое акционерное общество в октябре 1993 года такими компаниями как ОАО Московская городская телефонная сеть МГТС Deutsсhe Telecom DeTeMobil...