Исследование кипения водных растворов при повышенных давлениях и усовершенствование методики расчета испарителей кипящего типа при закритической минерализации

На правах рукописи

Лавриков александр владимирович

ИССЛЕДОВАНИЕ кипения водных растворов При повышенных давлениях и усовершенствование методики расчета испарителей кипящего типа при Закритической Минерализации

Специальность:01.04.14– " Теплофизика и теоретическая теплотехника "

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук

МОСКВА – 2008

Работа выполнена в Московском энергетическом институте (техническом университете) на кафедре инженерной теплофизики.

Научный руководитель: доктор технических наук,

профессор Кузма-Кичта Юрий Альфредович

Официальные оппоненты: доктор технических наук, профессор,

Ковалев Сергей Алексеевич

доктор технических наук, Бухаров Александр Васильевич

Ведущая организация:

Московский государственный университет инженерной экологии

Защита состоится 23 мая 2008 года в 10 часов 00 мин. на заседании диссертационного совета Д 212.157.04 при Московском энергетическом институте (техническом университете) по адресу: 111250, Москва, Красноказарменная ул., д.17, корп. Т, кафедра инженерной теплофизики, комн. Т-206.

С диссертацией можно ознакомиться в научно-технической библиотеке Московского энергетического института (технического университета).

Автореферат разослан апреля 2008

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

Ученый секретарь

диссертационного совета Д 212.157.04

к.ф.–м.н., доцент __________ Мика В.И.

Актуальность темы. В настоящее время термический метод водоподготовки на ТЭС реализуется с помощью испарителей кипящего типа с естественной циркуляцией. Исследование работы испарителей при глубоком концентрировании питательной воды выявило ряд особенностей их гидродинамического и теплового режима по сравнению со случаем маломинерализованной среды, не учитываемых моделями^[1],^[2]

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

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

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

Целью работы является:

Экспериментальное исследование внутренних характеристик кипения водных растворов в большом объеме при атмосферном и повышенном давлении применительно к условиям в испарителях кипящего типа, а также перепада давления при кипении водных растворов в трубе при низких массовых скоростях.

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

Научная новизна.

- Впервые обнаружено увеличение частоты отрыва паровых пузырей при кипении в большом объеме водного раствора Na2SO4 по сравнению с водой при атмосферном давлении;

- впервые получены экспериментальные данные по скорости всплытия паровых пузырей и фактору взаимодействия для водного раствора Na2SO4 в диапазонах давлений от 0.15 до 1 МПа и концентраций от 0 до 40 г/л;

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

- впервые определены составляющие полного перепада давлений в вертикальной трубе при кипении водного раствора Na2SO4. Обнаружено, что основной вклад в увеличение полного перепада давления по сравнению с водой вносят потери на трение;

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

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

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

Апробация работы. Основные положения и результаты работы были доложены: на одиннадцатой международной научно-технической конференции студентов и аспирантов “Радиоэлектроника, Электротехника и Энергетика”, МЭИ 2005 г.; XV международной школе-семинаре молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева, Калуга 2005 г; восьмой международной научно-технической конференции “Оптические Методы Исследования Потоков (ОИМП-2005)”, МЭИ 2005 г.; четвертой российской национальной конференции по теплообмену (РНКТ-4), Москва, МЭИ, октябрь 2006г.; XVI международной школе-семинаре молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева, Санкт Петербург 2007г.; девятой международной научно-технической конференции “Оптические Методы Исследования Потоков (ОИМП-2007)”, МЭИ 2007 г.; четырнадцатой международной научно-технической конференции студентов и аспирантов “Радиоэлектроника, Электротехника и Энергетика”, МЭИ 2008 г.; на заседании кафедры ИТФ МЭИ 5 марта 2008 г.

Публикации. Основное содержание работы отражено в 8 публикациях.

Структура и объем работы. Диссертация состоит из введения, четырех глав, выводов, списка использованной литературы и приложе­ний. Работа содержит 151 страниц основного машинописного текста, 89 рисунков, 5 таблиц, 17 страниц приложений, библиография содержит 59 наименований.

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

