Self-limiting and regenerative dynamics of perturbation growth on a vortex column

2013 ◽  
Vol 718 ◽  
pp. 39-88 ◽  
Author(s):  
Fazle Hussain ◽  
Eric Stout

AbstractWe study the mechanisms of centrifugal instability and its eventual self-limitation, as well as regenerative instability on a vortex column with a circulation overshoot (potentially unstable) via direct numerical simulations of the incompressible Navier–Stokes equations. The perturbation vorticity (${\boldsymbol{\omega} }^{\prime } $) dynamics are analysed in cylindrical ($r, \theta , z$) coordinates in the computationally accessible vortex Reynolds number, $\mathit{Re}({\equiv }\mathrm{circulation/viscosity} )$, range of 500–12 500, mostly for the axisymmetric mode (azimuthal wavenumber $m= 0$). Mean strain generates azimuthally oriented vorticity filaments (i.e. filaments with azimuthal vorticity, ${ \omega }_{\theta }^{\prime } $), producing positive Reynolds stress necessary for energy growth. This ${ \omega }_{\theta }^{\prime } $ in turn tilts negative mean axial vorticity, $- {\Omega }_{z} $ (associated with the overshoot), to amplify the filament, thus causing instability. (The initial energy growth rate (${\sigma }_{r} $), peak energy (${G}_{\mathit{max}} $) and time of peak energy (${T}_{p} $) are found to vary algebraically with $\mathit{Re}$.) Limitation of vorticity growth, also energy production, occurs as the filament moves the overshoot outward, hence lessening and shifting $\vert {- }{\Omega }_{z} \vert $, while also transporting the core $+ {\Omega }_{z} $, to the location of the filament. We discover that a basic change in overshoot decay behaviour from viscous to inviscid occurs at $Re\sim 5000$. We also find that the overshoot decay time has an asymptotic limit of 45 turnover times with increasing $\mathit{Re}$. After the limitation, the filament generates negative Reynolds stress, concomitant energy decay and hence self-limitation of growth; these inviscid effects are enhanced further by viscosity. In addition, the filament transports angular momentum radially inward, which can produce a new circulation overshoot and renewed instability. Energy decays at the $\mathit{Re}$ studied, but, at higher $\mathit{Re}$, regenerative growth of energy is likely due to the renewed mean shearing. New generation of overshoot and Reynolds stress is examined using a helical ($m= 1$) perturbation. Regenerative energy growth, possibly resulting in even vortex breakup, can be triggered by this new overshoot at practical $\mathit{Re}$ (${\sim }1{0}^{6} $ for trailing vortices), which are currently beyond the computational capability.

2019 ◽  
Vol 869 ◽  
pp. 553-586 ◽  
Author(s):  
Jinlong Wu ◽  
Heng Xiao ◽  
Rui Sun ◽  
Qiqi Wang

Reynolds-averaged Navier–Stokes (RANS) simulations with turbulence closure models continue to play important roles in industrial flow simulations. However, the commonly used linear eddy-viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g. flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness. Recent studies in plane channel flows found that even substituting Reynolds stresses with errors below 0.5 % from direct numerical simulation databases into RANS equations leads to velocities with large errors (up to 35 %). While such an observation may have only marginal relevance to traditional Reynolds stress models, it is disturbing for the recently emerging data-driven models that treat the Reynolds stress as an explicit source term in the RANS equations, as it suggests that the RANS equations with such models can be ill-conditioned. So far, a rigorous analysis of the condition of such models is still lacking. As such, in this work we propose a metric based on local condition number function for a priori evaluation of the conditioning of the RANS equations. We further show that the ill-conditioning cannot be explained by the global matrix condition number of the discretized RANS equations. Comprehensive numerical tests are performed on turbulent channel flows at various Reynolds numbers and additionally on two complex flows, i.e. flow over periodic hills, and flow in a square duct. Results suggest that the proposed metric can adequately explain observations in previous studies, i.e. deteriorated model conditioning with increasing Reynolds number and better conditioning of the implicit treatment of the Reynolds stress compared to the explicit treatment. This metric can play critical roles in the future development of data-driven turbulence models by enforcing the conditioning as a requirement on these models.


