автореферат диссертации по энергетике, 05.14.04, диссертация на тему:Повышение эффективности выпаривания отработанных варочных растворов целлюлозного производства

На правахржсописи

СУСЛОВ Вячеслав Александрович

«

ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ВЫПАРИВАНИЯ ОТРАБОТАННЫХ ВАРОЧНЫХ РАСТВОРОВ ЦЕЛЛЮЛОЗНОГО ПРОИЗВОДСТВА

Специальность: 05.14.04 - промышленная теплоэнергетика

Автореферат диссертации на соискание ученой степени доктора технических наук

1

САНКТ-ПЕТЕРБУРГ - 2006

Работа выполнена в ГОУ ВПО «Санкт-Петербургский государственный технологический университет растительных полимеров» на кафедре промышленной теплоэнергетики.

Официальные оппоненты:

— доктор техн. наук Готовский Михаил Абрамович

— доктор техн. наук, профессор Зысин Леонид Владимирович

— доктор техн. наук, профессор Иванов Олег Петрович

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

Защита диссертации состоится 28 февраля 2006 г. в 16-00 на заседании диссертационного совета Д 212.229.04 в ГОУ ВПО "Санкт-Петербургский государственный политехнический университет" по адресу: 195251, Санкт-Петербург, ул. Политехническая, 29 в аудитории 411 Ш К

С диссертацией можно ознакомиться в фундаментальной библиотеке ГОУ ВПО "Санкт-Петербургский государственный политехнический университет".

Автореферат разослан "25" января 2006 г.

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

Факс: (812)5521630 E-mail: kgl210@mail.ru

Ученый секретарь диссертационного совета

к А- Григорьев

А/96

3

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

В настоящее время в основном применяются сульфатный и сульфитный способы получения целлюлозы. Доля сульфатной целлюлозы от общего объема ее производства составляет около 90 %. В соответствии с этим в настоящей работе рассматриваются вопросы тепловых и гидродинамических процессов при выпаривании черного сульфатного щелока и водного пред-гидролизата на выпарных станциях сульфат-целлюлозных заводов.

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

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

библиотека |

СПтМрг А П ,

•ЗЬйР.1

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

Работа выполнялась в рамках региональной программы фундаментальных и прикладных исследований по проблеме комплексного использования и воспроизводства лесных ресурсов на 1965 - 1990 гг. АН СССР по теме «Исследование процессов и разработка новых видов оборудования для производства бумаги и картона», в рамках общесоюзной научно-технической программы Минлеспрома СССР по созданию и внедрению новой техники и технологии и государственной научно-технической программы "Комплексное использование и воспроизводство древесного сырья" 1991-1996 гг.

В работе приведены результаты исследования промышленных ВС, из 4

анализа которых получены данные по теплообмену при выпаривании сульфатных щелоков. Опыт, накопленный в период промышленных исследований, позволил спроектировать и осуществить монтаж стендовых и пилотных *> установок с геометрическими размерами кипятильных труб, равными размерам труб промышленных выпарных аппаратов (ВА). На стендах проведены экспериментальные исследования различных сторон выпаривания растворов в условиях, близких к условиям работы промышленных ВА во всем диапазоне режимов их работы. Приведенные зависимости позволяют определить значение любой составляющей общего термического сопротивления теплопередаче в ВА и оценить эффективность мероприятий по интенсификации теплообмена. В работе приводятся результаты исследований в области гидродинамики и теплообмена двухфазных потоков применительно к условиям работы ВА ЦБК. Исследования гидродинамики двухфазных потоков выполнялись для определения ее влияния на интенсивность теплообмена.

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

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

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

Научная новизна результатов.

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

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

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

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

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

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

Осуществлен теоретический анализ эффективности схем ВС и влияния режимных параметров процесса на интенсивность выпаривания.

На основе выполненных исследований разработаны конструкции перспективных выпарных аппаратов.

Создан программный метод теплового расчета ВС.

Основные методические положения работы. Разработка теоретических основ процесса выпаривания отработанных варочных растворов ЦБК базировалась на фундаментальных уравнениях тепло-и массопереноса с использованием теплофизических свойств пенообразующих растворов. Промышленные исследования теплообмена проведены на действующем оборудовании в производственных условиях на 14 ВС различных ЦБК. Позонное исследование теплообмена и гидродинамики при восходящем потоке паро-растворной двухфазной смеси проводилось на стенде с паровым обогревом кипятильной трубы с размерами, равными размерам труб промышленных ВА. Теплоотдача при гравитационном течении испаряющегося раствора исследовалась на стенде с электрообогреваемой трубой. Гидродинамика гравитационно стекающего раствора изучалась на стенде с трубой из кварцевого стекла. Динамика накипеобразования и диффузия пара при конденсации парогазовой смеси изучалась на работающих в промышленных условиях ВС. На пилотном ВА проверялись и отрабатывались эффективность внедряемых разработок: работоспособность распределительных устройств; интенсивность накипеобразования при выпаривании; оптимальные режимы выпаривания. Теоретические исследования проводились с использованием рабочих гипотез, а также с применением методов математического и цифрового моделирования.

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

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

ГУП «ЦНИИБуммаш» и ЗАО «ГИПРОБУМ». Полученные материалы в виде номограмм и расчетных зависимостей с учетом динамики накипеобразования используются как основа для квалифицированного управления работой ВС ЦБК. Результаты работы использовались при оптимизации режимов выпаривания и модернизации ВС Долинского, Архангельского и Котласского ЦБК. Основные положения научных исследований применяются в учебном процессе: при чтении лекций; при курсовом и дипломном проектировании; приведены в четырех учебных пособиях, два из них с грифом Министерства образования РФ.

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались на ежегодных научно-технических конференциях профессорско-преподавательского состава СПбГТУРП (ранее ЛТИ ЦБП), на Всесоюзной конференции «Теплофизика и гидрогазодинамика процессов кипения и конденсации» (Рига, 1985 г.), на VIII Всесоюзной конференции «Двухфазный поток в энергетических машинах и аппаратах» (Ленинград, 1990 г.), на Первой Всесоюзной национальной конференции по теплообмену «Кипение, кризисы кипения, закризисный теплообмен» (Москва, 1994 г.), на Российской научно-практической конференции «Организация природоохранной деятельности на предприятиях и пути ресурсосбережения» (С.-Петербург, 1996 г.), на научно-практической конференции «Опыт и перспективы развития ресурсосберегающих технологий и утилизации твердых отходов на предприятиях различных форм собственности» (С.-Петербург, 1998 г.), на Второй Российской национальной конференции по теплообмену «Кипение, кризисы кипения, закризисный теплообмен. Испарение, конденсация» (Москва, 1998 г.), на Третьей Российской национальной конференции по теплообмену «Кипение, кризисы кипения, закризисный теплообмен. Испарение, конденсация» (Москва, 2002 г.), на Международной научно-технической конференции «Современная наука и образование в решении проблем экономики европейского севера» (Архангельск, 2004 г.), на Международной научно-технической конференции «Ресурсо- и энергосбережение в целлюлозно-бумажной промышленности и городском коммунальном хозяйстве» (С.-Петербург, 2005 г.).

Публикации. По теме диссертации опубликовано 49 работ, в числе которых обзорная информация «Выпарные аппараты и схемы выпарных станций ЦБП», 4 учебных пособия, получено 3 авторских свидетельства.

Структура и объем диссертации. Диссертация состоит из введения, семи глав, выводов, списка литературы - 394 наименования, 137 рисунков, 20 таблиц. Диссертация включает 250 с. основного текста и приложений.

Автор защищает:

1 .Физическую модель движения двухфазного восходящего потока, образованного в результате кипения отработанных варочных пенообразующих растворов в длинных трубах ВА ЦБК.

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

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

4.Результаты теоретико-экспериментальных исследований теплоотдачи при конденсации парогазовой смеси на вертикальном пучке труб.

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

6. Рекомендации по управлению режимными параметрами с целью поддержания оптимальных режимов выпаривания, определяющие основу автоматизированных систем управления ВС ЦБП.

7.Модернизированную противоточную схему ВС как наиболее опти- « мальную для процесса выпаривания отработанных варочных растворов ЦБК.

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

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе приведены типовые варианты и основные конструктивно-технические характеристики ВС для выпаривания отработанных варочных растворов ЦБК. Дан критический анализ состояния теоретических и экспериментальных исследований по теплообмену при кипении и испарении жидкости в условиях восходящего и гравитационного ее течения в вертикальных трубах. Имеется весьма ограниченное количество исследований, посвященных данной проблеме, вследствие отсутствия достаточных и проверенных данных по теплообмену при выпаривании отработанных варочных растворов ЦБК. Большинство работ по исследованию теплообмена при кипении жидкостей выполнено применительно к чистой воде. Часть работ проведена на растворах, главным образом сахарных. В обоих случаях в основу обработки данных положены закономерности кипения чистых жидкостей. Полученные расчетные формулы для определения интенсивности кипения жидкостей имеют существенные рассогласования. Данные по теплообмену при конденсации водяного пара из парогазовых смесей на вертикальных пучках труб практически отсутствуют. На основе опубликованных данных вывести критерии моделирования процессов, протекающих в ВА, не представляется возможным. Поэтому наиболее достоверные экспериментальные данные в настоящее время могут быть получены только при исследовании процессов

выпаривания реальных отработанных варочных растворов ЦБП в трубах с размерами, аналогичными размерам труб промышленных ВА.

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

Вторая глава посвящена вопросам промышленных исследований режимов работы ВС ЦБК.

Представлена разработанная методика проведения теплотехнических исследований режимов работы ВС, в соответствии с которой изучалась теплопередача от насыщенного пара к кипящим растворам на 14 ВС различных ЦБК. ВС укомплектованы длиннотрубными вертикальными однопроходными аппаратами с восходящим движением выпариваемого раствора, длиной кипятильных труб 8+9 м, поверхностью теплообмена до 5500 м2. В результате исследований обнаружено, что зависимость интенсивности передачи теплоты от концентрации выпариваемых щелоков обратно пропорциональна. Значения коэффициентов теплопередачи в ВА увеличиваются с увеличением плотности теплового потока. ВС с небольшой поверхностью теплообмена, работающие по противоточным схемам питания, имели наилучший показатель тепловой эффективности.

В третьей главе приведено описание конструкций опытных стендов и методики эксперимента, полностью воспроизводящих процессы теплообмена и гидродинамики при выпаривании растворов ЦБК. Обоснованно приняты математические модели процессов теплообмена и в соответствии с этими моделями спланирован эксперимент. Приведены критерии для оценки соответствия модели результатам исследований. В соответствии с поставленной задачей исследования теплообмена при выпаривании отработанных растворов ЦБП были разработаны методики проведения исследований и обработки полученных опытных данных. Для проведения намеченных экспериментальных исследований были спроектированы и изготовлены:

- стенд с паровым обогревом для изучения теплопередачи при подъемном движении потока; основной рабочей частью стенда являлась кипятильная труба диаметром (1 = 51 х 1,5 мм и длиной Ь = 8,5 м, Ь/с1 =165; кипятильная труба разделена конденсатосборными стаканами на 19 измерительных

участков; гидродинамические и тепловые параметры стенда позволяли охватить весь диапазон работы промышленных установок;

- установка для исследования гидродинамики гравитационно стекающей пленки жидкости внутри кварцевой трубы, длиной 2 м, Ь/с1=95;

- установка для исследования теплоотдачи к стекающей пленке жидкости внутри стальной трубы с электрообогревом, длиной 2 м, Ь/с1=95;

- установка для исследования накипи и борьбы с ней, состоящая из однотрубного вертикального однопроходного ВА с принудительной циркуляцией раствора и паровым обогревом трубы длиной Ь = 8,5 м, ЫА = 163, аппарата магнитной обработки раствора, подогревателя, конденсатора, мерников конденсата, баков, щелоковой ловушки;

- экспериментальный выпарной аппарат с падающей пленкой. ,

Все установки были укомплектованы контрольно-измерительными

приборами и устройствами, обеспечивающими измерение температур стенки трубы и раствора, плотности теплового потока по длине трубы, параметров <

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

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

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

В четвертой главе приведены результаты экспериментальных исследований теплоотдачи и гидродинамики при кипении и испарении пенных растворов и воды при низких давлениях и движущихся с малыми скоростями в вертикальных трубах. Исследования проводились на установках в следующем диапазоне режимных параметров: при подъемном течении - концентрация черного сульфатного щелока на входе в трубу Ь = 0 + 50 % абсолютно сухих веществ (а.с.в.); водного предгидролизата - 0+30 % а.с.в.; давление греющею пара Ргр = 0,175 + 0,335 мПа; давление вторичного пара в сепараторе установки Рвт = 0,081+0,13 мПа; скорость циркуляции раствора м>п = 0,013+0,4 м/с; температура раствора на входе в трубу поддерживалась равной 96 °С; при этом плотность теплового потока изменялась в пределах ц = 2+100 кВт/м2; при гравитационном течении - <7 = 5+35 кВт/м ; объемная -

плотность орошения - Л, = 0,0001 +0,002 м3/(мхс). При этом выдерживалось равенство критериев кипения для стендовых установок и промышленных ВА ЦБК, диапазон изменения которых составлял: V

-=0,14+2,1; —2-=0,004+3,6.

Р -Г Уо Р-"о'Ын

Установлено, что интенсивность теплообмена а при всех гидродинамических режимах непостоянна по длине трубы. Отношение атах/ат]п достигало 50. Постоянные значения локальных а наблюдались только в случае течения некипящей жидкости. Во всех опытах, кроме режимов с кажущимся уровнем заполнения кипятильной трубы Ьур = 100 % при и>0 > 0,4 м/с, отмечалось падение теплоотдачи в верхней части трубы. При низких скоростях, менее 0,02 м/с и плотностях теплового потока более 40 кВт/м2, падение а отмечалось сразу же после входа раствора в трубу. Соотношение скорости цир-

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

Для выявления типовых режимов течения по данным настоящей работы нами были построены карты режимов течения двухфазного потока в координатах по - 1Ус1 и рп-х (массовая скорость - паросодержание). Из карт режимов следовало, что выпаривание черных сульфатных щелоков при тепловых режимах, в основном типовых для промышленных ВС, связано с наличием зоны ухудшенной теплоотдачи и, вследствие этого, интенсивным накипеоб-разованием. Экспериментальные данные, представленные в координатах рм> - х, показали, что возникновение пузырькового режима кипения воды ♦ происходит при более низких отрицательных паросодержаниях по сравне-

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

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

По полученным данным рассматривались зависимости: истинного ф, расходных - объемного /? и массового х паросодержаний от режимных параметров процесса по длине трубы Ь. Как показали эксперименты, влияние па-росодержания на теплоотдачу оказывается косвенным путем через изменение скоростей и структур двухфазного потока. Поскольку механизмы передачи теплоты, связанныес д и м>0, наилучшим способом учитываются

Рис. 1. Изменение интенсивности теплоотдачи в зависимости от паро-содержания па-рожидкостной смеси при К и- = 0,7:

1 - Ь=20 % а.с.в.;

2 - ¿=30 %;

3 - 5=40 %

критерием Кто влияние паросо-держания представлено на рис. 1 в функции а=А[(р) при сог^. По

этим кривым и ряду дополнительных признаков выделены четыре основных режима, с характерными соотношениями между q, wn, <р, определяющими доминирующее влияние того или иного механизма переноса теплоты на теплоотдачу. Область изменения паросодержаний с <р < 50 % и х < 0,005 определена как зона пузырькового кипения. При течении парорастворной смеси с низкими скоростями в режиме пузырькового кипения снижение а при увеличении wn или недогрева раствора до состояния насыщения АtH обусловлено уменьшением интенсивности процесса парообразования. Значительное повышение интенсивности теплоотдачи на участке пузырькового кипения связано с увеличением тепловой нагрузки. Отрицательное влияние на теплоотдачу увеличивающейся концентрации раствора в значительной мере компенсировалось пенообразованием, а независимость интенсивности теплоотдачи от паросодержания или скорости парожидкостной смеси объясняется определяющим влиянием механизма переноса теплоты, обусловленной процессом парообразования. Зависимость коэффициента теплоотдачи от режимных параметров процесса в области пузырькового кипения, представленная на рис. 2, с величиной достоверности аппроксимации R2 = 0,69 обобщается уравнением:

0^=111

А . 0.19 0.36

<4 '^0

(1)

где q= (7,0^86,5) кВт/м2; AtH = (0,1 *27) °С; w0= (0,01 +1) м/с.

1000

100

0,01 0,1 1 Скорость циркуляции, w[m/c]

♦ 1 Я2 A3

Рис. 2. Зависимость интенсивности теплоотдачи от скорости циркуляции:

1 - черный сульфатный щелок;

2 - вода; 3 - водный предгидроли-зат

Эмульсионный или пенный режим течения определялся в диапазоне параметров <р = (40^-70) %; х = (0,005 ^0,015); 1¥еэ < 11. На этом участке независимость коэффициента теплоотдачи от тепловой нагрузки объясняется низкой теплопроводностью пены. Увеличение скорости двухфазного потока, положительно сказывающееся на испарении междупленочной жидкости, способствовало разрушению пены и, следовательно, повышению ишенсивмости теплоотдачи. При возрастании скорости парорастворной смеси происходит наложение механизма турбулентного обмена в однофазной среде на механизм переноса, обусловленный процессом парообразования, в связи с чем а

увеличивается. Приведенная на рис. 3 зависимость интенсивности теплоотдачи от режимных параметров процесса в зоне пенного режима кипения с величиной достоверности аппроксимации Я2 = 0,54 описывается эмпирическим уравнением: апен = 4970 И^ • д ' 9 ; (2) где ^см = (0,045 -15,0) м/с; д = (10+110) кВт/м2.

10000

о

1000

100

10

100

1000

Тепловая нагрузка, кВт/м2

- Ряд1 ■ Ряд2 А РядЗ

Рис. 3. Зависимость интенсивности теплоотдачи от тепловой нагрузки: Ко - а\у~°

1 - черный сульфатный щелок;

2 - вода; 3 - водный предгидроли-зат

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

(80 ; 95) %; х > 0,015; \¥е0 >11, предшествующая резкому падению а, определялась как дисперсно-кольцевая. При д = (40*80) кВт/м2 тепловая нагрузка практически не влияла на теплообмен. С повышением концентрации раствора снижалась скорость движения его пристенной пленки, поэтому происходило уменьшение турбулентного обмена и, следовательно, интенсивности теплообмена. Приведенная на рис. 4 зависимость интенсивности теплоотдачи от режимных параметров процесса кипения в области дисперсно-кольцевого течения раствора с величиной достоверности аппроксимации

100000

£

10000

1000

10

Скорость смеси, м/с

-1 И2 АЗ

100

Рис. 4. Изменение интенсивности теплоотдачи в зависимости от изменения скорости парожидкостной смеси: 1 — черный сульфатный щелок; 2 - вода; 3 -водный предгид-ролизат; Ko = cu}-0■í)2blU,

