+78136960367 Дискретизация уравнений по пространственной координате осуществлена по методу контрольных объемов с аппроксимацией конвективных членов по схеме против потока. При этом по пространственной координате используется смещенная (шахматная) сетка, когда скалярные переменные потока (давление, энтальпии и объемные доли фаз) определяются относительно центров расчетных ячеек, а векторные переменные (скорости фаз) рассчитываются на границах последних.
Для численного решения систем конечно-разностных уравнений применена разработанная безытерационная полунеявная схема с автоматическим выбором шага интегрирования по времени. Более подробно используемая в коде КОРСАР численная схема представлена в [2, 10].
Для численного решения нестационарных уравнений теплопроводности в теплопроводящих конструкциях используется метод прогонки.
Для численного решения системы уравнений точечной модели кинетики используется неявный метод Эйлера. Для расчета динамики отдельных элементов оборудования интегрирование систем обыкновенных дифференциальных уравнений 1-го порядка осуществляется для каждого элемента автономно методами Эйлера или Рунге-Кутта-Мерсона с подключением к блоку расчета контурной теплогидравлики по явной схеме.