Все выпуски

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

Two algorithms of evaluation of the continuous wavelet transform with the Morlet wavelet are presented. The first one is the solution of PDE with transformed signal, which plays a role of the initial value. The second allows to explore the influence of central frequency variation via the diffusion smoothing of the data modulated by the harmonic functions. These approaches are illustrated by the analysis of the chaotic oscillations of the coupled Roessler systems.

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

The ring of Gaussian numbers is considered. The estimation of the module of some analogue of Weyl's trigonometrical sum with summation on Gaussian numbers is proved by methods of the analytical number theory. Multiplicative norm of Gaussian numbers is less than some integer.

Универсальные сценарии перехода к хаосу в динамических системах к настоящему моменту хорошо изучены. К типичным сценариям относятся каскад бифуркаций удвоения периода (сценарий Фейген-баума), разрушение тора малой размерности (сценарий Рюэля–Такенса) и переход через перемежаемость (сценарий Помо–Манневилля). В более сложных пространственно-распределенных динамических системах нарастающая с изменением параметра сложность поведения по времени тесно переплетается с формированием пространственных структур. Однако вопрос о том, могут ли в каком-то сценарии пространственная и временная оси полностью поменяться ролями, до сих пор остается открытым. В данной работе впервые предлагается математическая модель конвекции–реакции–диффузии, в рамках которой реализуется пространственный аналог перехода к хаосу через разрушение квазипериодического режима в рамках сценария Рюэля–Такенса. Исследуемая физическая система представляет собой два водных раствора кислоты (A) и основания (B), в начальный момент времени разделенных по пространству и помещенных в вертикальную ячейку Хеле–Шоу, находящуюся в статическом поле тяжести. При приведении растворов в контакт начинается фронтальная реакция нейтрализации второго порядка: A + B $\to$ C, которая сопровождается выделением соли (С). Процесс характеризуется сильной зависимостью коэффициентов диффузии реагентов от их концентрации, что приводит к возникновению двух локальных зон пониженной плотности, в которых независимо друг от друга возникают хемоконвективные движения жидкости. Слои, в которых развивается конвекция, все время остаются разделенными прослойкой неподвижной жидкости, но они могут влиять друг на друга посредством диффузии реагентов через прослойку. Формирующаяся хемо-конвективная структура представляет собой модулированную стоячую волну, постепенно разрушающуюся со временем, повторяя последовательность бифуркаций сценария разрушения двумерного тора. Показано, что в ходе эволюции системы пространственная ось, направленная вдоль фронта реакции, выполняет роль времени, а само время играет роль управляющего параметра.

In the last decades, universal scenarios of the transition to chaos in dynamic systems have been well studied. The scenario of the transition to chaos is defined as a sequence of bifurcations that occur in the system under the variation one of the governing parameters and lead to a qualitative change in dynamics, starting from the regular mode and ending with chaotic behavior. Typical scenarios include a cascade of period doubling bifurcations (Feigenbaum scenario), the breakup of a low-dimensional torus (Ruelle–Takens scenario), and the transition to chaos through the intermittency (Pomeau–Manneville scenario). In more complicated spatially distributed dynamic systems, the complexity of dynamic behavior growing with a parameter change is closely intertwined with the formation of spatial structures. However, the question of whether the spatial and temporal axes could completely exchange roles in some scenario still remains open. In this paper, for the first time, we propose a mathematical model of convection–diffusion–reaction, in which a spatial transition to chaos through the breakup of the quasi–periodic regime is realized in the framework of the Ruelle–Takens scenario. The physical system under consideration consists of two aqueous solutions of acid (A) and base (B), initially separated in space and placed in a vertically oriented Hele–Shaw cell subject to the gravity field. When the solutions are brought into contact, the frontal neutralization reaction of the second order A + B $\to$ C begins, which is accompanied by the production of salt (C). The process is characterized by a strong dependence of the diffusion coefficients of the reagents on their concentration, which leads to the appearance of two local zones of reduced density, in which chemoconvective fluid motions develop independently. Although the layers, in which convection develops, all the time remain separated by the interlayer of motionless fluid, they can influence each other via a diffusion of reagents through this interlayer. The emerging chemoconvective structure is the modulated standing wave that gradually breaks down over time, repeating the sequence of the bifurcation chain of the Ruelle–Takens scenario. We show that during the evolution of the system one of the spatial axes, directed along the reaction front, plays the role of time, and time itself starts to play the role of a control parameter.

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

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

