Все выпуски

В работе рассматривается эредитарный осциллятор, который характеризуется осцилляционным уравнением с производными дробных порядков $\beta$ и $\gamma$ в смысле Герасимова–Капуто. С помощью преобразования Лапласа были получены аналитические решения и функция Грина, которые определяются через специальные функции типа Миттаг-Леффлера и обобщенной функции Райта. Доказано, что при фиксированных значениях $\beta = 2$ и $\gamma = 1$ найденное решение переходит в классическое решение для гармонического осциллятора. Согласно полученным решениям были построены расчетные кривые и фазовые траектории эредитарного колебательного процесса. Установлено, что в случае внешнего периодического воздействия на эредитарный осциллятор могут возникать эффекты, присущие классическим нелинейным осцилляторам.

В приближении однородной намагниченности построена математическая модель трехслойной ячейки памяти MRAM c осью анизотропии, расположенной перпендикулярно запоминающему ферромагнитному слою ячейки (перпендикулярная анизотропия). Предполагается, что первоначально намагниченность свободного слоя ячейки ориентирована вдоль оси анизотропии и соответствует состоянию «нуль». Одновременное мгновенное включение спин-поляризованного тока и магнитного поля воздействует на намагниченность свободного слоя и может перевести ее в противоположное положение, соответствующее состоянию «единица». Математическое описание эффекта основано на классическом векторном уравнении Ландау–Лифшица с диссипативным членом в форме Гильберта. В нашей модели учтены взаимодействия намагниченности с внешним магнитным полем и эффективными полями анизотропии и размагничивания, а также с током инжекции в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от управляющих параметров: величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Установлено, что в данной системе, в отличие от системы с продольной анизотропией, дополнительные состояния равновесия отсутствуют. Проведен анализ устойчивости основных состояний равновесия по первому приближению. Построены бифуркационные диаграммы, характеризующие типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно, методом Рунге–Кутты, построены траектории переключения. Найдены комбинации управляющих параметров, при которых переключение невозможно. Найдены области существования устойчивых и неустойчивых предельных циклов системы. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Проведено сравнение значений порогового тока в моделях с продольной и перпендикулярной анизотропией при нулевом магнитном поле и показано, что в модели с перпендикулярной анизотропией ток переключения почти на порядок ниже, чем в модели с продольной анизотропией.

Астероидно-кометная опасность в течение последних десятилетий признана научными и правительственными кругами всех стран мира одной из самых существенных угроз развития и даже существования нашей цивилизации. Одним из аспектов деятельности по предотвращению этой опасности является изучение вторжения достаточно крупных метеорных тел в атмосферу и их движения в ней, сопровождаемых большим числом физическо-химических явлений. Особый интерес вызывает падение метеорных тел, для которых прослежены их траекторные и прочие характеристики, и найдены сами выпавшие метеориты или их фрагменты. В настоящей работе изучено падение именно такого тела. На основе комплексной физико-математической модели, определяющей движение и разрушение космических тел естественного происхождения в атмосфере Земли, рассмотрены движение и фрагментация очень яркого болида Бенешов (Benešov, EN070591), который был зарегистрирован в Чехии Европейской наблюдательной системой в 1991 г. Для этого болида были получены уникальные наблюдательные данные, включая спектры излучения. В настоящей работе проведено моделирование аэробаллистики метеороида Бенешов и его фрагментов с учетом их сложного характера разрушения под воздействием тепловых и силовых факторов. Скорость метеорного тела, унос массы под действием тепловых потоков определяются из решения системы уравнений классической физической теории метеоров. При этом учитывается переменность параметра уноса массы по траектории. Процесс фрагментации метеороида рассматривается в рамках модели последовательного дробления на основе статистической теории прочности, с учетом влияния масштабного фактора на предел прочности объекта. Проведены расчеты совместного обтекания системы тел (осколков метеорита) при проявлении эффекта интерференции. Для расчета обтекания конгломерата осколков метеороида разработан метод моделирования на системе сеток, который позволяет рассматривать фрагменты различных форм, размеров и масс, а также допускает достаточно произвольное их относительное положение в потоке. Из-за неточностей в расчете траектории ученые 23 года не могли найти осколки этого болида. Благодаря современным методикам и более точным расчетам ученые выявили место падения, которое оказалось существенно удаленным от ожидаемого. После этого были обнаружены четыре небольших обломка метеорита. Проведенные расчеты движения и разрушения болида Бенешов показывают, что на процессы его взаимодействия с атмосферой влияет множество факторов: массовые и прочностные характеристики болида, параметры движения, механизмы разрушения, процессы взаимодействия фрагментов, включая эффекты интерференции, и др.

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

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

