Все выпуски

В программном комплексе FlowVision проведено численное моделирование процесса перемешивания неизотермических потоков натриевого теплоносителя в тройнике для обоснования применимости различных подходов — URANS (Unsteady Reynolds Averaged Navier Stokers), LES (Large Eddy Simulation) и квази-DNS (Direct Numerical Simulation) — для предсказания осциллирующего характера течения в зоне смешения и получения температурных пульсаций. Одна из основных задач данной работы — выявление преимуществ и недостатков использования этих подходов.

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

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

Затем были проведены расчеты для трех режимов течения, отличающихся соотношением расходов и температур натрия во входных сечениях тройника. Для каждого режима выполнены расчеты с применением подходов URANS, LES и квази-DNS.

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

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

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

При программной реализации модели расчеты выполнялись в среде программирования MatLab. Для решения уравнений использовался метод Рунге–Кутты–Фехлберга. При этом предполагается, что модель будет представлена в виде компьютерной исследовательской программы, позволяющей реализовать различные гипотезы о механизме наблюдаемых процессов. Рассчитаны предполагаемые величины основных показателей газообмена в условиях гиперкапнии и гипоксии. Результаты расчетов, как по характеру, так и количественно, достаточно хорошо согласуются с данными, полученными в исследованиях на испытателях. Проведенная проверка на адекватность подтвердила, что погрешность вычислений находится в пределах погрешности данных медико-биологических экспериментов. Модель можно использовать при теоретическом прогнозировании динамики респираторных реакций организма человека в условиях измененной газовой среды.

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

Известно, что внутренняя подвижность молекул ДНК играет важную роль в функционировании этих молекул. Этим объясняется большой интерес исследователей к изучению особенностей внутренней динамики ДНК. Сложность, трудоемкость и дороговизна проведения исследований в этой области стимулируют поиск и создание более простых физических аналогов, удобных для имитации различных динамических режимов, возможных в ДНК. Одно из направлений такого поиска связано с использованием механического аналога ДНК — цепочки связанных маятников. В этой модели маятники имитируют азотистые основания, горизонтальная нить, на которой подвешены маятники, имитирует сахаро-фосфатную цепочку, а гравитационное поле имитирует поле, наводимое второй нитью ДНК. Простота и наглядность — основные достоинства механического аналога. Однако модель становится слишком громоздкой в тех случаях, когда необходимо моделировать длинные (более тысячи пар оснований) последовательности ДНК. Другое направление связано с использованием электронного аналога молекулы ДНК, который лишен недостатков механической модели. В данной работе мы исследуем возможность использования в качестве электронного аналога джозефсоновскую линию. Мы рассчитали коэффициенты прямых и непрямых преобразований для простого случая однородной, синтетической ДНК, последовательность которой содержит только аденины. Внутренняя подвижности молекулы ДНК моделировалась уравнением синус-Гордона для угловых колебаний азотистых оснований, принадлежащих одной из двух полинуклеотидных цепей ДНК. При этом вторая полинуклеотидная цепь моделировалась как некоторое усредненное поле, в котором происходят эти колебания. Преобразование, позволяющее перейти от ДНК к электронному аналогу, было получено двумя способами. Первый включает две стадии: (1) переход от ДНК к механическому аналогу (цепочке связанных маятников) и (2) переход от механического аналога к электронному (линии Джозефсона). Второй способ прямой. Он включает только одну стадию — прямой переход от ДНК к электронному аналогу.

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

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

