Теплоэнергетика, 2024, № 4, стр. 21-29

Неявный численный метод интегрирования уравнений сохранения двухжидкостной модели расчетного кода КОРСАР

Ю. В. Юдов a*, И. Г. Данилов a

a Научно-исследовательский технологический институт им. А.П. Александрова (НИТИ)
188540 г. Сосновый Бор, Копорское шоссе, д. 72, Россия

* E-mail: yudov@niti.ru

Поступила в редакцию 10.10.2023
После доработки 31.10.2023
Принята к публикации 01.11.2023

Аннотация

Представлен безытерационный неявный метод решения дискретных уравнений сохранения двухжидкостной модели расчетного кода (РК) КОРСАР/ГП (далее КОРСАР). Код КОРСАР/ГП является совместной разработкой сотрудников НИТИ и ОКБ “Гидропресс”. Он был аттестован в 2021 г. в Ростехнадзоре применительно к расчетному обоснованию безопасности реакторных установок с водо-водяным энергетическим реактором. В коде заложена полунеявная численная схема, которая ограничивает временной шаг интегрирования условием Куранта по скорости двухфазного потока. Для ускорения расчетов длительных переходных режимов реакторных установок разработан неявный численный метод без ограничения временного шага условием Куранта, основанный на методе SETS (stability-enhancing two-step). Он базируется на полунеявной схеме. Предварительно перед ее использованием на каждом шаге по времени решаются дискретные уравнения сохранения количества движения фаз для неявного выражения конвективных членов. После полунеявного этапа на новом временном слое вычисляются удельные (на единицу объема) массы и энергии фаз, которые являются донорными величинами в конвективных членах уравнений переноса. В отличие от SETS-метода, в неявном методе, разработанном для РК КОРСАР, используется полунеявная схема с линеаризацией нестационарных членов изменения удельных массы и энергии двухфазного потока, что позволяет проводить безытерационное решение дискретных уравнений. Однако реализация этой процедуры требует определения неизвестных скалярных переменных в расчетных ячейках: удельных энтальпий фаз, объемного паросодержания, давления – и согласования полей скоростей фаз с новым полем давления. Поэтому повторно применяется полунеявная схема с линеаризацией нестационарных членов в конце временного шага с перевычисленными донорными величинами. Работоспособность и эффективность разработанного неявного метода подтверждены путем решения по коду КОРСАР тестовой задачи с течением двухфазного потока под перепадом давления в обогреваемой горизонтальной трубе.

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

Список литературы

  1. Опыт создания и основные характеристики теплогидравлического расчетного кода нового поколения КОРСАР / В.А. Василенко, Ю.А. Мигров, С.Н. Волкова, Ю.В. Юдов, И.Г. Данилов, В.Г. Коротаев, В.В. Кутьин, Б.Р. Бондарчик, Д.В. Бенедиктов // Теплоэнергетика. 2002. № 11. С. 11–16.

  2. Юдов Ю.В. Двухжидкостная модель нестационарной контурной теплогидравлики и ее численная реализация в расчетном коде КОРСАР // Теплоэнергетика. 2002. № 11. С. 17–21.

  3. RELAP5/MOD2. Code manual, volume I: Code structure, system models and solution methods. NUREG/CR-4312-V1. Idaho, 1985.

  4. RELAP5/MOD3.3. Code manual, volume I: Code structure, system models and solution methods. NUREG/CR-5535/Rev1-V1. Idaho, 2001.

  5. TRAC-PF1/MOD2. Theory manual. NM 87545. Los Alamos, 1993.

  6. Hanna B. CATHENA: a thermalhydraulic code for CANDU analysis // Nucl. Eng. Des. 1998. V. 180. P. 113–131.

  7. TRACE V5.0. Theory manual: Field equations, solution methods and physical models. DC 20555-0001. Washington, 2012.

  8. Юдов Ю.В. Разработка двухжидкостной модели контурной теплогидравлики реакторных установок с водяным теплоносителем: дис. … канд. техн. наук. Сосновый Бор: НИТИ, 2001.

  9. Юдов Ю.В. Коррекция полунеявной численной схемы двухжидкостной модели кода КОРСАР // Теплоэнергетика. 2019. № 1. С. 75–84. https://doi.org/10.1134/S0040363619010090

Дополнительные материалы отсутствуют.