Элементы гидродинамики многофазной жидкости применительно к исследованию движения двухфазного теплоносителя в ядерных установках тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Бабицкий, Арнольд Фелицианович АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
1989 ГОД ЗАЩИТЫ
   
01.02.05 КОД ВАК РФ
Автореферат по механике на тему «Элементы гидродинамики многофазной жидкости применительно к исследованию движения двухфазного теплоносителя в ядерных установках»
 
Автореферат диссертации на тему "Элементы гидродинамики многофазной жидкости применительно к исследованию движения двухфазного теплоносителя в ядерных установках"

Московский ордена Ленина, ордена Трудового Красного Зншзгни, ордена Октябрьской революции государственный университет им. М.В.Ломоносова

&

/

с

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

Бабицкий Арнольд Фелицианович

УДК 532.525.6

ЭЛ$!ИШ ГИДР0ДИНА1МШ МНОГОФАЗНОЙ ЩЩСОСШ ГШШШШТЕЛЬНО К ИССЛЕДОВАНИЮ ДВИЖЕНИЯ ДВУХФАЗНОГО &ЮНОСИТЕЛЯ В ЯДЕРНЫХ УСТАНОВКАХ

(специальность 01.02.05-механика зхидкостей.гаэа и плаашзУ

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

Москва - 1989

Работа выполнена на кафедре вычислительной ыатематики Киевского ордена Ленина государственного университета им. Т.Г.Шевченко и в Киевском институте автоматики им.ХХУ съезда КПСС

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

О.В.ВОИНОВ;

.доктор технических наук, профессор Г.А.САЛТАНОВ;

доктор физико-математических наук, профессор Г.А.ТИРСКИЙ

Ведущая организация - Институт проблем механики АН СССР

Защита состоится " " ■_19 г. в _час.,

в вуд. __ но заседании специализированного совета Д.053.05.02

при Московской государственной университете им.!,!.В.Ломоносова.

I

Адрес: II9899,г.Москва, Ленинские горы, Главное здание ¡.¡ГУ.

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

Автореферат разослан " " -_ 19 г.

В.П.Карликов

Общая характеристика работы Актуальность проблем«. Развитие современной техники в значительной степени завися?'от развития фундаментальных науя. С пояз-ленком все новых и новых технических проблем аозннкаот и научанэ проблеет, одни из которых могут быть решены э рамках уке алогяв-0ЮССЯ научных представлений, другие же могут быть разрешены путем развитая фундаментальных наук. Так, например, тучно-тсхня-ческиа прсблскы, свазениые тем или икым образом с исследованием двняения спловннх сред, требуют дальнейаегб развития а такой дросноЯ науки, аак гидромеханика. И, если ранео тзхнпчеггоге задачи э этой "области ренадись на сскавэ шдронохмглкя пли газ одел д-»пкм (0днс£0!!т1сц?,чтпэй иеезкмаеыой или сгак-мизной гцдкоета), та • сэйпас необходаг:о решать с сдачи двкленяя кногокемюнентных ерзд, которцз доязнкм образом но гут быть разреяекц увэ а рэдках нового раздела, мзхенкхн - гидромеханики кногофазис1 гздкогтн. В настоящее время з гвязи с интенсификацией рабочих процоссоз п удеди-чгнмеи ^хорЪстл двиденил различных малин м механизмов зрут научно-технических проблей, связанных с ксслэдовашен течений етогофазных яидяоетей5*, быстро распиряется. К наиболее Ейдакя областям технического прилоязния гидромеханики кногофазкоИ жидкости молю отнести:

- дзияение сходных рабочих тел в различит иалмнах и квха-нязнах, включая атомную 'энергетику, реактпвнуя я рзяетнуэ тахн?~ну5

* Под гдаогофазной жидвостьа коми понимается Еедхоети лед газа п любого конечного чкег.л компонентов, зоторгз г-гогуг няходл??— гл в тзбом фазовом состоял;:,« : '¿зердом, гадком, г гиюсС; рознен.

- А -

- внйшео обтекание тел, шшочая обладеменмо сгшоязтов и об-^оггаша газом и частицами поверхности коо:шчсскюс аппаратов ;

- трубопроводный транспорт (гидротранспорт), движение сыпу-Щ1Х сред в погрузсчно-разгрузочных ьгрггатех ;

- внешние р внутренние кавитацконнио течения ;

- движение смесей жидкостей, газов н частиц в различных тех-

о

»¡алогических системах, включая барботая ;

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

- рад вопросов метеорология, особенно в области направленного ьоздоРстш на процессы выпадания осадков 1

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

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

тоэряз а бол1£кнст80 работ носит отпечаток узкотехничгской нглравяеняоети. а теоряэ движения вз-югоко&поиентных сред чщг гссрэ ^ассватрявеэг как некотороа спецггфичгсгоз допохиеико тех йхз кках рм~елсв обычной гвдродиндозш! шш газодкн&цлкя. По-еуцо2?зу а обяас«^ гвдродкнаикка многофазной кцдкоста заканчи--Сй5*ед пзрзад нааойдскЕЯ отдельных фактов и начинается период их егсяс^атаэацгл ¡а обобщения, о чеа свидетельствует гмхсд кз" печа-

га ряда нокорряфий как у пае в страна, чал и за рубеаси. Очзрздной этап развития гэдроданашнад шорофазной падкости - полученнэ еиа-датачасхях, полузшшрнчесинх н затем численных рсзенйй рада зласео! задач, енологячно тоиу, как это Сто оделено в классической гэдрс-пэхеипхз. Нал . таз вычислительной техники, существенно раетаряг^гЯ зозмотности реализации решений сгонных фнзико-иатематичзсют задач, осгестзенно, но нсклзочазт ¡юэташого развитая науки, на г/о-пд? едалать болез легким ргнгинв постоянно возшшащах научных пробяэа. Данная работа н посвящена ргзенгаз некоторых нэ актуальных проблем а области теоретической н прикладной Гйдродкнамшш иногсфагноЯ падкости.

Цель работы и предиэт исследований.

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

- Еияаленна общих закономерностей течения ¡многофазной падкости, характерных для плрокого круга задач. В деннс . случае рас-сиатртс^хед течения механических сносей {

- .цо.^енаэ первых интегралов уравнений движения и построение на Щ основе точных и приближенных аналитических резенай э&-дач ро ^еценир двухфазной падкости;

- разработка теоратичйских основ гидравлика гяогофазной. падкости {

- получение численного ргаения задач по течения двухфазно* кщкссул с учотом тепло!'лссообиена я динашпеского взаянодейст-сет комто.ченгоу смеси ;

Л тагов проктичасгсоэ пр-с:еме:пз полученных результатов тсо-рскгос:;:« иселодозаниЯ для рзгенгл ряда прикладных задач э со' -лаяти атогеной знергэтщгл д! приборостроения.

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

о

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

- получены первые интегралы уравнений движения типа интег -р&лов Бернулли, Громеко, Коши-Лагранжа, Эйлера;

- получены безразмерные параметры и критерии подобия (числа Эйлера, Рейнольдса, Струхаля, Фруда, Вебера и т.п.) ;

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

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

- на основе проинтегрированного уравнения движения получены полуэмпирические зависимости для определения потерь на трение в трубах, а также при резком и постепенном расширении двухфазного потока (формула Борда-Карно) ;

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

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

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

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

Практическая ценность исследований.

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

Результаты комплексных исследований, вклвчаяцих теоретические решения задач, иняенерно-конструкторские разработки и расчеты", с цайьк) создания приборов и тарировочных устройств для измерения объеыйых концентраций компонентов смеси по методу оа-

Н V -.у .

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

веса смеси й по методу кондуктометрии в 1974 г. были переданы на хоздоговорной основе и внедрены во ВНИИ аналитического при -боростроения Министерства приборостроения и средств автоматизации СССР. А образец прибора по измерении методом кондуктометрии объемных концентраций компонентов двухфазного потока в 1970-1971 годах "ак экспонат Института гидромеханики АН УСЙР (авторы А.Ф.

- в -

А.$. Бабицкий, А.М. Кравченко) экспонировался в павильоне "Наука на ВДНХ УССР.

Объем и структура работы.

Диссертация состоит из введения, Семи глав и заключения, содержание которое изложено на 240 страница* машинописного текста, иллюстрировано 61-м рисунком и одной таблицей. Библиография - 233 наименований.

Краткое содержание.

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

Что'касается истории вопроса, то одной из первых работ в области гидродинамики двухфазной-жидкости является статья В М.Иак-хавеева по теория турбулентного взвешивания насосов (1931 г.). В начале сороковых годов (1941) С.А.Христиановичем, а в конце сороковых годов (1949) Ы.Д.Мнллнонщиковым делаются первые постановки н решения а«дач по движение газированной жидкости в пористой среде. К этому же периоду относится работа Л.Д.Ландау "К гидродинамике гелия II", в которой на основе уравнений гидродинамики двухфазной (двухкомпонентной, состоящей из гелия I и гелия И) жидкости рассмотрены вопросы проницаемого движения гелия П.' Гидродинамические уравнения дгухфазного потока на основе осреднения пара-нэтров компонентов смеси по обгыу были записаны С.Г.Телетовш такзе в сороковых годах. Но теоретическое изучение течений сые-еэй в то врэыя носило одиночный и эпизодический характер, В пятидесяти* годах в связи с развитием гидротранспорта и строительствам гидросооружений появляется целый рпд работ по пщрод1шаыике сносей несЕимаеиой кидкоста и твердых частиц. 'Гак е

1951-1952 годах Н.А.Сквзкяшм опубликованы "Дифференциальны® уравнения двияення пулыш". Затем в 1953 году Г.И.Баренблатт получает на основе осреднения по времени уравнения для двм -пения взвешенных часткц в турбулентном потоке. А в 1956 году Х.А. Рахматулкн получил гидродинамические уравнения для потока сиеся идеальной сжимаемой жидкости н твердых частиц. Су -щественно новой чертой этой работы является то, что в ней при выводе уравнений было использовано понятно о материальном континуума.

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

Л.И.Седова и среди них такие, как "Введение в механику сплоя-

о

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

а

еплопных вред". К наиболее актуальным вопроса« в плане развития теории по исследована» течений многофазной аидкостп следует отнести:

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

- разре-ботяа кинематики ¡д > счэ г онл он гн тинх спап«:н1гх зрзц;

- ттецатнческое «писание кгхеннчвекокчь гепливого, ха-«гч» •.**го, а в ряде случаев я биологического всая&.эдэ1*с9эчя ;:скпснгнтов движущейся многофазной тадковг:, особенно

наличии фазовых превращений;

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

- создание' теоретических основ гидравлики многофазной жидкости , д

Наряду с развитием теории, естественно, развивались и экспериментальные исследования в области гидродинамики многофазной видкости. К наиболее существенным работам в зтой области, имеющим не только прикладное, но и большое научное значение', следует отнести экспериментальные работы Н.И.Семенова, С.И.Костерина ц А.И.Гуяова, В.Ф.Медведева по установлению критических явлений л определению скорости звука в потоке паро-газо-водяных смесей. В настоящее время имеется довольно большое число работ экспериментального плана, значительная часть которых обобщена в монографиях С.С.Кутателадзе, Н.А.Стыриковича, М.Е.Дейча, В.С.Данилина, Г.В.Цкклаури, Г.А.Салтанова, С.Соу« Однако следует отметить, что экспериментальные исследования еще не достигли такой полноты и совершенства, хан в гидродинамике однофазной жидкости, чем и объясняется полное отсутствие справочной литературы по гидродинамике к гидравлике многофазной жидкости. По-видимому, сдерживающим фактором в дальнейшем прогрессе экспериментальных работ является отсутствие аналитических и полуэмпирических решений ряда прикладных задач. В плане экспериментальных исследований, на наш взгляд, наиболее важными и актуальнчми вопросами являются:

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

- выявление г раскрытие закономерностей, определяющих как

- и - _

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

- создание экспериментальных основ гидравлики многофазной жидкости ;

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

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

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

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

- Г, 2 -

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

В данной работе для решения ряда задач применены как методы классической гидродинамики, обобщенны*-? на случай многокомпонентной среды, так и численные методы с использованием ЭЕ&. Такой подход позволяет сочетать универсальность и относительную простоту аналитических решений, дающих главным образом, качественный анализ, с конкретность» детализированных численных решений задач прикладного характера, дающих подробный и в ряде случаев более точный количественный результат. Рассмотренные в работе вопросы »сорви и полученные» решения как теоретических, так ы прикладных задач по течения паро-газо-жидкостных смесей безусловно в той или иной степени могут быть применены во шогих научно-технических областях, но общая направленность в прикладном еначении подчинена исследованию декесния рабочих тел в элементах иаакн и механизмов, .в основном, атомной энергетики.

'Ь первой главе существенное вниианне уделено формулировали» фаэмко-математнческой модели смеси я принятию рабочих гипотез, иа основе которых в дальнейшем строится математический аппарат ддш рвения различна а сдач. И сходя из понятия о многофазном кон-тпнууш и гмпитазы о раздельном действии сил вязкости в области г^мшодрРствия и ълутря самих компонентов принимается следующая

кодель многофазной аядкости:

- кйзэдый Компонент в отДЗАЬМости и смесь в цэлом представляют собой материальные коИ№нууш, т.о. предполагается, что каждый компонент з&йолйяЪ* Пространство занимаемого смесью объема непрерывно, йо сОхрМййт свон индивидуальные характеристики р. , С •

-V » Р » Т «С (плотность, объемную концентрацию, скоросТь, давление, температуру, удельную теплоемкость, j -го компонента и др.} ;

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

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

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

I

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

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

д:>;н£шш! системы к многокомпонентной сплошной среде.Ваедоннсэ по-

—♦

нятие о базовой скорости , с которой длилетсл наблюдатель

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

На основе принятой модгзли сиеси и обобщенного закона Ньэтенл дчя С2язи з отдельном г;схпакзнте стсл тзпэзра. напрлзеиий с т'.и-ззром скорости дьфор.шз';й зг-гшсаим з общзи е>*д» ¡гг.г.гхгмя иор?.';-швйсл'», двийенкя (з йераа уравнений Навье-С*сйеа) знзргнп (пря-тока ;епла/ для ] -х компонентов многофазной остл пра налички Лвдозих превращений

- н -

• * > (I)

(

(¿-о-

Я 17м- + а' (-° =

* йГ «у-4 •"< /

' "ЭР. < Ъ — <2)

* V + у/,

_ ,(3)

' 3'

' ^ -'удвлышй пэток, скорость отдедшзш

а шугранняя энергия уаасы фазовых про вращений.

Л - цензор напрякений, еоэнийшцшс при взаимодействии -го ссклгснснге е другкщ компонентам смеси.

В рэоузгьтаи? снзишеа гндродинаашеахих уравнений шогофаз-цгП авдкоета (I}, (2), (3) показано, что безргцшзрныо параметра и кратаран подобгая агграхиттся ч®раз слалогичнио параметры однофазной сздкостй к массовые концентрации кошонентов смеси

Чаааа С^рухсял, Зйлере, Пето :

Числа Фруда, Рейиольдса, Вебера для многофазной яядковтп

= ^ А. .

ЯЧ Гг. ' Яе Яе. '

Тепло а цассоперенос мэзду компонентам:! егдаси характеризуется

линейным разбором { I, ), чяслон Стэнтона( Зп )> параметром

1 ) определявши» удзлышй потоа фазовых прврк^енпй, а отасзрняеэ

поверхности раздола фаз к объему снэсн ( )

СН сС;

Эп = Х-ЙП , Зл. = -'- •

^ о; ' I

р".. 4 4

/ « $ ■ I • * ; 9 =

0 Ь <1 1

р. <у

В результат© интегрирования уравнения"двизенйя, записанного дяй смеси в ц«уюы в форме уравнения Гроиеко-ЛамЭе

х /v . .

(Ь (- +7> X *оЬ \ +-

(7)

¡7Т ■> •> * ¿77 Ьч 3 Ь (¿1 0) )

+ ТКДу'й + v•c^ и ^ 2_1е.уа-21 р -^-Оо-о ;

<1 * О) 1 у 4 ; 4 4-1 /> 4 > ' '

= 4}.-и. \ 4 =1.2,...,//)

4 М ,/ * / ' * /

получены интегралы Бэрнухлп, Гроаеко, Ксги-Латрад-м, Эйлера, ксторуе в обцек виде записи ггзтся так: -

- 16 -У

А

= сопзt

где введено обозначани© пготноет:: м9ханвч9СК.оД оцерзд на нищг массы к потенциала скорости ддщ ¡,коуяонентов еиася.

4 г. — Ъ-О-

I. - 4 ЛУ. + & + ср • ■ в ^(р .

Езлзз конкретизировать теруодкнсигч&скоо соатоянио j -х компонентов

А'

= аолв!

т^ = оопв£

р . С/>

*

4 »

Л

С = сх>п$Ь

= соп&1:

К = aonst .

то интегрМ Коаш^Еангранжа, например, будет иметь вид:

л— _

«О

-21 £' (к-1)~р

04:

ОЬ

(9)

О'-Л.г,/ , 6.'-i.fi.....п; 7 = л +

В связя с тен, что конкретизированная запись уравнений с учетва тепяо-иайсоебизнй Е гидродинамического взаимодействия п(Ш&11&нтов е^егл иагет быть осуществлена лишь на основе опиеа-нея ш: свойств для кенкргтекх фора двк&ения, во второй главе

.'Хг;:;.' нлябсяег ссзд-ые вопросы г?5Дро?ериодгз!Г^.з(ки паро-

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

В рээультате сок;эстного рассиотренпя 'дзигнгиия геометрических форд ксапонеитоз сн^си к проягходяцях пр*л атси флэкчзст кегх процессов вкязяшн егзз™ь осноеких стадий двухфазного потока:

Рис.{

Нулеидя стадия - однородной глд^ости с раствора п-""*!

в ной г*алой. 8та стадии длится до наступления момента К8<яжт?нея плдкостн гг-зом.

Перзад стадия - ссрдяя гозпегэгаяенкя пуз»!р.ъ»»г (ггяегкк зг-рпдоетгй) и ст&тичвска устсРлызая ,*х раширзнне. Втора I стадия - течений смеси пря статически иеуетойчягом рес^гг-рзння газовьх пузырькоц. 0

Т^-втья стадия - стадия пузырькового випекия.

Четвертая стадия - течение пара-газовой смеси с каплями кид-постн (капельная структур« потека).

Пятая стадия - стадия статически устойчивого сжатия парога-

оовых пузырьков,

* <.

Шестая стадия - захлопывание пузырьков. За пределами шестой етеднн потор вновь становится однофазным. 0

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

В это! са глав© рассиотран ряд вопросов по определению баро-«ропной зокгекмбатн для шогефавной ездиоетн о учетом фазовых превращений.

Доказано, что ееяз компонента баротропны, то и смесь в целом баретрогаа

/ г>.

= —ТГТТ

С, Сп. \Р. ) ' *

1 ^ \ 0, / 4 к »

1

А/ /v /р

П. ; т) —показатель политропы и удельный поток фаз' еих извращений ¿-го - компонента смеси.

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

где введены обозначения.

О) й ' / < О) •

у - V /3 _ ¿¿.

У

А р рЬ I ( -с;. •>; , 1.4 /

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

д ч йР ) < а>

а) /

* г с с* с.

о.2 = 2 ^а* + Р~йР~

4

( И )

ГД" « _ с/Р-

.--/—I

' ¿р.

■I 4

ЛР

а = 6* <1р

Если пренебречь влиянием сил поверхностного натяжения, то скорость с»укэ для двухф<з?!!ой жидкости будет равна

г

а =

а а2

I л

'« со

1 а2 + с а2 + а*аг(р- р)-^

1 г. I П'^Гп

сГР

oLP P ' . ' n-(i-î)

Пршар расчет«:скороети звука воздуховодяной смеси и сопоставление е экспериментальными данными С.Н. Костерина и Л^И. Гуно-ш представлен на рис. 2.

Зависимость скорости звука для поздуховодяной смеси от

- объзмной концентрации воздуха (С4 ) при Р = iû "кг/л.1

5 6*« 0 .}

— рлсчзт

д- экепаржел^алышо данино С.Н. Костерит с - вхсперхмэдъ'елыше даннш Л.И. Гурова

Рио. 2

Hos вавестно, доказательство равенства скорости звука скорости разгфсегрвнвния ехабых возмущений в потоко однородного газа дала EossiasHOCTb понучить ргзення ряда важных задач с газодинамике. Естественно, что не м-*н®о ваге« этот вопрос и для исследо-

нгн'.щ пя многофазной ™:д:£оз*П!. Пря анолопгсжа: прэдполозз-sraax, обдоно прздкмаггздс пр:з ибатэдовмсш точсгпкх огородного раза, бнли рассмотрены попроси рэдпрострзнсння сла&гс позг.^цо-Hiiíi в потока гзюгсфагнся гздкссти. Используя диис^пгсбсяо уево-зяя, кяшипяея гцдродпнмггеззяио урасненпя гщвалъной яцдЕовто пря отсутствии теазо-^гвсообхгна, и к&неюяшгссхвэ уз-ловкя езда

fv^J =Л..П ;

ЭЛ.. _i¿_

■Я. V •

i, F >

(¿ =i.2 j =-i,2(. ,.,y)

полезна система ypar.nc::::ü, спр?дзл.тггл сзорозть раепровтрсмэ-н!1я поверхности елабнх рззртгов яэ отнетгнхэ н j -ну ксггяэ-ненту с?:зе:«

Col п = QXC.a *Р.Я ;

V/ J J p¡ Jj a' '

lb ¡

Pj j J J vé

(12)

¿■да Ô ■ = // - ■O' i N - скорость пзсзыгягнкя пэ'зрхног-

J F ■ nj F - r _

ça слабоpo разрыва, общз!5 дгш ввех погока. п - одм-

изссгЯ Езггеор иоршя. СквФаа урагясякЯ (12) э об^ез елу4ао яэ~ ;1г:о?зя п дот 50 зоззгягняя исоахада!1 двгазшая гапззь

¡трсязстЯ, опсздвХйгг^а дпнеючссзсэ вэсаг:бдвае*кгэ веяссксятов Однсло е ncîîo^sa подегсноваа

л « я . - ■ (î3î

с j p¡ a Pj

коадрая следует па тоздестИ.

dC¿ = dCj_ dP dt ~ dP- dt саотсаа уравнений (12), (13) ршана в общзи вода.

А*

<

ML.ÁIl

■ 1 *

к. с.

4 4

d Р.

4

а.

3 d ci { -+■ —'-i

C.dP.

4 I

(14)

И для сиэеа в цзяогз

n

i'i

Л ¿C,

a:

4

K. dP; i ¿

= í

(15)

Q = d)

\ Pj '

Sl

a*

. 4

A

K. dP.

'i

(16)

8 -Q.+ U

Сi) 4 "J

- U.n U .= , ni Q

(i)

Ив еопоатавдегш; фарауя (II), (15), (Гб) ввдю, что в потоке ецдгозгл ( С. - I ) скорость звука н скорость раеирзотрзнв-

4

игл сзвэрзшоетк елаСш раэрегзов рш:и, с в кнегефазно г язтето она йяшкы ко ьначашга, т.о. В — й

В этей us глав© рассмотрено течении еш; пр« каше а скяь--тк уззр^тов. В частности иолучеки i вдродкна« чесуче; '.о^неяе-

ния на прсшом скачка для гдогофазной газдкоета, которые анагегач^ пи урэг-и&киш Рэшгкна-ГЪгсньо.

В четвертой глоео р^се::отрэ«о з одноксргеп прпб.тг^згнпп ета-штокарное двпгвняо дсухфаэной иадкозти, в результате чэго для частных случаев урашекш движения прокитограрзее'.ш а учете» фа~ эобнх превращений. Тая, используя барэтрогсгнэ оап::с:п:о5тл (9) пря , проинтегрировано урзшгнкэ двкаеняя слодуэ-

сС. = аопзЬ , цето взда.

¿Р

С01 Р<*

Р я

сС.

п АР

"ог ры I

(I

'ГО'

(17)

\

с - ^

поя

< i

пэ урзгпсткгд (17) ояо—

■и

<31

" Р0

1-Х.

2

Ро

_9х_

1-

•р ^

Р

(10)

'2

<У Я . Р

_— . _ —р_ + —а_ ^ —

О 1 л

Л

Л

А' .0 о

СЯО)

Г

Излученный Л. Прандтлоз: интеграл стационарюго течения сыесн ео-•ды и воздуха ,яшша®ея ooisiiai случаев интегралов £9) шш (19) к

!ii30st вод! '

у/ р

— + (l-s)- —— + 6-кСпР co^sf . 2 J ;

Здоеь <5 - отнегхнна мааса воздуха к кассе вода. Предполагается,

о

что К = ВТ ~ const , £ = const и иалс по бз-

¿itcsig. Интеграл (205 полностью совпадает е аналогичный интеграле:.: дай однородного газе, так каш np:i аго^ ко изменяете« состав' кяоя, «.«з. пс иаиеяюгея объемна ионцонтрацт-ш ее компонентов.

Пря еоог&з-гйгьукрх допущениях урздааасю движения возказяо щоагавррзрэва^ь .в бон наяоаенш условий о барзвропности компоне-юва u eiisssi в цааом. Ток, нелршэр, для однсиорного потока вда-сяьаоВ кногефаэмой зцдкости при отсутствия фазовых превращений , шзеоЕах сна а при постоянство площади поперечного еттт потека ( f я const ■ ) из уравнения (2) следует

и*; *<«>** ■ (an

В «syuia двукфзлшоП шщкоете шданоеть css&iaeuoro компонента и Escsaossfb ease: sxkskj ьирааить тгж вшшзозше кии объекы, ко-■ чзклз ©прадеякэдагя оОъкзшал кенцепграцшглз, скоростью и гвомотра-сй потока. Si»o обегоягельетво позволило прокнтогрярэвать уравне-газ £17) боа использования барогропных ваЕясшоетей для компонзн-i-га eyssa в целей. В чаетноето ураьнснио (17) проинтегрировано на оодения фуишри, олрздедгс^иЯ геомотрЕШ потека, ■

f « f(i--Urc); nir = " ;

J V 0. a / » о С -С " С

е наяозениеы условий

"1^ = 0 при j = i ] "^ = 4 при l/=f0. В результате интегрирования а поадвдукцшс преобразовашй получена аледующая формула

(22)

Пря j = i , т.е. joV = const и при ■О.-О' формулы (21) и (22) тоэдественны. При ft = const цз формулы (22) следует известный интеграл Р + = const . Зорцулы (21)

и (22) киек>т болыаое практическое значение для разработки теоретических основ гидравлики двухфазной падкости. Так проинтегрированное уравнение двигения (21) и чпаяа РэЙнояьдса (6) для многофазной. жидкости позволяет1 установить связь мацду скорость» движения ее компонентов н потерями на трзнкэ, например, о трубах тайно , как это сделано для однофазной шщкости на основе уравнения Бер)уллм.

d -г ы а р -у*''!

г т г '

п г

.V

he d3 J_ 2$~Г Re' '

.....*

О)

