Все выпуски
- 2024 Том 16
- 2023 Том 15
- 2022 Том 14
- 2021 Том 13
- 2020 Том 12
- 2019 Том 11
- 2018 Том 10
- 2017 Том 9
- 2016 Том 8
- 2015 Том 7
- 2014 Том 6
- 2013 Том 5
- 2012 Том 4
- 2011 Том 3
- 2010 Том 2
- 2009 Том 1
-
Моделирование неравномерного горения и напряженно-деформированного состояния пороховых элементов трубчатого заряда при выстреле
Компьютерные исследования и моделирование, 2022, т. 14, № 6, с. 1281-1300Врабо те представлена физико-математическая постановка задач внутренней баллистики артиллерийского выстрела для заряда, состоящего из совокупности пороховых трубок, и их напряженно-деформированного состояния. Горение и движение пучка пороховых трубок по каналу ствола моделируются эквивалентным трубчатым зарядом всестороннего горения. Предполагается, что эквивалентная трубка движется по оси канала ствола. Скорость движения эквивалентного трубчатого заряда и его текущее положение определяются из второго закона Ньютона. При расчете параметров течения использованы двумерные осесимметричные уравнения газовой динамики, для решения которых строится осесимметричная ортогонализированная разностная сетка, адаптирующаяся к условиям течения. Для численного решения системы газодинамических уравнений применяется метод контрольного объема. Параметры газа на границах контрольных объемов определяются с использованием автомодельного решения задачи о распаде произвольного разрыва С. К. Годунова. Напряженно-деформированное состояние моделируется для отдельной горящей пороховой трубки, находящейся в поле нестационарных газодинамических параметров. Расчет газодинамических параметров выстрела осуществляется без учета деформированного состояния пороховых элементов. При данных условиях рассмотрено поведение пороховых элементов при выстреле. Для решения нестационарной задачи упругости используется метод конечных элементов с разбиением области расчета на треугольные элементы. В процессе выгорания пороховой трубки расчетная сетка на каждом временном слое динамической задачи полностью обновляется в связи с изменением границ порохового элемента за счет горения. Представлены временные зависимости параметров внутрибаллистического процесса и напряженно-деформированного состояния пороховых элементов, а также распределения основных параметров течения продуктов горения в различные моменты времени. Установлено, что трубчатые пороховые элементы в процессе выстрела испытывают существенные деформации, которые необходимо учитывать при решении основной задачи внутренней баллистики. Полученные данные дают представления об уровне эквивалентных напряжений, действующих в различных точках порохового элемента. Представленные результаты говорят об актуальности сопряженной постановки задачи газовой динамики и напряженно-деформированного состояния для зарядов, состоящих из трубчатых порохов, поскольку это позволяет по-новому подойти к проектированию трубчатых зарядов и открывает возможность определения параметров, от которых существенно зависят физика процесса горения пороха и, следовательно, динамика процесса выстрела.
Ключевые слова: артиллерийская система, внутренняя баллистика, математическое моделирование, осесимметричная постановка, уравнения газовой динамики, трубчатые пороховые элементы, уравнения напряженно-деформированного состояния.
Simulation of uneven combustion and stress-strain state of powder elements of a tubular charge during firing
Computer Research and Modeling, 2022, v. 14, no. 6, pp. 1281-1300The paper presents the physical and mathematical formulation of the problems of internal ballistics of an artillery shot for a charge consisting of a set of powder tubes and their stress-strain state. Combustion and movement of a bundle of powder tubes along the barrel channel is modeled by an equivalent tubular charge of all-round combustion. It is assumed that the equivalent tube moves along the axis of the bore. The speed of movement of an equivalent tubular charge and its current position are determined from Newton’s second law. When calculating the flow parameters, two-dimensional axisymmetric equations of gas dynamics were used, for the solution of which an axisymmetric orthogonalized difference grid is constructed, which adapts to the flow conditions. The control volume method is used to numerically solve the system of gas-dynamic equations. The gas parameters at the boundaries of the control volumes are determined using a self-similar solution to the Godunov’s problem of the decay of an arbitrary discontinuity. The stress-strain state is modeled for a separate burning powder tube located in the field of gas-dynamic parameters. The calculation of the gas-dynamic parameters of the shot is carried out without taking into account the deformed state of the powder elements. The behavior of powder elements during firing is considered under these conditions. The finite element method with the division of the calculation area into triangular elements is used to solve the problem of elasticity. In the process of powder tube burnout, the computational grid on each time layer of the dynamic problem is completely updated due to a change in the boundaries of the powder element due to combustion. The paper shows the time dependences of the parameters of the internal ballistics process and the stress-strain state of powder elements, as well as the distribution of the main parameters of the flow of combustion products at different points in time. It has been established that the tubular powder elements during the shot experience significant deformations, which must be taken into account when solving the basic problem of internal ballistics. The data obtained give an idea of the level of equivalent stresses acting at various points of the powder element. The results obtained indicate the relevance of the conjugate formulation of the problem of gas dynamics and the stress-strain state for charges consisting of tubular powders, since this allows a new approach to the design of tubular charges and opens up the possibility of determining the parameters on which the physics of the combustion process of gunpowder significantly depends, therefore, and the dynamics of the shot process.
-
Разработка методического подхода и численное моделирование теплогидравлических процессов в промежуточном теплообменнике реактора БН
Компьютерные исследования и моделирование, 2023, т. 15, № 4, с. 877-894В работе представлены результаты трехмерного численного моделирования теплогидравлических процессов в промежуточном теплообменнике перспективного реактора на быстрых нейтронах с натриевым теплоносителем (БН) с учетом разработанного методического подхода.
Промежуточный теплообменник (ПТО) размещен в корпусе реактора и предназначен для передачи тепла от натрия первого контура, циркулирующего в межтрубном пространстве, натрию второго контура, циркулирующему внутри труб. Перед входными окнами ПТО при интегральной компоновке оборудования первого контура в реакторе БН имеет место температурное расслоение теплоносителя из-за неполного перемешивания разнотемпературных потоков на выходе из активной зоны. Внутри ПТО в районе входных и выходных окон теплообменника также реализуется сложное продольно-поперечное течение теплоносителя, которое приводит к неравномерному распределению расхода теплоносителя в межтрубном пространстве и, как следствие, к неравномерному распределению температуры и эффективности теплообмена по высоте и радиусу трубного пучка.
С целью подтверждения заложенных в проекте теплогидравлических параметров ПТО перспективного реактора БН был разработан методический подход для трехмерного численного моделирования теплообменника, размещенного в корпусе реактора, учитывающий трехмерную картину течения натрия на входе и внутри ПТО, а также обосновывающий рекомендации для упрощения геометрии расчетной модели ПТО. Численное моделирование теплогидравлических процессов в ПТО перспективного реактора БН проводилось с использованием программного комплекса FlowVision со стандартной $k-\varepsilon$-моделью турбулентности и моделью турбулентного теплопереноса LMS. Для повышения представительности численного моделирования трубного пучка ПТО выполнены верификационные расчеты однотрубного и многотрубного теплообменников «натрий – натрий» с соответствующими конструкции ПТО геометрическими характеристиками. Для определения входных граничных условий в модели ПТО выполнен дополнительный трехмерный расчет с учетом неравномерной картины течения в верхней смесительной камере реактора. Расчетная модель ПТО была оптимизирована за счет упрощения дистанционирующих поясов и выбора секторной модели. В результате численного моделирования ПТО получены распределения скорости натрия первого контура, температуры натрия первого и второго контуров. Удовлетворительное согласование результатов расчета с проектными данными по интегральным параметрам подтвердило принятые проектные теплогидравлические характеристики ПТО перспективного реактора БН.
Ключевые слова: методический подход, натрий, промежуточный теплообменник, реактор БН, сопряженный теплообмен, теплогидравлика, численное моделирование.
Development of a methodological approach and numerical simulation of thermal-hydraulic processes in the intermediate heat exchanger of a BN reactor
Computer Research and Modeling, 2023, v. 15, no. 4, pp. 877-894The paper presents the results of three-dimensional numerical simulation of thermal-hydraulic processes in the Intermediate Heat Exchanger of the advanced Sodium-Cooled Fast-Neutron (BN) Reactor considering a developed methodological approach.
The Intermediate Heat Exchanger (IHX) is located in the reactor vessel and intended to transfer heat from the primary sodium circulating on the shell side to the secondary sodium circulating on the tube side. In case of an integral layout of the primary equipment in the BN reactor, upstream the IHX inlet windows there is a temperature stratification of the coolant due to incomplete mixing of different temperature flows at the core outlet. Inside the IHX, in the area of the input and output windows, a complex longitudinal and transverse flow of the coolant also takes place resulting in an uneven distribution of the coolant flow rate on the tube side and, as a consequence, in an uneven temperature distribution and heat transfer efficiency along the height and radius of the tube bundle.
In order to confirm the thermal-hydraulic parameters of the IHX of the advanced BN reactor applied in the design, a methodological approach for three-dimensional numerical simulation of the heat exchanger located in the reactor vessel was developed, taking into account the three-dimensional sodium flow pattern at the IHX inlet and inside the IHX, as well as justifying the recommendations for simplifying the geometry of the computational model of the IHX.
Numerical simulation of thermal-hydraulic processes in the IHX of the advanced BN reactor was carried out using the FlowVision software package with the standard $k-\varepsilon$ turbulence model and the LMS turbulent heat transfer model.
To increase the representativeness of numerical simulation of the IHX tube bundle, verification calculations of singletube and multi-tube sodium-sodium heat exchangers were performed with the geometric characteristics corresponding to the IHX design.
To determine the input boundary conditions in the IHX model, an additional three-dimensional calculation was performed taking into account the uneven flow pattern in the upper mixing chamber of the reactor.
The IHX computational model was optimized by simplifying spacer belts and selecting a sector model.
As a result of numerical simulation of the IHX, the distributions of the primary sodium velocity and primary and secondary sodium temperature were obtained. Satisfactory agreement of the calculation results with the design data on integral parameters confirmed the adopted design thermal-hydraulic characteristics of the IHX of the advanced BN reactor.
-
Исследование путей распространения потенциала действия у высших растений с использованием модели ФитцХью-Нагумо
Компьютерные исследования и моделирование, 2011, т. 3, № 1, с. 77-84Двумерная система электрически связанных возбудимых клеток, описанных моделью ФитцХью-Нагумо, была использована для теоретического анализа механизмов распространения потенциала действия (ПД) в тканях высших растений. Показано, что в системе, состоящей из одинаковых элементов, увеличение электрической проводимости между клетками повышало как скорость распространения, так и порог возбуждения ПД. Для имитации симпласта проводящих пучков растений была использована система, состоящая из элементов со слабой электрической связью, которые описывали паренхимные клетки проводящих пучков, и элементов с сильной электрической связью, которые описывали ситовидные элементы. При этом показано,
что порог возбуждения ПД приближался к порогу системы, состоящей только из элементов со слабой электрической связью, а скорость распространения сигнала была значительно выше, нежели в такой системе.
An investigation of an action potential propagation in vascular plant using FitzHugh-Nagumo model
Computer Research and Modeling, 2011, v. 3, no. 1, pp. 77-84Просмотров за год: 3. Цитирований: 4 (РИНЦ).A two-dimensional system of excitable cells, describing by the FitzHugh-Nagumo model, has been used for a theoretical investigation of an action potential propagation (AP) in vascular plant tissues. It is shown that growth of electrical conductivity between cells increases the AP generation threshold and its propagation velocity in the homogeneous system, which has been formed by equal elements. The plant symplast has been
described by the heterogeneous system, including elements with low electrical conductivity, which simulate parenchyma cells, and elements with high electrical conductivity, which simulate sieve elements. Analysis of this system shows that the threshold of the AP generation is similar with this threshold in the homogeneous system
with low electrical conductivity; the velocity of the AP propagation is faster than one in this system. -
Влияние поверхности на особенности аморфизации системы Ni-Ag
Компьютерные исследования и моделирование, 2014, т. 6, № 2, с. 263-269В рамках метода молекулярной динамики, с использованием многочастичного потенциаламежатомного взаимодействия, на модели свободной частицы Ni60Ag40 диаметром 40 Å, исследована структурная эволюция в процессе закалки. Проведен сравнительный анализ структурных перестроек для моделей частицы и массивного образца. Отмечено снижение температуры начала и конца превращений для частицы. Показано, что в образовании перколяционного кластера из взаимопроникающих и контактирующих между собой икосаэдров, для модели свободной частицы, задействовано на 10 % атомов больше, чем для модели массивного тела.
Effect of the surface on characteristics of amorphization Ni-Ag system
Computer Research and Modeling, 2014, v. 6, no. 2, pp. 263-269Просмотров за год: 1. Цитирований: 1 (РИНЦ).Molecular dynamics simulation using the embedded-atom method is applied to study the structural evolution of the particle diameter of 40 Å during the quenching process. Was carried comparative analysis of the structural reconstruction for the particle and the bulk models. Was a reduction in temperature of the beginning and end of the transformation of the particle. In formation of a percolation cluster from interpenetrating and contacting icosahedrons, for model of the particle, it is involved for 10 percent of atoms more, than for model of a bulk.
-
Топологический анализ микроструктуры жидкой воды на примере модели TIP4P-EW
Компьютерные исследования и моделирование, 2014, т. 6, № 3, с. 415-426Проведено молекулярно-динамическое (МД) моделирование жесткой TIP4P-EW модели воды при нормальных условиях. В качестве структурных элементов системы рассматривались симплексы Делоне, а для выделения ее плотной части использовался топологический критерий, который позволяет идентифицировать на мгновенных снимках МД ячейки микроструктуру системы молекул воды. Геометрический анализ симплексов Делоне системы указывает на сильную их уплощенность по сравнению с правильным тетраэдром, что принципиально отличается от результатов для простых жидкостей. Статистика кластеров плотной части мгновенного снимка системы исследовалась в зависимости от их мощности и связности. Она схожа с таковой для простых жидкостей, а структура этой плотной части также представляет собой фрактальную поверхность, состоящую из свободных граней симплексов Делоне.
Topological microstructure analysis of the TIP4P-EW water model
Computer Research and Modeling, 2014, v. 6, no. 3, pp. 415-426Просмотров за год: 1. Цитирований: 1 (РИНЦ).Molecular dynamics (MD) simulations of rigid water model TIP4P-EW at ambient conditions were carried out. Delaunay’s simplexes were considered as structural elements of liquid water. Topological criterion which allows to identify the water microstructure in snapshot of MD cell was used to allocate its dense part. Geometrical analysis of water Delaunay’s simplexes indicates their strong flatness in comparison with a regular tetrahedron that is fundamentally different from the results for dense part of simple liquids. The statistics of TIP4P-EW water clusters was investigated depending on their cardinality and connectivity. It is similar to the statistics for simple liquids and the structure of this dense part is also a fractal surface consisting of the free edges of the Delaunay’s simplexes.
-
Перколяционное моделирование гидравлического гистерезиса в пористой среде
Компьютерные исследования и моделирование, 2014, т. 6, № 4, с. 543-558В работе рассматриваются различные модели гидравлического гистерезиса, возникающего при инвазивной ртутной порометрии. Для моделирования гидравлического гистерезиса используется изотропная перколяция узлов на трехмерных квадратных решетках с $(1,\,\pi)$-окрестностью. Феноменологически исследуется взаимосвязь данных инвазивной порометрии с параметрами перколяционной модели. Реализация перколяционной модели основана на библиотеках SPSL и SECP, выпущенных под лицензией GNU GPL-3 с использованием свободного языка программирования R.
Ключевые слова: инвазивная ртутная порометрия, гидравлический гистерезис, перколяция узлов, квадратная решетка, неметрическое расстояние Минковского, окрестность Мура, массовая фрактальная размерность, язык программирования R, библиотека SPSL, библиотека SECP.
Percolation modeling of hydraulic hysteresis in a porous media
Computer Research and Modeling, 2014, v. 6, no. 4, pp. 543-558Просмотров за год: 3. Цитирований: 1 (РИНЦ).In this paper we consider various models of hydraulic hysteresis in invasive mercury porosimetry. For simulating the hydraulic hysteresis is used isotropic site percolation on three-dimensional square lattices with $(1,\,\pi)$-neighborhood. The relationship between the percolation model parameters and invasive porosimetry data is studied phenomenologically. The implementation of the percolation model is based on libraries SPSL and SECP, released under license GNU GPL-3 using the free programming language R.
-
Модель анизотропной прочности со скалярным параметром поврежденности
Компьютерные исследования и моделирование, 2014, т. 6, № 6, с. 937-942В работе обсуждается возможность моделирования анизотропии прочности слоистой упругой среды с использованием скалярного параметра поврежденности. Сформулированы термодинамически согласованные определяющие уравнения. С помощью пакета SIMULIA/Abaqus моделируется процесс растяжения и сжатия образцов. Результаты расчета с использованием предложенной модели сравниваются с известными из литературы экспериментальными данными и предсказаниями традиционных моделей.
Ключевые слова: разрушение, теория поврежденности, анизотропия, термодинамическая согласованность, пользовательские подпрограммы.
Modeling of anisotropic strength using scalar damage parameter
Computer Research and Modeling, 2014, v. 6, no. 6, pp. 937-942Просмотров за год: 1.The paper discusses the possibility of modeling the strength anisotropy of layered elastic medium using a scalar damage parameter. Thermodynamically consistent constitutive equations are formulated. Using SIMULIA / Abaqus we numerically simulated the stretching and compression of the samples. The results of calculation using the proposed model are compared with the known experimental data from the literature and the predictions of traditional models.
-
Особенности численного моделирования поля модифицированного магнита спектрометра
Компьютерные исследования и моделирование, 2015, т. 7, № 1, с. 93-105В данной работе представлены результаты численного моделирования распределения магнитного поля некоторых модификаций спектрометрического магнита СП-40, используемого в экспериментальной установке НИС ОИЯИ. Основной особенностью такого магнита является прямоугольная апертура, следовательно, и область, в которой решается краевая задача, имеет всюду гладкую границу, за исключением угловой области ферромагнетика. В таких случаях решение задачи или производные решения могут иметь особенность. Изучено поведение магнитного поля в окрестности угловой области ферромагнетика и численным путем выбрана конфигурация магнита, для которой ширина области однородности (т. е. с $\Delta B/B < 1\,\%$) магнитного поля увеличилась с 0.5 м до 1.0 м, т. е. в два раза.
Ключевые слова: магнитные системы, математическое моделирование, краевая задача, эллиптические уравнения, поведение решения в угловой области ферромагнетика.
Peculiar features of numerical modeling of the modified spectrometer magnet field
Computer Research and Modeling, 2015, v. 7, no. 1, pp. 93-105Цитирований: 1 (РИНЦ).In this work, we propose a method of the numerical solution of the magnetostatic problem for domains with boundaries containing corners. With the help of this numerical method, the magnetic systems of rectangular configuration were simulated with high accuracy. In particular, the calculations of some modifications of the magnetic system SP-40 used in the NIS JINR experimental installation, are presented. The basic feature of such a magnet is a rectangular aperture, hence, the area in which the boundary-value problem is solved, has a smooth border everywhere, except for a finite number of angular points in the vicinity of which the border is formed by crossing two smooth curves. In such cases the solution to the problem or derivatives of the solution can have a special feature. A behavior of the magnetic field in the vicinity of an angular point is investigated, and the configuration of the magnet was chosen numerically. The width of the area of homogeneity of the magnetic field increased from 0.5 m up to 1.0 m, i. e. twice.
-
Построение адаптивной сетки в окрестности «угловой точки» ферромагнетика в численном моделировании магнитной системы
Компьютерные исследования и моделирование, 2016, т. 8, № 1, с. 89-99При численном решении краевой задачи магнитостатики в области с негладкой границей возникает вопрос о точности получаемого решения в окрестности угловой точки ферромагнетика [Zhidkov, Perepelkin, 2001]. В окрестности «угловой точки» возможен существенный рост модуля магнитного поля, что приводит к необходимости построения специальных численных алгоритмов при решении краевой задачи. В данной работе предложен алгоритм построения адаптивной сетки в окрестности угловой точки ферромагнетика, учитывающий характер поведения решения краевой задачи. Приводится пример расчета модельной задачи в области, содержащей угловую точку.
Ключевые слова: магнитные системы, математическое моделирование, краевая задача, эллиптические уравнения, поведение решения в угловой точке.
Construction of adaptive mesh in the domain with boundary «corner point» of ferromagnetic in the numerical simulation of magnetic systems
Computer Research and Modeling, 2016, v. 8, no. 1, pp. 89-99Просмотров за год: 2.At numerical solving of the boundary-value problem of magnetostatic in a domain with a boundary corner point, a question of accuracy of the obtained solution near the corner point of ferromagnetic arises [Zhidkov, Perepelkin, 2001]. Near the corner point an essential growth of the module of the magnetic field can take place, which leads to the necessity of constructing special numerical algorithms when solving the boundary-value problem. This work represents an algorithm of constructing an adaptive mesh in the domain with a boundary corner point of ferromagnetic taking into account the character of behaviour of the solution of the boundary-value problem. An example of calculating a model problem in the domain containing a corner point is given.
-
Работы А. С. Комарова по клеточно-автоматному моделированию популяционно-онтогенетических процессов у растений
Компьютерные исследования и моделирование, 2016, т. 8, № 2, с. 285-295Рассмотрены возможности моделирования в технике клеточных автоматов применительно к травянистым растениям и кустарничкам. Приводятся основные положения дискретного описания онтогенезов растений, на которых основывается математическое моделирование. В обзоре обсуждаются основные результаты, полученные с использованием моделей и раскрывающие закономерности функционирования ценопопуляций и сообществ. Описана модель CAMPUS и результаты компьютерного эксперимента по разрастанию двух клонов брусники с разной геометрией побегов. Публикация посвящена работам профессора А. С. Комарова, основоположника направления; дан список его основных публикаций по этой тематике.
Ключевые слова: компьютерные модели, индивидуально-ориентированный подход.
A.S. Komarov’s publications about cellular automata modelling of the population-ontogenetic development in plants: a review
Computer Research and Modeling, 2016, v. 8, no. 2, pp. 285-295The possibilities of cellular automata simulation applied to herbs and dwarf shrubs are described. Basicprinciples of discrete description of the ontogenesis of plants on which the mathematical modeling based are presents. The review discusses the main research results obtained with the use of models that revealing the patterns of functioning of populations and communities. The CAMPUS model and the results of computer experiment to study the growth of two clones of lingonberry with different geometry of the shoots are described. The paper is dedicated to the works of the founder of the direction of prof. A. S. Komarov. A list of his major publications on this subject is given.
Keywords: computer models, individual-based approach.Просмотров за год: 2. Цитирований: 6 (РИНЦ).
Журнал индексируется в Scopus
Полнотекстовая версия журнала доступна также на сайте научной электронной библиотеки eLIBRARY.RU
Журнал входит в систему Российского индекса научного цитирования.
Журнал включен в базу данных Russian Science Citation Index (RSCI) на платформе Web of Science
Международная Междисциплинарная Конференция "Математика. Компьютер. Образование"