


№ 4 (2023)
Статьи
Андрей Геннадьевич Куликовский (к 90-летию со дня рождения)



РАСЧЕТЫ СЖАТИЯ СФЕРИЧЕСКОЙ СЛОИСТОЙ СИСТЕМЫ УДАРНЫМИ ВОЛНАМИ С УЧЕТОМ ПЕРЕНОСА ТЕПЛОВОГО ИЗЛУЧЕНИЯ В КИНЕТИЧЕСКОЙ МОДЕЛИ
Аннотация
В настоящее время численное моделирование является основным, а зачастую и единственным инструментом для детального описания некоторых физических явлений при исследовании процессов сжатия веществ ударными волнами. Изучение поведения ударных волн и волн разряжения на простейших модельных тестах помогает при анализе более сложных расчетов, например, задач инерциального термоядерного синтеза на лазерных установках. В работе рассмотрены сравнительные расчеты тестовой задачи, моделирующей сжатие ударными волнами сферической слоистой системы, состоящей из двух веществ.



ДЕЙСТВИЕ ПУЛЬСИРУЮЩЕГО ИСТОЧНИКА В ЖИДКОСТИ ПРИ НАЛИЧИИ СДВИГОВОГО СЛОЯ
Аннотация
Решена двумерная нестационарная задача о развитии волнового движения в двухслойной жидкости конечной глубины, ограниченной сверху свободной поверхностью. Рассмотрены случаи, когда в невозмущенном состоянии один из слоев покоится, а в другом (приповерхностном или придонном) горизонтальная скорость потока линейно меняется по глубине. Определены дисперсионные зависимости и групповые скорости трех волновых мод, возникающих при наличии сдвигового потока. Вычислены вертикальные смещения свободной поверхности, вызванные включением пульсирующего источника, расположенного в изначально неподвижном слое жидкости. Задача рассматривается в линейной постановке, жидкость предполагается идеальной и несжимаемой.



ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ СТРУКТУРЫ ВОЗМУЩЕНИЙ ОТ ДВУХ ИМПУЛЬСНЫХ ИСТОЧНИКОВ В СВЕРХЗВУКОВОМ ПОГРАНИЧНОМ СЛОЕ ПЛАСТИНЫ
Аннотация
Описано развитие нового экспериментального метода введения в сверхзвуковой пограничный слой контролируемых возмущений с заданной частотно-волновой структурой. Представляются данные экспериментов по формированию возмущений от двух импульсных источников (импульсный тлеющий разряд) в ламинарном пограничном слое пластины при числе Маха, равном 2. Эксперименты выполнены в аэродинамической трубе Т-325 ИТПМ СО РАН. Локализованные источники располагались на одинаковом расстоянии от передней кромки пластины в 6 мм друг от друга по размаху. Пульсации потока измерялись с помощью однониточного датчика термоанемометра постоянного сопротивления, запись сигнала проводилась синхронно с зажиганием разрядов, что позволяло выделять возмущения от разрядов из фона случайных неконтролируемых “естественных” пульсаций пограничного слоя. Анализируются пространственно-временная структура и частотно-волновой состав генерируемых возмущений от одиночного и двух разрядов, работающих синхронно и с задержкой во времени. Получено, что наибольшие отличия в структуре возмущений от одного и двух источников наблюдаются в центральной области, тогда как на боковых границах возмущения пульсации близки во всех рассматриваемых случаях. В спектрах возмущений по поперечным волновым числам от двух разрядов формируются узлы и пучности, положение которых определяется расстоянием между источниками и временной задержкой в их работе.



РЕЗОНАНСНОЕ ВОЗДЕЙСТВИЕ ПЕРИОДИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ ПЛАЗМЕННЫХ АКТУАТОРОВ ПРИ УПРАВЛЕНИИ ТЕЧЕНИЕМ В ПОГРАНИЧНОМ СЛОЕ НА СТРЕЛОВИДНОМ КРЫЛЕ
Аннотация
Теоретически изучен процесс управления течением в пограничном слое на стреловидном крыле с помощью ряда плазменных актуаторов, установленных под углом к передней кромке. Исследован резонансный отклик управляемого течения в случае, когда угол наклона актуаторов близок углу наклона стационарных вихрей неустойчивости поперечного течения. Расчет управляемого течения произведен в рамках параболизованной системы уравнений Навье–Стокса. Вычислены характеристики устойчивости периодических течений, сформированных воздействием управляющей объемной силы.



