Simulation of a propelled wake with moderate excess momentum in a stratified fluid

2011 ◽  
Vol 692 ◽  
pp. 28-52 ◽  
Author(s):  
Matthew B. de Stadler ◽  
Sutanu Sarkar

AbstractDirect numerical simulation is used to simulate the turbulent wake behind an accelerating axisymmetric self-propelled body in a stratified fluid. Acceleration is modelled by adding a velocity profile corresponding to net thrust to a self-propelled velocity profile resulting in a wake with excess momentum. The effect of a small to moderate amount of excess momentum on the initially momentumless self-propelled wake is investigated to evaluate if the addition of excess momentum leads to a large qualitative change in wake dynamics. Both the amount and shape of excess momentum are varied. Increasing the amount of excess momentum and/or decreasing the radial extent of excess momentum was found to increase the defect velocity, mean kinetic energy, shear in the velocity gradient and the wake width. The increased shear in the mean profile resulted in increased production of turbulent kinetic energy leading to an increase in turbulent kinetic energy and its dissipation. Slightly larger vorticity structures were observed in the late wake with excess momentum although the differences between vorticity structures in the self-propelled and 40 % excess momentum cases was significantly smaller than suggested by previous experiments. Buoyancy was found to preserve the doubly inflected velocity profile in the vertical direction, and similarity for the mean velocity and turbulent kinetic energy was found to occur in both horizontal and vertical directions. While quantitative differences were observed between cases with and without excess momentum, qualitatively similar evolution was found to occur.

Author(s):  
G. D. MacIsaac ◽  
S. A. Sjolander ◽  
T. J. Praisner

Experimental measurements of the mean and turbulent flow field were preformed downstream of a low-speed linear turbine cascade. The influence of turbulence on the production of secondary losses is examined. Steady pressure measurements were collected using a seven-hole pressure probe and the turbulent flow quantities were measured using a rotatable x-type hotwire probe. Each probe was traversed downstream of the cascade along planes positioned at three axial locations: 100%, 120% and 140% of the axial chord (Cx) downstream of the leading edge. The seven-hole pressure probe was used to determine the local total and static pressure as well as the three mean velocity components. The rotatable x-type hotwire probe, in addition to the mean velocity components, provided the local Reynolds stresses and the turbulent kinetic energy. The axial development of the secondary losses is examined in relation to the rate at which mean kinetic energy is transferred to turbulent kinetic energy. In general, losses are generated as a result of the mean flow dissipating kinetic energy through the action of viscosity. The production of turbulence can be considered a preliminary step in this process. The measured total pressure contours from the three axial locations (1.00, 1.20 and 1.40Cx) demonstrate the development of the secondary losses. The peak loss core in each plane consists mainly of low momentum fluid that originates from the inlet endwall boundary layer. There are, however, additional losses generated as the flow mixes with downstream distance. These losses have been found to relate to the turbulent Reynolds stresses. An examination of the turbulent deformation work term demonstrates a mechanism of loss generation in the secondary flow region. The importance of the Reynolds shear stress to this process is explored in detail.


2018 ◽  
Vol 848 ◽  
pp. 788-820 ◽  
Author(s):  
Shaokai Zheng ◽  
P. J. K. Bruce ◽  
J. M. R. Graham ◽  
J. C. Vassilicos

A group of three multiscale inhomogeneous grids have been tested to generate different types of turbulent shear flows with different mean shear rate and turbulence intensity profiles. Cross hot-wire measurements were taken in a wind tunnel with Reynolds number$Re_{D}$of 6000–20 000, based on the width of the vertical bars of the grid and the incoming flow velocity. The effect of local drag coefficient$C_{D}$on the mean velocity profile is discussed first, and then by modifying the vertical bars to obtain a uniform aspect ratio the mean velocity profile is shown to be predictable using the local blockage ratio profile. It is also shown that, at a streamwise location$x=x_{m}$, the turbulence intensity profile along the vertical direction$u^{\prime }(y)$scales with the wake interaction length$x_{\ast ,n}^{peak}=0.21g_{n}^{2}/(\unicode[STIX]{x1D6FC}C_{D}w_{n})$($\unicode[STIX]{x1D6FC}$is a constant characterizing the incoming flow condition, and$g_{n}$,$w_{n}$are the gap and width of the vertical bars, respectively, at layer$n$) such that$(u^{\prime }/U_{n})^{2}\unicode[STIX]{x1D6FD}^{2}(C_{D}w_{n}/x_{\ast ,n}^{peak})^{-1}\sim (x_{m}/x_{\ast ,n}^{peak})^{b}$, where$\unicode[STIX]{x1D6FD}$is a constant determined by the free-stream turbulence level,$U_{n}$is the local mean velocity and$b$is a dimensionless power law constant. A general framework of grid design method based on these scalings is proposed and discussed. From the evolution of the shear stress coefficient$\unicode[STIX]{x1D70C}(x)$, integral length scale$L(x)$and the dissipation coefficient$C_{\unicode[STIX]{x1D716}}(x)$, a simple turbulent kinetic energy model is proposed that describes the evolution of our grid generated turbulence field using one centreline measurement and one vertical profile of$u^{\prime }(y)$at the beginning of the evolution. The results calculated from our model agree well with our measurements in the streamwise extent up to$x/H\approx 2.5$, where$H$is the height of the grid, suggesting that it might be possible to design some shear flows with desired mean velocity and turbulence intensity profiles by designing the geometry of a passive grid.


