skip to main content

Title: Coherence Analysis of Rotating Turbulent Pipe Flow
Rotating and swirling turbulence comprises an important class of flows, not only due to the complex physics that occur, but also due to their relevance to many engineering applications, such as combustion, cyclone separation, mixing, etc. In these types of flows, rotation strongly affects the characteristics and structure of turbulence. However, the underlying turbulent flow phenomena are complex and currently not well understood. The axially rotating pipe is an exemplary prototypical model problem that exhibits these complex turbulent flow physics. By examining the complex interaction of turbulent structures within rotating turbulent pipe flow, insight can be gained into the behavior of rotating flows relevant to engineering applications. Direct numerical simulations are conducted at a bulk Reynolds number up to Re_D = 19,000 with rotation numbers ranging from N = 0 to 3. Coherence analysis, including Proper Orthogonal Decomposition and Dynamic Mode Decomposition, are used to identify the relevant (highest energy) modes of the flow. Studying the influence of these modes on turbulent statistics (i.e. mean statistics, Reynolds stresses, turbulent kinetic energy, and turbulent kinetic energy budgets) allow for a deeper understanding of the effects of coherent turbulent flow structures in rotating flows.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
AIAA Scitech 2020 Forum
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In past experiments, simulations and theoretical analysis, rotation has been shown to dramatically effect the characteristics of turbulent flows, such as causing the mean velocity profile to appear laminar, leading to an overall drag reduction, as well as affecting the Reynolds stress tensor. The axially rotating pipe is an exemplary prototypical model problem that exhibits these complex turbulent flow physics. For this flow, the rotation of the pipe causes a region of turbulence suppression which is particularly sensitive to the rotation rate and Reynolds number. The physical mechanisms causing turbulence suppression are currently not well-understood, and a deeper understanding of these mechanisms is of great value for many practical examples involving swirling or rotating flows, such as swirl generators, wing-tip vortices, axial compressors, hurricanes, etc. In this work, Direct Numerical Simulations (DNS) of rotating turbulent pipe flows are conducted at moderate Reynolds numbers (Re=5300, 11,700, and 19,000) and rotation numbers of N=0 to 3. The main objectives of this work are to firstly quantify turbulence suppression for rotating turbulent pipe flows at different Reynolds numbers as well as study the effects of rotation on turbulence by analyzing the characteristics of the Reynolds stress tensor and the production and dissipation terms of the turbulence budgets. 
    more » « less
  2. New direct numerical simulation data of a fully-developed axially rotating pipe at Re = 5300 and Re = 19, 000 is used to examine the performance of the second-moment closure elliptic blending Reynolds stress model for a range of rotation rates from N=0 to N=3. In agreement with previous studies (using alternative second-moment closure models), the turbulence suppression observed by the DNS is over-predicted. This over-prediction is greatest at Re = 5, 300 and most noticeable in the poor prediction of the ut wt turbulent shear-stress component. At N=3 the flow is completely relaminarized in contrast to the DNS that is only partly relaminarized. The accuracy of the second-moment closure model is superior to the two-equation k − ω SST model which predicts pure solid-body rotation, however, both are equally poor at the highest rotation rates. The accuracy of each model is also assessed for the initial portion of a rotating pipe where in contrast to the fully- developed rotating pipe flow the turbulent suppression is under-predicted compared to the DNS. It is clear that greater work is required to understand the root cause of the poor prediction by these second-moment closure models and further DNS and experimental work is underway to assist this effort. 
    more » « less
  3. Turbulent shear flows are abundant in geophysical and astrophysical systems and in engineering-technology applications. They are often riddled with large-scale secondary flows that drastically modify the characteristics of the primary stream, preventing or enhancing mixing, mass and heat transfer. Using experiments and numerical simulations, we study the possibility of modifying these secondary flows by using superhydrophobic surface treatments that reduce the local shear. We focus on the canonical problem of Taylor–Couette flow, the flow between two coaxial and independently rotating cylinders, which has robust secondary structures called Taylor rolls that persist even at significant levels of turbulence. We generate these structures by rotating only the inner cylinder of the system, and show that an axially spaced superhydrophobic treatment can weaken the rolls through a mismatching surface heterogeneity, as long as the roll size can be fixed. The minimum hydrophobicity of the treatment required for this flow control is rationalized, and its effectiveness beyond the Reynolds numbers studied here is also discussed. 
    more » « less
  4. Highly-resolved, direct numerical simulations of turbulent channel flows with sub- Kolmogorov grid resolution are performed to investigate the characteristics of wall-bounded turbulent flows in the presence of sinusoidal wall waviness. The wall waviness serves as a simplified model to study the effects of well-defined geometric parameters of roughness on the characteristics of wall-bounded turbulent flows. In this study, a two-dimensional wave profile with steepness ranging from 0.06 to 0.25 and wave amplitudes ranging from 9 to 36 wall units were considered. For the smooth and wavy-wall simulations, the Reynolds number based on the friction velocity was kept constant. To study the effects of wave amplitude and wavelength on turbulence, two-dimensional time and spanwise averaged distributions of the mean flow, turbulent kinetic energy, and Reynolds stresses as well as turbulent kinetic energy production and dissipation are examined. Furthermore, in order to provide a more direct comparison with the smooth-wall turbulent channel flow one-dimensional pro- files of these quantities are computed by averaging them over one wavelength of the wave profile. A strong effect of the wall-waviness and, in particular, the wave amplitude and wavelength on the characteristics of the turbulence was obtained. Wall waviness mainly affected the inner flow region while all recorded turbulent statistics collapsed in the outer flow region. Significant reductions in turbulent kinetic energy, production and dissipation were obtained with increasing wave amplitudes when reported in inner scale. While production is lower for all wavy wall cases considered here in comparison to the smooth wall, reducing the wavelength caused an increase in production and a decrease in dissipation. 
    more » « less
  5. Direct numerical simulations are carried out to study the flow structure and transport properties in turbulent Rayleigh–Bénard convection in a vertical cylindrical cell of aspect ratio one with an imposed axial magnetic field. Flows at the Prandtl number $0.025$ and Rayleigh and Hartmann numbers up to $10^{9}$ and $1400$ , respectively, are considered. The results are consistent with those of earlier experimental and numerical data. As anticipated, the heat transfer rate and kinetic energy are suppressed by a strong magnetic field. At the same time, their growth with Rayleigh number is found to be faster in flows at high Hartmann numbers. This behaviour is attributed to the newly discovered flow regime characterized by prominent quasi-two-dimensional structures reminiscent of vortex sheets observed earlier in simulations of magnetohydrodynamic turbulence. Rotating wall modes similar to those in Rayleigh–Bénard convection with rotation are found in flows near the Chandrasekhar linear stability limit. A detailed analysis of the spatial structure of the flows and its effect on global transport properties is reported. 
    more » « less