Все выпуски

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

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

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.

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

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

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

We study the class of first order differential equations in partial derivatives of the Clairaut-type, which are a multidimensional generalization of the ordinary differential Clairaut equation to the case when the unknown function depends on many variables. It is known that the general solution of the Clairaut-type partial differential equation is a family of integral (hyper-) planes. In addition to the general solution, there can be particular solutions, and in some cases a special (singular) solution can be found.

The aim of the paper is to find a singular solution of the Clairaut-type equation in partial derivatives of the first order with a special right-hand side. In the paper, we formulate a criterion for the existence of a special solution of a differential equation of Clairaut type in partial derivatives for the case, when the function of the derivatives is a function of a linear combination of partial derivatives of unknown function. We obtain the singular solution for this type of differential equations with trigonometric functions of a linear combination of $n$-independent variables with arbitrary coefficients. It is shown that the task of finding a special solution is reduced to solving a system of transcendental equations containing initial trigonometric functions. The article describes the procedure for evaluation of a singular solution of Clairaut-type equation; the main idea is to find not partial derivatives of the unknown function, as functions of independent variables, but linear combinations of partial derivatives with some coefficients. This method can be used to find special solutions of Clairaut-type equations, for which this structure is preserved.

The work is organized as follows. The Introduction contains a brief review of some modern results related to the topic of the study of Clairaut-type equations. The Second part is the main one and it includes a formulation of the main task of the work and describes a method of evaluation of singular solutions for the Clairaut-type equations in partial derivatives with a special right-hand side. The main result of the work is to find singular solutions of the Clairaut-type equations containing trigonometric functions. These solutions are given in the main part of the work as an illustrating example for the method described earlier. In Conclusion, we formulate the results of the work and describe future directions of the research.

Представлен итерационный алгоритм, который численно решает нелинейные одномерные несингулярные интегральные уравнения Фредгольма и Вольтерры второго рода типа Урысона. Показано, что метод последовательных приближений Пикара может быть использован при численном решении такого типа уравнений. Сходимость числовой схемы гарантируется теоремами о неподвижной точке. При этом квадратурный алгоритм основан на явной форме встроенного правила Рунге–Кутты пятого порядка с адаптивным контролем размера шага. Возможность контроля локальных ошибок квадратур позволяет создавать очень точные автоматические числовые схемы и значительно уменьшить основной недостаток итераций Пикара, а именно чрезвычайно большое количество вычислений с увеличением глубины рекурсии. Наш алгоритм организован так, что по сравнению с большинством подходов нелинейность интегральных уравнений не вызывает каких-либо дополнительных вычислительных трудностей, его очень просто применять и реализовывать в программе. Наш алгоритм демонстрирует практически важные черты универсальности. Во-первых, следует подчеркнуть, что метод столь же прост в применении к нелинейным, как и к линейным уравнениям типа Фредгольма и Вольтерры. Во-вторых, алгоритм снабжен правилами останова, по которым вычисления могут в значительной степени контролироваться автоматически. Представлен компактный C++-код описанного алгоритма. Реализация нашей программы является самодостаточной: она не требует никаких предварительных вычислений, никаких внешних функций и библиотек и не требует дополнительной памяти. Приведены числовые примеры, показывающие применимость, эффективность, надежность и точность предложенного подхода.

We present the iterative algorithm that solves numerically both Urysohn type Fredholm and Volterra nonlinear one-dimensional nonsingular integral equations of the second kind to a specified, modest user-defined accuracy. The algorithm is based on descending recursive sequence of quadratures. Convergence of numerical scheme is guaranteed by fixed-point theorems. Picard’s method of integrating successive approximations is of great importance for the existence theory of integral equations but surprisingly very little appears on numerical algorithms for its direct implementation in the literature. We show that successive approximations method can be readily employed in numerical solution of integral equations. By that the quadrature algorithm is thoroughly designed. It is based on the explicit form of fifth-order embedded Runge–Kutta rule with adaptive step-size self-control. Since local error estimates may be cheaply obtained, continuous monitoring of the quadrature makes it possible to create very accurate automatic numerical schemes and to reduce considerably the main drawback of Picard iterations namely the extremely large amount of computations with increasing recursion depth. Our algorithm is organized so that as compared to most approaches the nonlinearity of integral equations does not induce any additional computational difficulties, it is very simple to apply and to make a program realization. Our algorithm exhibits some features of universality. First, it should be stressed that the method is as easy to apply to nonlinear as to linear equations of both Fredholm and Volterra kind. Second, the algorithm is equipped by stopping rules by which the calculations may to considerable extent be controlled automatically. A compact C++-code of described algorithm is presented. Our program realization is self-consistent: it demands no preliminary calculations, no external libraries and no additional memory is needed. Numerical examples are provided to show applicability, efficiency, robustness and accuracy of our approach.