1*2=0,71 описываются эмпирическим уравнением:

где М>. = (5-30) м/с; Ь = (0+50) % а.с.в.

СМ

Анализ опытных данных показал, что теплообмен при кипении пенных отработанных растворов в трубах выпарных аппаратов ЦБП определяется совместным воздействием режимных параметров (ц, \\>, А/„, Ь), которое проявляется через механизм формирования и перераспределения по высоте трубы зон с типовыми режимами кипения. Преимущественное развитие одной из зон может значительно изменить тенденцию теплообмена и быть причиной качественного различия данных, полученных при разных условиях. Длины участков с пузырьковым, эмульсионным и дисперсно-кольцевым режимами течения парощелокового потока увеличивались с повышением \у0 при <7=сопб1 и при юо > 0,07 м/с и д> 40 кВт/м2 участки могли занимать всю трубу. При повышении <7 до 100 кВт/м2 с одновременным уменьшением до ^0= 0,01 м/с на большей части трубы устанавливался дисперсный режим, и а резко падал. Наоборот, при увеличении до 0,1 м/с и снижении я до 5 кВт/м2 на большей части трубы преобладала экономайзерная зона. При повышении концентрации раствора режим ухудшенной теплоотдачи наступал при меньших <7, а снижение а проходило более плавно. В промышленных условиях подобное выделение указанных зон невозможно, в связи с чем управление технологическим процессом осуществляется изменением основных режимных параметров. Основным фактором, определяющим распределение указанных зон при кипении пенных растворов ЦБК, является соотношение между тепловым потоком и скоростью циркуляции раствора, поэтому критерий Км? = д/(гхр"х\\'о) принят в качестве меры, определяющей оптимальность процесса при выпаривании растворов ЦБК. Данные, представленные на рис. 5, дают возможность в зависимости от соотношения д и ум во вклад кипения выделить три характерных участка: развитого кипения при а =/(д); участка конвективного теплообмена при а = Ы'см); участка совместного влияния на теплоотдачу механизмов передачи теплоты за счет парообразования и конвективного обмена. Концентрация и величина недогрсва усиливают действие указанных механизмов передачи теплоты за счет зависимости от температуры вязкостных сил и сил межмолекулярного сцепления в растворе.

В области совместного влияния на теплоотдачу кипения раствора и его вынужденного движения для расчета коэффициента теплоотдачи может быть использована зависимость, предложенная А.А. Андреевским, В.М. Боришанским, В.Н. Фромзелем, Б.С. Фокиным:

а

= 1 + 7-10

-9

<*тр

Г-Р

Л3'2

апуз

канв У

(4)

где &тр -акоив + ^ '/а J ; апуз - коэффициент теплоотдачи при пузырьковом режиме кипения раствора, рассчитанный по зависимости (1); ™см - скорость парожидкостной смеси.

ф т г? =г

#§ £

£ Е £

о £

100000 10000 1000 100 10

шяшШ1 1 1 ♦

; ж ж жж! : ж ж хШР >

; ♦

0,01 0,1 1 10 Скорость потока, м/с

♦ 1 12 13 «4 ж5

100

Рис. 5. Теплоотдача при кипении воды и растворов в условиях вынужденного движения: 1 - расчетные данные для вынужденного движения однофазной жидкости; 2 - вода д=810 кВт/м2 [С.С.Кутателадзе]; 3 - вода д = 230 кВт/м2 [Д.А.Лабунцов]; 4, 5 - данные по кипящим растворам при д = 20 кВт/м2 и <7 = 80 кВт/м2 соответственно

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

м>/ V

V

Ч.УЯьР М Уя^У^'Л'

' g■l'a,r■fir■w'cp■it'cp■T,, I а ' уу'0' I

или в критериальной форме:

Яе, Рг, Рг, Kw, К, Кв, Ко, ^Д,-^,—Ей, Ыи .

\¥п I Р' X

(5)

(6)

По условиям эксперимента и независимости интенсивности теплообмена от отдельных характеристик исследуемого процесса из совокупности (6) были исключены Ёг, К, Кб, Ко, 1«А X /Л.' . Вместо них введены комплексы:

К1 =

р-Щ

Ргпл; а/ак

При этом аКОНв рассчитывался по формуле

для вынужденной конвекции в соответствии с режимом течения потока.

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

а Iаконе ~ '■> ^Г, РГ/1Л » ^/р'

(7)

Из результатов проведенного обобщения, представленного на рис. 6, видно, что устанавливаются две характерные области влияния безразмерных переменных в зависимости от критерия И. Для первой области, определяемой Ш < 0,05, при относительно низких тепловых нагрузках и высоких значениях недогрева характерно увеличение интенсивности теплоотдачи с уменьшением величины недогрева. Во второй области, характеризующейся высокими тепловыми нагрузками, малыми значениями недогрева и определяемой Ш > 0,05, интенсивность теплообмена возрастает с увеличением величины недогрева. Полученная зависимость с величиной достоверности аппроксимации И. = 0,59 описывается уравнением :

«/ = Рг°-6| Р / I /оч

где А = 12,5 • 10\ п = 0,214 при Ы < 0,05; А = 2910, п =-0,277 при ¡я > 0,05; к „ =0,1+2,0; к1 =0,005+3,5; Ргт= 1,5+25,0.

10000

о

1000

0,01

0.1

К!

♦1Я2л3

Рис. 6. Обобщение опытных данных по интенсивности теплоотдачи при кипении пенных растворов ЦБК:

Ко =-*-^

¿.«.453 р 0.6

1 - черный сульфатный щелок; 2 - водный пред-гидролизат; 3 - вода;

Опыты для изучения теплообмена при испарении стекающей пленки проводились при Ь = (0+60) % а.с.в. и температурах растворов, близких к температурам насыщения. Визуальные наблюдения показали, что процессы теплообмена находились в области нагрева и испарения с поверхности. В результате анализа данных эксперимента получены зависимости для средней и мак-

симальной скоростей пленки. По ним в исследуемом диапазоне значений чисел Рейнольдса идентифицированы характерные режимы течения пленки. Средние значения критических чисел Рейнольдса составили: Яе < 1000 - ла-минарно-волновое течение; Яе > 3000 - турбулентное течение. На основе данных по средней скорости пленки рассчитаны значения средней толщины пленки (рис.7). Из результата анализа профилей температур гравитационно

Рис. 7. Зависимость безразмерной средней толщины пленки жидкости 5' от числа Яе: 1 - данные настоящей работы, рассчитанные при Яе < 1600 по зависимости -= - *

( Л v

1/3

: 0,56 Re0-39.

при Re >1600 -8

S* =

< 2\,/3 у'

g

= 0,14 Re'

0.59

2 - по данным Г. Гимбути-

са; 3 - по данным И.

Федоткина; 4 - по данным К. Fiend; 5 - по данным

Е. Мальцева, В. Тимофеева;

6 - по данным Г. Гребера;

7 - по данным G. Folford

стекающей водяной пленки и локальных значений коэффициентов теплоотдачи при различных числах Рейнольдса следует, что в рассматриваемом диапазоне чисел Рейнольдса от 4х 10} до 14х 10 можно выделить два участка. Первый, на котором происходил существенный рост теплоотдачи до максимального значения на некотором расстоянии от входа в трубу до а™*, и второй - за точкой максимума, на котором наблюдалось незначительное снижение теплоотдачи вследствие продолжающегося роста турбулентного пограничного слоя. Затем фиксировалась относительная стабилизация. При более низких числах Рейнольдса проявлялся еще один участок - начальный, на котором теплоотдача уменьшалась сразу, начиная от входа раствора в трубу. Длина начального участка увеличивалась при уменьшении числа Ке (рис.8) и с увеличением числа Ре.

1900 а ! Вт/м2 °С

1 X

\ н

ООО V ♦

400 '-А-- и

Х -1

* -г

■ -э

о -4

£., см

Рис. 8. Локальные коэффициенты теплоотдачи по длине трубы при течении пленки щелока: 1 - Яе = 3045; 2 - Ле = 430; 3 - Яе = 255; 4 - Яе = 182

Данные по теплоотдаче к стабилизированным испаряющимся пленкам щелоков сульфатного производства с концентрациями до 60 % приведены на

Ш

рис. 9 в координатах Ко =/(Яе), где К„ = .

10000

«1 ■ 2 * 3 щ4 ж 5

Рис. 9. Зависимость интенсивности теплоотдачи от числа Яе:

„ _ ЛГи„ л„ — -

Рг0 75 '

1 - Ь = 60 % а.с.в.;

2 - Ь = 50 %;

3 - Ь = 45 %;

4 - Ь = 30 %;

5 - Ь = 20 %

Представленная на рис. 9 зависимость описывается уравнением:

= 0,0044Яе°-37Рг" %", (9)

где Ие = (0,2 * 10)х 103; Л- = 1,5+25,0.

Коэффициенты уравнения получены в результате процедуры регрессионного анализа. Величина достоверности аппроксимации составляет Я2= 0,83.

Анализ полученных данных показал неоднозначность влияния плотности теплового потока на теплоотдачу, отличающуюся от зависимостей, которые были получены для воды: для щелоков низкой концентрации и воды интенсивность теплоотдачи снижалась при увеличении плотности теплового потока; для щелоков с Ь = (35н-50) % интенсивность теплоотдачи практиче-

ски не зависела от плотности теплового потока; для щелоков высокой концентрации интенсивность теплоотдачи увеличивалась с возрастанием плотности теплового потока. Показатель степени «п» при критерии Кц изменялся от (-0,19) для случая испарения воды до 0,15 для щелока с Ь = 60 % и максимальном тепловом потоке.

В пятой главе представлены результаты исследований по динамике накипеобразования. Для определения влияния накипеобразования на тепловой режим работы станции, на изменение нагрузок поверхности испарения, на изменение концентраций щелоков по корпусам проведены годовые исследования режимов работы трех ВС Котласского ЦБК. Для характеристики процесса накипеообразования определялся термический коэффициент наки-

_ КН

пеобразования: V — ~—, (10)

где = — - условное термическое сопротивление накипи,

Л)„ Вт

8и - толщина слоя накипи, м; Хн - коэффициент теплопроводности накипи, Вт/(м. С); и = и,//г — массовое напряжение поверхности теплообмена, кг/(м2 хс); м> = 8(1 - Ь/В) - количество выпаренной в В А влаги , кг/с; 5 - количество раствора, подаваемого в ВА на выпаривание, кг/с; Ь, В- концентрации раствора на входе и выходе из В А, % а.с.в.; F — поверхность теплообмена В А, м2; г - время работы ВА после промывки, сутки.

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

Ки Кр

где Ки,Кр - коэффициент теплопередачи, полученный в результате исследования режимов работы ВС, и теоретический коэффициент теплопередачи, Вт/(м2х°С).

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

£ -£>о , (12)

где £ \{расч - расчетный температурный напор; - реальный темпе-

ратурный напор.

При достижении =0 обеспечение производительности ВС ста-

новится невозможным и, либо происходит очистка ее поверхности от накипи,

либо она будет работать с постепенно уменьшающейся производительностью.

Приведенные экспериментальные данные показывают, что условные термические сопротивления накипи составляют для большинства корпусов примерно 20 % (45 % - для отдельных) от полного термического сопротивления теплопередаче. На основании представленных данных появляются дополнительные возможности для применения наиболее рациональных мероприятий по экономии энергоресурсов, увеличения производительности ВС ЦБК, достигаемых при уменьшении накипи. По полученным данным определялся ежесуточный прирост термического сопротивления накипи для каждого корпуса по зависимости: =-, (13)

г

где - условные термические сопротивления накипи в ВА в конце и в начале эксперимента.

Постепенное увеличение Ян связано с увеличением содержания наки-пеообразователей в сульфатном щелоке и зависит главным образом от увеличения их концентрации в растворе. Среднесуточный прирост накипи (рис. 10) определяется зависимостью: /?„ = 1,658x10'3хихЬ037 . (14)

Рис. 10. Среднесуточный прирос! термического сопротивления накипи

Составлены номограммы для расчета оптимального периода работы между промывками ВС ЦБК.

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

при заданной производительности и, соответственно, повышение производительности всего комбината.

В шестой главе рассматриваются термические сопротивления передаче теплоты при конденсации = 1/асм: термическое сопротивление пленки конденсата Я^, термическое сопротивление фазового перехода Яф и диффузионное термическое сопротивление Яд. Гидродинамика течения пленки конденсата, определяющая интенсивность теплообмена при пленочной конденсации, имеет сходство вследствие постоянства плотности теплового потока с условиями течения и теплообменом при поверхностном испарении пленки, нагретой до температуры насыщения. Исходя из этого, а также считая, что физические параметры пленки конденсата по высоте трубы постоянны, коэффициент теплоотдачи при конденсации при указанных допущениях не будет зависеть от характера изменения температуры поверхности трубы. Следовательно, методика и результаты исследований гидродинамики и теплоотдачи от гравитационно стекающей и испаряющейся пленки воды внутри трубы могут быть использованы для определения термического сопротивления а пленки конденсата Я^1/а на наружной поверхности вертикальных труб. В исследуемом диапазоне значений чисел Рейнольдса нами было рекомендовано эмпирическое уравнение (9). На рис. 11 представлена эта зависимость в координатах Мим = /(Яе,^ и сопоставлена с данными, рассчитанными по уравнению Гимбутиса.

При продольном и направленном вниз движении пара силы трения,

0,4 0,35 0,3 0,25 0,2 0,15 0,1 0,05 0

- -- I 1 1

- — - - 1 " 1 1 ~ 1 . . . , 1

50

100

150 200

Явт ♦ 1 Я2 аз

250

300

350

Рис. 11. Теплоотдача при конденсации пара на поверхности вертикальной трубы: 1-данные эксперимента;

2 - теоретические данные, рассчитанные по уравнению Гимбутиса;

3 - данные, рассчитанные по уравнению (9)

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

Появление термического сопротивления фазового перехода обусловлено скачком температуры на границе раздела паровой и жидкой фаз. При малом коэффициенте конденсации скачок может быть значительным, особенно при низких давлениях. Сопротивление Яф может быть сопоставимо с термическим сопротивлением пленки конденсата Я^. Известно, что при конденсации чистого водяного пара при давлении Р > 104 Па в диапазоне работы ВС значением температурного скачка на границе раздела фаз пренебрегают.

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

де __~^ппов___ _ _^ ПО ~ 1Г пов___(15)

адОпО ~ ' п.те) + г ' Уп пов —'я.яое ) + г ' ~ тп ло«)[г + (' ис по« "'*«)]

где ал - конвективный коэффициент теплоотдачи от парогазовой смеси к пленке конденсата; г - теплота парообразования;]п.пов ~ поток массы пара у поверхности пленки конденсата; /? - коэффициент массоотдачи; т„о и тп т№ - относительные массовые концентрации пара в объеме и у поверхности.

Из уравнения (15) следует, что коэффициент теплоотдачи асм, куда входит Яд, зависит от интенсивности взаимосвязанных процессов тепло- и массообмена в парогазовой смеси и в пленке. Для расчета коэффициента асм предварительно определяют на основании аналогии процессов тепло- и массообмена коэффициент массоотдачи Д. Как следует из рис. 12, средний коэффициент массоотдачи для вертикального пучка труб в диапазоне изменения параметров работы ВС ЦБК определяется по зависимости В. А. Блюхина (16) с учетом полученного в результате исследований коэффициента с = 0,46:

к, п

О = 0.5 1/3 , (16)

е, 'п - тс

о

Рн1/

где - с = 0,46; л0 = (рпо - р„.пов); йэкв=4^П; £го = ур - количество неконденсирующихся газов.

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

t !

воды и растворов показывают, что, устанавливая правильные соотношения между количествами теплоты и раствора, подаваемых в аппарат, можно добиться высоких коэффициентов теплоотдачи, в том числе и в аппаратах, выпаривающих раствор высокой концентрации. Характеристикой этих соотношений является кажущийся уровень заполнения трубы, измеряемый наиболее просто по водомерному стеклу. Используя экстремальный характер зависимости а = /(Ьк), представленной на рис. 13, можно определить наиболее экономичные, с точки зрения передачи теплоты, режимы работы кипятильных труб. Для концентрации щелока до 35 % оптимальными являются

Рис. 12.

Массоотдача при конденсации пара из парогазовой смеси на трубах ВА:

о

0,01 0,1 1 Количество неконденсирующихся газов

*1 • 2 А 3 щ4

Я",

I/3

Re0'5 ' 1 - расчет по (16); 2,3,4 - данные эксперимента - ВС - 1, ВС - 2 и ВС - 3

уровни, составляющие 20-4-40 % высоты трубы, а для щелоков более высокой концентрации - уровни порядка 50^70 %. Таким образом, в аппаратах, выпаривающих плотные щелока, следует поддерживать повышенный кажущийся уровень.

о 6000 ел

80

«1 ■ 2 д 3

Рис. 13.

Изменение

среднего

коэффициента

теплоотдачи в

зависимости от

кажущегося

уровня

заполнения

кипятильной

трубы щелоком

концентрацией:

1 - 20 % а.с.в.;

2-30 %;

3-40 %

Из рис. 13 следует также, что в зависимости от Ьк для концентрированных щелоков значение а может быть выше, чем для средних и разбавленных. Кажущийся уровень в однопроходных ВА можно определить, кроме показаний водомерного стекла, по показанию специально установленного измерительного прибора. По этому показателю можно управлять режимами работы ВС с помощью средств автоматического регулирования.

Рис. 14. Значения необходимого (1) и реального (2) изменения скорости циркуляции раствора в ВА по мере его выпаривания: п - номер корпуса ВС по ходу движения раствора

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

Рис.15. Зависимость предельного влагосодержания от критериевАги Рг

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

Аг0,22

(1-^=025—^. (17)

рг».«> ■

Зависимость, характеризующая уравнение (17), представлена на рис. 15.

Механизм кризиса кипения рассматривают, как правило, в координатах Чкр~Хщ). Величиной, характеризующей начало кризиса теплообмена первого

f i

£ £

о-

