Все выпуски

В настоящее время ведутся активные разработки ядерных реакторов 4-го поколения с жидкометаллическими теплоносителями, в связи с чем актуальными являются расчеты их элементов и узлов с использованием программ трехмерного моделирования. Теплогидравлический анализ реакторных установок с жидкометаллическим теплоносителем признается одним из важнейших направлений комплекса взаимосвязанных задач по обоснованию параметров реакторных установок, включая обоснование безопасности. Сложность получения необходимой информации об условиях эксплуатации реакторного оборудования с жидкометаллическими теплоносителями на основе экспериментальных исследований требует привлечения численного моделирования. В качестве инструмента, описанного в статье исследования, использован отечественный CFD-код FlowVision, который имеет аттестат НТЦ ЯРБ для расчетного обеспечения безопасности ядерных реакторов. Ранее было доказано успешное применение данного расчетного кода для моделирования процессов в ядерных реакторах с натриевым теплоносителем. Поскольку на данный момент в ядерной отрасли в качестве перспективных реакторов рассматриваются установки со свинцово-висмутовым теплоносителем, необходимо обосновать пригодность кода FlowVision также и для моделирования течения такого теплоносителя, что и являлось целью данной работы. В статье приведены результаты численного моделирования потока свинцово-висмутовой эвтектики в пучке теплообменных труб парогенератора АЭС. В рамках CFD-моделирования процессов гидродинамики и теплообмена в пучке теплообменных труб произведены исследования сходимости по сетке, по шагу, выбрана модель турбулентности, определены коэффициенты гидравлического сопротивления решеток и проведено сравнение расчетов с использованием модели $k_\theta^{}$-$e_\theta^{}$ и без нее. По итогам исследования получено, что результаты расчета с использованием $k_\theta^{}$-$e_\theta^{}$-модели турбулентности более точно согласуются с корреляциями. В качестве дополнительной проверки точности результатов выполнена кросс-верификация с ПО STAR-CCM+, полученные результаты лежат в пределах погрешностей использованных для сравнения корреляций.

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.

В данной работе представлена модель пристеночных функций FlowVision, позволяющая моделировать турбулентные течения жидкости и газа около твердых непроницаемых поверхностей на разных сетках. Рассматриваются четыре модели турбулентности: $k-\varepsilon$ FlowVision, $k-\varepsilon$ Стандартная, SST $k-\omega$, SA. Обсуждаются особенности реализации моделей турбулентности в программном комплексе FlowVision. Демонстрируются результаты решения двух тестовых задач.

В настоящей работе приводится модификация опубликованной ранее модели турбулентности KEFV, позволяющая моделировать ламинарно-турбулентный переход в оторвавшемся потоке. Предлагаемая модель содержит настраиваемый механизм, с помощью которого можно задавать генерацию турбулентности в ламинаризованных зонах с большими градиентами скорости. Модель реализована в программном комплексе FlowVision. Демонстрируются результаты расчетов течения воздуха около низкорейнольдсового профиля Eppler-387. Условия течения характеризуются числом Рейнольдса Re = 200 000.

В работе представлены результаты трехмерного численного моделирования теплогидравлических процессов в промежуточном теплообменнике перспективного реактора на быстрых нейтронах с натриевым теплоносителем (БН) с учетом разработанного методического подхода.

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

С целью подтверждения заложенных в проекте теплогидравлических параметров ПТО перспективного реактора БН был разработан методический подход для трехмерного численного моделирования теплообменника, размещенного в корпусе реактора, учитывающий трехмерную картину течения натрия на входе и внутри ПТО, а также обосновывающий рекомендации для упрощения геометрии расчетной модели ПТО. Численное моделирование теплогидравлических процессов в ПТО перспективного реактора БН проводилось с использованием программного комплекса FlowVision со стандартной $k-\varepsilon$-моделью турбулентности и моделью турбулентного теплопереноса LMS. Для повышения представительности численного моделирования трубного пучка ПТО выполнены верификационные расчеты однотрубного и многотрубного теплообменников «натрий – натрий» с соответствующими конструкции ПТО геометрическими характеристиками. Для определения входных граничных условий в модели ПТО выполнен дополнительный трехмерный расчет с учетом неравномерной картины течения в верхней смесительной камере реактора. Расчетная модель ПТО была оптимизирована за счет упрощения дистанционирующих поясов и выбора секторной модели. В результате численного моделирования ПТО получены распределения скорости натрия первого контура, температуры натрия первого и второго контуров. Удовлетворительное согласование результатов расчета с проектными данными по интегральным параметрам подтвердило принятые проектные теплогидравлические характеристики ПТО перспективного реактора БН.