Во введении обоснована актуальность диссертационной работы.

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

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

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

- отсутствуют данные по скорости всплытия и эквивалентным диаметрам паровых пузырей, фактору взаимодействия для водных растворов при повышенных давлениях;

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

Во второй главе приводятся описания экспериментальных стендов, методик исследования.

Для исследования пузырькового кипения в большом объёме при атмосферном давлении усовершенствована установка, применявшаяся ранее для изучения механизма теплоотдачи при плёночном кипении. Установка (рис. 1) состоит из рабочего контура, оптической системы и измерительно-вычислительного комплекса. Оптическая система, состоит из гелий-неонового лазера (2), пучок излучения которого проходит через коллиматор (1), диафрагму и рабочую камеру (7). После рабочей камеры пучок излучения попадает на фотоприёмник (5), сигнал которого регистрируется цифровым быстродействующим вольтметром В7-43 (4) и измерительно-вычислительным комплексом (3). Рабочая камера представляет собой цилиндр, по торцам которого закреплены кварцевые стёкла. В качестве поверхности нагрева использовалась трубка (6) из стали 1Х18Н9Т наружным диаметром 3 мм, толщиной стенки 0,5 мм, длиной 80 мм. Нагрев рабочего участка осуществляется переменным электрическим током. Сигнал фотодатчика обрабатывается с помощью специально написанной программы, которая позволяет получать распределения частот отрыва, скоростей всплытия и отрывного диаметра парового пузырька для большого массива данных. Достоверность полученных данных проверялась в опытах на воде.

Исследование диаметров и скорости всплытия паровых пузырей для воды и водных растворов в диапазоне давлений 0.1-1 МПа проведено на установке, схема которой представлена на рис. 2. Рабочая камера (2) представляет собой толстостенный цилиндр, по торцам которого прикреплены кварцевые стекла. Кипение происходит на внутренней поверхности рабочей камеры, которая нагревается охранным нагревателем (1). Пар, образующийся в рабочей камере, поступает в холодильник (10). Давление в рабочей камере измеряется манометром (7). Сигналы датчиков поступают на измерительную карту(8) и компьютер (9).

Для измерения диаметров и скорости всплытия паровых пузырей применялась скоростная видеосъемка с помощью цифровой камеры (6). Съемка осуществляется в проходящем свете. Фокусное расстояние объектива (5) составляло 150 мм. Скорость съемки составляла 1000 кадров в секунду. Гильза (3) служила масштабом видеосъемки.

Распределения скорости всплытия паровых пузырей для воды при давлении 1 МПа, полученные в результате автоматической (а) и ручной обработки (б), представлены на рис. 3. В случае автоматической обработки распределение построено по данным для 300 паровых пузырей и при ручной обработке – 200 паровых пузырей. Наиболее вероятные значения, полученные при ручной и машинной обработке, отличаются не более 30%. Интересно отметить, что длительность ручной обработки составляет в данном случае ~10 часов и автоматической обработки ~2 часа..

Исследование гидродинамики при кипении водных растворов в трубе при пониженных массовых скоростях проведено на автоматизированном экспериментальном стенде (рис. 4). Установка обеспечивает проведение измерений в диапазонах давлений 0.1 – 7 МПа, массовых скоростей 10 – 650 кг/(м2с), тепловых нагрузок до 1.5 МВт/м2. Рабочий участок представляет собой трубу (1) с технически гладкой поверхностью из стали Х18Н10Т с внутренним диаметром 6.89 мм и толщиной стенки 0.55 мм. Нагреваемая длина трубы составила 1775 мм. Для разгрузки от давления канал помещен в толстостенный кожух (2). На входе в рабочий участок установлен сильфон (3). После холодильника (5) рабочая среда подается насосом (6) через ресивер (7) в рабочий участок. Расход жидкости измеряется датчиком Метран (4), перепад давления на рабочем участке - датчиком давления “Сапфир 22Д” (5). Данные измерений обрабатываются с помощью измерительно-вычислительного комплекса.