СТРУКТУРА ТЕЧЕНИЯ И ТЕПЛООБМЕН В СТАЦИОНАРНОМ ГАЗОКАПЕЛЬНОМ ПОТОКЕ ЗА ТОЧКОЙ ПАДЕНИЯ УДАРНОЙ ВОЛНЫ НА ПЛОСКУЮ СТЕНКУ
Аннотация
Рассматриваются задачи о структуре стационарного двумерного течения газокапельной смеси в пристеночной области за точкой падения косой и прямой волн уплотнения на плоскую стенку. В данном случае прямая волна соответствует ножке Маха при маховском отражении, а косая – регулярному режиму отражения падающей наклонной волны. Основной целью исследования является оценка влияния мелких капель жидкости, присутствующих в набегающем потоке, на равновесную температуру адиабатической стенки за точкой падения волны. Исследуется вопрос – насколько наличие падающей на стенку косой ударной волны может усилить эффект снижения равновесной температуры стенки присутствующими в потоке мелкими каплями. Область течения разбивается на внешнюю область “эффективно невязкого течения” и область асимптотического ламинарного пограничного слоя. Расчеты течения в каждой из областей проводятся в рамках двухконтинуальной модели газокапельной смеси с учетом фазового перехода (испарения) на поверхности капель. Исследованы наиболее интересные с точки зрения теплообмена волновые конфигурации, соответствующие неполному испарению капель за отраженной волной, а также “волнам с частичной и полной дисперсией” параметров. Используется упрощенная предельная схема формирования жидкой пленки оседающими на стенку каплями, эффектами неустойчивости и разбрызгивания пленки пренебрегается. На основании численных расчетов получены оценки возможного снижения равновесной температуры адиабатической стенки за точкой падения ударных волн в стационарном сверхзвуковом потоке газа, содержащем малую концентрацию примеси жидких капель.



ОБТЕКАНИЕ ТЕЛ ЗАПЫЛЕННЫМ ГАЗОМ ПРИ РАССЕЯНИИ ОТРАЖЕННЫХ ЧАСТИЦ
Аннотация
Рассмотрены обтекание плоской пластины конечной толщины в канале с большой дозвуковой скоростью и сверхзвуковое поперечное обтекание цилиндра двухфазным потоком газа с твердыми частицами. Передняя кромка пластины имеет форму клина или гладкое затупление постоянного радиуса. Поверхность клина и переднего затупления задается гладкой или шероховатой. Шероховатость моделируется двумерным профилем, который задается на основе эксперимента. Рассмотрены сферические частицы и смесь частиц в виде эллипсоидов вращения, прямоугольных призм, призм со срезанными вершинами и тетраэдров. Параметры каждой из форм варьируются. При определении поступательной и вращательной скоростей несферических частиц после отскока используется модель ударного взаимодействия, предложенная ранее и согласующаяся с экспериментальными данными по коэффициентам восстановления скорости центра масс. Наряду с монодисперсной примесью рассмотрена дисперсная фаза с разбросом частиц по размерам. На основе анализа численных результатов установлена роль исследованных факторов случайной природы на картину течения и параметры примеси.



