Prediction of mono-disperse gas–liquid turbulent flow in a vertical pipe

A two-fluid model in the Eulerian–Eulerian framework has been implemented for the prediction of gas volume fraction, mean phasic velocities, and the liquid phase turbulence properties for gas–liquid upward flow in a vertical pipe. The governing two-fluid transport equations are discretized using the finite volume method and a low Reynolds number $k-\varepsilon$ model is used to predict the turbulence field for the continuous liquid phase. In the present analysis, a fully developed one-dimensional flow is considered where the gas volume fraction profile is predicted using the radial force balance for the bubble phase. The current study investigates: (1) the turbulence modulation terms which represent the effect of bubbles on the liquid phase turbulence in the k–ε transport equations; (2) the role of the bubble induced turbulent viscosity compared to turbulence generated by shear; and (3) the effect of bubble size on the radial forces which results in either a center-peak or a wall-peak in the gas volume fraction profiles. The results obtained from the current simulation are generally in good agreement with the experimental data, and somewhat improved over the predictions of some previous numerical studies.     

Scaling of the near-wall layer beneath turbulent separated flow

This paper is a continuation of an earlier paper [P.E. Hancock, Velocity scales in the near-wall layer beneath reattaching turbulent separated and boundary layer flows, Eur. J. Mech. B Fluids 24 (2005) 425–438] in which it is proposed that each Reynolds stress has its own velocity scale. Two of these, $u_{\tau }^{\prime }$ and $w_{\tau }^{\prime }$, are directly related by definition to the r.m.s. of the wall-shear-stress fluctuations (${\tau }_x$ and ${\tau }_z$) in the streamwise and transverse directions. They are also velocity scales for the true dissipation of the turbulent kinetic energy and the Kolmogorov velocity and length scales at the surface. From asymptotic considerations it is shown that the other two scales are related to averages involving instantaneous gradients of wall-shear-stress fluctuations. The measurements, made using pulsed-wire anemometry into the viscous sublayer, show that $u_{\tau }^{\prime }$ and $w_{\tau }^{\prime }$ are also the velocity scales for the respective streamwise and transverse fourth-order velocity moments, together with the viscous velocity scale ($\nu /y$). Normalised, the fourth-order moments show an inner-layer-like behaviour independent of both position and direction, like that seen in the second-order moments [P.E. Hancock, Velocity scales in the near-wall layer beneath reattaching turbulent separated and boundary layer flows, Eur. J. Mech. B Fluids 24 (2005) 425–438]. However, not surprisingly, the third order moments exhibit an effect of mean shear, seen in the skewing of the probability distributions. Though not measured directly, the measurements imply the behaviour of the averaged products of fluctuations in wall-shear-stress and wall-pressure-gradient ($\overline{{\tau }_x\partial p/\partial x}$ and $\overline{{\tau }_z\partial p/\partial z}$). Normalised, they also are independent of position and direction. Some of the results presented apply more generally to the near-wall region beneath turbulent flow.     

Progress in the extension of a second-moment closure for turbulent environmental flows

An advanced second-moment closure for the double-averaged turbulence equations of porous medium and vegetation flows is proposed. It treats three kinds of second moments which appear in the double-averaged momentum equation. They are the dispersive covariance, the volume averaged (total) Reynolds stress and the micro-scale Reynolds stress. The two-component-limit pressure–strain correlation model is applied to model the total Reynolds stress equation whilst a novel scale-similarity non-linear $k''–''\epsilon$ two-equation eddy viscosity model is employed for the micro-scale turbulence. For the dispersive covariance, an algebraic relation is applied. Model validation in several fully developed homogeneous porous medium flows, porous channel flows and aquatic vegetation canopy flows is performed with satisfactory agreement with the data.     

Large eddy simulation of serration effects on an ultra-high-bypass-ratio engine exhaust jet

Serrated jet nozzles are considered to be an efficient and practical passive control approach for jet noise. However, some fundamental mechanisms of serration effects on jet noise are not fully understood, especially in terms of the sound source. In this paper, a high-fidelity simulation framework using large-eddy simulation (LES) is demonstrated to predict near-field turbulence and far-field acoustics from an ultra-high-bypass-ratio engine with round and serrated nozzles. Far-field sound is predicted using Ffowcs Willams–Hawkings (FWH) integration. The results show that the serrated nozzle increases mixing near the nozzle and hence the turbulence decay rate, reducing the turbulence level downstream. The serrations shift the energy from the low frequencies to the high frequencies and decrease overall sound pressure levels by about 3 dB over the low-frequency range. Sound sources are analysed based on fourth-order space–time correlations. There are six major source components ($R_{1111}$, $R_{2222}$, $R_{3333}$, $R_{1313}$, $R_{1212}$, and $R_{2323}$) inside the jet shear layers. The serrations are able to reduce the amplitude of these source terms, causing them to decay rapidly to a level below the round nozzle jet within 2D downstream of the nozzle.     

Response and wake patterns of two side-by-side elastically supported circular cylinders in uniform laminar cross-flow