- скорость базового компонента югзеи. По этим фориулгл !jcsho построить экспериментальные з&висшюста, гшалопгаш® известным графикам Никурадзе, для коэффициентов гидравлического

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

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

Е »той яе главе рассматривается также не менее важная задача гщфавлипи- определение гидравлических потерь пргл резко». и

и постепенном расширении двухфазного потока. Предполагается, поток является одномерным, давление скорости компонентов равны, причем давление при резком расширении потока изменяется скачком от £ до Р . Применяя уравнения сохранения массы, импульса и уравнение (22) с учетом потерь на "удар"

(25)

получена формула, аналогичная формула Бо^це-Карно

¿р л

+ £ •

I = и

и

(26)

¿„-О

А f

А

* -

А р

А р

при V < -Оа при 4У> V

Для звукового потока

чу

при р=рв- const.

Р>Ро ^

поато-

0

п потери капора на резкое раезярсние будут

б'оя'ьа®, чем для несжимаемой гэдкооти. Это объяснязтся тем, что в двухфазном потока пря переходе от сочснля f к сечении $ скорость уменыаается но только за ечет раеяирэния площади попе-рёчного сечения потока, но и за счет увзямче'ния плотности сиеея при повьсении даатснкл. Полученные формулы позволяет такгэ определять коэффициент смягчения удара Я> при постепенном раета-рении двухфазного потока. Результаты обработки экспериментальных

