Все выпуски

В данной работе рассматривается взаимодействие кинка уравнения $\varphi^4$ с двумя протяженными одинаковыми примесями. Протяженная примесь описывается с помощью функции прямоугольного вида. Анализируется случай притягивающей примеси. С помощью аналитических методов рассматривается случай малых амплитуд локализованных волн, когда возможно провести линеаризацию уравнений движения. Для численного решения использовался метод прямых для уравнений в частных производных. Для нахождения частот колебаний, локализованных на примесях волн, используется дискретное преобразование Фурье. Кинк запускался в направлении примесей с разными начальными скоростями. Изменялось также расстояние между двумя примесями. Показано, что при взаимодействии кинка с примесями на них возбуждаются долгоживущие локализованные волны бризерного типа. Исследована их структура и связанная динамика. Определено, как, изменяя параметры примесей и расстояние между ними, можно управлять типом и динамическими параметрами связанных колебаний, локализованных на примесях волн. Найдены возможные решения в виде синфазных, антифазных колебаний, в виде биений. Колебания локализованных волн происходят с излучением волн малой амплитуды. Спектр этих излучений состоит из двух частот. Первая приближенно равна $\sqrt{2}$, что соответствует величине частоты для хвоста воблингбризера уравнения $\varphi^4$. Вторая приближенно равна удвоенной частоте колебаний примесных мод. Найдено (как аналитически, так и численно) наличие двух возможных частот для связанных локализованных колебаний. Показано, что частоты сильно зависят от расстояния между примесями. С увеличением расстояния между примесями частоты сливаются в одну — частоту, полученную для случая одиночной примеси. Найденные численно и аналитически зависимости частот от расстояния между примесями хорошо совпадают для больших расстояний, когда взаимодействие между примесями слабое, и начинают заметно отличаться при малых расстояниях, когда взаимодействие между примесями сильное. Аналитическое значение величин полученных частот всегда больше численных. Показано, что зависимость амплитуды локализованных волн от начальной скорости кинка имеет несколько минимумов и максимумов.

In this paper, we consider the interaction of a kink of the $\varphi^4$ equation with two identical extended impurities. An extended impurity is described using a rectangular function. The case of an attractive impurity is analyzed. Using analytical methods, we consider the case of small amplitudes of localized waves, when it is possible to linearize the equations of motion. For the numerical solution, the method of lines for partial differential equations was used. To find the oscillation frequencies of waves localized on impurities, the discrete Fourier transform is used. The kink was launched in the direction of the impurities with different initial velocities. The distance between the two impurities was also varied. It is shown that when a kink interacts with impurities, long-lived localized breather-type waves are excited on them. Their structure and coupled dynamics are investigated. It is determined how, by changing the parameters of the impurities and the distance between them, it is possible to control the type and dynamic parameters of the coupled oscillations of the waves localized on the impurities. Possible solutions in the form of in-phase, antiphase oscillations, in the form of beats are found. The oscillations of localized waves occur with the emission of small-amplitude waves. The spectrum of these emissions consists of two frequencies. The first is approximately equal to $\sqrt{2}$, which corresponds to the frequency value for the wobbling breather tail of the $\varphi^4$ equation. The second is approximately equal to the doubled frequency of impurity mode oscillations. The presence of two possible frequencies for coupled localized oscillations is found both analytically and numerically. It is shown that the frequencies strongly depend on the distance between impurities. With increasing distance between impurities, the frequencies merge into one — frequency obtained for the case of a single impurity. The dependences of the frequencies on the distance between impurities found numerically and analytically coincide well for large distances, when the interaction between impurities is weak, and begin to differ noticeably at small distances, when the interaction between impurities is strong. The analytical value of the obtained frequencies is always greater than the numerical ones. It is shown that the dependence of the amplitude of localized waves on the initial kink velocity has several minima and maxima.

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

The mathematical model for aeroelastic oscillations of a narrow channel wall with a nonlinear-elastic suspension and interacting with a pulsating viscous gas layer is proposed. Within the framework of this model, the aeroelastic response of the channel wall and its phase response were determined and investigated. The authors simultaneously studied the influence of the nonlinear stiffness elastic suspension of the wall, compressibility and dissipative properties of gas, as well as the inertia of its motion on the wall oscillations. The model was elaborated based on the formulation and solution of the initial boundary-value plane problem of mathematical physics. The problem governing equations include the equations of dynamics for barotropic viscous gas, equation of dynamics for the rigid wall as the spring-mass nonlinear oscillator. Using the perturbation method, the asymptotic analysis of the problem was carried out. The solution of the equations of dynamics for the thin layer of viscous gas was obtained by the iteration method. As a result, the law of gas pressure distribution in the channel was determined and the initial problem of aeroelasticity was reduced to the study of the generalized Duffing equation. Its solution was realized by the harmonic balance method, which allowed us to determine the aeroelastic and phase responses of the channel wall in the form of implicit functions. The numerical study of these responses was carried out to evaluate the influence for inertia of gas motion and its compressibility, as well as a comparison of the results obtained with the special cases of creeping motion of viscous gas and incompressible viscous fluid. The results of this study have shown the importance of simultaneous consideration of compressibility and inertia of viscous gas motion when modeling aeroelastic oscillations of the considered channel wall.

