Краткие сообщения. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование
Метод исследования диссипативных свойств разностных схем в эйлеровых координатах
Краткое сообщение
В настоящее время численные методы расчета ударно-волновых течений жидкости и газа в эйлеровых координатах получили широкое распространение, поэтому исследование их характеристик является актуальной задачей. В работе представлен подход к оценке диссипативных свойств таких разностных схем на сильных разрывах. Идея метода заключается в построении уравнения производства энтропии, погрешность аппроксимации которого может быть выражена комбинацией погрешностей аппроксимации уравнений, составляющих разностную схему. В качестве критерия диссипативности разностной схемы используется уравнение производства энтропии на слабой ударной волне. В работе проведена оценка диссипативных свойств метода крупных частиц с использованием предложенного метода.
Бесплатно
Метод функций Грина в задаче о преобразовании случайного сигнала линейной динамической системой
Краткое сообщение
Рассматривается динамическая система, описываемая линейным дифференциальным уравнением высокого порядка с постоянными коэффициентами. Методом функций Грина установлена зависимость между числовыми характеристиками случайного сигнала на входе и выходе динамической системы, а именно между математическими ожиданиями и между корреляционными функциями. В отличие от известных результатов, не предполагается стационарность входного и выходного случайных сигналов.
Бесплатно
Методы исследования устойчивости и стабилизации некоторых систем с большим запаздыванием
Краткое сообщение
Статья посвящена исследованию свойств систем дифференциальных уравнений, содержащих большое (в частности, линейное) запаздывание. Системы с линейным запаздыванием имеют достаточно широкое применение в биологии, в частности, при моделировании распределения клеток в ткани организма; а также в теории нейронных сетей. Уравнения подобного типа встречаются также в задачах физики и механики, где важным моментом является асимптотическое поведение решения (в частности, асимптотическая устойчивость). При неустойчивости таких систем возникает задача стабилизации. Оптимальный алгоритм стабилизации основан на совокупности стабилизации систем обыкновенных дифференциальных уравнений и в дальнейшем разностных систем. Данный алгоритм достаточно просто реализуется с использованием численных методов решения систем дифференциальных уравнений с запаздыванием и решения матричных уравнений. Авторами составлена программа, позволяющая достаточно эффективно находить управляющее воздействие, осуществляющее стабилизацию некоторых систем.
Бесплатно
Моделирование влияния градиентов температуры на состояние свободной поверхности жидкой пленки
Краткое сообщение
В условиях межфазной неустойчивости в пленочных аппаратах протекают различные тепло-массообменные процессы (конденсация, испарение и др.). Исследуется межфазная неустойчивость, связанная с эффектом Марангони, стекающей под действием силы тяжести жидкой пленки при умеренных числах Рейнольдса. В рамках математической модели волнового течения неизотермической жидкой пленки, а именно, системы уравнений Навье - Стокса и уравнения неразрывности с граничными условиями, учитывающими процессы тепло и-массопереноса, выведено уравнение ее свободной поверхности, коэффициенты которого учитывают градиенты температуры. Градиенты температуры вызывают неоднородность поверхностного натяжения, что приводит к возникновению на межфазной поверхности термокапиллярных сил. Модельное уравнение состояния свободной поверхности жидкой пленки - нелинейное дифференциальное в частных производных четвертого порядка решалось методом конечных разностей. Результаты вычислительных экспериментов по влиянию градиентов температуры на нелинейное развитие возмущений на свободной поверхности жидкой пленки показали направленность эффекта Марангони как на усиление возмущений и возможность разрыва пленки, так и на гашение возмущений. В условиях возрастания скорости роста возмущений жидкие пленки с небольшими числами Рейнольдса наиболее устойчивы к воздействию градиентов температуры и выдерживают большие их значения.
Бесплатно
Моделирование динамики ударного импульса в трубе с внутренним слоем водной пены
Краткое сообщение
Численно исследована динамика сферической ударной волны, инициированной взрывом в центре трубы, содержащей газ и слой водной пены с объемным водосодержанием 0,2, расположенный около внутренней границы трубы. Система модельных уравнений водной пены включает законы сохранения массы, импульса, энергии каждой фазы и уравнение динамики водосодержания пены в однодавленческом, двухскоростном, двухтемпературном приближениях в трехмерной постановке и учитывает силы межфазного сопротивления и межфазный контактный теплообмен. Термодинамические свойства воздуха и воды, составляющих газокапельную смесь, описаны реалистическими уравнениями состояния. Численное решение поставленной задачи реализовано с использованием открытого программного комплекса OpenFOAM. Достоверность расчетов по предлагаемой модели подтверждена их согласованием с литературными экспериментальными данными. Для оценки эффективности пенной защиты решена аналогичная задача о распространении сферической ударной волны в воздухе при отсутствии пенного слоя. Показано, что наличие пенной преграды снижает скорость и амплитуду ударной волны, защищая стенки трубы от воздействия взрыва.
Бесплатно
Моделирование для наведения атмосферных черенковских телескопов в эксперименте TAIGA
Статья научная
Атмосферные черенковские телескопы TAIGA-IACT являются частью гибридного экспериментального комплекса TAIGA, расположенного в 50 км от озера Байкал в Тункинской долине, республика Бурятия. Для обеспечения работы телескопов требуется с высокой точностью определять направление телескопа и положение источников на камере. Система позиционирования телескопа состоит из шаговых двигателей, датчиков положения осей и CCD камеры, установленной на тарелке телескопа. Положение телескопа определяется с помощью датчиков положения осей и применения модели наведения, а также с помощью астрометрии звездного участка снимка CCD камеры. Положение источников на камере телескопа рассчитывается с помощью модели камеры и оцененному направлению телескопа. Угол поворота фокальной плоскости телескопа рассчитывается исходя из параметров модели наведения телескопа. В данной работе представлены разработанные модели и методы, используемые для наведения телескопов TAIGA-IACT.
Бесплатно
Моделирование колебаний уровня воды в дельте дона по прогнозам ветровой ситуации
Краткое сообщение
Представлена математическая модель, состоящая из модели гидродинамики всего Азовского моря, включая Таганрогский залив, и модели движения воды в русле Дона. Данная модель позволяет рассчитывать колебания уровня воды в устьевой области Дона в зависимости от ветровой ситуации над акваторией Азовского моря и Таганрогского залива. Такой подход дает возможность заблаговременно предсказать изменение уровенной поверхности в основных рукавах реки Дон, а также в его основном русле. Гидродинамика Азовского моря описывается уравнениями мелкой воды, а движение воды в основных рукавах дельты Дона - уравнением Сен-Венана. Задача решается конечно-разностными методами. Программа написана на языке FORTRAN и реализуется на высокопроизводительных вычислительных системах в среде MPI с применением пакета Aztec. Для прогнозов ветровой ситуации использовались данные Норвежского метеорологического института. Прогнозируемая ветровая нагрузка над акваторией Азовского моря задавалась по точкам, соответствующим расположению береговых гидрометеорологических станций с последующей интерполяцией на всю акваторию. Рассмотрены две задачи прогнозирования колебаний уровня воды. Первая задача заключалась в расчете колебаний уровня в зависимости от прогноза ветровой ситуации с заблаговременностью три дня. Во второй задаче изучалась динамика изменения качества прогноза в фиксированном промежутке времени по мере приближения к его началу. Качество прогнозов определялось сравнением прогнозируемых и наблюденных значений уровня, полученных на гидропостах. Результат сравнения показал, что представленная модель достаточно адекватно описывает гидродинамику в дельтовой области Дона в зависимости от ветровой ситуации над акваторией Азовского моря.
Бесплатно
Моделирование ромбоэдрической магнитострикции в сплавах Fe-Ga
Краткое сообщение
В данной работе представлены результаты моделирования ромбоэдрической магнитострикции сплавов Fe-Ga в кристаллических структурах кубической симметрии, полученные при помощи теории функционала плотности. Показано, что зависимость энергии магнитокристаллической анизотропии от степени малых деформаций является убывающей функцией в диапазоне концентраций от 3,125 до 25 ат.% и меняет знак при величине деформации более 1% сплавах с содержанием Ga 15,625, 21,875 и 25 ат.%. Константа ромбоэдрической магнитострикции в диапазоне концентрации Ga 12,5 - 18,75 ат.% хорошо согласуется с экспериментальными данными.
Бесплатно
Модель мемристоров Бернулли в виде полинома расщепленных сигналов
Краткое сообщение
В работе представлена поведенческая модель мемристора, в котором динамика тока описывается дифференциальным уравнением Бернулли. Поведенческая модель построена в виде двумерного полинома расщепленных сигналов для передаточной характеристики мемристора Бернулли при гармоническом входном сигнале. Расщепление входных сигналов обеспечивает однозначность соответствия вход-выход, адаптацию модели к заданному классу воздействий и, следовательно, ее простоту по сравнению с универсальными нелинейными моделями, например, рядом Вольтерры и нейронными сетями. Расщепление гармонического воздействия выполнено с помощью линии задержки. Показано, что векторный сигнал, содержащий воздействие и результат его задержки во времени на один шаг, обладает минимально возможной длиной согласно условиям расщепления. Двумерный полином третьей степени, построенный на элементах векторного сигнала, обеспечивает высокую точность моделирования передаточной характеристики мемристора Бернулли в среднеквадратичной метрике.
Бесплатно
Модель оптимизации усредненной скорости потока жидкости в трубе
Краткое сообщение
Исследуется задача условной оптимизации для математической модели двухлучевого ультразвукового расходомера, позволяющая получить точную оценку средней скорости турбулетного потока жидкости, проходящей вдоль трубопровода. На основе задачи условной оптимизации разработан способ нахождения количественных характеристик (усредненной скорости турбулентной жидкости вдоль трубопровода, усредненных скоростей ультразвуковых потоков вдоль каждого акустического пути) вектора весовых коэффициентов рассматриваемой модели. На основе метода конечных элементов построен алгоритм численного метода решения поставленной задачи. Приведены результаты вычислительных экспериментов распределения скорости турбулентного потока жидкости. Найдены эпюры и профиль скорости на основе экспериментальных данных.
Бесплатно
Краткое сообщение
Целью статьи является численное исследование решения начально-краевой задачи для модели термоконвекции ненулевого порядка. Рассматривается система, которая моделирует эволюцию скорости, градиента давления и температуры несжимаемой вязкоупругой жидкости Кельвина - Фойгта ненулевого порядка. Используя метод Галеркина, разработан алгоритм численного решения начально-краевой задачи для системы, моделирующей плоскопараллельную термоконвекцию несжимаемой жидкости ненулевого порядка, и реализована программа для персональных компьютеров нахождения численного решения указанной задачи. Получена графическая иллюстрация численного решения системы при заданных параметрах. Проведенное исследование основано на результатах теории полулинейных уравнений соболевского типа, поскольку начально-краевая задача для соответствующей системы дифференциальных уравнений в частных производных сводится к абстрактной задаче Коши для уравнения соболевского типа.
Бесплатно
Краткое сообщение
Авторами статьи был разработан неитерационный метод вычисления значений собственных функций возмущенных самосопряженных операторов, названный методом регуляризованных следов (РС). Он позволяет найти значения собственных функций возмущенных операторов, зная спектральные характеристики невозмущенного оператора и собственные числа возмущенного оператора. В отличие от известных методов нахождения собственных функций, метод РС не использует матрицы и значения собственных функций находятся по линейным формулам. Это значительно увеличивает его вычислительную эффективность по сравнению с классическими методами. Для применения метода РС на практике необходимо уметь суммировать функциональные ряды Релея - Шредингера возмущенных дискретных операторов. Ранее были получены формулы нахождения взвешенных поправок теории возмущений, что позволяло приближенно находить суммы функциональных рядов Релея - Шредингера, заменяя их частичными суммами, состоящими из этих поправок. В статье впервые получены формулы нахождения значений сумм функциональных рядов Релея - Шредингера возмущенных дискретных операторов в узловых точках. Проведены вычислительные эксперименты по нахождению значений собственных функций возмущенного одномерного оператора Лапласа. Результаты эксперимента показали высокую вычислительную эффективность разработанного метода суммирования рядов Релея - Шредингера.
Бесплатно
Начально-конечная задача для линейной стохастической модели Хоффа
Краткое сообщение
Линейная модель Хоффа, исследующая в линейном приближении динамику выпучивания двутавровых балок в конструкции, представляет собой множество линейных одномерных уравнений Хоффа, заданных на ребрах геометрического графа с условиями непрерывности и баланса потоков в его вершинах. Ранее детерминированная модель изучалась в разных аспектах многими специалистами. Стохастическая модель изучается впервые. В качестве метода исследования используется классический подход Ито - Стратоновича - Скорохода, распространенный на гильбертовы пространства и уравнения соболевского типа. Основной результат - теорема об однозначной разрешимости начально-конечной задачи с аддитивным белым шумом, под которым понимается обобщенная производная K-винеровского процесса. Решение представлено в виде формул, допускающих постановку вычислительных экспериментов.
Бесплатно
О конечном спектре трехточечных краевых задач
Краткое сообщение
Статья посвящена решению одной из проблем Джона Локкера, а именно вопросу, может ли краевая задача для дифференциального уравнения иметь конечный спектр. На задачи такого рода можно смотреть и как на обратные задачи определения коэффициентов дифференциального уравнения по заданному спектру. В работе показано, что если дифференциальное уравнение не имеет кратных корней характеристического уравнения, то тогда спектр соответствующей трехточечной краевой задачи не может быть конечным. Доказательство теоремы основано на том, что соответствующий характеристический определитель представляет собой целую функцию класса К, а так же результатах В.Б. Лидского и В.А. Садовничего, из которых следует, что количество корней характеристического уравнения (если они есть) бесконечно. Если же корни характеристического уравнения являются кратными, то спектр может быть конечным. Более того, существуют краевые задачи с наперед заданным конечным спектром.
Бесплатно
О моделировании деформаций пластин
Краткое сообщение
Рассматривается эллиптические краевые задачи четвертого порядка, лежащие в основе математических моделей деформаций пластин на упругих основаниях при смешанных краевых условиях четырех теоретически возможных типов. Предлагаются замещения этих задач в вариационной форме на их фиктивные продолжения. Решения последних задач с помощью модификаций методов фиктивных компонент сводятся к решениям задач в прямоугольной области. Приводятся оптимальные оценки сходимости итерационных процессов на непрерывном уровне. При простой дискретизации фиктивно продолженных задач по методу конечных элементов на параболических восполнениях получаются эффективные численные модификации методов фиктивных компонент простые при практической реализации на ЭВМ. Получаемые системы линейных алгебраических уравнений могут оптимально решаться с помощью методов итерационных факторизаций. В итоге предложенные численные методы являются логарифмически оптимальными или оптимальными по количеству арифметических операций, необходимых для достижения задаваемых относительных погрешностей.
Бесплатно
Об одной полулинейной математической модели соболевского типа высокого порядка
Краткое сообщение
В статье исследуется полулинейная математическая модель соболевского типа высокого порядка с относительно спектрально ограниченным оператором. Данная математическая модель строится на основе уравнения соболевского типа высокого порядка и условий Коши. В работе используются метод фазового пространства и теория относительно p-ограниченных операторов, разработанные Г.А. Свиридюком. При исследовании невырожденной математической модели используется подход, предложенный С. Ленгом; в статье он обобщается на дифференциальные уравнения высокого порядка. В работе рассмотрено два случая. В первом, когда оператор при старшей производной по времени является непрерывно обратимым, используются методы теории дифференцируемых банаховых многообразий и доказывается однозначная разрешимость задачи Коши. Во втором случае, когда оператор при старшей производной по времени имеет нетривиальное ядро. Как известно, задача Коши для уравнений соболевского типа принципиально не разрешима при произвольных начальных данных. В связи с этим возникает задача построения фазового пространства уравнения как множества допустимых начальных значений, содержащего решения уравнения, и изучения его морфологии. В данной работе для вырожденного уравнения строится локальное фазовое пространство.
Бесплатно
Обратная задача в управлении динамической системой
Краткое сообщение
Предложен синтез управлений движением манипуляционного робота (МР) в точку программной траектории прямым методом Ляпунова. При построении модели динамики МР представлен единой системой, фазовый вектор которой определяется данными как механизма, так и исполнительных приводов. Такой подход к синтезу управлений базируется на использовании в качестве функций Ляпунова первых интегралов движения системы, поэтому включаемый в рассмотрение привод должен иметь конкретный физический смысл, например, электропривод постоянного тока, гидропривод с дроссельным управлением и т.д., что позволит учесть его энергию при построении функции Ляпунова. Для определенности рассмотрен МР с электроприводами, с якорным управлением, что не ограничивает возможности перехода к другому типу привода. Вывод уравнений движения МР проведен на основе полной нелинейной модели МР с использованием тензорного анализа. Исследовано на устойчивость движение промышленного МР, рабочим органом которого является тело вращения.
Бесплатно
Ограниченные решения модели Баренблатта - Желтова - Кочиной в квазисоболевых пространствах
Краткое сообщение
Уравнения соболевского типа в банаховых пространствах изучены довольно полно. Квазисоболевы пространства - это квазинормируемые полные пространства последовательностей. Уравнения соболевского типа в таких пространствах начали изучаться совсем недавно. В данной статье рассматривается вопрос существования ограниченных на всей оси решений для модели Баренблатта - Желтова - Кочиной. Кроме введения и списка литературы, статья содержит две части. В первой содержатся предварительные сведения о свойствах операторов в квазибанаховых пространствах, а также об относительно ограниченных операторах. Во второй части приведен основной результат статьи о существовании ограниченных решений для модели Баренблатта - Желтова - Кочиной в квазисоболевых пространствах. Список литературы не претендует на полноту и отражает лишь вкусы и пристрастия авторов.
Бесплатно
Статья научная
В статье предлагается использование цифрового одномерного фильтра Калмана в реализации численных алгоритмов решения задачи оптимальных динамических измерений для восстановления динамически искаженного сигнала при наличии помех. Математическая модель сложного измерительного устройства построена как система леонтьевского типа, начальное состояние которой отражает условие Шоуолтера - Сидорова. Основным положением теории оптимальных динамических измерений является моделирование искомого входящего сигнала как решение задачи оптимального управления с минимизацией функционал штрафа, в котором оценивается расхождение моделируемого и наблюдаемого выходящего (или наблюдаемого) сигнала. Наличие помех на выходе измерительного устройства приводит к необходимости использования в численных алгоритмах цифровых фильтров. Сглаживающие фильтры, применяющиеся при неизвестных вероятностных параметрах помех, недостаточно эффективны при фильтрации пикообразных сигналов на малом временном промежутке. Кроме того, динамика измерений актуализирует рассмотрение фильтров, реагирующих на быстро меняющиеся данные. В статье предлагается включение процедуры фильтрации наблюдаемого сигнала в ранее разработанные численные алгоритмы, что позволяет либо расширить их применение, либо упростить функционал штрафа.
Бесплатно
Краткое сообщение
В работе исследована задача оптимального управления для неполного уравнения соболевского типа высокого порядка. Доказана теорема существования и единственности сильного решения задачи Коши для данного уравнения. Получены достаточные условия существования и единственности оптимального управления такими решениями. В работе используются идеи и методы, разработанные Г.А. Свиридюком и его учениками. Доказательство теоремы о существовании и единственности оптимального управления для исследуемой задачи опирается на теорию оптимального управления, развитую в работах Ж.-Л. Лионса.
Бесплатно
