The correspondence between drag enhancement and vortical structures in turbulent Taylor–Couette flows with polymer additives: a study of curvature dependence

2019 ◽  
Vol 881 ◽  
pp. 602-616 ◽  
Author(s):  
Jiaxing Song ◽  
Hao Teng ◽  
Nansheng Liu ◽  
Hang Ding ◽  
Xi-Yun Lu ◽  
...  

We report direct numerical simulation results that clearly elucidate the mechanism that leads to curvature dependence of drag enhancement (DE) in viscoelastic turbulent Taylor–Couette flow. Change in the angular momentum transport and its inherent link to transitions in vortical flow structures have been explored to depict the influence of the curvature of the flow geometry on DE. Specifically, it has been demonstrated that a transition in vortical structures with increasing radius ratio leads to weakening and elimination of the small-scale Görtler vortices and development and better organization (occupying the entire gap) of large-scale Taylor vortices as also evinced by the patterns of angular momentum current. The commensurate change in DE and its underlying mechanism are examined by contributions of convective flux and polymeric stress to the angular momentum current. The present finding paves the way for capturing highly localized elastic turbulence structures in direct numerical simulation by increasing geometry curvature in traditional turbulent curvilinear flows.

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Xingtuan Yang ◽  
Nan Gui ◽  
Gongnan Xie ◽  
Jie Yan ◽  
Jiyuan Tu ◽  
...  

This study investigates the anisotropic characteristics of turbulent energy dissipation rate in a rotating jet flow via direct numerical simulation. The turbulent energy dissipation tensor, including its eigenvalues in the swirling flows with different rotating velocities, is analyzed to investigate the anisotropic characteristics of turbulence and dissipation. In addition, the probability density function of the eigenvalues of turbulence dissipation tensor is presented. The isotropic subrange of PDF always exists in swirling flows relevant to small-scale vortex structure. Thus, with remarkable large-scale vortex breakdown, the isotropic subrange of PDF is reduced in strongly swirling flows, and anisotropic energy dissipation is proven to exist in the core region of the vortex breakdown. More specifically, strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region, whereas nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows.


2002 ◽  
Vol 450 ◽  
pp. 377-407 ◽  
Author(s):  
S. A. STANLEY ◽  
S. SARKAR ◽  
J. P. MELLADO

Turbulent plane jets are prototypical free shear flows of practical interest in propulsion, combustion and environmental flows. While considerable experimental research has been performed on planar jets, very few computational studies exist. To the authors' knowledge, this is the first computational study of spatially evolving three-dimensional planar turbulent jets utilizing direct numerical simulation. Jet growth rates as well as the mean velocity, mean scalar and Reynolds stress profiles compare well with experimental data. Coherency spectra, vorticity visualization and autospectra are obtained to identify inferred structures. The development of the initial shear layer instability, as well as the evolution into the jet column mode downstream is captured well.The large- and small-scale anisotropies in the jet are discussed in detail. It is shown that, while the large scales in the flow field adjust slowly to variations in the local mean velocity gradients, the small scales adjust rapidly. Near the centreline of the jet, the small scales of turbulence are more isotropic. The mixing process is studied through analysis of the probability density functions of a passive scalar. Immediately after the rollup of vortical structures in the shear layers, the mixing process is dominated by large-scale engulfing of fluid. However, small-scale mixing dominates further downstream in the turbulent core of the self-similar region of the jet and a change from non-marching to marching PDFs is observed. Near the jet edges, the effects of large-scale engulfing of coflow fluid continue to influence the PDFs and non-marching type behaviour is observed.


2021 ◽  
Vol 926 ◽  
Author(s):  
Jiaxing Song ◽  
Fenghui Lin ◽  
Nansheng Liu ◽  
Xi-Yun Lu ◽  
Bamin Khomami

