Direct numerical simulation of round turbulent jets in crossflow

2007 ◽  
Vol 574 ◽  
pp. 59-84 ◽  
Author(s):  
SUMAN MUPPIDI ◽  
KRISHNAN MAHESH

Direct numerical simulation is used to study a round turbulent jet in a laminar crossflow. The ratio of bulk jet velocity to free-stream crossflow velocity is 5.7 and the Reynolds number based on the bulk jet velocity and the jet exit diameter is 5000. The mean velocity and turbulent intensities from the simulations are compared to data from the experiments by Su & Mungal (2004) and good agreement is observed. Additional quantities, not available from experiments, are presented. Turbulent kinetic energy budgets are computed for this flow. Examination of the budgets shows that the near field is far from a state of turbulent equilibrium – especially along the jet edges. Also – in the near field – peak kinetic energy production is observed close to the leading edge, while peak dissipation is observed toward the trailing edge of the jet. The results are used to comment upon the difficulty involved in predicting this flow using RANS computations. There exist regions in this flow where the pressure transport term, neglected by some models and poorly modelled by others, is significant. And past the jet exit, the flow is not close to established canonical flows on which most models appear to be based.

1999 ◽  
Vol 387 ◽  
pp. 39-60 ◽  
Author(s):  
H. SHAN ◽  
B. MA ◽  
Z. ZHANG ◽  
F. T. M. NIEUWSTADT

A direct numerical simulation of transitional pipe flow is carried out with the help of a spectral element method and used to investigate the localized regions of ‘turbulent’ flow that are observed in experiments. Two types of such regions can be distinguished: the puff and the slug. The puff, which is generally found at low values of the Reynolds numbers, is simulated for Re = 2200 where the Reynolds number Re is based on the mean velocity UB and pipe diameter D. The slug occurs at a higher Reynolds number and it is simulated for Re = 5000. The computations start with a laminar pipe flow to which is added a prescribed velocity disturbance at a given axial position and for a finite time. The disturbance then evolves further into a puff or slug structure.The simulations confirm the experimentally observed fact that for a puff the velocity near the leading edge changes more gradually than for a slug where an almost discontinuous change is observed. The positions of the leading and trailing edges of the puff and slug are computed from the simulations as a function of time. The propagation velocity of the leading edge is found to be constant and equal to 1.56UB and 1.69UB for the puff and slug, respectively. For the trailing edge the velocity is found to be 0.73UB and 0.52UB, respectively. By rescaling the simulation results obtained at various times to a fixed length, we define an ensemble average. This method is used to compute the average characteristics of the puff and slug such as the spatial distribution of the mean velocity, the turbulent velocity fluctuations and also the wall shear stress. By computing particle trajectories we have investigated the entrainment and detrainment of fluid by a puff and slug. We find that the puff detrains through its trailing edge and entrains through its leading edge. The slug entrains fluid through its leading and through most of its trailing edge. As a consequence the fluid inside the puff is constantly exchanged with fluid outside whereas the fluid inside a slug remains there. These entrainment/detrainment properties which are in agreement with the measurements of Wygnanski & Champagne (1973) imply that the puff has the characteristics of a wave phenomenon while the slug can be characterized more as a material property which travels with the flow.Finally, we have investigated in more detail the velocity field within the puff. In a coordinate system that travels with the mean velocity we find recirculation regions both near the trailing and leading edges which agrees at least qualitatively with experimental data. We also find streamwise vortices, predominantly in the trailing-edge region which have been also observed in experiments and which are believed to play an important role in the dynamics of the transition process.


2011 ◽  
Vol 2011 ◽  
pp. 1-13 ◽  
Author(s):  
Nilanjan Chakraborty ◽  
Andrei N. Lipatnikov

The statistics of mean fluid velocity components conditional in unburned reactants and fully burned products in the context of Reynolds Averaged Navier Stokes (RANS) simulations have been studied using a Direct Numerical Simulation database of statistically planar turbulent premixed flame representing the corrugated flamelets regime combustion. Expressions for conditional mean velocity and conditional velocity correlations which are derived based on a presumed bimodal probability density function of reaction progress variable for unity Lewis number flames are assessed in this study with respect to the corresponding quantities extracted from DNS data. In particular, conditional surface averaged velocities(ui)¯Rsand the velocity correlations(uiu)j¯Rsin the unburned reactants are demonstrated to be effectively modelled by the unconditional velocities(ui)¯Rand velocity correlations(uiuj)¯R, respectively, for the major part of turbulent flame brush with the exception of the leading edge. By contrast, conditional surface averaged velocities(ui)¯Psand the velocity correlations(uiu)j¯Psin fully burned products are shown to be markedly different from the unconditional velocities(ui)¯Pand velocity correlations(uiuj)¯P, respectively.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