Работа посвящена исследованию возможности повышения эффективности решения задач внешней аэродинамики. Изучаются методические аспекты применения локально-адаптивных неструктурированных расчетных сеток и пристеночных функций для численного моделирования турбулентных течений около летательных аппаратов. Интегрируются осредненные по Рейнольдсу уравнения Навье–Стокса, которые замыкаются стандартной моделью турбулентности $k–\varepsilon$. Рассматривается обтекание крылового профиля RAE 2822 турбулентным дозвуковым потоком вязкого сжимаемого газа. Расчеты проводятся в программном ВГД-комплексе FlowVision. Анализируется эффективность применения технологии сглаживания диффузионных потоков и формулы Брэдшоу для турбулентной вязкости в качестве мер, повышающих точность решения аэродинамических задач на локально-адаптивных сетках. Результаты исследования показывают, что использование технологии сглаживания диффузионных потоков приводит к существенному уменьшению расхождений в величине коэффициента лобового сопротивления между результатами расчетов и экспериментальными данными. Кроме того, обеспечивается регуляризация распределения коэффициента поверхностного трения на криволинейной поверхности профиля. Эти результаты позволяют сделать вывод о том, что данная технология является эффективным способом повышения точности расчетов на локально-адаптивных сетках. Формула Брэдшоу для динамического коэффициента турбулентной вязкости традиционно используется в модели SST $k–\omega$. В настоящей работе исследуется возможность ее применения в стандартной $k–\varepsilon$-модели турбулентности. Результаты расчетов показывают, что, с одной стороны, данная формула обеспечивает хорошее согласование суммарных аэродинамических характеристик и распределения коэффициента давления по поверхности профиля с экспериментом. Помимо этого, она значительно повышает точность моделирования течения в пограничном слое и в следе. С другой стороны, использование формулы Брэдшоу при моделировании обтекания профиля RAE 2822 приводит к занижению коэффициента поверхностного трения. Поэтому в работе делается вывод о том, что практическое применение формулы Брэдшоу требует ее предварительной валидации и калибровки на надежных экспериментальных данных для рассматриваемого класса задач. Результаты работы в целом показывают, что при использовании рассмотренных технологий численное решение задач внешнего обтекания на локально-адаптивных сетках с применением пристеночных функций обеспечивает точность, приемлемую для оперативной оценки аэродинамических характеристик, а ПК FlowVision является эффективным инструментом решения задач предварительного аэродинамического проектирования, концептуального проектирования и оптимизации аэродинамических форм.

The work is dedicated to investigation of possibility to increase the efficiency of solving external aerodynamic problems. Methodical questions of using locally-adaptive grids and wall functions for numerical simulation of turbulent flows past flying vehicles are studied. Reynolds-averaged Navier–Stokes equations are integrated. The equations are closed by standard $k–\varepsilon$ turbulence model. Subsonic turbulent flow of perfect compressible viscous gas past airfoil RAE 2822 is considered. Calculations are performed in CFD software FlowVision. The efficiency of using the technology of smoothing diffusion fluxes and the Bradshaw formula for turbulent viscosity is analyzed. These techniques are regarded as means of increasing the accuracy of solving aerodynamic problems on locally-adaptive grids. The obtained results show that using the technology of smoothing diffusion fluxes essentially decreases the discrepancy between computed and experimental values of the drag coefficient. In addition, the distribution of the skin friction coefficient over the curvilinear surface of the airfoil becomes more regular. These results indicate that the given technology is an effective way to increase the accuracy of calculations on locally-adaptive grids. The Bradshaw formula for the dynamic coefficient of turbulent viscosity is traditionally used in the SST $k–\omega$ turbulence model. The possibility to implement it in the standard $k–\varepsilon$ turbulence model is investigated in the present article. The calculations show that this formula provides good agreement of integral aerodynamic characteristics and the distribution of the pressure coefficient over the airfoil surface with experimental data. Besides that, it essentially augments the accuracy of simulation of the flow in the boundary layer and in the wake. On the other hand, using the Bradshaw formula in the simulation of the air flow past airfoil RAE 2822 leads to under-prediction of the skin friction coefficient. For this reason, the conclusion is made that practical use of the Bradshaw formula requires its preliminary validation and calibration on reliable experimental data available for the considered flows. The results of the work as a whole show that using the technologies discussed in numerical solution of external aerodynamic problems on locally-adaptive grids together with wall functions provides the computational accuracy acceptable for quick assessment of the aerodynamic characteristics of a flying vehicle. So, one can deduce that the FlowVision software is an effective tool for preliminary design studies, for conceptual design, and for aerodynamic shape optimization.