- 2В -

данных Хенча и Джонстона по течению водо-воздушной смеси в пряыо-о

угольном диффузоре представлены на рис. 4

значения коэфф. "полноты удара" для потока водо-воздушной сыеси в прямоугольном д:|ффузо-рз

1 - С2 = 0

2 - С2 =0,33

3 - С2 = 0,55

4 - Сг = 0,75

Рис. 4

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

В пятой главе даны в одномерном приближении аналитические и пелузмпирические решения вадач по внутреннему течению двухфазной жидкости. Решения задач по истечению из солла идеальной д£у-„,ф«1з-ной еидкости без и при наличии фазовых превращений построен» на основе прокнтегри1\)ванного уравнения движения (16), в результата

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

постоянство плотности одного нз компонентов смеси шазт вид .х/п -1

П - 4

; (Р • * Сф) ■

Л =

с* л * Сагв

сип**си в*

* р п =4-4-

в «

п -1

4 -

при п *{ ; В "1п —— при п = 4

а.

в =

У

о-

Л*'

о о'-

К

С + с

н- i

при п= к

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

На основе уравнений сохранения массы, анергии н ур&енення двизения (25) получены решения задач дозвукового и сверхзвукового некзсантеропического истечения гатяцей падкости. Прячем давление даде идеальной еыесн (без учета днссипативных с:и) проке» ходит с повышением энтропии, что объясняется перераспределения

тепловой энергии.мваду компонентами смеси. Чем больше масса фа-о

эовых превращений, тем больше возрастает энтропия. Максимальное количествр-#ара;образуется при 0 = 0 и при этой увеличение энтропии (г«счегмаксимально. Скорость же истечения кипящей жидкости будет ыаксииадвдой.прн минимальном изменении энтропии смеси. Так как дляг.изолированной системы dS >> 0 , то изменения скорости движении «ИгЦНТропии смеси находятся в пределах:

(S =-const)><Ks) > V(SmaJ = 0 ;

s0 < s < з(о = о) (27)

(s - Sjb - S2p3 ; SSoipoi + So3pM)

Известно®-■реиение задачи иэоэнтропического истечения кипящей гид-кости, полученное на основе уравнений сохранения массы и энергии и накладываемого условия о постоянстве энтропии смеси ( S = S0-= const ), удовлетворяет только максимуму в пределе (27), а полученные решения задачи удовлетворяют всему диапазону возможных изменений параметров потока, ограничиваем« пределами (27). При этом рассмотрены три режима истечения кипящей жидкости.

- Дозвуковое истечение насыщенной кипящей жидкости при уело-внях Р0 = РН ; Р*Р46),>-Р*;