The flow physics of inertio-elastic turbulent Taylor–Couette flow for a radius ratio of $0.5$ in the Reynolds number ( $Re$ ) range of $500$ to $8000$ is investigated via direct numerical simulation. It is shown that as $Re$ is increased the turbulence dynamics can be subdivided into two distinct regimes: (i) a low $Re \leqslant 1000$ regime where the flow physics is essentially dominated by nonlinear elastic forces and the main contribution to transport and mixing of momentum, stress and energy comes from large-scale flow structures in the bulk region and (ii) a high $Re \geqslant 5000$ regime where inertial forces govern the flow physics and the flow dynamics is mainly governed by small-scale flow structures in the near-wall region. Flow–microstructure coupling analysis reveals that the elastic Görtler instability in the near-wall region is triggered via significant polymer extension and commensurately high hoop stresses. This instability gives rise to small-scale elastic vortical structures identified as elastic Görtler vortices which are present at all $Re$ considered. In fact, these vortices develop herringbone streaks near the inner wall that have a longer average life span than their Newtonian counterparts due to their elastic origin. Examination of the budgets of mean streamwise enstrophy, mean kinetic energy, turbulent kinetic energy and Reynolds shear stress demonstrates that increasing fluid inertia hinders the generation of elastic stresses, leading to a monotonic reduction of the elastic-related effects on the flow physics.


Author(s):  
Nilanjan Chakraborty ◽  
Alexander Herbert ◽  
Umair Ahmed ◽  
Hong G. Im ◽  
Markus Klein

AbstractA three-dimensional Direct Numerical Simulation (DNS) database of statistically planar $$H_{2} -$$ H 2 - air turbulent premixed flames with an equivalence ratio of 0.7 spanning a large range of Karlovitz number has been utilised to assess the performances of the extrapolation relations, which approximate the stretch rate and curvature dependences of density-weighted displacement speed $$S_{d}^{*}$$ S d ∗ . It has been found that the correlation between $$S_{d}^{*}$$ S d ∗ and curvature remains negative and a significantly non-linear interrelation between $$S_{d}^{*}$$ S d ∗ and stretch rate has been observed for all cases considered here. Thus, an extrapolation relation, which assumes a linear stretch rate dependence of density-weighted displacement speed has been found to be inadequate. However, an alternative extrapolation relation, which assumes a linear curvature dependence of $$S_{d}^{*}$$ S d ∗ but allows for a non-linear stretch rate dependence of $$S_{d}^{*}$$ S d ∗ , has been found to be more successful in capturing local behaviour of the density-weighted displacement speed. The extrapolation relations, which express $$S_{d}^{*}$$ S d ∗ as non-linear functions of either curvature or stretch rate, have been found to capture qualitatively the non-linear curvature and stretch rate dependences of $$S_{d}^{*}$$ S d ∗ more satisfactorily than the linear extrapolation relations. However, the improvement comes at the cost of additional tuning parameter. The Markstein lengths LM for all the extrapolation relations show dependence on the choice of reaction progress variable definition and for some extrapolation relations LM also varies with the value of reaction progress variable. The predictions of an extrapolation relation which involve solving a non-linear equation in terms of stretch rate have been found to be sensitive to the initial guess value, whereas a high order polynomial-based extrapolation relation may lead to overshoots and undershoots. Thus, a recently proposed extrapolation relation based on the analysis of simple chemistry DNS data, which explicitly accounts for the non-linear curvature dependence of the combined reaction and normal diffusion components of $$S_{d}^{*}$$ S d ∗ , has been shown to exhibit promising predictions of $$S_{d}^{*}$$ S d ∗ for all cases considered here.


2019 ◽  
Vol 15 (S359) ◽  
pp. 312-317
Author(s):  
Francoise Combes

AbstractGas fueling AGN (Active Galaxy Nuclei) is now traceable at high-resolution with ALMA (Atacama Large Millimeter Array) and NOEMA (NOrthern Extended Millimeter Array). Dynamical mechanisms are essential to exchange angular momentum and drive the gas to the super-massive black hole. While at 100pc scale, the gas is sometimes stalled in nuclear rings, recent observations reaching 10pc scale (50mas), may bring smoking gun evidence of fueling, within a randomly oriented nuclear gas disk. AGN feedback is also observed, in the form of narrow and collimated molecular outflows, which point towards the radio mode, or entrainment by a radio jet. Precession has been observed in a molecular outflow, indicating the precession of the radio jet. One of the best candidates for precession is the Bardeen-Petterson effect at small scale, which exerts a torque on the accreting material, and produces an extended disk warp. The misalignment between the inner and large-scale disk, enhances the coupling of the AGN feedback, since the jet sweeps a large part of the molecular disk.


