Моделирование, информатика и управление. Рубрика в журнале - Математическая физика и компьютерное моделирование

Развертывание региональной климатической модели для юга России на основе regcm 4.5
Статья научная
В работе описан первый этап развертывания региональной климатической модели для Юга России на основе RegCM 4.5, что требует актуализации набора схем параметризации подсеточных физических процессов для заданной территории. Моделирование общей циркуляции атмосферы требует высокой скорости счета радиационных характеристик, что предполагает высокую степень параметризации уравнения переноса излучения, в частности, характеристик молекулярного поглощения. Высокая селективность молекулярного поглощения в сравнении с аэрозольным и молекулярным рассеянием в видимом и инфракрасном диапазоне спектра сильно усложняет эту проблему. Подробно описаны схемы расчетов в радиационном блоке RegCM. Показано, что ряды экспонент обеспечивают высокую степень параметризации и высокую точность расчета, сопоставимую с методом line-by-line. Наши вычисления потоков длинноволнового излучения с различными версиями спектроскопической базы данных HITRAN показали, что различия между нисходящими потоками на нижней границе атмосферы для лета умеренных широт не превышали 0,1 Вт/м2, а для восходящих потоков на верхней границе атмосферы - не более 0,6 Вт/м2. Значения восходящего потока на верхней границе атмосферы при использовании спектроскопической базы данных HITRAN2012 дало величину 281,01 Вт/м2, а для нисходящего потока - 351,16 Вт/м2. Эти оценки позволяют сделать вывод о том, что современные спектроскопические базы данных обеспечивают высокую точность расчета широкополосных потоков. Для того чтобы оценить неопределенность в потоках излучения, мы использовали модель континуального поглощения паров Н2О, которая применялась во всем длинноволновом спектральном диапазоне, за исключением окна прозрачности 830-1130 см-1, где она заменялась на разработанную нами модель, основывающуюся на экспериментальных данных. Начальные и граничные условия определяются из глобальной климатической модели European Centre for Medium-Range Weather Forecast’s ERA-Interim (EIN15), а для задания температуры поверхности морей применяются данные Indian Ocean Sea Surface Temperature (IOSST). Этап верификации включает: анализ выбора вычислительной области, определение точности вертикальных профилей параметров, сравнение температурных полей с использованием данных наблюдений. Для оценки качества региональных климатических моделей мы предлагаем анализировать возможность появления особых метеорологических событий, прежде всего, связанных с сильными осадками. Такой подход представляется более надежным по сравнению с традиционным рассмотрением распределений температуры.
Бесплатно

Разработка классификатора фотоизображений патологий для ультрамалого набора данных
Статья научная
Цель работы: создать алгоритм и реализовать его в программном средстве для классификации фотоизображений патологий центральной области глазного дна человека, выявляемых при помощи исследования аутофлюоресценции, по 8 типам-паттернам: нормальный, минимальные изменения, фокальный, пятнистый, линейный, кружевоподобный, ретикулярный, крапчатый. Методы: алгоритмы машинного обучения (сверточные нейронные сети) и компьютерного зрения (гистораммные методы, перцептивные хэш-алгоритмы). Главная особенность задачи: ультрамалый набор уникальных фотоизображений с точно диагностируемым типом патологии (18 штук). Точность прогнозов при решении задачи с помощью нейросети 12,5 %. Точность прогнозов разработанного алгоритма с использованием комбинации гистограмм, перцептивного хэша и 1 опорного фото нормального состояния глазного дна равна 60% при подборе параметров классификатора из набора 1 фото на 1 патологию. При использовании 3 опорных фото нормы - 85 %. Предложенное решение может использоваться в медицине, офтальмологии, фотонике и оптике биотканей, машинном обучении как в научно-исследовательских, так и учебных целях.
Бесплатно

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

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

Статья научная
Моделирование процесса зарастания мелководий высшей водной растительностью имеет важное практическое значение для рыбохозяйственной отрасли нашей страны и является неотъемлемой частью мониторинга водных объектов. В статье представлены результаты разработки системы компьютерного зрения средствами машинного обучения на базе архитектур SegNet и U-Net для оценки зарастания высшей водной растительностью Цимлянского водохранилища. Для обучения и тестирования моделей применялся набор данных, состоящий из 200 пар снимков Landsat, охватывающих 24 различных участка Цимлянского водохранилища за разные годы, а также соответствующих им разметок зарастания. Процесс обучения SegNet продолжался в течение 50 эпох, U-Net обучалась в течение 30 эпох. Каждая эпоха обучения включала итерации по обучающим данным, вычисление функции потерь, обратное распространение градиентов и обновление весов с использованием оптимизатора. После каждой эпохи производилась валидация модели на валидационной выборке для оценки ее производительности. Точность модели SegNet составила 0,869, U-Net – 0,881. Для оценки качества сегментации зарастания на тестовой выборке были измерены коэффициенты Жаккара (IoU). Модель U-Net показала IoU на уровне 0,665, у SegNet этот показатель составил 0,633.
Бесплатно

