NUMERICAL RESEARCH ON THE EFFECT OF FIBERS ON THE PROPERTY OF TURBULENT FIBER SUSPENSION IN A CHANNEL

2009 ◽  
Vol 23 (03) ◽  
pp. 509-512 ◽  
Author(s):  
SUHUA SHEN ◽  
JIANZHONG LIN

To explore the rheological property in turbulent channel flow of fiber suspensions, the equation of probability distribution function for mean fiber orientation and the Reynolds averaged Navier-Stokes equation with the term of additional stress resulted from fibers were solved with numerical methods to get the distributions of the mean velocity and turbulent kinetic energy. The simulation results show that the effect of fibers on turbulent channel flow is equivalent to an additional viscosity. The turbulent velocity profiles of fiber suspension become gradually sharper by increasing the fiber concentration and/or decreasing the Reynolds number. The turbulent kinetic energy will increase with increasing Reynolds number and fiber concentration.

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.


Author(s):  
Zacharie M. J. Durand ◽  
Shawn P. Clark ◽  
Mark F. Tachie ◽  
Jarrod Malenchak ◽  
Getnet Muluye

The effect of Reynolds number on three-dimensional offset jets was investigated in this study. An acoustic Doppler velocimeter simultaneously measured all three components of velocity, U, V and W, and turbulence intensity, urms, vrms, and wrms, and all three Reynolds shear stresses, uv, uw, and vw. Turbulent kinetic energy, k, was calculated with all three values of turbulence intensities. Flow measurements were performed at Reynolds numbers of 34,000, 53,000 and 86,000. Results of this experimental study indicate the wall-normal location of maximum mean velocity and jet spread to be independent of Reynolds number. The effects on maximum mean velocity decay are reduced with increasing Reynolds number. Profiles of mean velocities, U, V and W, turbulence intensities, urms, vrms, and wrms, and turbulent kinetic energy, k, show independence of Reynolds number. Reynolds shear stress uv was independent of Reynolds number while the magnitude of uw was reduced at higher Reynolds number.


2009 ◽  
Vol 131 (10) ◽  
Author(s):  
J. D. Schwarzkopf ◽  
C. T. Crowe ◽  
P. Dutta

A closed set of volume averaged equations for modeling turbulence in the carrier phase of particle-laden flows is presented. The equations incorporate a recently developed dissipation transport equation that contains an additional production of dissipation term due to particle surfaces. In the development, it was assumed that each coefficient was the sum of the coefficient for single phase flow and a coefficient quantifying the contribution of the particulate phase. To assess the effects of this additional production term, a numerical model was developed and applied to particles falling in a channel of downward turbulent air flow. Boundary conditions were developed to ensure that the production of turbulent kinetic energy due to mean velocity gradients and particle surfaces balanced with the turbulent dissipation near the wall. The coefficients associated with the production of dissipation due to mean velocity gradients and particle surfaces were varied to assess the effects of the dispersed phase on the carrier phase turbulent kinetic energy across the channel. The results show that the model predicts a decrease in turbulent kinetic energy near the wall with increased particle loading, and that the dissipation coefficients play a critical role in predicting the turbulent kinetic energy in particle-laden turbulent flows.


Volume 1 ◽  
2004 ◽  
Author(s):  
K. T. Christensen ◽  
Y. Wu

Stereo particle-image velocimetry (PIV) has become a widely-used method for studying complex flows because it allows one to acquire instantaneous, three-component velocity data on a planar domain with high spatial resolution. However, the accuracy of such measurements must be carefully evaluated before stereo PIV data can be faithfully used in the development of sophisticated turbulence models, assessment of appropriate computational boundary conditions, and in the validation of advanced computations. To this end, the accuracy of stereo PIV is assessed directly in an actual turbulent environment: two-dimensional turbulent channel flow. This flow is a challenging test of stereo PIV because the turbulent velocity fluctuations are quite small compared to the mean (typically less than ten percent of the mean velocity) and strong velocity gradients exist in the near-wall region. Measurements are made in the streamwise–wall-normal plane along the channel’s spanwise centerline using both stereoscopic and conventional 2-D PIV. A large ensemble of statistically independent velocity realizations are acquired with each method at a friction Reynolds number Reτ = u*h/ν = 934. Single-point statistics are computed from the experimental data and compared to statistics determined from a direct numerical simulation (DNS) of turbulent channel flow at a nearly-identical friction Reynolds number of 940 [5]. Excellent agreement is found in the outer region of the flow (y/h > 0.15, where h is the half-height of the channel). For y/h < 0.15, both the conventional and stereo PIV results differ from the DNS data. These differences are most-likely a manifestation of errors associated with strong velocity gradients and intense turbulent events present in this region of the flow.


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.


2013 ◽  
Vol 8 (1) ◽  
pp. 155892501300800
Author(s):  
Jianzhong Lin ◽  
Suhua Shen ◽  
Xiaoke Ku

The equations of averaged momentum, turbulence kinetic energy, turbulence dissipation rate with the additional term of the fibers, and the equation of probability distribution function for mean fiber orientation are derived and solved numerically for fiber suspension flowing in a turbulent boundary layer. The mathematical model and numerical code are verified by comparing the numerical results with the experimental ones in a turbulent channel flow. The effects of Reynolds number, fiber concentration and fiber aspect-ratio on the mean velocity profile, turbulent kinetic energy, Reynolds stress, turbulent dissipation rate and eddy viscosity coefficient are analyzed. The results show that the velocity profiles become full, and the turbulent kinetic energy, Reynolds stress and eddy viscosity coefficient increase, while turbulent dissipation rate decreases, as the Reynolds number, fiber concentration and fiber aspect-ratio increase. The effect of the fiber aspect-ratio on the turbulent properties is larger than that of the Reynolds number, but smaller than that of the fiber concentration in the range of parameters considered in this paper.


Entropy ◽  
2019 ◽  
Vol 21 (7) ◽  
pp. 675 ◽  
Author(s):  
T.-W. Lee

There are two components in this work that allow for solutions of the turbulent channel flow problem: One is the Galilean-transformed Navier-Stokes equation which gives a theoretical expression for the Reynolds stress (u′v′); and the second the maximum entropy principle which provides the spatial distribution of turbulent kinetic energy. The first concept transforms the momentum balance for a control volume moving at the local mean velocity, breaking the momentum exchange down to its basic components, u′v′, u′2, pressure and viscous forces. The Reynolds stress gradient budget confirms this alternative interpretation of the turbulence momentum balance, as validated with DNS data. The second concept of maximum entropy principle states that turbulent kinetic energy in fully-developed flows will distribute itself until the maximum entropy is attained while conforming to the physical constraints. By equating the maximum entropy state with maximum allowable (viscous) dissipation at a given Reynolds number, along with other constraints, we arrive at function forms (inner and outer) for the turbulent kinetic energy. This allows us to compute the Reynolds stress, then integrate it to obtain the velocity profiles in channel flows. The results agree well with direct numerical simulation (DNS) data at Reτ = 400 and 1000.


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.


Sign in / Sign up

Export Citation Format

Share Document