Все выпуски

Производство высокочистого водорода необходимо для экологически чистой энергетики и различных химико-технологических процессов. Значительная часть водорода будет производиться за счет конверсии метана. Методом измерения удельной водородопроницаемости исследуются различные сплавы, перспективные для использования в газоразделительных установках. Требуется оценить параметры диффузии и сорбции, чтобы иметь возможность численно моделировать различные сценарии и условия эксплуатации материала (включая экстремальные), выделять лимитирующие факторы. В статье представлены нелинейная модель водородопроницаемости в соответствии со спецификой эксперимента, численный метод решения краевой задачи и результаты параметрической идентификации модели для сплава V₈₅Ni₁₅.

High-purity hydrogen is required for clean energy and a variety of chemical technology processes. A considerable part of hydrogen is to be obtained by methane conversion. Different alloys, which may be wellsuited for use in gas-separation plants, were investigated by measuring specific hydrogen permeability. One had to estimate the parameters of diffusion and sorption to numerically model the different scenarios and experimental conditions of the material usage (including extreme ones), and identify the limiting factors. This paper presents a nonlinear model of hydrogen permeability in accordance with the specifics of the experiment, the numerical method for solving the boundary-value problem, and the results of parametric identification for the alloy V₈₅Ni₁₅.

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

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.

Теоретические и экспериментальные исследования взаимодействия водяного пара с пористыми материалами проводятся как на макро-, так и на микроуровне. На макроуровне исследуется влияние структуры расположения индивидуальных пор на процессы взаимодействия водяного пара с пористым материалом как сплошной средой. На микроуровне исследуется зависимость характеристик взаимодействия водяного пара с пористой средой от геометрии и размеров индивидуальной поры.

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

Theoretical and experimental investigations of water vapor interaction with porous materials are carried out both at the macro level and at the micro level. At the macro level, the influence of the arrangement structure of individual pores on the processes of water vapor interaction with porous material as a continuous medium is studied. At the micro level, it is very interesting to investigate the dependence of the characteristics of the water vapor interaction with porous media on the geometry and dimensions of the individual pore.

In this paper, a study was carried out by means of mathematical modelling of the processes of water vapor interaction with suffering pore of the cylindrical type. The calculations were performed using a model of a hybrid type combining a molecular-dynamic and a macro-diffusion approach for describing water vapor interaction with an individual pore. The processes of evolution to the state of thermodynamic equilibrium of macroscopic characteristics of the system such as temperature, density, and pressure, depending on external conditions with respect to pore, were explored. The dependence of the evolution parameters on the distribution of the diffusion coefficient in the pore, obtained as a result of molecular dynamics modelling, is examined. The relevance of these studies is due to the fact that all methods and programs used for the modelling of the moisture and heat conductivity are based on the use of transport equations in a porous material as a continuous medium with known values of the transport coefficients, which are usually obtained experimentally.

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

The background social strain of a society can be quantitatively estimated using various statistical indicators. Mathematical models, allowing to forecast the dynamics of social strain, are successful in describing various social processes. If the number of interacting groups is small, the dynamics of the corresponding indicators can be modelled with a system of ordinary differential equations. The increase in the number of interacting components leads to the growth of complexity, which makes the analysis of such models a challenging task. A continuous social stratification model can be considered as a result of the transition from a discrete number of interacting social groups to their continuous distribution in some finite interval. In such a model, social strain naturally spreads locally between neighbouring groups, while in reality, the social elite influences the whole society via news media, and the Internet allows non-local interaction between social groups. These factors, however, can be taken into account to some extent using the term of the model, describing negative external influence on the society. In this paper, we develop a continuous social stratification model, describing the dynamics of two societies connected through migration. We assume that people migrate from the social group of donor society with the highest strain level to poorer social layers of the acceptor society, transferring the social strain at the same time. We assume that all model parameters are constants, which is a realistic assumption for small societies only. By using the finite volume method, we construct the spatial discretization for the problem, capable of reproducing finite propagation speed of social strain. We verify the discretization by comparing the results of numerical simulations with the exact solutions of the auxiliary non-linear diffusion equation. We perform the numerical analysis of the proposed model for different values of model parameters, study the impact of migration intensity on the stability of acceptor society, and find the destabilization conditions. The results, obtained in this work, can be used in further analysis of the model in the more realistic case of inhomogeneous coefficients.

Рассматривается пространственно-временная модель тритрофической системы, описывающая взаимодействие жертвы, хищника и суперхищника в среде с неоднородным распределением ресурса. Учитываются всеядность суперхищника (Intraguild Predation, IGP), диффузия и направленная миграция (таксис), который моделируется с помощью логарифмической функции от ресурса и плотности жертвы. Основное внимание уделено анализу мультистабильности системы и роли косимметрии в формировании континуальных семейств стационарных решений. С использованием численно-аналитического подхода изучаются пространственно-однородные и неоднородные стационарные решения. Установлено, что при выполнении дополнительных соотношений между параметрами, характеризующими локальное взаимодействие хищников, и коэффициентами диффузии система обладает косимметрией, что приводит к возникновению семейства устойчивых стационарных решений, пропорциональных функции ресурса. Показано, что косимметрия не зависит от функции ресурса в случае неоднородной среды. Проведено исследование устойчивости стационарных распределений с помощью спектрального метода. Нарушение условий косимметрии приводит к разрушению семейства и появлению изолированных стационарных состояний, а также к длительным переходным процессам, отражающим память системы об исчезнувшем семействе. В зависимости от начальных условий и параметров в системе реализуются переходы к режимам с одним хищником (выживание хищника или суперхищника) или к сосуществованию хищников. Численные эксперименты на основе метода прямых (разностная схема по пространственной переменной и метод Рунге – Кутты для интегрирования по времени) подтверждают мультистабильность системы и иллюстрируют исчезновение семейства решений при разрушении косимметрии.