В предлагаемой статье приводятся результаты аналитического и компьютерного исследования хаотической эволюции регулярного поля скорости, возникающего под действием крупномасштабной гармонической вынуждающей силы. Авторами получено аналитическое решение для функции тока течения и ее производных величин (скорости, завихренности, кинетической энергии, энстрофии и палинстрофии). Проведено численное моделирование эволюции течения с помощью пакета программ OpenFOAM (на основе модели несжимаемой среды), а также двух собственных реализаций, использующих приближение слабой сжимаемости (схемы КАБАРЕ и схемы МакКормака). Расчеты проводились на последовательности вложенных сеток с 64², 128², 256², 512², 1024² ячейками для двух характерных (асимптотических) чисел Рейнольдса Rea, характеризующих ламинарную и турбулентную эволюцию течения соответственно. Моделирование показало, что разрушение аналитического решения происходит в обоих случаях. Энергетические характеристики течения обсуждаются на основе кривых энергии, а также скоростей диссипации. Для самой подробной сетки эта величина оказывается на несколько порядков меньше своего гидродинамического (вязкого) аналога. Разрушение регулярной структуры течения наблюдается для любого из численных методов, в том числе на поздних стадиях ламинарной эволюции, когда полученные распределения близки к аналитическим значениям. Можно предположить, что предпосылкой к развитию неустойчивости выступает ошибка, накапливаемая в процессе счета. Эта ошибка приводит к неравномерностям в распределении завихренности и, как следствие, к появлению вихрей различной интенсивности, взаимодействие которых приводит к хаотизации течения. Для исследования процессов производства завихренности мы использовали две интегральные величины, определяемые на ее основе, — интегральные энстрофию ($\zeta$) и палинстрофию $(P)$. Постановка задачи с периодическими граничными условиями позволяет установить простую связь между этими величинами. Кроме того, $\zeta$ может выступать в качестве меры вихреразрешающей способности численного метода, а палинстрофия определяет степень производства мелкомасштабной завихренности.

Данная работа посвящена численному исследованию высокоскоростных слоев смешения сжимаемых потоков. Рассматриваемая задача имеет широкий спектр применения в практических задачах и, несмотря на кажущуюся простоту, является достаточно сложной в плане моделирования, потому что в слое смешения в результате неустойчивости тангенциального разрыва скоростей поток от ламинарного течения переходит к турбулентному режиму. Поэтому полученные численные результаты рассмотренной задачи сильно зависят от адекватности используемых моделей турбулентности. В представленной работе данная задача исследуется на основе двухжидкостного подхода к проблеме турбулентности. Данный подход возник сравнительно недавно и достаточно быстро развивается. Главное преимущество двухжидкостного подхода — в том, что он ведет к замкнутой системе уравнений, тогда как известно, что давний подход Рейнольдса ведет к незамкнутой системе. В работе представлены суть двухжидкостного подхода для моделирования турбулентной сжимаемой среды и методика численной реализации предлагаемой модели. Для получения стационарного решения поставленной задачи применен метод установления и использована теория пограничного слоя Прандтля, которая ведет к упрощенной системе уравнений. В рассматриваемой задаче происходит смешение высокоскоростных потоков. Следовательно, необходимо моделировать также перенос тепла и давление нельзя считать постоянным, как это делается для несжимаемых потоков. При численной реализации конвективные члены в гидродинамических уравнениях аппроксимировались против потока вторым порядка точности в явном виде, а диффузионные члены в правых частях уравнений аппроксимировались центральной разностью в неявном виде. Для реализации полученных уравнений использовался метод прогонки. Для коррекции скорости через давления использован метод SIMPLE. В работе проведено исследование двухжидкостной модели турбулентности при различных начальных возмущениях потока. Полученные численные результаты показали, что хорошее соответствие с известными опытными данными наблюдается при интенсивности турбулентности на входе $0,1 < I < 1 \%$. Для демонстрации эффективности предлагаемой модели турбулентности представлены также данные известных экспериментов, а также результаты моделей $k − kL + J$ и LES. Показано, что двухжидкостная модель по точности не уступает известным современным моделям, а по затрате вычислительных ресурсов является более экономичной.

