


Vol 62, No 4 (2024)
Articles
К юбилею Игоря Владимировича Ломоносова



Plasma Investigations
Структурный переход в сильнонеидеальных кулоновских кластерах
Abstract
Методом молекулярной динамики исследованы равновесные состояния кулоновских кластеров, содержащих до 5000 частиц, при высоких значениях кулоновского параметра неидеальности. Показано, что при числе частиц более 2000 в кулоновских кластерах появляется кристаллизованное ядро с доминирующей структурой ГПУ, которое плавится при величине параметра неидеальности, близкой к характерной для бесконечной однокомпонентной плазмы. В целом большой кластер представляет собой кристаллическое ядро, окруженное сферическими оболочками, число которых не зависит от полного числа частиц, но зависит от параметра неидеальности. Большая величина порога структурного перехода объясняется близостью потенциальных энергий кристаллической структуры и системы вложенных сферических оболочек. Фактор сжимаемости частиц, определенный как по вириалу сил, так и по энергии, оказывается близким к нулю в широком диапазоне чисел частиц, что позволяет экстраполировать этот результат к классической однокомпонентной плазме. Предложено объяснение данного результата, основанное на модели ячеек Вигнера–Зейтца.



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



Thermophysical Properties of Materials
Теплопроводность жидких сплавов системы Na–Pb
Abstract
Методом лазерной вспышки измерена теплопроводность λ жидких сплавов натрий–свинец (10, 21, 31, 41, 50 и 63 ат. % Pb) в интервале температуры от линии ликвидуса до 1070 К с неопределенностью 4–6%. Построены температурные и концентрационные зависимости λ. Установлено, что для большинства исследованных сплавов теплопроводность с температурой растет монотонно и имеет довольно низкие значения по сравнению с чистыми расплавами Na и Pb: от 3.7 до 12.0 Вт/ (мК). На концентрационной зависимости λ выявлен пологий минимум в интервале концентрации свинца XPb≈ 20–50 ат. %. Отмечены корреляции между полученными результатами по теплопроводности и другими теплофизическими свойствами, которые косвенно подтверждают существующие в литературе представления о формировании в жидкой системе Na–Pb ионных комплексов.



Расчет фазовой диаграммы индия на основе законов подобия, связанных с линией единичного фактора сжимаемости
Abstract
Кривая сосуществования жидкость–газ (бинодаль) для индия на плоскости плотность–температура рассчитана с помощью законов подобия, связанных с линией единичного фактора сжимаемости, и численного моделирования методом погруженного атома. Для этого использовалась ранее разработанная модель, учитывающая общие асимптотические свойства бинодали, следующие из теории критических явлений и используемых законов подобия. Эта модель ранее была проверена для веществ с известной бинодалью и успешно применялась для ряда металлов, для которых положение бинодали и критических точек не может быть пока получено из измерений и более строгих методов расчета. Кроме этого, методом погруженного атома рассчитана плотность жидкости на линии кристаллизации индия при известных из измерений давлениях и температурах. Данные об этой плотности до сих пор отсутствовали в литературе.



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



Влияние процессов закалки и отжига на высокотемпературные кинетические свойства сплавов NiV
Abstract
Приведены результаты экспериментального исследования температуропроводности и удельного электрического сопротивления твердых растворов и интерметаллических соединений сплавов системы NiV. Эксперимент показал, что процесс закалки–отжига наиболее сильно влияет на свойства упорядоченных интерметаллических соединений на основе Ni2V, Ni3V и NiV3. Объяснение полученных результатов дано на основании теории явлений переноса, рассматривающей два механизма рассеяния носителей: примесный и фононный.



Heat and Mass Transfer and Physical Gasdynamics
Местное ускорение ламинарного потока и аномальная интенсификация теплообмена в канале с двумя рядами наклонных канавок
Abstract
Проведено экспериментальное и численное исследование аномальной интенсификации ламинарного отрывного течения и теплообмена в канале с двумя рядами из 26 плотно расположенных канавок под углом наклона ±135° при равномерном потоке на входе и изменении чисел Re от 1000 до 5500. Обоснованы локальное ускорение потока с достижением максимальной скорости порядка 1.8 и утончение пристеночного слоя с ростом продольной скорости, доходящей до величины 1.4 на расстоянии y = 0.005 от стенки при Re = 2500, над входными сферическими сегментами канавок. Установлена корреляция ускорения и аномальной интенсификации ламинарного отрывного течения и теплообмена. Отмечено, что при Re = 2500 минимальная величина относительного отрицательного трения доходит до –7, а относительная теплоотдача от структурированного участка канала при Re = 5500 достигает 4.6.



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



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



О некоторых особенностях осаждения одиночных капель на модель с полусферическим торцем при различных температурах
Abstract
В работе предпринята попытка экспериментального изучения некоторых эффектов, наблюдаемых при гравитационном осаждении одной или двух капель на модель с полусферическим торцем при различных температурах. Путем покадрового анализа видеозаписей проанализированы особенности взаимодействия капель с поверхностью модели. Эксперименты показали, что в случае комнатных температур капли (Td = 25–27 °C) и поверхности модели (TW = 25–27 °C) процессы растекания и собирания второй капли проходят не на поверхности модели, а на поверхности первой капли, и сопровождаются повторным растеканием первой капли и ее повторным собиранием. В случае пониженных температур капли (Td = 7–9 °C) и поверхности (TW = –3...+0.5 °C) процесс собирания первой капли отсутствовал, а процесс растекания второй капли отличался от растекания первой капли. Полученные результаты призваны улучшить понимание механизма взаимодействия между каплей и поверхностью в процессе ударного замораживания и совершенствовать технологии защиты от обледенения.



