All issues
- 2026 Vol. 18
- 2025 Vol. 17
- 2024 Vol. 16
- 2023 Vol. 15
- 2022 Vol. 14
- 2021 Vol. 13
- 2020 Vol. 12
- 2019 Vol. 11
- 2018 Vol. 10
- 2017 Vol. 9
- 2016 Vol. 8
- 2015 Vol. 7
- 2014 Vol. 6
- 2013 Vol. 5
- 2012 Vol. 4
- 2011 Vol. 3
- 2010 Vol. 2
- 2009 Vol. 1
-
Исследование достижимости цели в медицинском квесте
Компьютерные исследования и моделирование, 2025, т. 17, № 6, с. 1149-1179В работе представлено экспериментальное исследование древовидной структуры, возникающей при медицинском обследовании. При каждой встрече с медицинским специалистом пациент получает некоторое количество направлений на консультации других специалистов или на анализы. Возникает дерево направлений, каждую ветвь которого должен пройти пациент. В зависимости от разветвленности дерева оно может быть как конечным (и в этом случае обследование может быть завершено), так и бесконечным, когда цель пациента не может быть достигнута. В работе как экспериментально, так и теоретически изучаются критические свойства перехода системы из леса конечных деревьев в лес бесконечных в зависимости от вероятностных характеристик дерева.
Для описания предлагается модель, в которой дискретная функция вероятности числа ветвей на узле повторяет динамику непрерывного гауссового распределения. Характеристики распределения Гаусса (математическое ожидание $x_0$, среднеквадратичное отклонение $\sigma$) являются параметрами модели. В выбранной постановке задача относится к проблематике ветвящихся случайных процессов (ВСП) в неоднородной модели Гальтона – Ватсона.
Экспериментальное изучение проводится путем численного моделирования на конечных решетках. Построена фазовая диаграмма, определены границы областей различных фаз. Проведено сравнение с фазовой диаграммой, полученной из теоретических критериев для макросистем, установлено адекватное соответствие. Показано, что на конечных решетках переход является размытым.
Описание размытого фазового перехода проведено с помощью двух подходов. В первом (стандартном) подходе переход описывается с помощью так называемой функции включения, имеющей смысл доли одной из фаз в общем множестве. Установлено, что такой подход в данной системе неэффективен, поскольку найденное положение условной границы размытого перехода определяется только размером выбранной экспериментальной решетки и не несет объективного смысла.
Предлагается второй (оригинальный) подход, основанный на введении в рассмотрение параметра порядка, равного обратной средней высоте дерева, и анализа его поведения. Установлено, что динамика такого параметра порядка в сечениях $\sigma = \text{const}$ с очень небольшими отличиями имеет вид распределения Ферми – Дирака ($\sigma$ выполняет ту же функцию, что и температура для распределения Ферми – Дирака, $x_0$ — функцию энергии). Для параметра порядка подобрано эмпирическое выражение, введен и рассчитан аналог химического потенциала, который и имеет смысл характерного масштаба параметра порядка, то есть тех значений $x_0$, при которых условно можно считать, что порядок сменяется беспорядком. Этот критерий положен в основу определе- ния границы условного перехода в данном подходе. Установлено, что эта граница соответствует средней высоте дерева, равной двум поколениям. На основании обнаруженных свойств предложены рекомендации для медицинских учреждений, позволяющие контролировать обеспечение конечности траектории пациентов.
Рассмотренная модель и метод ее описания с помощью условно-бесконечных деревьев имеют приложение ко многим иерархическим системам. К таким системам можно отнести сети маршрутизации интернет-соединений, бюрократические сети, торговые, логистические сети, сети цитирования, игровые стратегии, задачи популяционной динамики и пр.
Ключевые слова: медицинское обследование, ветвящийся случайный процесс, модель Гальтона – Ватсона, размытые фазовые переходы, конечные системы, условно-бесконечные траектории, макросистема, функция включения, области почти чистых фаз, параметр порядка, химический потенциал, фазовая диаграмма, критическое поведение.
Research on the achievability of a goal in a medical quest
Computer Research and Modeling, 2025, v. 17, no. 6, pp. 1149-1179The work presents an experimental study of the tree structure that occurs during a medical examination. At each meeting with a medical specialist, the patient receives a certain number of areas for consulting other specialists or for tests. A tree of directions arises, each branch of which the patient should pass. Depending on the branching of the tree, it can be as final — and in this case the examination can be completed — and endless when the patient’s goal cannot be achieved. In the work both experimentally and theoretically studied the critical properties of the transition of the system from the forest of the final trees to the forest endless, depending on the probabilistic characteristics of the tree.
For the description, a model is proposed in which a discrete function of the probability of the number of branches on the node repeats the dynamics of a continuous gaussian distribution. The characteristics of the distribution of the Gauss (mathematical expectation of $x_0$, the average quadratic deviation of $\sigma$) are model parameters. In the selected setting, the task refers to the problems of branching random processes (BRP) in the heterogeneous model of Galton – Watson.
Experimental study is carried out by numerical modeling on the final grilles. A phase diagram was built, the boundaries of areas of various phases are determined. A comparison was made with the phase diagram obtained from theoretical criteria for macrosystems, and an adequate correspondence was established. It is shown that on the final grilles the transition is blurry.
The description of the blurry phase transition was carried out using two approaches. In the first, standard approach, the transition is described using the so-called inclusion function, which makes the meaning of the share of one of the phases in the general set. It was established that such an approach in this system is ineffective, since the found position of the conditional boundary of the blurred transition is determined only by the size of the chosen experimental lattice and does not bear objective meaning.
The second, original approach is proposed, based on the introduction of an parameter of order equal to the reverse average tree height, and the analysis of its behavior. It was established that the dynamics of such an order parameter in the $\sigma = \text{const}$ section with very small differences has the type of distribution of Fermi – Dirac ($\sigma$ performs the same function as the temperature for the distribution of Fermi – Dirac, $x_0$ — energy function). An empirical expression has been selected for the order parameter, an analogue of the chemical potential is introduced and calculated, which makes sense of the characteristic scale of the order parameter — that is, the values of $x_0$, in which the order can be considered a disorder. This criterion is the basis for determining the boundary of the conditional transition in this approach. It was established that this boundary corresponds to the average height of a tree equal to two generations. Based on the found properties, recommendations for medical institutions are proposed to control the provision of limb of the path of patients.
The model discussed and its description using conditionally-infinite trees have applications to many hierarchical systems. These systems include: internet routing networks, bureaucratic networks, trade and logistics networks, citation networks, game strategies, population dynamics problems, and others.
-
Численное моделирование возникновения пика напряжения при отражении ударно-волнового импульса от зернистой пористой среды
Компьютерные исследования и моделирование, 2026, т. 18, № 2, с. 359-375Исследование упругих волн в пористых средах актуально для задач поиска полезных ископаемых, при использовании пористых экранов для демпфирования ударно-волновых воздействий, для изучения строения земной коры и т. д. Упругие свойства пористой среды, о которых можно судить по скорости распространения разного рода волн, зависят от степени консолидации пористой среды. Например, насыпные среды (песок, стеклянные шарики, гранулированные материалы) обладают низкой скоростью звука (порядка 100 м/с), уплотнение таких сред сопровождается некоторым увеличением скорости, а их консолидация (песчаник, цементация газогидратом) приводит к многократному возрас- танию скорости акустических волн, порядка 2000–3000 м/с. В данной работе теоретически исследуется динамика волнового импульса в ударной трубе, содержащей слой насыпной среды. Численное моделирование проведено для условий эксперимента. Приводится описание экспериментальной установки типа «ударная труба». Установка состоит из камеры высокого давления (КВД), камеры низкого давления (КНД) и секции насыпной среды. Ударно-волновой импульс (УВИ) создается из-за разрыва диафрагмы между КВД и КНД. Динамика УВИ регистрируется пьезоэлектрическими датчиками, расположенными заподлицо с внутренней стороны трубы. В ударной трубе, оснащенной секцией насыпных сред, волна испытывает многократные отражения от поверхности изучаемой пористой среды и верхнего торца трубы. Переотраженные сигналы используются в качестве зондирующих импульсов для изучения изменений в пористой среде, вызванных повторными прохождениями ударно-волнового импульса, с периодом около 10 мс. Используется математическая модель, включающая уравнения сохранения массы, уравнения сохранения импульсов и энергии для газовой фазы и твердых частиц с замыкающими соотношениями. Описание процесса проводится для одномерного плоского движения газовой и дисперсной фаз. Для численного решения используется аппроксимация уравнений, основанная на методе контрольного объема. Численные результаты показали, что предложенная модель качественно и количественно правильно описывает появление резкого кратковременного увеличения полного напряжения (пика) при повторном прохождении импульса через слой насыпной среды, наблюдаемого в экспериментах.
Ключевые слова: насыпная среда, ударная труба, волновой импульс, математическая модель, многофазная среда, метод контрольного объема, численные расчеты, образование пика.
Numerical modeling of the occurrence of a stress peak during the reflection of a shock wave pulse from a granular porous medium
Computer Research and Modeling, 2026, v. 18, no. 2, pp. 359-375The study of elastic waves in porous media is relevant for mineral exploration, the use of porous screens for shock wave damping, and the study of the structure of the earth’s crust. The elastic properties of a porous medium, which can be judged by the propagation velocity of various types of waves, depend on the degree of consolidation of the porous medium. For example, bulk media (sand, glass beads, granular materials) have a low sound velocity (about 100 m/s); compaction of such media is accompanied by a slight increase in velocity, while their consolidation (sandstone, gas hydrate cementation) leads to a multiple increase in the acoustic wave velocity, on the order of 2000–3000 m/s. This paper theoretically investigates the dynamics of a wave pulse in a shock tube containing a layer of a bulk medium. Numerical modeling was performed under experimental conditions. A description of a shock tube experimental setup is provided. The setup consists of a high-pressure volume (HPV), a low-pressure volume (LPV), and a bulk medium section. A shock wave pulse (SWP) is generated by the rupture of a diaphragm between the HPV and LPV. The SWP dynamics are recorded by piezoelectric sensors located flush on the inside of the tube. In the shock tube, equipped with a bulk medium section, the wave experiences multiple reflections from the surface of the porous medium under study and the upper end of the tube. The reflected signals are used as probe pulses to study changes in the porous medium caused by repeated passages of the shock wave pulse, with a period of approximately 10 ms. A mathematical model is used that includes the equations of conservation of mass, momentum, and energy for the gas phase and solid particles with closure relations. The process is described for one-dimensional planar motion of the gas and dispersed phases. The numerical solution utilizes an approximation of the equations based on the control volume method. Numerical results have shown that the proposed model accurately describes, qualitatively and quantitatively, the occurrence of a sharp, short-term increase in the total voltage (peak) during repeated pulse passage through a layer of bulk material, as observed in experiments.
-
Классификатор на базе оптимизированной архитектуры YOLO11 с применением инструментов CBAM и Grad-CAM для надежной диагностики больных растений
Компьютерные исследования и моделирование, 2026, т. 18, № 4, с. 871-889Надежная диагностика болезней растений компьютерным зрением требует не только высокой точности классификации, но и стабильной обобщающей способности и интерпретируемого принятия решений. При этом механизмы внимания безусловно полезны для моделей глубокого обучения, но эффективность решения задачи зависит от места и способа интеграции этих механизмов в архитектуру нейронной сети. В данном исследовании представлен классификатор на базе оптимизированной архитектуры YOLOv11m, а также результаты анализа результатов диагностики болезней растений в зависимости от сверточных блоков CBAM (Convolutional Block Attention Module), внедренных на разных уровнях архитектуры. Выполнен сравнительный анализ трех архитектурных вариантов модели: базовой архитектуры YOLO11m-Cls без механизма внимания, варианта с добавлением CBAM только в магистральную часть backbone и гибридной архитектуры, сочетающей сокращенное внимание в магистральной части backbone с дополнительным блоком CBAM на уровне классификационного слоя head. Все модели были обучены и протестированы в одинаковых экспериментальных условиях с использованием обширного набора данных, содержащего около 90 тысяч изображений 38 типов болезней растений. Экспериментальные результаты ясно показывают, что равномерное внедрение CBAM в магистральную часть backbone снижает устойчивость модели и ухудшает ее способность к обобщению, что проявляется в повышении валидационной функции потерь и понижении точности (Top-1 ≈ 90,5%) по сравнению с базовой архитектурой. Напротив, гибридная архитектура обеспечивает баланс между устойчивостью и способностью детектировать с повышением параметров точности Top-1 до 99,71% и Top-5 до 99,99% на этапе валидации, близким по значениям к параметрам базовой модели. Анализ интерпретируемости на основе метода Grad-CAM также положительно характеризует гибридную модель, которая формирует улучшенные и биологически обоснованные карты активации, четко фокусирующиеся именно на признаках заболевания и практически исключающие отвлечение внимания на фоновые элементы изображения. Важно отметить, что целью данной работы является не достижение наилучшей производительности среди всех классификаторов, а предоставление доказательств на уровне архитектурного проектирования относительно того, каким образом размещение механизма внимания влияет на устойчивость и интерпретируемость в классификационных моделях на базе YOLO. Полученные результаты предоставляют практические рекомендации по проектированию надежных и интерпретируемых систем искусственного интеллекта для агропромышленного сектора.
Ключевые слова: диагностика болезней растений, компьютерное зрение, сверточная нейронная сеть, YOLO11, Convolutional Block Attention Module (CBAM), Explainable Artificial Intelligence (XAI), Grad-CAM.
Design-optimized YOLO11 classification via strategic CBAM attention injection and Grad-CAM explainability for reliable plant disease diagnosis
Computer Research and Modeling, 2026, v. 18, no. 4, pp. 871-889Reliable plant disease diagnosis requires not only high classification accuracy, but also stable generalization and interpretable decision-making. Although attention mechanisms are useful to these deep learning models, the performance also depends on where and how they are integrated into the network architecture. This study presents a design-optimized YOLO11m-based classification framework that systematically investigates the impact of Convolutional Block Attention Module (CBAM) injection at different architectural levels for plant disease diagnosis. We perform a comparative modelcontrolled analysis of three model architectures: (i) the baseline YOLO11m-Cls architecture lacking attention, (ii) only adding the backbone block along with CBAM and (iii) a hybrid architecture that includes reduced backbone attention combined with CBAM added at classification head. All models are trained and tested under the same experimental settings using a largescale dataset with around 90 000 images of 38 types of plant diseases. Experimental results clearly show that the uniform injection of CBAM into the backbone reduces stability but causes generalization to worsen with a higher validation loss and significantly lower Top-1 accuracy (≈ 90.5%), while the hybrid attention design balances stability and discrimination, with Top-1 accuracy up to 99.71%, Top-5 accuracy up to 99.99% and near-baseline validation behavior respectively Grad-CAMbased interpretability analysis also demonstrates that the hybrid model generates enhanced and biologically interpretable activation maps, which are more disease-specific with less distraction from background. Notably, the aim of our work is not for achieving superior performance over all classifiers but instead only to provide design-level evidence on how placing attention modulates rigidity and interpretability in YOLO-based classification models. The results provide practical architectural considerations for building dependable and interpretable AI systems in agriculture.
-
Исследование образования комплекса флаводоксина и фотосистемы 1 методами прямого многочастичного компьютерного моделирования
Компьютерные исследования и моделирование, 2009, т. 1, № 1, с. 85-91С помощью компьютерной модели, основанной на методах многочастичного прямого моделирования и броуновской динамики, изучается кинетика образования комплекса между компонентами фотосинтетической электронтранспортной цепи — белком флаводоксином и мембранным комплексом фотосистемы 1. Моделируется броуновское движение нескольких сотен молекул флаводоксина, учитываются электростатические взаимодействия и сложная форма молекул. С помощью данной модели удалось воспроизвести экспериментальную немонотонную зависимость константы связывания флаводоксина с фотосистемой 1. Это говорит о том, что для описания такого вида зависимости достаточно учета только электростатических взаимодействий.
Investigation of complex formation of flavodoxin and photosystem 1 by means of direct multiparticle computer simulation
Computer Research and Modeling, 2009, v. 1, no. 1, pp. 85-91Views (last year): 4. Citations: 2 (RSCI).Kinetics of complex formation between components of the photosynthetic electron transport chain — flavodoxin and membrane complex photosystem I has been studied using computer model based on methods of multiparticle simulation and Brownian dynamics. We simulated Brownian motion of several hundreds of flavodoxin molecules, taking into account electrostatic interactions and complex shape of the molecules. Our model could describe experimental nonmonotonic dependence of the association rate constant for flavodoxin and photosystem I. This lets us conclude that electrostatic interactions are sufficient to form such kind of nonmonotonic dependence.
-
Стехиометрия метаболических путей в динамике клеточных популяций
Компьютерные исследования и моделирование, 2011, т. 3, № 4, с. 455-475Проанализированы проблемы соответствия кинетических моделей клеточного метаболизма описываемому ими объекту. Изложены основы стехиометрии полного метаболизма и его больших частей. Описана биоэнергетическая форма стехиометрии, основанная на универсальной единице восстановленности химических соединений (редоксон). Выведены уравнения материально-энергетического баланса (биоэнергетической стехиометрии) метаболических потоков, в том числе баланса протонов с высоким электрохимическим потенциалом μH+ и макроэргических соединений. Получены соотношения, выражающие выход биомассы, скорость потребления источника энергии для роста и другие физиологически важные величины через биохимические характеристики клеточной энергетики. Вычислены значения максимального энергетического выхода биомассы при использовании клетками различных источников энергии. Эти значения совпадают с экспериментальными данными.
Ключевые слова: кинетические модели, законы сохранения, биоэнергетика, выход биомассы, скорость потребления источника энергии, поддержание клеток.
The stoichiometry of metabolic pathways in the dynamics of cellular populations
Computer Research and Modeling, 2011, v. 3, no. 4, pp. 455-475Views (last year): 5. Citations: 1 (RSCI).The problem has been considered, to what extent the kinetic models of cellular metabolism fit the matter which they describe. Foundations of stoichiometry of the whole metabolism and its large regions have been stated. A bioenergetic representation of stoichiometry based on a universal unit of chemical compound reductivity, viz., redoxon, has been described. Equations of mass-energy balance (bioenergetic variant of stoichiometry) have been derived for metabolic flows including those of protons possessing high electrochemical potential μH+, and high-energy compounds. Interrelations have been obtained which determine the biomass yield, rate of uptake of energy source for cell growth and other important physiological quantities as functions of biochemical characteristics of cellular energetics. The maximum biomass energy yield values have been calculated for different energy sources utilized by cells. These values coincide with those measured experimentally.
-
Математическая модель гидридного фазового перехода в частице порошка симметричной формы
Компьютерные исследования и моделирование, 2012, т. 4, № 3, с. 569-584В статье предложена математическая модель фазового перехода на примере гидрирования/дегидрирования порошка металла. Рассматривается одна частица, форма которой обладает некоторой симметрией. Шар, цилиндр и плоская пластина являются частными случаями симметричных форм. Модель описывает как сценарий «сжимающегося ядра» (формирование слоя новой фазы на поверхности частицы с его последующим утолщением), так и сценарий «образования и роста зародышей», при которых сплошной слой не формируется до полного исчезновения старой фазы. Модель представляет собой неклассическую диффузионную краевую задачу со свободной границей и нелинейными граничными условими III рода. Предположения симметрии позволяют свести задачу к одной пространственной переменной. Модель апробирована на серии экспериментальных данных. Показано, что влияние формы частиц на кинетику несущественно. Также показано, что ансамбль частиц различных форм с распределением по размерам может быть аппроксимирован одной частицей «среднего» размера простой формы, что оправдывает использование в моделях упрощающих предположений.
Ключевые слова: гидрирование, дегидрирование, фазовый переход, математическое моделирование, симметрия формы.
Mathematical model of hydride phase change in a symmetrical powder particle
Computer Research and Modeling, 2012, v. 4, no. 3, pp. 569-584Views (last year): 2. Citations: 2 (RSCI).In the paper we construct the model of phase change. Process of hydriding / dehydriding is taken as an example. A single powder particle is considered under the assumption about its symmetry. A ball, a cylinder, and a flat plate are examples of such symmetrical shapes. The model desribes both the "shrinking core"(when the skin of the new phase appears on the surface of the particle) and the "nucleation and growth"(when the skin does not appear till complete vanishing of the old phase) scenarios. The model is the non-classical boundary-value problem with the free boundary and nonlinear Neumann boundary condition. The symmetry assumptions allow to reduce the problem to the single spatial variable. The model was tested on the series of experimental data. We show that the particle shape’s influence on the kinetics is insignificant. We also show that a set of particles of different shapes with size distribution can be approxomated by the single particle of the "average" size and of a simple shape; this justifies using single particle approximation and simple shapes in mathematical models.
-
Об исследовании плазменной турбулентности на основе анализа спектров
Компьютерные исследования и моделирование, 2012, т. 4, № 4, с. 793-802В статье рассмотрены примеры анализа спектров экспериментальных данных для выявления характерных структур процессов, формирующих турбулентность в плазме. Основу метода составляет использование оригинального алгоритма, идеологически близкого к бутстреппроцедуре для одновыборочной задачи. В качестве базовой модели для описания тонкой структуры стохастических процессов предлагаются конечные сдвиг-масштабные смеси нормальных законов. Для отыскания статистических оценок (максимального правдоподобия) предполагается использование широко известного EM-алгоритма. Для нескольких серий спектров, полученных для разных режимов низкочастотной плазменной турбулентности, демонстрируется эффективность использования предложенного метода исследования.
On the investigation of plasma turbulence by the analysis of the spectra
Computer Research and Modeling, 2012, v. 4, no. 4, pp. 793-802Views (last year): 2. Citations: 4 (RSCI).The article describes the examples of the analysis of the experimental data spectra for identifying typical structures of processes forming plasma turbulence. The method is based on the original algorithm which is close to the one-sample bootstrap. The base model for description of the fine structure of stochastic processes is finite local-scale normal mixtures. For finding the statistical estimates (maximum likelihood estimates) well known EM algorithm is used. The efficiency of the proposed research technique is demonstrated for a number of spectra’s set obtained in different modes of low-frequency plasma turbulence.
-
Модель деформации полимерных нанокомпозитов на основе клеточных автоматов
Компьютерные исследования и моделирование, 2014, т. 6, № 1, с. 131-136Данная статья посвящена моделированию процесса деформации полимерных нанокомпозитов, содержащих «жесткие» и «мягкие» включения, с использованием клеточных автоматов и параллельных вычислений. В статье описан алгоритм расчета по модели, приведены сравнения с экспериментальными данными и описан программный комплекс для проведения численного эксперимента.
Deformation model of polymer nanocomposites based on cellular automata
Computer Research and Modeling, 2014, v. 6, no. 1, pp. 131-136Views (last year): 3. Citations: 2 (RSCI).This paper discusses the modeling of the deformation of polymer nanocomposites containing "hard" and "soft" inclusions, using cellular automata and parallel computing. The paper describes an algorithm based on the model, a comparison with experimental data is shown, software for the numerical experiment is described.
-
Компонентный анализ бинарных сред методом акустической рефлектоимпедансометрии
Компьютерные исследования и моделирование, 2015, т. 7, № 2, с. 301-313Рассматривается компьютерная модель компонентного анализа бинарных сред, основанная на применении нового метода — акустической рефлектоимпедансометрии и реализованная в среде графического программирования LabVIEW. Обсуждаются перспективы метрологического и приборного обеспечения экспериментальных приложений модели.
Ключевые слова: ультразвук, диагностика, неразрушающий контроль, коэффициент отражения, акустический импеданс.
Component analysis of binary media using acoustic reflectoimpedancemetry
Computer Research and Modeling, 2015, v. 7, no. 2, pp. 301-313Citations: 4 (RSCI).A computer model of component analysis of binary media, based on application of a new method acoustic reflecto-impedancemetry and realized in graphic programming environment LabVIEW is considered. Prospects of metrological and instrumental provisions of experimental applications of the model are discussed.
-
Расчетные исследования процесса перемешивания неизотермических потоков натриевого теплоносителя в тройнике
Компьютерные исследования и моделирование, 2017, т. 9, № 1, с. 95-110В программном комплексе FlowVision проведено численное моделирование процесса перемешивания неизотермических потоков натриевого теплоносителя в тройнике для обоснования применимости различных подходов — URANS (Unsteady Reynolds Averaged Navier Stokers), LES (Large Eddy Simulation) и квази-DNS (Direct Numerical Simulation) — для предсказания осциллирующего характера течения в зоне смешения и получения температурных пульсаций. Одна из основных задач данной работы — выявление преимуществ и недостатков использования этих подходов.
Численное исследование пульсаций температуры, возникающих в жидкости и в стенках тройника в процессе перемешивания неизотермических потоков натриевого теплоносителя, проведено в рамках математической модели, предполагающей, что рассматриваемое течение турбулентное, плотность жидкости не зависит от давления и что между теплоносителем и стенками тройника происходит теплообмен. При моделировании турбулентного теплопереноса в рамках подхода URANS применялась модель турбулентного теплопереноса LMS.
Исследование было проведено в два этапа. На предварительном этапе были определены влияние расчетной сетки на формирование осциллирующего течения и характер температурных пульсаций в рамках указанных выше подходов к моделированию турбулентности. В результате этого исследования были выработаны критерии построения расчетных сеток для каждого из подходов и произведена оценка потребных вычислительных ресурсов.
Затем были проведены расчеты для трех режимов течения, отличающихся соотношением расходов и температур натрия во входных сечениях тройника. Для каждого режима выполнены расчеты с применением подходов URANS, LES и квази-DNS.
На заключительном этапе работы был проведен сравнительный анализ численных и экспериментальных данных. Определены и сформулированы преимущества и недостатки использования каждого из указанных подходов к моделированию процесса перемешивания неизотермических потоков натриевого теплоносителя в тройнике.
Ключевые слова: подходы URANS, LES, DNS к моделированию турбулентности, натриевый теплоноситель, неизотермические потоки, пульсации температуры.
Numerical investigations of mixing non-isothermal streams of sodium coolant in T-branch
Computer Research and Modeling, 2017, v. 9, no. 1, pp. 95-110Views (last year): 3.Numerical investigation of mixing non-isothermal streams of sodium coolant in a T-branch is carried out in the FlowVision CFD software. This study is aimed at argumentation of applicability of different approaches to prediction of oscillating behavior of the flow in the mixing zone and simulation of temperature pulsations. The following approaches are considered: URANS (Unsteady Reynolds Averaged Navier Stokers), LES (Large Eddy Simulation) and quasi-DNS (Direct Numerical Simulation). One of the main tasks of the work is detection of the advantages and drawbacks of the aforementioned approaches.
Numerical investigation of temperature pulsations, arising in the liquid and T-branch walls from the mixing of non-isothermal streams of sodium coolant was carried out within a mathematical model assuming that the flow is turbulent, the fluid density does not depend on pressure, and that heat exchange proceeds between the coolant and T-branch walls. Model LMS designed for modeling turbulent heat transfer was used in the calculations within URANS approach. The model allows calculation of the Prandtl number distribution over the computational domain.
Preliminary study was dedicated to estimation of the influence of computational grid on the development of oscillating flow and character of temperature pulsation within the aforementioned approaches. The study resulted in formulation of criteria for grid generation for each approach.
Then, calculations of three flow regimes have been carried out. The regimes differ by the ratios of the sodium mass flow rates and temperatures at the T-branch inlets. Each regime was calculated with use of the URANS, LES and quasi-DNS approaches.
At the final stage of the work analytical comparison of numerical and experimental data was performed. Advantages and drawbacks of each approach to simulation of mixing non-isothermal streams of sodium coolant in the T-branch are revealed and formulated.
It is shown that the URANS approach predicts the mean temperature distribution with a reasonable accuracy. It requires essentially less computational and time resources compared to the LES and DNS approaches. The drawback of this approach is that it does not reproduce pulsations of velocity, pressure and temperature.
The LES and DNS approaches also predict the mean temperature with a reasonable accuracy. They provide oscillating solutions. The obtained amplitudes of the temperature pulsations exceed the experimental ones. The spectral power densities in the check points inside the sodium flow agree well with the experimental data. However, the expenses of the computational and time resources essentially exceed those for the URANS approach in the performed numerical experiments: 350 times for LES and 1500 times for ·DNS.
Indexed in Scopus
Full-text version of the journal is also available on the web site of the scientific electronic library eLIBRARY.RU
The journal is included in the Russian Science Citation Index
The journal is included in the RSCI
International Interdisciplinary Conference "Mathematics. Computing. Education"




