Все выпуски

В данной работе представлены новые результаты исследований, которые с 2011 года проводятся в Институте экономики РАН под руководством академика РАН В.И. Маевского. Эти работы направлены на развитие теории переключающегося воспроизводства и соответствующих математических моделей, особенностью которых является то, что в них в явном виде моделируется взаимодействие финансового и реального секторов экономики, а сама экономика страны дезагрегируется не по отраслевому принципу (машиностроение, сельское хозяйство, услуги и т. д.), а на производственные подсистемы, отличающиеся друг от друга возрастом основного капитала. Одна их математических сложностей работы с подобными моделями, называемыми моделями переключающегося режима воспроизводства (ПРВ), заключается в сложности моделирования конкурентных отношений между подсистемами разного возраста. Поэтому до сих пор в моделях ПРВ рассматривалось взаимодействие конечного числа производственных подсистем, сами модели имели дискретно-непрерывный характер, расчеты делались исключительно на ЭВМ, получение аналитических зависимостей было затруднено. В настоящей работе показано, что для частного случая сбалансированного экономического роста и континуального множества производственных подсистем возможно получение аналитических выражений, позволяющих лучше понять особенности влияния денежно-кредитной политики на экономическую динамику. Кроме чисто научного интереса это имеет важное практическое значение, поскольку позволяет оценивать возможную реакцию реального сектора экономики на изменения в денежно-кредитной сфере без проведения сложных имитационных расчетов.

This paper presents new research results that have been conducted at the Institute of Economics of the Russian Academy of Sciences since 2011 under the leadership of Academician of the Russian Academy of Sciences V. I.Mayevsky. These works are aimed at developing the theory of switching mode of reproduction and corresponding mathematical models, the peculiarity of which is that they explicitly model the interaction of the financial and real sectors of the economy, and the country’s economy itself is not disaggregated according to the sectoral principle (engineering, agriculture, services, etc.), but by production subsystems that differ from each other by the age of the fixed capital. One of the mathematical difficulties of working with such models, called models of switching mode of reproduction (SMR), is the difficulty of modeling competitive relationships between subsystems of different “ages”. Therefore, until now, the interaction of a finite number of production subsystems has been considered in the SMR models, the models themselves were of a discrete-continuous nature, calculations were done exclusively on computers, and obtaining analytical dependencies was difficult. This paper shows that for the special case of balanced economic growth and a continuum of production subsystems, it is possible to obtain analytical expressions that allow a better understanding of the impact of monetary policy on economic dynamics. In addition to purely scientific interest, this is of great practical importance, since it allows us to assess the possible reaction of the real sector of the economy to changes in the monetary sphere without conducting complex simulation calculations.

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

Для описания предлагается модель, в которой дискретная функция вероятности числа ветвей на узле повторяет динамику непрерывного гауссового распределения. Характеристики распределения Гаусса (математическое ожидание $x_0$, среднеквадратичное отклонение $\sigma$) являются параметрами модели. В выбранной постановке задача относится к проблематике ветвящихся случайных процессов (ВСП) в неоднородной модели Гальтона – Ватсона.

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

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

Предлагается второй (оригинальный) подход, основанный на введении в рассмотрение параметра порядка, равного обратной средней высоте дерева, и анализа его поведения. Установлено, что динамика такого параметра порядка в сечениях $\sigma = \text{const}$ с очень небольшими отличиями имеет вид распределения Ферми – Дирака ($\sigma$ выполняет ту же функцию, что и температура для распределения Ферми – Дирака, $x_0$ — функцию энергии). Для параметра порядка подобрано эмпирическое выражение, введен и рассчитан аналог химического потенциала, который и имеет смысл характерного масштаба параметра порядка, то есть тех значений $x_0$, при которых условно можно считать, что порядок сменяется беспорядком. Этот критерий положен в основу определе- ния границы условного перехода в данном подходе. Установлено, что эта граница соответствует средней высоте дерева, равной двум поколениям. На основании обнаруженных свойств предложены рекомендации для медицинских учреждений, позволяющие контролировать обеспечение конечности траектории пациентов.

Рассмотренная модель и метод ее описания с помощью условно-бесконечных деревьев имеют приложение ко многим иерархическим системам. К таким системам можно отнести сети маршрутизации интернет-соединений, бюрократические сети, торговые, логистические сети, сети цитирования, игровые стратегии, задачи популяционной динамики и пр.

