Все выпуски

Репрессилятором называют первую в синтетической биологии генную регуляторную сеть, искусственно сконструированную в 2000 году. Он представляет собой замкнутую цепь из трех генетических элементов — $lacI$, $\lambda cI$ и $tetR$, — которые имеют естественное происхождение, но в такой комбинации в природе не встречаются. Промотор каждого гена контролирует следующий за ним цистрон по принципу отрицательной обратной связи, подавляя экспрессию соседнего гена. В данной работе впервые рассматривается нелинейная динамика модифицированного репрессилятора, у которого имеются запаздывания по времени во всех звеньях регуляторной цепи. Запаздывание может быть как естественным, т. е. возникать во время транскрипции/трансляции генов в силу многоступенчатого характера этих процессов, так и искусственным, т. е. специально вноситься в работу регуляторной сети с помощью методов синтетической биологии. Предполагается, что регуляция осуществляется протеинами в димерной форме. Рассмотренный репрессилятор имеет еще две важные модификации: расположение на той же плазмиде гена $gfp$, кодирующего флуоресцентный белок, а также наличие в системе накопителя для белка, кодируемого геном $tetR$. В рамках детерминистского описания методом разложения на быстрые и медленные движения получена система нелинейных дифференциальных уравнений с запаздыванием на медленном многообразии. Показано, что при определенных значениях управляющих параметров единственное состояние равновесия теряет устойчивость колебательным образом. Для симметричного репрессилятора, у которого все три гена идентичны, получено аналитическое решение для нейтральной кривой бифуркации Андронова–Хопфа. Для общего случая асимметричного репрессилятора нейтральные кривые построены численно. Показано, что асимметричный репрессилятор является более устойчивым, так как система ориентируется на поведение наиболее стабильного элемента в цепи. Изучены нелинейные динамические режимы, возникающие в репрессиляторе при увеличении надкритических значений управляющих параметров. Кроме предельного цикла, отвечающего поочередным релаксационным пульсациям белковых концентраций элементов, в системе обнаружено существование медленного многообразия, не связанного с этим циклом. Долгоживущий переходный режим, который отвечает многообразию, отражает процесс длительной синхронизации пульсаций в работе отдельных генов. Производится сравнение полученных результатов с известными из литературы экспериментальными данными. Обсуждается место предложенной в работе модели среди других теоретических моделей репрессилятора.

В работе исследуются особенности групповой динамики особей-агентов в компьютерной модели популяции животных, взаимодействующих между собой и с возобновимым ресурсом. Такого типа динамика были ранее обнаружены в работе [Белотелов, Коноваленко, 2016]. Модельная популяция состоит из совокупности особей. Каждая особь характеризуется своей массой, которая отождествляется с энергией. В ней подробно описана динамика энергетического баланса особи. Ареал обитания моделируемой популяции представляет собой прямоугольную область, на которой равномерно произрастает ресурс (трава).

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

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

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

В работе проведен расчет 1/2-волнового преобразователя продольных ультразвуковых колебаний и волнового преобразователя продольно-крутильных ультразвуковых колебаний методом конечных элементов в среде программного пакета Abaqus. В результате расчета определены размеры концентратора продольно-крутильных колебаний и частотно-амплитудные характеристики преобразователей продольных и продольно-крутильных колебаний. Наложение продольно-крутильных колебаний на инструмент в процессах поверхностно-упрочняющей обработки деталей из титановых сплавов позволит уменьшить адгезионную составляющую трения в зоне контакта. Проведено сравнение результатов конечноэлементного расчета частотно-амплитудных характеристик преобразователей с экспериментальными данными, погрешность расчета не превышает 2.5 %.

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

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

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