Задача демографической динамики ставится в рамках модели Мак-Кендрика – фон Ферстера на основе уравнения динамики возрастного состава. Вид функций распределения рождений, смертности и миграции населения определяется на основе имеющейся статистической информации. Приводится численное решение задачи. Представлены анализ и прогноз демографических показателей. На основе задачи демографической динамики формулируется экономико-математическая модель динамики человеческого капитала. В задаче моделирования динамики человеческого капитала рассматриваются три составляющие: образовательная, составляющая здоровья и культурная (духовная) составляющая. Для описания эволюции составляющих человеческого капитала используется двумерное уравнение типа уравнения переноса. Объемы инвестиций в составляющие человеческого капитала определяются на основе расходных статей бюджета и частных расходов с учетом характерного временного жизненного цикла демографических элементов. Для прогнозирования динамики суммарной величины человеческого капитала используется одномерное кинетическое уравнение. Приводится методика расчета динамики данного фактора как функции времени. Представлены расчетные данные по динамике человеческого капитала для Российской Федерации. Как показали исследования, величина человеческого капитала интенсивно нарастала до 2008 года, в дальнейшем наступил период стабилизации, но после 2014 года имеет место отрицательная динамика данной величины.

Теоретические и экспериментальные исследования взаимодействия водяного пара с пористыми материалами проводятся как на макро-, так и на микроуровне. На макроуровне исследуется влияние структуры расположения индивидуальных пор на процессы взаимодействия водяного пара с пористым материалом как сплошной средой. На микроуровне исследуется зависимость характеристик взаимодействия водяного пара с пористой средой от геометрии и размеров индивидуальной поры.

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

В работе рассмотрено приложение методов кинетической теории к задачам гемодинамики. Для моделирования выбраны решеточные уравнения Больцмана. Данные модели описывают дискретизированную по пространственной и временной координате динамику движения частиц на одномерной решетке. Хорошо известно, что в пределе малых длин свободного пробега решеточные уравнения Больцмана описывают уравнения гидродинамики. Если течение достаточно медленное (мало число Маха), то данные уравнения гидродинамики переходят в уравнения Навье – Стокса для сжимаемого газа. Если в получающихся гидродинамических уравнениях переменные, отвечающие плотности и скорости звука, считать площадью поперечного сечения сосуда и скоростью распространения пульсовой волны давления, то выводятся хорошо известные в биомеханике нелинейные уравнения распространения несжимаемой вязкой жидкости (крови) в эластичном сосуде для частного случая постоянной пульсовой скорости.

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

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

В статье рассмотрены режимы с обострением в социальной и биологической истории. Проведен анализ возможных причин резкого ускорения биологических и социальных процессов в определенные исторические эпохи. С использованием математического моделирования показано, что гиперболические тренды в социальной и биологической эволюции могут быть следствием переходных процессов в периоды расширения экологических ниш. Ускорение биологического видообразования связано с тем, что более ранние виды своей жизнедеятельностью изменяют среду обитания, делая ее более разнообразной, насыщая органикой, тем самым создавая благоприятные условия для появления новых видов. В социальной истории расширение экологических ниш связано с технологическими революциями, важнейшими из которых были: неолитическая революция — переход от присваивающего хозяйства к производящему (10 тыс. лет назад), «городская революция» — переход от неолита к бронзовому веку (5 тыс. лет назад), «осевое время» — переход к массовому освоению железных орудий (2.5 тыс. лет назад), промышленная революция — переход от ручного труда к машинному (200 лет назад). Все эти технологические революции сопровождались резким демографическим ростом, изменениями в социальной и политической сфе- рах. Так, наблюдаемый в последние столетия гиперболический характер роста некоторых демографических, экономических и других показателей мировой динамики — это следствие переходных процессов, начавшихся вследствие промышленной революции (замены ручного труда машинным) и предваряющих переход общества на новую стадию своего развития. Точка сингулярности гиперболического тренда характеризует окончание начального этапа этого процесса и переход к завершающей его стадии. Предложена математическая модель, описывающая демографические и экономические изменения в эпохи перемен. Показано, что прямым аналогом современной ситуации в этом смысле является «осевое время» (период с 8 века до нашей эры до начала нашей эры). Наличие такой аналогии позволяет заглянуть в будущее, изучая прошлое.

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

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

Тенсегрити-структуры отличаются рядом свойств, важных для коллаборативной робототехники, задач разведывания и движения в недетерминированных средах: естественной податливостью, компактностью при транспортировке, малым весом при значительной удароустойчивости и жесткости. При этом открытыми остаются многие вопросы управления такими структурами, что в свою очередь связано со сложностью описания их динамики.

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

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