Все выпуски

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

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

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

Функциональное моделирование процессов свертывания крови, в частности возникновения фибрин–полимерных сгустков, имеет большое значение для прикладных вопросов медицинской биофизики. Несмотря на некоторые неточности в математических моделях, качественные результаты представляют огромный интерес для экспериментаторов как средство анализа возможных вариантов развития их работ. При достижении хорошего количественного совпадения с экспериментальными результатами такие модели могут быть использованы для технологических применений. Целью данной работы является моделирование процесса многоступенчатой полимеризации фибрина и сопряженного с ними золь-гель-перехода — возникновения фибрин-полимерной сетки в точечной системе. Для программной реализации и численных экспериментов используется неявный метод Розенброка второго порядка с комплексными коэффициентами (CROS). В работе представлены результаты моделирования и проведен анализ чувствительности численных решений к коэффициентам математической модели методами вариации. Показано, что в физиологическом диапазоне параметров констант модели существует лаг-период 20 секунд между началом реакции и возникновением зародышей фибрин-полимерной сетки, что хорошо соответствует экспериментальным наблюдениям подобных систем. Показана возможность появления нескольких $(n = 1–3)$ последовательных золь-гель-переходов. Такое необычное поведение системы является прямым следствием наличия нескольких фаз в процессе полимеризации фибрина. На последнем этапе раствор олигомеров фибрина длины 10 может достичь полуразбавленного состояния. Это, в свою очередь, приведет к исключительно быстрой кинетике формирования фибрин-полимерной сетки, управляемой вращательной диффузией олигомеров. Если же состояние полуразбавленного раствора не достигается, то образование фибрин-полимерной сетки контролируется трансляционной диффузией, которая является существенно более медленным процессом. Такой дуализм в процессе золь-гель-перехода привел к необходимости введения функции переключения в уравнения для кинетики образования фибрин-полимера. Ситуация с последовательными золь-гель-переходами соответствует экспериментальным системам, где вследствие физических процессов, таких как пресипитация, фибрин-полимерная сетка может быть быстро удалена из объема.

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

Предложена трехкомпонентная модельпланк тонного сообщества с дискретным временем. Сообщество представлено зоопланктоном и двумя конкурирующими за ресурсы видами фитопланктона: токсичным и нетоксичным. Модельдв ух связанных уравнений Рикера, ориентированная на описание динамики конкурентного сообщества, используется для описания взаимодействия двух видов фитопланктона и позволяет неявно учитывать ограничение роста биомассы каждого из видов-конкурентов доступностью внешних ресурсов. Изъятие фитопланктона за счет питания зоопланктоном описывается трофической функцией Холлинга II типа с учетом насыщения хищника. Способность фитопланктона защищаться от хищничества и избирательность питания хищника учтены в виде ограничения потребления: зоопланктон питается только нетоксичным фитопланктоном.

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

Работа посвящена изучению волновых и релаксационных эффектов при импульсном истечении смеси газа с большим содержанием твердых частиц из цилиндрического канала при его начальном частичном заполнении. Задача сформулирована в двухскоростной двухтемпературной постановке и решалась численно гибридным методом крупных частиц второго порядка аппроксимации. Численный алгоритм реализован в виде параллельных вычислений с использованием базовых языковых средств Free Pascal. Применимость и точность метода для волновых потоков концентрированных газовзвесей подтверждены сопоставлением с тестовыми асимптотически точными решениями. Погрешность расчета на сетке невысокой детализации вх арактерных зонах течения двухфазной среды составила 10⁻⁶ . . . 10⁻⁵.

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

Изучено влияние релаксационных свойств при изменении размеров частиц на динамику ограниченного слоя газодисперсной среды. Сопоставление истечения ограниченного слоя газовзвеси с различными размерами частиц показывает, что для мелких частиц (число Стокса меньше 0,001) наблюдается аномальное явление одновременного существования ударно-волновых структур в сверх- и дозвуковом потоке газа и взвеси. С увеличением размеров дисперсных включений скачки уплотнения в области двухфазной смеси сглаживаются, а для частиц (число Стокса больше 0,1) — практически исчезают. При этом ударно-волновая конфигурация сверхзвукового газового потока на выходе из канала сохраняется, а положения и границы энергонесущих объемов газовзвеси при изменении размеров частиц близки.

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

В данной работе рассматривается одна из самых значимых моделей популяционной генетики — модель Кроу–Кимуры. В последнее десятилетие были исследованы модели с ландшафтами приспособленности малой размерности. Цель статьи состоит в анализе модели Кроу–Кимуры c многомерным ландшафтом приспособленности в рамках формализма Гамильтона–Якоби. Для случая однопикового ландшафта приспособленности выводятся точные аналитические выражения, которые подтверждаются численно.

В представленной работе обсуждается возможность упрощенного трехмерного нестационарного моделирования процесса плазменно-стимулированного горения газообразного топлива в сверхзвуковом потоке воздуха. Расчеты проводились в программном комплексе FlowVision. В работе выполнен анализ геометрии эксперимента и сделан вывод о ее существенной трехмерности, связанной как с дискретностью подачи топлива в поток, так и с наличием локализованных плазменных образований. Предложен вариант упрощения расчетной геометрии, основанный на симметрии аэродинамического канала и периодичности пространственных неоднородностей. Выполнено тестирование модифицированной $k–\varepsilon$ модели турбулентности FlowVision (KEFV) в условиях сверхзвукового потока. В этих расчетах в области источников тепла и инжекции топлива использовалась подробная сетка без пристеночных функций, а на удаленных от ключевой области поверхностях пристеночные функции были включены. Это позволило существенно уменьшить количество ячеек расчетной сетки. Сложная задача моделирования воспламенения углеводородного топлива при воздействии плазмы была существенно упрощена путем представления плазменных образований как источников тепла и использования одной брутто-реакции для описания горения топлива. На базе геометрии аэродинамического стенда ИАДТ-50 ОИВТ РАН с помощью моделирования в программном комплексе ПК FlowVision проведены калибровка и параметрическая оптимизация подачи газообразного топлива в сверхзвуковой поток. Продемонстрировано хорошее совпадение экспериментальной и синтетической теневой картины потока при инжекции топлива. Проведено моделирование потока для геометрии камеры сгорания Т131 ЦАГИ с инжекцией топлива и генерацией плазмы. В результате моделирования для заданного набора параметров продемонстрировано воспламенение топлива, что совпало с результатами эксперимента. Отмечена важность адаптации расчетной сетки с повышением пространственного разрешения в области объемных источников тепла, моделирующих зону электрического разряда. Достигнуто удовлетворительное качественное совпадение распределений давления, полученных в моделировании и эксперименте.