Численное моделирование турбулентного течения и теплопереноса в канале методом спектральных элементов: влияние пространственного разрешения

Abstract

Актуальность: турбулентные течения жидкости и газа реализуются во многих инженерных приложениях. В частности, транспортировка нефти и газа по трубопроводам, которая достигает 12 % в структуре энергопотребления нефтедобычи, осуществляется при высоких числах Рейнольдса, т. е. в турбулентном режиме. Задача контролировать и управлять характеристиками турбулентных потоков является крайне актуальной как с научной, так и с экономической точки зрения. Развивающиеся компьютерные технологии делают методы численного моделирования перспективными для решения подобных оптимизационных задач. Цель работы. Для аккуратного численного расчета уравнений Навье-Стокса особенно при высоких числах Рейнольдса необходимо использовать разностные схемы с низкой диссипацией, т. е. высокого порядка точности аппроксимации. Хорошим методом для пространственной дискретизации является метод спектральных элементов. В работе вышеобозначенным методом рассчитывается каноническое турбулентное течение в прямоугольном канале. Цель работы состоит в исследовании точности метода спектральных элементов при расчете отфильтрованных уравнений Навье-Стокса с замыканием подсеточных напряжений при помощи распространенной динамической модели Смагоринского. Методы. Для численного исследования используется открытый пакет программ Nek5000 для расчета дифференциальных уравнений в частных производных при помощи метода спектральных элементов. Этот код использует неструктурированные гексагедральные расчетные сетки и протокол MPI для ускорения вычислений в мультипроцессорном режиме. Результаты. Проведено исследование турбулентного течения в канале при числах Рейнольдса 2800 и 6800. Результаты прямого численного моделирования (DNS) уравнений Навье-Стокса и моделирования методом крупных вихрей (LES) очень хорошо согласуются с профилями средней скорости, температуры и их пульсаций из литературных данных. Даже для самых грубых расчетных сеток с количеством узлов в 25-30 раз меньшим, чем у DNS, удалось очень точно воспроизвести пик пульсаций скорости около стенки, соответствующий полосчатым когерентным вихревым структурам. Как следствие, скорость трения предсказывается с погрешностью <3,0 %, что указывает на перспективность использования метода спектральных элементов.

Relevance. Turbulent flows are readily met in a variety of engineering applications. In particular, oil and gas transport through pipelines reaching 12 % of the oil-related power consumption is performed at very high Reynolds numbers, i. e. in the turbulent regime. Control of turbulent flow characteristics is the problem of high importance both from scientific and economical point of view. Rapidly developing computing capabilities make the numerical modeling a promising tool to solve various engineering optimization tasks. The aim of the research. For accurate numerical solution of the Navier-Stokes equations especially at high Reynolds numbers it is necessary to use numerical schemes with low dissipation, i.e. with a high order of approximation. A good method for spatial discretization is the spectral element method (SEM). In the work a canonical turbulent flow in a plane channel is simulated with SEM. The main aim of the research is to investigate the accuracy of SEM when applied to the filtered Navier-Stokes equations closed with the use of the dynamic Smagorinsky model. Methods. For numerical simulations of partial differential equations with SEM the authors have used the open-source code Nek5000 which employs unstructured hexahedral meshes and MPI for parallel computing. Results. The authors carried out the numerical computations of the turbulent channel flow at the Reynolds numbers 2800 and 6800. The results of the direct numerical simulation (DNS) of the Navier-Stokes equations and Large eddy simulations (LES) are in good agreement with the data from the literature, i. e. the profiles of the time-averaged velocity, temperature and their fluctuations. Even for very coarse LES mesh with 25-30 times less grid points than for DNS the accuracy of the fluctuations profile near the wall representing streaky structures is very high. As a result the friction velocity is predicted within the <3,0 % error indicating promising prospects for SEM.