Все выпуски

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

The influence of colored noise on the equilibrium regimes of nonlinear dynamical systems is investigated. To study the response of the system to small perturbations, we use an asymptotic approach that develops the stochastic sensitivity function technique. The stochastic sensitivity of equilibrium in a general multidimensional dynamical system is defined by some matrix. For this stochastic sensitivity matrix, we obtain a matrix algebraic equation. An exact solution of this equation is given for an important class of nonlinear oscillators with perturbations in the form of colored noises. This theory is applied to the parametric study of the response of the electronic generator with hard excitation to colored noises with various correlation times. The dependence of the dispersion of random states on the characteristic correlation time is investigated. It is shown that this dependence can be nonmonotonic and have maxima corresponding to the resonances. The paper discusses the probabilistic mechanism of the stochastic generation of large-amplitude oscillations caused by color noise.

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

A mathematical model of competition under conditions of biological invasion, written in the form of a system of nonlinear parabolic equations, is considered. The competition of two closely related species — resident and invader — is studied. The dynamics of populations in a heterogeneous area is determined by local interaction and diffusion. For the invader population, interspecific taxis and directed migration caused by heterogeneity of living conditions are taken into account. In computational experiments, sets of migration parameters corresponding to various invasion scenarios are determined. An analysis of the influence of initial distributions on competitive exclusion and coexistence of species is given.

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

We numerically investigate the orbits of star cluster moving in the plane of the Galaxy in a circular orbit (Bok’s problem). Shuster–Plummer model is used as a model of the potential of the cluster. We examine the impact of initial conditions on the star movements, in particular on the star opportunity to fly out of the cluster. The mass calculation of star orbits for various initial values of energy and momentum with respect to the cluster is carried out. The probabilities for a star to leave a cluster are estimated.

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

The formation process of a spherical shock impulse in gas and its interaction with a protective aqueous foam barrier, accompanied by formation of vortex flows, are numerically investigated. The problem is solved in a two-dimensional axisymmetric formulation using a two-phase model of a gas-liquid mixture, which includes the laws of conservation of mass, momentum and energy of the mixture and an equation for the dynamics of volume content of phases.The numerical implementation is carried out on the basis of the OpenFOAM package using the standard compressibleMultiphaseInterFoam solver, modified in accordance with the conditions of the problem and model representations. The discretization of the system of equations in the chosen solver is carried out by the method of finite volumes using the computational Pimple algorithm. A significant decrease in the intensity of the shock wave in its interaction with the aqueous foam barrier is shown and the causes leading to vortex formation in the gas region are revealed. The reliability of the results obtained is estimated by comparison with solutions of a similar problem by other numerical methods.

В статье рассмотрены основные принципы постановок задач в механике твердого тела при наличии связей (с сухим трением и без). Основное внимание уделено предыстории начальных условий задачи, которая должна быть корректно определена таким образом, чтобы не требовалось введения дополнительных гипотез и допущений, выводящих исследование за рамки динамики твердого тела без ударов. Тогда динамика движения (и/или равновесия) твердых тел может быть описана однозначно и без каких-либо парадоксальных ситуаций (парадоксов Пэнлеве). Эта методика иллюстрируется на трех известных задачах механики: опирание твердого тела на одну точку при наличии сухого трения, движение стержня с ползунами в направляющих с сухим трением, опирание твердого тела на две точки с сухим трением («скамейка»).

