Direct numerical simulation of turbulence in injection-driven three-dimensional cylindrical flows

2011 ◽  
Vol 670 ◽  
pp. 176-203 ◽  
Author(s):  
JU ZHANG ◽  
THOMAS L. JACKSON

Incompressible turbulent flow in a periodic circular pipe with strong injection is studied as a simplified model for the core flow in a solid-propellant rocket motor and other injection-driven internal flows. The model is based on a multi-scale asymptotic approach. The intended application of the current study is erosive burning of solid propellants. Relevant analysis for easily accessible parameters for this application, such as the magnitudes, main frequencies and wavelengths associated with the near-wall shear, and the assessment of near-wall turbulence viscosity is focused on. It is found that, unlike flows with weak or no injection, the near-wall shear is dominated by the root mean square of the streamwise velocity which is a function of the Reynolds number, while the mean streamwise velocity is only weakly dependent on the Reynolds number. As a result, a new wall-friction velocity $\(u_\tau{\,=\,}\sqrt{\tau_w/\rho}\)$, based on the shear stress derived from the sum of the mean and the root mean square, i.e. $\(\tau_{w,inj} {\,=\,} \mu |{\partial (\bar{u}+u_{rms})}/{\partial r}|_w\)$, is proposed for the scaling of turbulent viscosity for turbulent flows with strong injection. We also show that the mean streamwise velocity profile has an inflection point near the injecting surface.

1995 ◽  
Vol 297 ◽  
pp. 101-122 ◽  
Author(s):  
Masahito Asai ◽  
Michio Nishioka

Subcritical transition in a flat-plate boundary layer is examined experimentally through observing its nonlinear response to energetic hairpin eddies acoustically excited at the leading edge of the boundary-layer plate. When disturbed by the hairpin eddies convecting from the leading edge, the near-wall flow develops local three-dimensional wall shear layers with streamwise vortices. Such local wall shear layers also evolve into hairpin eddies in succession to lead to the subcritical transition beyond the x-Reynolds number Rx = 3.9 × 104, where the momentum thickness Reynolds number Rθ is 127 for laminar Blasius flow without excitation, and is about 150 under the excitation of energetic hairpin eddies. It is found that in terms of u- and v-fluctuations, the intensity of the near-wall activity at this critical station is of almost the same order as or slightly less than that of the developed wall turbulence. The development of wall turbulence structure in this transition is also examined.


2019 ◽  
Vol 864 ◽  
pp. 708-745 ◽  
Author(s):  
Kenzo Sasaki ◽  
Ricardo Vinuesa ◽  
André V. G. Cavalieri ◽  
Philipp Schlatter ◽  
Dan S. Henningson

Three methods are evaluated to estimate the streamwise velocity fluctuations of a zero-pressure-gradient turbulent boundary layer of momentum-thickness-based Reynolds number up to $Re_{\unicode[STIX]{x1D703}}\simeq 8200$, using as input velocity fluctuations at different wall-normal positions. A system identification approach is considered where large-eddy simulation data are used to build single and multiple-input linear and nonlinear transfer functions. Such transfer functions are then treated as convolution kernels and may be used as models for the prediction of the fluctuations. Good agreement between predicted and reference data is observed when the streamwise velocity in the near-wall region is estimated from fluctuations in the outer region. Both the unsteady behaviour of the fluctuations and the spectral content of the data are properly predicted. It is shown that approximately 45 % of the energy in the near-wall peak is linearly correlated with the outer-layer structures, for the reference case $Re_{\unicode[STIX]{x1D703}}=4430$. These identified transfer functions allow insight into the causality between the different wall-normal locations in a turbulent boundary layer along with an estimation of the tilting angle of the large-scale structures. Differences in accuracy of the methods (single- and multiple-input linear and nonlinear) are assessed by evaluating the coherence of the structures between wall-normally separated positions. It is shown that the large-scale fluctuations are coherent between the outer and inner layers, by means of an interactions which strengthens with increasing Reynolds number, whereas the finer-scale fluctuations are only coherent within the near-wall region. This enables the possibility of considering the wall-shear stress as an input measurement, which would more easily allow the implementation of these methods in experimental applications. A parametric study was also performed by evaluating the effect of the Reynolds number, wall-normal positions and input quantities considered in the model. Since the methods vary in terms of their complexity for implementation, computational expense and accuracy, the technique of choice will depend on the application under consideration. We also assessed the possibility of designing and testing the models at different Reynolds numbers, where it is shown that the prediction of the near-wall peak from wall-shear-stress measurements is practically unaffected even for a one order of magnitude change in the corresponding Reynolds number of the design and test, indicating that the interaction between the near-wall peak fluctuations and the wall is approximately Reynolds-number independent. Furthermore, given the performance of such methods in the prediction of flow features in turbulent boundary layers, they have a good potential for implementation in experiments and realistic flow control applications, where the prediction of the near-wall peak led to correlations above 0.80 when wall-shear stress was used in a multiple-input or nonlinear scheme. Errors of the order of 20 % were also observed in the determination of the near-wall spectral peak, depending on the employed method.


