Все выпуски

В работе представлена новая версия программного комплекса FlowVision, предназначенного для автоматизации инженерных расчетов в области вычислительной гидродинамики: FlowVision 3.09.05. Программный комплекс (ПК) FlowVision используется для решения различных прикладных задач в различных областях промышленности. Его популярность основана на том, что он позволяет решать сложные нетрадиционные задачи, находящиеся на стыке различных дисциплин, с одной стороны, и, с другой стороны, на парадигме полной автоматизации таких трудоемких для инженера процессов, как построение расчетной сетки. FlowVision — это программный комплекс, полностью отчуждаемый от разработчиков. Он имеет развитый графический интерфейс, систему задания расчетного проекта и систему визуализации течений различными методами — от построения контуров (для скалярных переменных) и векторов (для векторных переменных) на плоскостях и поверхностях до объемной визуализации расчетных данных. Кроме этого, ПК FlowVision предоставляет пользователю возможность вычислять интегральные характеристики на поверхностях и в ограниченных объемах.

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

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

В данной работе предложены три возможных механизма формирования сегментированных волн и спиралей. Структуры такого рода были обнаружены в реакции Белоусова–Жаботинского, диспергированной в обращенной микроэмульсии аэрозоля OT. Первый механизм обусловлен взаимодействием двух подсистем, одна из которых возбудима, а другая неустойчива по Тьюрингу. Показано, как под воздействием поперечной неустойчивости из однородной гладкой спиральной волны формируется сегментированная спираль. В зависимости от свойств подсистем мы демонстрируем несколько различных по виду и форме сегментированных спиральных волн. В качестве второго механизма мы предлагаем «дробление» бегущей волны в окрестности бифуркационной точки коразмерности два, в которой пересекаются границы тьюринговской и волновой неустойчивостей. Наконец, мы показываем, что сегментированные волны могут возникать в некоторых простых двухкомпонентных моделях типа «реакция–диффузия», имеющих более одного стационарного состояния, в частности, в модели ФитцХью–Нагумо.

Предложен формальный модельный механизм формирования осциллонов, которые были обнаружены во множестве физических систем, а также в химической реакции Белоусова–Жаботинского, протекающей в обращенной масляной микроэмульсии аэрозоля ОТ. В предложенном механизме возникновение осциллонов происходит в результате взаимодействия двух подсистем. В первой подсистеме при подходящем наборе параметров в результате жесткого локального возбуждения возможно образование уединенных стационарных структур, которые определяют пространственное распределение параметра второй подсистемы, изменение которого вызывают в ней локальные осцилляции.

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

В приближении однородной намагниченности построена математическая модель трехслойной ячейки памяти MRAM c осью анизотропии, расположенной перпендикулярно запоминающему ферромагнитному слою ячейки (перпендикулярная анизотропия). Предполагается, что первоначально намагниченность свободного слоя ячейки ориентирована вдоль оси анизотропии и соответствует состоянию «нуль». Одновременное мгновенное включение спин-поляризованного тока и магнитного поля воздействует на намагниченность свободного слоя и может перевести ее в противоположное положение, соответствующее состоянию «единица». Математическое описание эффекта основано на классическом векторном уравнении Ландау–Лифшица с диссипативным членом в форме Гильберта. В нашей модели учтены взаимодействия намагниченности с внешним магнитным полем и эффективными полями анизотропии и размагничивания, а также с током инжекции в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от управляющих параметров: величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Установлено, что в данной системе, в отличие от системы с продольной анизотропией, дополнительные состояния равновесия отсутствуют. Проведен анализ устойчивости основных состояний равновесия по первому приближению. Построены бифуркационные диаграммы, характеризующие типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно, методом Рунге–Кутты, построены траектории переключения. Найдены комбинации управляющих параметров, при которых переключение невозможно. Найдены области существования устойчивых и неустойчивых предельных циклов системы. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Проведено сравнение значений порогового тока в моделях с продольной и перпендикулярной анизотропией при нулевом магнитном поле и показано, что в модели с перпендикулярной анизотропией ток переключения почти на порядок ниже, чем в модели с продольной анизотропией.