The study of complex processes in various spheres of human activity is traditionally based on the use of mathematical models. In modern conditions, the development and application of such models is greatly simplified by the presence of high-speed computer equipment and specialized tools that allow, in fact, designing models from pre-prepared modules. Despite this, the known problems associated with ensuring the adequacy of the model, the reliability of the original data, the implementation in practice of the simulation results, the excessively large dimension of the original data, the joint application of sufficiency heterogeneous mathematical models in terms of complexity and integration of the simulated processes are becoming increasingly important. The more critical may be the external constraints imposed on the value of the optimized functional, and often unattainable within the framework of the constructed model. It is logical to assume that in order to fulfill these restrictions, a purposeful transformation of the original model is necessary, that is, the transition to a mathematical model with a deliberately improved solution. The new model will obviously have a different internal structure (a set of parameters and their interrelations), as well as other formats (areas of definition) of the source data. The possibilities of purposeful change of the initial model investigated by the authors are based on the realization of the idea of strategic reflection. The most difficult in mathematical terms practical implementation of the author's idea is the use of simulation models, for which the algorithms for finding optimal solutions have known limitations, and the study of sensitivity in most cases is very difficult. On the example of consideration of rather standard discrete- event simulation model the article presents typical methodological techniques that allow ranking variable parameters by sensitivity and, in the future, to expand the scope of definition of variable parameter to which the simulation model is most sensitive. In the transition to the “improved” model, it is also possible to simultaneously exclude parameters from it, the influence of which on the optimized functional is insignificant, and vice versa — the introduction of new parameters corresponding to real processes into the model.

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

The article presents the results of numerical modeling of the vortex spatial non-stationary motion of the medium arising near the lateral and bottom surfaces of the descent module during its movement in the atmosphere of Mars. The numerical study was performed for the high-speed streamline regime at various angles of attack. Mathematical modeling was carried out on the basis of the Navier – Stokes model and the model of equilibrium chemical reactions for the Martian atmosphere gas. The simulation results showed that under the considered conditions of the descent module motion, a non-stationary flow with a pronounced vortex character is realized near its lateral and bottom surfaces. Numerical calculations indicate that, depending on the angle of attack, the nonstationarity and vortex nature of the flow can manifest itself both on the entire lateral and bottom surfaces of the module, and, partially, on their leeward side. For various angles of attack, pictures of the vortex structure of the flow near the surface of the descent vehicle and in its near wake are presented, as well as pictures of the gas-dynamic parameters fields. The non-stationary nature of the flow is confirmed by the presented time dependences of the gas-dynamic parameters of the flow at various points on the module surface. The carried out parametric calculations made it possible to determine the dependence of the aerodynamic characteristics of the descent module on the angle of attack. Mathematical modeling is carried out on the basis of the conservative numerical method of fluxes, which is a finitevolume method based on a finite-difference writing of the conservation laws of additive characteristics of the medium using «upwind» approximations of stream variables. To simulate the complex vortex structure of the flow over descent module, the nonuniform computational grids are used, including up to 30 million finite volumes with exponential thickening to the surface, which made it possible to reveal small-scale vortex formations. Numerical investigations were carried out on the basis of the developed software package based on parallel algorithms of the used numerical method and implemented on modern multiprocessor computer systems. The results of numerical simulation presented in the article were obtained using up to two thousand computing cores of a multiprocessor complex.