The 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.

В статье представлена технология построения систем поддержки принятия врачебных решений (СППВР), основанная на сервис-контейнерах с использованием Docker и веб-интерфейсе, работающем непосредственно в браузере без установки специализированного программного обеспечения на рабочую станцию врача. Предложена модульная архитектура, где для каждого прикладного модуля формируется самостоятельный сервис-контейнер, объединяющий мини-веб-сервер, пользовательский интерфейс и вычислительные компоненты обработки медицинских изображений. Взаимодействие между браузером и серверной частью реализовано через постоянное двунаправленное соединение по веб-сокетам с бинарной сериализацией сообщений в формате MessagePack, что обеспечивает малые задержки и эффективную передачу больших объeмов данных. Для локального хранения изображений и результатов анализа применены средства браузера (IndexedDB) с оболочкой Dexie.js, что ускоряет повторный доступ к данным. Трeхмерная визуализация и базовые операции с DICOM-данными реализованы с использованием библиотек Three.js и AMI.js: такая связка обеспечивает интеграцию интерактивных элементов, возникающих в контексте задачи (аннотации, ориентиры, метки, 3D-модели), в трeхмерные медицинские изображения.

Серверные компоненты и функциональные модули собраны в виде набора взаимодействующих контейнеров, управляемых средствами Docker. Рассмотрены выбор базовых образов, подходы к минимизации контейнеров до уровня исполняемых файлов без внешних зависимостей, организация многоступенчатой сборки, включающей отдельный «сборочный» контейнер. Описан сервис-хаб, выполняющий запуск прикладных контейнеров по обращению пользователя, проксирование запросов, управление сессиями и перевод контейнера из общего режима в монопольный при начале вычислений. Приведены примеры прикладных модулей (оценка фракционного резерва кровотока, расчeт количественного отношения потока, расчет смыкания створок аортального клапана), показана интеграция React-интерфейса и трeхмерной сцены, а также представлены политика версионирования, автоматизированные проверки воспроизводимости результатов и порядок развeртывания на целевой площадке.

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

The article presents a technology for building clinical decision support systems (CDSS) based on service containers using Docker and a web interface that runs directly in the browser without installing specialized software on workstation of a clinician. A modular architecture is proposed in which each application module is packaged as an independent service container combining a lightweight web server, a user interface, and computational components for medical image processing. Communication between the browser and the server side is implemented via a persistent bidirectional WebSocket connection with binary message serialization (MessagePack), which provides low latency and efficient transfer of large data. For local storage of images and analysis of results, browser facilities (IndexedDB with the Dexie.js wrapper) are used to speed up repeated data access. Three-dimensional visualization and basic operations with DICOM data are implemented with Three.js and AMI.js: this toolchain supports the integration of interactive elements arising from the task context (annotations, landmarks, markers, 3D models) into volumetric medical images.

Server components and functional modules are assembled as a set of interacting containers managed by Docker. The paper discusses the choice of base images, approaches to minimizing containers down to runtime-only executables without external utilities, and the organization of multi-stage builds with a dedicated build container. It describes a hub service that launches application containers on user request, performs request proxying, manages sessions, and switches a container from shared to exclusive mode at the start of computations. Examples of application modules are provided (fractional flow reserve estimation, quantitative flow ratio computation, aortic valve closure modeling), along with the integration of a React-based interface with a three-dimensional scene, a versioning policy, automated reproducibility checks, and the deployment procedure on the target platform.

It is demonstrated that containerization ensures portability and reproducibility of the software environment, dependency isolation and scalability, while the browser-based interface provides accessibility, reduced infrastructure requirements, and interactive real-time visualization of medical data. Technical limitations are noted (dependence on versions of visualization libraries and data formats) together with practical mitigation measures.