Здесь Р0 , Р , - давление в начальном, критическом и в

конечном сечении потока. Рн - давление насыщенного пара.

- Сверхзвуковое истечение недогретой жидкости при условиях

Р' >Р

О Ьы X •

Решение задачи значительно упрощается при допущении, что Р - Р^ Такое допущение не вносит существенной погрешности в решение за-

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

- Сверхзвуковое истечение насыщенной яидкости, т.е. когда вы-

Это наиболее сложный случай, так как нужно определять параметры в критическом сечении потока. Использование уравнения движения (22) и формулы (II) для скорости звука позволяет найти приближенное, но достаточно точное и сравнительно простое решение задачи. Сделанное в работе сопоставление с опытными данными других авторов (Корасвв Э.К., Миронов D.B., Перемыщев В.В., Фуске Р., Генри X; и др.) по критическому расходу кипящей яидкости подтверждает вполне . достаточную точность приближенных решений для прикладных задач.

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

с

полняются уело г ш: PQ = Рн

Ьых '

нений.

Уравнения неразрывности и движения для смеси

г/ ¿Л 1 * О- с1 X з ах

Уравнения неразрывности к движения для второго компонента Смеси Сг ¿к Кг ¿х + 4)л ¿к * } ¿х 1якл >

(30)

¿X • ¿Л ¿х