В работе развивается иерархический метод математического и компьютерного моделирования интервально-стохастических тепловых процессов в сложных электронных системах различного назначения. Разработанная концепция иерархического структурирования отражает как конструктивную иерархию сложной электронной системы, так и иерархию математических моделей процессов теплообмена. Тепловые процессы, учитывающие разнообразные физические явления в сложных электронных системах, описываются системами стохастических, нестационарных и нелинейных дифференциальных уравнений в частных производных, и в силу этого их компьютерное моделирование наталкивается на значительные вычислительные трудности даже с применением суперкомпьютеров. Иерархический метод позволяет избежать указанных трудностей. Иерархическая структура конструкции электронной системы в общем случае характеризуется пятью уровнями: 1 уровень — активные элементы ЭС (микросхемы, электро-, радиоэлементы); 2 уровень — электронный модуль; 3 уровень — панель, объединяющая множество электронных модулей; 4 уровень — блок панелей; 5 уровень — стойка, установленная в стационарном или подвижном помещении. Иерархия моделей и моделирования стохастических тепловых процессов строится в порядке, обратном иерархической структуре конструкции электронной системы, при этом моделирование интервально-стохастических тепловых процессов осуществляется посредством получения уравнений для статистических мер. Разработанный в статье иерархический метод позволяет учитывать принципиальные особенности тепловых процессов, такие как стохастический характер тепловых, электрических и конструктивных факторов при производстве, сборке и монтаже электронных систем, стохастический разброс условий функционирования и окружающей среды, нелинейные зависимости от температуры факторов теплообмена, нестационарный характер тепловых процессов. Полученные в статье уравнения для статистических мер стохастических тепловых процессов представляют собой систему 14-ти нестационарных нелинейных дифференциальных уравнений первого порядка в обыкновенных производных, решение которых легко реализуется на современных компьютерах существующими численными методами. Рассмотрены результаты применения метода при компьютерном моделировании стохастических тепловых процессов в электронной системе. Иерархический метод применяется на практике при тепловом проектировании реальных электронных систем и создании современных конкурентоспособных устройств.

A hierarchical method of mathematical and computer modeling of interval-stochastic thermal processes in complex electronic systems for various purposes is developed. The developed concept of hierarchical structuring reflects both the constructive hierarchy of a complex electronic system and the hierarchy of mathematical models of heat exchange processes. Thermal processes that take into account various physical phenomena in complex electronic systems are described by systems of stochastic, unsteady, and nonlinear partial differential equations and, therefore, their computer simulation encounters considerable computational difficulties even with the use of supercomputers. The hierarchical method avoids these difficulties. The hierarchical structure of the electronic system design, in general, is characterized by five levels: Level 1 — the active elements of the ES (microcircuits, electro-radio-elements); Level 2 — electronic module; Level 3 — a panel that combines a variety of electronic modules; Level 4 — a block of panels; Level 5 — stand installed in a stationary or mobile room. The hierarchy of models and modeling of stochastic thermal processes is constructed in the reverse order of the hierarchical structure of the electronic system design, while the modeling of interval-stochastic thermal processes is carried out by obtaining equations for statistical measures. The hierarchical method developed in the article allows to take into account the principal features of thermal processes, such as the stochastic nature of thermal, electrical and design factors in the production, assembly and installation of electronic systems, stochastic scatter of operating conditions and the environment, non-linear temperature dependencies of heat exchange factors, unsteady nature of thermal processes. The equations obtained in the article for statistical measures of stochastic thermal processes are a system of 14 non-stationary nonlinear differential equations of the first order in ordinary derivatives, whose solution is easily implemented on modern computers by existing numerical methods. The results of applying the method for computer simulation of stochastic thermal processes in electron systems are considered. The hierarchical method is applied in practice for the thermal design of real electronic systems and the creation of modern competitive devices.