Author(s):  
Bohua Sun

This study revisits the Reynolds-averaged Navier--Stokes equations (RANS) and finds that the existing literature is erroneous regarding the primary unknowns and the number of independent unknowns in the RANS. The literature claims that the Reynolds stress tensor has six independent unknowns, but in fact the six unknowns can be reduced to three that are functions of the three velocity fluctuation components, because the Reynolds stress tensor is simply an integration of a second-order dyadic tensor of flow velocity fluctuations rather than a general symmetric tensor. This difficult situation is resolved by returning to the time of Reynolds in 1895 and revisiting Reynolds' averaging formulation of turbulence. The study of turbulence modeling could focus on the velocity fluctuations instead of on the Reynolds stress. An advantage of modeling the velocity fluctuations is, from both physical and experimental perspectives, that the velocity fluctuation components are observable whereas the Reynolds stress tensor is not.


Author(s):  
Bohua Sun

This paper attempts to clarify an long-standing issue about the number of unknowns in the Reynolds-Averaged Navier-Stokes equations (RANS). This study shows that all perspectives regarding the numbers of unknowns in the RANS stem from the misinterpretation of the Reynolds stress tensor. The current literature consider that the Reynolds stress tensor has six unknown components; however, this study shows that the Reynolds stress tensor actually has only three unknown components, namely the three components of fluctuation velocity. This understanding might shed a light to understand the well-known closure problem of turbulence.


1999 ◽  
Vol 122 (1) ◽  
pp. 179-183 ◽  
Author(s):  
Robert E. Spall ◽  
Blake M. Ashby

Solutions to the incompressible Reynolds-averaged Navier–Stokes equations have been obtained for turbulent vortex breakdown within a slightly diverging tube. Inlet boundary conditions were derived from available experimental data for the mean flow and turbulence kinetic energy. The performance of both two-equation and full differential Reynolds stress models was evaluated. Axisymmetric results revealed that the initiation of vortex breakdown was reasonably well predicted by the differential Reynolds stress model. However, the standard K-ε model failed to predict the occurrence of breakdown. The differential Reynolds stress model also predicted satisfactorily the mean azimuthal and axial velocity profiles downstream of the breakdown, whereas results using the K-ε model were unsatisfactory. [S0098-2202(00)01601-1]


Author(s):  
H. Naji ◽  
O. El Yahyaoui ◽  
G. Mompean

The ability of two explicit algebraic Reynolds stress models (EARSMs) to accurately predict the problem of fully turbulent flow in a straight square duct is studied. The first model is devised by Gatski and Rumsey (2001) and the second is the one derived by Wallin and Johansson (2000). These models are studied using a priori procedure based on data resulting from direct numerical simulation (DNS) of the Navier-Stokes equations, which is available for this problem. For this case, we show that the equilibrium assumption for the anisotropy tensor is found to be correct. The analysis leans on the maps of the second and third invariants of the Reynolds stress tensor. In order to handle wall-proximity effects in the near-wall region, damping functions are implemented in the two models. The predictions and DNS obtained for a Reynolds number of 4800 both agree well and show that these models are able to predict such flows.


1989 ◽  
Vol 208 ◽  
pp. 639-656 ◽  
Author(s):  
Fei Li ◽  
Sheila E. Widnall

A model is constructed for perturbations created in the surrounding laminar flow by a turbulent spot in plane Poiseuille flow. The turbulent spot is represented as a distribution of increased Reynolds stress, which travels steadily through the surrounding laminar flow. The Navier-Stokes equations are linearized and are solved by using Fourier transforms in the plane parallel to the channel walls and a finite-difference method in the direction perpendicular to the walls. The travelling Reynolds stress distribution acts as a forcing term in the equations.Numerical results show that a packet of oblique waves are generated around the disturbance when the force is antisymmetric with respect to the channel centreline, whereas no identifiable wave crests are found when the forcing is symmetric. Furthermore, wavelengths of the typical waves composing the packet are insensitive to the size of the region of Reynolds stress. The dependencies of the flow field on Reynolds number and spot speed are investigated. In the case of symmetric forcing, the flow is forced around the disturbance, causing distortions to the basic velocity profiles. These results are in qualitative agreement with experimental observations.