Репрессилятором называют первую в синтетической биологии генную регуляторную сеть, искусственно сконструированную в 2000 году. Он представляет собой замкнутую малоразмерную цепь из трех генов: $lacI$, $\lambda cI$ и $tetR$, которые в такой комбинации в природе не встречаются. Промотор каждого гена контролирует следующий за ним цистрон по принципу отрицательной обратной связи, подавляя экспрессию соседнего гена. Ранее в работе [Брацун и др., 2018] была предложена математическая модель запаздывающего репрессилятора и изучены ее свойства в рамках детерминистского описания. Предполагается, что запаздывание может быть как естественным, т. е. возникать во время процессов транскрипции/трансляции в силу многоступенчатого характера этих процессов, так и искусственным, т. е. специально вноситься в работу регуляторной сети с помощью методов генной инженерии. Данная работа посвящена стохастическому описанию динамических процессов в запаздывающем репрессиляторе, которое является важным дополнением детерминистского анализа из-за сильных флуктуаций и небольшого числа молекул, принимающих обычно участие в генной регуляции. Стохастическое исследование было проведено численно с помощью алгоритма Гиллеспи, модифицированного для систем с запаздыванием. Приводятся описание алгоритма, его программная реализация и результаты тестовых расчетов для одногенного авторепрессора с запаздыванием. При исследовании репрессилятора обнаружено, что стохастическое описание в ряде случаев дает новую информацию о поведении системы, которая не сводится к детерминистской динамике даже при усреднении по большому числу реализаций. В подкритической области, где детерминистский анализ предсказывает абсолютную устойчивость системы, было обнаружено возбуждение квазирегулярных колебаний, вызываемых нелинейным взаимодействием шума и запаздывания. Выше порога возникновения неустойчивости обнаружено спонтанное изменение фазы колебаний из-за внезапной временной деградации этих колебаний. Ранее в детерминистском анализе был обнаружен долгоживущий переходный режим, который отвечает движению фазовой траектории по медленному многообразию и отражает процесс длительной синхронизации пульсаций в работе отдельных генов. Показано, что в стохастическом случае переход к кооперативному режиму работы генов репрессилятора происходит в среднем на два порядка быстрее. Построено распределение вероятности соскока фазовой траектории с медленного многообразия и определено наиболее вероятное время такого перехода. Обсуждается влияние внутреннего шума химических реакций на динамические свойства репрессилятора.

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

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

С использованием разработанного алгоритма проведено численное исследование распространения устойчивой детонации в среде с переменной плотностью. Исследован режим с относительно небольшой амплитудой колебаний плотности, при котором колебания параметров за фронтом детонационной волны происходят с частотой колебаний плотности среды. Показана связь периода колебаний параметров со временем прохождения характеристик C₊ и C₀ по области, которую условно можно считать зоной индукции. Сдвиг по фазе между колебаниями скорости детонационной волны и плотности газа перед волной оценен как максимальное время прохождения характеристики C₊ по зоне индукции.

Предлагается метод анализа тремора земной поверхности, измеряемого средствами космической геодезии с целью выделения прогностических эффектов активизации сейсмичности. Метод иллюстрируется на примере совместного анализа совокупности синхронных временных рядов ежесуточных вертикальных смещений земной поверхности на Японских островах для интервала времени 2009–2023 гг. Анализ основан на разбиении исходных данных (1047 временных рядов) на блоки (кластеры станций) и последовательном применении метода главных компонент. Разбиение сети станций на кластеры производится методом k-средних из критерия максимума псевдо-статистики. Для Японии оптимальное число кластеров было выбрано равным 15. К временным рядам главных компонент от блоков станций применяется метод разложения Хуанга на последовательность независимых эмпирических мод колебаний (Empirical Mode Decomposition, EMD). Для обеспечения устойчивости оценок волновых форм EMD-разложения производилось усреднение 1000 независимых аддитивных реализаций белого шума ограниченной амплитуды. С помощью разложения Холецкого ковариационной матрицы волновых форм первых трех EMD-компонент в скользящем временном окне определены индикаторы аномального поведения тремора. Путем вычисления корреляционной функции между средними индикаторами аномального поведения и выде- лившейся сейсмической энергии в окрестности Японских островов установлено, что всплески меры ано- мального поведения тремора предшествуют выбросам сейсмической энергии. Целью статьи является про- яснение распространенных гипотез о том, что движения земной коры, регистрируемые средствами космической геодезии, могут содержать прогностическую информацию. То, что смещения, регистрируемые геодезическими методами, реагируют на последствия землетрясений, широко известно и многократно демонстрировалось. Но выделение геодезических эффектов, предвещающих сейсмические события, является значительно более сложной задачей. В нашей статье мы предлагаем один из методов обнаружения прогностических эффектов в данных космической геодезии.

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

В работе Бурмина (2019 г.) предложена сферически-симметричная модель Земли. Она характеризуется тем, что в ней внешнее ядро обладает вязкостью и, следовательно, эффективным модулем сдвига, отличным от нуля. Для этой модели Земли проведен расчет с высокой детальностью с несущей частотой в 1 герц. В результате аналитического расчета обнаружено, что впереди PKP-волн возникают высокочастотные колебания небольшой амплитуды, так называемые предвестники. Аналитический расчет показал, что теоретические сейсмограммы для этой модели Земли во многом похожи на экспериментальные данные. При этом ключевым моментом сравнения является возникновение предвестников впереди PKP-волн. Это подтверждает правильность идей, положенных в основу ее построения.