Все выпуски

На основе гидродинамического 3D-моделирования с использованием уравнения состояния газа твердых сфер Карнахана – Старлинга выполнено исследование развития периодических возмущений сходящейся сферической ударной волны, приводящих к ограничению кумуляции. Метод решения системы уравнений Эйлера на подвижных (сжимающихся) сетках позволяет с высокой точностью проследить эволюцию фронта сходящейся ударной волны в широком диапазоне изменения ее радиуса. Скорость сжатия расчетной сетки адаптируется к движению фронта ударной волны, при этом движение границ расчетной области выбирается из условия сверхзвуковой скорости ее движения относительно среды. Это приводит к тому, что решение на этапе сжатия определяется только начальными данными. Применена схема TVD второго порядка аппроксимации для реконструкции вектора консервативных переменных на границах расчетных ячеек в сочетании со схемой Русанова для расчета численного вектора потоков. Выбор обусловлен сильной тенденцией к проявлению в расчетах численной неустойчивости типа «карбункул», известной для других классов течений. Использование сжимающихся сеток позволило исследовать детальную картину течения на масштабе прекращения кумуляции, что невозможно в рамках метода геометрической динамики ударных волн Уизема (Whitham), применявшегося ранее другими авторами для расчета сходящихся ударных волн. Исследование показало, что ограничение кумуляции связанно с переходом от маховского взаимодействия сегментов сходящейся ударной волны к регулярному вследствие прогрессирующего роста отношения азимутальной скорости на фронте ударной волны к радиальной при уменьшении ее радиуса. Установлено, что это отношение представляется в виде произведения ограниченной осциллирующей функции радиуса и степенной функции радиуса с показателем степени, зависящим от начальной плотности упаковки в модели твердых сфер. Показано, что увеличение параметра плотности упаковки в модели твердых сфер приводит к значительному увеличению давлений, достигаемых в ударной волне с нарушенной симметрией. Впервые в расчете показано, что на масштабе прекращения кумуляции течение сопровождается формированием высокоэнергетичных вихрей, в которые вовлечено вещество, подвергшееся наибольшему ударно-волновому сжатию. Оказывая влияние на процессы тепло- и массопереноса в области наибольшего сжатия, это обстоятельство является важным для актуальных практических применений сходящихся ударных волн в целях инициирования реакций (детонации, фазовых переходов, управляемого термоядерного синтеза).

В данной работе рассматривается взаимодействие кинка уравнения $\varphi^4$ с двумя протяженными одинаковыми примесями. Протяженная примесь описывается с помощью функции прямоугольного вида. Анализируется случай притягивающей примеси. С помощью аналитических методов рассматривается случай малых амплитуд локализованных волн, когда возможно провести линеаризацию уравнений движения. Для численного решения использовался метод прямых для уравнений в частных производных. Для нахождения частот колебаний, локализованных на примесях волн, используется дискретное преобразование Фурье. Кинк запускался в направлении примесей с разными начальными скоростями. Изменялось также расстояние между двумя примесями. Показано, что при взаимодействии кинка с примесями на них возбуждаются долгоживущие локализованные волны бризерного типа. Исследована их структура и связанная динамика. Определено, как, изменяя параметры примесей и расстояние между ними, можно управлять типом и динамическими параметрами связанных колебаний, локализованных на примесях волн. Найдены возможные решения в виде синфазных, антифазных колебаний, в виде биений. Колебания локализованных волн происходят с излучением волн малой амплитуды. Спектр этих излучений состоит из двух частот. Первая приближенно равна $\sqrt{2}$, что соответствует величине частоты для хвоста воблингбризера уравнения $\varphi^4$. Вторая приближенно равна удвоенной частоте колебаний примесных мод. Найдено (как аналитически, так и численно) наличие двух возможных частот для связанных локализованных колебаний. Показано, что частоты сильно зависят от расстояния между примесями. С увеличением расстояния между примесями частоты сливаются в одну — частоту, полученную для случая одиночной примеси. Найденные численно и аналитически зависимости частот от расстояния между примесями хорошо совпадают для больших расстояний, когда взаимодействие между примесями слабое, и начинают заметно отличаться при малых расстояниях, когда взаимодействие между примесями сильное. Аналитическое значение величин полученных частот всегда больше численных. Показано, что зависимость амплитуды локализованных волн от начальной скорости кинка имеет несколько минимумов и максимумов.

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

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

The paper is an example of computer simulation in mechanical engineering. Velocity field in a hydrocyclone are determined numerically, because for direct measurements it is difficult to achieve them. The numerical simulation of 3D fluid dynamics based on the k-eps RNG model of turbulence in the hydrocyclone with the injector, containing 5 tangentially directed nozzles is considered. It is shown that the direction of movement of
injected fluid in the hydrocyclone depends on the water flow rate through the injector. The calculations show in accordance with the experiments that the dependence of the Split-parameter on the injected water flow rate has a non-monotone character associated with the ratio of power of the main flow and the injected fluid.

В данной работе решается задача поиска конструкции магнитной системы для создания магнитного поля с требуемыми характеристиками в заданной области. На основе анализа математической модели магнитной системы предлагается достаточно общий подход к решению нелинейной обратной задачи, которая описывается уравнением Фредгольма H(z) = ∫_SIJ(s)G(z, s)ds, z ∈ S _H, s ∈ S _I . Необходимо определить распределение плотности тока J(s), а также расстановку источников тока для создания поля H(z). В работе предлагается метод решения этих задачс помощью регуляризованных итерационных процессов. На примере конкретной магнитной системы проводится численное исследование влияния различных факторов на характер создаваемого магнитного поля.

Проведено математическое моделирование нестационарных режимов естественной конвекции в замкнутой пористой цилиндрической полости с теплопроводной оболочкой конечной толщины в условиях конвективного теплообмена с внешней средой. Краевая задача математической физики, сформулированная на основе модели Дарси–Буссинеска в безразмерных переменных «функция тока – температура», реализована численно методом конечных разностей. Детально проанализировано влияние проницаемости пористой среды 10^–5≤Da<∞, отношения толщины твердой оболочки к внутреннему радиусу цилиндра 0.1≤h/L≤0.3, относительного коэффициента теплопроводности 1≤λ_1,2≤20 и безразмерного времени 0≤τ≤1000 как на локальные распределения изолиний функции тока и температуры, так и на интегральные комплексы, отражающие интенсивность конвективного течения и теплопереноса.

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

В статье приводятся результаты моделирования деформации метеорита при попадании в плотные слои атмосферы разрывным методом Галёркина на неструктурированных треугольных сетках и методом сглаженных частиц. В качестве исходных данных брались материалы по Челябинскому метеориту. Проводилась серия расчётов, где варьировались характеристики материала метеорита и угол входа в плотные слои атмосферы.

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