Все выпуски

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

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

В ограниченной по переменной $z$ области, имеющей слабо горизонтальную неоднородность, исследуется задача определения сверточного ядра $k(t,x)$, $t>0$, $x\in {\Bbb R}$, входящего в гиперболическое интегро-дифференциальное уравнение второго порядка. Предполагается, что это ядро слабо зависит от переменной $x$ и разлагается в степенной ряд по степеням малого параметра $\varepsilon$. Построен метод нахождения первых двух коэффициентов $k_{0}(t)$, $k_{1}(t)$ этого разложения по заданным первым двум моментам по переменной $x$ решения прямой задачи при $z=0$.

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

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

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

В данной работе рассматриваются системы материальных точек в евклидовом пространстве, взаимодействующих как друг с другом, так и с внешним полем. В частности, рассматриваются системы частиц, взаимодействие между которыми описывается однородным потенциалом степени однородности α=-2. Для этих систем существует дополнительная скрытая симметрия, которой соответствует первый интеграл движения, называемый нами интегралом Якоби. Данный интеграл указывался ранее в различных работах, начиная с Якоби, однако мы приводим его в более общем виде.

Исследуется эволюция угла наклона оси вращения планеты в поле притяжения звезды и внешних планет, входящих в планетную систему. Считаем, что исследуемая планета является динамически-симметричным твердым телом $(A = B)$. Полагаем также, что сама планета и внешние планеты движутся по кеплеровским эллипсам вокруг звезды со средними движениями $\omega$ и $\omega_2,\ldots ,\omega_N$, где $N$ - число небесных тел, воздействующих на планету. В переменных Депри-Андуайе получена функция Гамильтона задачи в рамках спутникова приближения. Проведено осреднение функции Гамильтона по быстрым переменным вращательного и орбитального движений при условии отсутствия резонансов между быстрыми частотами указанных движений. Показано, что осредненная функция Гамильтона содержит, помимо классических параметров, параметры $D_i$, являющиеся функционалами на семействе орбит исследуемой планеты и внешних планет. Показано, что осредненная функция Гамильтона допускает разделение переменных и, как следствие, существует три первых интеграла в инволюции. При рассмотрении гравитационных моментов от внешних планет как малых возмущений, получены, с помощью интеграла энергии осредненных уравнений, явные приближенные формулы для угла нутации исследуемой планеты. Получены также приближенные формулы для возмущенного периода прецессии планеты. Проведены расчеты размаха колебаний по углу нутации планеты, возмущенного периода ее прецессии для частного случая планетной системы, состоящей из звезды, самой планеты и массивной внешней планеты (подобной Юпитеру) с симметрично расположенными орбитами, плоскости которых пересекаются под углом $\gamma$.

В статье разработано приближенно-аналитическое решение задачи конформного отображения внутренних точек произвольного многоугольника на единичный круг. На предварительном этапе задача конформного отображения сформулирована в виде краевой задачи (задача Шварца). Последняя сведена к решению интегрального уравнения Фредгольма второго рода с ядром типа Коши относительно неизвестной комплексной функции плотности на границе области с последующим вычислением интеграла Коши. Разработанное приближенно-аналитическое решение основано на разложении ядра Коши в системе многочленов Лежандра первого и второго рода. Выполнена априорная и апостериорная оценки сходимости и точности заданного решения. Определены экспоненциальная сходимость решения в $L_2\left([0,1]\right)$ и полиномиальная в $C\left([0,1]\right)$. Для наглядного сравнения результативности разработанного решения приведены расчеты на тестовых примерах.

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