Спектральная зависимость оптической плотности среды в УФ и видимой областях при образовании сажи и полиароматических углеводородов при пиролизе углеводородов
Abstract
В работе измерены температурные зависимости оптической плотности, отражающей образование полиароматических углеводородов и конденсированных частиц сажи при ударно-волновом пиролизе смесей метана, этилена, ацетилена и бензола с аргоном, на длинах волн 313, 405 и 633 нм в диапазонах температур 1500–2600 К и давлений 3.5–6.4 бар. Произведена оценка вклада конденсированной фазы в поглощение на длинах волн 313 и 405 нм. Определены температурные зоны наиболее интенсивного формирования полиароматических углеводородов и сажи при пиролизе исследуемых смесей. Наибольший выход полиароматических углеводородов и сажи наблюдался при пиролизе бензола, наименьший – при пиролизе метана и этилена. Отмечено, что максимум выхода полиароматических углеводородов и сажи в случае пиролиза ацетилена существенно сдвинут по температуре относительно других углеводородов, что свидетельствует об отличии кинетического механизма сажеобразования в случае пиролиза ацетилена.



Особенности тепло- и массообмена при горении гранулированной смеси Zr + 0.5C в спутном потоке аргона
Abstract
В работе впервые исследованы закономерности синтеза керамики ZrC0.5 из гранулированной смеси в спутном потоке аргона. В предварительных экспериментах без продува аргоном показано, что примесный газ не влияет на скорость горения смесей Zr + 0.5C, а передачу горения от гранулы к грануле обеспечивает движение расплава циркония. В спутном потоке аргона обнаружен новый режим горения гранулированных смесей, характеризующийся высокой скоростью распространения фронта горения и сопровождающийся уменьшением как продольного, так и поперечного размеров сгоревшего образца с образованием так называемого фингера. Совокупность наблюдаемых явлений можно объяснить в предположении, что из-за значительной усадки образца в продольном направлении проницаемость продуктов уменьшается, что препятствует фильтрации аргона через продукты, направляя поток газа к боковой поверхности цилиндрического реактора в образующийся фингер. Увеличение потока газа в узком канале вдоль боковой поверхности реактора в свою очередь вызывает увеличение видимой скорости и повышение температуры горения смеси на периферии образца. Проведенные эксперименты и расчеты показывают, что уменьшение размеров образца в продольном направлении обеспечивается перепадом давления фильтрующегося газа, а в поперечном сечении – действием сил поверхностного натяжения из-за градиента температуры.



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



New energy and modern technologies
Современное состояние углеродной энергетики
Abstract
На основе данных научных центров и результатов долгосрочных научных проектов проанализированы тепловое состояние Земли и глобальные процессы, протекающие в атмосфере и на поверхности. Анализ выполнен в форме числовых оценок, и на их основе построена физическая картина глобальных процессов с участием углерода и углекислого газа. Показано, что парниковый эффект за счет антропогенного углекислого газа, инжектируемого в атмосферу, не является главной причиной наблюдаемого в последние десятилетия роста глобальной температуры, который в свою очередь мал по сравнению с изменениями в прошлом. Рост концентрации углекислого газа в атмосфере при ее современных значениях не влияет на здоровье человека, но увеличивает эффективность фотосинтеза, так что урожаи сельскохозяйственных культур за индустриальный период выросли в полтора раза. Истощение ресурсов ископаемых газа и нефти в этом веке делает актуальной задачу производства синтетического жидкого и газообразного топлива из угля. Показано, что замена метана как топлива водородом наряду с низкой эффективностью ведет также к большим рискам при массовом использовании.



Динамика генерации свободных носителей в кремнии с различным типом легирования при воздействии терагерцевыми импульсами
Abstract
Представлены результаты экспериментальных исследований генерации свободных носителей в легированном кремнии n- и р-типа при воздействии переднего фронта терагерцевого импульса пикосекундной длительности с амплитудой электрического поля до 20 МВ/см. Экспериментально показано, что коэффициент пропускания пробного фемтосекундного лазерного импульса при увеличении напряженности поля от 10 до 20 МВ/см во время действия первого периода терагерцевого импульса одинаков для образцов кремния с различным типом легирования. Проведено численное моделирование динамики заполнения носителями зоны проводимости в кремнии n- и р-типа. Показано, что в ходе увеличения напряженности поля до ~10 МВ/см, когда концентрация электрон-дырочных пар становится соизмеримой с концентрацией примесных электронов (дырок), доминирует электронная (n-тип) или дырочная (р-тип) ударная ионизация, что необходимо учитывать при расчетах.



Short Communications
Подавление волн неустойчивости в стохастически возбужденной турбулентной низкоскоростной струе с использованием упреждающего управления
Abstract
Впервые продемонстрирована возможность упреждающего управления когерентными структурами (волнами неустойчивости) в низкоскоростной стохастически возбужденной турбулентной струе при числах Маха M = 0.11 и Рейнольдса Re = 1.2 × 105. Управляющие возмущения вводились в струю с помощью плазменного актуатора на основе барьерного разряда, размещенного на кромке сопла. Возбуждение струи осуществлялось естественными широкополосными возмущениями. Сигнал для управления был получен с помощью термоанемометра, установленного внутри сопла. Получено снижение пульсаций скорости, ассоциированных с волнами неустойчивости, на 2 дБ по всей длине струи в области чисел Струхаля 0.2 ≤ Sh ≤ 2. Показано, что подавление возмущений приводит к некоторому уменьшению толщины сдвигового слоя на границе струи.



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