В работе в целом повторен подход из известной работы Д. Армбрустера, Д. Гукенхеймера, Ф.Холмса, где использован подход, основанный на применении метода Галёркина. Результаты такого анализа расширены и развиты. Использование возможностей современных компьютеров помогло существенно дополнить анализ этой задачи. В частности, найти все решения в четырех- и пятичленных аппроксимациях Галёркина, которые для изучаемой краевой задачи следует интерпретировать как состояния равновесия второго рода. Также дан анализ их устойчивости в смысле определения А. М. Ляпунова.

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

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

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

Статья посвящена численному исследованию распространения возмущений внутри плотной засыпки песка, вызванных последовательным воздействием ударной волны, падающей по нормали к поверхности засыпки из воздуха. Постановка задачи следует натурным опытам А.Т. Ахметова с соавторами. Целью работы является объяснение возможных причин усиления сигнала на датчике давления внутри засыпки, которое наблюдается в опытах при некоторых условиях. Математическая модель основана на одномерной системе уравнений Баера – Нунциато для описания плотных течений двухфазных сред с учетом межгранулярных напряжений в фазе частиц. Вычислительный алгоритм основан на методе Годунова для уравнений Баера – Нунциато.

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

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

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

Модель динамики численности популяций лесных насекомых использована для моделирования взаимодействий «лес–насекомые» и оценки возможных повреждений лесных насаждений насекомыми-вредителями. Согласно этой модели популяция рассматривалась как система автоматической регуляции, в которой входные переменные характеризуют влияние модифицирующих (прежде всего климатических) факторов, а цепи обратной связи описывают влияние регулирующих факторов (паразитов и хищников, внутрипопуляционных взаимодействий). На основе этой модели популяционной динамики предложена методика стресс-тестирования — оценки рисков повреждений и гибели лесных насаждений по отношению к вспышкам массового размножения насекомых. Такой опасный вид лесных вредителей, как сосновая пяденица (Bupalus piniarius L.), рассматривался в качестве объекта анализа; проводились компьютерные эксперименты по оценке рисков возникновения вспышек массового размножения при возможных климатических изменениях на территории Средней Сибири. Модельные эксперименты по- казали, что при достаточно умеренном потеплении (не более 4 °С в летний период) риск воздействия насекомых на лес существенно не возрастает. Однако более сильное потепление на территории Средней Сибири в сочетании с уменьшением количества осадков в летний период может вызвать существенное увеличение частоты вспышек массового размножения основного вредителя сосновых лесов — сосновой пяденицы.

В работе представлены результаты моделирования расчетных случаев приводнения возвращаемого аппарата (ВА) пилотируемого транспортного корабля нового поколения в условиях штиля. Рассмотрены случаи посадки ВА с работающими и с выключенными двигательными установками.

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

Для моделирования переноса межфазных границ применяется метод VOF (Volume of fluid), где перенос контактной поверхности описывается конвективным уравнением, а поверхностное натяжение на межфазной границе учитывается давлением Лапласа. Ключевой особенностью метода является расщепление поверхностных ячеек, куда заносятся данные соответствующей фазы. Уравнения для обеих фаз (уравнения неразрывности, импульса, энергии и другие) в поверхностных ячейках решаются совместно.

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

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