, (Я - Я)

Уравнения притока тепла для с&есч и баланса тепла на границе раздела ое компонентов •

* * «у. сллУ* к ; «>

' * ^ сзз,

Уравньаия терюдннаитеескогс состояния и динамического вэаимо-двйстшя ко&понентой смзси (уравнение Рело'я)

+ А* Пл. + А п . к

,, б,

+ . (35)

Здесь £ - отношение динамических к статическим силам, действующим на пузырзк

Уравнения (28) - (36) записаны в безразмерно:.: виде: безразмерные констакты. К' основной сиСТеыз уравнений добавляется Уравнения, определяющие' коэффициент гидродинамического еопро-ШЪлеИМ пузырьков С^ Сдавленно насыщенных лоров Рн , реометр.® потока $ (для' прямой задачи) и распределение давления (для обратной задачи). Причем линия распределения давления при решении обратной задачи'долина пройти через точки задаваемого давления на входе и выходе * йэ' сопла и через точку критического сечения, в котором &--Р и (1} = 0 .0 отличне от Мвестиого решения аналогичной задача истечения однородного газа,. 0И гИтЬ^сЬо Р обычно является известной константой (для во-з'духйР = 0,528) в данной' случае значение критического давления неизвйетй'о и определяется' й результате решения задачи.

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

£Эб) выполнялось как линейным (методом Эйлера, Рунге-Кутта, Адамса), тая и нелинейным экспоненциальным методом. Практика выполнения расчетов показала, что при решении таких задач более выгодно применять нелинейный метод, причем с пос-

