Часто просматриваемые статьи

111. Говорухин В.Н., Филимонова А.М.
     Расчет плоских геофизических течений невязкой несжимаемой жидкости бессеточно-спектральным методом
     Компьютерные исследования и моделирование, 2019, т. 11, № 3, с. 413-426

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

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

     Ключевые слова: идеальная жидкость, вихревые методы, вихревые структуры, геофизические течения.
     Просмотров за год: 16.
112. Мадера А.Г., Кандалов П.И.
     Математическое моделирование интервально стохастических тепловых процессов в технических системах при интервальной неопределенности определяющих параметров
     Компьютерные исследования и моделирование, 2016, т. 8, № 3, с. 501-520

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

     Ключевые слова: математическое моделирование, тепловой процесс, техническая система, интервальный, стохастический, нелинейный, нестационарный, статистические меры, математическое ожидание, дисперсия, ковариация.
     Просмотров за год: 15. Цитирований: 6 (РИНЦ).
113. Борисов А.В., Розенблат Г.М.
     Новые алгоритмы составления дифференциальных уравнений движения экзоскелета с переменной длиной звеньев и управлением в шарнирах-суставах
     Компьютерные исследования и моделирование, 2017, т. 9, № 2, с. 201-210

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

     Ключевые слова: экзоскелет, звено переменной длины, шарнир-сустав, управление, дифференциальные уравнения движения, рекурсия, матрица, метод, интеграл энергии, численное решение.
     Просмотров за год: 15. Цитирований: 2 (РИНЦ).
114. Фомин А.А., Фомина Л.Н.
     О сходимости неявного итерационного полинейного рекуррентного метода решения систем разностных эллиптических уравнений
     Компьютерные исследования и моделирование, 2017, т. 9, № 6, с. 857-880

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

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

     Ключевые слова: система линейных алгебраических уравнений, итерационный метод решения, сходимость метода.
     Просмотров за год: 15. Цитирований: 1 (РИНЦ).
115. Белотелов Н.В., Апалькова Т.Г., Мамкин В.В., Курбатова Ю.А., Ольчев А.В.
     Некоторые особенности взаимосвязи термодинамических характеристик земной поверхности с потоками водяного пара и диоксида углерода на сплошной свежей вырубке
     Компьютерные исследования и моделирование, 2017, т. 9, № 6, с. 965-980

     В статье рассматриваются некоторые особенности временной изменчивости и взаимосвязь эксергии коротковолнового и длинноволнового излучения с потоками тепла, водяного пара (Н₂О) и диоксида углерода (СО₂) на сплошной свежей вырубке в смешанном хвойно-мелколиственном лесу в Тверской области. На основе проведенного сравнения коэффициентов радиационной эффективности и эффективности эксергии, введенных Ю. М. Свирежевым, было показано, что в течение первых восьми месяцев после проведения сплошной рубки лесная экосистема функционирует как «тепловая машина», то есть на вырубке доминируют процессы диссипации энергии над продукционными процессами. Для проверки справедливости полученного вывода был выполнен статистический анализ временной изменчивости метеорологических рядов наблюдений и данных по среднесуточным значениям потоков явного тепла, затрат тепла на испарение и потоков СО₂ с применением тригонометрических многочленов, который подтвердил полученный ранее вывод. Для среднесуточных значений потоков СО₂, валовой первичной продуктивности растительного покрова на вырубке, а также потоков явного тепла удалось построить статистические модели, линейно зависящие от эксергии коротковолнового и длинноволнового излучения. Анализ этих зависимостей также подтвердил вывод, полученный на основе рассмотрения временной изменчивости коэффициентов радиационной эффективности и эффективности эксергии. Используя разбиение временных рядов на отдельные временные интервалы «весна–лето» и «лето–осень», удалось выявить, что в процессе зарастания вырубки травянистой растительностью в летние месяцы связь между этими параметрами и величиной эксергии усиливается. Анализ линейной связи временных рядов затрат тепла на испарение и эксергии показал ее статистическую незначимость. В свою очередь, линейная связь между затратами тепла на испарение и температурой оказалась статистически значимой. Температура выступала в роли ключевого фактора, повышающего точность модели, а эксергия оказывалась несущественным фактором. Это указывает на то, что межсуточная временная изменчивость испарения активно зарастающей вырубки определялась главным образом температурой.

     Ключевые слова: сплошная вырубка, термодинамические характеристики, эксергия, коротковолновая радиация, длинноволновое излучение, временные ряды, потоки тепла, водяного пара и диоксида углерода.
     Просмотров за год: 15. Цитирований: 1 (РИНЦ).