Преобразование Фурье является основным инструментом для оценки спектральных характеристик турбулентного течения. Применение преобразования Фурье (как правило, дискретного) к первой степени пульсационной продольной или поперечной компоненты вектора скорости позволяет оценить спектральную плотность энергии (ESD) или мощности (PSD). Для оценки ESD и PSD турбулентного сигнала, полученного по результатам численного моделирования обтекания тела, создается массив значений сигнала с дискретизацией во временной или пространственной области. Референтное распределение спектральной характеристики ESD (закон масштабирования) в области волновых чисел инерциального поддиапазона следует из двух гипотез А.Н. Колмогорова и определяется законом $−\frac{5}{3}$. Закон $−\frac{5}{3}$ используется также в большинстве работ для оценки распределения ESD и PSD в частотной области. При этом распределение спектра мощности PSD получается из распределения спектра энергии ESD нормировкой к времени сканирования сигнала. Альтернативная спектральная характеристика энергии (ESS) пульсаций скорости может быть определена преобразованием Фурье для квадрата пульсаций скорости. Размерность ESS в пространственной области совпадает с размерностью закона $−\frac{5}{3}$ А.Н. Колмогорова в области волновых чисел. При дискретизации сигнала во временной области для ESS ранее получен закон масштабирования $−2$ в частотной области. В настоящей работе вводится альтернативная оценка спектра мощности сигнала (PSS), полученная с использованием преобразования Фурье для третьей степени пульсаций скорости. Из гипотез А.Н. Колмогорова следует, что в частотной области закон масштабирования спектра PSS определяется степенью $−\frac{5}{2}$. В качестве приложения рассматривается нестационарное обтекание отсека цилиндрической 3D-поверхности при числе Рейнольдса 3900. Численное моделирование выполнено с использованием пакета ANSYS Fluent на базе решения уравнений Навье – Стокса в несжимаемой постановке. Пространственно-временные характеристики вектора скорости турбулентного потока анализируются с применением правильной ортогональной декомпозиции (POD). Для оценки ESS и PSS используется преобразование Фурье с дискретизацией сигнала по времени.

Рассматриваются вопросы моделирования течения в дальних областях осесимметричного безымпульсного турбулентного следа, возникающего при обтекании тела. Основными особенностями исследуемого течения являются его практически бессдвиговый характер и близость к автомодельному режиму, что подтверждается имеющимися экспериментальными данными. Для описания течения в дальнем следе привлекается полуэмпирическая $e-\varepsilon$-модель турбулентности, основанная на алгебраической модели рейнольдсовых напряжений Роди. С помощью теоретико-группового анализа получена автомодельная редукция привлекаемой модели турбулентности к системе обыкновенных дифференциальных уравнений. Проведено построение приближенного автомодельного решения второго рода соответствующей краевой задачи для редуцированной системы уравнений с помощью метода, основанного на сшивании усечений асимптотических разложений ее решения, найденных в окрестностях граничных точек. Привлекаемый подход также позволил осуществить выделение приближенного значения показателя автомодельности и других неизвестных параметров задачи в процессе выполнения процедуры сшивания асимптотических разложений. Сопоставление сконструированного в виде рядов Пюизё приближенного решения поставленной краевой задачи, включающего 16 аппроксимирующих слагаемых, с численным решением, найденным с помощью метода стрельбы, показало, что максимальная относительная погрешность между решениями не превосходит 2%. Установленное в данной работе приближенное значение показателя автомодельности задачи хорошо согласуется с точным значением, полученным ранее другими авторами в результате аналитического решения соответствующей нелинейной задачи на собственные значения. Проведено сопоставление приближенного решения с доступными экспериментальными данными по осесимметричному безымпульсному турбулентному следу, полученными в ИГиЛ СО РАН, имеет место удовлетворительное соответствие. Построенные в работе приближенные формулы имеют практическое значение и могут быть использованы для быстрой оценки значений гидродинамических величин в дальних областях безымпульсного осесимметричного турбулентного следа. Развиваемый метод построения приближенных решений краевых задач свободной турбулентности может найти применение и в других прикладных областях, которые не связаны с теорией турбулентности.

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

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

В данной работе представлены результаты расчетов аэродинамических характеристик тематической модели летательного аппарата схемы «летающее крыло» в безграничном потоке при разных углах атаки, полученные в рамках первого этапа работы по созданию математической модели рабочей части аэродинамической трубы Т-102 ЦАГИ. Расчеты выполнялись с использованием двухпараметрической k–ε модели турбулентности со специальными пристеночными функциями, приспособленными для расчета отрывных течений. В рамках данной работы исследовались основные продольные аэродинамические характеристики, было выполнено сравнение с результатами экспериментальных исследований в аэродинамической трубе Т-102 ЦАГИ с учетом погрешностей.