НЕИЗОТЕРМИЧЕСКОЕ ТЕЧЕНИЕ ЖИДКОСТИ В СКВАЖИНЕ ПРИ ИНДУКЦИОННОМ НАГРЕВЕ ОБСАДНОЙ КОЛОННЫ
Аннотация
Исследованы особенности полей скорости течения и температуры восходящего потока жидкости в металлической круглой трубе (обсадной колонне, установленной в добывающей скважине) в условиях ее локального индукционного нагрева. Результаты исследований основаны на численном решении уравнений Навье–Стокса в приближении Буссинеска–Обербека, расчеты выполнены в программном пакете Ansys Fluent (Лицензия ANSYS Academic Research CFD в рамках договора с Башкирским государственным университетом от 15.06.2020). Рассмотрены расходы жидкости 10 и 50 кубических метров в сутки, соответствующие ламинарному и переходному режимам течения в обсадной трубе. Установлено наличие локальных возмущений поля скорости и температуры в пристеночной области нагретой обсадной трубы. Возмущения температуры в жидкости достигают нескольких градусов Кельвина, причем локальная скорость потока в пристеночной области обсадной колонны, возрастающая за счет естественной тепловой конвекции, в несколько раз превышает среднюю по сечению скорость потока. Показано возникновение областей вихревого движения потока над интервалом индукционного нагрева, обусловленного естественной тепловой конвекцией.



ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ СТРУКТУРЫ СИСТЕМЫ ТРЕЩИН НА ФИЛЬТРАЦИЮ ЖИДКОСТИ В ПОРОУПРУГОЙ СРЕДЕ
Аннотация
Рассматривается двумерная однофазная фильтрация слабосжимаемой жидкости в деформируемой трещиновато-пористой среде. Для совместного моделирования процессов фильтрации и связанного с ними изменения напряженно-деформированного состояния среды использована модель пороупругой среды, моделирование трещиноватости выполнено с помощью модели дискретных трещин. Трещины в рассматриваемой области имели случайное положение и ориентацию, распределение трещин по длинам подчинялось степенному закону. Исследовалась зависимость фильтрационных свойств трещиновато-пористой среды от ее напряженно-деформированного состояния и структуры системы трещин. Численное исследование выполнено для вариантов систем трещин, полученных путем множественной случайной генерации. Установлено, что фильтрационные свойства трещиновато-пористой среды определяются главным образом структурой системы трещин, характеризуемой параметром перколяции. Показано, что существенное влияние напряженно-деформированного состояния среды на ее фильтрационные свойства наблюдается для только связных систем трещин. Предложена формула для аппроксимации зависимости эквивалентной проницаемости трещиновато-пористой среды от параметров характеризующих связность системы трещин, напряженно-деформированное состояние среды, деформационные и фильтрационные свойства трещин.



ИССЛЕДОВАНИЕ МЕХАНИЗМА ВЗАИМОДЕЙСТВИЯ УГЛЕВОДОРОДНОГО ПЛАМЕНИ С ЭЛЕКТРИЧЕСКИМ ПОЛЕМ
Аннотация
Рассматривается один из возможных механизмов взаимодействия ограниченной области горения (факела пламени) с приложенным электрическим полем. Исследование основано на применении методов электрогидродинамики (ЭГД) для описания химически реагирующих многокомпонентных неравновесных частично ионизованных газовых смесей. Показано, что в окрестности границ пламени при определенных условиях образуются зоны объемного электрического заряда, на которые можно воздействовать полем. Исследуется характер изменения этих зон под влиянием поля.



ФОРМИРОВАНИЕ ПРОСТРАНСТВЕННЫХ ВНУТРЕННИХ ВОЛН ЗА ТЕЛОМ, ДВИГАЮЩИМСЯ В СТРАТИФИЦИРОВАННОЙ ВЯЗКОЙ ЖИДКОСТИ
Аннотация
Рассмотрено равномерное движение диска в горизонтальном направлении вдоль оси его симметрии в покоящейся стратифицированной вязкой жидкости. Диск генерирует пространственные гравитационные внутренние волны, занимающие все пространство между диском и местом его старта. Волны наблюдаются при помощи двухцветной бета-плюс-визуализации вихревой структуры течения, рассчитанного при помощи системы уравнений Навье–Стокса в приближении Буссинеска. В настоящей работе существенно дополнен опубликованный ранее механизм формирования полуволн над осью симметрии диска, где основное внимание уделялось периодическому процессу зарождения деформированных вихревых колец над местом старта диска, происходящему в силу гравитационной и сдвиговой неустойчивостей; левое полукольцо трансформируется в полуволну впадин или гребней, а правое полукольцо исчезает со временем. В настоящей работе установлено, что левые части правых нечетных полуколец превращаются в осевые части полуволн гребней.



