Все выпуски

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

Проводится исследование динамической эволюции шести моделей рассеянных звездных скоплений по данным о фазовых координатах звезд, полученных при численном интегрировании уравнений движения звезд. Для этой цели используются фазовые координаты звезд для 100 равноотстоящих моментов времени от начального t=0 до t_m≅5.1τ_vr (τ_vr - начальное время бурной релаксации скопления). На этом интервале времени ошибки, связанные с округлением и экспоненциальным нарастанием возмущений в исходных координатах звезд, существенно не сказываются на статистических выводах о характере движения звезд скопления. Метод исследования основан на вычислениях взаимных корреляционных функций C_1,2=C_1,2(τ,r) (τ - временная задержка, r - расстояние между точками) для флуктуаций фазовой плотности и применении Фурье-преобразования функций C_1,2 для расчета спектра частот и дисперсионных соотношений. Анализ графиков функций C_1,2, спектров частот и дисперсионных кривых подтверждает существование в моделях волн фазовой плотности, позволяет установить полный спектр радиальных колебаний фазовой плотности, отделить устойчивые колебания от неустойчивых, рассчитать периоды колебаний фазовой плотности и инкременты нарастания неустойчивых колебаний фазовой плотности. Подтверждены теоретические оценки периодов известных неустойчивых гомологических колебаний ядер моделей скоплений. Указываются некоторые астрофизические приложения полученных результатов: возникновение иррегулярных структур в рассеянных скоплениях, слабая турбулентность в движениях звезд скоплений.

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

В работе представлены результаты расчетного исследования локальной структуры восходящего газожидкостного потока в вертикальной трубе. Математическая модель основана на использовании двухжидкостного эйлерова подхода с учетом обратного влияния пузырьков на осредненные характеристики и турбулентность несущей фазы. Турбулентная кинетическая энергия жидкости рассчитывается с применением двухпараметрической изотропной модели турбулентности $k - \varepsilon$, модифицированной для двухфазных сред. Для описания динамики распределения пузырьков по размерам используются уравнения сохранения количества частиц для отдельных групп пузырьков с различными диаметрами для каждой фракции с учетом процессов дробления и коалесценции. Численно исследовано влияние изменения степени дисперсности газовой фазы, объемного расходного газосодержания, скорости дисперсной фазы на локальную структуру и поверхностное трение в двухфазном потоке. Сравнение результатов моделирования с экспериментальными данными показало, что разработанный подход позволяет адекватно описывать турбулентные газожидкостные течения в широком диапазоне изменения газосодержания и начальных размеров пузырьков.

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

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

Рассматриваются два подхода к решению задачи математического моделирования обтекания метаемых тел: численное решение уравнений движения сплошной среды Навье-Стокса, осредненных по Рейнольдсу (RANS - Reynolds-averaged Navier–Stokes), с использованием модели турбулентности и прямое численное моделирование (DNS - Direct Numerical Simulation). Тестирование рассматриваемых подходов проводится при решении задачи обтекания тел вращения с простой геометрией: сферы и цилиндра с конической головной частью, для которых известны значения коэффициентов сопротивления при различных числах Маха. Проведено качественное и количественное сравнение результатов обтекания рассматриваемых тел сверхзвуковым потоком, полученным по методикам RANS и DNS. Апробация методики численного моделирования проводится для метаемого тела (снаряда) характерной формы. Представлены результаты численного моделирования обтекания снаряда для широкого диапазона параметров: чисел Маха и углов нутации. Выполнено сравнение расчетных значений коэффициентов сопротивления с эмпирическими эталонными зависимостями по законам 1943 и 1958 годов.

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