В третьей главе приводятся результаты экспериментального исследования. Распределения отрывных диаметров паровых пузырей при различных тепловых нагрузках для воды и водного раствора Na2SO4 с концентрациями 8 г/л и 30г/л представлены на рис.5 и рис.6. В кружках на рисунках показано количество обработанных паровых пузырей. По вертикали отложены доли паровых пузырьков данного диаметра в общем количестве пузырей на распределении. Из анализа распределений установлено влияние тепловой нагрузки и концентрации раствора на дисперсию. С ростом тепловой нагрузки дисперсия отрывных диаметров паровых пузырей уменьшается (рис. 5). Распределение отрывных диаметров паровых пузырей для раствора Na2SO4 с концентрацией 8 г/л имеет меньшую дисперсию, чем для воды и водного раствора Na2SO4 с концентрацией 30 г/л при примерно одинаковой тепловой нагрузке.

Полученные распределения позволяют надежно определить наиболее вероятные значения измеряемых характеристик процесса кипения. Сравнение опытных данных по отрывным диаметрам паровых пузырей при атмосферном давлении для воды с расчетными (рис.7) показало, что формула Ягова В.В. описывает результаты с отклонением до 30%. Каждая точка на графике является наиболее вероятной величиной, полученной при обработке распределений, построенных по массивам из 600-5000 значений.

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

Видеокадры всплытия паровых пузырей для воды и водного раствора Na2SO4 с концентрацией 20 г/л и 40 г/л показаны на рис. 9. На одном из кадров (рис. 9в) отмечена пена вблизи зеркала испарения. В работе впервые проведены измерения скорости всплытия и диаметра всплывающих паровых пузырей для водного раствора Na2SO4 с концентрацией до 40 г/л при давлениях 0.1-1 МПа.

Рис. 5. Распределения отрывных диаметров паровых пузырей для водного раствора Na2SO4, C= 8 г/л, при различных тепловых нагрузках, р=0.1 МПа. а) q=50 кВт/ м2; б) q=96 кВт/ м2 ; в) q=150 кВт/м2

Рис. 6. Распределения отрывных диаметров паровых пузырей для воды и водного раствора Na2SO4, 8 г/л и 30 г/л соответственно. р=0.1 МПа. а) q=98 кВт/м2; б) q=96 кВт/ м2; в) q=100 кВт/м2

При кипении водных растворов в случае повышенных давлений выделены три группы паровых пузырьков:

-с диаметром d>1 мм, их скорость всплытия меньше, чем для воды и зависит от концентрации раствора;

-с диаметром 0.5 мм<d<1 мм, для этой группы паровых пузырей влияние концентрации водного раствора на скорость всплытия не проявляется;

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

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

Рис.9. Кадры видеосъемки. а) Вода, p=0.8 МПа; б) Раствор Na2SO4, 40 г/л, р=0.5 МПа; в) Раствор Na2SO4, 20 г/л, р= 1 МПа

Распределения диаметров всплывающих паровых пузырей для воды и водного раствора Na2SO4 c концентрацией 20 г/л при различных давлениях представлены на рис 10а, б, в, г. Обращает на себя внимание то, что максимум функции распределения диаметров всплывающих паровых пузырей для водного раствора смещен в сторону меньших значений. Данные видеосъемки позволяют предположить, что это вызвано пенообразованием при кипении водного раствора сульфата натрия.

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

Для изучения взаимодействия паровых пузырей при всплытии в большом объеме были построены трехмерные распределения. Пример такого экзотического распределения для водного раствора Na2SO4 с концентрацией 8 г/л при р=0.1 МПа представлен на рис. 12. Полученное распределение иллюстрирует взаимодействие паровых пузырей при их всплытии.

Для определения значений фактора взаимодействия для водных растворов на основе полученных распределений (рис. 10) были найдены относительные доли крупных паровых пузырей и сравнены с аналогичными величинами для воды. Результаты такой обработки представлены на рис. 13, где показаны экспериментальные данные по доле крупных паровых пузырей (d>0.4 мм) по отношению к аналогичной величине при атмосферном давлении для воды при кипении в большом объеме. Поскольку в исследованных условиях возможно слияние паровых пузырей, то эти данные можно применить для оценки фактора взаимодействия при кипении водных растворов в трубе.

