Математическое моделирование. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование
Модификация метода крупных частиц для исследования течений газовзвесей
Статья научная
В данной работе приводится модификация метода крупных частиц в приложении к исследованиям течений газовзвесей. Показано, что предложенная модификация метода крупных частиц позволяет проводить расчеты поведения ударных волн в газовзвесях без введения в явном виде искусственной вязкости. Это позволило избежать искажения физической картины течения газовзвеси, связанной с наличием осцилляций, имеющих место при распространении ударных волн в неоднородных средах. В данной работе было установлено, что для проведения расчетов распространения ударных волн в газовзвесях с большими числами Куранта может быть использована явная модификация метода крупных частиц. Это позволило значительно сократить время расчета задачи и избежать проведения сложных итерационных процедур, присущих неявным разностным схемам. Было показано, что предложенная в данной работе модификация метода крупных частиц является эффективной и позволяет проводить расчеты даже сильных ударных волн в газовзвесях.
Бесплатно
Модификация метода крупных частиц для решения задач распространения ударных волн и волн разрежения
Статья научная
В данной работе предложена модификация метода крупных частиц. Проведен численный анализ различных модификаций метода крупных частиц применительно к задачам волновой динамики (газовой динамики). Решены задачи расчета распада произвольного разрыва, а также задачи о распространении стационарных ударных волн с отражением от жесткой стенки. Было показано, что предложенная модификация метода крупных частиц наилучшим образом совпадает с аналитическим решением задачи об отражении плоской ударной волны от жесткой стенки. Проведенный численный анализ показал, что данная модификация позволяет проводить устойчивые расчеты течений с большими градиентами изменения параметров. Значительным достоинством предложенной модификации является тот факт, что рассмотренные в работе задачи могут быть решены без введения в законы сохранения "искусственной" вязкости, а также при больших числах Куранта.
Бесплатно
Статья научная
С помощью модифицированного метода Галеркина доказывается однозначная регулярная разрешимость краевой задачи для уравнения смешанного типа второго порядка с произвольным многообразием изменения типа. Теория таких уравнений исходит из прикладных задач, в частности, трансзвуковой газовой динамики. Исследование проведено для случаев, когда вблизи нижнего основания цилиндрической области уравнение имеет эллиптический тип, а вблизи верхнего основания цилиндра - гиперболический или эллиптический тип. В последнем случае разрешимость этой краевой задачи ранее была изучена авторами с помощью другой методики, там впервые была сформулирована ее постановка. Кроме того, в настоящей работе получена оценка погрешности приближенных решений краевой задачи относительно точного решения через параметр регуляризации и собственные значения спектральной задачи Дирихле для уравнения Лапласа.
Бесплатно
Молекулярно-динамическое исследование формирования границы раздела между металлическими нанопленками
Статья научная
В настоящее время изучение совокупности механизмов, ответственных за формирование наноструктур на подложке при осаждении металлов из пучка или из газовой фазы, не завершено. Для выявления этих механизмов необходимы методы исследования процессов на атомарном уровне, среди которых одним из самых мощных является метод молекулярной динамики, что обусловливает его актуальность при решении данной задачи. В работе исследуется влияние ориентации и температуры подложки на формирование границы раздела Ag/Cu.
Бесплатно
Нахождение значений первых собственных функций возмущенных дискретных операторов с простым спектром
Статья научная
В работе получены аналитические формулы для нахождения перых «взвешенных» поправок теории возмущений возмущенных самосопряженных операторов в случае, когда собственные значения невозмущенных операторов простые. Получены оценки остатков сумм функциональных рядов Рэлея-Шредингера. Разработан метод нахождения значений собственных функций возмущенных дискретных операторов с простым спектром.
Бесплатно
Статья научная
В работе получены аналитические формулы для вычисления первых четырех поправок теории возмущений дискретных полуограниченных снизу операторов, когда собственные значения невозмущенных операторов имеют произвольную кратность.
Бесплатно
Начально-конечная задача для неоднородного уравнения Буссинеска - Лява
Статья научная
Рассматривается начально-конечная задача для неоднородного уравнения Буссинеска - Лява. Проводится редукция к абстрактной начально-конечной задаче для уравнения соболевского типа второго порядка. Получены достаточные условия для однозначной разрешимости исходной и абстрактной задач.
Бесплатно
Статья научная
Объектом исследования в работе являются нелинейные обратные коэффициентные задачи для нестационарных дифференциальных уравнений высокого порядка типа псевдогиперболических. Более точно, изучаются задачи определения вместе с решением соответствующего уравнения также неизвестного коэффициента при решении или же при производной решения по временной переменной. Отличительной особенностью рассматриваемых задач является то, что неизвестный коэффициент является функцией лишь от времени. В качестве дополнительного условия в работе используется условие интегрального переопределения. Доказываются теоремы существования регулярных (имеющих все обобщенные по С.Л. Соболеву производные, входящие в уравнение) решений. Техника доказательства основана на переходе от исходной обратной задачи к новой, уже прямой, задаче для вспомогательного интегро-дифференциального уравнения, доказательстве ее разрешимости и построении по решению вспомогательной задачи решения исходной обратной задачи.
Бесплатно
Нелинейный метод проекционной регуляризации
Статья научная
В статье рассмотрен метод проекционной регуляризации, в котором параметр регуляризации выбран из принципа невязки. Получена оценка погрешности этого метода на классе корректности Мг.
Бесплатно
Новые гиперболические модели многокомпонентных гетерогенных сред
Статья научная
Разработка математически корректных и физически непротиворечивых моделей много-фазных сред является актуальной задачей, поскольку не все существующие к настоящему времени модели гетерогенных сред являются таковыми. В данной работе для многокомпонентной среды предлагаются две новые модели - в одно- и многоскоростном приближениях. Модели основаны на законах сохранения. Учитываются вязкие и теплопроводящие свойства смеси. Для приведенных моделей строятся автомодельные решения типа бегущей волны. На примере бинарной смеси расчеты, произведенные в одно- и многоскоростном приближениях. Показывается, что при использовании релаксационных законов для диссипативных процессов системы уравнений относятся к гиперболическому типу.
Бесплатно
О декомпозиции разностных схем при численном решении дифференциально-алгебраических уравнений
Статья научная
Рассматриваются квазилинейные системы обыкновенных дифференциальных уравнений (ОДУ), с тождественно вырожденной матрицей перед производной искомой вектор-функции и разностные схемы, применяемые для их решения. В работе обсуждаются условия, обеспечивающие на каждом шаге вычислительного процесса возможность последовательного решения алгебраических (конечных) уравнений и подстановки этих решений в уравнения динамики. Приведены результаты численных экспериментов для систем ОДУ, описывающих прямоточную котельную установку.
Бесплатно
О дискретизации линейных дифференциальных уравнений
Статья научная
Рассмотрены некоторые вопросы получения дискретного описания дифференциальной системы (ДС) на равномерной сетке. Рассматриваются ДС в виде системы n линейных обыкновенных дифференциальных уравнений первого порядка с постоянными коэффициентами или одно уравнение n-го порядка для наблюдаемого функционала состояния ДС. Изучаемые вопросы дискретизации важны для задач вариационной идентификации и аппроксимации динамических процессов моделями этого типа в конечном интервале. Дано сравнение аналитического равномерного (на основе теоремы Гамильтона-Кэли) и локальных методов дискретизации: на основе разделенных разностей и с помощью интерполяции выборок из n +1 отсчетов многочленами Тейлора степени n. Получена общая формула локальной дискретизации, прозволяющая сравнивать ее разностный и интерполяционные методы. Показано с использованием свойств обратных матриц Вандермонда, что в полученной общей формуле локальной дискретизации ее интерполяционному методу соответствуют (n + 1)-матрицы Тейлора (из коэффициентов многочленов Тейлора), а разностному — (n + 1)-матрицы Паскаля (из чисел треугольников Паскаля). Показано, что невырожденность матрицы наблюдаемости ДС на сетке есть необходимое и достаточное условие как для аналитической дискретизируемости, так и для приведения дискретной системы (описания ДС сетке) к каноническому фробениусовскому виду. Он эквивалентен одному обыкновенному разностному уравнению для наблюдаемой переменной с постоянными коэффициентами. Это уравнение есть основа известного вариационного метода идентификации. Показано, что интерполяционный метод локальной дискретизации есть первое (линейное) приближение формулы равномерной аналитической дискретизации. Показано, что нулевое приближение ее не зависит от коэфффициентов ДС и есть вектор коэффициентов n-й разности. Показано также, что нулевое приближение матрицы наблюдаемости ДС н и матрицы наблюдамости полиномиальной системы y (n) =0 на сетке есть n-матрица Тейлора.
Бесплатно
О задаче минимальной реализации
Статья научная
Предполагается, что для линейной конечномерной стационарной динамической системы с дискретным временем известна степень МакМиллана и конечная последовательность ее марковскиx параметров , . Рассматриваются задачи восстановления по этим данным переходной матрицы-функции системы, минимальных индексов и взаимно простых дробных факторизаций , минимальных решений соответствующих уравнений Безу, минимальной реализации . Для каждой из них существует отдельный алгоритм решения. В данной работе предлагается единый подход к исследованию этих проблем. Он основан на методе индексов и существенных многочленов конечной последовательности матриц. Этот метод был ранее разработан для явного решения задачи факторизации Винера - Хопфа мероморфных матриц-функций. Показано, что решение всех вышеуказанных задач может быть получено, как только будут найдены индексы и существенные многочлены последовательности . Вычисление индексов и существенных многочленов можно осуществить средствами линейной алгебры. Для матриц с элементами из поля рациональных чисел алгоритм реализован в среде Maple в виде процедуры ExactEssPoly.
Бесплатно
О качественном анализе семейства дифференциальных уравнений с первыми интегралами выше 2-й степени
Статья научная
Исследуется семейство дифференциальных уравнений, возникшее в результате обобщения классических интегрируемых случаев динамики твердого тела. Исследуемая система допускает полиномиальные первые интегралы 4 и 6 степени. При определенных ограничениях на параметры семейства дифференциальные уравнения интерпретируются как уравнения движения твердого тела в центральном поле сил, идеальной жидкости, электрически заряженного тела. Проводится качественный анализ уравнений: находятся особые инвариантные множества различной размерности и исследуется их устойчивость по Ляпунову. Для анализа задачи используются обобщения метода Рауса - Ляпунова и программные средства компьютерной алгебры.
Бесплатно
О корректной разрешимости задачи Коши для обобщенного телеграфного уравнения
Статья научная
В работе устанавливается равномерно корректная разрешимость задачи Коши для обобщенного телеграфного уравнения с переменными коэффициентами, частным случаем которого является классическое телеграфное уравнение. Установление корректной разрешимости математических задач является одним из основных условий при их численной реализации. Как известно, для классического телеграфного уравнения решение задачи Коши находится в классе дважды непрервно дифференцируемой функции и с помощью метода Римана выписывается в явном виде. Однако, при этом вопрос устойчивости решения в зависимости от начальных данных, требующий использования соответствующих метрических пространств в этих работах не обсуждается. Между тем этот вопрос является наиболее важным при корректной численной реализации решения задачи, когда его существование и единственность доказаны. В настоящей заметке методами теории полугрупп линейных преобразований, устанавливается равномерно корректная разрешимость задачи Коши в пространствах функций интегрируемых с экспоненциальным весом для некоторого класса дифференциальных уравнений с переменными коэффициентами. Получено точное решение задачи Коши и указаны условия на коэффициенты, при которых задача раномерно корректна в некоторых функциональных пространствах. Следствием из этих результатов является равномерная корректность задачи Коши для классического телеграфного уравнения с постоянными коэффициентами.
Бесплатно
О корректной разрешимости некоторых задач фильтрации в пористой среде
Статья научная
В работе методом теории полугрупп линейных преобразований устанавливается равномерно корректная разрешимость начально-краевых задач для одного класса интегрально-дифференциальных уравнений, рассматриваемых в ограниченной и полуограниченной областях, которые описывают процессы нестационарной фильтрации сжимающей жидкости в пористой среде. Частный случай таких уравнений на полубесконечной прямой с условием Дирихле на границе рассматривался в работе Ю.И. Бабенко. В этой работе требовалось найти градиент давления на границе области. Здесь ответ получен формальным применением дробного интегро-дифференцирования, не затрагивая вопроса о корректной разрешимости и устойчивости решения к погрешностям по исходным данным. При этом решение задачи представляется в виде формального ряда с неограниченным оператором, сходимость которого также не обсуждается. Метод теории сильно непрерывных полугрупп преобразований позволяет установить равномерно корректную разрешимость задач Дирихле и Неймана как для конечных так и бесконечных областей. Это дает возможность в случае задачи Дирихле корректно вычислить градиент давления на границе и значение решения на границе в случае условий Неймана. Здесь же доказана устойчивость решения по начальным данным.
Бесплатно
О криптоанализе системы BBCRS на двоичных кодах Рида - Маллера
Статья научная
В работе рассматривается система система BBCRS - модификация криптосистемы Мак-Элиса, предложенная М. Балди и др. В модификации матрица публичного ключа представляет собой произведение трех матриц: невырожденной -матрицы , порождающей матрицы секретного -кода и невырожденной -матрицы специального вида. Отличие системы BBCRS от системы, предложенной Р. Мак-Элисом, состоит в том, что подстановочная матрица, используемая в системе Мак-Элиса, заменена матрицей , представляющей сумму подстановочной матрицы и матрицы малого ранга . Позже В. Готье и др. построили атаку, позволяющую дешифровать сообщения в случае, когда - обобщенный код Рида - Соломона (ОРС-код) и . Ключевыми этапами построенной атаки являются, во-первых, нахождение пересечения линейных оболочек и , натянутых соответственно на строки матриц и , а во-вторых, нахождение кода по подкоду . В настоящей работе строится атака в случае, когда - двоичный код Рида - Маллера порядка и длины при . В построенной в настоящей работе атаке этапы нахождения кодов и полностью отличаются от соответствующих этапов для ОРС-кодов, а остальные шаги атаки адаптируют известные результаты криптоанализа системы BBCRS на ОРС-кодах.
Бесплатно
О моделировании с использованием дифференциально-алгебраических уравнений в частных производных
Статья научная
Рассматриваются эволюционные системы дифференциальных уравнений в частных производных, зависящие от одной пространственной переменной. Предполагается, что матрицы перед производными искомой вектор-функции вырожденные во всей области определения. Такие системы принято называть дифференциально-алгебраическими уравнениями (ДАУ) в частных производных. Свойства ДАУ существенно отличаются от свойств невырожденных систем. В частности, невозможно судить о типе систем по виду корней характеристических уравнений. В работе вводится понятие расщепляемых систем. Под такими уравнениями понимаются системы, допускающие существование невырожденных преобразований, расщепляющих исходный объект на подсистемы с единственным решением, функциональным произволом от одной из переменных и собственно невырожденную подсистему уравнений в частных производных. Этот прием позволяет исследовать структуру общих решений ДАУ и в ряде случаев установить разрешимость начально краевых задач.
Бесплатно
О модельных движениях в задаче управления при функциональных ограничениях на помеху
Статья научная
Рассматривается задача управления системой, описываемой обыкновенным дифференциальным уравнением. Предполагается, что значения управления и помехи в каждый момент времени содержатся в некоторых компактных множествах. Предполагается также, что помехи удовлетворяют некоторым дополнительным ограничениям функционального характера, отражающим природу рассматриваемой задачи. Качество управления оценивается функционалом, заданым на множестве фазовых траекторий рассматриваемой системы, и непрерывным в метрике равномерной сходимости. Ранее установлено, что стратегия с полной памятью разрешает данную задачу управления при компактных ограничениях на помеху и при других функциональных ограничениях, которые к ним сводятся. Вместе с тем, построенные для этих случаев стратегии не являлись универсальными, то есть они зависели от начальной позиции движения системы. Также оставался открытым вопрос о возможности разрешения задач управления с функциональными ограничениями в более узком (классическом) множестве стратегий - позиционных стратегий. В данной статье приводится конструкция оптимальной стратегии, использующая в цепи обратной связи вспомогательную модель управляемой системы и обладающая свойством универсальности. Даны примеры, мотивирующие расширение класса разрешающих стратегий до стратегий с полной памятью.
Бесплатно
О некоторых обратных задачах для математических моделей тепломассопереноса
Статья научная
В настоящей работе рассмотрены вопросы корректности некоторых обратных задач для математических моделей, возникающих при описании процессов тепломассопереноса. По данным первой начально-краевой задачи и условию Неймана на боковой поверхности цилиндра (таким образом, на боковой поверхности цилиндра заданы данные Коши) восстанавливаются решение параболического уравнения второго порядка и коэффициент этого уравнения, принадлежащий ядру некоторого дифференциального уравнения первого порядка и характеризующий параметры среды. Неизвестный коэффициент может в том числе входить и в главную часть дифференциального оператора. Решение уравнения ищется в пространствах Соболева с достаточно большим показателем суммируемости, а неизвестный коэффициент в классе непрерывных функций. Показано, что локально по времени задача имеет единственное устойчивое решение.
Бесплатно
