Все выпуски

С помощью численных методов исследованы пространственно-временная структура и спектральные свойства излучения долгоживущих нелинейных локализованных волн бризерного типа для двух нелинейных дифференциальных уравнений типа Клейна – Гордона — уравнения синус-Гордона и уравнения $\varphi^4$ — при наличии одномерной пространственной модуляции потенциала, играющей роль притягивающей примеси. Рассмотрен случай протяженной примеси прямоугольного профиля, которая выступает в качестве потенциальной ямы для движущегося топологического солитона (кинка) и генератора локализованных колебаний. Для решения начально-краевых задач применялись метод конечных разностей и метод линий. В целях повышения вычислительной производительности и сокращения времени моделирования реализована векторная оптимизация численных алгоритмов на языке Python с использованием библиотеки NumPy, что позволило сохранить требуемую точность расчетов на длительных временных интервалах. Установлено, что возбуждаемые после прохождения примеси локализованные волны сопровождаются постепенным затуханием их амплитуды из-за потери энергии на излучение волн. Временная зависимость этого излучения разделена на зоны сильной и слабой нелинейности. В первой из них излучения прошедшего кинка накладываются на излучение возбужденного бризера. Спектральный анализ Фурье излучаемых волн выявил наличие двух основных частот во всех исследованных точках наблюдения: первая соответствует нижней границе линейного фононного спектра ($\omega \approx 1$ для синус-Гордона и $\omega \approx \sqrt{2}$ для $\varphi^4$), а вторая связана с собственной частотой колебаний бризера. Для уравнения синус-Гордона вторая частота равна утроенной частоте колебаний бризера, тогда как для уравнения $\varphi^4$ — удвоенной частоте. Выявлено, что по мере удаления от примеси вклад частоты фононного спектра становится доминирующим. Показано, что амплитуда излучения медленно убывает со временем, следуя за уменьшением амплитуды самого бризера, а ее абсолютная величина вдали от примеси для случая уравнения синус-Гордона значительно превосходит аналогичный показатель для уравнения $\varphi^4$. Полученные количественные закономерности показывают, что варьированием параметров примеси можно изменять характеристики (спектр частот и амплитуду) излучаемых бризером волн.

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

На основе метода расщепления для нестационарного уравнения Шредингера предложена разностная схема численного решения нестационарной системы двух уравнений Шредингера с оператором спин-орбитального взаимодействия для двухкомпонентной спинорной волновой функции. Выполнено компьютерное моделирование эволюции волновых функций внешних нейтронов с различными проекциями полного момента на межъядерную ось и вероятности их передачи при лобовых столкновениях ядер ¹⁸O и ⁵⁸Ni.

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

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

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

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

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

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

В данной работе рассматриваются формулы Рэлея, полученные из интегральных формул Кирхгофа, которые в дальнейшем могут быть применены для получения миграционных изображений. Актуальность проведенных в работе исследований обусловлена распространенностью применения миграции в интересах сейсмической разведки нефти и газа. Предлагаемый подход позволит существенно повысить качество сейсмической разведки в сложных случаях, таких как вечная мерзлота и шельфовые зоны южных и северных морей. Особенностью работы является использование упругого приближения для описания динамического поведения геологической среды, в отличие от широко распространенного акустического приближения. Сложность применения системы уравнений, описывающей состояние линейно-упругой среды, для получения формул Рэлея и алгоритмов на их основе возникает из-за значительного роста количества вычислений, математической и аналитической сложности итоговых алгоритмов по сравнению со случаем акустической среды. Поэтому в промышленной сейсморазведке в настоящий момент не используют алгоритмы миграции для случая упругих волн, что создает определенные трудности, так как акустическое приближение описывает только продольные сейсмические волны в геологических средах. В данной статье представлены итоговые аналитические выражения, которые можно использовать для разработки программных комплексов, используя описание упругих сейсмических волн (продольных и поперечных), тем самым охватывая весь диапазон сейсмических волн (продольных отраженных PP-волн, продольных отраженных SP-волн, поперечных отраженных PS-волн и поперечных отраженных SS-волн). Также в работе приведены результаты сравнения численных решений, полученных на основе формул Рэлея, с численными решениями, полученными сеточно-характеристическим методом. Ценность такого сравнения обусловлена тем, что метод на основе интегралов Рэлея основан на аналитических выражениях, в то время как сеточно-характеристический метод является методом численного интегрирования решения по расчетной сетке. В проведенном сравнении рассматривались различные типы источников: модель точечного источника, широко используемого в морской и наземной сейсморазведке, и модель плоской волны, которую также иногда применяют в полевых исследованиях.