2003 ◽  
Vol 125 (6) ◽  
pp. 988-998 ◽  
Author(s):  
Chun-Ho Liu

The turbulence structure and passive scalar (heat) transport in plane Couette flow at Reynolds number equal to 3000 (based on the relative speed and distance between the walls) are studied using direct numerical simulation (DNS). The numerical model is a three-dimensional trilinear Galerkin finite element code. It is found that the structures of the mean velocity and temperature in plane Couette flow are similar to those in forced channel flow, but the empirical coefficients are different. The total (turbulent and viscous) shear stress and total (turbulent and conductive) heat flux are constant throughout the channel. The locations of maximum root-mean-square streamwise velocity and temperature fluctuations are close to the walls, while the location of maximum root-mean-square spanwise and vertical velocity fluctuations are at the channel center. The correlation coefficients between velocities and temperature are fairly constant in the center core of the channel. In particular, the streamwise velocity is highly correlated with temperature (correlation coefficient ≈−0.9). At the channel center, the turbulence production is unable to counterbalance the dissipation, in which the diffusion terms (both turbulent and viscous) bring turbulent kinetic energy from the near-wall regions toward the channel center. The snapshots of the DNS database help explain the nature of the correlation coefficients. The elongated wall streaks for both streamwise velocity and temperature in the viscous sublayer are well simulated. Moreover, the current DNS shows organized large-scale eddies (secondary rotations) perpendicular to the direction of mean flow at the channel center.


2017 ◽  
Vol 816 ◽  
pp. 306-330 ◽  
Author(s):  
C. McKenna ◽  
M. Bross ◽  
D. Rockwell

Impingement of a streamwise-oriented vortex upon a fin, tail, blade or wing represents a fundamental class of flow–structure interaction that extends across a range of applications. It can give rise to unsteady loading known as buffeting and to changes of the lift to drag ratio. These consequences are sensitive to parameters of the incident vortex as well as the location of vortex impingement on the downstream aerodynamic surface, generically designated as a wing. Particle image velocimetry is employed to determine patterns of velocity and vorticity on successive cross-flow planes along the vortex, which lead to volume representations and thereby characterization of the streamwise evolution of the vortex structure as it approaches the downstream wing. This evolution of the incident vortex is affected by the upstream influence of the downstream wing, and is highly dependent on the spanwise location of vortex impingement. Even at spanwise locations of impingement well outboard of the wing tip, a substantial influence on the structure of the incident vortex at locations significantly upstream of the leading edge of the wing was observed. For spanwise locations close to or intersecting the vortex core, the effects of upstream influence of the wing on the vortex are to: decrease the swirl ratio; increase the streamwise velocity deficit; decrease the streamwise vorticity; increase the azimuthal vorticity; increase the upwash; decrease the downwash; and increase the root-mean-square fluctuations of both streamwise velocity and vorticity. The interrelationship between these effects is addressed, including the rapid attenuation of axial vorticity in presence of an enhanced defect of axial velocity in the central region of the vortex. When the incident vortex is aligned with, or inboard of, the tip of the wing, the swirl ratio decreases to values associated with instability of the vortex, thereby giving rise to enhanced values of azimuthal vorticity relative to the streamwise (axial) vorticity, as well as relatively large root-mean-square values of streamwise velocity and vorticity.


2013 ◽  
Vol 718 ◽  
pp. 1-4 ◽  
Author(s):  
B. J. McKeon

AbstractMarusic et al. (J. Fluid Mech., vol. 716, 2013, R3) show the first clear evidence of universal logarithmic scaling emerging naturally (and simultaneously) in the mean velocity and the intensity of the streamwise velocity fluctuations about that mean in canonical turbulent flows near walls. These observations represent a significant advance in understanding of the behaviour of wall turbulence at high Reynolds number, but perhaps the most exciting implication of the experimental results lies in the agreement with the predictions of such scaling from a model introduced by Townsend (J. Fluid Mech., vol. 11, 1961, pp. 97–120), commonly termed the attached eddy hypothesis. The elegantly simple, yet powerful, study by Marusic et al. should spark further investigation of the behaviour of all fluctuating velocity components at high Reynolds numbers and the outstanding predictions of the attached eddy hypothesis.


