Математическое моделирование. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование

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

Обобщенная модель курьера с дополнительными ограничениями
Статья научная
Конструируется математическая модель процесса последовательного выбора вариантов перемещений и выполнения комплекса работ, осложненных взаимным влиянием действий на различных временных промежутках и условиями предшествования. Исследуется задача маршрутизации с ограничениями и функциями стоимости, включающими зависимость от списка заданий. Постановка ориентирована на решение инженерных задач, возникающих в атомной энергетике и машиностроении. В первом случае допускаются ограничения, зависящие от списка заданий, не выполненных на текущий момент и касающихся демонтирования излучающих элементов оборудования. Во втором случае возможны ограничения, связанные с обеспечением жесткости листа при резке деталей на станках с числовым программным управлением (ЧПУ); в этом случае возникает зависимость от списка уже выполненных работ. Метод решения, связанный с использованием широко понимаемого динамического программирования, излагается в форме алгоритма на функциональном уровне. При наличии условий предшествования не предусматривается построение всего массива значений функции Беллмана. Для конкретного варианта задачи, связанного с листовой резкой на машинах с ЧПУ, предлагаемый (оптимальный) алгоритм реализован на ПЭВМ; приведены результаты вычислительного эксперимента.
Бесплатно

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

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

Статья научная
Исследуется система уравнений с, вообще говоря, переменными коэффициентами, описывающая функционально-инвариантные решения волнового уравнения в пространстве R 3(t, x, y). Хорошо известно, что для единичной матрицы коэффициентов все функционально-инвариантные решения описываются формулой Соболева. В работе доказано, что если решение рассматриваемой системы имеет максимальный произвол (который понимается в смысле теории совместности переопределенных систем дифференциальных уравнений в частных производных), то коэффициенты волнового уравнения связаны алгебраическим соотношением второго порядка (гиперболическим или эллиптическим) и, кроме того, дифференциальным соотношением второго порядка. На множестве дифференциальных уравнений естественно действует группа преобразований, индуцированных заменами пространственных переменных. Получена полная классификация рассматриваемых систем относительно этой группы. Доказано, что есть ровно три класса эквивалентности. В работе используются классические методы теории Рикье исследования переопределенных систем дифференциальных уравнений в частных производных.
Бесплатно

Обратные задачи определения граничных режимов для некоторых уравнений соболевского типа
Статья научная
В работе изучается разрешимость обратных задач нахождения вместе с решением некоторых уравнений соболевского типа также неизвестных коэффициентов специального вида, определяющих граничные режимы (граничные данные) в первой или соответственно третьей начально-краевых задачах. Наличие в подобных задачах неизвестного коэффициента предполагает, что наряду с краевыми и начальными условиями, характерными для соответствующего класса дифференциальных уравнений, задается также дополнительное условие - условие переопределения. В настоящей работе условие переопределения есть условие интегрального переопределения - условие равенства нулю некоторых интегралов по сечениям цилиндрической области плоскостями t=const. Цель работы - доказательство существования регулярных (имеющих все обобщенные по С.Л. Соболеву производные, входящие в уравнение) решений. Наряду с конкретными результатами приведены некоторые возможные их обобщения.
Бесплатно

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

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

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

Оптимальное управление решениями начально-конечной задачи для уравнения Буссинеска - Лява
Статья научная
В работе исследована задача оптимального управления для уравнения соболевского типа второго порядка с относительно полиномиально ограниченным пучком операторов. Доказана теорема существования и единственности сильного решения начальноконечной задачи для данного уравнения. Получены достаточные, а в случае когда бесконечность является устранимой особой точкой A-резольвенты пучка операторов, и необходимые условия существования и единственности оптимального управления такими решениями. Исследована начально-конечная задача для уравнения Буссинеска - Лява, моделирующего продольные колебания упругого стержня. В работе используются идеи и методы, разработанные Г.А. Свиридюком и его учениками. Доказательство теоремы о существовании и единственности оптимального управления для исследуемой задачи опирается на теорию оптимального управления, развитую в работах Ж.-Л. Лионса.
Бесплатно

Оптимизация полигармонического импульса
Статья научная
В теории и практике создания некоторых технических устройств имеется необходимость оптимизации тригонометрических полиномов. В статье изложено решение задачи оптимизации тригонометрического полинома (полигармонического импульса). f (t) := Σ fk cos(kt) with the asymmetry coefficient k=1 k :=f max / |fmin|, f max : f (t, λ), f min := min t f (t, λ).Вычислены оптимальные значения главных амплитуд. В основу представленного в статье анализа положено понятие «минимального страта Максвелла», под которым подразумевается модмножество многочленов фиксированной степени с максимально возможным количеством минимумов при условии, что все минимумы расположены на одном уровне (значения многочлена во всех точках минимума равны между собой). Многочлен f (t) при выполнении данного условия называется максвелловским. Отправной точкой проведенного исследования послужил экспериментально найденный авторами оптимальный набор значений коэффициентов fk для произвольного п. Позже появилось доказательство единственности оптимального многочлена с максимальным количеством минимумов на отрезке [0, π] и найдена общая формула масквелловского многочлена степени п, связанная с ядром Фейера, для которого коэффициент несимметрии равен п. Возникла естественная гипотеза о том, что ядро Фейера задает оптимальный многочлен. В настоящей статье дано обоснование справедливости этой гипотезы.
Бесплатно

Статья научная
Рассматривается задача маршрутизации перемещений с ограничениями и функциями стоимости, допускающими зависмость от списка заданий. Предполагается, что начальное условие процесса с дискретным временем может выбираться в пределах метрического пространства, удовлетворяющего условию полной ограниченности. По постановке задачи предполагается посещение конечной системы мегаполисов (непустых конечных множеств) с выполнением тех или иных работ, стоимости которых зависят всякий раз от пункта прибытия и пункта отправления. Стоимости перемещений и выполняемых работ агрегируются аддитивно. Для решения используется вариант широко понимаемого динамического программирования, обеспечивающий нахождение ε-оптимального решения при любом значении ε>0.
Бесплатно

Ослабление сферических ударных волн в гетерогенных средах
Статья научная
В данной работе приводится анализ численного исследования взаимодействия сферических ударных волн (УВ) с механическими и тепловыми неоднородностями, расположенными на плоской поверхности. Результаты расчетов показали, что наличие на плоской поверхности гетерогенного слоя приводит к ослаблению падающей УВ внутри этого слоя в зависимости от объемной доли конденсированной фазы. Достаточно большая объемная доля твердой фазы в гетерогенном слое приводит к деформации головного фронта УВ и сильное ее торможение, в результате чего отражения от стенки практически не происходит и ножка Маха не образуется, однако при a2>=5*10-2 - она образуется над поверхностью гетерогенного слоя. Увеличение температуры газа в слое приводит к возрастанию скорости (УВ), в то время как амплитуда ее снижается. Приведенные в работе расчеты показали, что увеличение температуры газа в слое способствует увеличению момента силы, действующей на решетку из твердых частиц. При высоких температурах газа в гетерогенном слое, несмотря на силовое и тепловое взаимодействие газа и частиц, происходит образование предвестника.
Бесплатно

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

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

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

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

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

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

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