ОСОБЕННОСТИ ФОРМИРОВАНИЯ ФАКЕЛА ПРИ СТОЛКНОВЕНИИ ДВУХ ЛАМИНАРНЫХ ГАЗОВЫХ СТРУЙ
Аннотация
Представлены результаты исследований процесса взаимодействия двух соударяющихся осесимметричных ламинарных микроструй пропан/бутановой смеси при диффузионном горении и без горения. Истечение газовой смеси осуществлялось через круглые трубки с равными скоростями. В процессе эксперимента варьировалось поперечное положение трубок относительно друг друга с сохранением угла между ними. Выявлены особенности формирования результирующей струи в зависимости от поперечного положения трубок. Если трубки находятся в одной плоскости, то результирующая струя формируется в ортогональной к ней плоскости. Данный процесс наблюдается при взаимодействии горящих и не горящих струй. При увеличении скорости истечения струй обнаружена область локального разрыва фронта пламени.



ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ СВЕРХЗВУКОВОГО ПОТОКА С ОБЛАСТЬЮ ТЕПЛОВЫДЕЛЕНИЯ ПРОДОЛЬНО-ПОПЕРЕЧНЫМ РАЗРЯДОМ
Аннотация
Проведено численное моделирование сверхзвукового стационарного воздушного потока внутри расширяющегося аэродинамического канала с прямоугольным сечением – лабораторной модели прямоточного воздушно-реактивного двигателя. С помощью экспериментальных данных была проведена валидация аэродинамической модели в случае отсутствия зоны объемного тепловыделения. После валидации модели было проведено численное моделирование сверхзвукового потока с включенной зоной объемного тепловыделения. Получены трехмерные распределения скорости, температуры и давления в сверхзвуковом стационарном воздушном потоке. Показано, что при объемной плотности тепловой мощности источника эквивалентной средней общей мощности разряда W = 10 кВт разряд нагревает газ до температуры T = 1700–4200 К, что приводит к ускорению потока без его теплового запирания. При плотности тепловой мощности источника эквивалентной средней общей мощности разряда W = 20 кВт газ нагревается сильнее до 6700 К, но начинается локальное тепловое запирание потока.



ИССЛЕДОВАНИЕ РЕЖИМОВ ТЕПЛООБМЕНА В ДОЗВУКОВЫХ СТРУЯХ ДИССОЦИИРОВАННОГО АЗОТА ВЫСОКОЧАСТОТНОГО ИНДУКЦИОННОГО ПЛАЗМОТРОНА ПРИ ДОПОЛНИТЕЛЬНОМ НАГРЕВЕ ПОВЕРХНОСТИ ЛАЗЕРНЫМ ИЗЛУЧЕНИЕМ
Аннотация
На индукционном ВЧ-плазмотроне ВГУ-4 (ИПМех РАН) проведены эксперименты по исследованию теплообмена цилиндрической водоохлаждаемой медной модели диаметром 30 мм, оснащенной калориметрическим датчиком с тепловоспринимающей поверхностью диаметром 13.8 мм из графита, при комбинированном режиме нагрева поверхности плазмой азота и лазерным излучением, а также для случаев нагрева поверхности только лазерным излучением или только струей плазмы азота. Эксперименты в струях ВЧ-плазмотрона проведены при давлении в барокамере установки p = 1 × 104 Па, массовом расходе азота G = 2.4 г/c, мощности ВЧ-генератора плазмотрона по анодному питанию Na.p. = 22 кВт. Для выбранных экспериментальных режимов установлено, что струя диссоциированного азота и высокочастотный индукционный разряд не оказывают заметного влияния на проходящий через них лазерный луч. Получены значения плотности теплового потока в зависимости от подведенной мощности лазерного излучения. Для рассмотренных условий экспериментов проведено численное моделирование дозвукового течения плазмы азота в кварцевом разрядном канале и в пространстве барокамеры установки ВГУ-4, основанное на решении полных уравнений Навье–Стокса методом Патанкара–Сполдинга.



АНАЛИЗ АВТОКОЛЕБАТЕЛЬНЫХ ПРОЦЕССОВ В КАВЕРНЕ С ОТКРЫТЫМ ТИПОМ ТЕЧЕНИЯ НА ОСНОВЕ ДАННЫХ ВИХРЕРАЗРЕШАЮЩИХ РАСЧЕТОВ
Аннотация
Рассматриваются механизмы автоколебательных процессов, возникающих в каверне с открытым типом течения, и дается их обоснование, полученное в результате детального изучения взаимодействия явлений гидродинамической, расходно-волновой и резонансной природы. Теоретические выводы подкрепляются анализом данных вычислительных экспериментов, проведенных различными авторами.