- з -ч

тоянным шагом интегрирования ( А Х -10 -т- iO j

О íX éi ), так как выигрыш времени за счет переменного шага обычно меньше затрат машинного времени на проце -дуру выбора величины шага.

М ; а ; О

Г : Т i '■Í

i,5

1,0

. 5

/

з i

2

i ^

' S • 74 •

_ У

0,95

O.S 0.9 а

i-а: г-i: з-м ю"х; ч-i\ - s-т \а = —

а

0,5

1-Pt. :2-PH ; з-¡:ч-'>а:5-<?:б-сг.

з г i о -i -г

\

\

\

J

о о.г

0,6

Рис. 5

Расчет истечения из сопла паро-воздухо-водяной смеси

Значения функции £ определяется из уравнения (36) о запаздыва-ниеы' на один шаг и в случае необходимости с повторным уточнением значений всех функций. Один из примеров расчета показан на рис. 5, из которого видно, что в зоне резкого возрастания объемной концентрации газа (пара) динамические силы в несколько раз превышают статические силы, действующие на паро-газовые пузырьки, причем в этой зоне функция § имеет экстремальные значения обоих знаков. При резких расширениях потока, а также при резких пиках тепловой нагрузки ^ х значения функции £ изменяются в пределах нескольких порядков. По рассмотренной схеме решены прямые и обратные задачи течения в соплах и в контуре атомных реакторов двухкомпонентной и трехкомпонентной двухфазной жидкости при пузырьковой, капельной и пузырьково-капельной структуре потока:

