Исследование кризисов кипения и смены режимов многофазных течений в условиях, характерных для аварий на реакторах АЭС

Актуальность темы

.

В связи с авариями на атомных реакторах TMI-2 (США, 1979) и в Чернобыле (СССР, 1986) в последние два десятилетия повысились требования к безопасности АЭС, и вырос интерес к исследованию возможных аварийных процессов. Современный подход к анализу аварий заключается в создании и использовании специальных численных кодов. Это вызвано большими техническими и финансовыми трудностями при проведении полномасштабных экспериментов с воспроизведением условий аварии. Построение кодов основано на разработке физических моделей, которые сначала валидируются на мелкомасштабных экспериментах, а затем используются для описания реальных процессов. В связи с постоянным обновлением экспериментальной базы, а также с растущей необходимостью описывать процессы во все более широком диапазоне параметров возникает потребность в непрерывном совершенствовании кодов. Это делает крайне актуальным уточнение старых и разработку новых физически обоснованных моделей.

Одним из важнейших элементов аварии, требующих адекватного моделирования, является контакт нагретых конструкций или нагретых материалов с теплоносителем. Как правило, такие контакты сопровождаются бурным кипением. Протекание аварии в целом и тяжесть ее последствий в существенной степени зависят от характера кипения и реализующегося при этом уровня теплоотдачи. В зависимости от условий возможны различные режимы кипения: пузырьковое, пленочное или переходное. Со сменой режимов кипения, или иначе — с кризисами, связано резкое изменение теплоотдачи. Серьезную опасность представляет, например, кризис пузырькового кипения при импульсном энерговыделении в твэл в ходе реактивностной аварии [1]. Ухудшение теплоотдачи при наступлении кризиса может повлечь плавление оболочки и разрушение топливных стержней. Следует отметить, что с изменением внешних условий также могут количественно и качественно меняться характеристики каждого из режимов кипения [2]. Это сказывается в том числе и на уровне теплоотдачи.

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

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

Цель работы.

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

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

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

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

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

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

Диссертация состоит из введения, четырех глав, выводов и списка литературы. Каждая глава начинается кратким обзором по соответствующей проблеме и заканчивается выводами, наиболее важные из которых сформулированы в конце работы. Диссертация изложена на 139 страницах текста и содержит 2 таблицы, 55 рисунков, 117 наименований цитированной литературы.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ.

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

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

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

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