All issues

В работе рассматриваются интегро-дифференциальные уравнения влагопереноса дробного порядка с оператором Бесселя. Изучаемые уравнения содержат оператор Бесселя, два оператора дробного дифференцирования Герасимова – Капуто с разными порядками $\alpha$ и $\beta$. Рассмотрены два вида интегро-дифференциальных уравнений: в первом случае уравнение содержит нелокальный источник, т.е. интеграл от неизвестной функции по переменной интегрирования $x$, а во втором — случае интеграл по временной переменной $\tau$, обозначающий эффект памяти. Подобные задачи возникают при изучении процессов с предысторией. Для решения дифференциальных задач при различных соотношениях $\alpha$ и $\beta$ получены априорные оценки в дифференциальной форме, откуда следуют единственность и устойчивость решения по правой части и начальным данным. Для приближенного решения поставленных задач построены разностные схемы с порядком аппроксимации $O(h^2+\tau^2)$ при $\alpha=\beta$ и $O(h^2+\tau^{2-\max\{\alpha,\beta\}})$ при $\alpha\neq\beta$. Исследование единственности, устойчивости и сходимости решения проводится с помощью метода энергетических неравенств. Получены априорные оценки решений разностных задач при различных соотношениях $\alpha$ и $\beta$, откуда следуют единственность и устойчивость, а также сходимость решения разностной схемы к решению исходной дифференциальной задачи со скоростью равной порядку аппроксимации разностной схемы.

The paper considers integro-differential equations of fractional order moisture transfer with the Bessel operator. The studied equations contain the Bessel operator, two Gerasimov – Caputo fractional differentiation operators with different orders $\alpha$ and $\beta$. Two types of integro-differential equations are considered: in the first case, the equation contains a non-local source, i.e. the integral of the unknown function over the integration variable $x$, and in the second case, the integral over the time variable τ, denoting the memory effect. Similar problems arise in the study of processes with prehistory. To solve differential problems for different ratios of $\alpha$ and $\beta$, a priori estimates in differential form are obtained, from which the uniqueness and stability of the solution with respect to the right-hand side and initial data follow. For the approximate solution of the problems posed, difference schemes are constructed with the order of approximation $O(h^2+\tau^2)$ for $\alpha=\beta$ and $O(h^2+\tau^{2-\max\{\alpha,\beta\}})$ for $\alpha\neq\beta$. The study of the uniqueness, stability and convergence of the solution is carried out using the method of energy inequalities. A priori estimates for solutions of difference problems are obtained for different ratios of $\alpha$ and $\beta$, from which the uniqueness and stability follow, as well as the convergence of the solution of the difference scheme to the solution of the original differential problem at a rate equal to the order of approximation of the difference scheme.

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

The study of the structural characteristics of composites and nanostructures is of fundamental importance in materials science. Theoretical and numerical modeling and simulation of the mechanical properties of nanostructures is the main tool that allows for complex studies that are difficult to conduct only experimentally. One example of nanostructures considered in this work are carbon nanotubes (CNTs), which have good thermal and electrical properties, as well as low density and high Young’s modulus, making them the most suitable reinforcement element for composites, for potential applications in aerospace, automotive, metallurgical and biomedical industries. In this review, we reviewed the modeling methods, mechanical properties, and applications of CNT-reinforced metal matrix composites. Some modeling methods applicable in the study of composites with polymer and metal matrices are also considered. Methods such as the gradient descent method, the Monte Carlo method, methods of molecular statics and molecular dynamics are considered. Molecular dynamics simulations have been shown to be excellent for creating various composite material systems and studying the properties of metal matrix composites reinforced with carbon nanomaterials under various conditions. This paper briefly presents the most commonly used potentials that describe the interactions of composite modeling systems. The correct choice of interaction potentials between parts of composites directly affects the description of the phenomenon being studied. The dependence of the mechanical properties of composites on the volume fraction of the diameter, orientation, and number of CNTs is detailed and discussed. It has been shown that the volume fraction of carbon nanotubes has a significant effect on the tensile strength and Young’s modulus. The CNT diameter has a greater impact on the tensile strength than on the elastic modulus. An example of works is also given in which the effect of CNT length on the mechanical properties of composites is studied. In conclusion, we offer perspectives on the direction of development of molecular dynamics modeling in relation to metal matrix composites reinforced with carbon nanomaterials.

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