Рис. 10. Распределения всплывающих паровых пузырей по диаметрам.

Рис.13. Отношение доли крупных паровых пузырей для воды и водных растворов Na2SO4 по сравнению с аналогичной величиной для воды при атмосферном давлении. Экспериментальные данные: 1) вода; 2) раствор Na2SO4, 20 г/л; 3) раствор Na2SO4, 40 г/л; Расчетные данные: 4) формула (2) при b=0; 5) отношение факторов взаимодействия для воды при повышенном и атмосферном давлениях [1] формула (3)

Формула, описывающая экспериментальные данные по фактору взаимодействия, представлена в виде

, (1)

где

, (2)

, (3)

b– коэффициент, выбранный на основе сопоставления с имеющимися опытными данными по фактору взаимодействия для воды; С - концентрация, % ; p – давление среды, МПа.

В случае кипения в большом объеме b =0, при кипении в трубе b =0.026. Формула (1) получена для диапазона концентраций Na2SO4 от 0 до 40 г/л и давлениях от 0.1 до 1 МПа.

На рис. 14 представлены данные по полному перепаду давления на рабочем участке при кипении водных растворов в трубе и отдельным составляющим потерь давления (рис. 15).

Данные получены при давлении p=0.7 МПа и скоростях циркуляции: w=0.19 м/с, w=0.18 м/с, w=0.16 м/с. Вклады нивелирной составляющей и потерь на ускорение потока не превышают 10% от полного перепада в случае воды и водного раствора. Местные гидравлические сопротивления, напротив, составляют значительную часть от общего перепада давления в трубе. Потери на трение для водного раствора больше, чем для воды, и сильно зависят от паросодержания. Из полученных данных можно заключить, что основной причиной повышения гидравлического сопротивления при кипении водных растворов в трубе являются потери на трение.

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

Фактор взаимодействия для опускной щели в случае растворов закритической минерализации:

.

(4)

Для труб греющей секции для фактора взаимодействия предложено использовать формулу (1).

Как показали измерения Карцева А.С., скорость циркуляции при закритической минерализации в испарителе И-600 составляет w0~0.1-0.2 м/c. Условие захвата пара в опускную щель имеет вид:

w0>wпуз.

(5)

и в рассмотренных условиях оно выполняется, так как wпуз~ 0.03-0.1 м/c

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

.

(6)

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

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

, ,

(7)

где- площадь сечения труб греющей секции;

- площадь сечения, через которое идет слив из греющей секции в опускную щель;

- массовое расходное паросодержание на выходе из труб греющей секции.

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

Коэффициенты теплопередачи, рассчитанные с учетом особенностей гидродинамики испарителя и границ ухудшения теплообмена для водных растворов, сопоставлены на рис. 16 с известными опытными данными для испарителя И-600. Расхождение между экспериментальными и рассчитанными коэффициентами теплопередачи не превышает 20% для 90% точек.

Выводы

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

2. Проведено исследование кипения воды и водного раствора Na2SO4 с концентрациями 8 г/л и 30 г/л на горизонтальной трубке в большом объеме при атмосферном давлении. Получены данные по теплоотдаче и внутренним характеристикам кипения. Построены распределения скоростей всплытия паровых пузырей. Впервые получены распределения частот отрыва и отрывных диаметров паровых пузырей для водного раствора Na2SO4. Получены трехмерные распределения отрывных диаметров паровых пузырей, вид которых существенно меняется для последовательных интервалов времени за период наблюдения. Обнаружено, что для каждого диаметра всплывающего парового пузыря существует свое распределение скоростей всплытия. Установлено, что при кипении водных растворов происходит увеличение частоты отрыва паровых пузырей по сравнению с водой.

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