2010 ◽  
Vol 652 ◽  
pp. 373-404 ◽  
Author(s):  
KYLE A. BRUCKER ◽  
SUTANU SARKAR

Direct numerical simulations (DNS) of axisymmetric wakes with canonical towed and self-propelled velocity profiles are performed atRe= 50 000 on a grid with approximately 2 billion grid points. The present study focuses on a comparison between towed and self-propelled wakes and on the elucidation of buoyancy effects. The development of the wake is characterized by the evolution of maxima, area integrals and spatial distributions of mean and turbulence statistics. Transport equations for mean and turbulent energies are utilized to help understand the observations. The mean velocity in the self-propelled wake decays more rapidly than the towed case due to higher shear and consequently a faster rate of energy transfer to turbulence. Buoyancy allows a wake to survive longer in a stratified fluid by reducing the 〈u1′u3′〉 correlation responsible for the mean-to-turbulence energy transfer in the vertical direction. This buoyancy effect is especially important in the self-propelled case because it allows regions of positive and negative momentum to become decoupled in the vertical direction and decay with different rates. The vertical wake thickness is found to be larger in self-propelled wakes. The role of internal waves in the energetics is determined and it is found that, later in the evolution, they can become a dominant term in the balance of turbulent kinetic energy. The non-equilibrium stage, known to exist for towed wakes, is also shown to exist for self-propelled wakes. Both the towed and self-propelled wakes, atRe= 50000, are found to exhibit a time span when, although the turbulence is strongly stratified as indicated by small Froude number, the turbulent dissipation rate decays according to inertial scaling.


1984 ◽  
Vol 8 (3) ◽  
pp. 165-170
Author(s):  
L.P. Hackman ◽  
A.B. Strong ◽  
G.D. Raithby

This paper reports predictions of the mean velocity, the turbulent kinetic energy and the pressure and skin friction coefficients for turbulent flow over a backward facing step based on the standard k – ε closure for the turbulence shear stresses. In previous publications, errors due to the numerical algorithm as distinct from the turbulence model have been carefully assessed using different numerical schemes and finite volume geometries and it is argued that the current results are numerically accurate. Thus one can now assess the accuracy of the k – ε model of turbulence independently of numerical error. The results predicted herein were found to be in reasonable agreement with relevant experimental data.


2012 ◽  
Vol 134 (6) ◽  
Author(s):  
G. D. MacIsaac ◽  
S. A. Sjolander ◽  
T. J. Praisner

Experimental measurements of the mean and turbulent flow field were preformed downstream of a low-speed linear turbine cascade. The influence of turbulence on the production of secondary losses is examined. Steady pressure measurements were collected using a seven-hole pressure probe and the turbulent flow quantities were measured using a rotatable x-type hotwire probe. Each probe was traversed downstream of the cascade along planes positioned at three axial locations: 100%, 120%, and 140% of the axial chord (Cx) downstream of the leading edge. The seven-hole pressure probe was used to determine the local total and static pressure as well as the three mean velocity components. The rotatable x-type hotwire probe, in addition to the mean velocity components, provided the local Reynolds stresses and the turbulent kinetic energy. The axial development of the secondary losses is examined in relation to the rate at which mean kinetic energy is transferred to turbulent kinetic energy. In general, losses are generated as a result of the mean flow dissipating kinetic energy through the action of viscosity. The production of turbulence can be considered a preliminary step in this process. The measured total pressure contours from the three axial locations (1.00, 1.20, and 1.40Cx) demonstrate the development of the secondary losses. The peak loss core in each plane consists mainly of low momentum fluid that originates from the inlet endwall boundary layer. There are, however, additional losses generated as the flow mixes with downstream distance. These losses have been found to relate to the turbulent Reynolds stresses. An examination of the turbulent deformation work term demonstrates a mechanism of loss generation in the secondary flow region. The importance of the Reynolds shear stresses to this process is explored in detail.