Исследование упругих волн в пористых средах актуально для задач поиска полезных ископаемых, при использовании пористых экранов для демпфирования ударно-волновых воздействий, для изучения строения земной коры и т. д. Упругие свойства пористой среды, о которых можно судить по скорости распространения разного рода волн, зависят от степени консолидации пористой среды. Например, насыпные среды (песок, стеклянные шарики, гранулированные материалы) обладают низкой скоростью звука (порядка 100 м/с), уплотнение таких сред сопровождается некоторым увеличением скорости, а их консолидация (песчаник, цементация газогидратом) приводит к многократному возрас- танию скорости акустических волн, порядка 2000–3000 м/с. В данной работе теоретически исследуется динамика волнового импульса в ударной трубе, содержащей слой насыпной среды. Численное моделирование проведено для условий эксперимента. Приводится описание экспериментальной установки типа «ударная труба». Установка состоит из камеры высокого давления (КВД), камеры низкого давления (КНД) и секции насыпной среды. Ударно-волновой импульс (УВИ) создается из-за разрыва диафрагмы между КВД и КНД. Динамика УВИ регистрируется пьезоэлектрическими датчиками, расположенными заподлицо с внутренней стороны трубы. В ударной трубе, оснащенной секцией насыпных сред, волна испытывает многократные отражения от поверхности изучаемой пористой среды и верхнего торца трубы. Переотраженные сигналы используются в качестве зондирующих импульсов для изучения изменений в пористой среде, вызванных повторными прохождениями ударно-волнового импульса, с периодом около 10 мс. Используется математическая модель, включающая уравнения сохранения массы, уравнения сохранения импульсов и энергии для газовой фазы и твердых частиц с замыкающими соотношениями. Описание процесса проводится для одномерного плоского движения газовой и дисперсной фаз. Для численного решения используется аппроксимация уравнений, основанная на методе контрольного объема. Численные результаты показали, что предложенная модель качественно и количественно правильно описывает появление резкого кратковременного увеличения полного напряжения (пика) при повторном прохождении импульса через слой насыпной среды, наблюдаемого в экспериментах.

The 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.

Надежная диагностика болезней растений компьютерным зрением требует не только высокой точности классификации, но и стабильной обобщающей способности и интерпретируемого принятия решений. При этом механизмы внимания безусловно полезны для моделей глубокого обучения, но эффективность решения задачи зависит от места и способа интеграции этих механизмов в архитектуру нейронной сети. В данном исследовании представлен классификатор на базе оптимизированной архитектуры 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. Полученные результаты предоставляют практические рекомендации по проектированию надежных и интерпретируемых систем искусственного интеллекта для агропромышленного сектора.

Reliable 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. Моделируется броуновское движение нескольких сотен молекул флаводоксина, учитываются электростатические взаимодействия и сложная форма молекул. С помощью данной модели удалось воспроизвести экспериментальную немонотонную зависимость константы связывания флаводоксина с фотосистемой 1. Это говорит о том, что для описания такого вида зависимости достаточно учета только электростатических взаимодействий.

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.

В работе обсуждаются методические основы и проблемы моделирования мировой динамики. Излагаются подходы к построению новой имитационной модели глобального развития и первичные результаты моделирования. В основу построения модели положен эмпирический подход, основанный на анализе статистики основных социально-экономических показателей. На основании этого анализа выделены основные переменные. Для этих переменных составлены динамические уравнения (в непрерывно-дифференциальной форме). Связи между переменными подбирались исходя из динамики соответствующих показателей в прошлом и на основании экспертных оценок, при этом использовались эконометрические методы, основанные на регрессионном анализе. Были проведены расчеты по полученной системе динамических уравнений, результаты представлены в виде пучка траекторий для тех показателей, которые непосредственно наблюдаемы и по которым имеется статистика. Таким образом, имеется возможность оценить разброс траекторий и понять прогнозные возможности представленной модели.

The work discusses the methodological basis and problems of modeling of world dynamics. Outlines approaches to the construction of a new simulation model of global development and the results of the simulation. The basis of the model building is laid empirical approach which based on the statistical analysis of the main socio-economic indicators. On the basis of this analysis identified the main variables. Dynamic equations (in continuous differential form) were written for these variables. Dependencies between variables were selected based on the dynamics of indicators in the past and on the basis of expert assessments, while econometric techniques were used, based on regression analysis. Calculations have been performed for the resulting dynamic equations system, the results are presented in the form of a trajectories beam for those indicators that are directly observable, and for which statistics are available. Thus, it is possible to assess the scatter of the trajectories and understand the predictive capability of this model.