2013 ◽  
Vol 444-445 ◽  
pp. 462-467
Author(s):  
Dang Guo Yang ◽  
Yong Hang Wu ◽  
Jin Min Liang ◽  
Jun Liu

A numerical simulation method on noise prediction, which incorporates aerodynamics and sound wave equations based on acoustic analogy, is presented in the paper. Near-field unsteady aerodynamic characteristic can be obtain by large eddy simulation (LES), and far-field propagation of sound waves and spatial sound-field can be obtain by solving the time-domain integral equations of Ffowcs Williams and Hawings (FW-H). Based on the method, a numerical simulation was done on a two-dimension cylinder and a three-dimension flat plate with blunt leading edge. The agreement of numerical results with experiment data validated the Feasibility of the method. The results also indicate that LES can describe vortex generation and shedding in the flow-fields, and FW-H formulation, which has taken time-lag between sound emission and reception times into account, can simulate time-effect of sound propagation toward far-fields.


2019 ◽  
Vol 862 ◽  
pp. 781-815 ◽  
Author(s):  
Y. Kuwata ◽  
Y. Kawaguchi

Lattice Boltzmann direct numerical simulation of turbulent open-channel flows over randomly distributed hemispheres at $Re_{\unicode[STIX]{x1D70F}}=600$ is carried out to reveal the influence of roughness parameters related to a probability density function of rough-surface elevation on turbulence by analysing the spatial and Reynolds- (double-) averaged Navier–Stokes equation. This study specifically concentrates on the influence of the root-mean-square roughness and the skewness, and profiles of turbulence statistics are compared by introducing an effective wall-normal distance defined as a wall-normal integrated plane porosity. The effective distance can completely collapse the total shear stress outside the roughness sublayer, and thus the similarity of the streamwise mean velocity is clearer by introducing the effective distance. In order to examine the influence of the root-mean-square roughness and the skewness on dynamical effects that contribute to an increase in the skin friction coefficient, the triple-integrated double-averaged Navier–Stokes equation is analysed. The main contributors to the skin friction coefficient are found to be turbulence and drag force. The turbulence contribution increases with the root-mean-square roughness and/or the skewness. The drag force contribution, on the other hand, increases in particular with the root-mean-square roughness whereas an increase in the skewness does not increase the drag force contribution because it does not necessarily increase the surface area of the roughness elements. The contribution of the mean velocity dispersion induced by spatial inhomogeneity of the rough surfaces substantially increases with the root-mean-square roughness. A linear correlation is confirmed between the root-mean-square roughness and the equivalent roughness while the equivalent roughness monotonically increases with the skewness. A new correlation function based on the root-mean-square roughness and the skewness is developed with the available experimental and direct numerical simulation data, and it is confirmed that the developed correlation reasonably predicts the equivalent roughness of various types of real rough surfaces.


1997 ◽  
Vol 334 ◽  
pp. 353-379 ◽  
Author(s):  
KRISHNAN MAHESH ◽  
SANJIVA K. LELE ◽  
PARVIZ MOIN

Direct numerical simulation and inviscid linear analysis are used to study the interaction of a normal shock wave with an isotropic turbulent field of vorticity and entropy fluctuations. The role of the upstream entropy fluctuations is emphasized. The upstream correlation between the vorticity and entropy fluctuations is shown to strongly influence the evolution of the turbulence across the shock. Negative upstream correlation between u′ and T′ is seen to enhance the amplification of the turbulence kinetic energy, vorticity and thermodynamic fluctuations across the shock wave. Positive upstream correlation has a suppressing effect. An explanation based on the relative effects of bulk compression and baroclinic torque is proposed, and a scaling law is derived for the evolution of vorticity fluctuations across the shock. The validity of Morkovin's hypothesis across a shock wave is examined. Linear analysis is used to suggest that shock-front oscillation would invalidate the relation between urms and Trms, as expressed by the hypothesis.


2005 ◽  
Vol os-14 (3) ◽  
pp. 1558925005os-14
Author(s):  
Eric M. Moore ◽  
Dimitrios V. Papavassiliou ◽  
Robert L. Shambaugh

An unconventional melt blowing die was analyzed using computational fluid dynamics (CFD). This die has an annular configuration wherein the jet inlet is tapered (the cross-sectional area decreases) as the air approaches the die face. It was found that the flow characteristics of this die are different from conventional slot and annular dies. In particular, for the tapered die the near-field normalized turbulent kinetic energy was found to be lower at shallow die angles. Also, it was found that the peak mean velocity behavior was intermediate between that of conventional annular and slot dies. The centerline turbulence profiles were found to be qualitatively similar to those of annular dies; quantitatively, higher values were present for tapered dies.


Sign in / Sign up

Export Citation Format

Share Document