2017 ◽  
Vol 822 ◽  
pp. 848-879 ◽  
Author(s):  
J. Singh ◽  
M. Rudman ◽  
H. M. Blackburn

Direct numerical simulations of turbulent pipe flow of power-law fluids at $Re_{\unicode[STIX]{x1D70F}}=323$ are analysed in order to understand the way in which shear thinning or thickening affects first- and second-order flow statistics including turbulent kinetic energy production, transport and dissipation in such flows. The results show that with shear thinning, near-wall streaks become weaker and the axial and azimuthal correlation lengths of axial velocity fluctuations increase. Viscosity fluctuations give rise to an additional shear stress term in the mean momentum equation which is negative for shear-thinning fluids and which increases in magnitude as the fluid becomes more shear thinning: for an equal mean wall shear stress, this term increases the mean velocity gradient in shear-thinning fluids when compared to a Newtonian fluid. Consequently, the mean velocity profile in power-law fluids deviates from the law of the wall $U_{z}^{+}=y^{+}$ in the viscous sublayer when traditional near-wall scaling is used. Consideration is briefly given to an alternative scaling that allows the law of wall to be recovered but which results in loss of a common mean stress profile. With shear thinning, the mean viscosity increases slightly at the wall and its profile appears to be approximately logarithmic in the velocity log layer. Through analysis of the turbulent kinetic energy budget, undertaken here for the first time for generalised Newtonian fluids, it is shown that shear thinning decreases the overall turbulent kinetic energy production but widens the wall-normal region where it is generated. Additional dissipation terms in the mean flow and turbulent kinetic energy budget equations arise from viscosity fluctuations; with shear thinning, these result in a net decrease in the total viscous dissipation. The overall effect of shear thinning on the turbulent kinetic energy budget is found to be largely confined to the inner layers, $y^{+}\lesssim 60$.


2016 ◽  
Vol 790 ◽  
pp. 57-70 ◽  
Author(s):  
L. Djenidi ◽  
R. A. Antonia ◽  
N. Lefeuvre ◽  
J. Lemay

Self-preservation (SP) solutions on the axis of a turbulent round jet are derived for the transport equation of the second-order structure function of the turbulent kinetic energy ($k$), which may be interpreted as a scale-by-scale (s.b.s.) energy budget. The analysis shows that the mean turbulent energy dissipation rate, $\overline{{\it\epsilon}}$, evolves like $x^{-4}$ ($x$ is the streamwise direction). It is important to stress that this derivation does not use the constancy of the non-dimensional dissipation rate parameter $C_{{\it\epsilon}}=\overline{{\it\epsilon}}u^{\prime 3}/L_{u}$ ($L_{u}$ and $u^{\prime }$ are the integral length scale and root mean square of the longitudinal velocity fluctuation respectively). We show, in fact, that the constancy of $C_{{\it\epsilon}}$ is simply a consequence of complete SP (i.e. SP at all scales of motion). The significance of the analysis relates to the fact that the SP requirements for the mean velocity and mean turbulent kinetic energy (i.e. $U\sim x^{-1}$ and $k\sim x^{-2}$ respectively) are derived without invoking the transport equations for $U$ and $k$. Experimental hot-wire data along the axis of a turbulent round jet show that, after a transient downstream distance which increases with Reynolds number, the turbulence statistics comply with complete SP. For example, the measured $\overline{{\it\epsilon}}$ agrees well with the SP prediction, i.e. $\overline{{\it\epsilon}}\sim x^{-4}$, while the Taylor microscale Reynolds number $Re_{{\it\lambda}}$ remains constant. The analytical expression for the prefactor $A_{{\it\epsilon}}$ for $\overline{{\it\epsilon}}\sim (x-x_{o})^{-4}$ (where $x_{o}$ is a virtual origin), first developed by Thiesset et al. (J. Fluid Mech., vol. 748, 2014, R2) and rederived here solely from the SP analysis of the s.b.s. energy budget, is validated and provides a relatively simple and accurate method for estimating $\overline{{\it\epsilon}}$ along the axis of a turbulent round jet.