We consider basic concepts for setting the problems of motion of a rigid body with constraints (with and without dry friction). The main accent is placed upon the prehistory of initial condition of a problems, which should be formulated in a correct manner which would not require introducing additional hypothesis and assumptions which make one to leave the frames of the rigid body dynamics without impacts. With such correct formulation, the dynamics of motion (or equilibrium) of rigid bodies can be described without occurence of some paradoxic situations (Painlev'e paradoxes). The presented methodology is illustrated by three well-known problems in mechanics: 1) rigid body with a single contact point with a surface in the presence of dry friction, 2) sliding bar in the sliding ways with dry friction, 3) rigid body with two point contact in the presence of dry friction («bench»).

Критически обсуждаются различные способы определения иррегулярных и регулярных сил в звездных системах. Наиболее удовлетворительным кажется определение Эддингтона, согласно которому регулярная сила - это сила притяжения сплошной гравитирующей среды, получающейся «размешиванием» вещества по системе. Интерес представляет также определение регулярной силы как математического ожидания случайной силы. Подчеркивается, что время пересечения τ_c, характерное время действия регулярных сил, определяет темп коллективных процессов в системе. Существенно, что регулярные силы могут приводить и, как правило, приводят к бесстолкновительной стохастизации. В этой связи рассматривается квазиэнтропия, среднее по фазовому пространству значение произвольной выпуклой функции от крупнозернистой функции распределения. Максимум квазиэнтропии для невращающихся систем возможен только при изотропном распределении скоростей. Приводятся найденные Антоновым выражения для ее первой и второй вариаций. Если вторая вариация положительна хотя бы для некоторого изменения плотности, то это означает, что данное состояние системы не является наивероятнейшим. Отсюда следует, что эволюция вдоль последовательности политропных шаров невозможна без поступления в систему дополнительной энергии. Напоминается классификация видов фазового размешивания в бесстолкновительных системах.

Кратко рассматривается проблема столкновительной релаксации в гравитирующих системах. Излагается подход к ее решению с точки зрения теории геодезических потоков с последующим усреднением по ансамблю, что требует знания закона распределения случайной силы. Чтобы избежать обрезания распределения Хольцмарка на малых прицельных расстояниях, использовано распределение случайной силы, найденное Петровской. В этом случае оказывается, что отношение эффективного времени стохастизации к времени пересечения пропорционально N^⅓/(ln N)^½, где N>>1 - число тел в системе. Полученная временная шкала столкновительной эволюции практически совпадает с шкалой, ранее предложенной Генкиным.

Various ways of definition of irregular (random) and regular (smoothed) forces in stellar systems are critically discussed. The most satisfactory is Eddington's one according to which the regular force is an attraction force of a continuous fluid resulting from spreading a stellar mass over a system. Also, a definition of the regular force as a mathematical expectation of a random force is of interest. It is emphasized that the crossing time, τ_c, a time scale of regular forces, characterizes the rate of collective processes in the system, including collisionless relaxation, that (as a rule) occurs in gravitating systems. The quasi-entropy, i.e., a result of averaging of an arbitrary convex function of a coarse-grained distribution function over the phase space, is discussed as a measure of collisionless stochastization. For non-rotating systems the maximum of quasi-entropy can be reached only for isotropic velocity distributions. Formulas for the first and second variations of quasi-entropy, found by Antonov, are given. If there exists the density variation so that the second variation of quasi-entropy is positive, then the present state of the system is not the most probable. It follows from this assertion that evolution along a sequence of polytropic spheres is not possible without some energy input to the system. We recall the classification of forms of the phase mixing in collisionless systems.

The problem of collisional relaxation in gravitating systems is briefly discussed. We state the approach to its analysis on the basis of studying geodesic flows and the ensemble averaging as the next step, which requires the knowledge of distribution of a random force. To avoid truncation of Holtsmark's distribution at small impact parameters the distribution of random force by Petrovskaya was used. In that case the ratio of the effective stochastization time to the crossing time is proportional to N^⅓/(ln N)^½, where N>>1 is the number of stars in the system. This evolutionary time scale is close to the one found earlier by Genkin.

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

Mathematical models were developed and the nonlinear boundary value problem of dynamics thinwalled shells of the arbitrary form under action shock pulse is formulated. Dependence of processes of deformation on speed loading, compressibility of a material, finite deformations and large displacements of a shell middle surface, formation and kinetic of plasticity zones of a material during action of a shock wave are considered. Parameterization of a shell surface is carried out by bi-cubic splines. For the description of nonlinear, time and speed dependents of a shell material behavior with anisotropic hardening the generalized model of microplasticity is developed on the account of viscosity of deformation, hysteresis losses and Baushinger's effect. The solution of boundary value problems on the basis of difference schemes is constructed. Results of modeling of nonlinear wave processes in a assemble shell under action of explosion also are presented.

Рассмотрена динамика вращения твердого тела (ротатора) вокруг неглавной оси Oz, проходящей через его центр масс, с учетом диссипативных моментов: сухого трения M_fr, возникающего в опорах из-за поперечных динамических реакций, и квадратичного по угловой скорости ω аэродинамического сопротивления M_R=-c|ω|ω. Показано, что уравнение динамики и вытекающие из него кинетики вращения тела качественно различны в общем и частном случаях инерционных и диссипативных параметров: осевого момента инерции J_zz, коэффициентов c и α=M_fr/√ε²+ω⁴ (ε - угловое ускорение). В частном случае равенства J_zz=c=α обнаружено отсутствие физически возможного решения для вращения по инерции в рамках динамики абсолютно твердого тела. Парадокс разрешается через нормализующее причинно-следственные связи введение запаздывающих величин ε(t-τ) и ω(t-τ), определяющих в согласии с принципом Даламбера поперечные реакции в опорах оси M_x,y(t-τ) и пару M_fr(t-τ). Последняя же определяла темп потери кинетического момента dK_z(t)/dt в момент времени t. Кинетика вращения при этом имеет импульсивный характер так называемого фрикционно-аэродинамического удара. Также путем численного интегрирования продемонстрирована необычная угловая кинетика φ(t) затухающих колебаний ротатора под действием упругого момента M_e=-κφ, характеризующаяся наличием двух фаз: кратковременного стартового участка, зависящего от начальных условий, затем резко переходящего в фазу почти синусоидальных колебаний с медленно убывающей амплитудой.

The article studies the rotational dynamics of a rigid body (rotator) around the central but non-principal axis Oz passing through its center of mass under the action of dry frictional torque M_fr=α√ε²+ω⁴ caused by inertia forces in the axis's supports and the drag momentum M_R=-c|ω|ω quadratic in angular speed ω. It has been shown that the dynamical equations and the equations of the kinetics of the body's rotation, which follow from the dynamical equations, are qualitatively different in general and particular cases of the inertial and dissipative parameters involved: the axial moment of inertia J_zz and the coefficients c and α=Mfr/√ε²+ω⁴ (where ε is the angular acceleration). It is found that in the particular case of the equality J_zz=c=α a physical feasible solution for the inertial rotation within the dynamics of a perfectly rigid body is absent. The paradox is resolved by the introduction of the lagged angular velocity ω(t-τ) and acceleration ε(t-τ) as factors defining due to D'Alembert principle the supports' transversal reactions M_x,y(t-τ) and hence the value of M_fr(t-τ). The last one determines the loss rate of kinetic momentum, i.e. the dK_z(t)/dt at time t. The rotational kinetics had a type of frictional-aerodynamic impact. Also, by numerical integration, there was shown the unusual angular kinetics φ(t) of the damping oscillations of the rotator under the action of the elastic torque M_e=-κφ. The kinetics was characterized by the presence of two phases: the short starting part strongly depending on initial conditions followed by the phase of almost sine wave oscillations with extremely slow damping.