Математическое моделирование. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование
Статья научная
В работе исследуется задача о магнитогидродинамическом течении несжимаемой проводящей полимерной жидкости в плоском канале, помещенном в магнитное поле. По стенкам канала пропущен электрический ток проводимости, а на самих стенках поддерживается разная температура. За основу математической модели магнитной гидродинамики жидких полимеров, рассмотренной в работе, берется обобщенная реологическая модель Покровского - Виноградова с привлеченными к ней уравнениями Максвелла. Для полученной краевой задачи изучаются стационарные решения специального вида, являющиеся аналогами известных вязких течений Пуазейля и Куэтта. Задача для таких решений сводится к краевой задаче для системы нелинейных обыкновенных дифференциальных уравнений второго порядка. Эту задачу мы преобразуем в систему интегральных уравнений, решения которой находим методом простой итерации. Проводится анализ решений для различных параметров задачи, и изучается характер влияния этих параметров на режим течения. Результаты работы демонстрируют возможность контроля за течением жидкого полимера в плоском канале при помощи внешнего магнитного поля и неравномерного нагрева.
Бесплатно
Стохастическая задача Коши в гильбертовом пространстве: модели, примеры, решения
Статья научная
Работа посвящена стохастической задаче Коши для нелинейного уравнения первого порядка со значениями в сепарабельном гильбертовом пространстве с мультипликативным шумом в некотором другом гильбертовом пространстве. В первую очередь в работе представлена модель временной структуры процентной ставки, которая является мерой текущего рынка облигаций. Случайность процесса, описывающего временные структуру цены облигации, обусловлена тем, что экономические показатели изменяются во времени и неизвестны заранее. Рассмотрены методы вычисления форвардной кривой, описывающей временную структуры цены облигации, и переход от них к решению задачи Коши указанного вида. Приведены условия на исходные отображения, необходимые для существования и единственности решения, и построены примеры отображений, удовлетворяющих этим условиям. Рассмотрены слабое и мягкое решения задачи Коши, приведены результаты существования и единственности мягкого решения, показана связь мягкого и слабого решений, из которой следует существование и единственность слабого решения задачи Коши.
Бесплатно
Статья научная
Рассмотрены стохастические модели температурного поля в зоне резания и радиальной составляющей силы резания при шлифовании. Стохастичность взаимодействия инструмента и заготовки вызвана вероятностным строением шлифовального круга, что обусловлено технологией его производства. Рассмотрены вспомогательные модели, необходимые для правильного функционирования основных моделей. Построена информационная схема программного комплекса, включающая основные модули и связи между ними. Показаны масштабы проводимых расчетов и необходимость применения параллельных технологий. Структурно программа разбита на три блока, каждый из данных блоков представляет собой набор вложенных циклов (глубина вложенности от 2 до 4). При этом циклы, начиная с глубины 2, допускают эквивалентное преобразование к виду, содержащему независимые итерации. Цикл верхнего уровня допускает распараллеливание с условием синхронизации входных данных в начале каждой итерации. В параллельной реализации программного комплекса использована комбинация технологий MPI и OpenMP.
Бесплатно
Статья научная
Теория уравнений соболевского типа переживает эпоху своего расцвета. Большое число исследований посвящено детерминированным уравнениям и системам. Однако в натурных экспериментах возникают математические модели, содержащие случайные возмущения, например, в виде белого шума. Поэтому в последнее время все чаще появляются исследования, посвященные стохастическим дифференциальным уравнениям. В данной работе в рамках теории уравнений соболевского типа рассмотрена математическая модель Буссинеска - Лява с аддитивным белым шумом. При изучении модели полезными оказались методы и результаты теории уравнений соболевского типа с относительно p-ограниченными операторами. Поскольку модель представлена вырожденным уравнением математической физики, то к ней трудно применимы существующие ныне подходы Ито - Стратоновича - Скорохода. Мы используем уже хорошо зарекомендовавший себя при решении уравнений соболевского типа метод фазового пространства, заключающийся в редукции сингулярного уравнения к регулярному, определенному на некотором подпространстве исходного пространства. В первой части статьи собраны основные факты теории (L,р)-ограниченных операторов. Во второй -рассмотрена задача Коши для стохастического линейного уравнения соболевского типа высокого порядка. В качестве примера приведена математическая модель Буссинеска - Лява.
Бесплатно
Стохастическое моделирование замкнутых кривых на плоскости
Статья научная
Наиболее универсальный метод имитационного моделирования - стохастическое моделирование. Первоначально Энрико Ферми в 1930-х годах в Италии, а затем Джон фон Нейман и Станислав Улам в 1940-х в Лос-Аламосе предложили использовать стохастический подход для аппроксимации многомерных интегралов в уравнениях переноса, возникших в связи с задачей о движении нейтрона в изотропной среде. После начала использования компьютеров произошeл большой прорыв, и этот метод стал применяться в самых разных задачах, для решения которых стохастический подход оказался более эффективным, чем другие математические методы. В данной работе изучается форма случайного выпуклого овала на плоскости и более общая задача форма случайной замкнутой кривой на плоскости, исследуется изопериметрическое отношение - отношение квадрата длины кривой к площади ограниченной кривой. Величина этого отношения в силу изопериметрического неравенства ограниченна и характеризует отклонение кривой от окружности. Определяется конечномерное многообразие замкнутых регулярных кривых на плоскости и его бесконечномерный аналог. Изучается вероятностные распределения изопериметрического отношения на них. Основной результат состоит в установлении аналитического закона вероятностного распределения отношения - как распределения Фреше являющиеся частным случаем обобщенного распределения экстремальных значений. Основным используемым методом является разложение Фурье опорной функции множества на плоскости и применение математических пакетов Mathematica и Matlab при стохастическом моделировании.
Бесплатно
Точные решения в нелинейной модели теплопередачи
Статья научная
Настоящая работа продолжает большой цикл публикаций авторов, посвященных решениям нелинейного уравнения теплопроводности, которые имеют тип тепловой волны, распространяющейся по нулевому фону с конечной скоростью. Изучается проблема построения точных решений искомого типа для нелинейного уравнения теплопроводности с источником (стоком) и определения их свойств. Особенностью подобных решений является то, что на фронте тепловой волны параболический тип уравнения имеет вырождение, из-за чего оно приобретает необычные для параболических уравнений свойства. Рассмотрены два типа решений: простая волна, которая движется с постоянной скоростью и имеет вид уединенной волны (солитона); волна с экспоненциальным законом движения фронта. В обоих случаях построение редуцируется к задачам Коши для обыкновенных дифференциальных уравнений (ОДУ) второго порядка, которые наследуют особенность от исходной задачи. Построены фазовые портреты ОДУ, установлены свойства траекторий, проходящих через особые точки. Также получены разложения искомых решений в степенные ряды, для которых найдены оценки радиуса сходимости.
Бесплатно
Трехмерная симметричная задача протекания для уравнений Навье - Стокса
Статья научная
В работе Ж. Лерэ (1933) доказана рассуждением от противного разрешимость краевой задачи для уравнений Навье - Стокса при дополнительном условии нулевого потока через каждую связную компоненту границы области течения. При этом же условии Э. Хопф (1941) получил априорную оценку решения. Остается открытым вопрос: имеет ли эта задача решение при выполнении лишь необходимого условия суммарного нулевого потока? Ранее разрешимость трехмерной задачи протекания установлена при малых значениях потоков (Х. Фуджита, 1961; Р. Финн, 1961), либо при условии близости течения к потенциальному (Х. Фуджита и Х. Моримото, 1995). В серии работ М.В. Коробкова, К. Пилецкаса и Р. Руссо (2011 - 2015) положительный ответ на этот вопрос получен для плоских и осесимметричных течений без ограничений на величину потоков. В данной работе задача протекания для уравнений Навье - Стокса рассматривается в трехмерной области типа сферического слоя. Получена априорная оценка решения этой задачи при следующих дополнительных условиях: течение имеет плоскость симметрии; поток через внутреннюю границу области положителен. Из этой оценки вытекает разрешимость указанной задачи.
Бесплатно
Статья научная
Математическое моделирование широко применяется для исследований во всех естественных науках, в отраслях промышленности, в экономике, биологии и других областях. Для решения конкретных задач используются уже существующие или создаются новые модели и численные методы. Наиболее надежным способом проверки качества разностной схемы является сравнение численного решения, где это возможно, с точным решением задачи. В качестве такого эталонного решения построено точное решение задачи о сходящейся ударной волне и о динамическом сжатии газа, находящегося в сферическом сосуде с непроницаемой стенкой. В начальный момент времени наружная граница газа скачком начинает двигаться с отрицательной скоростью, и в газ от границы начинает распространяться ударная волна. Ускорение границы и сферичность определяют движение ударной волны и структуру течения газа между фронтом ударной волны и границей. Изложенная постановка задачи принципиально отличается от ранее известных постановок задачи о схождении автомодельной ударной волны к центру симметрии и ее отражении от центра, в которых отсутствует граница газа.
Бесплатно
Управление в бинарных моделях с разладкой
Статья научная
В статье рассматривается задача динамического управления портфелем для бинарной модели с разладкой. Рассматривается апостериорный подход, то есть в процессе решения задачи производится обнаружение разладки с последующей кластеризацией вершин дерева. На основе этой кластеризации выстраивается алгоритм вычисления оптимального динамического портфеля, применимый для бинарных моделей с разладкой. Используются как симметричный, так и несимметричный штрафы за недостижение поставленной цели управления. Далее анализируется возможность использования бинарной модели в качестве средства приближения непрерывной модели Блэка - Шоулса с разладкой, причем исследуется возможность редукции NP-полной задачи к P-полной задаче с потерей информации.
Бесплатно
Уравнение Хоффа как модель упругой оболочки
Статья научная
Исследуется разрешимость задачи Коши для уравнения Хоффа, моделирующего процесс выпучивания двутавровой балки при постоянной нагрузке и при высоких температурах. Это уравнение относится к классу полулинейных (у оператора действующего на исходную функцию можно выделить линейную часть и нелинейную) уравнений соболевского типа. Разрешимость абстрактных уравнений соболевского типа в банаховых пространствах исследовалась в работах Г.А. Свиридюка и его учеников с помощью метода фазового пространства. Уравнение Хоффа задается на гладком компактном ориентированном римановом многообразии без края. Многобразие в данном случае понимается, как упругая двухсторонняя оболочка. Удается редуцировать исходную задачу к задаче Коши для абстрактного уравнения соболевского типа и применить общую теорию. Редукция основана на теории Свиридюка относительно р-ограниченных операторов и теории Ходжа -Кодаиры о расщеплении пространств дифференциальных форм в прямые суммы подпространств. В результате получена теорема о простоте фазового пространства уравнения Хоффа в случае попадания или нет параметра, характеризующего нагрузку, в спектр оператора Лапласа - Бельтрами.
Бесплатно
Устойчивость стационарного решения неавтономной линеаризованной модели Хоффа на геометрическом графе
Статья научная
Статья посвящена исследованию устойчивости стационарного решения для неавтономной линеаризованной модели Хоффа на геометрическом графе. Такая модель позволяет описывать конструкцию из двутавровых балок, находящуюся под внешним давлением и воздействием высоких температур. Используя условия устойчивости стационарного решения для такой модели, можно описать условия стабильности конструкции, описываемой данной моделью на геометрическом графе. Отметим, что для линеаризованной модели Хоффа нельзя применить метод экспоненциальных дихотомий, так как относительный спектр оператора уравнения может пересекаться с мнимой осью. Поэтому для исследования устойчивости мы будем применять второй метод Ляпунова. Статья кроме введения и списка литературы содержит две части. В первой из них приводятся условия разрешимости неавтономной линеаризованной модели Хоффа на геометрическом графе, а во второй исследуется устойчивость стационарного решения этой модели.
Бесплатно
Устойчивые методы восстановления зашумленных изображений
Статья научная
Рассматривается задача восстановления зашумленных изображений. Для решения используются два регуляризирующих алгоритма, основу которых составляет тихоновская регуляризация с использованием двух различных недифференцируемых стабилизаторов. Для решения задачи негладкой минимизации привлекается проксимальный метод и субградиентный процесс. Приводятся результаты расчетов на суперкомпьютере «Уран».
Бесплатно
Фазовое пространство модифицированного уравнения Буссинеска
Статья научная
В статье доказана однозначная разрешимость задачи Коши для полулинейного уравнения соболевского типа второго порядка. В работе используются идеи и техника, разработанные Свиридюком Г.А. при исследовании задачи Коши для полулинейного уравнения соболевского типа первого порядка, и Замышляевой А.А. при решении задачи Коши для линейного уравнения соболевского типа высокого порядка. В работе так же используется теория дифференцируемых банаховых многообразий, которая окончательно оформилась в работах С. Ленга. В качестве приложения приведена начально-краевая задача для модифицированного уравнения Буссинеска. Рассмотрено два случая - первый, когда оператор L при старшей производной по времени непрерывно обратим, тогда для любой точки из касательного расслоения исходного банахова пространства существует единственное решение, лежащее в этом пространстве как траектория. Особое внимание было уделено второму случаю, когда оператор L не является непрерывно обратимым, тогда уравнение Буссинеска является вырожденным, и было построено для него локальное фазовое пространство. Приводятся условия, при которых фазовое пространство данного уравнения является простым банаховым многообразием.
Бесплатно
Фазовое пространство первой начально-краевой задачи для системы Осколкова высшего порядка
Статья научная
В последние десятилетия теория уравнений соболевского типа активно изучается в различных аспектах. Применение полугруппового подхода к теории уравнений соболевского типа получило плодотворное развитие в работах научного направления, которое возглавляет Г.А. Свиридюк. Данная работа примыкает к этому научному направлению. В работе исследуется первая начально-краевая задача для системы Осколкова. В нашем случае система моделирует динамику плоскопараллельного течения несжимаемой вязкоупругой жидкости Кельвина - Фойгта высшего порядка. Эта задача имеет преимущество, так как фазовое пространство для вышеуказанной системы может быть описано полностью при различных значениях параметра, который характеризует упругие свойства жидкости. Изложению этого факта и посвящена данная статья. Исследование проводится в русле теории полулинейных автономных уравнений соболевского типа на основе понятий квазистационарной траектории и относительно спектрально ограниченного оператора.
Бесплатно
Статья научная
В статье изложена допускающая алгоритмизацию методика приближенного вычисления нормализованных ключевых функций в задаче о ветвлении периодических экстремалей гладкого функционала действия вблизи его точки минимума. Периодические экстремали такого функционала служат прототипами периодических колебаний динамических систем, сегнетоэлектрических фаз кристаллов, нелинейных периодических волн и т.д. Изучение бифуркации циклов динамических систем посредством ключевых уравнений и ключевых функций было недавно проведено в работах А.П. Карповой, Н.А. Копытина, Е.В. Деруновой и Ю.И. Сапронова в случаях двойных резонансов p 1 : p 2 : p 3, p 12 3. В настоящей статье рассмотрен мало изученный случай p 1 = p 2 = p 3 = 1. Предложенная в статье исследовательская схема опирается на вариационную версию метода Ляпунова - Шмидта, в соответствии с которой численное и качественное описание бифуркаций циклов сводится к анализу ветвления критических точек так называемой ключевой функции. В качестве демонстрационной модели рассмотрен функционал действия, соответстваующий обыкновенному дифференциальному уравнению шестого порядка. Приведены примеры раскладов ветвей критических точек и описан подход к классификации таких раскладов, основанный на разбиении бифурцирующих ветвей на подмножества с одинаковыми индексами Морса и на описании взаимных примыканий бифурцирующих критических точек.
Бесплатно
Формирование кольцевых структур в высыхающей под шаблоном пленке коллоидного раствора
Статья научная
В работе предложен метод формирования концентрических кольцевых структур из микро- и наночастиц на непроницаемой поверхности в результате испарения пленки коллоидного раствора под специальным шаблоном. Описанный способ изучается теоретически. Математическая модель начальной стадии процесса решена численно. Модель включает уравнение неразрывности, уравнение движения и уравнение конвекции-диффузии для случая системы с переменной массой. Результаты расчетов показали, что в заданных шаблоном местах массовая доля частиц, движимых радиальным течением, увеличивается значительно. Для подтверждения метода необходимо проведение эксперимента.
Бесплатно
Численное исследование дисперсионных волн, возникающих при движении подводного оползня
Статья научная
Исследуются поверхностные волны, возникающие при сходе подводного оползня по криволинейному склону дна глубокого водоема. Для изучения таких волн используются модели мелкой воды первого и второго приближения. Оползень описывается в рамках модели движения квазидеформируемого тела по криволинейной поверхности под действием внешних сил. Численный алгоритм решения нелинейно-дисперсионных уравнений основан на конечно-разностной аппроксимации системы уравнений гиперболического типа, аналогичной системе уравнений мелкой воды первого гидродинамического приближения и отличающейся от последней лишь правой частью, и уравнения эллиптического типа для осредненной по глубине дисперсионной составляющей давления. Выполнено сопоставление численных результатов, полученных в рамках бездисперсионной модели мелкой воды и нелинейно-дисперсионной модели.
Бесплатно
Численное моделирование внутрикамерных нестационарных турбулентных течений. Часть 1
Статья научная
В работе представлена методика моделирования трехмерных внутренних нестационарных турбулентных течений, в частности, течения продуктов сгорания в твердотопливном ракетном двигателе. Приведена система определяющих уравнений, описывающая поток сжимаемого вязкого газа, записанная в цилиндрической системе координат. Предложен вычислительный алгоритм, разработанный на основе модифицированной схемы расщепления векторов потоков, относящийся к классу методов, использующих подход Годунова. Данный алгоритм пригоден для сквозного расчета внутреннего течения по всему тракту ракетного двигателя, включающего как зоны дозвукового течения в камере, так и зону сверхзвукового течения в сопле. Полученные результаты численного моделирования внутреннего турбулентного течения газа в модельном ракетном двигателе показывают осциллирующий ударно-волновой характер процессов, протекающих в камере двигателя на начальном этапе работы. Определено время выхода двигателя на стационарный режим.
Бесплатно
Численное моделирование внутрикамерных нестационарных турбулентных течений. Часть 2
Статья научная
Течение потока газа в твердотопливном ракетном двигателе определяется особенностями физико-химических процессов, протекающих в камере сгорания, и процесса истечения газа из сопла. В работе предложена методика моделирования внутренних нестационарных турбулентных течений в ракетном двигателе с зарядом твердого топлива телескопического типа. Алгоритм разработан на основе системы уравнений сохранения гидромеханических параметров, описывающих поток сжимаемого вязкого газа. Система уравнений была записана с использованием цилиндрической системы координат. Представленная численная методика относится к классу алгоритмов, использующих подход Годунова. Основой разработанного алгоритма является схема расщепления векторов потоков, модифицированная для вязких течений. Предложенный алгоритм позволяет производить сквозной расчет течения продуктов сгорания твердого топлива по всему тракту ракетного двигателя. Результаты, полученные в ходе численного моделирования потока в ракетном двигателе, позволяют сделать анализ зависимости температуры газа на стенке двигателя от скорости горения низкотемпературной внешней шашки телескопического заряда.
Бесплатно
Численное моделирование решений обратной граничной задачи теплопроводности
Статья научная
В работе рассмотрена обратная граничная задача теплопроводности. Для ее решения предложены различные подходы, основанные на использовании преобразований Лапласа и Фурье. Применение преобразования Лапласа позволило получить операторное уравнение, характеризующее явную зависимость искомой граничной функции от исходных данных на другой границе. Метод, основанный на использовании прямого и обратного преобразований Фурье по переменной, характеризующей время, позволяет получать устойчивые решения, погрешность которых является неулучшаемой по порядку. Предложенные подходы послужили основой для разработки алгоритмов численного решения рассматриваемой задачи и для проведения вычислительного эксперимента, в результате которого были решены обратные задачи для некоторых модельных функций.
Бесплатно