2014 ◽  
Vol 742 ◽  
pp. 171-191 ◽  
Author(s):  
Matteo Bernardini ◽  
Sergio Pirozzoli ◽  
Paolo Orlandi

AbstractThe high-Reynolds-number behaviour of the canonical incompressible turbulent channel flow is investigated through large-scale direct numerical simulation (DNS). A Reynolds number is achieved ($Re_{\tau } = h/\delta _v \approx 4000$, where $h$ is the channel half-height, and $\delta _v$ is the viscous length scale) at which theory predicts the onset of phenomena typical of the asymptotic Reynolds number regime, namely a sensible layer with logarithmic variation of the mean velocity profile, and Kolmogorov scaling of the velocity spectra. Although higher Reynolds numbers can be achieved in experiments, the main advantage of the present DNS study is access to the full three-dimensional flow field. Consistent with refined overlap arguments (Afzal & Yajnik, J. Fluid Mech. vol. 61, 1973, pp. 23–31; Jiménez & Moser, Phil. Trans. R. Soc. Lond. A, vol. 365, 2007, pp. 715–732), our results suggest that the mean velocity profile never achieves a truly logarithmic profile, and the logarithmic diagnostic function instead exhibits a linear variation in the outer layer whose slope decreases with the Reynolds number. The extrapolated value of the von Kármán constant is $k \approx 0.41$. A near logarithmic layer is observed in the spanwise velocity variance, as predicted by Townsend’s attached eddy hypothesis, whereas the streamwise variance seems to exhibit a shoulder, perhaps being still affected by low-Reynolds-number effects. Comparison with previous DNS data at lower Reynolds number suggests enhancement of the imprinting effect of outer-layer eddies onto the near-wall region. This mechanisms is associated with excess turbulence kinetic energy production in the outer layer, and it reflects in flow visualizations and in the streamwise velocity spectra, which exhibit sharp peaks in the outer layer. Associated with the outer energy production site, we find evidence of a Kolmogorov-like inertial range, limited to the spanwise spectral density of $u$, whereas power laws with different exponents are found for the other spectra. Finally, arguments are given to explain the ‘odd’ scaling of the streamwise velocity variances, based on the analysis of the kinetic energy production term.


Author(s):  
Chung-Chu Chen ◽  
Tong-Miin Liou

Laser-Doppler velocimetry (LDV) measurements are presented of turbulent flow in a two-pass square-sectioned duct simulating the coolant passages employed in gas turbine blades under rotating and non-rotating conditions. For all cases studied, the Reynolds number characterized by duct hydraulic diameter (Dh) and bulk mean velocity (Ub) was fixed at 1 × 104. The rotating case had a range of rotation number (Ro = ΩDh/Ub) from 0 to 0.2. It is found that both the skewness of streamwise mean velocity and magnitude of secondary-flow velocity increase linearly, and the magnitude of turbulence intensity level increases non-linearly with increasing Ro. As Ro is increased, the curvature induced symmetric Dean vortices in the turn for Ro = 0 is gradually dominated by a single vortex most of which impinges directly on the outer part of leading wall. The high turbulent kinetic energy is closely related to the dominant vortex prevailing inside the 180-deg sharp turn. For the first time, the measured flow characteristics account for the reported spanwise heat transfer distributions in the rotating channels, especially the high heat transfer enhancement on the leading wall in the turn. For both rotating and non-rotating cases, the direction and strength of the secondary flow with respect to the wall are the most important fluid dynamic factors affecting local heat transfer distributions inside a 180-deg sharp turn. The role of the turbulent kinetic energy in affecting the overall enhancement of heat transfer is well addressed.


Author(s):  
Shinji Honami ◽  
Wataru Tsuboi ◽  
Takaaki Shizawa

This paper presents the effect of flame dome depth on the total pressure performance and flow behavior in a sudden expansion region of the combustor diffuser without flow entering the dome head. The mean velocity and turbulent Reynolds stress profiles in the sudden expansion region were measured by a Laser Doppler Velocitmetry (LDV) system. The experiments show that total pressure loss is increased, when flame dome depth is increased. Installation of an inclined combuster wall in the sudden expansion region is suggested from the viewpoint of a control of the reattaching flow. The inclined combustor wall is found to be effective in improvement of the diffuser performance. Better characteristics of the flow rate distribution into the branched channels are obtained in the inclined wall configuration, even if the distorted velocity profile is provided at the diffuser inlet.


Sign in / Sign up

Export Citation Format

Share Document