2015 ◽  
Vol 764 ◽  
pp. 362-394 ◽  
Author(s):  
T. Dairay ◽  
V. Fortuné ◽  
E. Lamballais ◽  
L.-E. Brizzi

AbstractDirect numerical simulation (DNS) of an impinging jet flow with a nozzle-to-plate distance of two jet diameters and a Reynolds number of 10 000 is carried out at high spatial resolution using high-order numerical methods. The flow configuration is designed to enable the development of a fully turbulent regime with the appearance of a well-marked secondary maximum in the radial distribution of the mean heat transfer. The velocity and temperature statistics are validated with documented experiments. The DNS database is then analysed focusing on the role of unsteady processes to explain the spatial distribution of the heat transfer coefficient at the wall. A phenomenological scenario is proposed on the basis of instantaneous flow visualisations in order to explain the non-monotonic radial evolution of the Nusselt number in the stagnation region. This scenario is then assessed by analysing the wall temperature and the wall shear stress distributions and also through the use of conditional averaging of velocity and temperature fields. On one hand, the heat transfer is primarily driven by the large-scale toroidal primary and secondary vortices emitted periodically. On the other hand, these vortices are subjected to azimuthal distortions associated with the production of radially elongated structures at small scale. These distortions are responsible for the appearance of very high heat transfer zones organised as cold fluid spots on the heated wall. These cold spots are shaped by the radial structures through a filament propagation of the heat transfer. The analysis of probability density functions shows that these strong events are highly intermittent in time and space while contributing essentially to the secondary peak observed in the radial evolution of the Nusselt number.


2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.


2018 ◽  
Vol 859 ◽  
pp. 819-838 ◽  
Author(s):  
Josef Hasslberger ◽  
Sebastian Ketterl ◽  
Markus Klein ◽  
Nilanjan Chakraborty

The local flow topology analysis of the primary atomization of liquid jets has been conducted using the invariants of the velocity-gradient tensor. All possible small-scale flow structures are categorized into two focal and two nodal topologies for incompressible flows in both liquid and gaseous phases. The underlying direct numerical simulation database was generated by the one-fluid formulation of the two-phase flow governing equations including a high-fidelity volume-of-fluid method for accurate interface propagation. The ratio of liquid-to-gas fluid properties corresponds to a diesel jet exhausting into air. Variation of the inflow-based Reynolds number as well as Weber number showed that both these non-dimensional numbers play a pivotal role in determining the nature of the jet break-up, but the flow topology behaviour appears to be dominated by the Reynolds number. Furthermore, the flow dynamics in the gaseous phase is generally less homogeneous than in the liquid phase because some flow regions resemble a laminar-to-turbulent transition state rather than fully developed turbulence. Two theoretical models are proposed to estimate the topology volume fractions and to describe the size distribution of the flow structures, respectively. In the latter case, a simple power law seems to be a reasonable approximation of the measured topology spectrum. According to that observation, only the integral turbulent length scale would be required as an input for the a priori prediction of the topology size spectrum.


2007 ◽  
Vol 594 ◽  
pp. 59-69 ◽  
Author(s):  
MATTHEW J. RINGUETTE ◽  
MINWEI WU ◽  
M. PINO MARTÍN

We demonstrate that data from direct numerical simulation of turbulent boundary layers at Mach 3 exhibit the same large-scale coherent structures that are found in supersonic and subsonic experiments, namely elongated, low-speed features in the logarithmic region and hairpin vortex packets. Contour plots of the streamwise mass flux show very long low-momentum structures in the logarithmic layer. These low-momentum features carry about one-third of the turbulent kinetic energy. Using Taylor's hypothesis, we find that these structures prevail and meander for very long streamwise distances. Structure lengths on the order of 100 boundary layer thicknesses are observed. Length scales obtained from correlations of the streamwise mass flux severely underpredict the extent of these structures, most likely because of their significant meandering in the spanwise direction. A hairpin-packet-finding algorithm is employed to determine the average packet properties, and we find that the Mach 3 packets are similar to those observed at subsonic conditions. A connection between the wall shear stress and hairpin packets is observed. Visualization of the instantaneous turbulence structure shows that groups of hairpin packets are frequently located above the long low-momentum structures. This finding is consistent with the very large-scale motion model of Kim & Adrian (1999).


Sign in / Sign up

Export Citation Format

Share Document