В данной работе описана методика генерации синтетической базы данных стенозов, состоящей из 1620 записей. Каждая запись представляет собой результаты численного эксперимента по моделированию трехмерного течения вязкой несжимаемой жидкости через трубку с переменным сечением: перепад давлений, средний поток, усредненная по сечению скорость кровотока на входе в трубку, максимальная степень сужения стеноза, длина стеноза, асимметрия стеноза, радиус трубки, число Рейнольдса. База данных валидировалась путем сравнения с другими моделями (с эластичными стенками) и стендовыми экспериментами и показала отклонение перепадов давлений не выше 4%. База данных синтетических стенозов использовалась для обучения физически информированной нейронной сети для быстрой оценки перепада давления по четырем ключевым входным показателям: число Рейнольдса, длина стеноза, степень стеноза, степень асимметрии стеноза. Физическая информированность достигалась за счет введения штрафов в функцию потерь за отсутствие положительного перепада давления и за отсутствие монотонности перепада давления по входным параметрам. Физически информированная нейронная сеть показала более высокую точность на гемодинамических значимых стенозах при тестировании на валидационной выборке и на новых стенозах, не представленных в базе данных. Средняя относительная ошибка на стенозах длиной в 8 радиусов здорового сосуда составила 6% для физически информированной сети и 13% для классической нейронной сети. Ошибки на коротких стенозах длиной в 4 радиуса оказались почти одинаковы: 9,5% для физически информированной сети и 10% для классической нейронной сети. Разработанный метод функциональной оценки гемодинамической значимости стенозов может использоваться как самостоятельный инструмент для клинической оценки стенозов и как компонент сетевых моделей кровотока. Наибольшую актуальность подход приобретает при моделировании многососудистых поражений, которые доминируют в клинической практике. Ключевое преимущество метода заключается в физической корректности результатов и точности, сопоставимой с классическим моделированием, при значительно меньших вычислительных затратах.

This paper describes a method for generating a synthetic database of stenoses, consisting of 1620 entries. Each entry represents the results of a numerical experiment simulating the three-dimensional flow of a viscous incompressible fluid through a tube with a variable cross-section: pressure drop, mean flow rate, cross-sectionally averaged inlet blood flow velocity, maximum stenosis severity, stenosis length, stenosis asymmetry, tube radius, and Reynolds number. The database was validated by comparison with other models (with elastic walls) and bench experiments, showing a deviation in pressure drops of no more than 4%. The synthetic stenosis database was used to train a physics-informed neural network for the rapid estimation of pressure drop based on four key input parameters: Reynolds number, stenosis length, stenosis severity, and stenosis asymmetry coefficient. The physics-informed aspect was achieved by introducing penalties into the loss function for the absence of a positive pressure drop and for the lack of monotonicity of the pressure drop with respect to the input parameters. The physics-informed neural network demonstrated higher accuracy on hemodynamically significant stenoses when tested on a validation set and on new stenoses not represented in the database. The mean relative error for stenoses with a length of 8 healthy vessel radii was 6% for the physics-informed network and 13% for a classical neural network. The errors for short stenoses with a length of 4 radii were nearly identical: 9.5% for the physics-informed network and 10% for the classical neural network. The developed method for the functional assessment of the hemodynamic significance of stenoses can be used both as a standalone tool for clinical stenosis evaluation and as a component of network blood flow models. The approach becomes most relevant when modeling multi-vessel disease, which is predominant in clinical practice. The key advantage of the method lies in the physical correctness of the results and accuracy comparable to classical modeling, but with significantly lower computational costs.

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

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

В ходе исследования были изучены волновые и структурные характеристики разлета концентрированных газовых взвесей в зависимости от объемной доли частиц, их размеров и начального отношения давлений. Было установлено, что показатель политропы и скорость звука в газовой взвеси зависят не только от температуры, но и от давления и концентрации частиц. С увеличением давления значения показателя политропы и скорости звука возрастают, а с увеличением объемной доли дисперсной фазы — уменьшаются. При распаде произвольного разрыва с увеличением начального давления наблюдается необычное, по сравнению с «чистой» газодинамикой, уменьшение относительной скорости смеси в зоне постоянного течения. Это объясняется нелинейной зависимостью скорости звука в газодисперсной смеси от давления.

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

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

