Все выпуски

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

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

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

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

Перспективным методом повышения количества осадков в засушливом климате является способ создания вертикальной высокотемпературной струи, насыщенной гигроскопическим аэрозолем. Такая установка позволяет создавать искусственные облака с возможностью образования осадков в безоблачной атмосфере, в отличие от традиционных способов искусственного увеличения осадков, в которых предусматривается повышение эффективности осадко-образования только в естественных облаках путем их засева ядрами кристаллизации и конденсации. Для увеличения мощности струи добавляются хлорид кальция, карбамид, пищевая соль в виде грубодисперсного аэрозоля, а также нанопорошок NaCl/TiO₂, который способен конденсировать значительно больше водяного пара, чем перечисленные типы аэрозолей. Дисперсные включения в струе также являются центрами кристаллизации и конденсации в создаваемом облаке для повышения возможности осадкообразования. Для моделирования конвективных течений в атмосфере применяется математическая модель атмосферных течений большого масштаба FlowVision, решение уравнений движения, энергии и массопереноса проводится в относительных переменных. Рассматриваемая постановка задачи разделена на две части: модель начальной струи и постановка атмосферных течений большого масштаба FlowVision. Нижняя область, где происходит течение начальной высокоскоростной струи, моделируется в сжимаемой постановке с решением уравнения энергии относительно полной энтальпии. Данное разделение задачи на две отдельные подобласти необходимо, чтобы корректно провести численный расчет начальной турбулентной струи при высокой скорости (M > 0,3). Приводятся основные математические зависимости модели. С использованием представленной модели проведены численные эксперименты, для исходных данных взяты экспериментальные данные из натурных испытаний установки по созданию искусственных облаков, проведенные в Объединенных Арабских Эмиратах. Получено хорошее согласие с экспериментом: в 55% проведенных расчетов значение вертикальной скорости на высоте 400 м (более 2 м/с) и высота подъема струи (более 600 м) находятся в пределах погрешности 30% от экспериментальных характеристик, а в 30% расчетах — полностью согласуются с экспериментом. Результаты численного моделирования позволяют оценить возможность использования метода высокоскоростной струи для стимулирования искусственной конвекции и, в конечном итоге, для создания осадков. Расчеты проведены с использованием программного комплекса FlowVision на суперкомпьютере «Торнадо ЮУрГУ».

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

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

В дальнейшем планируется использовать разработанный инструментарий для анализа конкретных ситуаций, в том числе для анализа закономерностей цивилизационных циклов: взлета и падения Римской империи, СССР, современной западной цивилизации.

Рассматривается пространственно-временная модель тритрофической системы, описывающая взаимодействие жертвы, хищника и суперхищника в среде с неоднородным распределением ресурса. Учитываются всеядность суперхищника (Intraguild Predation, IGP), диффузия и направленная миграция (таксис), который моделируется с помощью логарифмической функции от ресурса и плотности жертвы. Основное внимание уделено анализу мультистабильности системы и роли косимметрии в формировании континуальных семейств стационарных решений. С использованием численно-аналитического подхода изучаются пространственно-однородные и неоднородные стационарные решения. Установлено, что при выполнении дополнительных соотношений между параметрами, характеризующими локальное взаимодействие хищников, и коэффициентами диффузии система обладает косимметрией, что приводит к возникновению семейства устойчивых стационарных решений, пропорциональных функции ресурса. Показано, что косимметрия не зависит от функции ресурса в случае неоднородной среды. Проведено исследование устойчивости стационарных распределений с помощью спектрального метода. Нарушение условий косимметрии приводит к разрушению семейства и появлению изолированных стационарных состояний, а также к длительным переходным процессам, отражающим память системы об исчезнувшем семействе. В зависимости от начальных условий и параметров в системе реализуются переходы к режимам с одним хищником (выживание хищника или суперхищника) или к сосуществованию хищников. Численные эксперименты на основе метода прямых (разностная схема по пространственной переменной и метод Рунге – Кутты для интегрирования по времени) подтверждают мультистабильность системы и иллюстрируют исчезновение семейства решений при разрушении косимметрии.

Представлена математическая модель, описывающая динамику инфекции ВИЧ-1 в отдельно взятом лимфоузле в начальный период развития инфекции. В рамках модели инфицирование индивидуума задается неотрицательной финитной функцией, описывающей скорость поступления первоначальных вирусных частиц в лимфоузел. Уравнения модели построены с учетом следующих факторов: 1) взаимодействие вирусных частиц с наивными Т-лимфоцитами CD4+, находящимися в различных фазах клеточного цикла; 2) контактное взаимодействие между размножающимися наивными Т-лимфоцитами CD4+ и инфицированными Т-лимфоцитами CD4+, производящими вирусные частицы. Спецификой контактных межклеточных взаимодействий является образование комплексов, состоящих из пар указанных клеток. Длительности существования комплексов задаются функциями распределения на конечных промежутках времени. Модель записана в форме высокоразмерной системы нелинейных дифференциальных уравнений с запаздыванием, включая два уравнения с распределенным запаздыванием, и дополнена неотрицательными начальными данными. При отсутствии инфекции ВИЧ-1 модель сводится к четырем дифференциальным уравнениям с запаздыванием, описывающим численность наивных Т-лимфоцитов CD4+ в различных фазах клеточного цикла. Показана глобальная разрешимость модели (существование и единственность решения на полуоси) и установлена неотрицательность компонент решения. Для проведения вычислительных экспериментов с моделью разработан алгоритм численного решения используемой системы дифференциальных уравнений на основе полунеявной схемы Эйлера для случая равномерного распределения длительностей существования комплексов. Представлены результаты вычислительных экспериментов, направленных на приближение численного решения модели к описанию кинетики развития инфекции ВИЧ-1 в ее острой фазе, включая фазу эклипса. В качестве наблюдаемой использована переменная, описывающая количество вирусных частиц на один миллилитр крови на 10–12-е сутки после начала острой инфекции. Численно исследована динамика наблюдаемой переменной в зависимости от вариации параметров модели, отражающих закономерности формирования комплексов и образования клеток, производящих вирусные частицы. Показана возможность затухания инфекции ВИЧ-1 в лимфоузле при определенных значениях некоторых из параметров модели.

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

На основе математическоймо дели динамики биомасс фитопланктона построен способ оценки содержания хлорофилла в районе моря с учетом его распределения по глубине. Модель построена на основе уравнения «реакция-диффузия», учитывает основные влияющие факторы: минеральное питание, освещенность и температуру. Используется спутниковая информация о поверхностном слое моря. Приведен пример расчетов для залива Петра Великого (Японское море).