Файл: Лебедев Д.П. Тепло- и массообмен в процессах сублимации в вакууме.pdf

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

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

и коэффициента теплоотдачи замораживания льда в по­ ристой пластине

(рис. 5-2).

Из формул (5-8) и (5-9) следует, что коэффициент теплоотдачи при сублимации сильно зависит от тепловой нагрузки, в то время как коэффициент теплоотдачи за­ мораживания льда главным образом зависит от вакуума в сублиматоре.

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

Принятая обработка экспериментальных данных мо­ жет использоваться только для толстых крупнопористых пластин, в которых зоны льда, жидкости и пара экспе-

1 Практическое применение полученных зависимостей ограничи­ вается тепловыми нагрузками от 1 000 до 5 000 ккал/м2 (критической нагрузкой промерзания и критической нагрузкой выброса).

221

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

5-5. МАТЕМ АТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА СУБЛИМАЦИИ В ПРОНИЦАЕМОЙ ПЛАСТИНЕ

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

Детальный анализ механизма процесса сублимации в одиночном цилиндрическом капилляре, «пленочной мо­ дели» с сеткой позволяет при соответствующих допуще­ ниях записать уравнения переноса тепла и массы для всего пористого тела. Как показали исследования изме­

222

условиях размеры зон и температурные поля по тол­ щине пластины во времени непрерывно меняются, оста­ ваясь в то же время постоянными относительно опреде­ ленного среднего значения. Это среднее значение опре­ деляется режимными параметрами: давлением воды рв, интенсивностью тепловой нагрузки q и в меньшей сте­ пени вакуумом в сублиматоре рк (при рк< 67 н/мг), а также структурными и теплофизическими характери-

Рис. 5-21. Физическаямодель процесса сублимации.

а — зонная модель для капиллярно-пористого тела; б — зонная модель для мелкопористых пластин.

стиками пластины: размером пор, пористостью, теплопроводностью и т. д. Изменение размеров зон и темпе­ ратурных полей имеет пульсирующий характер. Цикл

и опускания поверхности сублимации Дтг. В [Л. 5-6] сделана общая постановка рассматриваемой задачи. При выводе уравнений, описывающих рассматриваемый процесс, делаются следующие допущения:

1.Процесс теплообмена является одномерным.

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

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

3. На одной и той же глубине основной материал (пористая пластина) и наполнитель (жидкость, лед, пар) имеют одинаковую температуру [Л. 5-7].

223

4. Теплопроводность паров воды пренебрежимо мала по сравнению с теплопроводностью основного каркаса.

5. Движение границ льда происходит равномерно по всей пластине, т. е. предполагается равномерность пористости и одинаковые условия подвода тепла на нижней поверхности для всей пластины.

6 . Конвективный и лучистый теплообмен на наруж­ ной поверхности и внутри пластины пренебрежимо мал.

7. Теплофизические свойства металлокерамики: су­ хой, смоченной и замороженной — полагаются известны­ ми и постоянными.

8 . Движение жидкости в пористой пластине в зоне

с силами внешнего давления. Предложенная физическая модель процесса подробно разобрана в [Л. 5-8].

Рассмотрим безразмерные системы уравнений, опи­ сывающие поле температур в пористой пластине, для периода микропрорыва и _периода_ промерзания влаги.

Период микропрорыва (ті^ т^ тг):

(5-10)

F O2 dz

224

fi(x)\ fz(x); f3(x) получаются из решения системы урав­ нений (5-10) при граничных условиях (5-11) —(5-17).

В результате решения системы уравнений (5-18) мо­ гут быть получены зависимости:

прежних уравнений и граничных условий, описывающих процесс, и т. д.

К' :/С7;