Моделирование конвекции с использованием одномерных кодов
Показана возможность физически корректного моделирования процессов тепломассопереноса в однофазных средах с помощью одномерных компьютерных кодов. Приведены рекомендации по выбору оптимальных нодализационных схем. Описано можливість фізично коректного моделювання процесів тепломасоперенесення в о...
Saved in:
| Published in: | Ядерна та радіаційна безпека |
|---|---|
| Date: | 2013 |
| Main Author: | |
| Format: | Article |
| Language: | Russian |
| Published: |
Державне підприємство "Державний науково-технічний центр з ядерної та радіаційної безпеки" Держатомрегулювання України та НАН України
2013
|
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/97484 |
| 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: | Моделирование конвекции с использованием одномерных кодов / А.А. Чуклин // Ядерна та радіаційна безпека. — 2013. — № 4. — С. 10-15. — Бібліогр.: 8 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859636645997838336 |
|---|---|
| author | Чуклин, А.А. |
| author_facet | Чуклин, А.А. |
| citation_txt | Моделирование конвекции с использованием одномерных кодов / А.А. Чуклин // Ядерна та радіаційна безпека. — 2013. — № 4. — С. 10-15. — Бібліогр.: 8 назв. — рос. |
| collection | DSpace DC |
| container_title | Ядерна та радіаційна безпека |
| description | Показана возможность физически корректного моделирования процессов
тепломассопереноса в однофазных средах с помощью одномерных компьютерных кодов.
Приведены рекомендации по выбору оптимальных нодализационных схем.
Описано можливість фізично коректного моделювання процесів
тепломасоперенесення в однофазних середовищах за допомогою одновимірних
комп’ютерних кодів. Наведено рекомендації з вибору оптимальних нодалізаційних схем.
This paper describes the possibility of physically correct simulation of heat and mass
transfer in single-phase flows with help of one-dimensional computer codes that are widely used
to develop computer models of nuclear power plants for analyzing their emergency regimes. For
solving the problem of convection modeling by means of computer codes, the results of the bench
experiment are compared with results of the experiment performed using the system of
differential equations and with results of the experiment performed using the MELCOR
computer code. Recommendations for the selection of optimal nodalization schemes for
developing computer models with convection are provided.
|
| first_indexed | 2025-12-07T13:16:58Z |
| format | Article |
| fulltext |
10 ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013
УДК 621.036.24.004.415.3
А. А. Чуклин
Севастопольская служба научно-технической поддержки,
ОП НТЦ НАЭК «Энергоатом», г. Севастополь, Украина
Моделирование конвекции
с использованием
одномерных кодов
Показана возможность физически корректного моделирования
процессов тепломассопереноса в однофазных средах с помощью од-
номерных компьютерных кодов. Приведены рекомендации по выбору
оптимальных нодализационных схем.
К л ю ч е в ы е с л о в а: конвекция, компьютерное моделирование,
одномерные коды, MELCOR.
О. О. Чуклін
Моделювання конвекції з використанням одновимірніx
кодів
Описано можливість фізично коректного моделювання процесів
тепломасоперенесення в однофазних середовищах за допомогою од-
новимірних комп’ютерних кодів. Наведено рекомендації з вибору опти-
мальних нодалізаційних схем.
К л ю ч о в і с л о в а: конвекція, комп’ютерне моделювання, однови-
мірні коди, MELCOR.
© А. А. Чуклин, 2013
Р
ассматриваемый в данной статье компьютерный
код MELCOR является одномерным с «сосредото-
ченными» параметрами. Параметры по всему мо-
делируемому контрольному объему эквивалентны
друг другу, что значительно усложняет модели-
рование трехмерных систем. Существенная проблема, воз-
никающая при моделировании неодномерных процессов
тепломассообмена и гидродинамики, — учет рециркуляции,
влияющей на их динамику.
Рециркуляция в объемах связана с возникновением
конвективного перемешивания. Конвекция — это явление
переноса теплоты в жидкостях, газах или сыпучих средах
потоками вещества за счет переноса массы. Естественная
конвекция возникает в веществе самопроизвольно при его
неравномерном нагревании в поле тяготения.
Для создания квазитрехмерных моделей приходится
прибегать к более детальной нодализации моделируемых
элементов. Одна из задач, стоящих перед разработчи-
ком компьютерной модели энергоустановки, заключается
в выборе оптимального количества элементов нодализа-
ционной схемы.
В случае упрощенного моделирования программа может
не учесть некоторые физические процессы, в частности:
конвекцию (разогрев недогретого теплоносистеля
в компенсаторе давления и т. п.);
перетоки сред между участками моделируемого объекта
(например, верхняя камера смешения реактора, заданная
одним контрольным объемом, не позволит учесть перето-
ки теплоносителя при отключении одного из ГЦН);
температурные изменения по длине и высоте объекта
(например, разбивка трубчатки парогенератора на участки
по длине и высоте).
С другой стороны, использование неоправданно дета-
лизированной нодализации приведет как к увеличению
расчетного времени, так и к накоплению вычислительной
погрешности [1].
Для рассмотрения задачи учета конвекции одномерны-
ми кодами сравним результаты стендового эксперимента
с результатами расчетного эксперимента, выполненного
с использованием системы дифференциальных уравнений,
а также с результатами расчетного эксперимента, выпол-
ненного с использованием компьютерного кода MELCOR.
Стендовый эксперимент. Экспериментальный стенд
представляет собой заполненную водой металлическую
емкость, к одной из стенок которой подводится тепло
от нагретой до 100 °С жидкости, которая, в свою очередь,
получает тепло от трубчатого электронагревателя мощ-
ностью 2000 Вт. В качестве измерительных приборов ис-
пользуются видеорегистратор, двухмерная измерительная
плоскость и термопары (ТП), расположенные согласно
схеме, изображенной на рис. 1. Для визуализации поля
скоростей в емкость вводится краситель — раствор фу-
корцина (лат. fucorcinum, также жидкость Кастеллани или
краска Кастеллани), который обладает характерным ярко-
малиновым цветом.
В процессе эксперимента вблизи греющей стенки на-
блюдается устойчивый поток жидкости, устремленный
вверх. Толщина слоя находится в пределах 0,019 м, ско-
рость подъема слоя — 0,05 м/с. Движение жидкости в верх-
ней части объема происходит в направлении «от» греющей
стенки. Скорость и толщина слоя жидкости, движущейся
вдоль поверхности, изменяются со временем в зависимости
от отдаленности потока от греющей стенки. Средняя ско-
рость потока на поверхности составляет 0,027 м/с. Ввиду
отсутствия «холодного» источника (на противоположной
от греющей стенки стороне), движение слоя жидкости
ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013 11
Моделирование конвекции с использованием одномерных кодов
вдоль стенки вниз (подобно слою жидкости, движущему-
ся вверх вблизи «горячей» стенки) почти не происходит.
Согласно показаниям термопар, в верхней части объема
выделяется зона с нагретой жидкостью, толщина которой
с течением времени увеличивается. Начинается медленное
закручивание всего нагретого потока (рис. 2).
Рис. 2. Графическое изображение результатов
стендового эксперимента
Поскольку «холодный» источник отсутствует, процесс
является нестационарным: в связи с постепенным нагре-
ванием жидкости изменяются поля температур и скоро-
стей. Для сравнения результатов стендового эксперимента
с результатами расчетов рассматривается интервал вре-
мени от начала введения «горячего» источника в систему
до повторного достижения «языком» нагретой жидкости
греющей стенки.
Расчетный эксперимент с использованием матема-
тической модели. Для расчетного эксперимента создана
математическая модель экспериментальной установки.
Наиболее популярной моделью для описания конвекции
в жидкостях и газах являются уравнения тепловой конвек-
ции в приближении Буссинеска — Обербека (А. Oberbeck,
1879; J. Boussinesq, 1903). Модель включает уравнение
Навье — Стокса, уравнение теплопроводности и уравне-
ние несжимаемости. Основная идея приближения состоит
в особенности учета зависимости плотности от темпера-
туры [2].
Примем, что характер течения жидкости в рассматрива-
емом объеме одинаков во всех вертикальных плоскостях,
перпендикулярных к греющей поверхности. В таком слу-
чае для исследования кинематики и динамики потока до-
статочно рассмотреть плоскую задачу, т. е. когда скорость
и давление жидкости зависят только от двух координат (х; y).
Уравнением движения, описывающим изменение ско-
рости во времени и пространстве, является уравнение
Навье — Стокса [2], которое в переменных «скорость —
давление» имеет следующий вид:
( )dw p
w w w
d
∇
+ ∇ = − + µ∇
τ ρ
, (1)
где w — скорость; τ — время; ρ — плотность; p — давление;
µ — динамический коэффициент вязкости; Δ — оператор
Лапласа; ∇ — оператор Набла.
Уравнение неразрывности для несжимаемой жидкости
в двухмерной задаче имеет вид
0yx ww
x y
∂∂
+ =
∂ ∂
, (2)
где wx — проекция скорости на координатную ось х; wу —
проекция скорости на координатную ось у; у — у-координата;
х — х-координата.
Чтобы сделать систему уравнений замкнутой, необхо-
димо добавить дифференциальное уравнение, описыва-
ющее температурное поле в движущейся жидкости, т. е.
уравнение энергии. При записи дифференциального урав-
нения энергии принимается, что жидкость однородна
и изотропна, энергия деформации мала по сравнению
с изменением внутренней энергии. Уравнение энергии
имеет вид [3]
2 2
2 2x y
t t t t t
w w a
x y x y
∂ ∂ ∂ ∂ ∂
+ + = − + ∂τ ∂ ∂ ∂ ∂
, (3)
где t — температура; а — коэффициент теплопроводности.
Для моделирования плоских течений широко исполь-
зуются переменные «функция тока — вихрь скорости»
вместо физических (естественных) переменных «ско-
рость — давление» [4, 5]. Проекции скоростей на коорди-
натные оси х и у, представленные через функцию тока ψ,
имеют вид
; x yw w
y x
∂ψ ∂ψ
= = −
∂ ∂
, (4)
где ψ — функция тока.
Вихрь скорости двухмерных течений представлен в виде
yx ww
x y
∂∂
ω = −
∂ ∂
, (5)
где ω — вихрь скорости.
Приняв во внимание (4), получим уравнение [5, 7]
.∆ψ = −ω (6)
Рис. 1. Схема
экспери ментального стенда
Двухмерная
измерительная
плоскость
Емкоскть с водой
Источник тепла,
заполненный водой
Нагревательный
элемент
—Место
расположения ТП
12 ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013
А. А. Чуклин
Для двухмерного случая уравнение Навье — Стокса (1)
в переменных «функция тока — вихрь скорости» имеет
вид [4]
( ) .
w
w
∂
+ ⋅ ∇ ω = µ∆ω
∂τ
(7)
Система уравнений (3), (6) и (7) является искомой си-
стемой дифференциальных уравнений, используемой
для задачи расчета конвекции в объеме. Полученная си-
стема уравнений решалась при помощи программного па-
кета Flexpde [6]. В качестве граничных условий для рас-
сматриваемой системы уравнений использовались данные,
соответствующие параметрам экспериментального стенда.
Результаты расчета представлены на рис. 3 и 4.
Результаты расчета показывают, что вблизи грею-
щей стенки наблюдается устойчивый поток жидкости,
устремленный вверх. Толщина слоя в среднем составляет
0,02 м, средняя скорость подъема слоя — 0,06 м/с (рис. 3).
Движение жидкости в верхней части объема происходит
в направлении «от» греющей стенки. Скорость и толщина
слоя жидкости, движущейся вдоль поверхности, изменя-
ется со временем в зависимости от отдаленности потока
от греющей стенки. Средняя скорость потока на поверх-
ности составляет 0,03 м/с. Подобно результатам стендового
эксперимента, начинается медленное закручивание всего
нагретого потока (рис. 3). Рассогласования в полученных
значениях обусловлены как некоторыми неточностями за-
данных граничных условий расчета, так и погрешностью
измерительных приборов стендовой установки.
Расчетный эксперимент с использованием кода MELCOR.
Модель экспериментального стенда, созданная с помощью
кода MELCOR, состоит из шести контрольных объемов,
шести путей соединения и трех тепловых структур (рис. 5),
исходные параметры которых соответствуют реальным па-
раметрам рассмотренного выше стенда.
Контрольные объемы Cv001 и Cv002 представляют со-
бой вертикальные объемы, которые моделируют «подъ-
емный» и «опускной» участки нагреваемой жидкости,
соответственно. Контрольные объемы Cv011 и Cv012 —
горизонтальные объемы, моделирующие «верхний»
и «нижний» участки нагреваемой жидкости. Размеры
Рис. 3. Результаты расчетного эксперимента:
распределение скорости потока в объеме
Рис. 4. Результаты расчетного эксперимента:
распределение температуры в объеме
ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013 13
Моделирование конвекции с использованием одномерных кодов
контрольных объемов принимались исходя из результатов,
полученных расчетными и экспериментальными путями.
Наличие контрольного объема Cv012 позволяет учесть
возможное присутствие «холодного» источника, от воздей-
ствия которого жидкость будет двигаться вниз. Источник
тепла модулируется контрольным объемом Cv100, взаимо-
действующим с конвективным объемом через теплопере-
дающую поверхность.
Теплопередающая стенка моделируется с использо-
ванием пакета Heat Structures [1]. Тепловые структуры
(HS00101, HS00301 и HS00401) с «левой» стороны взаимо-
действуют с контрольным объемом Cv100, а с «правой» —
с контрольными объемами Cv001, Cv011 и Cv012 (в зави-
симости от высотной отметки).
Контрольный объем Cv500 модулирует окружающую
среду.
Результаты расчета показали, что скорость движения
потока вблизи греющей стенки (Cv001) равна 0,015 м/с.
Полученная величина меньше экспериментального
значения и значения, рассчитанного с использованием
математической модели (0,05…0,06 м/с). Заниженное
значение скорости движения потока жидкости вблизи
греющей стенки снижает теплоотдачу от стенки к жидкости,
что негативно влияет на корректность результатов расчета.
Было принято решение произвести более детальную
нодализацию рассматриваемой установки (рис. 6).
В результате применения более детальной разбивки
рассматриваемого объема скорость вблизи греющей стенки
(контрольный объем Cv001) увеличилась до 0,051 м/с
(рис. 7).
Так как рассогласования параметров контрольных
объемов Cv002—Cv010 относительно невелики, можно
принять допущение об их эквивалентности и, без влияния
на результат, снизить количество расчётных объемов и пу-
тей соединения, объединив объемы Cv002—Cv010 в один
(рис. 8).
Для выяснения влияния соотношения сечений «подъ-
емного» и «опускного» участков модели на результаты рас-
чета были проведены эксперименты, результаты которых
Рис. 5. Нодализационная
схема установки
Рис. 6. Нодализационная схема установки с детальной
разбивкой «подъемного» и «опускного» участков
Рис. 7. Скорость движения жидкости
в контрольных объемах Сv001-Cv010
Рис. 8. Нодализационная схема установки
с объединенным опускным участком
14 ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013
А. А. Чуклин
показали, что оптимальная толщина подъемного участка
находится в районе 0,02 м, что соответствует толщине подъ-
емного участка экспериментального стенда. Аналогичным
образом получена оптимальная высота верхнего контроль-
ного объема, которая с некоторой зоной нечувствительно-
сти соответствует толщине верхнего слоя эксперименталь-
ной установки. Результаты расчета представлены на рис. 9.
Рис. 9. Скорость движения жидкости
в объемах Cv001 и Cv002
В результате произведенных изменений скорость
на подъемном участке осталась равной 0,051 м/с. Скорость
в опускном участке соответствует средней скорости не-
объединенных опускных объемов (Cv002—Cv010) и равна
0,006 м/c. Скорость движения жидкости на верхнем уча-
стке равна 0,04 м/c.
Параметры установки, полученные эксперименталь-
ным путем, путем расчёта системы дифференциальных
уравнений и при моделировании компьютерным кодом
MELCOR, сведены в табл. 1.
Таблица 1. Сравнительная таблица
экспериментальных и расчетных результатов
Параметр
Стендо-
вый экс-
перимент
Матема-
тическая
модель
Компью-
терная
модель
Скорость движения
жидкос ти подъемного
участка, м/с
0,05 0,06 0,051
Толщина подъемного
участ ка, м 0,019 0,02 0,02
Скорость движения
жидкос ти вдоль
поверхности, м/с
0,027 0,03 0,038
Скорость движения
жидкости опускного
участка, м/с
0,005 0,009 0,006
Данные, приведенные в табл. 1, физически корректны
и имеют достаточную сходимость. Рассогласования в зна-
чениях обусловлены:
неточностью заданных граничных условий расчетов;
погрешностью измерительных приборов стендовой
установки;
использованием различных подходов для реализации
расчета.
Основываясь на полученных результатах, можно ут-
верждать, что для физически корректного моделирования
конвективных систем с использованием одномерных кодов
необходимо, наряду с основными параметрами системы,
обладать информацией о толщине подъемного и верхнего
участков. Такую информацию можно получить исходя
из температуры греющей поверхности, температуры нагре-
ваемой среды, количества теряемой теплоты и т. п.
В качестве примера реализации предложенной мето-
дики по моделированию конвективного перемешивания
с использованием одномерных компьютерных кодов рас-
смотрим моделирование ламинарного участка бассейна
с граничными условиями, эквивалентными рассматривае-
мому экспериментальному стенду.
Для определения параметров ламинарного подъемного
участка воспользуемся формулой [8]
( )0 250 5 0 253 93 0 952 ,, , , Pr , Pr Grx − −δ = ⋅ + , (8)
где δ — ширина пограничного слоя; x — х-координата
рассматриваемого элемента; Pr — число Прандтля; Gr —
число Грасгофа.
Локальные числа Грасгофа и Прандтля определяются
по выражениям [8]
( )3
0
2Gr wgx T Tβ −
=
ν
, (9)
где β — коэффициент объемного расшения; g — ускорение
свободного падения; ν — кинематический коэффициент
вязкости; Тw — температура «горячей» стенки; Т0 — тем-
пература жидкости вдали от греющей стенки;
Pr ,
a
ν
= (10)
где a — коэффициент теплопередачи.
Определив ширину пограничного слоя, можно найти
скорость свободной конвекции, т. е. скорость жидкости,
движущейся вблизи греющей поверхности вверх:
( ) 2 3 40
0 2
1 1 1
4 2 3 12
x w
g
w T T y y y y
ρ β δ = − − + − µ δ δ
, (11)
где ρ0 — начальная плотность жидкости; µ — динамиче-
ский коэффициент вязкости; у — у-координата рассматри-
ваемого элемента.
Начальным условиям: температура стенки 90 °С, тем-
пература воды 20 °С, высота подъемного участка 0,15 м —
соответствуют (рис. 10) ширина пограничного слоя
0,021 м (8) и скорость конвекции 0,052 м/c (11). При мо-
делировании рассматриваемой системы с использованием
одномерных кодов для получения физически корректных
результатов необходимо производить разбиение на нода-
лизационные элементы с учетом рассчитанных значений
толщины и скорости подъемного слоя. Аналогично моде-
лируется конвекция с турбулентным подъемным участком.
Если получить расчётные значения параметров подъемно-
го участка невозможно, следует проводить стендовый экс-
перимент, результаты которого, согласно предложенной
методике, должны быть использованы при моделировании
рассматриваемой системы.
ISSN 2073-6237. Ядерна та радіаційна безпека 4(60).2013 15
Моделирование конвекции с использованием одномерных кодов
Выводы
В используемой модели стенда для получения кор-
ректных результатов моделирования конвекции в объеме
вполне достаточно использования четырех контрольных
объемов, моделирующих «подъемный», «опускной», «верх-
ний» и «нижний» участки конвективной емкости. При этом
параметры нодализационных элементов (толщина участка
и скорость движения среды) должны соответствовать рас-
четным или экспериментальным данным. Применение
большего количества элементов в рассматриваемой схеме
не приводит к получению более корректных результатов.
В представленных расчетах количество элементов
в смоделированной установке практически не влияет
на расчетное время. Однако в случае построения слож-
ной модели, такой как модель энергоустановки, следует
стремиться к уменьшению количества объемов, если это
не влияет на корректность результатов. Увеличение ко-
личества элементов, используемых для моделирования
сложных систем, ведет как к значительному увеличению
расчетного времени, так и к накоплению вычислительной
погрешности.
Согласно полученным результатам, моделирование про-
странственных систем с перетоками сред внутри объема
(например, конвективное перемешивание) с использова-
нием одномерных кодов вполне допустимо и при правиль-
ной нодализации дает физически корректные результаты.
Список использованной литературы
1. MELCOR Computer Code Manuals. Reference Manuals. Ver-
sion 1.8.5 October 2000, NUREG/CR-6119, Vol. 1, Rev. 2 SAND2000-
2417/1, Sandia National Laboratories, Albuquerque, NM 87185–0739.
2. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. Теоретическая
физика. — Т. VI. — 3-е изд., перераб. — М.: Наука, 1986. —736 с.
3. http://www.fast-const.ru/articles.php?article_id = 21
4. Никонов В. В. Развитие вихревых методов расчета обте-
кания тел несжимаемыми невязким и вязким потоками: дис. …
канд. техн. наук: 01.02.05 — механика жидкости, газа и плазмы /
Никонов Валерий Владимирович. — Самара, 2007. — 174 с.
5. Самарский А. А. Вычислительная теплопередача / А. А. Са-
марский, П. Н. Вабищевич. — М.: Едиториал УРСС, 2003. — 784 с.
6. FlexPDE_User Guide_Version 2.19_10/30/00_Copyright ©2000
PDE Solutions Inc
7. Григорьев Ю. Н. Метод вихрей в ячейках для плоских те-
чений сжимаемых, неоднородных и релаксирующих газов /
Ю. Н. Григорьев; Ин-т вычислительных технологий (Новоси-
бирск) // Всерос. конф. по вычислительной математике КВМ-
2007. 18–20 июня 2007 г. : Тезисы докладов : Численное решение
диф. и интеграл. уравнений. — Новосибирск: Ин-т вычислитель-
ной математики и математической геофизики СОРАН. — 2007.
8. Дульнев Г. Н. Теория тепло- и массообмена / Г. Н. Дуль-
нев. — СПб: НИУ ИТМО, 2012. — 195 с.
References
1. MELCOR Computer Code Manuals. Reference Manuals. Ver-
sion 1.8.5 October 2000, NUREG/CR-6119, Vol. 1, Rev. 2 SAND2000-
2417/1, Sandia National Laboratories, Albuquerque, NM 87185–0739.
2. L. Landau, E. Lifshitz. Hydrodynamics. Theoretical physics. —
V. VI — 3rd publication. — M.: Science, 1986. — 736 p.) (Rus)
3. http://www.fast-const.ru/articles.php?article_id = 214.
4. Nikonov V.. Development of methods for calculating the vortex
flow of bodies and non-viscous incompressible viscous flows: Dis. Sci-
ences: 01.02.05 — Mechanics of liquid, gas and plasma / Valery Niko-
nov. — Samara, 2007. — 174 p. (Rus)
5. Samarsky A. Computational Heat Transfer / A. Samarsky,
P. Vabishchevich. — Moscow: Editorial URSS, 2003. — 784. (Rus)
6. FlexPDE_User Guide _Version 2.19_10/30/00_Copyright ©
2000 PDE Solutions Inc.
7. Grigoryev Yu. Method of vortices in the cells for plane flows
of compressible, non-uniform and relaxing gases / Yu. Grigoryev, In-
stitute of Computing Technology (Novosibirsk) / Conf. on Computa-
tional Mathematics ICCM-2007. 18–20 June 2007: Abstracts: The nu-
merical solution of differential. and integral. equations. —Novosibirsk:
Institute of Computational Mathematics and Mathematical Geophysics
SB RAS. — 2007. (Rus)
8. Dulnev G. Theory of Heat and Mass Transfer / G. Dulnev. —
St. Petersburg: ITMO, 2012. — 195 p. (Rus)
Получено 16.07.2013.
Рис. 10. Распределение скорости подъемного участка
в зависимости от отдаления от греющей стенки
|
| id | nasplib_isofts_kiev_ua-123456789-97484 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 2073-6231 |
| language | Russian |
| last_indexed | 2025-12-07T13:16:58Z |
| publishDate | 2013 |
| publisher | Державне підприємство "Державний науково-технічний центр з ядерної та радіаційної безпеки" Держатомрегулювання України та НАН України |
| record_format | dspace |
| spelling | Чуклин, А.А. 2016-03-28T18:33:26Z 2016-03-28T18:33:26Z 2013 Моделирование конвекции с использованием одномерных кодов / А.А. Чуклин // Ядерна та радіаційна безпека. — 2013. — № 4. — С. 10-15. — Бібліогр.: 8 назв. — рос. 2073-6231 https://nasplib.isofts.kiev.ua/handle/123456789/97484 621.036.24.004.415.3 Показана возможность физически корректного моделирования процессов тепломассопереноса в однофазных средах с помощью одномерных компьютерных кодов. Приведены рекомендации по выбору оптимальных нодализационных схем. Описано можливість фізично коректного моделювання процесів тепломасоперенесення в однофазних середовищах за допомогою одновимірних комп’ютерних кодів. Наведено рекомендації з вибору оптимальних нодалізаційних схем. This paper describes the possibility of physically correct simulation of heat and mass transfer in single-phase flows with help of one-dimensional computer codes that are widely used to develop computer models of nuclear power plants for analyzing their emergency regimes. For solving the problem of convection modeling by means of computer codes, the results of the bench experiment are compared with results of the experiment performed using the system of differential equations and with results of the experiment performed using the MELCOR computer code. Recommendations for the selection of optimal nodalization schemes for developing computer models with convection are provided. ru Державне підприємство "Державний науково-технічний центр з ядерної та радіаційної безпеки" Держатомрегулювання України та НАН України Ядерна та радіаційна безпека Моделирование конвекции с использованием одномерных кодов Моделювання конвекції з використанням одномірніx кодів Simulation of Convection Using One-Dimensional Codes Article published earlier |
| spellingShingle | Моделирование конвекции с использованием одномерных кодов Чуклин, А.А. |
| title | Моделирование конвекции с использованием одномерных кодов |
| title_alt | Моделювання конвекції з використанням одномірніx кодів Simulation of Convection Using One-Dimensional Codes |
| title_full | Моделирование конвекции с использованием одномерных кодов |
| title_fullStr | Моделирование конвекции с использованием одномерных кодов |
| title_full_unstemmed | Моделирование конвекции с использованием одномерных кодов |
| title_short | Моделирование конвекции с использованием одномерных кодов |
| title_sort | моделирование конвекции с использованием одномерных кодов |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/97484 |
| work_keys_str_mv | AT čuklinaa modelirovaniekonvekciisispolʹzovaniemodnomernyhkodov AT čuklinaa modelûvannâkonvekcíízvikoristannâmodnomírníxkodív AT čuklinaa simulationofconvectionusingonedimensionalcodes |