Методом молекулярной динамики в работе исследуется взаимодействие нелинейной локализованной колебательной моды с бивакансией Pt в кристалле Pt₃Al. Выявлены зависимости времени жизни нелинейной локализованной колебательной моды от начальной температуры модельного кристалла, от начального отклонения атома Al от положения равновесия, а также расстояния до внедренной бивакансии Pt в плоскости (111) кристалла.

By method of molecular dynamics investigated the interaction of nonlinear localized modes with bivacancies Pt crystal Pt₃Al. Identified dependences of the lifetime of the nonlinear localized modes from the initial temperature of the crystal model, the initial atom Al deviation from its equilibrium position, as well as the distance to the introduced bivacancy Pt in (111) plane of the crystal.

Рассматривается вертикально-распределенная трехкомпонентная модель морской экосистемы. Состояние планктонного сообщества с учетом питательных веществ анализируется в условиях активных перемещений зоопланктона в вертикальном столбе воды. Аналитически получены условия ДС-неустойчивости системы в окрестности пространственно-однородного равновесия. Численно определены области параметров, при которых пространственнооднородное равновесие устойчиво к небольшим пространственно-неоднородным возмущениям, неустойчиво по Тьюрингу и колебательно неустойчиво. Исследовано влияние параметров, определяющих биологические характеристики зоопланктона и пространственные перемещения планктона, на возможность образования пространственных структур. Показано, что при малой скорости потребления фитопланктона на пространственную неустойчивость существенно влияет убыль зоопланктона, а при больших значениях этого параметра имеют значение перемешивание фитопланктона и пространственные перемещения зоопланктона.

A vertically distributed three-component model of marine ecosystem is considered. State of the plankton community with nutrients is analyzed under the active movement of zooplankton in a vertical column of water. The necessary conditions of the Turing instability in the vicinity of the spatially homogeneous equilibrium are obtained. Stability of the spatially homogeneous equilibrium, the Turing instability and the oscillatory instability are examined depending on the biological characteristics of zooplankton and spatial movement of plankton. It is shown that at low values of zooplankton grazing rate and intratrophic interaction rate the system is Turing instable when the taxis rate is low. Stabilization occurs either through increased decline of zooplankton either by increasing the phytoplankton diffusion. With the increasing rate of consumption of phytoplankton range of parameters that determine the stability is reduced. A type of instability depends on the phytoplankton diffusion. For large values of diffusion oscillatory instability is observed, with a decrease in the phytoplankton diffusion zone of Turing instability is increases. In general, if zooplankton grazing rate is faster than phytoplankton growth rate the spatially homogeneous equilibrium is Turing instable or oscillatory instable. Stability is observed only at high speeds of zooplankton departure or its active movements. With the increase in zooplankton search activity spatial distribution of populations becomes more uniform, increasing the rate of diffusion leads to non-uniform spatial distribution. However, under diffusion the total number of the population is stabilized when the zooplankton grazing rate above the rate of phytoplankton growth. In general, at low rate of phytoplankton consumption the spatial structures formation is possible at low rates of zooplankton decline and diffusion of all the plankton community. With the increase in phytoplankton predation rate the phytoplankton diffusion and zooplankton spatial movement has essential effect on the spatial instability.

В программном комплексе FlowVision проведено численное моделирование процесса перемешивания неизотермических потоков натриевого теплоносителя в тройнике для обоснования применимости различных подходов — URANS (Unsteady Reynolds Averaged Navier Stokers), LES (Large Eddy Simulation) и квази-DNS (Direct Numerical Simulation) — для предсказания осциллирующего характера течения в зоне смешения и получения температурных пульсаций. Одна из основных задач данной работы — выявление преимуществ и недостатков использования этих подходов.

Численное исследование пульсаций температуры, возникающих в жидкости и в стенках тройника в процессе перемешивания неизотермических потоков натриевого теплоносителя, проведено в рамках математической модели, предполагающей, что рассматриваемое течение турбулентное, плотность жидкости не зависит от давления и что между теплоносителем и стенками тройника происходит теплообмен. При моделировании турбулентного теплопереноса в рамках подхода URANS применялась модель турбулентного теплопереноса LMS.

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

Затем были проведены расчеты для трех режимов течения, отличающихся соотношением расходов и температур натрия во входных сечениях тройника. Для каждого режима выполнены расчеты с применением подходов URANS, LES и квази-DNS.

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

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.