Modeling of gas-particle suspensions with large pressure and density gradients is of practical interest in the study of volcanic phenomena, explosions at different altitudes, as well as in technogenic problems related to the operation of space technology and the formation of space debris. This work presents numerical and analytical investigations of the expansion of gas suspensions with a high ratio (up to six orders of magnitude) of pressures and densities. For numerical modeling, a high-resolution hybrid large-particle method was employed. Under the conditions considered, the accuracy of the method was confirmed by comparison with asymptotically exact solutions. The study examined the wave and structural characteristics of concentrated gas suspension expansion depending on particle volume fraction, particle size, and initial pressure ratio. It was found that the polytropic index and sound speed in the gas suspension depend not only on temperature but also on pressure and particle concentration. With increasing pressure, both the polytropic index and sound speed rise, while with increasing particle volume fraction they decrease. In the case of an arbitrary discontinuity decay, an unusual effect is observed compared with “pure” gas dynamics: the relative velocity of the mixture in the uniform flow region decreases as the initial pressure increases. This is explained by the nonlinear dependence of the sound speed in a gas-dispersed mixture on pressure. With increasing particle size (Stokes number), the mixture flow splits into gaseous and dispersed components. At the initial moment, the contact discontinuity separating the mixture from the rarefied gas region splits into two contact boundaries: gaseous and dispersed. A practical conclusion is that when the particle size changes by two orders of magnitude, the gas-dynamic parameters of the mixture in the rarefaction wave region and up to the medium interface remain close to each other. During spatial expansion, the initial cylindrical shape of the dispersed medium successively transforms into a cross-section resembling a hexagon. At the next stage of expansion, the particles redistribute to form a bilateral conical structure. Eventually, a dispersed formation close to a spherical shape emerges.

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

In this paper a numerical model for blood flow through a venous bifurcation with an obliterating clot is investigated. We studied propagation of perturbations of blood flow velocity and perturbations of pressure inside the vein. The model is built in acoustic (linear) approximation. Computational results reveal conditions for clot resonance oscillation, which can cause its detachment and thromboembolism.

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

The aim of research is to develop software system for solving gas dynamic problem in multiply connected integration domains of regular shape by high-performance computing system. Comparison of the various technologies of parallel computing has been done. The program complex is implemented using multithreaded parallel systems to organize both multi-core and massively parallel calculation. The comparison of numerical results with known model problems solutions has been done. Research of performance of different computing platforms has been done.

Методом математического моделирования на основе метода крупных частиц исследуется динамика распространения пожара в однородном степном массиве.

We consider the two-dimensional mathematical model of wildfire. Numerical solution algorithm based on the method of large particles was developed for this model.

The paper is an example of computer simulation in mechanical engineering. Velocity field in a hydrocyclone are determined numerically, because for direct measurements it is difficult to achieve them. The numerical simulation of 3D fluid dynamics based on the k-eps RNG model of turbulence in the hydrocyclone with the injector, containing 5 tangentially directed nozzles is considered. It is shown that the direction of movement of
injected fluid in the hydrocyclone depends on the water flow rate through the injector. The calculations show in accordance with the experiments that the dependence of the Split-parameter on the injected water flow rate has a non-monotone character associated with the ratio of power of the main flow and the injected fluid.

Статья представляет собой пример компьютерного моделирования в области инженерной механики. Численным методом находятся поля скорости в гидроциклоне, которые недоступны прямому измерению. Рассматривается численное моделирование трехмерной гидродинамики на основе k-ε RNG модели турбулентности в гидроциклоне со встроенным инжектором, содержащим 5 тангенциально направленных сопла. Показано, что направление движения инжектируемой жидкости зависит от расхода жидкости через инжектор. Расчеты показывают в соответствии с экспериментами, что зависимость сплит-параметра от расхода инжектируемой жидкости имеет немонотонный характер, связанный с отношением мощности основного потока и инжектируемой жидкости.

В данной работе решается задача поиска конструкции магнитной системы для создания магнитного поля с требуемыми характеристиками в заданной области. На основе анализа математической модели магнитной системы предлагается достаточно общий подход к решению нелинейной обратной задачи, которая описывается уравнением Фредгольма H(z) = ∫_SIJ(s)G(z, s)ds, z ∈ S _H, s ∈ S _I . Необходимо определить распределение плотности тока J(s), а также расстановку источников тока для создания поля H(z). В работе предлагается метод решения этих задачс помощью регуляризованных итерационных процессов. На примере конкретной магнитной системы проводится численное исследование влияния различных факторов на характер создаваемого магнитного поля.

In this paper the problem of searching for the design of the magnetic system for creation a magnetic field with the required characteristics in the given area is solved. On the basis of analysis of the mathematical model of the magnetic system rather a general approach is proposed to the solving of the inverse problem, which is written by the Fredgolm equation H(z) = ∫_SIJ(s)G(z, s)ds, z ∈ S _H, s ∈ S _I . It was necessary to define the current density distribution function J(s) and the existing winding geometry for creation of a required magnetic field H(z). In the paper a method of solving those by means of regularized iterative processes is proposed. On the base of the concrete magnetic system we perform the numerical study of influence of different factors on the character of the magnetic field being designed.

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

The purpose of this paper is to generalize the macroscopic hydrodynamic vehicular traffic models by using the algorithm for constructing the adequate state equation — dependence the pressure from traffic density by taking into account the real experimental data (possibly using the parametric solutions for model equations). It is proved that this kind of state equation which closed model equations system and obtained from the experimentally observed form of the fundamental diagram — dependence the traffic intensity from its density, completely determines the all properties of the used phenomenological model.