Все выпуски

Математическое моделирование композиционных материалов играет важную роль в современной технике, а решение и исследование обратных граничных задач теплообмена невозможно без использования систем собственных функций задачи Штурма-Лиувилля для дифференциального уравнения с разрывными коэффициентами. Одним из важнейших свойств таких систем является их полнота в соответствующих пространствах. Это свойство систем позволяет доказать теоремы существования и единственности как для прямых задач, так и обратных граничных задач теплопроводности, а также обосновать численные методы решения таких задач. В настоящей статье доказана полнота в пространстве $L_2[r_0,r_2]$ задачи Штурма-Лиувилля для дифференциального оператора второго порядка с разрывным коэффициентом. Эта задача возникает при исследовании и решении обратной граничной задачи теплопроводности для полого шара, состоящего из двух шаров с различными коэффициентами температуропроводности. Доказана самосопряженность, инъективность, а также положительная определенность этого оператора.

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

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

Численно исследуется явление возникновения скоростного бафтинга при обтекании профиля NACA0012 трансзвуковым потоком. Формулируется математическая модель, основанная на алгоритмах высокого порядка аппроксимации, позволяющая рассчитывать нестационарные отрывные течения. Модель базируется на интегрировании квазигидродинамических уравнений. Проводится параметрическое исследование обтекания профиля высокоскоростным потоком вязкого газа в зависимости от угла атаки. Анализируются как мгновенные, так и осредненные картины течения. Получены распределения пульсационных характеристик течений при различных углах атаки. Выявляются закономерности возникновения отрыва пограничного слоя, определено влияние скачков уплотнения на характер течения вблизи поверхности профиля. Определяется критический угол атаки, при котором начинает иметь место скоростной бафтинг.

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

Рассматриваются два подхода к решению задачи математического моделирования обтекания метаемых тел: численное решение уравнений движения сплошной среды Навье-Стокса, осредненных по Рейнольдсу (RANS - Reynolds-averaged Navier–Stokes), с использованием модели турбулентности и прямое численное моделирование (DNS - Direct Numerical Simulation). Тестирование рассматриваемых подходов проводится при решении задачи обтекания тел вращения с простой геометрией: сферы и цилиндра с конической головной частью, для которых известны значения коэффициентов сопротивления при различных числах Маха. Проведено качественное и количественное сравнение результатов обтекания рассматриваемых тел сверхзвуковым потоком, полученным по методикам RANS и DNS. Апробация методики численного моделирования проводится для метаемого тела (снаряда) характерной формы. Представлены результаты численного моделирования обтекания снаряда для широкого диапазона параметров: чисел Маха и углов нутации. Выполнено сравнение расчетных значений коэффициентов сопротивления с эмпирическими эталонными зависимостями по законам 1943 и 1958 годов.