- Течение двухфазной андкости при капельной структуре потока

- Течение двухфазной аидкости при однородной пузырьковой структуре потока

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

- Течение двухфазной жидкости при переменной пузырьково-капельной структуре потока .

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

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

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

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

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

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

величины гидравлических потерь, габаритов и мощности реактора

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

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

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

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

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

о

гих параметров теплоноснтедл при разрыве п&тюпроеода первого кон-

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

- Получена система гидродинамических уравнений для многофазной жидкости при наличии переноса массы, импульса и тепла с записью уравнений движения в форме уравнений Громеко и Навье-Сто.-.са.

- Получены бсэраэмерше параметры и критерии подооня для кно-

гофазиой жидкости.

-' Получены в обобщенном и конкретизированном виде первые интегралы уравнений движения, аналогичные интегралам Бернулли, Гро-меко, Коши-Лагранжа, Эйлера. Установлено, что первые интегралы уравнений движения существуют только при постоянстве массовых (весовых) концентрациях компонентов смеси.

- 3 результате рассмотрения вопросов динамики многофазной жидкости получены зависимости для определения скорости звука и чисел Маха в многофазном потоке, определены скорости распространения слабых разрывов, гидродинамические соотношения на поверхности сильаьх • разрывов и условия на прямом скачг.е (уравнения Рзнкина-Гюгоньо).

П. Разработаны теоретические основы гидравлики двухфазной жидкости:

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

- Решена задача по определению потерь "на удар" при резИом и постепенном расширении двухфазного потока. Показало, что гидравлические пот«фи "на удар" значительно больше в двухф&чном, чем в однофазном сотоке.

Ш. Получено в одномерном приближении решение ряда задач по стационарному течению двухфазной жидкости:

- Истечение из сопла идеальной двухфазной жидкости без и при наличии фазовых превращений.

- Неизоэнтропическое дозвукозое и сверхзвуковое истечение кипящей жидкости.

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

- -

превращений для капельной, пузырьковой и пузырьково-капельной неоднородной структуры потока.

I

V. Разработанные вопросы теории и полученные - решения задач использованы при наследования движения двухфазного теплоносителя в ядершх энергетических установках.

о

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

- Разработан комплекс программ для ЭВМ и выполнены систематизированные расчеты (численные эксперименты) различных вариантов праектируешх ядерных рзакторов большой единичной мощности.

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

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

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

- В области приборостроения для расчета параметров двухфаз-иого потока в измерительно-таркровочных устройствах. Кроме того, разрабо^йя прийор для измерения объелшых концентраций смеси с 7,сухфйомо:1 потоке.

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

экономический эффект.

Работа выполнялась с 1964 г. по 1973 г. в Институте гидромеханики АН УССР, а с 1973 г. - в Киевском госуниверситете на кафедре вычислительной математики факультета кибернетики.