Расчет 3D-формы гиперупругого тела для моделей нелинейной теории упругости методом Ньютона
Статья научная
Рассматриваются способы расчета деформаций объектов с гиперупругими материалами в рамках нелинейной теории упругости. В статье приводится алгоритм решения данной задачи, как минимизация функционала запасенной энергии, методом Ньютона, дается его подробное описание и способы реализации. Было разработано программное обеспечение, реализующее данный метод и позволяющее производить численные эксперименты и компьютерное моделирование деформаций гиперупругих тел.
Бесплатно

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

Симуляционное моделирование движения конечности человека
Статья научная
Как известно, моделирование любых процессов основывается на некоторых закономерностях, происходящих как внутри объекта моделирования, так и снаружи (изменение окружающей среды, в которой находится объект). При изучении сложных биосистем выявление таких закономерностей осложняется тем фактом, что подобные системы обладают хаотической структурой. В таких системах невозможно произвольно повторить начальное состояние 𝑥𝑖, любое промежуточное и конечное 𝑥𝑘. Соответственно и моделирование сложных биосистем должно производиться на основе случайных закономерностей. Созданная симуляционная модель позволяет максимально точно моделировать движение конечности человека в двумерной плоскости. Выделяются три принципиальные основополагающие идеи в работе модели: во-первых, включение регуляторных механизмов модели осуществляется за счет нахождения Ф-решений; во-вторых, траектория удержания конечности в пространстве (моделирование тремора человека) изменяет свое направление случайным образом; в-третьих, все изменения системы происходят на основе случайных чисел (в модели нет констант). Основополагающая часть симуляционной модели базируется на поиске Ф-решений. Созданная симуляционная модель является масштабируемой, что позволит в дальнейшем исследовании перевести ее к трехмерному варианту с возможностью изменять количество мышечных пучков, задействованных в движении конечности человека.
Бесплатно

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

Сценарии осцилляций и гибели в новой непрерывной модели эруптивной фазы инвазии чужеродного вида
Статья научная
Рассматривается проблема моделирования развития особых популяционных процессов, которые включают прохождение эруптивной фазы динамики. Подобные непродолжительные, но ураганные режимы часто связаны с последствиями инвазий нежелательных биологических видов. Процессы при вселении вида часто могут развиваться через отложенную во времени фазу стремительного увеличения его численности. Завершение фазы зависит от многих факторов. Вспышки отдельных биологических видов оказывают столь сильное давление на среду, что достижение ненулевого балансового равновесия проблематично. Подобные явления трактуются нами как переходный процесс с эруптивной фазой к неопределенному заранее состоянию биотической среды. В зависимости от противодействия, что ярко видно на примерах динамики насекомых-вредителей, сценарии подобных явлений могут развиваться различным образом, в том числе с разрушением среды обитания. Разработана новая модель на основе уравнения с отклоняющимся аргументом, где описан вариант развития повторной вспышки катастрофического характера. Сценарий реализуется при возникновении негармонического цикла 𝑁*(𝑟, 𝑡), который при специфических условиях не может быть орбитально устойчивым, но становится переходным. Цикл завершается тривиальным значением. Моделируемый сценарий наиболее резкой формы эруптивной фазы заканчивается в вычислительном эксперименте полной гибелью инвазионной популяции, но без образования неограниченной сверху траектории из колебаний. Гибель происходит при разрушении релаксационного цикла экстремальной амплитуды в рассмотренном нами уравнении популяционной вспышки в предыдущей нашей работе [2].
Бесплатно

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

Цепной алгоритм сжатия 3D-моделей
Статья научная
В статье подробно излагается алгоритм сжатия информации о геометрическом строении трехмерных пространственных моделей и многомерных триангуляций, основанный на использовании смежности граней. Этот алгоритм преобразует набор граней 3D-модели в список цепочек (list of chains), последовательно расположенных в пространстве и смежных между собой. Сжатие информации происходит за счет отсутствия дублирования номеров вершин, образующих грани модели. Описанный в статье алгоритм состоит из трех основных частей. В первой части по множеству граней модели строится специальный граф граней – ребра графа соответствуют смежным граням. Используя алгоритм обхода вершин графа в глубину, этот граф разбивается на простые цепи. Вторая часть алгоритма преобразует каждую цепь графа в последовательность номеров вершин, участвующих в формировании граней этой цепочки. Третья часть алгоритма призвана выполнять обратное действие – переводить построенную последовательность номеров вершин обратно в наборы кортежей, состоящих из номеров вершин, соответствующих граням 3D-модели. Указанный алгоритм распространен и на случай пространственных триангуляций полигональных областей. Программная реализация алгоритма для частного случая 3D-моделей выполнена в виде встраиваемых модулей в программу Blender. Архивы модулей свободно доступны в репозитории автора статьи по адресу: https://github.com/KlyachinVA/LocFile.
Бесплатно

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

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

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