Amid the ongoing trend of rapid urbanization and the intensive development of megacities and large cities worldwide, the impact of man-made vibrations on residential structures and infrastructure is increasing. The operation of subway systems, construction using pile-driving and drilling equipment, and heavy traffic have become active sources of wave disturbances, which can be a decisive factor in reducing the structural stability of buildings and, consequently, their long-term reliability. This paper proposes a numerical calculation using the grid-characteristic method to model elastic waves propagating through soil layers and load-bearing structures from various sources. By solving the direct problem of numerical pulse simulation and varying its location, the values of velocity vector projections and components of the Cauchy stress tensor were obtained at each time step. Two scenarios were examined: the first simulates the impact of noise generated by construction work or nearby traffic, while the second demonstrates how a subway running through an underground tunnel affects multi-story residential buildings. Wave propagation patterns from these sources were visualized in terms of the parameters of interest, enabling a quick and convenient comprehensive analysis of the problem. The analysis of the obtained data will help adjust the timing and types of repair work, identify structural weak points, and develop innovative methods for preserving historical buildings that are cultural heritage sites. Additionally, it will allow for the most economically optimal construction of modern buildings near architectural landmarks, provide an efficient and safe action plan in emergencies, and modernize existing construction technologies to enhance the comfort of residential buildings, office structures, and other socially significant facilities. It will also aid in selecting the most suitable locations for modern high-precision manufacturing plants.

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

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

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

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

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

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

A mathematical model of fluid movement in the portal vein is considered. The fundamental circumstance determining the whole variety of portal hemodynamic phenomena is the absence of a valve apparatus. The flow direction is solely a function of the pressure gradient, and, therefore, it is fundamentally reversible when the boundary conditions of the system change.

Calculations were performed in the area representing a CT image fragment of the portal vein of a particular patient, which does not contain vascular bifurcations. The interpolation of patient Dopplerography data published in the press was used as boundary conditions for the flow. The calculations were carried out using the FlowVision industrial hydrodynamics software package. A comparison of calculations using the ideal fluid model and the Kuemada model of viscoplastic flow is carried out.

Calculations were performed for different values of the resistance coefficient corresponding to the physiological norm and with an increased value of the resistance coefficient.

At a normal value of the resistance coefficient, the flow in the portal vein is characterized by strong convective mixing.

As a result of calculations, it was found that when using the Kuemada model, the flow in the portal vein is stratified. The nature of the stratification depends on the hematocrit. At a normal value of the resistance coefficient, a plastic core of the flow is formed as the velocity decreases. With an increased value of the resistance coefficient during flow inversion, the flow core is also formed. When the flow is inverted, the plastic core continues to move in the forward direction, while the wall layers of the liquid begin to move in the opposite direction.

It is noted that in order to obtain correct modeling results, it is necessary to refine the blood composition of the portal vein and, possibly, refine the rheological model. Such information can be obtained both from comparing simulation data with clinical data, and by laboratory examination of portal vein blood.

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

Isothermal electroconvection in a dielectric liquid arising in a plane-parallel electrode system due to unipolar injection of charges from the cathode is considered. Spatially periodic rolls structures stability is investigated.

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

The paper describes the results of computer simulation of time-dependent temperature fields arising in polar dielectrics irradiated by focused electron bunches with average electron energy when analyzing with electron microscopy techniques. The mathematical model was based on solving several-dimensional nonstationary heat conduction equation with use of numerical finite element method. The approximation of thermal source was performed taking into account the estimation of initial electron distribution determined by Monte-Carlo simulation of electron trajectories. The simulation program was designed in Matlab. The geometrical modeling and calculation results demonstrated the main features of model sample heating by electron beam were presented at the given experimental parameters as well as source approximation.