We consider a spatiotemporal model of a tritrophic system describing the interaction between prey, predator, and superpredator in an environment with nonuniform resource distribution. The model incorporates superpredator omnivory (Intraguild Predation, IGP), diffusion, and directed migration (taxis), the latter modeled using a logarithmic function of resource availability and prey density. The primary focus is on analyzing the multistability of the system and the role of cosymmetry in the formation of continuous families of steady-state solutions. Using a numerical-analytical approach, we study both spatially homogeneous and inhomogeneous steady-state solutions. It is established that under additional relations between the parameters governing local predator interactions and diffusion coefficients, the system exhibits cosymmetry, leading to the emergence of a family of stable steady-state solutions proportional to the resource function. We demonstrate that the cosymmetry is independent of the resource function in the case of a heterogeneous environment. The stability of stationary distributions is investigated using spectral methods. Violation of the cosymmetry conditions results in the breakdown of the solution family and the emergence of isolated equilibria, as well as prolonged transient dynamics reflecting the system’s “memory” of the vanished states. Depending on initial conditions and parameters, the system exhibits transitions to single-predator regimes (survival of either the predator or superpredator) or predator coexistence. Numerical experiments based on the method of lines, which involves finite difference discretization in space and Runge –Kutta integration in time, confirm the system’s multistability and illustrate the disappearance of solution families when cosymmetry is broken.

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

We consider a model of populations that are closely related and share a common areal. System of nonlinear parabolic equations is formulated that incorporates nonlinear diffusion and migration flows induced by nonuniform densities of population and carrying capacity. We employ the method of lines and study the impact of migration on scenarios of local competition and coexistence of species. Conditions on system parameters are determined when a nontrivial family of steady states is formed.

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

Population dynamics is considered in a modified chemostat model including diffusion, chemotaxis, and nonlocal competitive losses. To account for influence of the external environment on the population of the ecosystem, a random parameter is included into the model equations. Computer simulations reveal three dynamic modes depending on system parameters: the transition from initial state to a spatially homogeneous steady state, to a spatially inhomogeneous distribution of population density, and elimination of population density.

На основе математическоймо дели динамики биомасс фитопланктона построен способ оценки содержания хлорофилла в районе моря с учетом его распределения по глубине. Модель построена на основе уравнения «реакция-диффузия», учитывает основные влияющие факторы: минеральное питание, освещенность и температуру. Используется спутниковая информация о поверхностном слое моря. Приведен пример расчетов для залива Петра Великого (Японское море).

A model, describing the influence of external factors on temporal evolution of phytoplankton distribution in a horizontally-homogenous water layer, is presented. This model is based upon the reactiondiffusion equation and takes into account the main factors of influence: mineral nutrients, insolation and temperature. The mineral nutrients and insolation act oppositely on spatial phytoplankton distribution. The results of numerical modeling are presented and the prospect of applying this model to reconstruction of phytoplankton distribution from sea-surface satellite data is discussed. The model was used to estimate the chlorophyll content of the Peter the Great Bay (Sea of Japan).

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

The article contains results of modeling of bi-modal microstructure evolution of VT6 alloy during isothermal forging. The model of recrystallization based on dislocation approach to nucleation and growth of nuclear was made to calculate recrystallization processes of the secondary alpha phase. The globularization process of lamellar alpha phase was calculated with assumption of diffusion-control migration of beta-phase boundary due to grain boundary diffusion of vanadium atom. The theoretical results were compared with experimental one.

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

Metal hydrides are an interesting class of chemical compounds that can reversibly bind a large amount of hydrogen and are, therefore, of interest for energy applications. Understanding the factors affecting the kinetics of hydride formation and decomposition is especially important. Features of the material, experimental setup and conditions affect the mathematical description of the processes, which can undergo significant changes during the processing of experimental data. The article proposes a general approach to numerical modeling of the formation and decomposition of metal hydrides and solving inverse problems of estimating material parameters from measurement data. The models are divided into two classes: diffusive ones, that take into account the gradient of hydrogen concentration in the metal lattice, and models with fast diffusion. The former are more complex and take the form of non-classical boundary value problems of parabolic type. A rather general approach to the grid solution of such problems is described. The second ones are solved relatively simply, but can change greatly when model assumptions change. Our experience in processing experimental data shows that a flexible software tool is needed; a tool that allows, on the one hand, building models from standard blocks, freely changing them if necessary, and, on the other hand, avoiding the implementation of routine algorithms. It also should be adapted for high-performance systems of different paradigms. These conditions are satisfied by the HIMICOS library presented in the paper, which has been tested on a large number of experimental data. It allows simulating the kinetics of formation and decomposition of metal hydrides, as well as related tasks, at three levels of abstraction. At the low level, the user defines the interface procedures, such as calculating the time layer based on the previous layer or the entire history, calculating the observed value and the independent variable from the task variables, comparing the curve with the reference. Special algorithms can be used for solving quite general parabolic-type boundary value problems with free boundaries and with various quasilinear (i.e., linear with respect to the derivative only) boundary conditions, as well as calculating the distance between the curves in different metric spaces and with different normalization. This is the middle level of abstraction. At the high level, it is enough to choose a ready tested model for a particular material and modify it in relation to the experimental conditions.