Основные разделы диссертации докладывались и обсуидались на Ш и 1У съездах по теоретической и прикладной механика , на семинарах и конференциях по аэрогидродинамике, теплообмену и массообмену, на соответствующих семинарах Институтов гидромеханики, механики, иата-матихи, технической теплофизики, электродинамики АН УССР, Института проблем механики АН УССР, НИИ механики при ЬТ, Института тепло-п массообмена БССР, ЭНИН им. Г.М. Кржижановского, на кафедра аэродинамики н теплообмена КРУ, на кафедре газовой и волновой динамики МГУ и в других научных организациях.

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

1. Бабицкий А.Ф. Гидродинамические уравнения для парогазовой смеси. В кн.: Гидродинамика больших скоростей, выл. 3, Киев 1967,

с. 63 - 70.

2. Бабицки^ А.Ф. Возникновение паровых пузырькоз при кавитации и кипении - В кн. Теплообмен и гидродинамика в двухфазных средах. Киев I9C7, с. 102 - 106.

3. Баб1цькйй А.Ф., 1вченко В.М. Перш1 1нтеграли р1внянь .руху багатофаэно! £1дйня з прихованою к1нетикою - ДАН Укра1нсько1 PCP, сер1я А. ф1зйко-техн1чн1 та математичн1 науки 1967, № 12, с.1091 -- 1092.

4. Бабицкий А.Ф. Модель и гидродинамические уравнения многофазной лидкости при наличии фазовых превращений. - В кн. Гидродинамика, вып. 15, Киев, 1969, с. 80 - 89.

5. Бабицкий А.Ф. Одномерное движение двухфазной жидкости по

трубе переменного сечения. - Прикладная механика, 1970 т.6,вып.7, с. 107 - 112

6. БабЗдький А.Ф. Перш1 1нтеграли одновим1рних р1внянь нерозрив-hoctI I руху для двофазно! р1днини. - ДАН УССР, сер. ф1зико-тех-н1чн1 та математичн1 науки. .., 1971, JP7, с.615-618.

7. Бабицкий А.Ф. Одномерная задача истечения двухфазной жидкости из сопла. - Прикладная механика, 1972, т.8, вып.10,с.88-93.

8. Баб1цький А.Ф. 1нтегрування р!внянь руху баротропно1 двофазно! р1дини,- ДАН Укра1нсько1 PCP, сер. А, ф1зико-техн!чн1 та математичн1 науки, 1972, № 8, с.724-726.

9. Бабицкий А.Ф. Слабые разрывы в потоке многофазной жидкости,- В кн- Аналитические и качественные методы теории дифференциальных уравнений. Издание Института математики АН УССР, Киев,

1972, с.29-37.

10. Бабицкий А.Ф., Кравченко А.М. Измерение методом кондукто-метрии объемных концентраций компонентов в двухфазном потоке. -В кн.: Гидромеханика, вып. 20, Киев, 1972, с.92-95.

11. Баб1цькнй А.Ф. _ОдноЕш1рна задача вит1кання 1з сопла дво-'Ьазно! р1дин!' при наявност! фазових перетворень.- В кн.: Край о в! эадач1 г1дроаеромехан1ки, 1нститут математики АН УССР, Ки1в,

1973, с.75-85.

12. Бабицкий А.Ф. Соотношения на прямом скачке уплотнения (уравнения Рэнкина-Гюгоньо) для многофазной жидкости.- В кн.: Аналитические методы теории фильтрации и уеплопроводности. Издание шститута математики АН УССР, Киев, 1973, с.61-69.

13. Бабицкий А.Ф. К вопросу об истечении кипящей жидкости. -

В кн.: Теплофизика и теплотехника, вып. 25, Киев, 1973, c.II5-iI7.

14. Бабицкий А.Ф. Формы и стадии движения паро-газо-Еодпных смесей,- В кн.: Гидромеханика, вып.26, Киев, 1974, с.40—49.

15. Бабицкий А.Ф. Превые интегралы уравнений движения для многофазной жидкости - В кн. Вычислительная и прикладная математика, вып. 26, Киев, 1975, с.48 - 55.

' 16. Бабицкий А.Ф. Определение сопротивления движению двухфазной • жидкости - В кн. Гидравлика и гидротехника, Киев 1975, с. 27 - 32.

17. Бабицкий А.Ф., Белоусов Б.Н. Численное рещение задачи истечения из сопла капельной двухфазной жидкости - В кн. Вычислительная и прикладная математика, вып. 27, Киев 1975, с.5£ - 58.

18. Бабицкий А.Ф., Белоусов Б.Н., Численное решение задачи истечения из сопла пузырьковой двухфазной жидкости при наличии фазовых превращений. Институт математики АН УССР, Киев 1976, с. 14 - 23.

19. Баб1цький А.Ф. 0собливост1 кТнематики багатскомпонентного суц1льного середовища. В1сник Ки1вського ун1верси!ету, сер. математика та механ1ка, № 18, 1976, с. 81 - 85.

20. Ляшко I.1.,Баб1цький А.Ф. Про розвиток теоретично1 г1дродина-м!ки багатофазноI рГдини. В1еник АН УРСР, 1? 4, 1976, с. 5 - 13.

21. Бабицкий А.Ф., Кохановский И.Н., Демяненко С.Н. Численное решение задачи истечения из сопла пузырьково-капельной двухфазной жидкости - В кн. Вычислительная и прикладная математика, вып. 32, Киев, 1977, с. 10 - 17.

22. Бабицкий А.Ф., Белоусов Б.Н., Высочанский В.В. Численное решение задачи истечения трехкомпонентной двухфазной жидкости при наличии фазовых превращений - В кн. Вычислительная и прикладная математика, вып. 35, Киев, 1978, с. 7 - 15.

а)