Представлена гиперболическая модель вязкого теплопроводного газа, в которой для гиперболизации уравнений использован подход Максвелла–Каттанео, обеспечивающий распространение волн с конечными скоростями. В модифицированной модели вместо оригинальных законов Стокса и Фурье использовались их релаксационные аналоги и показано, что при стремлении времен релаксации $\tau_\sigma^{}$ и $\tau_w^{}$ к нулю гиперболизированные уравнения приводятся к классической системе Навье–Стокса негиперболического типа с бесконечными скоростями перемещения вязких и тепловых волн. Отмечено, что рассматриваемая в работе гиперболизированная система уравнений движения вязкого теплопроводного газа инвариантна не только по отношению к преобразованиям Галилея, но и к повороту, поскольку при дифференцировании по времени компонентов тензора вязких напряжений использована производная Яуманна. Для интегрирования уравнений модели применены гибридный метод Годунова (ГМГ) и многомерный узловой метод характеристик. ГМГ предназначен для интегрирования гиперболических систем, в которых имеются как уравнения, записанные в дивергентном виде, так и уравнения, не приводящиеся к таковому (оригинальный метод Годунова применяется только для систем уравнений, представленных в дивергентной форме). При вычислении потоковых переменных на гранях смежных ячеек использован линеаризованный римановский решатель. Для дивергентных уравнений применена конечно-объемная, а для недивергентных — конечноразностная аппроксимация. Для расчета ряда задач в работе также использовался неконсервативный многомерный узловой метод характеристик, который базируется на расщеплении исходной системы уравнений на ряд одномерных подсистем, для решения которых использован одномерный узловой метод характеристик. С помощью описанных численных методов решен ряд модельных одномерных задач о распаде произвольного разрыва, а также рассчитано двумерное течение вязкого газа при взаимодействии ударного скачка с прямоугольной ступенькой, непроницаемой для газа.

A hyperbolic model of a viscous heat-conducting gas is presented, in which the Maxwell – Cattaneo approach is used to hyperbolize the equations, which provides finite wave propagation velocities. In the modified model, instead of the original Stokes and Fourier laws, their relaxation analogues were used and it is shown that when the relaxation times $\tau_\sigma^{}$ и $\tau_w^{}$ tend to The hyperbolized equations are reduced to zero to the classical Navier – Stokes system of non-hyperbolic type with infinite velocities of viscous and heat waves. It is noted that the hyperbolized system of equations of motion of a viscous heat-conducting gas considered in this paper is invariant not only with respect to the Galilean transformations, but also with respect to rotation, since the Yaumann derivative is used when differentiating the components of the viscous stress tensor in time. To integrate the equations of the model, the hybrid Godunov method (HGM) and the multidimensional nodal method of characteristics were used. The HGM is intended for the integration of hyperbolic systems in which there are equations written both in divergent form and not resulting in such (the original Godunov method is used only for systems of equations presented in divergent form). A linearized solver’s Riemann is used to calculate flow variables on the faces of adjacent cells. For divergent equations, a finitevolume approximation is applied, and for non-divergent equations, a finite-difference approximation is applied. To calculate a number of problems, we also used a non-conservative multidimensional nodal method of characteristics, which is based on splitting the original system of equations into a number of one-dimensional subsystems, for solving which a one-dimensional nodal method of characteristics was used. Using the described numerical methods, a number of one-dimensional problems on the decay of an arbitrary rupture are solved, and a two-dimensional flow of a viscous gas is calculated when a shock jump interacts with a rectangular step that is impermeable to gas.

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

A compact finite difference scheme has been developed for modeling convection in a porous medium saturated with a fluid. We consider the problem for a rectangular domain with anisotropic permeability and thermal conductivity properties in terms of stream function and temperature deviation, taking into account Darcy's law. Boundary conditions of impenetrability and a linear distribution of temperature are set. This model is cosymmetric when certain conditions are imposed on the permeability and thermal conductivities. One parametric family of stationary convection regimes arises when mechanical equilibrium loses stability. A numerical method with a fourth-order finite difference approximation for spatial variables and a Runge – Kutta integrator for time has been developed. It has been proved that this scheme preserves cosymmetry. Numerical results for evaluating the critical Rayleigh number have been presented. We compare them with results obtained using a second-order finite-difference method. We show that critical Rayleigh numbers are repeated twice with very high accuracy, which proves cosymmetry preservation. Numerical evaluation of convective regimes and spectral properties are presented. The efficiency of the developed compact finite difference scheme on a nine-point stencil is assessed.