4. Разработаны установка и методика исследования кипения в большом объеме при повышенных давлениях с помощью высокоскоростной цифровой видеокамеры. Впервые проведены измерения скорости всплытия и диаметра всплывающих паровых пузырей для водного раствора Na2SO4 с концентрацией от 0 до 40 г/л при давлениях 0.1-1 МПа. При кипении водных растворов в области повышенных давлений выделены три группы паровых пузырей: с диаметром d>1 мм, их скорость всплытия меньше, чем для воды и зависит от концентрации раствора; с диаметром 0.5 мм<d<1 мм, для этой группы паровых пузырей влияние концентрации водного раствора на скорость всплытия не проявляется; пена, состоящая из паровых пузырей диаметром d<0.4 мм, движется с жидкостью или увлекается большими паровыми пузырями.

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

6. Получены данные по перепаду давления при кипении в трубе раствора Na2SO4 при концентрациях до 10 г/л для w=0,160,19 м/с в диапазоне х=0,10,87. Обнаружено, что с ростом минерализации среды гидравлическое сопротивление при кипении водных растворов в трубе возрастает до 30%. Для исследованных условий определены вклады отдельных составляющих полного перепада давления при кипении в трубе воды и водного раствора.

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

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

1. Седлов А.С., Кузма-Кичта Ю.А., Коньков Е.О., Лавриков А.В., Хромов А.Н. Влияние концентраций питательной воды на режимы работы водоподготовительных испарительных установок // Журнал Энергосбережение и Водоподготовка, декабрь, 2006. N 6(44), c.73-74.
2. Лавриков А.В., Кузма-Кичта Ю.А. Исследование характеристик кипения водных растворов с использованием лазерной диагностики // Одиннадцатая международная научно-техническая конференция студентов и аспирантов “Радиоэлектроника, Электротехника и Энергетика”. Тезисы докладов. М.: Издательство МЭИ, 2005, Т.3. с.54-55.
3. Кузма-Кичта Ю.А., Седлов А.С, Карцев А.С., Коньков Е.О., Лавриков А.В. Исследование кипения водных растворов солей с помощью лазерной диагностики // Труды XV международной Школы-семинара молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева. М.: Издательство МЭИ, 2005, Т.1. c.220-224.
4. Кузма-Кичта Ю.А., Седлов А.С, Карцев А.С., Коньков Е.О., Лавриков А.В. Лазерная диагностика кипения водных растворов // Восьмая международная научно-техническая конференция “Оптические Методы Исследования Потоков (ОИМП-2005)”. Труды конференции. М.: 2005, с.432-435.
5. Кузма-Кичта Ю.А., Седлов А.С.,. Коньков Е.О, Лавриков А.В. Исследование гидравлического сопротивления при кипении водных растворов в канале // Четвертая Российская национальная конференция по теплообмену (РНКТ-4). Труды конференции. М.: Издательский дом МЭИ, 2006. Т.5. с.159-162
6. Лавриков А.В., Кузма-Кичта Ю.А., Стенина. Н.А. Лазерная диагностика характеристик кипения воды и водных растворов // Труды XVI международной Школы-семинара молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева. 21–25 мая 2007, г. Cанкт-Петербург.Т2, с 473-426
7. Кузма-Кичта Ю.А., Лавриков А.В., Шустов М.В. Исследование кипения водных растворов с помощью лазерной диагностики // Девятая международная научно-техническая конференция “Оптические Методы Исследования Потоков (ОИМП-2007)”. Труды конференции. М.: 2007, с.432-435.
8. Лавриков А.В., Шустов М.В. Афонин С.Ю., Кузма-Кичта Ю.А. Исследование кипения водных растворов с помощью лазерной диагностики // Четырнадцатая международная научно-техническая конференция студентов и аспирантов “Радиоэлектроника, Электротехника и Энергетика”. Тезисы докладов. М.: Издательство МЭИ, 2008, Т.3. с.16-17.

Подписано в печать Л-

Печ.л. Тираж Заказ.

Типография МЭИ (ТУ), Красноказарменная, 13.

---

^[1] Кутепов А.М., Стерман Л.С., Стюшин Н.Г. Гидродинамика и теплообмен при парообразовании. – М.: Высшая школа, 1986. – 392 с.

^[2] Седлов А.С., Кузма-Кичта Ю.А. Гидродинамика и теплообмен при кипении водных растворов. Москва. Изд МЭИ, 2007. – 164 c.