(В

I

¡I

i К

■ ---

■

•

г » ■

• Л

-в-

-0.1

-0,05

0,05

0,1

Паросодержание, X

♦ 1 ■ 2

Рис. 16. Зависимость Чкр=/(х) при кипении сульфатного щелока:

1-Ь < 20 %;

2-Ъ > 30 %

рода в области отрицательных и малых значений паросодержаний, может являться паросодержание Хкр. На основании опытных данных, полученных на стендовой установке с разной концентрацией раствора и w0 = 0,015+0,025 м/с, pw = 25 кг/(м2хс), на рис. 16 построена зависимость q = f(x) для кризисных условий при кипении пенообразующих растворов ЦБК. Кризис фиксировался в точке роста температуры стенки трубы и одновременного падения коэффициента теплоотдачи. lía рис. 16 показано, что изменения величины критической плотности теплового потока находятся в очень малом интервале значений паросодержания х вследствие пенообразующих свойств растворов.

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

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

Перечисленные механизмы образования кризиса теплообмена при кипении пенообразующих растворов наблюдаются до значений паросодержа-ний х « 0,05. Для первой рассмотренной области кризиса теплообмена первого рода при х < 0,001 получено уравнение

дКр - 0,56(1 + 6)1Л (ан + 0,5)ол , [кВт/м2 ] , (18) где Ь = 10^60 % а.с.в. - концентрация раствора; = 1 + 30 °С - недогрев раствора до температуры кипения.

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

Критическая плотность теплового потока при кипении пенообразующих растворов ЦБК на участке х - 0,001 +0,005, где кризис теплообмена первого рода возникает в случае образования в парогенерирующей трубе микропленки в условиях неравновесного изменения расхода жидкости в пристенной пленке и ее высыхания, описывается уравнением:

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

Поскольку исследовались режимы кипения растворов в диапазоне изменения параметров работы ВА ЦБК, область кризисных явлений, соответствующая дальнейшему увеличению паросодержания не изучалась.

ОСНОВНЫЕ ВЫВОДЫ

По результатам проведенных промышленных, стендовых и теоретических исследований режимов работы ВС сделаны следующие выводы:

1. Предложены пути совершенствования энергоиспользования при применении прот ивоточной схемы ВС на ЦБК. Экспериментально показана возможность и даны параметры безнакипного выпаривания пенных отработан-

ных растворов ЦБП в длиннотрубных аппаратах с гравитационным течением выпариваемого раствора.

2. Разработаны средства, методики для экспериментального исследования процессов теплообмена, создан экспериментальный комплекс для проведения стендовых испытаний и проведены комплексные исследования теплопередачи при выпаривании восходящих и гравитационно-стекающих пенных растворов ЦБК на стендовых установках и промышленных ВС.

3. На основе экспериментальных данных разработаны: физическая модель движения многокомпонентных пенных двухфазных парорастворных потоков в кипятильных трубах ВА в условиях восходящего и гравитационного течения; определено влияние физико-химических характеристик черного сульфатного щелока, водного предгидролизата и режимных параметров работы ВА на интенсивность выпаривания.

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

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

6. Предложена методика расчета для определения плановых промывок

ВС.

7. Экспериментально установлены особенности теплоотдачи к кипящим сильнопенящимся растворам в трубе: отмечено быстрое развитие кризиса теплоотдачи как при высоких (д < 100 кВт/м2), так и при низких (<7 < 10 кВт/м2) тепловых потоках при определенных режимных параметрах; в качестве параметров, характеризующих предельный режим работы испарительных каналов с подъемным движением двухфазного потока, рекомендовано принимать предельное влагосодержание (1-<р)пр и критические тепловые нагрузки qKp•, для каналов с гравитационно-стекающей пленкой жидкости режим начала ухудшения теплоотдачи характеризуется предельной плотностью орошения ^пред I предложены расчетные зависимости для определения величин (1-<р), Чкр> ГПред •

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

9. Результаты исследований по динамике накипеобразования, приведенные в настоящей работе в виде номограмм и расчетных зависимостей, являются основой для управления работой ВУ ЦБК, обеспечивают максимальную продолжительность их эксплуатации между стадиями промывки, в результате чего повышается производительность ВС и всего комбината.

10. Обобщенные методы расчета теплогидродинамических процессов, основные концепции комплексной оптимизации энергоемких технологических звеньев теплотехнологических схем, предложения по снижению наки-пеобразования на кипятильных поверхностях труб в виде программного метода расчета ВС используются при модернизации, разработке, внедрении и промышленной отладке ВС ЦБК, ЗАО «ГИПРОБУМ», ГУП «ЦНИИБум-маш». Результаты исследований используются в учебном процессе СПбГТУРП.

Основные положения и результаты диссертационной работы изложены в следующих публикациях:

1.Суслов В.А.Сопоставление тепловой эффективности выпарных станций сульфатного производства, поставляемых на предприятия ЦБП зарубежными фирмами [Текст] / В.А. Суслов, C.B. Антуфьев, B.JI. Мовсесян //Машины и оборудование ЦБП: межвуз. сб. научн. тр., вып. 5.-JI., JITA, 1977.-С. 131-133.

2.Суслов В.А. Результаты сопоставления опытных и расчетных данных по теплообмену при кипении сульфатного щелока[Текст] / В.А. Суслов //Машины, конструирование, расчеты и оборудование ЦБП: межвуз.сб. научн. тр., вып. 6.- Л., ЛТА, 1978.-С. 96-101.

3.Суслов В.А. Стенд для исследования механизма теплообмена и гидродинамики при пленочном кипении сульфатного щелока [Текст] / В.А. Суслов, В.И. Куш-нир //Машины, конструирование, расчеты и оборудование ЦБП: межвуз. сб. научн. тр., вып. 6.-Л., ЛТА, 1978.-С. 102-106.

4.Мовсесян В.Л. Выпарные аппараты и схемы выпарных станций ЦБП [Текст]: обзорн. информ. ЦИНТИхимнефтемаш / Мовсесян В.Л., Суслов В.А.- М., 1979,50 с.

5.Суслов В.А. Анализ процесса теплоотдачи при кипении в вертикальных трубах пленочных выпарных аппаратов [Текст] / В.А. Суслов, В.А. Веденеев, Г.В. Ри-жинашвили // Машины, конструирование, расчеты и оборудование целлюлозно-бумажных производств: межвуз. сб. научн. тр., вып. 7.- Л., ЛТА, 1979.-С. 127-130.

6.Суслов В.А. Исследование теплообмена при кипении слабого черного сульфатного щелока в вертикальных трубах выпарных аппаратов [Текст] / В.А. Суслов // Машины и аппараты ЦБП: межвуз. сб. научн. тр., вып. 8.- Л., ЛТИ ЦБП. 1980.-С. 102-105.

7.Суслов В.А. Исследование процессов теплообмена при испарении щелоко-вых растворов в производстве бумаги [Текст] / В.А. Суслов // Химия и технология бумаги: межвуз. сб. научн. тр.,вып. 9.-Л., ЛТИ ЦБП, 1981.-С. 137-141.

8.Суслов В.А.Исследование теплоотдачи при выпарке черного сульфатного щелока [Текст] / В.А. Суслов, Г.В. Рижинашвили //Бумажная промышленность,-1981.-№ 12.-С. 20-21.

9.Суслов В.А. Результаты промышленных исследований режимов работы выпарных станций ЦБП [Текст] / В.А. Суслов // Машины и аппараты ЦБП : межвуз. сб. научн. тр., вып. 10.-Л., ЛТИ ЦБП, 1982.-С.133-135.

Ю.Суслов В.А. Анализ теплообмена при кипении сульфатного щелока по усредненным опытным данным [Текст] / В.А. Суслов //Машины и аппараты ЦБП: межвуз.сб. научн.тр., вып. ll.-JI., 1983.-С. 113-115.

11. Суслов В.А. Позонное исследование теплоотдачи при кипении водного предгидролизата в вертикальной трубе [Текст] / В.А. Суслов, В.И. Кушнир; ВНИ-ПИЭИлсспром. - М., 1985.- Деп. № 1649-лб.

12.Суслов В.А. Исследование теплообмена при выпаривании черных сульфатных щелоков [Текст] / В.А. Суслов //Теплофизика и гидродинамика процессов кипения и конденсации: материалы Всесоюзн. конф.,т.1,ч.4- Рига, РПИ, 1985.-С. 117-119.

В.Суслов В.А. Проектирование выпарных установок ЦБП [Текст]: учеб. пособие / В .А. Суслов, В .Л. Мовсесян, А.Ф. Мурзич, В.А. Ганичев.-Л.: ЛТА., 1987.- 83 е.- (В надзаг.: ЛТИ ЦБП).

14.Суслов В.А. Интенсификация работы выпарных установок глиноземного производства [Текст] / В.А. Суслов, Е.М. Милин, А.М. Насопов, A.A. Яновский, Д.Р. Сандлер // Цветная металлургия.- 1988.- № 6.- С.31-35.

15. A.c. 1432122 СССР. Устройство для регенерации тепла и химикатов и парогазовых выбросов растворителя плава содорегенерационного котлоагрегата [Текст] / Л.В. Романова, A.B. Наумов, В.А. Суслов и др.(СССР).- 1988.

16.Суслов В.А. Ступенчатый метод расчета выпарных аппаратов [Текст] / В.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на промпред-приятиях и ТЭС: межвуз. сб. научн. тр.- Л., ЛТИ ЦБП, 1988.-С. 8-13.

17. Мурзич А.Ф. Промышленные тепломассообменные процессы и установки в ЦБП [Текст]: учеб. пособие / А.Ф. Мурзич, В.Л. Мовсесян, В.А. Суслов, Е.К. Гу-сев.-Л.: ЛТИ ЦБП, 1989.- 80 с.

18.Суслов В.А. Исследование теплоотдачи при кипении черного сульфатного щелока в экономайзерной зоне выпарных аппаратов [Текст] / В.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС : межвуз. сб. научн. тр.- Л., ЛТИ ЦБП, 1989.-С. 10-13.

19. Суслов В.А. Расчет коэффициента теплопередачи в выпарных аппаратах отрасли [Текст] / В.А. Суслов, Е. В. Алексеев // Бумажная промышленность.- 1989.-№ 2.-С. 30-31.

20.Суслов В.А. Стенд для исследования процесса накипеобразования в кипятильных трубах выпарных аппаратов ЦБП [Текст] / В.А. Суслов, В.В. Семенов // Машины и аппараты целлюлозно-бумажного производства: межвуз. сб. научн. тр.-Л„ ЛТИ ЦБП, 1989, с. 91-92.

21.Суслов В.А. Теплоотдача при кипении черного сульфатного щелока в вертикальных трубах выпарных аппаратов целлюлозно-бумажной промышленности [Текст] / В.А. Суслов //Двухфазный поток в энергетических машинах и аппаратах : тез. докл. 8-й Всесоюзн. конф. - Л., 1990 .-С. 352.

22.Суслов В.А. Исследование процесса накипеобразования в трубах выпарных аппаратов батарей Котласского ЦБК [Текст] / В.А. Суслов, В.В. Семенов, Л.А. Михайлов //Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС : межвуз. сб. научн. тр. - Л.,ЛТИ ЦБП, 1990. -С. 6-9.

23.Романова Л.В. Установка для конденсации парогазовой смеси из бака плава содорегенерационных хотлоагрегатов [Текст] / Л.В. Романова, В.А. Суслов // V Минский международный форум по тепло-и массообмену : тез. докл. и сообщ.-Минск, 2004.-С. 119-120.

24.Суслов В.А. Обобщение усредненных по поверхности теплообмена опытных данных, полученных при выпаривании водного предгидролизата [Текст] / В.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС: межвуз. сб. научн. тр. - Л., ЛТИ ЦБП, 1990.-С. 9-12.

25. A.c. 1646087 СССР. Выпарной аппарат [Текст] / В.А. Суслов и др.(СССР) .1991.

26.Суслов В.А. Оптимальные скорости подачи щелока в трубы выпарных аппаратов сульфат-целлюлозного производства [Текст] / В.А. Суслов, Г.В. Рижина-

швили // Проблемы экономии топливно-энергетических ресурсов на промпред-прияиях и ТЭС : межвуз. сб. научн. тр.- СПб., СПбТИ ЦБП, 1992.- С. 19-21.

27. A.c. 1829154. Выпарной аппарат для кристаллизующихся растворов [Текст] / В .А. Суслов и др. -1992.

28.Суслов В.А. Параметрическое уравнение для расчета коэффициента теплопередачи при выпаривании черных сульфатных щелоков [Текст] / В.А. Суслов, Г.В. Рижинашвили // Проблемы экономии топливно-энергетических ресурсов на пром-предприяиях и ТЭС : межвуз. сб. научн. тр. - СПб., СПбГТУ РП, 1993. - С. 14-18.

29.Суслов В.А. Теплообмен и гидродинамика пленочного течения в выпарных аппаратах целлюлозного производства [Текст] / В.А. Суслов, Г.В. Рижинашвили // Кипение, кризисы кипения, закризисный теплообмен : тр. 1-й Рос. нац. конф. по теплообмену. Т. 4,- М„ МЭИ, 1994,- С. 225 -229.

30.Суслов В.А. Анализ технико-экономических показателей различных схем выЛарных установок целлюлозного производства [Текст] / В. А. Суслов, Г.В. Рижинашвили // Целлюлоза. Бумага. Картон. -1994.- № 11-12,- С. 21-22.

31.Суслов В.А. Способ перевода выпарных аппаратов на режим работы с нисходящей пленкой раствора [Текст] / В.А. Суслов, Г.В. Рижинашвили // Организация природоохранной деятельности на предприятиях и пути ресурсосбережения: тез. докл. Рос. научно-практ. конф.- СПб., 1996.-С. 23-25.

32.Суслов В А. Влияние термокапиллярного эффект на теплопередачу в гравитационно стекающей пленке жидкости [Текст] / В.А. Суслов, В.Ю. Бойков, Г.В. Рижинашвили // Проблемы экономии топливно-энергетических ресурсов на пром-предприятиях и ТЭС :межвуз. сб. научн. тр.- СПб., СПбГТУ РП, 1997.- С. 82-86.

33.Суслов В.А. Особенности работы распределительных устройств типа спиральных вставок для выпарных аппаратов с падающей пленкой [Текст] / В.А. Суслов, Г.В. Рижинашвили, В.Ю. Бойков // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС :межвуз. сб. научн. тр.-СПб., СПбГТУ РП, 1997,- С. 51-55.

34.Суслов В.А. Исследование теплоотдачи при кипении сульфатного щелока в трубах выпарных аппаратов в зоне кольцевого течения [Текст] /В.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на предприятиях и ТЭС : межвуз. сб. научн. тр.- СПб., СПбГТУ РП, 1997.-С. 10-13.

35.Суслов В.А. Основные процессы при выпаривании щелоков целлюлозного производства [Текст] : учеб. пособие для вузов /В.А. Суслов. - СПб.: СПбГТУ РП, 1998,- 92 с.

36.Суслов В.А Повышение эффективности сжигания сульфатных щелоков целлюлозно-бумажного производства [Текст] / В.А. Суслов, Г.В. Рижинашвили, В.Ю. Бойков // Опыт и перспективы развития ресурсосберегающих технологий и утилизации твердых отходов на предприятиях различных форм собственности: материалы научно-практ. конф.- СПб.,СПбГТУ РП, 1998.-С. 94-96.

37.Суслов В.А. Теплообмен на начальном участке трубы при нагревании падающей пленки жидкости [Текст] / В.А. Суслов, Г.В. Рижинашвили, В.Ю. Бойков // Кипение, кризисы кипения, закризисный теплообмен. Испарение. Конденсация : тр. 2-й Рос. нац. конф. по теплообмену . Т.4.- М., МЭИ, 1998.-С. 362-363.

38.Суслов В.А. Оптимизация работы выпарных станций ОАО «Котласский ЦБК»[Текет] / В.А. Суслов, М.А. Чайка // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС: межвуз. сб. научн. тр.-СПб., СПбГТУ РП, 2002.-С. 64-69.

39.Суслов В.А. Исследование паросодсржания двухфазного потока при кипении щелока целлюлозно-бумажной промышленности в длинных трубах [Текст] /

В.А. Суслов, М.И. Чайка // 3-я Российская национальная конференция по теплообмену: программа. - М., МЭИ, 2002.-С. 48.

40.Романова Л.В. Установка для конденсации парогазовой смеси из бака плава содорегенерационных котлоагрегатов [Текст] / Л.В. Романова, В.А. Суслов // V Минский международный форум по тепло-и массообмену : тез. докл. и сообщ.-Минск, 2004.- С. 119-120.

41.Суслов В.А. Кризисы теплообмена при кипении щелочных растворов в трубах выпарных аппаратов ЦБП [Текст] / В.А. Суслов // Современная наука и образование в решении проблем экономики европейского севера: материалы Международной научно-техн. конф., посвященной 75-летию АЛТИ-АГТУ. Т. 1. - Архангельск, 2004,- С. 371-373.

42.Суслов В.А. Повышение надежности выпарных аппаратов при выпаривании черных сульфатных щелоков [Текст] / В Л. Суслов // Целлюлоза. Бумага. Кар-тон,-2004.-№ 10,- С. 58-61.

43 .Суслов В.А. Основные процессы при выпаривании щелоков целлюлозного производства [Текст]: учеб. пособие для вузов / В.А. Суслов. - СПб.: СПбГТУ РП, 2004,-102 с.

44.Суслов В.А. Результаты промышленных исследований режимов работы выпарных станций ЦБП и их сопоставление по эффективности выпаривания [Текст] / В.А. Суслов // Изв. ВУЗов. Лесной журнал,- 2005.- № 1-2,- С. 168-172.

45.Суслов В.А. Определение коэффициентов накипеобразования на выпарных батареях целлюлозно-бумажных комбинатов [Текст] / В.А. Суслов // Изв. ВУЗов. Лесной журнал.- 2005.- № 1-2.- С. 86-99.

46.Суслов В.А. Исследование теплообмена при кипении отработанных щелочных растворов в трубах выпарных аппаратов [Текст] / В.А. Суслов // Теплоэнергетика.- 2005.- № 8.- С. 68-71.

47.Су слов В.А.Зависимость локальных коэффициентов теплоотдачи при кипении от параметров и структур двухфазных пенообразующих водных растворов [Текст] / В .А. Суслов // Теплоэнергетика.- 2005.- № 9.- С. 74-77.

48.Суслов В.А. Исследование гидродинамики и теплоотдачи при гравитационном отекании пленки сульфатного щелока в вертикальной трубе выпарного аппарата [Текст] / В.А. Суслов, Г.В. Рижинашвили // Целлюлоза. Бумага. Картон.-2005.-№5.-С. 68-71.

49.Суслов В.А. Определение коэффициента теплоотдачи при конденсации парогазовой смеси на вертикальном пучке труб выпарных аппаратов [Текст] / В.А. Суслов // Ресурсо- и энергосбережение в целлюлозно-бумажной промышленности и городском коммунальном хозяйстве: материалы Международной научно-практической конференции: сборник трудов, СПб., 2005.-С. 295-300.

Условные обозначения

а - плотность теплового потока, Вт/м2; г - скрытая теплота парообразования, кДж/кг; р"- плотность пара, кг/м3; - скорость циркуляции жидкости, м/с; а - коэффициент теплоотдачи при кипении раствора в трубе; а,™, - коэффициент теплоотдача при вынужденном движении однофазной жидкости, Вт/(м2°С); а - коэффициент поверхностного натяжения, н/м; р - плотность жидкости, кг/м3; X - коэффициент теплопроводности, Вт/(м°С); V - коэффициент кинематической вязкости, м /с; а - коэффициент температуропроводности, м2/с; ¿1н- недогрев жидкости до

11-2 190

состояния насыщения, °С; сР - теплоемкость жидкости, кДж/(кг*°С); <р - истинное объемное паросодержание;

р - расходное объемное паросодержание; <1 - диаметр трубы, м; Ь - длина кипятильной трубы; х - весовое паросодержание потока; - скорость парожидкостной смеси, м/с; В - концентрация раствора, % а.с.в.; pw- массовая скорость потока, кг/(м2с); % - ускорение свободного падения, м/с; Гу- объемная плотность орошения, м3/(м.с); 5* - средняя толщина пленки при стабилизированном течении, м; - кажущийся уровень заполнения кипятильной трубы, м.

Критерии подобия

Ы (у1}"*

Ыи = — - критерий Нуссельта; Ыия я =— — - модифицированное число X )

