Все выпуски

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

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

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

С помощью характеристик уравнения Гамильтона-Якоби-Беллмана описана структура множества точек, в которых минимаксное решение недифференцируемо.

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

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

Изучаются статистические характеристики множества достижимости управляемой системы, которая параметризована с помощью топологической динамической системы. Получены оценки снизу характеристик, связанных с инвариантностью заданного множества на конечном промежутке времени. Рассматривается также следующая задача, возникающая во многих приложениях. Пусть заданы числа λ₀ ∈ (0, 1] и θ > 0. Необходимо найти условия, которым должны удовлетворять управляемая система и множество X, чтобы для заданного σ ∈ Σ относительная частота поглощения множества достижимости A(t,σ,X) системы заданным множеством M на любом отрезке времени длины θ была бы не менее λ₀. Отметим, что характеристика θ предполагается заданной в зависимости от прикладной задачи. В частности, если управляемый процесс имеет периодический характер, то θ является периодом данного процесса. Результаты работы иллюстрируются на примерах управляемых систем, которые описывают различные модели роста популяции.

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

В работе рассматриваются результаты решения задачи стационарного течения вязкой несжимаемой жидкости в плоском канале с обратным уступом и прогреваемой нижней стенкой в широком диапазоне числа Рейнольдса $100\leqslant \text{Re}\leqslant 1000$ и параметра расширения потока $1.11 \leqslant ER \leqslant 10$. Исследование выполнено путем численного интегрирования системы двумерных уравнений Навье-Стокса в переменных «скорость-давление» на равномерных сетках с шагом 1/300. Достоверность полученных результатов подтверждается их сравнением с литературными данными. Приводятся подробные картины течения и перегрева жидкости в зависимости от двух основных параметров задачи: $\text{Re}$ и $ER$. Показывается, что с одновременным ростом параметров $\text{Re}$ и $ER$ существенно усложняется структура течения - увеличиваются количество вихрей и их размеры вплоть до образования вихря за уступом с двумя центрами вращения. Также показывается, что характерная высота зоны прогрева течения слабо зависит от $\text{Re}$ и $ER$ и в конечном счете ближе к выходу из канала составляет приблизительно половину его высоты. Для всех центров вихрей определяются их основные характеристики: координаты, экстремумы функции тока, завихренности. Анализируется сложное немонотонное поведение профилей коэффициентов трения, сопротивления и теплоотдачи (числа Нуссельта) по длине канала. Показывается, что эти коэффициенты в одинаковой степени сильно зависят как от числа Рейнольдса, так и от параметра расширения канала, достигая своих максимальных значений при максимальных значениях $\text{Re}$ и $ER$.

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