|
Реферат: Модель теплового состояния аппарата сепарации (Технология)
Модель теплового состояния аппарата сепарации
Ставится задача определения времени, необходимого для окончания
процесса сепарации аппарата восстановления титана, и теплового состояния
сепарируемой массы во время процесса.
Нагрев аппарата происходит в три стадии:
Прогрев реакционной массы. Оканчивается, когда на внутренней поверхности
стенки аппарата достигается температура кипения магния, соответствующая
поддерживаемому в аппарате давлению.
Кипение летучих. Будем полагать, что фронты кипения Mg и MgCl2 движутся
поступательно внутрь аппарата от стенки, образуя коаксиальные
цилиндрические поверхности (см. рис. ниже).
Стадия конечного прогрева после выкипания большей части летучих. Тепловые
свойства аппарата определяются свойствами титановой губки.
Аппарат находится в печи сепарации. Тепло к нему подводится вследствие
теплообмена излучением с нагревателями печи и конвективного теплообмена с
воздухом, заполняющим печь. В первом приближении, суммарный тепловой поток
[pic], воспринимаемый аппаратом, можно выразить как
[pic],
где [pic]– температура наружной поверхности стенки аппарата,
[pic] – температура нагревателей,
[pic] – температура воздуха в печи,
[pic] – коэффициент теплоотдачи при свободной конвекции,
[pic] – интегральный коэффициент теплопередачи излучением, зависящий от
степени черноты тел и углового коэффициента облучения.
Для системы печь – аппарат можно принять коэффициент теплопередачи
излучением
[pic],
где [pic] Вт/(м2К4),
[pic],
[pic] м2 – площадь излучателей (нагревателей печи),
[pic] м2 – площадь поверхности аппарата,
[pic]м – радиус аппарата,
[pic] м – высота аппарата.
Коэффициент лучистого теплообмена для системы воздух – аппарат
определяется, исходя из парциальных давлений паров воды и углекислого газа.
В рассматриваемых условиях наличие паров воды маловероятно, а коэффициент
теплового излучения CO2 в интервале температур 600...1200 оС близок к 0.3,
то есть
[pic],
где [pic]– коэффициент излучения углекислого газа заданной толщины (0,25 м)
и парциального давления (0.005 атм),
[pic] – то же для условно бесконечного слоя.
Аналогичным образом можно записать тепловые потоки для нагревателей
печи и воздуха, заполняющего печь.
Расход тепла из аппарата происходит излучением через крышку в реторту-
конденсатор:
[pic]
где [pic] – коэффициент теплопередачи от аппарата к реторте-конденсатору,
[pic] – температура реторты-конденсатора (( 390 К).
Согласно записанным выше соображениям, порядок [pic] можно оценить
следующим образом:
[pic].
Множитель 1/2 принят из-за того, что между аппаратом и конденсатором
находится тепловой экран, как минимум вдвое снижающий лучистый тепловой
поток.
Кроме того, происходит унос тепла вместе с продуктами возгонки.
Оценить его можно, только достоверно зная массовый поток и температуру
сублимированных продуктов. Этот вопрос выходит за рамки настоящего
исследования.
На первой стадии можно рассматривать аппарат как сплошное
цилиндрическое тело. Задача нагрева бесконечного цилиндра, помещенного в
подогревающую среду, имеет аналитическое решение
[pic],
где [pic] – относительная температура,
[pic] – температура цилиндра на радиусе [pic],
[pic] – начальная температура цилиндра (до нагрева),
[pic] – температура подогревающей среды (воздуха в печи),
[pic] – n-й корень характеристического уравнения [pic],
[pic] – критерий Био,
[pic] – коэффициент теплоотдачи от подогревающей среды,
[pic] – радиус аппарата, [pic] м,
[pic] – коэффициент теплопроводоности материала цилиндра,
[pic], [pic] – коэффициенты,
[pic] – критерий Фурье,
[pic] – коэффициент температуропроводности материала цилиндра,
[pic] – плотность материала цилиндра,
[pic] – теплоемкость материала цилиндра,
[pic] – время прогрева,
[pic] – функция Бесселя k-го порядка, являющаяся решением уравнения
[pic];
[pic].
Оценим порядок критериев, входящих в это уравнение.
По окончании процесса восстановления в аппарате содержится порядка 60%
(массовых) Ti, 20...30% Mg и 10...20% MgCl2. Плотность титана 4.35, магния
1.8, MgCl2 порядка 2.7, следовательно, средняя плотность реакционной массы
[pic]кг/м3.
Примем следующие зависимости от температуры теплофизических
параметров:
[pic] Дж/(кг К) в твердом состоянии,
[pic] Дж/(кг К) в жидком состоянии,
[pic] Дж/(кг К) в твердом состоянии
[pic] Дж/(кг К) в жидком состоянии,
[pic] Дж/(кг К),
[pic] Вт/(м К) в твердом состоянии
[pic] Вт/(м К) в жидком состоянии
[pic] Вт/(м К)
[pic] Вт/(м К)
[pic] кДж/кг
[pic] кДж/кг
Таким образом, при температуре 800 К, средней в рассматриваемом
диапазоне, средняя теплоемкость аппарата
[pic]Дж/м3.
и средняя теплота парообразования
[pic] кДж/моль.
Теплопроводность титановой губки можно оценить по соотношению
[pic],
где [pic] – теплопроводность титана, [pic] Вт/(м2К),
[pic]– пористость титановой губки (отношение объема пор к общему объему).
|Пористость губки|Теплопроводность, |
| |Вт/(м2К), |
|0.2 |14.0 |
|0.3 |11.0 |
|0.4 |8.0 |
При средней пористости блока 0.2...0.3 можно принять [pic].
Порядок средней теплопроводности аппарата в целом можно оценить,
исходя из массового состава блока:
[pic] Вт/(м2К).
Тогда коэффициент температуропроводности изменяется в пределах [pic],
а число Фурье [pic].
Так как время прогрева порядка нескольких часов или даже суток, то
величина [pic]. При больших значениях [pic] достаточно одного члена ряда:
[pic].
Оценим порядок критерия Био. Коэффициент теплопередачи при свободной
конвекции (cм. ниже) можно принять порядка 8...10. Тогда для аппарата
диаметром 1.5 метра среднее значение [pic]. При малых [pic] можно считать,
что температура прогреваемого цилиндра по всему сечению близка к
одинаковой, то есть равна температуре на поверхности аппарата.
Величина [pic] может быть найдена из соответствующих критериальных
зависимостей, например, для ламинарного режима
[pic],
где [pic] – критерий Нуссельта,
[pic]– теплопроводность воздуха в печи,
[pic] – критерий Прандтля,
[pic] – кинематическая вязкость воздуха, зависящая от температуры,
[pic] – температуропроводность воздуха, зависящая от температуры,
[pic] – критерий Грасгофа,
[pic] – ускорение свободного падения,
[pic],
[pic] – средняя температура воздуха в печи.
Для оценочных расчетов коэффициента теплоотдачи при свободной
конвекции на вертикальной поверхности иногда используется формула [pic].
Оценим порядок: пусть [pic] и [pic]. Примем температуру воздуха в печи
порядка средней между температурой нагревателей и аппарата: [pic]. тогда
[pic]. В выкладках, приведенных ниже, используется именно эта зависимость,
хотя, поскольку интегрирование приведенной ниже системы обыкновенных
дифференциальных уравнений возможно только численными методами, можно
применять и более точные методы определения [pic]. Можно также определить
[pic] экспериментальным путем и затем использовать эмпирическую
зависимость.
Зависимость теплоемкости вещества от температуры, как правило, хорошо
приближается полиномами 1...2 степеней: [pic].
Запишем уравнения теплового баланса для нагревателей печи, воздуха и
аппарата:
[pic]
[pic]
[pic]
где[pic] – средняя теплоемкость нагревателей, воздуха или аппарата,
зависящая от температуры реакционной смеси,
[pic] – масса нагревателей (1), воздуха(f) или аппарата(o),
[pic] – мощность, подводимая к печи,
[pic] – тепловой КПД печи, определяемый экспериментально.
Полученная система дифференциальных уравнений легко разрешается любым
численным методом. Она достаточно устойчива, что позволяет поддерживать шаг
интегрирования явным методом порядка 100 секунд.
[pic]
На графике показана зависимость температуры аппарата от времени (при
условии отсутствия кипения), полученная решением системы (О) при
эффективной подводимой мощности [pic] кВт, и по эмпирическому уравнению
[pic],
где [pic] – заданная температура, [pic] = 0.4, аппроксимирующему
экспериментальные данные.
Зависимость давления паров Mg и MgCl2 от температуры хорошо
приближается формулой [pic], где [pic] – эмпирические коэффициенты, или
[pic]. Для магния [pic], [pic]; для MgCl2 [pic], [pic]. Задаваясь давлением
в аппарате, находим температуру кипения. Время прогрева получим, интегрируя
приведенную выше систему до момента, когда температура сравняется с
температурой кипения магния при заданном давлении в аппарате.
Максимальная температура, до которой можно прогревать аппарат,
определяется исходя из скорости взаимодействия губки с материалом реторты и
лежит в пределах [pic]К. Имеет смысл во время прогрева препятствовать
началу кипения, поддерживая в аппарате некоторое давление. Так, в
указанном интервале температур магний вскипает при давлении ниже
0.026...0.05 МПа. Затем производится сброс давления до давления сепарации
(10...100 Па), в результате чего начинается интенсивное объемное кипение по
всему объему аппарата. При указанном выше содержании летучих компонентов
можно ожидать, что силы поверхностного натяжения будут препятствовать
появлению гидростатической составляющей давления в объеме аппарата. Можно
оценить количество летучих, которые выкипят за время сброса давления,
исходя из баланса тепла:
[pic],
где [pic] – удельная теплота парообразования. Окончательно имеем
[pic].
Это уравнение имеет аналитическое решение
[pic],
где [pic] – начальная масса летучих, или при постоянной теплоемкости
аппарата
[pic].
где [pic] – изменение температуры, [pic]. По окончании объемного кипения
вся реакционная масса будет иметь одинаковую температуру [pic].
Расчет по программе, реализующей приведенную выше методику, дает
следующую зависимость массы летучих в аппарате от количества циклов нагрев
– сброс давления:
[pic]
Таким образом, за 7 циклов теоретически возможно полностью удалить все
летучие из аппарата.
Неясным является следующий вопрос: как распределяется тепло,
аккумулированное аппаратом в процессе нагрева под давлением, между
летучими? В приведенном расчете предполагалось, что количество поглощаемой
на испарение теплоты пропорционально теплопроводности и количеству магния
или хлористого магния в аппарате.
[pic]
Рассмотрим стадию сепарации, на которой происходит кипение остатков
летучих. Будем полагать, что при достижении реакционной массой температуры
кипения, соответствующей поддерживаемому давлению, вся подводимая теплота
уходит на испарение. Тогда баланс тепла на фронте кипения позволяет найти
скорость его движения. Очевидно, что
[pic],
где [pic] – тепловой поток, подводимый к фронту,
[pic] – тепло, расходуемое на кипение,
[pic] – тепловой поток, пропускаемый фронтом.
Тепловой поток [pic] считаем полностью расходуемым на испарение:
[pic]
или
[pic].
Будем предполагать квазистационарное распределение температур в
реакционной массе. Это справедливо, если прогрев происходит достаточно
медленно ([pic]), что доказано выше. В этом случае зависимость температуры
от радиуса цилиндра при заданных граничных условиях 1 рода (температура на
внешней и внутренней поверхности) выражается формулой
[pic],
где [pic] – текущий радиус,
[pic] – температура наружной (горячей) поверхности,
[pic] – температура внутренней (холодной) поверхности,
[pic] – радиус наружной (горячей) поверхности,
[pic] – радиус внутренней (холодной) поверхности,
а удельный линейный тепловой поток в цилиндр формулой [pic] Вт/м.
В принятой модели не рассматривается теплообмен блока с внешней
средой, кроме подвода теплоты от печи, поэтому вся теплота [pic] (см. рис.)
расходуется на испарение магния, а теплота [pic], расходуемая на испарение
MgCl2, находится из приведенных выше зависимостей.
Окончательно, получим систему уравнений для определения скоростей
движения фронтов кипения:
[pic];
[pic];
[pic];
[pic];
[pic];
[pic];
[pic];
[pic].
Интегрируя эту систему до тех пор, пока [pic], находим время,
необходимое для испарения всех летучих.
Возможен и другой подход к расчету второй стадии. Можно предположить,
что после интенсивного объемного кипения летучие не будут заполнять
протяженные сплошные области, а распределятся по блоку “каплями”, не имея
друг с другом хорошего теплового контакта. Тогда теплопроводность блока
будет определяться исключительно теплопроводностью губки, и
[pic]
и, соответственно, [pic]; [pic]; [pic]
Следовательно, нельзя принимать допущение, что температура блока по
сечению постоянна. Используя формулу
[pic],
справедливую при больших [pic], можно получить зависимость температуры на
оси цилиндра от времени. Температура подогревающего воздуха, в общем
случае, неизвестна, но принимая во внимание большую продолжительность
процесса и малое изменение температуры поверхности аппарата на 2 стадии,
можно предположить, что она приблизительно равна температуре нагревателей
печи. Недостаток такого предположения частично покрывается тем, что в
дальнейшем выводе опущен теплоподвод излучением от нагревателей печи и
подогревающего воздуха.
На внешней поверхности цилиндра ([pic]) имеем [pic]
или
[pic].
Отсюда можно получить время прогрева поверхности до величины [pic] как
[pic]
или
[pic];
[pic];
[pic];
[pic].
При [pic] имеем [pic] и
После этого, по условиям процесса, температура у стенки аппарата
поддерживается постоянной и равной [pic]. Изменение температуры на оси
цилиндра при граничных условиях 1 рода определяется формулой
[pic],
или
[pic];
[pic];
[pic];
[pic];
[pic].
Таким образом, получены все необходимые и достаточные аналитические
выражения для определения продолжительности процесса сепарации на 1 стадии
и максимальной продолжительности 2 стадии при различных допущениях о
характере процесса. Продолжительность 3 стадии (высокотемпературной
выдержки) определяется по условию малости остатка летучих в объеме блока.
Она зависит от объемного коэффициента теплопередачи от блока к каплям
летучих и скорости их кипения (оба параметра неизвестны, и для них
отсутствуют эмпирические зависимости), поэтому должна определяться
экспериментально.
Реферат на тему: Модернизация cсудовой энергетической установки грузового судна
1. Анализ показателей судна и его энергетической установки
Для выполнения последующих обоснований и расчётов необходимо привести
краткую характеристику указанного в задании проекта судна (табл. 1.1) и его
энергетической установки (табл. 1.2).
Таблица 1.1
Характеристика судна проекта 559Б
|Параметры, единицы измерения |Значения |
|Класс |“О” |
|Размерения корпуса, м: | |
|длина |80 |
|ширина |15 |
|Водоизмещение, т |1632 |
|Грузоподъёмность, т |1200 |
|Мощность, кВт |588 |
|Осадка, м |1,7 |
|Скорость в полном грузу, км/ч |16,6 |
|Автономность, сут. |6 |
|Число мест для экипажа |12 |
|Количество движителей |2 |
|Тип движителя |ГВН |
|Сухая масса СЭУ, т |35 |
|Габариты машинного отделения, | |
|м: |18 |
|длина |14 |
|ширина | |
Таблица 1.2
Характеристика основных элементов ЭУ судна проекта 559Б
|Элементы ЭУ и их параметры, единицы измерения |Значения |
|Главные двигатели: | |
|количество |2 |
|марка |8NVD36 |
|номинальная эффективная мощность, кВт |294 |
|номинальная частота вращения коленчатого вала, |500 |
|мин-1 | |
|род топлива |дизельное |
|удельный эффективный расход, кг/(кВт(ч): | |
| топлива |0,228 |
| масла |0,00163 |
|Главная передача мощности: | |
|тип |прямая |
|Судовая электростанция: | |
|количество дизель-генераторов |2 |
|марка дизель-генераторов (дизеля) |ДГ-50/1 |
| |(6Ч12/14) |
Продолжение табл. 1.2
|Элементы ЭУ и их параметры, единицы измерения |Значения |
|номинальная эффективная мощность (дизеля), кВт |50 (58,8) |
|номинальная частота вращения коленчатого вала, |1500 |
|мин-1 | |
|удельный эффективный расход топлива, кг/(кВт(ч) |0,269 |
|КОТЕЛЬНАЯ УСТАНОВКА: | |
|марка автономного котла |КОАВ-68 |
|количество |1 |
|теплопроизводительность, кДж/ч |285600 |
|расход топлива, кг/ч |8,2 |
На основе данных табл. 1.1 и 1.2 выполняется расчёт показателей
установки судна.
Эффективная мощность главной ЭУ:
[pic] кВт,
где [pic] – количество главных двигателей СЭУ;
[pic] кВт – номинальная эффективная мощность главного двигателя СЭУ.
Энергооснащённость судна:
[pic] кВт/т,
где [pic] т – водоизмещение судна в полном грузу.
Энергонасыщенность по отношению к длине машинного отделения:
[pic] кВт/м,
где [pic] – количество вспомогательных двигателей;
[pic] кВт – номинальная эффективная мощность вспомогательного
двигателя;
[pic] м – длина машинного отделения.
Энергонасыщенность по отношению к площади машинного отделения:
[pic] кВт/м2,
где [pic] м2 – площадь машинного отделения.
Энергоёмкость работы судна:
[pic] кДж/(т(км),
где [pic] т – грузоподъёмность судна;
[pic] км/ч – скорость судна в полном грузу.
Удельная масса ЭУ:
[pic] кг/кВт,
где [pic] т – сухая масса ЭУ.
Абсолютный КПД установки:
[pic],
где [pic] кг/ч – общий расход топлива на СЭУ;
[pic] кДж/кг – низшая удельная теплота сгорания дизельного топлива.
Эффективный КПД установки:
[pic]
[pic],
где [pic] кг/(кВт(ч) – удельный эффективный расход топлива главного
двигателя;
[pic] кг/(кВт(ч) – удельный эффективный расход топлива
вспомогательного двигателя;
[pic] – количество автономных котлов;
[pic] – КПД вспомогательного автономного котла;
[pic] кДж/ч – теплопроизводительность автономного котла.
КПД судового (пропульсивного) комплекса:
[pic],
где [pic] – эффективный КПД главного двигателя:
[pic];
[pic] – КПД главной передачи установки;
[pic] – КПД валопровода;
[pic] – пропульсивный КПД.
КПД энергетического комплекса:
[pic]
[pic]
[pic],
где [pic] – эффективный КПД вспомогательного двигателя:
[pic];
[pic] – КПД электрогенератора;
[pic] кг/ч – расход топлива автономным котлом.
Результаты расчёта представлены в табл. 1.3.
Таблица 1.3
Показатели энергетической установки судна проекта 559Б
|Наименование показателей, единицы |Значения |
|измерения | |
|Эффективная мощность главной ЭУ, кВт |588 |
|Энергооснащённость, кВт/т |0,36 |
|Энергонасыщенность по отношению к: | |
|длине МО, кВт/м |39,2 |
|площади МО, кВт/м2 |2,8 |
|Энергоёмкость работы судна, кДж/(т(км) |106 |
Продолжение табл. 1.3
|Наименование показателей, единицы |Значения |
|измерения | |
|Удельная масса ЭУ, кг/кВт |59,5 |
|Эффективный КПД установки |0,286 |
|Абсолютный КПД установки |0,286 |
|КПД судового комплекса |0,214 |
|КПД энергетического комплекса |0,362 |
| |