Author(s):  
Soshi Kawai

This paper addresses the error in large-eddy simulation with wall-modeling (i.e., when the wall shear stress is modeled and the viscous near-wall layer is not resolved): the error in estimating the wall shear stress from a given outer-layer velocity field using auxiliary near-wall RANS equations where convection is not neglected. By considering the behavior of turbulence length scales near a wall, the cause of the errors is diagnosed and solutions that remove the errors are proposed based solidly on physical reasoning. The resulting method is shown to accurately predict equilibrium boundary layers at very high Reynolds number, with both realistic instantaneous fields (without overly elongated unphysical near-wall structures) and accurate statistics (both skin friction and turbulence quantities).


2001 ◽  
Vol 439 ◽  
pp. 217-253 ◽  
Author(s):  
SEDAT F. TARDU

The effect of time-periodical blowing through a spanwise slot on the near-wall turbulence characteristics is investigated. The blowing velocity changes in a cyclic manner from 0 to 5 wall units. The frequency of the oscillations is nearly equal to the median frequency of the near-wall turbulence. The measurements of the wall shear stress and the streamwise velocity are reported and discussed. The flow field near the blowing slot is partly relaminarized during the acceleration phase of the injection velocity which extends 40 wall units downstream. The imposed unsteadiness is confined to the buffer layer, and the time-mean structural parameters under unsteady blowing are found to be close to those of isotropic turbulence in this region. The relaminarized phase is unstable and gives way to a coherent spanwise structure that increases the shear from 80 to 300 wall units downstream of the slot in a predictable way. This phenomenon is strongly imposed-frequency dependent.


2015 ◽  
Vol 783 ◽  
pp. 379-411 ◽  
Author(s):  
I. Marusic ◽  
K. A. Chauhan ◽  
V. Kulandaivelu ◽  
N. Hutchins

In this paper we study the spatial evolution of zero-pressure-gradient (ZPG) turbulent boundary layers from their origin to a canonical high-Reynolds-number state. A prime motivation is to better understand under what conditions reliable scaling behaviour comparisons can be made between different experimental studies at matched local Reynolds numbers. This is achieved here through detailed streamwise velocity measurements using hot wires in the large University of Melbourne wind tunnel. By keeping the unit Reynolds number constant, the flow conditioning, contraction and trip can be considered unaltered for a given boundary layer’s development and hence its evolution can be studied in isolation from the influence of inflow conditions by moving to different streamwise locations. Careful attention was given to the experimental design in order to make comparisons between flows with three different trips while keeping all other parameters nominally constant, including keeping the measurement sensor size nominally fixed in viscous wall units. The three trips consist of a standard trip and two deliberately ‘over-tripped’ cases, where the initial boundary layers are over-stimulated with additional large-scale energy. Comparisons of the mean flow, normal Reynolds stress, spectra and higher-order turbulence statistics reveal that the effects of the trip are seen to be significant, with the remnants of the ‘over-tripped’ conditions persisting at least until streamwise stations corresponding to $Re_{x}=1.7\times 10^{7}$ and $x=O(2000)$ trip heights are reached (which is specific to the trips used here), at which position the non-canonical boundary layers exhibit a weak memory of their initial conditions at the largest scales $O(10{\it\delta})$, where ${\it\delta}$ is the boundary layer thickness. At closer streamwise stations, no one-to-one correspondence is observed between the local Reynolds numbers ($Re_{{\it\tau}}$, $Re_{{\it\theta}}$ or $Re_{x}$ etc.), and these differences are likely to be the cause of disparities between previous studies where a given Reynolds number is matched but without account of the trip conditions and the actual evolution of the boundary layer. In previous literature such variations have commonly been referred to as low-Reynolds-number effects, while here we show that it is more likely that these differences are due to an evolution effect resulting from the initial conditions set up by the trip and/or the initial inflow conditions. Generally, the mean velocity profiles were found to approach a constant wake parameter ${\it\Pi}$ as the three boundary layers developed along the test section, and agreement of the mean flow parameters was found to coincide with the location where other statistics also converged, including higher-order moments up to tenth order. This result therefore implies that it may be sufficient to document the mean flow parameters alone in order to ascertain whether the ZPG flow, as described by the streamwise velocity statistics, has reached a canonical state, and a computational approach is outlined to do this. The computational scheme is shown to agree well with available experimental data.


Sign in / Sign up

Export Citation Format

Share Document