В данной работе решается динамическая задача теории упругости о совместном нормальном и касательном воздействии на полупространство. С помощью этой задачи моделируется процесс наклонного падения метеорита на земную поверхность. Проведены исследования и расчеты поверхностной волны Рэлея. Полученное решение использовано в качестве внешнего воздействия на высотное здание, находящееся на некотором расстоянии от места падения для оценки безопасности и устойчивости его конструкции. Проведены численные эксперименты на основе конечно-элементного программного комплекса STARK ES. Рассчитаны амплитуды колебаний верхних этажей выбранного объекта при таком динамическом воздействии. Также проведено ихсистема тическое сравнение с результатами расчета при колебаниях основания, соответствующихст андартной акселерограмме 8-балльного землетрясения.

In this paper the dynamic elasticity problem of the simultaneous normal and tangential impact on the half-space is solved. This problem simulates the oblique incidence of meteorite on the Earth’s surface. The surface Rayleigh wave is investigated. The resulting solution is used as an external effect on the high-rise building, located at some distance from the spot of falling for the safety and stability assessment of its structure. Numerical experiments were made based on the finite element software package STARK ES. Upper floors amplitudes of the selected object were calculated under such dynamic effects. Also a systematic comparison with the results at the foundation vibrations, relevant to  standard a 8-point earthquake accelerograms, was made.

Построена статическая трехуровневая теоретико-игровая модель системы управления судовыми балластными водами. Используются методы иерархического управления при одновременном учете условий поддержания экологической системы в заданном состоянии. Проводится сравнение результатов исследования модели с точки зрения игр Гермейера $\Gamma_1$ и $\Gamma_2$. Приведены примеры численных расчетов в ряде характерных случаев.

The static three-level game-theoretic model of a control system of the ship’s water ballast is built. The methods of hierarchical control in view of requirements of keeping the system in the given state are used. A comparison of the results of study of the model in terms of $\Gamma_1$ and $\Gamma_2$. Germeyer’s games is conducted. Numerical calculations for some typical cases are given.

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

The aim of the present study is to show how engineering approaches and ideas work in clinical restorative dentistry, in particular, how they affect the restoration design and durability of restored endodontically treated teeth. For these purposes a 3D-computational model of a first incisor including the elements of hard tooth tissues, periodontal ligament, surrounding bone structures and restoration itself has been constructed and numerically simulated for a variety of restoration designs under normal chewing loadings. It has been researched the effect of different adhesive technologies in root canal on the functional characteristics of a restored tooth. The 3D model designed could be applied for preclinical diagnostics to determine the areas of possible fractures of a restored tooth and prognosticate its longevity.

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

The physical and mathematical model of combustion of the polydisperse suspension of coal dust was developed. The formulation of the problem takes into account the evaporation of particle volatile components during the heating, the particle emitting and the gas heat transfer to a surrounding area via the sphere volume side surface, heat transfer coefficient as a function of temperature. The polydisperse of coal-dust is taken into consideration. N — the number of fraction. Fractions are subdivided into inert and reacting particles. The oxidizer mass balance equation takes into consideration the oxidizer consumption per each reaction (heterogeneous on the particle surface and homogenous in the gas). Exothermic chemical reactions in gas are determined by Arrhenius equation with second-order kinetics. The heterogeneous reaction on the particle surface was first-order reaction. The numerical simulation was solved by Runge–Kutta–Merson method. Reliability of the calculations was verified by solving the partial problems. During the numerical calculation the percentage composition of inert and reacting particles in coal-dust and their total mass were changed for each simulation. We have determined the influence of the percentage composition of inert and reacting particles on burning characteristics of polydisperse coal-dust methane-air mixture. The results showed that the percent increase of volatile components in the mixture lead to the increase of total pressure in the volume. The value of total pressure decreases with the increasing of the inert components in the mixture. It has been determined that there is the extremism radius value of coarse particles by which the maximum pressure reaches the highest value.