This paper is devoted to the problem of edge criticality identification and ranking in complex networks, which is a part of a modern research direction in the novel network science. The diffusion importance belongs to the set of acknowledged methods that help to identify the significant connections in the graph that are critical to retaining structural integrity. In the present work, we develop the Iterative Diffusion Importance algorithm that is based on the re-estimation of critical topological features at each step of the graph deconstruction. The Iterative Diffusion Importance has been compared with methods such as diffusion importance and degree product, which are two very well-known benchmark algorithms. As for benchmark networks, we tested the Iterative Diffusion Importance on three standard networks, such as Zachary’s Karate Club, the American Football Network, and the Dolphins Network, which are often used for algorithm efficiency evaluation and are different in size and density. Also, we proposed a new benchmark network representing the airplane communication between Japan and the US. The numerical experiment on finding the ranking of critical edges and the following network decomposition demonstrated that the proposed Iterative Diffusion Importance exceeds the conventional diffusion importance by the efficiency for 2–35% depending on the network complexity, the number of nodes, and the number of edges. The only drawback of the Iterative Diffusion Importance is an increase in computation complexity and hencely in the runtime, but this drawback can be easily compensated for by the preliminary planning of the network deconstruction or protection and by reducing the re-evaluation frequency of the iterative process.

This paper is devoted to the problem of edge criticality identification and ranking in complex networks, which is a part of a modern research direction in the novel network science. The diffusion importance belongs to the set of acknowledged methods that help to identify the significant connections in the graph that are critical to retaining structural integrity. In the present work, we develop the Iterative Diffusion Importance algorithm that is based on the re-estimation of critical topological features at each step of the graph deconstruction. The Iterative Diffusion Importance has been compared with methods such as diffusion importance and degree product, which are two very well-known benchmark algorithms. As for benchmark networks, we tested the Iterative Diffusion Importance on three standard networks, such as Zachary’s Karate Club, the American Football Network, and the Dolphins Network, which are often used for algorithm efficiency evaluation and are different in size and density. Also, we proposed a new benchmark network representing the airplane communication between Japan and the US. The numerical experiment on finding the ranking of critical edges and the following network decomposition demonstrated that the proposed Iterative Diffusion Importance exceeds the conventional diffusion importance by the efficiency for 2–35% depending on the network complexity, the number of nodes, and the number of edges. The only drawback of the Iterative Diffusion Importance is an increase in computation complexity and hencely in the runtime, but this drawback can be easily compensated for by the preliminary planning of the network deconstruction or protection and by reducing the re-evaluation frequency of the iterative process.

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

The dual task of integral geometry – to define for a given averaging operator the function class where inversion of that operator is possible – is solved. Those classes are defined ambiguously. Full description of those classes is given in the form of minimal complimentary information necessary to know about the function. The possible to give a constructive description of the class is researched and in the case of a finite averaging system the inversion formulas are given.

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

An approach to the decrease of norm of the correction in Newton’s methods of optimization, based on the Cholesky’s factorization is presented, which is based on the integration with the technique of the choice of leading element of algorithm of linear programming as a method of solving the system of equations. We investigate the issues of increasing of the numerical stability of the Cholesky’s decomposition and the Gauss’ method of exception.

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

The different types of the reduced equations are known in the dynamics a heavy rigid body with a fixed point. Since the Euler−Poisson’s equations admit the three first integrals, then for the first approach the obtaining new forms of equations are usually based on these integrals. The system of six scalar equations can be transformed to a third-order system with them. However, in indicated approach the reduced system will have a feature as in the form of radical expressions a relatively the components of the angular velocity vector. This fact prevents the effective the effective application of numerical and asymptotic methods of solutions research. In the second approach the different types of variables in a problem are used: Euler’s angles, Hamilton’s variables and other variables. In this approach the Euler−Poisson’s equations are reduced to either the system of second-order differential equations, or the system for which the special methods are effective. In the article the method of finding the reduced system based on the introduction of an auxiliary variable is applied. This variable characterizes the mixed product of the angular momentum vector, the vector of vertical and the unit vector barycentric axis of the body. The system of four differential equations, two of which are linear differential equations was obtained. This system has no analog and does not contain the features that allows to apply to it the analytical and numerical methods. Received form of equations is applied for the analysis of a special class of solutions in the case when the center of mass of the body belongs to the barycentric axis. The variant in which the sum of the squares of the two components of the angular momentum vector with respect to not barycentric axes is constant. It is proved that this variant exists only in the Steklov’s solution. The obtained form of Euler−Poisson’s equations can be used to the investigation of the conditions of existence of other classes of solutions. Certain perspectives obtained equations consists a record of all solutions for which the center of mass is on barycentric axis in the variables of this article. It allows to carry out a classification solutions of Euler−Poisson’s equations depending on the order of invariant relations. Since the equations system specified in the article has no singularities, it can be considered in computer modeling using numerical methods.