This study investigates the influence of a time-periodic (modulation) magnetic field upon the development of ferroconvection in a densely packed medium saturated with couple stress ferromagnetic fluid. The Darcy model is used to describe the flow in porous medium. The research is important from practical and theoretical point of view. A time-periodic magnetic field is essential in circumscribing channels where the effect of gravity is less or nonexistent to generate circulation. There are numerous engineering uses for this in the manufacturing of magnetic field sensors, charged particle electrode materials, modulators, magnetic resonators, and optical devices. The resulting physical eigenvalue problem is dealt with by using isothermal boundary conditions and the regular perturbation technique with a small time-periodic amplitude. The onset criteria were defined on the supposition that the exchange of stability principle holds. The shift in the thermal Rayleigh number is dependent on the associated parameters: magnetic parameter, Vadasz number, couple stress parameter, porosity, and frequency of the time-periodic function. The results in this case indicate that the onset of ferroconvection can be enhanced or reduced by appropriate changes in the governing parameters.

This study investigates the influence of a time-periodic (modulation) magnetic field upon the development of ferroconvection in a densely packed medium saturated with couple stress ferromagnetic fluid. The Darcy model is used to describe the flow in porous medium. The research is important from practical and theoretical point of view. A time-periodic magnetic field is essential in circumscribing channels where the effect of gravity is less or nonexistent to generate circulation. There are numerous engineering uses for this in the manufacturing of magnetic field sensors, charged particle electrode materials, modulators, magnetic resonators, and optical devices. The resulting physical eigenvalue problem is dealt with by using isothermal boundary conditions and the regular perturbation technique with a small time-periodic amplitude. The onset criteria were defined on the supposition that the exchange of stability principle holds. The shift in the thermal Rayleigh number is dependent on the associated parameters: magnetic parameter, Vadasz number, couple stress parameter, porosity, and frequency of the time-periodic function. The results in this case indicate that the onset of ferroconvection can be enhanced or reduced by appropriate changes in the governing parameters.

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

We investigated numerically the dynamics of kinks of modified sine-Gordon equation in the model with localized spatial modulation of a periodic potential (or impurity). We considered the case of two identical impurities. We showed the possibility of collective effects of the influence of impurities, which are heavily dependent on the distance between them. We demonstrated the existence of a certain critical value of the distance between impurities, which has two qualitatively different scenarios of the dynamic behavior of kink.

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

In this paper we consider software implementation of three-dimensional modeling of thermal processes in multilayer integrated circuits based on low-temperature co-veneering ceramic. The results obtained by the software implemented by the example of the radio frequency receiver module based on low-temperature ceramics for autonomous navigation systems. And also provides a comparison with the results of certified software product.

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

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

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

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

The results of verification carried out for investigations of hydrodynamic effect on reentry conicalsegmental space vehicle are presented in the paper. The program complex Flow Vision is used for this analysis. The purpose of the study is verification of using Flow Vision program complex for problem solving mentioned above on the base of comparison between calculated and experimental data, obtained on the Apollo landing models and new development reentry spacecraft of manned transporting spaceship designed by RSC Energia. The comparison was carried out through the data of pressure values on spacecraft model surfaces during its water landing and inertia center motion parameters.

The results of study show good agreement between experimental and calculated data of force effects on vehicle construction during water landing and its motion parameters in the water medium. Computer simulation sufficiently well reproduces influence of initial velocities & water entry angles variations on water landing process.

Using of computer simulation provides simultaneous acquisition of all data information needed for investigation of water landing peculiarities during construction design, notably, hydrodynamic effects for structural strength calculations, parameters and dynamics of center mass motion and vehicle revolution around center mass for estimation water landing conditions, as well as vehicle stability after landing.

Obtained results confirm suitability of using Flow Vision program complex for water landing vehicle investigations and investigations of influence of different landing regimes through wide initial condition change range, that permits considerably decrease extent of expensive experimental tests and realize landing conditions which are sufficiently complicated for realizing in model physical experiments.