Vortex-induced vibrations (VIV) of two side-by-side elastically supported circular cylinders in a uniform flow with the Reynolds number of 100 are numerically investigated by using the immersed boundary method. The cylinders are constrained to oscillate in the cross-flow direction with a center-to-center spacing ratio $T/D$ ranging from 2 to 5. The structural damping is set to zero to enable large vibration amplitudes in the range of reduced velocity $U_r=3-10$. It is found that the proximity of the cylinders does not have a significant impact to the lock-in region and cylinder responses, except at a small spacing ratio of $T/D=2$. The critical spacing ratio is determined as $T/D=4$ and beyond that the interaction between the cylinders is negligible. The following six near-wake patterns are observed; the irregular pattern, in-phase flip-flopping pattern, out-of-phase flip-flopping pattern, in-phase-synchronized pattern, antiphase-synchronized pattern and the biased antiphase-synchronized pattern. These patterns are plotted in a plane of $U_r$ and $T/D$, together with approximate borderlines to distinguish one region from the others. The time histories, spectral features and wavelet transform contours of drag and lift forces are presented to elucidate the mechanisms of the in-phase and out-of-phase flip-flopping phenomena. It is established that the in-phase flip-flopping stems from the long-short near-wake pattern and its low-frequency flip-over, whereas the out-of-phase pattern originates from the large vortex shedding from the fictitious bluff-body with an augmented characteristic length.     

Numerical verification of the weak turbulent model for swell evolution

The purpose of this article is to numerically verify the theory of weak turbulence. We have performed numerical simulations of an ensemble of nonlinearly interacting free gravity waves (a swell) by two different methods: by solving the primordial dynamical equations describing the potential flow of an ideal fluid with a free surface, and by solving the kinetic Hasselmann equation, describing the wave ensemble in the framework of the theory of weak turbulence. In both cases we have observed effects predicted by this theory: frequency downshift, angular spreading and formation of a Zakharov–Filonenko spectrum $I_{\omega }\sim {\omega }^{\mathrm{-4}}$. To achieve quantitative coincidence of the results obtained by different methods, we have to augment the Hasselmann kinetic equation by an empirical dissipation term $S_{\mathrm{diss}}$ modeling the coherent effects of white-capping. Using the standard dissipation terms from the operational wave predicting model (WAM) leads to a significant improvement on short times, but does not resolve the discrepancy completely, leaving the question about the optimal choice of $S_{\mathrm{diss}}$ open. In the long run, WAM dissipative terms essentially overestimate dissipation.     

A priori study for the modelling of velocity–interface correlations in the stratified air–water flows

This work concerns the modelling of stratified two-phase turbulent flows with interfaces. We consider an equation for an intermittency function $\alpha (\mathbf{x}\text{,}t)$ which denotes the probability of finding an interface at a given time t and a given point $\mathbf{x}$. In Wacławczyk and Oberlack (2011) a model for the unclosed terms in this equation was proposed. Here, we investigate the performance of this model by a priori tests, and finally, based on the a priori data discuss its possible modification and improvements.     

Experimental investigation on wake characteristics behind a yawed square cylinder

The wake vortical structures of a square cylinder at different yaw angles to the incoming flow ($\alpha =$0°, 15°, 30° and 45°) are studied using a one-dimensional (1D) hot-wire vorticity probe at a Reynolds number (Re) of about 3600. The results are compared with those obtained in a yawed circular cylinder wake. The Strouhal number (St_N) as well as the mean drag coefficient ($C_{D_N}$), normalized by the velocity component normal to the cylinder axis, follow the independent principle (IP) satisfactorily up to $\alpha =$40°. Using the phase-averaging analysis, both the coherent and the remaining contributions of velocity and vorticity are quantified. The flow patterns of the coherent spanwise vorticity (${\omega }_z$) display obvious Kármán vortex streets and their maximum concentrations decrease as $\alpha$ increases. Similar phenomena are also shown in the coherent contours of the streamwise ($u$) and transverse ($v$) velocities as well as the Reynolds shear stress ($uv$). The contours of the spanwise velocity ($w$) and Reynolds shear stress ($uw$), however, experience an increasing trend for the maximum concentrations with increasing yaw angle. These results indicate an enhancement of the three-dimensionality of the wake and the reduction of vortex shedding strength as $\alpha$ increases. While general similarities to the wake behind a yawed circular cylinder are found in terms of flow features, some differences between the two wakes at different yaw angles are highlighted.     

Instability of Navier slip flow of liquidsInstabilité des écoulements glissant de liquides.

We investigate the stability problem related to the basic slip flows of liquids in plane microchannels by using the Navier slip concept. We found that if the Navier slip parameter ($N_s$) equals 0.06, the critical Reynolds number (${\mathit{Re}}_{\mathrm{cr}}$) becomes 213.6. There are short-wave instabilities, however, when we further increase $N_s$ to 0.07 or 0.08. ${\mathit{Re}}_{\mathrm{cr}}$ becomes 132.9 for $N_s=0.08$ if we neglect the short-wave instability. To cite this article: A.K.-H. Chu, C. R. Mecanique 332 (2004).