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

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

Saved in:
Bibliographic Details
Date:2010
Main Authors: Шараевский, И.Г., Фиалко, Н.М., Зимин, Л.Б., Меранова, Н.О.
Format: Article
Language:Russian
Published: Інститут проблем безпеки атомних електростанцій НАН України 2010
Series:Проблеми безпеки атомних електростанцій і Чорнобиля
Subjects:
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/58200
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока / И.Г. Шараевский, Н.М. Фиалко, Л.Б. Зимин, Н.О. Меранова // Проблеми безпеки атомних електростанцій і Чорнобиля: наук.-техн. зб. — 2010. — Вип. 13. — С. 69–75. — Бібліогр.: 20 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-58200
record_format dspace
spelling nasplib_isofts_kiev_ua-123456789-582002025-02-10T01:41:03Z Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока Проблеми прогнозування кризи тепловіддачі в дисперсно-кільцевому режимі течії двофазного парорідинного потоку Problems of forecasting of the heat transfer crisis in annulad-disperse mode of the diphasic vapour-liquid flow Шараевский, И.Г. Фиалко, Н.М. Зимин, Л.Б. Меранова, Н.О. Проблеми безпеки атомних електростанцій Рассмотрены пути и доказана возможность прогнозирования кризиса теплоотдачи 2-го рода. Розглянуто шляхи й доведено можливість прогнозування кризи тепловіддачі 2-го роду. Ways are considered and possibility of forecasting for 2 type heat transfer crisis is proved. 2010 Article Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока / И.Г. Шараевский, Н.М. Фиалко, Л.Б. Зимин, Н.О. Меранова // Проблеми безпеки атомних електростанцій і Чорнобиля: наук.-техн. зб. — 2010. — Вип. 13. — С. 69–75. — Бібліогр.: 20 назв. — рос. 1813-3584 https://nasplib.isofts.kiev.ua/handle/123456789/58200 532.529 ru Проблеми безпеки атомних електростанцій і Чорнобиля application/pdf Інститут проблем безпеки атомних електростанцій НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
topic Проблеми безпеки атомних електростанцій
Проблеми безпеки атомних електростанцій
spellingShingle Проблеми безпеки атомних електростанцій
Проблеми безпеки атомних електростанцій
Шараевский, И.Г.
Фиалко, Н.М.
Зимин, Л.Б.
Меранова, Н.О.
Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
Проблеми безпеки атомних електростанцій і Чорнобиля
description Рассмотрены пути и доказана возможность прогнозирования кризиса теплоотдачи 2-го рода.
format Article
author Шараевский, И.Г.
Фиалко, Н.М.
Зимин, Л.Б.
Меранова, Н.О.
author_facet Шараевский, И.Г.
Фиалко, Н.М.
Зимин, Л.Б.
Меранова, Н.О.
author_sort Шараевский, И.Г.
title Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
title_short Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
title_full Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
title_fullStr Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
title_full_unstemmed Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
title_sort проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока
publisher Інститут проблем безпеки атомних електростанцій НАН України
publishDate 2010
topic_facet Проблеми безпеки атомних електростанцій
url https://nasplib.isofts.kiev.ua/handle/123456789/58200
citation_txt Проблемы прогнозирования кризиса теплоотдачи в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока / И.Г. Шараевский, Н.М. Фиалко, Л.Б. Зимин, Н.О. Меранова // Проблеми безпеки атомних електростанцій і Чорнобиля: наук.-техн. зб. — 2010. — Вип. 13. — С. 69–75. — Бібліогр.: 20 назв. — рос.
series Проблеми безпеки атомних електростанцій і Чорнобиля
work_keys_str_mv AT šaraevskiiig problemyprognozirovaniâkrizisateplootdačivdispersnokolʹcevomrežimetečeniâdvuhfaznogoparožidkostnogopotoka
AT fialkonm problemyprognozirovaniâkrizisateplootdačivdispersnokolʹcevomrežimetečeniâdvuhfaznogoparožidkostnogopotoka
AT ziminlb problemyprognozirovaniâkrizisateplootdačivdispersnokolʹcevomrežimetečeniâdvuhfaznogoparožidkostnogopotoka
AT meranovano problemyprognozirovaniâkrizisateplootdačivdispersnokolʹcevomrežimetečeniâdvuhfaznogoparožidkostnogopotoka
AT šaraevskiiig problemiprognozuvannâkriziteplovíddačívdispersnokílʹcevomurežimítečíídvofaznogoparorídinnogopotoku
AT fialkonm problemiprognozuvannâkriziteplovíddačívdispersnokílʹcevomurežimítečíídvofaznogoparorídinnogopotoku
AT ziminlb problemiprognozuvannâkriziteplovíddačívdispersnokílʹcevomurežimítečíídvofaznogoparorídinnogopotoku
AT meranovano problemiprognozuvannâkriziteplovíddačívdispersnokílʹcevomurežimítečíídvofaznogoparorídinnogopotoku
AT šaraevskiiig problemsofforecastingoftheheattransfercrisisinannuladdispersemodeofthediphasicvapourliquidflow
AT fialkonm problemsofforecastingoftheheattransfercrisisinannuladdispersemodeofthediphasicvapourliquidflow
AT ziminlb problemsofforecastingoftheheattransfercrisisinannuladdispersemodeofthediphasicvapourliquidflow
AT meranovano problemsofforecastingoftheheattransfercrisisinannuladdispersemodeofthediphasicvapourliquidflow
first_indexed 2025-12-02T13:32:06Z
last_indexed 2025-12-02T13:32:06Z
_version_ 1850403535559065600
fulltext ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 69 УДК 532.529 ПРОБЛЕМЫ ПРОГНОЗИРОВАНИЯ КРИЗИСА ТЕПЛООТДАЧИ В ДИСПЕРСНО-КОЛЬЦЕВОМ РЕЖИМЕ ТЕЧЕНИЯ ДВУХФАЗНОГО ПАРОЖИДКОСТНОГО ПОТОКА И. Г. Шараевский, Н. М. Фиалко, Л. Б. Зимин, Н. О. Меранова Институт проблем безопасности АЭС НАН Украины, Киев Рассмотрены пути и доказана возможность прогнозирования кризиса теплоотдачи 2-го рода. Ключевые слова: кризис теплоотдачи 2-го рода, двухфазный поток, дисперсно-кольцевой режим, тепловыделяющие сборки, надежность и безопасность АЭС. Для реакторов кипящего типа (РБМК, BWR) наибольшую опасность представляет кризис теплоотдачи 2-го рода, соответствующий началу перехода от дисперсно-кольцевой структуры потока (рис. 1, в) к дисперсному (рис. 1, г). Физика кризисных явлений в дисперс- но-кольцевом режиме течения двухфазного потока принципиально отличается от теплогид- равлических процессов, сопровождающих трансформацию структуры двухфазного пристен- ного слоя при переходе пузырькового кипения к его пленочной форме в потоке недогретой и насыщенной жидкости (рис. 1, а, б). В дисперсно-кольцевом потоке жидкость течет в виде пленки по стенке канала и в виде капель, распределенных в потоке пара в центре канала. При давлениях менее 5–7 МПа капли воды в основном сосредоточиваются в пристенной зоне, а при более высоких давлениях находятся в центральной области потока [1]. Рис. 1. Схематическое изображение основных физических закономерностей механизма кри- зисных явлений в цилиндрическом парогене- рирующем канале при различных режимах восходящего течения двухфазного потока, а также профилей истинного паросодержания ϕ при наступлении кризиса теплоотдачи (К – сечение кризиса): а – недогретое пузырьковое кипение; б – на- сыщенное пузырьковое кипение; в – дисперс- но-кольцевой режим течения; г – дисперсный режим течения. В отличие от кризиса теплоотдачи при кипении недогретой и насыщенной жидкости (см. рис. 1, а, б), физика кризисных явлений 2-го рода в дисперсно-кольцевом режиме течения двухфазного парожидкостного потока, согласно терминологии [2], связана с пре- кращением расхода жидкости в пристенной жидкостной пленке (см. рис. 1, в). Величина этого расхода в соответствующем сечении парогенерирующего канала определяется интен- сивностью следующих тепловых и гидродинамических процессов (схематически изобра- женных на рис. 2): испарения жидкости – B = q/r; уноса влаги из пристенной пленки в результате пузырькового кипения в ее объеме ПE ; механического срыва капель с гребней волн МE ; выпадения капель из ядра потока на поверхность пленки D . При этом МП EEE += . Явление кризиса гидравлического сопротивления в парогенерирующем канале соот- ветствует области трансформации пристенной жидкостной пленки от состояния с мелкими волнами (рябью) на ее поверхности к микропленке с уже подавленным в ней поверхностным кипением и сопровождается характерными изменениями акустического спектра. И. Г. ШАРАЕВСКИЙ, Н. М. ФИАЛКО, Л. Б. ЗИМИН, Н. О. МЕРАНОВА ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 70 Рис. 2. Теплогидравлические процессы в дисперсно– кольцевом режиме течения при возникновении кризиса теплоотдачи 2-го рода по данным [2]: 1 – начальный расход (GПЛ. 0); 2 – испарение пленки (B=q/r); 3 – механический унос капель (EМ); 4 – пузырь- ковый унос капель (EП); 5 – выпадение капель (D); 6 – сечение кризиса теплоотдачи 2-го рода (хгр); 7 – область мелких волн (ряби); 8 – жидкостная микропленка с подавленным кипением; 9 – сечение кризиса гидрав- лического сопротивления (х∆∆∆∆Р). Для аксиальной координаты Z парогенерирующего канала с периметром Р и сечением А дифференциальное уравнение расхода жидкости в рассмотренной пристенной пленке имеет вид [1] )( BED A P dZ dGпл −−= . (1) Аналитическое решение уравнения (1) наталкивается на значительные трудности, связанные с отсутствием необходимой информации о деталях механизма тепловых и гидродинамических процессов, формирующих расход жидкости в пристенной пленке и, таким образом, непосредственно обусловливающих конкретные особенности кризисных явлений. С целью получения хотя бы качественного согласования полученных решений с экспериментом предпринимаются попытки (такие, например, как [3]) эмпирического под- бора констант, необходимых для интегрирования уравнения (1). В настоящее время аналитические исследования кризиса теплоотдачи в условиях дисперсно-кольцевого режима восходящего течения двухфазного потока ведутся в контексте двух основных подходов. В соответствии с первым из них главным объектом аналитического исследования является тот вид тепломассообменного процесса (см. рис. 2), который приме- нительно к соотвествующему участку и конкретным теплогидравлическим условиям в парогенерирующем канале оказывает доминирующее влияние на формирование пристенной жидкостной пленки. Поскольку кризис теплоотдачи второго рода обусловлен высыханием пристенной жидкостной пленки, то уравнение (1) подлежит интегрированию от того сечения парогенерирующего канала, которое соответствует началу дисперсно-кольцевого режима течения, до той его аксиальной координаты, в которой достигается GПЛ. 0 = 0, т.е. отсутствие расхода в ней. Известен ряд методов решения уравнения (1), таких, как [4], различающихся подходами к выбору граничных условий, а также моделями для определения параметров D и E. Последние требуют надежной информации о таких трудноопределимых параметрах, как текущая Ci и равновесная Cp концентрации капель, а также связывающий эти величины коэффициент пропорциональности K – коэффициент массообмена для каждого из процессов: iKCD = ; pKCE = . (2) Согласно данным [4], равновесная концентрация капель Cp определяющим образом влияет на расход жидкости в пленке, хотя эта взаимосвязь установлена авторами только для адиабатного потока. В рамках гипотезы остается также и предположение о равенстве коэф- ПРОБЛЕМЫ ПРОГНОЗИРОВАНИЯ КРИЗИСА ТЕПЛООТДАЧИ ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 71 фициентов массообмена К в соотношении (2). Кроме того, нет необходимой ясности и с определением двух составляющих параметра Е, который характеризует интенсивность уноса жидкости из пристенной пленки (см. рис. 2), а именно: пузырькового уноса ЕП, обуслов- ленного выносом влаги в паровое ядро генерируемыми в пленке паровыми пузырями, а также механического уноса ЕМ, связанного со срывом капель с гребней волн потоком пара. В настоящее время отсутствуют экспериментальные методы контроля указанных составляю- щих уноса жидкости из пленки. В силу этой причины не удается разделить влияние каждого из этих факторов уноса, что в свою очередь, не дает возможности уточнить и скоррек- тировать рассматриваемую математическую модель. Кроме того, модель [4] не учитывает возможности перехода при определенных тепло- гидравлических условиях пузырькового кипения в пристенном жидкостном слое к его пленочному режиму (факт такого перехода подтверждается экспериментальными результа- тами [5]), что может инициировать резкое снижение интенсивности теплоотдачи еще до возникновения условий пересыхания жидкостной пленки на стенке парогенерирующего канала. В этой связи приходится констатировать, что экспериментальные методы, позво- ляющие надежно идентифицировать в промежутках между гребнями волн (рис. 2) момент перехода пузырькового кипения к его пленочной форме в настоящее время фактически отсутствуют. Вполне очевидно также и то, что интегрирование уравнения (1) в контексте рассматриваемого детерминированного аналитического подхода предполагает наличие надежной информации о точке начала дисперсно-кольцевого режима течения. Вместе с тем, ввиду вероятностного характера дисперсно-кольцевого режима, равно, как и других типов структурных форм двухфазного парожидкостного потока, границы между ними являются физически нечеткими, а экспериментальные методы надежной идентификации этих границ до настоящего времени не разработаны. В соответствии со вторым подходом к аналитическим исследованиям кризиса тепло- отдачи в дисперсно-кольцевом режиме восходящего течения двухфазного парожидкостного потока решению подлежит замкнутая система дифференциальных уравнений в частных производных, включающая уравнения сохранения массы, импульса и энергии. Основное препятствие в решении указанной системы уравнений состоит в наличии внутри выде- ленного дифференциального объема границы раздела фаз, изменяющейся в соответствии с вероятностными закономерностями, характерными для процесса кипения в целом [6]. Такое противоречие между аналитическим детерминированным подходом к описанию данной структурной формы двухфазного парожидкостного потока и стохастической природой кипе- ния органически присуще и любому другому аналитическому подходу такого рода. Это относится и к наиболее совершенным математическим моделям данного режима течения, таким, например, как [7–9], в которых двухфазный поток представляется в виде взаимопроникающих и взаимодействующих континуумов, а дисперсно-кольцевая структура потока подразделяется на три характерных зоны: стенку канала, жидкостную пленку на ней и паровое ядро, причем тепломассообмен между этими зонами (процессы испарения пленки, осаждения и уноса капель, и др.) учитывается на их границах через условия сопряжения. Ввиду высокой сложности указанной задачи ее решение, а именно: определение расходов жидкости в пленке, полей температур, коэффициентов теплоотдачи, требует введения ряда упрощающих предположений. Последние связаны с введением комплекса замыкающих соот- ношений, которые должны приближенно охарактеризовать долю жидкости в ядре потока на нижней границе дисперсно-кольцевой структуры, коэффициенты диффузии и турбулентной вязкости, интенсивность капельного уноса и ряд других параметров. Естественно, что ука- занные упрощения существенно снижают точность получаемых аналитических решений. При этом с позиций общей оценки эффективности направлений аналитических подхо- дов к расчету кризиса теплоотдачи применительно к рассматриваемым условиям необходимо обратить внимание на следующее принципиальное обстоятельство. Достигнутое в ряде ра- бот, таких, например, как [10–12], удовлетворительное совпадение теоретически рассчитан- И. Г. ШАРАЕВСКИЙ, Н. М. ФИАЛКО, Л. Б. ЗИМИН, Н. О. МЕРАНОВА ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 72 ных значений критического теплового потока (КТП) с экспериментально полученными ве- личинами, в целом, не должно давать оснований для оптимизма при оценке продуктивности математических моделей кризисных явлений. Анализ этих работ свидетельствует о реша- ющей зависимости конечных результатов аналитического расчета КТП от выбора граничных условий, который, в силу ряда упрощающих допущений, сделанных авторами, далеко не всегда может быть признан физически обоснованным. Еще более сложны не только для аналитического расчета КТП, но также и для экс- периментального исследования этого важнейшего эксплуатационного параметра процессы теплообмена и гидродинамики в тепловыделяющих сборках (ТВС) энергетических ядерных реакторов. Как известно [7], Активные зоны (АкЗ) этих реакторных установок сформиро- ваны из ТВС с продольноомываемыми цилиндрическими ТВЭЛ и парогенерирующими кана- лами сложной геометрической формы. Для таких ТВС характерен интенсивный тепломас- сообмен между отдельными теплогидравлическими ячейками [13]. На начальном этапе экспериментальных исследований определение КТП в ТВС осуществлялось по среднепуч- ковым параметрам. Так был получен ряд эмпирических формул, каждая из которых, однако, была пригодна лишь для одной конкретной геометрии ТВС. Такова, например, получившая широкое применение в практике эксплуатации отечественных ядерных реакторов водо- водяного типа эмпирическая формула [14], которая, однако, позволяет рассчитать КТП только для одного значения давления и пригодна только для ТВС ВВЭР. Результаты дальнейших измерений температур и расходов по сечению пучка показа- ли, что энтальпия, массовая скорость, паросодержание в различных ячейках могут суще- ственно отличаться друг от друга. Поэтому было признано, что получение надежных зависи- мостей, построенных по осредненным по сечению пучка параметрам, даже в ограниченных диапазонах невозможно. Следующим этапом исследований стало рассмотрение стержневой сборки ТВЭЛ в виде совокупности отдельных теплогидравлически неравноценных подка- налов (теплогидравлических ячеек). Сущность такого поканального аналитического подхода состоит в определении (на основе анализа соответствующих уравнений сохранения) наиболее теплонапряженной ячейки, в которой КТП достигается раньше, чем в других. Наибольшая сложность подобного расчетного анализа состоит в необходимости адекватного учета взаимодействия между отдельными теплогидравлическими ячейками ТВС. При этом указанное взаимодействие обусловлено не только турбулентным обменом и конвективным переносом, но также и дрейфом пара в те подканалы сборки, в которых массовые скорости теплоносителя выше. Необходимо также отметить, что именно в двухфазном парожидкостном потоке перенос массы, количества движения и энтальпии между ячейками особенно существенен, причем вскипание теплоносителя в наиболее теплонапряженных ячейках значительно интен- сифицирует эти процессы. В этом контексте весьма важным обстоятельством следует считать и факт невозможности надежного обнаружения начала генерации паровой фазы на теплоотдающей поверхности ТВЭЛ существующими измерительными средствами [15]. Данные работы [16], приведенные на рис. 3, которые получены на основе гомогенной модели двухфазного парожидкостного потока, дают представление об интенсивности межка- нального обмена в функции от его массового паросодержания. Эти результаты свидетель- ствуют о том, что характерной особенностью межканального обмена в двухфазном потоке является наличие максимума его интенсивности в области перехода от пузырькового к дисперсно-кольцевому потоку при балансовом массовом паросодержании х ≈ 0,05÷0,15. Наличие этого достаточно широкого максимума является вполне очевидным свидетельством также и физически нечеткого (вероятностного) характера границы между двумя подлежа- щими идентификации диагностическими объектами, а именно: пузырьковым режимом течения в ТВС и его дисперсно-кольцевой структурной формой. В этой связи необходимо отметить, что, ввиду значительных трудностей, связанных и идентификацией основных классов структурных форм двухфазного парожидкостного потока в парогенерирующих ПРОБЛЕМЫ ПРОГНОЗИРОВАНИЯ КРИЗИСА ТЕПЛООТДАЧИ ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 73 каналах и, в особенности, в стержневых сборках ТВЭЛ, карты режимов течения в ТВС до настоящего времени практически не созданы. Рис. 3. Зависимость интенсивности межка- нального обмена в ТВС (Sth – безразмерный коэффициент обмена энтальпий) от массового паросодержания двухфазного парожидкостного потока по данным [16] (сплошная линия – при наличии дистанционирующей решетки; пунк- тирная линия – при ее отсутствии). Следует также подчеркнуть, что в широко известных теплогидравлических кодах «Совrа» и «Пучок» [17] теплогидравлические параметры ТВС с учетом вышеизложенного взаимодействия подканалов стержневой сборки осреднены отдельно для каждого из трех основных видов отдельных теплогидравлических ячеек (центральных, боковых и угловых). Дальнейшее развитие расчетных методик, основанных на учете межканального взаимо- действия в ТВС, связано с детальным анализом течения жидкостной пленки на поверхности ТВЭЛ, а также массообменных процессов между отдельными теплогидравлическими ячейками стержневой сборки, причем работой [18] начало таким исследованиям уже положено. Надежная информация об этих процессах создаст необходимую основу для аналитического решения уравнения (1) уже для каждого поперечного сечения ТВС с учетом изменения интенсивности массообменных процессов по периметру ТВЭЛ. Не снижая ценности вышеупомянутых теоретических исследований кризисных явле- ний в области высоких паросодержаний, следует подчеркнуть, что в настоящее время эти работы не позволяют с отвечающей требованиям практики точностью предсказать значения КТП, необходимые для оценки безопасных условий эксплуатации АкЗ ядерных реакторов. Таким образом, экспериментальные исследования кризисов теплоотдачи продолжают оставаться единственным источником информации об этом сложном и многоплановом теплогидравлическом процессе. Кроме того, с позиций важной для обеспечения надежности и безопасности совре- менных АЭС проблемы прогнозирования кризисных явлений на поверхности ТВЭЛ в реальных условиях эксплуатации АкЗ водоохлаждаемых энергетических ядерных реакторов необходимо отметить следующее. Во-первых, измерение ряда важнейших технологических параметров, например, массовой скорости, в каждой конкретной ТВС (не говоря уже о распределении этого пара- метра по отдельным теплогидравлическим ячейкам стержневой сборки) в настоящее время не может быть обеспечено с приемлемой для надежного расчета КТП точностью. Как изве- стно, контроль величины реального расхода в каждой ТВС в настоящее время не может быть обеспечен. Отсутствие информации об этом и ряде других теплогидравлических параметров АкЗ (например, в реакторах ВВЭР расход теплоносителя через АкЗ определяется исклю- чительно на основе приближенных оценок этого параметра и только интегрально, т.е. в це- лом по циркуляционной петле) не позволяет надежно оценить величину реального запаса до кризиса в каждой ТВС на основе используемых для этих расчетов зависимостей по КТП. Во-вторых, выполненные в ряде работ, таких, например, как [19], исследования влияния ряда физических аномалий в АкЗ ЯР, связанных как с отклонениями геометрии парогенерирующего канала от проектной, так и с появлением локальных всплесков тепло- выделения на поверхности ТВЭЛ («горячие пятна», связанные с распуханием топливного сердечника), свидетельствуют о существенном влиянии указанных факторов на величину И. Г. ШАРАЕВСКИЙ, Н. М. ФИАЛКО, Л. Б. ЗИМИН, Н. О. МЕРАНОВА ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 74 КТП, рассчитанную по известным зависимостям. Вполне очевидно, что указанные факторы реально не могут быть проконтролированы в процессе эксплуатации АкЗ ядерных реакторов. Так, в [19] исследовано влияние местного сужения кольцевого канала на КТП при давлениях 7...12 МПа, массовых скоростях (0,7…2,2)·103 кг/(м2·с), паросодержаниях от – 0,15 до 0,35. На основании выполненных экспериментов сделан вывод о том, что местное сужение канала приведет к существенному снижению КТП, и такое снижение тем больше, чем меньше массовая скорость. В [20] показано, что даже незначительный эксцентриситет стержневой сборки ТВЭЛ в пределах (0,6...1)·10–3 м может снизить величину КТП на 20...30 %. Кроме того, касаясь правомерности переноса результатов исследований КТП, которые получены на теплогидравлических стендах с экспериментальными электрообогреваемыми стержневыми сборками, имеющими обычно относительно небольшое (по сравнению с реаль- ными ТВС энергетических ядерных реакторов) число ТВЭЛ, необходимо отметить следу- ющее. Известно [13], что семистержневой модельный пучок имеет три основных типа эле- ментарных ячеек (центральная, боковая, угловая), которые существенно отличаются друг от друга в теплогидравлическом отношении. При переносе экспериментальных данных, напри- мер, с семистержневого пучка на реальные многостержневые ТВС необходимо учитывать, что для ячеек каждого из указанных типов в такой кассете распределение теплогидравли- ческих параметров иное, чем в экспериментальной семистержневой сборке. Выводы 1. Имеющиеся в настоящее время эмпирические расчетные зависимости для опре- деления критического теплового потока в области высоких паросодержаний (соответ- ствующей дисперсно-кольцевой структуре потока) не позволяют предсказать наступление кризиса теплоотдачи 2-го рода с точностью, необходимой для обеспечения надежности водоохлаждаемых ядерных реакторов. 2. В процессе эксплуатации ядерных реакторов оперативные данные о реальных значениях комплекса важнейших теплогидравлических параметров ТВС отсутствуют. Это также не позволяет на основе существующих расчетных зависимостей по критическому тепловому потоку надежно оценить величину запаса мощности до кризиса теплоотдачи на поверхности ТВЭЛ. 3. Налицо актуальная необходимость разработки принципиально новых методов оперативной шумовой диагностики предаварийных и аномальных режимов генерации паровой фазы на теплоотдающей поверхности ТВЭЛ в процессе эксплуатации ядерных реакторов. Эти методы следует рассматривать не только как дополнение к реализуемым в настоящее время расчетным методикам определения располагаемого запаса до кризиса теплоотдачи в ТВС, но также и как фактическую альтернативу использованию существу- ющих методик расчета критического теплового потока. 4. Возможность распознавания кризиса гидравлического сопротивления в парогенери- рующих каналах следует считать главным определяющим условием возможности прогнози- рования кризиса теплоотдачи 2-го рода, происходящего в области перехода дисперсно- кольцевого режима течения в дисперсный. СПИСОК ЛИТЕРАТУРЫ 1. Investigation of the liquid phase distribution in two-phase flow / V. M. Kascheev, P. L. Kirilov, Y. V. Muranov e.a. // Proc. Fourth Int. Topical Meeting on Nuclear Reactor Thermal-Hydraulics. - Karlsruhe, F.R.G. Oct. 10 - 13. – 1997. - Vol. 2. - P. 1157 - 1162. 2. Дорощук В.Е. Кризисы теплообмена при кипении воды в трубах. – 2-е изд. – М.: Энергоатом- издат, 1983. – 120 с. 3. Levy S., Healzer J.M., Adollahian D. Prediction of critical heat flux in vertical pipe flow // Nuclear Eng. Design. - 1981. - Vol. 65(1). - P. 131 - 140. 4. Hewitt G.F. Critical Heat Flux in Flow boiling // 6-th Int. Heat Transfer Conf. Toronto. - 1978. - Vol. VI. - P. 143 - 179. ПРОБЛЕМЫ ПРОГНОЗИРОВАНИЯ КРИЗИСА ТЕПЛООТДАЧИ ________________________________________________________________________________________________________________________ ПРОБЛЕМИ БЕЗПЕКИ АТОМНИХ ЕЛЕКТРОСТАНЦІЙ І ЧОРНОБИЛЯ ВИП. 13 2010 75 5. Козловский С.П., Пометько Р.С., Пашичев В.В. Кризис теплоотдачи и распределение жидкости в дисперсно-кольцевом режиме течения // Теплофизика высок. температур. – 1978.- Т. 16, № 1. - С. 94 - 98. 6. Толубинский В.И. Теплообмен при кипении. – К.: Наук. думка, 1980. - 315 с. 7. Каращев В.А., Муранов Ю.М. Расчет поля скорости в двухфазном дисперсно-кольцевом потоке // ТВТ.- 1976.- Т. 14, № 5. - С. 1015 - 1021. 8. Atwo-dimensional mathematical model of annular-dispersed and dispersed flows / P. L. Kirilov, V. M. Kascheev, Yu. V. Muranov a.e. // Int. Jour. Heat Mass Transfer. - 1987. - Vol. 30, No. 4. - P. I. - Р. 801. 9. Кащеев В.М., Муранов Ю.В., Юрьев Ю.С. Влияние замыкающих коэффициентов в двумерной математической модели парокапельного потока на результаты вычислительного эксперимента // Вопросы атомной науки. Сер. Физика и техника ядерных реакторов. – 1987. – Вып. 6 – С. 56 - 59. 10. Method of Critical Power Prediction Based on Film Flow Model COUPLED With Subchannel Analysis / Tomiyama A.Yokomizo, Yoshimoto Y a.e. // Jour. Nucl. Sci. Technology. - 1988. - Vol. 25. - № 42. - Р. 137 - 145. 11. Multi-fluid modeling of annular two-phase flow / T. Saito Hughes E.G., Carbon M.W. // Nuclear Eng. Design. - 1978. - Vol. 50. - P. 225 - 271. 12. The calculation of CHF in forced convection boiling / Whalley P.B., Hutehinson P. Hewitt G.F. a.e. // Proc. 5-th Int. Heat Transfer Conf., Tokyo. - 1974. - Vol. IY. - P. 290 - 294. 13. Шараевский И.Г., Домашев Е.Д., Архипов А.П. Метод локальной верификации карт режимов течения в тепловыделяющих сборках // Пром. теплотехника. – 2002. – Т. 24. - № 2-3. – С. 127 - 134. 14. Экспериментальное исследование и статический анализ данных по кризису теплообмена в пучках стержней для реактора ВВЭР / Ю. А. Безруков, В. И. Астахов, В. Г. Брантов и др. // Теплоэнергетика. – 1976. - № 2. - С. 83 - 88. 15. Делайе Дж., Гио М., Ритмюллер М. Теплообмен и гидродинамика двухфазных потоков в атомной и тепловой энергетике / Пер. с англ. – М.: Энергоатомиздат, 1984. – 424 с. 16. Молекулярный и турбулентный межканальный обмен в двухфазных потоках / А. П. Сорокин, А. В. Жуков, Ю. Н. Корниенко и др. // Препринт ФЭИ. А 1870. - Обнинск, 1987. 17. Нигматулин Б.И., Мелихов О.И., Соловьев С.Л. Состояние и развитие отечественных системных теплогидравлических кодов для моделирования аварийных и нестационарных процессов на АЭС с ВВЭР // Теплоэнергетика. – 2001. - № 3. – С. 17 - 20. 18. Subchannel analysis by the FIDAS Code based on the three-fluid model / S. Sugsvars, T. Sakai, Watanabe K., Rumens H. E. C. // Proc. IV Intern. Topical Meeting on Nuclear Thermal-Hydraulics, Karlsruhe. FRG. - 1989. - Vol. 1. - P. 560 - 567. 19. Исследование кризиса теплообмена при движении кипящей воды в кольцевом канале / В. С. Осмачнин и др. // ТФ-74. Исслед. критических тепловых потоков в пучках стержней. СЭВ. – М., 1974. – С. 251 - 266. 20. Влияние геометрического фактора на критический тепловой поток в пучке стержней // Тепло- массообмен / Т. Коберн и др. - Минск, 1976. - Т. 3, ч. 2. – С. 90 - 95. ПРОБЛЕМИ ПРОГНОЗУВАННЯ КРИЗИ ТЕПЛОВІДДАЧІ В ДИСПЕРСНО-КІЛЬЦЕВОМУ РЕЖИМІ ТЕЧІЇ ДВОФАЗНОГО ПАРОРІДИННОГО ПОТОКУ І. Г. Шараєвський, Н. М. Фіалко, Л. Б. Зімін, Н. О. Меранова Розглянуто шляхи й доведено можливість прогнозування кризи тепловіддачі 2-го роду. Ключові слова: криза тепловіддачі 2-го роду, двофазний потік, дисперсно-кільцевий режим, тепловиділяючі збірки, надійність і безпека АЕС. PROBLEMS OF FORECASTING OF THE HEAT TRANSFER CRISIS IN ANNULAD-DISPERSE MODE OF THE DIPHASIC VAPOUR-LIQUID FLOW I. G. Sharaevsky, N. M. Fialko, L. B. Zimin, N. O. Meranova Ways are considered and possibility of forecasting for 2 type heat transfer crisis is proved. Key words: сrisis of 2 type heat transfer, two-phase flow, dispersed-annular mode, fuel assemblies, reliability and safety of nuclear power plant. Поступила в редакцию 26.08.08