116. Зыза А.В.
     Компьютерное исследование полиномиальных решений уравнений динамики гиростата
     Компьютерные исследования и моделирование, 2018, т. 10, № 1, с. 7-25

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

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

     Третье решение построено для уравнений движения гиростата в магнитном поле с учетом эффекта Барнетта–Лондона. Для него структура такова: первая и вторая компоненты вектора угловой скорости—многочлены второй степени, квадрат третьей компоненты—многочлен четвертой степени по вспомогательной переменной, которая находится обращением эллиптического интеграла Лежандра.

     Все построенные решения не имеют аналогов в динамике твердого тела с неподвижной точкой.

     Ключевые слова: полиномиальные решения, уравнения Кирхгофа–Пуассона, гиростат, потенциальные и гироскопические силы, эффект Барнетта–Лондона, эллиптические интегралы Лежандра.
     Просмотров за год: 15.
117. Леонов А.В., Колтовская Е.В., Чичерина О.В.
     Биогидрохимический портрет Белого моря
     Компьютерные исследования и моделирование, 2018, т. 10, № 1, с. 125-160

     Биогидрохимический портрет Белого моря построен с помощью расчетов на CNPSi-модели по систематизированным среднемноголетним наблюдениям (среднемесячные гидрометеорологические, гидрохимические и гидробиологические параметры морской среды). Также в расчетах использована уточненная информация о выносе в морские акватории биогенных веществ со стоком основных рекритоков (Нива, Онега, Северная Двина, Мезень, Кемь, Кереть). Параметры морской среды — значения температуры, освещенности, прозрачности, биогенной нагрузки. Для девяти районов моря (заливы Кандалакшский, Онежский, Двинский, Мезенский, Соловецкие о-ва, Бассейн, Горло, Воронка, губа Чупа) характеристики портрета моря включают: изменение в течение года концентраций органических и минеральных соединений биогенных элементов (С, N, P, Si), биомассы организмов низших трофических звеньев (гетеротрофные бактерии, диатомовый фитопланктон, растительноядный и хищный зоопланктон) и другие показатели (скорости изменения концентраций веществ и биомасс организмов, внутренние и внешние потоки веществ, балансы отдельных веществ и биогенных элементов в целом). Расчетные по среднемноголетним данным показатели состояния морской среды (температура воды, соотношения минеральных фракций N < P) и доминирующего диатомового фитопланктона в море (обилие, продукция, биомасса, содержание хлорофилла а) сравнивали с результатами отдельных съемок (за 1972–1991 и 2007–2012 гг.) по районам моря. При очевидных отличиях способов оценки значений показателей (по наблюдениям — аналитические методы, а при расчетах на модели — вычисления по соответствующим уравнениям) отмечена близость расчетных показателей состояния фитопланктона приведенным в литературе данным по фитопланктону Белого моря. Так, литературные оценки годовой продукции диатомовых водорослей в Белом море находятся в пределах 1.5–3 млн т С (при продолжительности вегетации 180 сут), а по расчетам она составляет ~2 и 3.5 млн т С при принимаемых периодах вегетации в 150 и 180 сут соответственно.

     Ключевые слова: экосистема Белого моря, биогенные вещества (БВ), гетеротрофный бактериопланктон, диатомовый фитопланктон, растительноядный и хищный зоопланктон, детрит, трофическая цепь, CNPSi-модель биотрансформации БВ, экологический портрет Белого моря, сравнение наблюдаемых и расчетных показателей диатомовых водорослей (обилие, продукция, биомасса, хлорофилл а).
     Просмотров за год: 15. Цитирований: 1 (РИНЦ).
118. Шепелев В.В., Фортова С.В., Опарина Е.И.
     Использование программного комплекса Turbulence Problem Solver (TPS) для численного моделирования взаимодействия лазерного излучения с металлами
     Компьютерные исследования и моделирование, 2018, т. 10, № 5, с. 619-630

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

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

     Ключевые слова: пакет прикладных программ, уравнения газовой динамики, фемтосекундные лазеры, уравнения состояния.
     Просмотров за год: 15.
119. Акопов А.С., Бекларян Л.А., Бекларян А.Л., Сагателян А.К.
     Укрупненная модель эколого-экономической системы на примере Республики Армения
     Компьютерные исследования и моделирование, 2014, т. 6, № 4, с. 621-631

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

     Ключевые слова: экологическое моделирование, системная динамика, многопараметрическая оптимизация, имитационное моделирование.
     Просмотров за год: 14. Цитирований: 7 (РИНЦ).
120. Башкирцева И.А., Бояршинова П.В., Рязанова Т.В., Ряшко Л.Б.
     Анализ индуцированного шумом разрушения режимов сосуществования в популяционной системе «хищник–жертва»
     Компьютерные исследования и моделирование, 2016, т. 8, № 4, с. 647-660

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

     Ключевые слова: популяционная динамика, случайные возмущения, функция стохастической чувствительности, доверительные области.
     Просмотров за год: 14. Цитирований: 4 (РИНЦ).

• Общая информация
• Политика открытого доступа
• Публикационная этика
• Редакционная коллегия
• Авторский указатель
• Для авторов
• Для рецензентов
• Контакты
• Часто просматриваемые статьи
• Наиболее цитируемые статьи
• Подача статей ONLINE