Нуссельта; К1 = . ^ ,. ; Км* = —%-; Ко = ———; К = ~—критерии фазо-

Р'^-Л!, Г-р гр-у с-Д/

г „ „ Ы 4Г 4 Г

во го превращения; Кн ---кретерий иедогрева; Яе = — = — = —- критерий

Ы, у ру V

и>2 _

Рейнольдса; /т = — - критерий Фруда; Рг = у/а - критерий Прандгля; Wэ —

{"омнУр .

Jcr-g-p

эмульсионное число Вебера.

СУСЛОВ Вячеслав Александрович

ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ВЫПАРИВАНИЯ ОТРАБОТАННЫХ ВАТОЧНЫХ РАСТВОРОВ ЦЕЛЛЮЛОЗНОГО ПРОИЗВОДСТВА

Автореферат

Подписано к печати 13.01.2006. Сдано в производство 16.01.2006. Формат 60x84/16 _Уч--изд. л. 2,0. Тираж 100 экз. Заказ № 961._

Ризограф ГОУ ВПО СПбГТУРП. 198095, Санкт-Петербург, ул. Ивана Черных, 4.

ВВЕДЕНИЕ

Глава 1 СОСТОЯНИЕ И РАЗВИТИЕ ИССЛЕДОВАНИЙ ПО ИЗУ- 14 ЧЕНИЮ ТЕПЛОПЕРЕДАЧИ ПРИ ВЫПАРИВАНИИ ВОДЫ И РАСТВОРОВ В ТРУБАХ ВЫПАРНЫХ АППАРАТОВ

1.1 Аппаратурное исполнение станций для выпаривания отра- 14 ботанных растворов целлюлозного производства

1.2 Теплообмен при кипении восходящего потока воды в тру- 21 бах

1.2.1 Механизмы кипения в условиях восходящего потока и ин- 21 тенсивность теплообмена, усредненная по поверхности трубы

1.2.2 Классификация восходящих режимов течения

1.2.2.1 Теплообмен в зоне пузырькового кипения

1.2.2.2 Теплообмен в зоне дисперсно-кольцевого течения жидко- 29 сти

1.3 Нисходящее течение жидкости в трубе. Классификация 31 режимов течения

1.3.1 Гидродинамические режимы течения падающей пленки

1.3.1.1 Ламинарный режим

1.3.1.2 Волновое течение

1.3.1.3 Турбулентное течение

1.3.2 Средняя толщина квазистационарной пленки

1.3.3 Влияние теплового потока на режимы парообразования в 37 пленке

1.3.3.1 Теплоотдача при нагреве пленки в условиях ее течения по 38 вертикальной поверхности

1.3.3.2 Теплоотдача при поверхностном испарении пленки, нагре- 40 той до Тн

1.3.3.3 Теплоотдача при пузырьковом кипении в пленке 43 жидкости

1.3.3.4 Кризисные явления в пленке жидкости при теплообмене

1.3.3.5 Другие факторы, влияющие на теплообмен при гравита- 50 ционном течении пленки

1.4 Теплообмен при кипении растворов

1.4.1 Теплообмен при кипении отработанных растворов ЦБП

1.4.1.1 Физические свойства черных сульфатных щелоков

1.4.1.2 Характеристика и состав предгидролизата

1.4.1.3 Накипеобразование в трубах выпарных аппаратов ЦБП

1.4.2 Анализ работ по теплообмену при кипении отработанных 63 растворов ЦБП

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

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

Глава 2 ПРОМЫШЛЕННЫЕ ИССЛЕДОВАНИЯ РЕЖИМОВ РА- 72 БОТЫ ВЫПАРНЫХ УСТАНОВОК ЦБП

2.1 Методика проведения теплотехнических исследований те- 72 плообмена на выпарных станциях

2.2 Исследования теплообмена на выпарных станциях ЦБП 73 Выводы

Глава 3 КОНСТРУКЦИИ ОПЫТНЫХ СТЕНДОВ И МЕТОДИКА 78 ЭКСПЕРИМЕНТА

3.1 Анализ процесса теплообмена; планирование теплотехни- 78 ческого эксперимента

3.2 Стендовая установка для исследования теплоотдачи и гид- 85 родинамики восходящего потока кипящей жидкости

3.2.1 Описание стендовой установки для восходящего потока 85 кипящей жидкости

3.2.2 Методика проведения опытов и обработки опытных дан- 87 ных при исследовании теплоотдачи к восходящему потоку

3.2.3 Тарировка стендовой установки на воде для восходящего 93 потока

3.3 Экспериментальная установка для исследования гидроди- 95 намики гравитационно стекающей пленки жидкости; методика эксперимента

3.3.1 Описание установки для исследования гидродинамики 95 стекающей пленки

3.3.2 Метод меток для исследования скоростей в пленке воды

3.4 Экспериментальная установка для исследования теплоот- 106 дачи к стекающему потоку

3.4.1 Описание установки для исследования теплоотдачи к сте- 106 кающему потоку

3.4.2 Методика определения локальной теплоотдачи и длины 106 начального участка

3.4.3 Определение температуры стенки и вычисление коэффи- 110 циентов теплоотдачи

3.4.4 Тарировка установки для исследования теплоотдачи к сте- 110 кающему потоку

3.4.5 Оценка погрешностей эксперимента

3.5 Стенд для исследования процесса накипеобразования в ки- 115 пятильных трубах выпарных аппаратов ЦБК

3.6 Экспериментальный выпарной аппарат с падающей плен

Выводы

Глава 4 ИССЛЕДОВАНИЕ ТЕПЛООТДАЧИ И ГИДРОДИНАМИКИ ПРИ КИПЕНИИ МНОГОКОМПОНЕНТНЫХ ПЕННЫХ РАСТВОРОВ В ТРУБАХ ПОД НИЗКИМ ДАВЛЕНИЕМ

4.1 Восходящий поток кипящей жидкости в трубах.

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

4.1.2 Основные режимы течения двухфазных потоков при кипении пенных растворов

4.1.3 Теплоотдача к кипящим растворам при пузырьковом режиме течения

4.1.4 Теплоотдача к кипящим растворам при эмульсионном (пенном) режиме кипения

4.1.5 Теплоотдача к кипящим растворам при дисперсно-кольцевом режиме течения

4.1.6 Анализ результатов исследований теплообмена при восходящем парожидкостном потоке

4.2 Гравитационное течение жидкости

4.2.1 Визуальные наблюдения

4.2.1.1 Опыты на воде

4.2.1.2 Опыты на черном сульфатном щелоке

4.2.2 Исследования гидродинамики пленки

4.2.3 Исследование теплообмена в пленке жидкости

4.2.3.1 Конвективный теплообмен при нагревании воды

4.2.3.2 Теплообмен при поверхностном испарении воды

4.2.3.3 Теплообмен при нагревании и испарении щелока

4.2.3.4 Влияние плотности теплового потока на теплоотдачу к пленке жидкости

4.2.3.5 Локальная теплоотдача и участок стабилизации Выводы

Глава 5 ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА НАКИПЕОБРАЗО-ВАНИЯ НА ВЫПАРНЫХ СТАНЦИЯХ ЦБК

5.1 Методика проведения опытов

5.2 Исследование режимов работы выпарных батарей для определения термического сопротивления накипи

5.2.1 Исследование режимов работы выпарных батарей картон-но-бумажного производства Котласского ЦБК

5.2.2 Исследование режимов работы выпарной батареи TEC-III Котласского ЦБК

5.3 Интенсивность изменений режимов работы выпарных станции

5.4 Продолжительность работы выпарных батарей между промывками

Выводы

Глава 6 КОНДЕНСАЦИЯ ПАРОГАЗОВОЙ СМЕСИ НА ПУЧКАХ 184 ВЕРТИКАЛЬНЫХ ТРУБ

6.1 Термические сопротивления передаче теплоты при конденсации

6.1.1 Термическое сопротивление при пленочной конденсации 186 неподвижного пара

6.1.1.1 Результаты исследования гидродинамики и теплообмена 188 при течении водяной пленки

6.1.1.2 Влияние касательных напряжений на границе раздела фаз

6.1.2. Термическое сопротивление фазового перехода

6.1.3. Диффузионное термическое сопротивление 194 Выводы

Глава 7 ОПТИМАЛЬНЫЕ ЗНАЧЕНИЯ РЕЖИМНЫХ ПАРАМЕТ- 199 РОВ ПРИ ВЫПАРИВАНИИ

7.1 Зависимость оптимального кажущегося уровня заполнения 199 кипятильной трубы от скорости питания и тепловой нагрузки

7.2 Зависимость интенсивности теплообмена от истинного па- 205 росодержания и влагосодержания кипящего потока

7.3 Зависимость интенсивности теплообмена от массового 210 расходного паросодержания

7.3.1 Механизм кризиса теплообмена первого рода до точки ин- 211 версии

7.3.2 Механизм кризиса теплообмена первого рода за точкой 215 инверсии

7.4 Влияние плотности теплового потока на интенсивность на- 220 кипеобразования в условиях пленочного течения жидкости

7.5 Зависимость интенсивности теплообмена от условий вы- 222 паривания

7.6 Перспективы конструктивного развития выпарных 223 аппаратов

7.7 Выбор оптимального количества ступеней выпаривания и 225 распределение по ним поверхности испарения

7.8 Тепловой расчет выпарных станций 227 Выводы 229 ОСНОВНЫЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ 232 ПРОМЫШЛЕННРСТИ

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

О 20 40 60 80 100 120

Производство продукци отрасли, %

РисЛ. Производство промышленной продукции по отраслям в 2004 г.: 1.Легкая. 2.Электроэнергетика. 3 .Химическая и нефтехимическая. 4.Лесная, деревообрабатывающая и целлюлозно-бумажная. 5.Машино- строение и металлообработка. 6.Черная металлургия. 7.Промышленность строительных материалов. 8.Пищевая. 9.Топливная. 10.Промышленность в целом.

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

Несмотря на то, что ЦБП - некрупная отрасль в составе всей промышленности России, она является лидером российского лесного комплекса. Ее индекс роста производства составил в 2004 г. 105,4 %, при этом наблюдается рост производства всех основных видов целлюлозно-бумажной продукции. Прирост по товарной целлюлозе за 2004 г. составил - 4,5 %, по производству бумаги

6,1 %, по тарному картону - 6,7 %. Наибольший прирост достигнут по производству ящиков из картона - 17,9 %. Данные по динамике производства целлюлозно-бумажной продукции представлены в табл. 1.

Таблица 1.

Вид продукции 2003 2004 Темп роста, % к 2003

Целлюлоза по варке, тыс. т. 5751,81 5915,62 102,8

Целлюлоза товарная, тыс. т. 2301,01 2404,22 104,5

Бумага всех видов, тыс. т. 3654,82 3878,94 106,1

Бумага газетная, тыс. т. 1814,21 1978,99 109,1

Бумага офсетная, тыс. т. 449,21 475,95 106,0

Картон всех видов, тыс. т. 2520,23 2774,21 110,1

Картон тарный, тыс. т. 1872,18 1007,81 106,7

Ящики из картона, млн. кв. м. 1420,00 1674,99 117,9

Тетради школьные, млн. шт. 1434,26 1525,80 92,4

ЦБП имеет мощные теплоэлектростанции, крупные и дорогие водоподго-товительные сооружения, а также установки для очистки и обезвреживания сточных вод [161]. По потреблению тепла ЦБП занимает 4 место среди других отраслей промышленности, электроэнергии - 6, а по потреблению воды - 3. Расходы на цели энергетики составляют до 20 %.

Собственные ТЭЦ и котельные предприятий вырабатывают в настоящее время более 80 % всей потребляемой отраслью тепловой энергии и около 40 % всей потребляемой электроэнергии. Установленная мощность ТЭЦ достигла 20 млн. кВт. ТЭЦ ЦБП в виде топлива используют угли - 27 %, мазуты - 37 %, природный газ - 15 %, биотопливо в виде утилизируемых собственных отходов - 21 %: частично выпаренные черные сульфатные щелока; щепу; кору; опилки и другие вторичные энергоресурсы, за счет чего частично решаются задачи, возникшие после подписания Киотского протокола. Основной расход топлива направлен на производство тепловой (75 %) и электрической (17 %) энергии. На технологические цели расходуется до 6 % потребляемого топлива. До 2 % топлива расходуется на выработку тепловой энергии для коммунально-бытового потребления. Удельный расход топлива на выработку единицы тепловой и электрической энергии на энергетических ТЭЦ и котельных ЦБП вследствие большого количества устаревшего оборудования составляет «312 г/кВт*ч или 782 кг/ГДж, что в два раза выше удельного расхода топлива на современных энергетических ТЭЦ.

Усть-Илимский; Братский; Сыктывкарский лесопромышленные комплексы; Котласский; Кондопожский; Сясьский; Архангельский; Соликамский; Све-тогорский, Балахнинский целлюлозно-бумажные комбинаты (ЦБК) потребляют до 50 % расходуемого тепла в отрасли, а 8 из них потребляют 50 % топлива и электроэнергии.

В составе различных производств ЦБК наиболее крупными потребителями электроэнергии являются следующие: бумаги -19,6 %; целлюлозы — 15,3 %; древесной массы - 13,6 % и разных видов картона - 7,1 %. Наиболее крупными потребителями тепловой энергии являются производства: целлюлозы - 36,1 %; бумаги - 15,9 %; картона - 8,5 %. На их долю приходится около 60,5 % всего расхода потребляемой тепловой энергии в отрасли. Данные, представленные в табл. 2 по производству основных видов продукции на примере ОАО «Архангельский ЦБК», подтверждают данное положение. Поэтому искать резервы повышения эффективности производства, интенсификации производственных процессов, экономии, в том числе электрической и тепловой энергии, устранять ее потери, лишние расходы необходимо, прежде всего, в этих производствах.

Таблица 2

Показатель Ед. изм. 2003 2004 % к 2003

Варка целлюлозы т. 770745 788200 102,3

Бумага т. 80138 82079 102,4

Картон т. млн. м2. 250456 1671,6 255949 1712,8 102,2 102,5

Гофробумага и картон универсальный т. млн.м2. 183949 1441,5 183239 1448,5 99,5 102,6

Картонно-транспортная тара тыс. м2. 124279,5 128431,2 103,3

Плиты древесноволокнистые тыс. м2. 8223,4 7901,4 96,1

Тетради ученические тыс.шт. 331597 358924 108,2

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

В реальных эксплуатационных условиях все технологические и энергетические процессы являются взаимосвязанными. Поэтому совершенствование взаимосвязей между энергетикой и технологией производственных процессов должно в значительной мере способствовать качественному ведению технологического процесса, рациональному использованию топливно-энергетических и вторичных ресурсов, улучшению экологической обстановки и снижению себестоимости продукции. С помощью этих взаимосвязей на базе глубокого изучения тепловых и технологических процессов в ЦБП должно создаваться новое, экономически эффективное высокопроизводительное технологическое оборудование и осуществляться модернизация существующего. В соответствии с этим, работа выполнялась в рамках региональной программы фундаментальных и прикладных исследований по проблеме комплексного использования и воспроизводства лесных ресурсов на 1965 -s-1990 г.г. АН СССР по теме «Исследование процессов и разработка новых видов оборудования для производства бумаги и картона»; в рамках общесоюзной научно-технической программы Минлеспрома СССР по созданию и внедрению новой техники и технологии и государственной научно-технической программы "Комплексное использование и воспроизводство древесного сырья" 1991 * 1996 г.г.

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

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

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

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

Производительность ВС по выпаренной влаге определяется из уравнения

F -At п где R - термическое сопротивление теплопередатеплового баланса: W =

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

Состав омпибиою рвааскпь nooh*NO7S3 '

СлоБыи NtpHblU

Белый целок

Реокиив

Барм Промы&кя Отстой иым

OKVCAlxuC

Peel 6

Укр«ЩЙЕни»м мерный mtJWt, «диц.гО-»^ 1»

Н^СО^Он^гНоОУЬСоСО,^ Коусти в! оЬжиг

Зеленый ш,елок Пдв|

Выпорке

Эобсб ко

OpUlUb ояоба

Содорггене-рационный агрегат

1^бсио(ные реакции^ ZNaOvWlO'rtajCO'HjO но^о1>+гс-наг&-»гсог

С ост с Б гиова:

Ка-С0 -Ь5т707«"

Ко 2 5,-20+ 3Q7. Ho^SO,,- ! » 64» jfiaptwwciQ чериыД"" йИшпшйыЗ mt*0Kf кони,.60* 6S% испарит еаь J

Черныи UtflQK на «рорсуцци. кони,. 60» 65V,

Рис. 2. Оборот химикатов и щелоков в производстве сульфатной целлюлозы выпаривания. Ее термическое сопротивление теплопередаче достигает 50 % и более от общего [57, 239]. Соотношение между отдельными составляющими общего термического сопротивления теплообмену имеет большое значение для анализа влияния различных факторов на интенсивность теплопередачи. Эти соотношения постоянно меняются для всех корпусов ВС в зависимости от продолжительности их работы вследствие постепенного увеличения термического сопротивления накипи. Периодическая химическая очистка поверхности теплообмена выпарных аппаратов (ВА) обуславливает необходимость поочередного их выведения из технологического режима, что снижает производительность ВС. Стремление увеличить продолжительность работы ВА между промывками при относительно малых температурных напорах определяет применение больших поверхностей теплообмена ВА и их работу при относительно низких значениях тепловых нагрузок (5*15)-103 кВт. Малые полезные температурные разности 1ST = t?p — ввт и температурные напоры определяются не параметрами греющего пара, которые могут быть легко повышены, а ограничивается допустимой по технологическим условиям температурой кипения раствора. Ее превышение приводит к разложению щелоков с выделением летучих серосодержащих элементов. Поэтому локальные ухудшения теплообмена и, связанные с этим, повышения средних и локальных значений температур раствора, допустимых с точки зрения общей интенсивности теплопередачи, могут оказаться крайне нежелательными. Температура вторичного пара из последнего корпуса определяется технико-экономическими возможностями конденсатора. При этом полезная разность температур, приходящаяся на один корпус, составляет 7 -fl5°C.

В настоящее время практически на всех предприятиях сульфат-целлюлозного производства установлены ВС, укомплектованные пятью -девятью вертикальными длиннотрубными ВА. На отдельных станциях установлены концентраторы различных конструкций, в том числе с гравитационным течением выпариваемого раствора и пластинчатой поверхностью теплообмена. Размеры суммарной поверхности теплообмена ВС в зависимости от конструкций аппаратов и тепловых схем лежат в пределах от 10000 до 30000м2. К рациональным конструкциям (ВА) предъявляются высокие требования с точки зрения стоимости их изготовления, технологической надежности и интенсивности процесса теплопередачи. Удовлетворение этим требованиям возможно только при наличии надежных данных и закономерностей по теплообмену в ВА.

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

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

Для условий оптимального режима при кипении в трубах, рядом исследователей на основании теории подобия получены полуэмпирические критериальные зависимости. Анализ литературных данных [221] показал, что большинство работ по исследованию теплообмена при кипении жидкостей выполнено применительно к чистой воде. Часть работ проведена на растворах, главным образом сахарных. В обоих случаях в основу обработки данных положены закономерности кипения чистых жидкостей. Распространение зависимостей, полученных для конкретных сред и условий, на область кипения сульфатных щелоков и других подобных растворов не правомерно. В весьма ограниченном количестве работ по теплообмену при кипении отработанных пенных растворов ЦБП [227] отсутствуют достаточно надежные рекомендации для расчета и регулирования режимов работы выпарной аппаратуры ЦБП. На основе опубликованных данных вывести критерии моделирования ВА не представляется возможным, поэтому наиболее достоверные экспериментальные данные в настоящее время могут быть получены только при исследовании процессов выпаривания отработанных растворов ЦБП в трубах натуральных, по отношению к промышленным ВА, геометрических размеров и в рабочем диапазоне режимных параметров.

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

В соответствии с этим сформулированы основные цели и задачи работы:

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

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

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

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

5. Изучить процесс конденсации пара из парогазовой смеси на вертикальных трубах ВА и разработать на основе экспериментальных данных и их анализа расчетные зависимости для определения интенсивности теплообмена.

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

7. Создать программный метод теплового расчета выпарных станций.

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

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

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

1. СОСТОЯНИЕ И РАЗВИТИЕ ИССЛЕДОВАНИЙ ПО ИЗУЧЕНИЮ ТЕПЛОПЕРЕДАЧИ ПРИ ВЫПАРИВАНИИ ВОДЫ И РАСТВОРОВ В ТРУБАХ ВЫПАРНЫХ АППАРАТОВ

ОСНОВНЫЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ ПРОМЫШЛЕННРСТИ

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

На основе результатов проеденных исследований сделаны следующие выводы:

1. Предложены пути совершенствования энергоиспользования при применении противоточной схемы ВС на ЦБК. Экспериментально показана возможность и даны параметры безнакипного выпаривания пенных отработанных растворов ЦБП в длиннотрубных аппаратах с гравитационным течением выпариваемого раствора.

2. Разработаны средства, методики для экспериментального исследования процессов теплообмена, создан экспериментальный комплекс для проведения стендовых испытаний и проведены комплексные исследования теплопередачи при выпаривании восходящих и гравитационно-стекающих пенных растворов ЦБК на стендовых установках и промышленных ВС.

3. На основе экспериментальных данных разработаны: физическая модель движения многокомпонентных пенных двухфазных парорастворных потоков в кипятильных трубах ВА в условиях восходящего и гравитационного течения; определено влияние физико-химических характеристик черного сульфатного щелока, водного предгидролизата и режимных параметров работы ВА на интенсивность выпаривания.

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

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

6. Предложена методика расчета для определения плановых промывок

• ВС.

7. Экспериментально установлены особенности теплоотдачи к кипящим сильнопенящимся растворам в трубе: отмечено быстрое развитие кризиса теплоотдачи как при высоких (д < 100 квт/м2), так и при низких (д < 10 кВт/м2) тепловых потоках при определенных режимных параметрах; в качестве параметров, характеризующих предельный режим работы испарительных каналов с подъемным движением двухфазного потока, рекомендовано принимать предельное влагосодержание (1-<р)пр и критические тепловые нагрузки qKp; для каналов с гравитационно-стекающей пленкой жидкости режим начала ухудшения теплоотдачи характеризуется предельной плотностью орошения

Гпред.> предложены расчетные зависимости для определения величин (1~(р), а Якр.> ГПред.

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

9. Результаты исследований по динамике накипеобразования, приве-Ц* денные в настоящей работе в виде номограмм и расчетных зависимостей, являются основой для управления работой ВУ ЦБК, обеспечивают максимальную продолжительность их эксплуатации между стадиями промывки, в результате чего повышается производительность ВС и всего комбината.

10. Обобщенные методы расчета теплогидродинамических процессов, основные концепции комплексной оптимизации энергоемких технологических звеньев теплотехнологических схем, предложения по снижению накипеобразования на кипятильных поверхностях труб в виде программного метода расчета ВС используются при модернизации, разработке, внедрении и промышленной отладке ВС ЦБК, ЗАО «ГИПРОБУМ», ГУП «ЦНИИБум

• маш». Результаты исследований используются в учебном процессе СПбГТУРП.

1.Авксентюк Б.П. Кризис теплоотдачи первого рода при кипении жидкостей в трубах //Тезисы докладов Всесоюзной конференции: Теплообмен в парогенераторах.- Новосибирск, 1986.-С. 14-15.

2. Адлер Ю.П. Планирование эксперимента при поиске оптимальных условий Текст. /Ю.П. Адлер, Е.П. Маркова, Ю.В. Грановский. М.: Наука, 1976.- 279 с.

3. Аладьев И.Т. Теплоотдача при пузырчатом кипении Текст. / И.Т. Аладьев //Конвективный и лучистый теплообмен: сб.науч.тр.- М.: Изд. АН СССР, I960.- С. 233-255.

4. Алексеенко С. В. Трение при отекании пленки по вертикальной стенке Текст. / С.В. Алексеенко, В.Е. Накоряков, Б.Г. Покусаев, В.В. Христофоров //ИФЖ.- 1973. т. 24, № 5. С 824-830.

5. Аксененко М.Д. Приемники оптического излучения Текст.: Справочник / М.Д. Аксененко, М.Л. Бараночников.- М.: Радио и связь, 1987.-296с.

6. Андреев А. Ф. Об устойчивости ламинарного течения тонких слоев жидкости Текст. /А.Ф. Андреев //ЖТФ.- 1963. -т. 45, вып. ЗОТ.- С. 755-759.

7. Андреевский А.А. Волновое течение тонких слоев вязкой жидкости Текст. / А.А. Андреевский // Температурный режим и гидравлика парогенераторов.- М.: Изд-во АН СССР, 1978.- вып. 1.- С. 181-230.

8. Андреевский А.А. Охлаждение поверхности нагрева двухфазным пароводяным потоком Текст. / А.А. Андреевский //Конвективный теплообмен в элементах парогенераторов и теплообменников: сб. науч.тр. ЦКТИ, вып.86.- Л., 1968.- С. 3-25.

9. Н.Ардашев В.А. Исследование теплообмена при выпаривании гравитационно-стекающей пленки жидкости в вертикальных трубах Текст.: авто-реф. дис. . канд. техн. наук / Ардашев В.А.- Киев, 1983.

10. Архипов Г.А. Экспериментальное исследование теплоотдачи при конденсации водяного пара на вертикальной трубе Текст. / Г.А.Архипов, Г. Фенек // Расчет и конструирование энергооборудования с конденсацией пара: сб. научн. тр.- Л.: ЛКИ ,1990.- С. 131-159.

11. З.Байдаков В.Г. Кинетика вскипания газированной жидкости Текст. / В.Г. Байдаков, A.M. Каверин, Г.Ш. Болтачев : тез. докл. и сообщ. V Минского международного форума по тепло-и массообмену.- Минск, 2004 .- С. 6.

12. Бартоломей Г.Г. Истинные объемные паросодержания при кризисе теплообмена Текст. / Г.Г. Бартоломей, Ю.А. Дунаев, Ю.В. Харитонов // Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесозн. конф., т. 1.- J1, 1990.- С. 213-214.

13. Бартоломей Г.Г. Истинное объемное паросодержание при кипении с недогревом для различных законов тепловыделения по длине канала Текст. /Г.Г. Бартоломей, J1.C. Саботинов //Ядерная энергия .- 1976.- №3.-С. 33-41.

14. Бартоломей Г.Г. Экспериментальное исследование истинного объемного паросодержания при малых массовых скоростях Текст. / Г.Г. Бартоломей // Теплообмен в парогенераторах : тез. докл. Всесозн. конф. — Новосибирск, 1986.- С. 250-251.

15. Барулин Ю.Д. Исследование теплогидравлических характеристик труб применительно к работе парогенераторов ТЭС Текст. / Ю.Д. Барулин //Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесоюзн. конф., т.1.- Л., 1990.- С. 13-15.

16. Беглес А. Исследование режимов течения кипящей воды Текст. / А. Беглес, Н. Сю // Достижения в области теплообмена: кн. / под ред. В.М. Боришанского.- М.: Мир, 1970.- С. 30.

17. Безродный М.К. Кризисы теплопереноса в замкнутых двухфазных термосифонах Текст.: автореф. дис. . докт. техн. наук/ Безродный М.К.Киев, КПИ, 1984.

18. Безродный М.К. Влияние скорости пара на толщину пленки жидкости при спутном восходящем течении Текст. / М.К Безродный, Ю.В. Ан-тошко : тр. Первой Российской нац. конф. по теплообмену, т. 6.-М.: МЭИ, 1994.- С. 35-41.

19. Берман Л.Д. О критериях подобия для совместно протекающих процессов тепло- и массообмена в гетерогенных системах Текст. /Л.Д. Берман //ЖТФ.- 1958, т. 28, № 11.-С. 2617-2629.

20. Берман Л.Д. К определению коэффициента массоотдачи при расчете конденсации пара, содержащего примесь воздуха Текст. / Л.Д. Берман //Теплоэнергетика.- 1969.- № 10.- С. 68-71.

21. Блюхин В. А. Исследование гидродинамики и тепломассообмена нестационарных течений жидких пленок Текст. / В.А. Блюхин, Л.А. Кали-муляна // Тепломассообмен-5.- Минск, 1976, т. 4.- С. 33-40.

22. Бобе Л.С. Тепло и массообмен при конденсации пара из парогазовой смеси при турбулентном течении внутри трубы Текст. / JI.C. Бобе, В.А. Солоухин // Теплоэнергетика.- 1972.- № 9,- С. 27-30.

23. Бобе JI.C. К расчету конденсации пара при поперечном омывании труб парогазовой смесью Текст. /JI.C. Бобе, Д.Д. Малышев // Теплоэнергетика.- 1971.- № 12.- С.84-86.

24. Бойков В.Ю. Теплообмен при выпаривании в стекающей пленке черного щелока сульфат-целлюлозного производства Текст. : автореф. дис. . канд. техн.наук / Бойков Вячеслав Юрьевич.- СПб., 1999.-15с.

25. Боришанский В.М. Теплоотдача к двухфазному потоку Текст. / В.М. Боришанский //Теплоэнергетика .- 1969.-№ 5.- С. 58-61.

26. Боришанский В.М. Теоретическое обоснование теплового расчета парогенераторов в докризисном режиме Текст. / В.М. Боришанский : сб. на-уч.тр. ЦКТИ, вып. 108.-М., 1971.-С. 25.

27. Бошняк Л.Л. Измерения при теплотехнических исследованиях Текст. / Л.Л. Бошняк.- Л.: Машиностроение, 1974.- 448 с.

28. Буглаев В.Т. Теплоотдача при поперечном обтекании горизонтального трубного пучка конденсирующимся паром Текст. /В.Т.Буглаев, М.М Андреев //Энергомашиностроение.- 1973.-№ 1.- С. 36-38.

29. Буглаев В.Т. Теплообмен при конденсации пара атмосферного давления на горизонтальном трубном пучке Текст. /В.Т. Буглаев, М.М. Андреев, М.М. Казаков //Теплоэнергетика.- 1970.- № 11.- С. 66-67.

30. Веденеев В.А. Исследование теплообмена при вынужденном движении растворов в трубах Текст. /В.А. Веденеев // Материалы научно-технической конференции, вып. 2.-JI.: ЛТИЦБП, 1974.-С. 132-134.

31. Вишнев И.П. Обобщенные зависимости для расчета теплоотдачи при кипении жидкостей Текст. /И.П. Вишнев //Химическое и нефтяное машиностроение.-1978.- № 7.- С. 17-19.

32. Воинов Н.А. Теплоотдача в пленке жидкости, стекающей по гладкой и шероховатой поверхности Текст. / Н.А. Воинов, А.Н. Николаев, Н.А. Николаев // ТОХТ.- 1998.- № 1, т.32.- С. 28-32.

33. Воронцов Е.Г. Теплообмен в жидкостных пленках Текст. / Е.Г. Воронцов, Ю.М. Тананайко. Киев, Техника.- 1972.- 194 е., ил.

34. Воронцов Е.Г. О минимальной плотности орошения вертикальных пленочных аппаратов Текст. / Е.Г. Воронцов // ИФЖ. -1968.- № 6, т. 14.- С. 1075-1078.

35. Воронцов Е.Г. Минимальная плотность орошения и равномерность орошения при эксцентричном щелевом распределении жидкости в пленку Текст. / Е.Г. Воронцов, Е.Д. Малимов // ЖПХ 47.- 1974.- № 7.- С. 15371542.

36. Гаврилин А.В. Исследование теплообмена при сгущении молока в тонком слое Текст. / А.В. Гаврилин // Молочная промышленность. -1968.-№ 8.-С. 28-31.

37. Ганчев Б.Г. Расчет локальных значений средней толщины турбулентной пленки жидкости, стекающей по вертикальной поверхности Текст./ Б.Г. Ганчев, В.М Козлов, В.В. Лозовецкий // Изв. вузов СССР. Машиностроение.- 1970.-№ 1.-С. 112-116.

38. Ганчев Б.Г. Стекание пленки жидкости в вертикальном канале Текст. / Б.Г. Ганчев, В.М. Козлов, В.В. Лозовецкий // Исследование процессов в энергетических установках: сб. научн. тр. МВТУ, вып.2 1975.- № 207.- С. 40-45.

39. Ганчев Б.Г. Исследование предельных условий пленочного охлаждения вертикальной поверхности Текст. / Б.Г. Ганчев, А.Е. Боков // Судовые энергетические установки.-Владивосток, 1981.-С. 100-108.

40. Ганчев Б.Г. Исследование гравитационного течения пленки жидкости по стенкам вертикального канала большой длины Текст. / Б.Г. Ганчев, В.М. Козлов//ПМТФ.- 1973. -№ 1. -С. 128-135.

41. Ганчев Б.Г. Исследование термокапиллярной устойчивости при гравитационном стекании пленки жидкости.Текст. / Б.Г. Ганчев, А.Е. Боков //ИФЖ.- 1980. № 4, Т.39.- С. 581-591.

42. Ганчев Б.Г. Исследование местных коэффициентов теплоотдачи в пленке жидкости, стекающей по вертикальной поверхности. Текст. / Б.Г. Ганчев, В.М. Козлов, В.В. Лозовецкий, В.М. Никитин // Изв. вузов СССР. Машиностроение.- 1970.- № 9.- С. 114-117.

43. Ганчев Б.Г. Исследование скоростей в стекающих пленках жидкости в условиях развивающегося волнового движения. Текст. / Б.Г. Ганчев, В.М. Козлов // Исследование процессов в энергетических установках: сб.науч.тр. МВТУ.- 1975.- № 207, вып. 2.-С.2-61.

44. Гельперин Н.И. Рациональная форма корреляции коэффициентов теплоотдачи к органическим жидкостям, кипящим при атмосферном давлении в вертикальных трубчатых аппаратах Текст. / Н.И. Гельперин // Химическое и нефтяное машиностроение.-1972.-№6.-С.9-11.

45. Гимбутис Г.И.Теплообмен при гравитационном течении пленки жидкости Текст. / Г.И. Гимбутис. Вильнюс: Мокслас, 1988.-233 е.: ил.

46. Гимбутис Г.И. Теплоотдача при турбулентном течении гравитационной пленки жидкости на вертикальной поверхности Текст. / Г.И. Гимбутис // Теплообмен. -1974:сб. науч.тр.- М.: 1975.-С. 44-50.

47. Гимбутис Г.И.Теплоотдача при поверхностном кипении гравитационной турбулентной пленки воды Текст. / Г.И. Гимбутис, А.Ю. Дробавичус //Механика-7. :сб .науч.тр.-Каунас: 1977.-С.90-92.

48. Глоба В.З. Разработка выпарной установки мыльно-щелочных растворов Текст./ автореф. дис. . канд.техн. наук / Глоба Виктор Зиновьевич-Л., 1985.-15с.:ил.

49. Гогонин И.И. К обобщению опытных данных по критическим тепловым потокам в стекающих пленках жидкости Текст. / И.И. Гогонин, А.Р. Дорохов // Изв. СО АН СССР, сер. техн. наук. -1980.- № 8 , вып.2.-С. 100102.

50. Гогонин И.И. К вопросу образования «сухих пятен» в стекающих тонких пленках жидкости Текст. / И.И. Гогонин, А.Р. Дорохов, В.Н. Бога-чев //Изв. СО АН СССР, сер. техн. наук.- 1977.- № 13, вып З.-С. 46-51.

51. Гогонин И.И. Расчет локального тепломассообмена при конденсации на пакетах труб Текст. / И.И. Гогонин //Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесоюзн.конф. , т.2,- Л., 1990.-С. 134136.

52. Горбовский Б.Г. Устройство и обслуживание выпарных станций сульфат-целлюлозного производства Текст.: пособие для мастеров и выпарщиков / Б.Г. Горбовский М.; Л.: Гослесбумиздат,1953.- 191с.:ил.

53. Горянов Ю.А. Гидродинамика вспененных потоков в трубах различной ориентации и гибах Текст. / Ю.А. Горянов //Двухфазный поток в энергетических машинах и аппаратах: тез. докл.8-й Всесоюзн. конф., т.2.- Л., 1990.-С. 191.

54. Гребер Г.Основы учения о теплообмене Текст. /Г. Гребер, С. Эрк, У. Григулль .-М.: Иностран. лит., 1958.- 566 е.: ил.

55. Грибаненков А.В. Исследование теплоотдачи при кипении многокомпонентных и бинарных водных растворов в вертикальной трубке в условиях вынужденного движения Текст. :автореф. дис. . канд. техн. наук/ Грибаненков А.В.- М.: МИХМ, 1970.

56. Дас Каперу Р. Экспериментальное исследование гидродинамики и теплоотдачи в стекающих пленках воды, водных растворов солей и поверхностно-активных веществ Текст. : автореф. дис. канд.физ. наук /Дас Каперу Р. Киев, 1970.

57. Демехин Е. А. О нестационарных волнах в слое вязкой жидкости Текст. / Е.А. Демехин, В.Я. Цикадов // Изв. АН СССР. Механика жидкости и газа, №3.- 1981.-С. 151-154.

58. Джиджи Л.М. Пузырьковый пограничный слой и профили температуры при кипении жидкости, движущейся принудительно в канале Текст. / Л.М. Джиджи, И.А. Кларк // Теплопередача .- М.:Мир, 1969.- № 1.-С.66.

59. Доманский И.В. Определение режимов устойчивой работы выпарных аппаратов с падающей пленкой Текст. / И.В. Доманский, В.Н. Соколов //ЖПХ, т.40, вып.2.-1967.-С. 365-370.

60. Доманский И.В. Теплоотдача к падающей пленке жидкости, предварительно нагретой до температуры кипения Текст. / И.В. Доманский, В.Н. Соколов//ЖПХ, т.40, вып.1.- 1967.- С. 66-71.

61. Дорощук В.Е. Кризисы теплообмена при кипении воды в трубах Текст. / В.Е.Дорощук. М.: Энергоатомиздат, 1983.- 167 с.

62. Дрейцер Г.А. Анализ современного состояния теплообмена при конденсации теплоносителей Текст. / Г.А. Дрейцер // Двухфазный поток в энергетических машинах и аппаратах : тез. докл. 8-й Всесоюзн. конф., т.2.- Л., 1990.-С. 140-142.

63. Дымов М.И Теплообмен в тонком слое жидкости при регулярно-волновом режиме течения Текст. / М.И. Дымов, В.Х. Кириллов, И.Л. Чертков // Холодильная техника и технология. -1975.- № 21.-С. 50-54.

64. Дымов М.И. Устойчивость течения тонкого жидкого слоя в теплообменниках пленочного типа Текст. / М.И. Дымов, В.Х. Кириллов // Холодильная техника и технология.- 1977.- № 25.-С. 59-63.

65. Жмай Л.А. Минимальная плотность орошения в пленочных трубчатых аппаратах Текст./ Л.А. Жмай , В.М. Олевский //Химия и технология продуктов органического синтеза. Процессы и аппараты, вып.1., ч.2: сб.науч.тр. ГИАП.- М., ОНТИ, 1969.- С. 5-12

66. Живайкин Л.Я. Течение пленки на вертикальной поверхности Текст. / Л.Я. Живайкин, Б.В. Волгин //ЖПХ, т. 34, вып. 6.- 1961.-С. 1236-1242.

67. Живайкин Л.Я. Теплообмен между вертикальной стенкой и стекающей пленкой жидкости при ее охлаждении Текст. / Л.Я. Живайкин, В.И. Холостых, И.Г. Бляхер // Изв. ВУЗов СССР. Энергетика.- 1978.- № 4.-С. 132135.

68. Живайкин JI.Я. Исследование гидродинамики и теплообмена при пленочном течении жидкости Текст. / Л.Я. Живайкин, А.А. Шехтман , И.Г. Бляхер , В.И. Холостых: тр. Уральского научн. исслед. хим. ин-та, вып.41.-1976.-С. 47-52.

69. Иванов О.П. Исследование теплообмена при кипении смесей фреон-фреон, фреон-масло Текст.: автореф. дис. . канд.техн. наук / Иванов Олег Петрович. Л., ЛТИХП, 1966.

70. Иванов О.П. Исследование и интенсификация теплообмена в конденсаторах холодильных машин Текст. : автореф. дис. . докт. техн. наук / Иванов Олег Петрович. — Л., ЛТИХП, 1975.

71. Ивановский М.Н. Особенности пленочной и капельной конденсации паров металлов Текст. / М.Н. Ивановский // Тепло- массоперенос при физико-химических превращениях / Ивановский М.Н., Милованов Ю.В., Сорокин В.П.- Минск: Наука и техника, 1968.- С. 387-394.

72. Иванский А.Л. О нелинейных волнах на вертикальной пленке жидкости Текст./ А.Л. Иванский // ПМТФ.- 1980.- № 2.-С. 52-58.

73. Ивашкевич А.А. О граничном паросодержании при кипении воды в трубах Текст. / А.А.Ивашкевич //Двухфазный поток в энергетических машинах и аппаратах, т. 1:тез. докл. 8-й Всесоюзн.конф.- Л., 1990.-С. 402-404.

74. ЮО.Ивашкевич А.А. Обобщение опытных данных по кризису кипения при вынужденном течении воды в трубах Текст. / А.А.Ивашкевич //Кипение, кризисы кипения, закризисный теплообмен, т.4: тр. Первой Рос. нац. конф. по теплообмену.- М., МЭИ, 1994.- С. 134-138.

75. Ивашкевич А.А. К вопросу обобщения опытных данных по кризису кипения воды в трубах Текст. / А.А. Ивашкевич // Кипение, кризисы кипения, закризисный теплообмен, т 4: тр. Второй Рос.нац. конф. по теплообмену.- М., МЭИ, 1998.- С. 145-148.

76. Ильин И.Н. Начало кипения жидкостей Текст. / И.Н. Ильин, В.П. Гривцов, У.О. Саука // Кипение и конденсация: сб. науч.тр. Рига, РПИ, 1981.- С. 93-117.

77. Исаченко В.П. Теплопередача Текст. / В.П. Исаченко, В.А. Осипо-ва, А.С. Сукомел. М.: Энергия, 1975.- 488 с.

78. Исаченко В.П. Теплообмен при конденсации Текст. / В.П. Исаченко. М.: Энергия, 1977.- 239 е.: ил.

79. Исаченко В.П. Теплоотдача при испарении воды из пористой стенки, омываемой воздухом Текст. / В.П. Исаченко, В.В. Взоров, В.А. Верто-градский // Теплоэнергетика. -1961.- №1. С. 65-72.

80. Юб.Ионайтис P.P. Патентно-техническое исследование способов образования пленки жидкости Текст. / P.P. Ионайтис, H.JI. Шведов // Атомная техника за рубежом. -1980.- № 1.- С. 8-15.

81. Ю7.Кабаньков О.Н. Исследование теплообмена и границ устойчивости стекающих пленок жидкости по гладкой и структурированной поверхности Труды. / О.Н. Кабаньков, JI.A. Сукомел , В.В. Ягов // Двухфазные течения, т. 6 .- М.: МЭИ, 1994.- С. 63-69.

82. Ю9.Калинин Н.Н. Химическая технология древесины. Расчет выпарной станции Текст. / Н.Н. Калинин Д.: JITA, 1971.- 19 с.

83. О.Капица JI.JL Волновое течение тонких слоев вязкой жидкости Текст. / Л.Л. Капица//ЖЭТФ, т. 18, вып.1, 1948.- С. 3-28.

84. Ш.Капица Л.Л. Волновое течение тонких слоев вязкой жидкости Текст. / Л.Л. Капица, С.П. Капица // ЖЭТФ, т. 19, вып. 2, 1949.-С. 105-120.

85. Кашинский О.Н. Структура восходящего снарядного течения в вертикальной трубе Текст. / О.Н. Кашинский, Р.С. Горезин , В.В. Рандин : тр. 3-й Рос. нац. конф. по теплообмену,т. 4.- М., МЭИ, 2002.-С. 68-72.

86. Кашинский О.Н. Эволюция локальных характеристик восходящего снарядного течения в вертикальной трубе Текст. / О.Н. Кашинский, В.В.

87. Рандин, А.С. Курдюмов // V Минский международный форум по тепло-и массообмену: тез. докл. и сообщ.-. Минск, 2004.- С. 43.

88. Пб.Квурт Ю. П. Влияние частоты переменного тока на результаты измерения толщины пленки жидкости методом электропроводности Текст. / Ю.П. Квурт, JI.M. Холпанов, В.А. Малюсов, Н.М. Жаворонков // ЖПХ.-1981.- № 5, т. 54.-С. 1068-1072.

89. П.Кириллов В.Х. Волновое течение тонкого слоя вязкой жидкости по наклонной поверхности Текст. / В.Х. Кириллов, В.И. Логачевский // Холодильная техника и технология.- Киев, 1982.- №3.- С. 117-122.

90. Кириллов П.Л. Справочник по теплогидравлическим расчетам (ядерные реакторы, теплообменники, парогенераторы) Текст. / П.Л. Кириллов, Ю.С. Юрьев, В.П. Бобков.- М.: Энергоатомиздат, 1984.- 294 с.

91. Козлов В.М. Исследование гравитационного течения пленки жидкости методом нейтронной диагностики Текст. / В.М. Козлов, В.В. Гусев, М.Г. Месропов // ТОХТ.- 1976.- № 1,т.10.-С. 69- 73.

92. Копп И.З. Влияние реальной структуры поверхности теплообмена на генерацию зародышей паровой фазы Текст./ И.З. Копп // Двухфазный поток в энергетических машинах и аппаратах: тез.докл. 8-й Всесоюзн. конф. Т.1.-Л., 1990.-С. 265-266.

93. Комшилов Н.Ф. Сульфатный черный щелок и его использование Текст. /Н.Ф. Комшилов.- М.:Лесная промышленность, 1969,- 184 с.

94. Кремлевский П.П. Расходомеры и счетчики количества Текст. /П.П. Кремлевский.-Л.: Машиностроение, 1975.- 776 с.

95. Кремлевский П.П. Расходомеры Текст. / П.П. Кремлевский.- Л.: Машгиз, 1963.- 656 с.

96. Кружилин Г.Н. Уточнение нуссельтовской теории теплообмена при конденсации Текст. / Г.Н. Кружилин // ЖТФ.- 1937. -Т. 7, вып. 20/21.- с. 2017-2022.

97. Кузма-Кичта Ю.А.Исследование ухудшения теплообмена при кипении трехкомпонентных водных растворов при низких массовых скоростях и давлениях Текст. / Ю.А. Кузма-Кичта: тр. 3-й Рос. нац. конф. по теплообмену, т. 4.- М., МЭИ, 2002.- С. 123-127.

98. Кулов Н.Н. Массоотдача в стекающих пленках жидкости Текст. / Н.Н. Кулов, В.В. Максимов, В.А. Малюсов, Н.М. Жаворонков // ТОХТ. -1983.-№3,т. 17.-С. 291-306.

99. Кутателадзе С.С. Теплопередача при конденсации и кипении Текст. / С.С. Кутателадзе .-М.; Л.: Машгиз, 1952.- 231с.: ил.

100. Кутателадзе С.С. Тепломассообмен и трение в турбулентном пограничном слое Текст. / С.С. Кутателадзе, А.И. Леонтьев.-2-е изд., перераб.-М.: Энергоатомиздат, 1985.- 320 с.

101. Кутателадзе С.С. Влияние скорости циркуляции на коэффициент теплоотдачи при кипении в трубах Текст. / С.С. Кутателадзе //Энергомашиностроение.-1961.-№1.- С. 12-15.

102. Кутателадзе С.С. Тепломассообмен и волны в газожидкостных системах Текст. / С.С. Кутателадзе, В.Е. Накоряков. Новосибирск: Наука, 1984.-301 е.: ил.

103. Кутателадзе С.С. Гидродинамика газожидкостных систем Текст. / С.С. Кутателадзе, М.А. Стырикович .-М.: Энергия, 1976.- 296 е.: ил.

104. Кутателадзе С.С. Основы теории теплообмена Текст. / С.С. Кутателадзе. Новосибирск: Наука, 1970.- 658 е.: ил.

105. Кутепов А.М.Гидродинамика и теплообмен при парообразовании Текст.: учеб. пособие для ВТУЗов/ A.M. Кутепов, Л.С. Стерман, Н.Г. Стю-шин. -3-е изд., испр. -М.: Высшая школа, 1986.- 352 с.

106. Лабунцов Д.А. Обобщенные зависимости для теплоотдачи при пузырьковом кипении жидкостей Текст. /Д.А. Лабунцов // Теплоэнергетика.-№5.- 1960.- С. 76-81.

107. Лабунцов Д.А. Исследование паросодержаний неравновесных двухфазных потоков Текст. / Д.А. Лабунцов// Сб. трудов НИЭИ им. Г.М. Кржижановского, вып. 35, 1976.-С. 88-98.

108. Лабунцов Д.А. Теплоотдача при пленочной конденсации чистых паров на вертикальных поверхностях и горизонтальных трубах Текст. / Д.А. Лабунцов // Теплоэнергетика.- 1957.- № 7.-С. 72-80.

109. Лабунцов Д.А. Теплообмен при конденсации пара на вертикальной поверхности в условиях турбулентного стекания пленки конденсата Текст. / Д.А. Лабунцов // ИФЖ.- I960.- № 8,т.З.-С. 3-12.

110. Лабунцов Д.А. Вопросы теплообмена при пузырьковом кипении жидкости Текст. / Д.А. Лабунцов //Теплоэнергетика.- 1972.- № 9.- С. 14-19.

111. Лабунцов Д.А. О влиянии конвективного переноса тепла и сил инерции на теплообмен при ламинарном течении конденсатной пленки Текст. / Д.А. Лабунцов // Теплоэнергетика.- 1956.- № 12.-С. 47-50.

112. Ландау Л.Д. Теоретическая физика Текст.: учеб. пособие. В 10 т. Т. VI. Гидродинамика /Л.Д. Ландау, Е.М. Лившиц. М.: Наука, 1988. - 736 с.

113. ИЗ.Левеншпиль О. Инженерное оформление химических процессов Текст. / О. Левеншпиль.- М.:Химия,1969.- 624 с.

114. Липсман B.C. Исследование конвективной теплоотдачи при поступательном и закрученном пленочном течении воды и сахарных раствороввнутри вертикальных труб Текст.: автореф. дис. . канд. техн. наук / Липс-ман B.C.- Киев, КТИПП ,1972.

115. Лобачев А.Г.Исследование истинных объемных паросодержаний в обогреваемой трубе при подъемном и опускном движении двухфазного потока Текст. /А.Г. Лобачев// Теплоэнергетика,- 1973.-№ 5.-С. 75-77.

116. Мальцев Е.Д. Исследование гидродинамики и теплообмена при пленочном течении морской воды Текст. / Е.Д. Мальцев, B.C. Тимофеев // Сб. трудов МИСИ.-М.: МИСИ.- 1972.-№ 89.-С. 159-167.

117. Милашенко В.И. Исследование гидродинамических характеристик дисперсно-кольцевых парожидкостных потоков Текст.: автореф. дис. . канд. техн. наук / Милашенко В.И. М.,1977.

118. Минухин Л.А. Тепломассообмен при конденсации пара из парогазовой смеси, поперечно обтекающей вертикальные трубки Текст. /Л.А.Минухин, С.В. Григоренко // Кипение и конденсация: сб. науч. тр. РПИ.- Рига, 1981.-С. 83-93.

119. Миропольский З.Л. Объемные паросодержания при опускном движении пароводяной смеси в трубах и межтрубном пространстве Текст. /З.Л. Миропольский, Р.И. Шнеерова, А.И. Карамышева//Теплоэнергетика.- 1978.-№4.-С. 60-62.

120. Миропольский З.Л. Теплоотдача при конденсации перегретого и насыщенного пара внутри труб Текст. / З.Л. Миропольский, Р.И. Шнеерова // Двухфазный поток в энергетических машинах и аппаратах : тез. докл.8-й Всесоюзн. конф.Т.2- Л., 1990.-С. 49-51.

121. Митрофанов С.М. Моделирование процесса взрывного пристеночного парообразования Текст. / С.М. Митрофанов, П.А. Павлов // Труды 3-й Рос. нац. конф. по теплообмену. Т. 4.- М., МЭИ, 2002.-С. 144-147.

122. Михеев М.А. Основы теплопередачи Текст. / М.А. Михеев.- М.-Л.: Госэнергоиздат, 1956.- 392 е.: ил.

123. Мовсесян В.Л. Выпарные аппараты и схемы выпарных станций ЦБП Текст.: обзорн. информ. ЦИНТИхимнефтемаш.- М., 1979.- 50 с.

124. Исследование механизма теплообмена при пленочном кипении сульфатного щелока с целью выявления путей повышения производительности выпарных аппаратов Текст.: отчет о НИР: 64.2 / НИС ЛТИ ЦБП.- Л., 1981.- Исполн.: Мовсесян В.Л., Суслов В.А.-ГР № 81052256.

125. Исследование процесса выпаривания водных предгидролизатов Текст.: отчет о НИР:62.8, разд.8 / НИС ЛТИ ЦБП.- Л., 1984,-Исполн.: Мов-сесян В.Л., Суслов В.А.- ГР № 01.890091732.

126. Разработка рекомендаций по интенсификации и стабилизации режимов работы выпарных установок сульфат-целлюлозного производства Текст.: отчет о НИР 316 (доп. согл. № 1) / НИС ЛТИ ЦБП.-Л, 1987.- Ис-полн.: Мовсесян В.Л., Суслов В.А.- ГР № 01.0017821.

127. Молочников Ю.С. Обобщение экспериментальных данных по истинным объемным паросодержаниям при кипении воды с недогревом Текст. / Ю.С. Молочников, Г.Н. Баташова, В.Н. Михайлов, В.А. Солодкий // Теплоэнергетика.- 1982.-№ 7.-С. 47-50.

128. Моркущин С.А. Параметры волнового течения тонкого слоя вязкой жидкости по вертикальной плоскости. Текст. / С.А. Моркущин, А.Г. Шейк-ман // Гидродинамика и теплообмен в энергетических установках: сб. на-уч.тр.- Свердловск, 1975.-С. 22-25.

129. Мурзич А.Ф. Промышленные тепломассообменные процессы и установки в ЦБП Текст.: учеб. пособие / А.Ф. Мурзич, В.Л. Мовсесян, В.А. Суслов, Е.К. Гусев.-Л.: ЛТИ ЦБП, 1989.- 80 с.

130. Муратова Т.М. Кинетический анализ процессов испарения и конденсации Текст. / Т.М. Муратова, Д.А. Лабунцов // ТВТ АН СССР, т. 5, вып. 5-6: сб.науч .тр.- 1969.-С. 959-967.

131. Накоряков В.Б. Стационарные двумерные катящиеся волны на вертикальной пленке жидкости Текст. / В.Б. Накоряков, Б.Г. Покусаев, С.В. Алексеенко //ИФЖ.-1976.- № 5,т.30.-С. 780-785

132. Накоряков В.Е. Экспериментальное исследование неизотермической абсорбции стекающей пленкой Текст. / В.Е.Накоряков, А.П. Бурдуков, Н.С. Буфетов, Н.И. Григорьева, А.Р. Дорохов // ТОХТ.- 1980.- № 5,т.14.-С. 755-758.

133. Накоряков В.Е. Мгновенный профиль скорости в волновой пленке жидкости Текст. / В.Е. Накоряков, Б.Г. Покусаев, С.В. Алексеенко, В.В. Орлов//ИФЖ.- 1977.-№ 33,Т.З.-С. 399-404.

134. Невструева Е.И. Тепло-и массообмен в атомных энергетических установках с водоохлаждающими реакторами Текст. / Е.И. Невструева .-М., 1978.-112с.-( Итоги науки и техники.Тепло- и массообмен. т.1)

135. Невструева Е.И. Взаимосвязь тепловых и гидродинамических характеристик в двухфазном неравновесном потоке Текст. / Е.И. Невструева, В.В. Тютяев// Тепломассообмен V: сб.науч.тр.Т. 3, ч. 2.- Минск, 1976.-С. 13-20.

136. Непенин Ю.Н. Технология целлюлозы Текст. Т. II. : учеб. пособие для вузов / Ю.Н. Непенин. М.: Лесная промышленность, 1990.- 600 с.

137. Непенин Ю.Н. Технология целлюлозы Текст. Т. II: учеб. пособие для вузов/ Ю.Н. Непенин .- М.:Гослесбумиздат, 1963.- 935 с.

138. Несис Е.И. Основные проблемы теории теплообмена при кипении Текст. /Е.И. Несис //Двухфазный поток в энергетических машинах и аппаратах : тез. докл. 8-й Всесоюзн. конф. Т.1.- Л., 1990.-С. 233-234.

139. Нигматулин Б.И. Динамика многофазных сред Т. 2 Текст. / Б.И. Нигматулин.- М.:Наука, 1987.- 359 с.

140. Нигматулин Б.И. Исследование распределения жидкости между ядром и пленкой в дисперсно-кольцевом потоке Текст. / Б.И. Нигматулин,

141. B.И. Милашенко, Ю.З. Шугаев //Теплоэнергетика.- 1976.- №5.-С. 77-79.

142. Джон Г. Перри. Справочник инженера-химика.Т. 1 Текст. / Джон Г. Перри. Л.:Химия, 1969.- 639 с.

143. C. Шпак, В.Ф. Юхименко, Ю.В. Якубовский // Кипение и конденсация: сб. науч. тр.- Рига,1986.-С. 46-54.

144. Поляков Ю.А. Элементарный состав и поверхностное натяжение черного сульфатного щелока в процессе выпаривания Текст. / Ю.А. Поляков, А.Б. Маршак //Труды ЛТИЦБП, вып. 27.- Л., 1970.-С. 21-27.

145. Преображенский В.П. Теплотехнические измерения и приборы Текст. / В.П. Преображенский. М.- Л.: Госэнергоиздат, 1953.- 384 с.

146. Промышленная теплоэнергетика и теплотехника Текст.: Справоч-ник/А.М. Бакластов, В.М. Бродянский, Б.П. Голубев и др.; Под общ. ред. В.А. Григорьева и В.М. Зорина.- М.: Энергоатомиздат, 1983.- 552 е.: ил.

147. Прядко Н.А. Совершенствование теплообменных аппаратов и установок пищевой промышленности Текст.: автореф. дис. .докт. техн. наук / Прядко Н.А. -Киев, КТИПП, 1986.

148. Ракицкий В.А. Гидродинамика и теплообмен при пленочном течении жидкости по гладкой гидрофобной поверхности Текст.: автореф. дис. . канд. техн. наук / Ракицкий В.А.-. Киев, 1990.

149. Рассохин Н.Г. Расчет теплоотдачи при кипении Текст. / Н.Г. Рассохин //Теплоэнергетика.- 1970. № 9.-С. 58-61.

150. Рассохин Н.Г. Парогенераторные установки атомных электростанций Текст. / Н.Г. Рассохин .- М.:Атомиздат, 1980.- 360 с.

151. Раушер Миле Денни. Экспериментальное исследование пленочной конденсации при обтекании горизонтальной трубы нисходящим потоком па-ро-воздушной смеси Текст. / Миле Дени Раушер //Теплопередача .- 1974,- № 1, т. 96.-С. 86-92.

152. Ривкин C.JI. Теплофизические свойства воды и водяного пара Текст. / C.JI. Ривкин, А.А. Александров.- М.: Энергия, 1980.- 424 с.

153. Роговая И.А. Влияние физических свойств некоторых жидкостей на параметры волнового течения пленки Текст. / И.А.Роговая, В.М. Олевский, Н.Г. Рунова // ТОХТ.- 1973.- № 5,т.7-С. 734-739.

154. Розенов В.М. Теплообмен при кипении Текст. / В.М. Розенов //Современные проблемы теплообмен: сб. науч. тр.- М.; JL,Энергия.- 1966.-С. 212-260.

155. А.С. 1432122 СССР. Устройство для регенерации тепла и химикатов и парогазовых выбросов растворителя плава содорегенерационного котлоаг-регата Текст. / Л.В. Романова, А.В. Наумов, В.А. Суслов и др.(СССР).- 1988.

156. Романова Л.В. Установка для конденсации парогазовой смеси из бака плава содорегенерационных котлоагрегатов Текст. / Л.В. Романова, В.А. Суслов // V Минский международный форум по тепло-и массообмену : тез. докл. и сообщ.- Минск, 2004.-С. 119-120.

157. Сагань И.И. Исследование течения двухфазных потоков в вертикальной круглой трубе Текст. / И.И. Сагань, С.И. Ткаченко, Н.Ю. Тобилевич, Ю.К. Пинчук // Изв. вузов СССР. Пищевая технология.- 1970.-№ 3.-С. 162-165.

158. Синкевич А.Е. Обобщенный метод вычисления коэффициента теплоотдачи при пленочной конденсации пара в вертикальном канале Текст. /

159. A.Е Синкевич //Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесоюзн. конф. Т.1.- Л., 1990.-С. 67-69.

160. Слесаренко В. Н. Дистилляционные опреснительные установки Текст. / В.Н. Слесаренко.- М.: Энергия, 1980.- 244 е.: ил.

161. Слесаренко В.Н. Исследование теплообмена при термической дистилляции морской воды в тонкопленочных опреснительных установках Текст. / В.Н. Слесаренко // Тепломассообмен V: сб.науч.тр. Т. 3, ч. 1.-Минск, 1976.-С. 372-374.

162. Слесаренко В.Н. Исследование преснения морской воды Текст. /

163. B.Н. Слесаренко, Ю.В. Якубовский, Б.Я. Карастлев// Труды ДВПИ, т.71, вып. 2.- 1971.-С. 15-25

164. Смолин В.Н. Истинное объемное паросодержание в парогенери-рующих каналах Текст. / В.Н. Смолин // Двухфазные течения : тр. 1-й Рос. нац. конф. по теплообмену. Т.6.-М., МЭИ.- 1994.-С. 216-225.

165. Соловьева Г.И. Влияние конструкции исходного распределителя жидкости на гидродинамику в аппаратах со стекающей пленкой Текст.: ав-тореф. дис. . канд. техн. наук /Соловьева Г.И.- М., МИХМ.- 1984.

166. Соколов В.Н. Газожидкостные реакторы Текст. / В.Н. Соколов, ИВ. Доманский,- Л.: Машиностроение, 1976.- 214 е.: ил.

167. Сорокин Д.Н. К методике измерения локального паросодержания при кипении в большом объеме Текст. / Д.Н. Сорокин: препринт.- Об-нинск:ФЭИ-626,1975.- 19 с.

168. Сорокин А.Ф. Тепло-и массообмен при кипении растворов в выпарных установках Текст.: автореф. дис. . докт. техн. наук / Сорокин А.Ф.-М., МЭИ.- 1956.

169. Справочник по теплообменникам Текст.В 2 т. Т.1/ Под ред. Б.С. Петухова, B.C. Шикова ; пер. с англ. М.: Энергоатомиздат, 1987. - 560 е.: ил.

170. Стерман Л.С. К теории теплоотдачи при кипении жидкости Текст. / Л.С. Стерман //ЖТХ, т. XXIII, вып. 2,- 1953.-С. 341-351.

171. Стерман JI.С. Исследование теплообмена при кипении воды и этилового спирта в трубах Текст. / Л.С. Стерман, В.Г. Морозов, С.А. Ковалев // ИФЖ,- 1959 .- № 10.-С. 40-46.

172. Стырикович М.А. Методы экспериментального изучения процессов генерации пара Текст. / М.А. Стырикович, М.И. Резников.- М.: Энергия, 1977.- 278 с.

173. Стюшин Н.Г. К расчету интенсивности теплообмена при кипении в условиях вынужденного движения жидкости Текст. / Н.Г. Стюшин //Труды МИХМ, вып. 42.- М., 1972.-С. 26-31.

174. Стюшин Н.Г. К вопросу о теплообмене при кипении растворов Текст. / Н.Г. Стюшин, В.И. Астафьев //Теоретические основы химической технологии.-1975.- т. IX, № 4.-С.555-562.

175. Стюшин Н.Г. О влиянии концентрации на интенсивность теплообмена и критическую плотность теплового потока при кипении смесей и растворов Текст. / Н.Г. Стюшин, В.И. Астафьев //Теоретические основы химической технологии.- 1977.- т. XI, № 1.-С. 59-65.

176. Субботин В.И. Временные и структурные характеристики газожидкостного потока при снарядном течении Текст. /В.И.Субботин //Теплоэнергетика. 1976.-№1.-С. 67-70.

177. Суслов В.А. Результаты сопоставления опытных и расчетных данных по теплообмену при кипении сульфатного щелока Текст. / В.А. Суслов //Машины, конструирование, расчеты и оборудование ЦБП: межвуз. сб. научн. тр., вып. 6.- Л., ЛТА, 1978.-С. 96-101.

178. Суслов В.А. Исследование теплообмена при кипении слабого черного сульфатного щелока в вертикальных трубах выпарных аппаратов

179. Текст. / В.А. Суслов // Машины и аппараты ЦБП: межвуз. сб. научн. тр., вып. 8.- Л., ЛТИ ЦБП, 1980.-С. 102-105.

180. Суслов В.А. Исследование теплоотдачи при выпарке черного сульфатного щелока Текст. / В.А. Суслов, Г.В. Рижинашвили // Бумажная промышленность. 1981 .-№ 12.- С. 20-21.

181. Суслов В.А. Исследование процессов теплообмена при испарении щелоковых растворов в производстве бумаги Текст. / В.А. Суслов // Химия и технология бумаги: межвуз. сб. научн. тр., вып. 9.-Л., ЛТИ ЦБП, 1981.-С. 137-141.

182. Суслов В.А. Результаты промышленных исследований режимов работы выпарных станций ЦБП Текст. / В.А. Суслов // Машины и аппараты ЦБП : межвуз. сб. научн. тр., вып. 10.-Л., ЛТИ ЦБП, 1982.-С. 133-135.

183. Суслов В.А. Влияние режимных параметров на интенсивность теплообмена при кипении сульфатного щелока в длиннотрубных выпарных аппаратах Текст.: автореф. дис. . канд. техн. наук / Суслов Вячеслав Александрович. Л., ЛТИЦБП, 1982.-19 с.

184. Суслов В.А. Анализ теплообмена при кипении сульфатного щелока по усредненным опытным данным Текст. / В.А. Суслов //Машины и аппараты ЦБП: межвуз. сб. нучн. тр., вып. 11.-Л., 1983.-С. 113-115.

185. Суслов В.А. Исследование теплообмена при выпаривании черных сульфатных щелоков Текст. / В.А. Суслов //Теплофизика и гидродинамика процессов кипения и конденсации : материалы Всесоюзн. конф., т.1,ч.4- Рига, РПИ, 1985.-С. 117-119.

186. Суслов В.А. Позонное исследование теплоотдачи при кипении водного предгидролизата в вертикальной трубе Текст. / В.А. Суслов, В.И. Куш-нир; ВНИПИЭИлеспром. М., 1985.- Деп. № 1649-лб.

187. Суслов В.А. Проектирование выпарных установок ЦБП Текст. : учеб. пособие / В.А. Суслов, В.Л. Мовсесян, А.Ф. Мурзич, В.А. Ганичев.-Л.: ЛТА., 1987.- 83 е.: ил. (В надзаг.: ЛТИ ЦБП).

188. Суслов В.А. Ступенчатый метод расчета выпарных аппаратов Текст. / В.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС: межвуз. сб. научн. тр.- Л., ЛТИ ЦБП,1988.-С. 8-13.

189. Суслов В.А. Интенсификация работы выпарных установок глиноземного производства Текст. / В.А. Суслов, Е.М. Милин, A.M. Насопов, А.А. Яновский, Д.Р. Сандлер // Цветная металлургия.- 1988.- № 6.- С.31-35.

190. Суслов В.А. Стенд для исследования процесса накипеобразования в кипятильных трубах выпарных аппаратов ЦБП Текст. / В.А. Суслов, В.В. Семенов // Машины и аппараты ЦБП: межвуз. сб. научн. тр. Л., ЛТИ ЦБП,1989.- С. 91-92.

191. Суслов В.А. Расчет коэффициента теплопередачи в выпарных аппаратах отрасли Текст. / В.А. Суслов, Е. В. Алексеев // Бумажная промышленность.- 1989.- № 2.-С. 30-31.

192. А.С. 1646087 СССР. Выпарной аппарат Текст. / В.А. Суслов и др. (СССР) .-1991.

193. А.с. 1829154. Выпарной аппарат для кристаллизующихся растворов Текст. / В.А. Суслов и др. -1992.

194. Суслов В.А. Анализ технико-экономических показателей различных схем выпарных установок целлюлозного производства Текст. / В. А. Суслов, Г.В. Рижинашвили //Целлюлоза. Бумага. Картон. 1994.- № 11-12.- С. 21-22.

195. Суслов В.А. Исследование теплоотдачи при кипении сульфатного щелока в трубах выпарных аппаратов в зоне кольцевого течения Текст. /

196. B.А. Суслов // Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС : межвуз. сб. научн. тр.- СПб., СПбГТУ РП, 1997.1. C. 10-13.

197. Суслов В.А. Основные процессы при выпаривании щелоков целлюлозного производства Текст. : учеб. пособие для вузов /В.А. Суслов. — СПб.: СПбГТУ РП, 1998.- 92 с.

198. Суслов В.А. Оптимизация работы выпарных станций ОАО «Котласский ЦБК» Текст. / В.А. Суслов, М.А. Чайка // Проблемы экономии топливно-энергетических ресурсов на предприятиях и ТЭС: межвуз. сб. научн. тр.-СПб., СПбГТУ РП, 2002.- С. 64-69.

199. Суслов В.А. Основные процессы при выпаривании щелоков целлюлозного производства Текст.: учеб. пособие для вузов / В.А. Суслов. — СПб.: СПбГТУ РП, 2004.- 102 с.

200. Суслов В.А. Повышение надежности выпарных аппаратов при выпаривании черных сульфатных щелоков Текст. / В.А. Суслов // Целлюлоза. Бумага. Картон.- 2004.- № ю.- С. 58-61.

201. Суслов В.А. Результаты промышленных исследований режимов работы выпарных станций ЦБП и их сопоставление по эффективности выпаривания Текст. / В.А. Суслов // Изв. ВУЗов. Лесной журнал.- 2005.- № 1-2.- С. 168-172.

202. Суслов В.А. Определение коэффициентов накипеобразования на выпарных батареях целлюлозно-бумажных комбинатов Текст. / В.А. Суслов // Изв. ВУЗов. Лесной журнал.- 2005.- № 1-2.- С. 86-99.

203. Суслов В.А. Исследование теплообмена при кипении отработанных щелочных растворов в трубах выпарных аппаратов Текст. / В.А. Суслов // Теплоэнергетика.- 2005.- № 8.- С. 68-71.

204. Суслов В.А.Зависимость локальных коэффициентов теплоотдачи при кипении от параметров и структур двухфазных пенообразующих водных растворов Текст. / В.А. Суслов // Теплоэнергетика .- 2005.- № 9.- С. 74-77.

205. Суслов В.А. Исследование гидродинамики и теплоотдачи при гравитационном стекании пленки сульфатного щелока в вертикальной трубе выпарного аппарата Текст. / В.А. Суслов, Г.В. Рижинашвили // Целлюлоза. Бумага. Картон.-2005.- № 5.-С. 68-71.

206. Тананайко Ю. М. Методы расчета и исследования пленочных процессов Текст. / Ю.М. Тананайко, Е.Г. Воронцов.- Киев: Техника, 1975.- 311 е.: ил.

207. Таубман Е.И. Экспериментальное и теоретическое исследование теплообмена при испарении жидкости в стекающей пленке Текст. / Е.И. Та-убман, Ю.И. Калишевич // Тепломассообмен V : сб. науч. тр., т. IV.- Минск, 1976.-С. 437-439.

208. Гортышев Ю.Ф.Теория и техника теплофизического эксперимента Текст.: учеб. пособие для вузов / Ю.Ф. Гортышев, Ф.Н. Дресвянников, Н.С. Идиатуллин и др.; под ред. В.К. Щукина.- М.:Энергоатомиздат, 1985.-360 е.: ил.

209. Тепловой расчет котельных агрегатов (Нормативный метод) Текст.- М.:Энергия, 1973.- 358 с.

210. Накоряков В.Е.Теплообмен при пленочной конденсации и пленочном кипении в элементах оборудования АЭС Текст. / В.Е.Накоряков. М.: Энергоатомиздат,1993. - 207 е.: ил.

211. Теплопередача в двухфазном потоке Текст. / Под ред. Д.Баттерворса и Г.Хьюитта.- М.:Энергия, 1980.- 328 с.

212. Тихомиров В.К. Пены. Теория и практика их получения и разрушения Текст. / В.К. Тихомиров. М.: Химия, 1975.- 262 с.

213. Ткаченко С.И. Обобщенные методы расчета теплогидродинамиче-ских процессов и применение их для оптимизации выпарных устано-вокТекст. : автореф. дис. . докт. техн. наук / Ткаченко С.И. -М., МЭИ, 1986.

214. Ткаченко С.И.Двухфазные течения в элементах теплотехнического оборудования Текст. / С.И. Ткаченко, А.Ю. Пинчук, С.К. Кузьменко // Двухфазный поток в энергетических машинах и аппаратах : тез. докл. 8-й Всесоюзн. конф. Т.1.- Л., 1990.-С. 206-207.

215. Тобилевич Н.Ю. Исследование рабочих процессов при выпаривании сахарных растворов Текст. : автореф. дис. . докт. техн. наук / Тобилевич Н.Ю. Киев, КТИ ПП, 1969.

216. Товажнянский Л.Л. Расчет пластинчатых конденсаторов многокомпонентных парогазовых смесей Текст. / Л.Л. Товажнянский // Двухфазный поток в энергетических машинах и аппаратах. Т.2 : тез. докл. 8-й Всесоюзн. конф.-Л., 1990.-С. 29-30.

217. Толубинский В.И. Механизм теплообмена при кипении и его интенсивность Текст. / В.И. Толубинский // Тепло-и массоперенос при фазовых превращениях. Минск, 1974 -Ч. 1- С. 63-77.

218. Толубинский В.И. Влияние давления на интенсивность теплоотдачи при кипении водоспиртовых смесей Текст. / В.И. Толубинский // Теплофизика и теплотехника. Вып. 20.- Киев, 1971.-С. 3-5.

219. Тонг Л. Теплоотдача при кипении и двухфазное течение Текст. / Л. Тонг; пер. с англ. В.Я. Сидорова.- М.: Мир, 1969.- 340 с.

220. Уоллис Г. Одномерные двухфазные течения Текст. / Г. Уоллис.-М.: Мир, 1972.-440 с.

221. Федоров Г.С. Исследование распределения жидкости в трубчатых пленочных аппаратах Текст.: автореф. дис. . канд. техн. наук / Федоров Г.С.- Л.,ЛТИ, 1972.

222. Федоткин И.М. Исследование толщины пленки в аппаратах с нисходящим прямотоком Текст. / И.М. Федоткин, М.М. Чепурной // Пищевая промышленность.- 1976.- Вып. 22.- С. 36-70.

223. Федоткин И.М. Пленочные теплообменные аппараты и пути интенсификации теплообмена в них Текст. / И.М. Федоткин, В.Ф. Фирисюк.- Киев, 1969.- 92 е.: ил.

224. Филиппов Б.С. Сравнительная характеристика плотности черных щелоков, полученных при сульфатной варке хвойной и лиственной древесины Текст. / Б.С. Филиппов, С.И. Третьяков //Труды АЛТИ.- Вып. XXXII, 1972.-С. 71-76.

225. Филоненко В.И. Исследование гидродинамики и теплообмена при кипении жидкости в кольцевых каналах с двухсторонним обогревом Текст. : автореф. дис. канд. техн. наук / Филоненко В.Н. -Киев, КТИПП, 1982.

226. Френкель Б.А. Измерение малых и микрорасходов продуктов нефтехимических производств Текст.: тем.обзор / Б.А. Френкель.-М.: ЦНИИТЭ-нефтехим,1973.- 120 с.

227. Фукс С.И. Тепло-и массообмен при конденсации чистого пара и содержащего примесь воздуха при боковой подаче в трубный пучок Текст. / С.И. Фукс, Э.П. Зернова// Теплоэнергетика.- 1970.- № 3.- С. 59-63.

228. Хана О.Т. Теплоотдача при турбулентном течении в трубе жидкостей, вязкость которых зависит от температуры Текст. / О.Т. Хана, О.С. Сандал // Теплопередача.- М.,Мир, 1978.- № 2.-С. 55-61.- (Серия «С» )

229. Хартман К. Планирование эксперимента в исследовании технологических процессов Текст. / К. Хартман.- М.: Мир, 1977.- 447 с.

230. Холостых В. И. Течение пленки жидкости по вертикальной поверхности Текст. / В.И. Холостых, И.Г. Бляхер, А.А. Шехтман // ИФЖ.- 1972.-№ 3 , т.2-С. 494-498.

231. Холпанов JL П. О массообмене в пленке жидкости при волнообразовании Текст. / Л.П. Холпанов, В.Я. Шкадов, В.А. Малюсов, Н.М. Жаворонков // ТОХТ.- 1967.- № 1, т.1.- С. 73-79.

232. Холпанов Л. П. Исследование гидродинамики и массообмена в пленке жидкости с учетом входного участка Текст. / Л.П. Холпанов, В.Я. Шкадов, В.А. Малюсов, Н.М. Жаворонков // ТОХТ. -1976.- № 5, т. 10.- С. 659-669.

233. Цвелодуб О. Ю. Пространственные волновые режимы на поверхности тонкой вязкой пленки жидкости Текст. / О.Ю. Цвелодуб.-Новосибирск: ИТФ, 1991. -31 е.: ил.

234. ЗОЗ.Чепурной М.Н. Гидродинамика и теплообмен в опускных дисперсно-кольцевых потоках Текст. / М.Н. Чепурной, Е.П. Ларюшкин // Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесо-юзн.конф. Т.1.-Л, 1990.-С. 36-38.

235. Чернобыльский И.И. Гидродинамика и теплоотдача к орошающей пленке жидкости при ее гравитационном течении по вертикальной поверхности теплообмена Текст. / И.И. Чернобыльский, Е.Г. Воронцов // Тепло-и массоперенос.- М.: Энергия, 1968. Т.1.- С. 259-266.

236. Чистяков С.Ф. Теплотехнические измерения и приборы Текст. / С.Ф. Чистяков, Д.В. Радун .-М.: Энергия, 1972.- 392 с.

237. ЗОб.Чистяков В.А. Позонный теплогидравлический расчет конденсатора водяного пара Текст. / В.А. Чистяков, Г.А. Архипов, Н.И. Тищенко //Двухфазный поток в энергетических машинах и аппаратах : тез. докл. 8-й Всесоюзн. конф. Т.2.- Л., 1990.-С. 132-133.

238. Шейнкман А.Г. Особенности гравитационного волнового течения пленок жидкости по вертикальной и наклонной пластине Текст. / А.Г. Шейнкман, В.М. Пахалуев // Гидродинамика и теплообмен.- Свердловск, 1974.-С. 40-44.

239. Шкадов В. Я. К пленочной теории волновых течений пленки жидкости Текст. / В.Я. Шкадов, Л.П. Холпанов, В.А. Малюсов, Н.М. Жаворонков // ТОХТ. -1970.- Т. 4.-С. 859-867.

240. ЗЮ.Шлихтинг Г. Теория пограничного слоя Текст. / Г. Шлихтинг.-М.: Наука, 1969.- 742 с.

241. ЗИ.Элконин Н.В. Пенообразование сульфатного мыла и черного щело-каТекст. / Н.В. Элконин // Бумажная промышленность.- 1941.- № 2,- С. 2126.

242. Юнсал М. Нелинейный анализ устойчивости пленки конденсата Текст. / М. Юнсал , В. Томас // Теплопередача.- 1980.- № З.-С. 108-115.

243. Ягов В.В. Механизм теплообмена при развитом пузырьковом кипении Текст. / В.В. Ягов // Двухфазный поток в энергетических машинах и аппаратах: тез. докл. 8-й Всесоюзн. конф. Т.1.- Л., 1990.-С. 235-237.

244. Заявка № 2554071, Германия.332.3аявка № 50-35506, Япония.

245. ЗЗЗ.Заявка № 52-14713, Япония.334.Патент № 4040898, США.335.Заявка № 50-7564, Япония.336.Патент № 3097389, США.

246. Bays G.S., McAdams W.H. Heat transfer coefficients in falling film heaters, streamline flow//Industr. Engin. Chem. 1937. Vol. 29,N 11. p. 1240-1246.

247. Benjamin Т. B. Wave formation, in laminar flow down an inclined plate //J. Fluid. Mechanics. 1957. Vol. 2. P. 554-574.

248. Bodenheimer V.B. Viscosity variation on high solids black liquor. -"Southernplag and paper manufacturer", v.32, № 6, 1969, p. 54-55.

249. Brainier N., Maron D. M. Characteristics of inclined thin films, waviness and the associated mass transfer a Int. J. Heat and Mass Transfer. 1982. Vol. 25, № l,p. 99-100.

250. Brauer H. Forschungsheft, 457, 1956, 53 p.

251. Burnout and nucleation in climbing film flow / G.F.Hewitt, H.A. Kearsey, P.M.C. Lacey, D.J.Pulling Int. J. Heat Mass Trans. 1965, v. 8, p. 793.

252. Butterworth D. Simplified Methods for Condensation on a Vertical Surface with Vapor Shear. UKAEA Rept. AERE-R9683, 1981.

253. Cathro K.J. Tait W.F. Heat transfer to liquids boiling inside tubes. — "Australian journal of applied science", v. 8, № 4, 1957, p. 279-304.

254. Chen J.C. Correlation for boiling heat transfer to saturated liquids in convective flow. "Industrial Engineering Chemical Process Design and Development", № 5,1966, p. 322.

255. Chun K. R., Seban R. A. Heat transfer to evaporating liquid films //J. Heat Transfer. 1971. Vol. 93, N 4. p. 391-396.

256. Colburn A.P. Method of Correlation Forced Convection Heat-Transfer Date and a Comparison with Fluid Friction. Trans. AIChE, vol.30, pp. 187-193, 1934.

257. Dukler A.E. Fluid mechanics and heat transfer in vertical falling-film systems. // Chem. Engin. Progr. Symposium. 1960. Vol. 56, № 30, p. 1-10.

258. Esmail M. N. Wave profiles on inclined falling film // Can. J. Chem. Engin. 1980. Vol. 58, N2. P. 145-150.

259. Garwin L., Kehy B.W. Inclined falling films.// Industr. Engin. Chem. 1955. Vol. 47. № 3, p. 392-395.

260. Gerhart K. Stoff- und Warmeubergang bei der Kondesation von Dampfen aus in ringspaltstromenden Gemischen mit Luft. "VDI-Forschungsheft 539", 1970, s. 25-48.

261. Gudmundson G. Heat transfer in industrial black liquor evaporation plants. "University of technology", Geteborg, 1973, p. 161.

262. Gupta A.S., Holland F.A. Heat transfer studies in a climbing film evaporator. "Canadian journal of chemical engineering", v. 44, № 4, 1966, p. 7781.

263. Gupta A.S., Holland F.A. Heat transfer studies in a climbing film evaporator, part II: Heat transfer regions. "Canadian journal of chemical engineering", v. 44, №4, 1966, p.326-329.

264. Fiend K. Strmungsuntersuchungen bie Gregenstorm von Rieseffilmen und Cras in lotrevhten Rohren. VDI Forshungsheft, p. 481, 1960. 35

265. Folford G.D. The flow of liquid in thin films.//Advansed in Chemical Engineering Bd. 5-Y/London, Academic. 1988. Vol. 62, N 10, p. 627-631.

266. Fuita T. Theory of laminar film condensation. -New York etc.: Springer-Verl.-1991.-XYIII, p. 213

267. Fuita Т., Ueda T. Heat transfer to falling liquid films and film breakdown. Saturated liquid films with nucleate boiling // Int. J. Heat and Mass Transfer. 1978. Vol. 28, p. 109-118.

268. Fujii T. Theory of laminar film condensation. New York etc.: Springer Verl. 1991. XYIII, 213 p.

269. Haase В., Die Warmeubergang am siedenden Riesenfilm. Chem. Techn. 22, 1970, №5, s. 283-287.

270. Herbert L. S., Sterns U. J. An experimental investigation of heat transfer to water in film flow. Part 1. Non-boiling runs with and without induced swirl // Can. J. Chem. Engin. 1968. Vol. 46, p. 401-407.

271. Hoffman M. A., Potts W. W. Experimental behavior of falling liquid films at high surface tension numbers fl Industr. Engin. Chem. Fundamentals. 1979. Vol. 18, N 1, p. 27-33.

272. НоЫег Т., Lu Sin Zu. Badania urzanzen zraszajacych rury pionowe od zewnatrz. Chem. Stosow. 5. 1961. № 2, s. 153-168.

273. Kirkbride G.G. Heat transfer by condensing in vertical tube. // Ind. And Eng. Chem. 1934. V. 26. № 4, p. 432-434.

274. Javdani K., Goren S. L. Finite amplitude wavy flow of thin films // Progr. Heat and Mass Transfer. 1971-72. Vol. 6, p. 253-262.

275. McAdams W. H., Drew Т. В., Bays G. 8. Heat transfer to falling-water films // Transactions of the ASME. 1940. Vol. 62, N 10, p. 627-631.

276. Mikielewicz J., Moszynski J.R. Minimum thickness of a liquid film flowing vertically down solid surface. / Int. J. Heat and Mass Transfer. 1976. Vol. 19. № 7, p. 771-776.

277. Moles F.D. and Show J.H.G. Boiling heat transfer to subcooled liquidds under condition of forced convection/ "Trans. Institute Chemical Ingineers", v. 50, 1972, p. 76-84.

278. Nakaya C. Waves of large amplitude on a fluid film down a vertical wal I jl J. Phys. Soc. Jap. 1977. Vol. 43, N 5, p. 1821-1822.

279. Nakaya C. Solitary waves on a viscous fluid film down a vertical wall. // J. Phys. Soc. Jap. 1983. Vol. 52. N. 2, p. 359-360.

280. Nusselt W. Die Oberflachenkondensation des Wasserdamptes // Zeitschrift Vol. 1916. Bd. 60. s. 541-546.

281. Nusselt W. Delattoeustausch am Berieselungsk Uhler. // Zeitschrift VDI. 1923. Bd. 67. №9, s. 206-210.373.0shinovo Т. and Charles M.E. Vertical two-phase flow. Part 1. Flow pattern correlations. Can. J. Chem. Eng. 1947, V. 52, p. 531-540.

282. Pouter A.B., Aswald K.M. Minimum thickness of a liquid film flowing down a vertical surface-validity of Mkielewicz and Moszynski's equation. // Fot. J. Heat and Mass transfer. 1977. Vol. 20.№ 5, p. 575-576.

283. Rohsenov W.M., Webber J.H., Ling A.T. Effect of Velocity on Turbulent Film Condensation, Trans. ASME, vol.78, pp. 1637-1643, 1956.

284. Ruckenstein E. On the break up of thin liquid layers flowing along a surface. // Int. J. Heat and Mass Transfer. 1971. Vol. 14. № 1, p. 165-169.

285. Saveanu T. Chim Acad. RPR/ fil. Jasi, Rev chim. Bucharest, 21, 1960 №1; p. 149-157

286. Sexaher T. Ocr Wantleubergang am senkrechten Rohr. // Forsch. lug. Wes. 1939. Bd. 30. № 6, s. 286-296.

287. Sinek J. Raveling E. H. Heat transfer in falling long-tube vertical evaporators Chem. Engin. Progr. 1962. Vol. 58, N 12, p. 74-80.

288. Slesarenko V. Hydrodynamics and heat transfer during seawater boiling in thin films desalination plants // Desalination. V. 21., № 3, 1977.

289. Solomon T. R., Armstrong R. C, Brown. R. A. Travelling waves on vertical films: Numerical analysis using finite element method // Physics Fluid, 6(6),-1994, pp. 2202-2220.

290. Spindler B. Linear stability of liquid films with interfacial phase change // Int. J. Heat and Mass Transfer. 1982. Vol. 25, N 2, p. 161-173.

291. Staler P. J., Krantz W. B. Spatially growing three-dimensional waves on falling film flow // Intern. J. Multiphase Flow. 1977. Vol. 3, N 6, p. 609-614.

292. Streit D.E. Free Falling Evaporation Sistem. // Chemical Engineering Progress 1983. - V. 79. - № 9., pp. 41-45.

293. Strove H. Der Warmeubergang an einem verdamptenden Rieselfilm // VDI-Forschungsheft 534. Dusseldorf, 1969. 36 s.

294. Takabama H., Kato S. Longitudinal flow characteristics of vertically falling liquid films without concurrent gas flow a Int. J. Multiphase Flow. 1980. Vol. 6, N 3, p. 203-215.

295. Takeda H., Hayakawa Т., Fujita S. Vapor hold up boiling heat transfer coefficient of natural circulation vertical tube evaporator. "Kagaku kagaku", v. 36, №7, 1972, p. 760-768.

296. Thomas D.G., Young G. Ind. and Eng. Chem. Process Des. and Develop. 9, 1970, №2, p. 317-323.

297. Tougou H. Long waves on a film flow of a viscous fluid down the surface of a vertical cylinder a J. Phys. Soc. Jap. 1977. Vol. 43, N 1, p. 318-325.

298. Ueda Т., Inane M., Nagatome S. Critical heat flux and droplet entrapment rate in boilling of falling liquid films // Int. J. Heat and Mass Transfer. 1981. Vol. 24, N7, p. 1257-1266.

299. Unterberg W., Edwards O.K. Evaporation from falling saline water films in laminar transitional flow.// AIChE Journal. 1965. Vol. 11, № 6, p. 1073-1080.

300. Wang С. Y. Liquid film flowing slowly down a wavy incline // AIChE Journal. 1981. Vol. 27, N2, p. 207-212.

301. Wilke W. Warmelibergang an Rieselfilme 11 VDI-Forschungsheft 490. 1962. 36 s.

302. Zuber N., Staub F. W. Stability of dry patches forming in liquid films flowing over heated surfaces // Int. J. Heat and Mass Transfer. 1966. Vol. 9, p. 897-905.