Численное моделирование высотомера малых высот с двойной частотной модуляцией
Статья научная
В работе рассматривается задача повышения точности высотомера малых высот, работающего по принципу частотной модуляции. Рассмотрен вариант подавления «дискретной ошибки» высотомера путем использования дополнительной «медленной» частотной модуляции несущего колебания и усреднения результатов счета. Достоинством такого метода является простота технической реализации, выражающаяся в минимальной доработке СВЧ-тракта и регистрирующего устройства, которое должно работать в режиме усредняющего счета. Показано, что при ограничении на занимаемую Полосу частот линейная форма модуляции не является оптимальной. Приведены результаты расчета погрешности, полученные математическим моделированием работы высотомера при различных формах дополнительной частотной модуляции. Показано, что использование сложных видов «медленной» модуляции с положительной третьей производной и оптимизация в заданном диапазоне высот позволяет дополнительно в 2-3 раза снизить среднюю погрешность измерения относительно линейной формы модуляции, без расширения занимаемой полосы частот.
Бесплатно

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

Численное моделирование двумерных газодинамических течений в многокомпонентных неравновесных средах
Статья научная
Рассмотрена динамика двумерных течений многокомпонентного неравновесного газа с учетом релаксационных процессов, вязкости, теплопроводности, химических реакций, внешних источников нагрева и охлаждения. На основе численного газодинамического метода MUSCL построена вычислительная модель, позволяющая с высоким пространственным разрешением изучать нелинейные волновые структуры (слабые ударные волны, ударно-волновые импульсы, детонационные волны), возникающие в неравновесной среде из-за развития газодинамических неустойчивостей (акустической, тепловой, тангенциального разрыва скорости течения). Разработана параллельная версия численного алгоритма с использованием технологий OpenMP – CUDA – GPUDirect для гибридных вычислительных систем с несколькими графическими процессорами (multi-GPU), позволяющая существенно повысить производительность вычислений и ускорить расчеты в сотни раз по сравнению с версиями кода для CPU. Проведено численное моделирование ударно-волновых структур в плоском двумерном канале с дозвуковым потоком неравновесного газа при наличии внешних источников нагрева и охлаждения. Показано, что в результате взаимодействия дозвукового потока неравновесного акустически активного газа с твердой стенкой формируются ударно-волновые импульсы (УВИ) как с плоским фронтом в поперечном потоку направлении, так и с изгибно деформированными УВИ, на фронте которых возникают локальные максимумы газодинамических величин из-за нарастания акустически неустойчивых поперечных возмущений. Эти возмущения на нелинейной стадии эволюции образуют в канале систему косых ударных волн, которые, многократно отражаясь от стенок канала и взаимодействуя между собой, приводят к формированию сложной нерегулярной структуры течения. С увеличением ширины канала количество локальных максимумов в распределении параметров течения на фронте ударно-волновых импульсов возрастает. Исследованы нелинейные волновые структуры в плоских сверхзвуковых струях неравновесного колебательно-возбужденного газа, возникающие в результате развития неустойчивостей Кельвина — Гельмгольца и отражательных резонансных гармоник симметричных и антисимметричных мод струи. Показано, что учет колебательной неравновесности среды усиливает неустойчивость тангенциального разрыва скорости в струе и увеличивает интенсивность ударно-волновых и вихревых структур, возникающих на нелинейной стадии развития этих неустойчивостей. В численных моделях сверхзвуковых сдвиговых течений неравновесного колебательно-возбужденного газа обнаружены новые ударно-вихревые структуры высокой интенсивности, которые могут иметь важное значение в практических приложениях.
Бесплатно

Статья научная
Астрономические наблюдения показывают, что остатки сверхновых, даже с близкой к сферической формой, обычно имеют разномасштабные искажения в виде ряби или вытянутых волокон газа. Например, в остатке 0509-67.5, на ударной волне которого отчетливо видно волнистое возмущение с 15 изгибами. В качестве одного из возможных механизмов генерации таких структур рассматривается глобальная неустойчивость течения. В работе [26] в рамках линейного анализа было показано, что эта неустойчивость имеет резонансный характер, то есть возмущения определенной длины волны растут быстрее, а потому рябь в оболочке остатка проявится преимущественно в определенном диапазоне масштабов. В данной работе мы представляем результаты численного моделирования нелинейной стадии этой неустойчивости, вызванной малыми возмущениями во внешней среде в зависимости от их масштабов и интенсивности. В частности, показано, что внешние возмущения вызывают турбулизацию течения за фронтом ударной и формирование внутри оболочек радиально вытянутых филаментов с вихревой структурой,количество которых определяется номером гармоники возмущения .
Бесплатно

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