1995 ◽  
Vol 39 (04) ◽  
pp. 263-283 ◽  
Author(s):  
F. Sotiropoulos ◽  
V. C. Patel

ABSTRACT The Reynolds-averaged Navier-Stokes equations are solved to assess the importance of the turbulence model in the prediction of ship stern and wake flows. Solutions are obtained with a two-equation scalar turbulence model and a seven-equation Reynolds-stress tensor model, both of which resolve the flow up to the wall, holding invariant all aspects of the numerical method, including solution domain, initial and boundary conditions, and grid topology and density. Calculations are carried out for two tanker forms used as test cases at recent workshops, and solutions are compared with each other and with experimental data. The comparisons reveal that the Reynolds-stress model accurately predicts most of the experimentally observed flow features in the stern and near-wake regions whereas the two-equation model predicts only the overall qualitative trends. In particular, solutions with the Reynolds-stress model clarify the origin of the stern vortex.


Author(s):  
Charles G. Speziale

The high-Reynolds-number turbulent flows of technological importance contain such a wide range of excited length and time scales that the application of direct or large-eddy simulations is all but impossible for the foreseeable future. Reynolds stress models remain the only viable means for the solution of these complex turbulent flows. It is widely believed that Reynolds stress models are completely ad hoc, having no formal connection with solutions of the full Navier-Stokes equations for turbulent flows. While this belief is largely warranted for the older eddy viscosity models of turbulence, it constitutes a far too pessimistic assessment of the current generation of Reynolds stress closures. It will be shown how secondorder closure models and two-equation models with an anisotropic eddy viscosity can be systematically derived from the Navier-Stokes equations when one overriding assumption is made: the turbulence is locally homogeneous and in equilibrium. A brief review of zero equation models and one equation models based on the Boussinesq eddy viscosity hypothesis will first be provided in order to gain a perspective on the earlier approaches to Reynolds stress modeling. It will, however, be argued that since turbulent flows contain length and time scales that change dramatically from one flow configuration to the next, two-equation models constitute the minimum level of closure that is physically acceptable. Typically, modeled transport equations are solved for the turbulent kinetic energy and dissipation rate from which the turbulent length and time scales are built up; this obviates the need to specify these scales in an ad hoc fashion. While two-equation models represent the minimum acceptable closure, second-order closure models constitute the most complex level of closure that is currently feasible from a computational standpoint. It will be shown how the former models follow from the latter in the equilibrium limit of homogeneous turbulence. However, the two-equation models that are formally consistent with second-order closures have an anisotropic eddy viscosity with strain-dependent coefficients - a feature that most of the commonly used models do not possess.


2015 ◽  
Vol 798 ◽  
pp. 546-550 ◽  
Author(s):  
Asel Beketaeva ◽  
Amr H. Abdalla ◽  
Yekaterina Moisseyeva

The three-dimensional supersonic turbulent flow in presence of symmetric transverse injection of round jet is simulated numerically. The simulation is based on the Favre-averaged Navier-Stokes equations coupled with Wilcox’s turbulence model. The numerical solution is performed using ENO scheme and is validated with the experimental data that include the pressure distribution on the wall in front of the jet in the plane symmetry. The numerical simulation is used to investigate in detail the flow physics for a range of the pressure ratio . The well-known primary shock formations are observed (a barrel shock, a bow shock, and the system of λ-shock waves), and the vortices are identified (horseshoe vortex, an upper vortex, two trailing vortices formed in the separation region and aft of the bow shock wave, two trailing vortices that merge together into one single rotational motion). During the experiment the presence of the new vortices near the wall behind the jet for the pressure ratio is revealed.


Sign in / Sign up

Export Citation Format

Share Document