DNS of viscoelastic turbulent channel flow with rectangular orifice at low Reynolds number

Direct numerical simulations of turbulent viscoelastic-fluid flow in a channel with a rectangular orifice were performed to investigate the influence of viscoelasticity on turbulence statistics and turbulent structures downstream of the orifice. The geometry considered is periodic rectangular orifices with 1:2 expansion. The constitutive equation follows the Giesekus model, valid for polymer (or surfactant) solutions, which are generally capable of reducing the turbulent frictional drag in a smooth channel. The friction Reynolds number and the Weissenberg number were set to 100 and 20-30, respectively. A drag reduction of about 20% was achieved in the viscoelastic flows. The onset Reynolds number for the transition from a symmetric to an asymmetric state was found to be shifted to higher values than that for the Newtonian flow. In the viscoelastic flow, the turbulent kinetic energy was decreased and fewer turbulent eddies were observed, as the Kelvin-Helmholtz vortices were quickly damped. Away from the orifice, quasi-streamwise vortices in the viscoelastic flow were sustained for a longer period, accompanied by energy exchange from elastic energy of the viscoelastic fluid to kinetic energy.     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

Flow structure of microbubble-laden turbulent channel flow measured by PIV combined with the shadow image technique

The turbulence structure of a horizontal channel flow with microbubbles is experimentally investigated using combined particle image velocimetry (PIV) in order to clarify the mechanism of drag reduction caused by microbubbles. A new system which simultaneously measures the liquid phase and the dispersed bubbles is proposed, based on a combination of particle tracking velocimetry (PTV), laser-induced fluorescence (LIF) and the shadow image technique (SIT). To accurately obtain the velocity of the liquid phase, tracer particles which overlap with the bubble shadow images are almost entirely eliminated in the post-processing. Finally, the turbulence characteristics of the flow field are presented, including measurements for both phases, and the bubble effect on the turbulence is quantified.     

Extensional behavior influence on viscoelastic turbulent channel flow

In this work we use in the simulation of a viscoelastic turbulent channel flow a modification of the finitely extensible of non-linear elastic dumbbells with the Peterlin approximation (FENE-P) constitutive model for dilute polymer solutions, applicable to high extensional deformations. The new feature introduced by this modification is that the free energy of the polymer (since it is assumed to be entirely entropically driven) remains always bounded (FENE-PB). The characteristics of the model under steady shear flow, pure elongational flow and transient extensional behavior are presented. It is found that the FENE-PB model is more shear thinning than FENE-P. Most importantly, it also shows a higher extensional viscosity than the FENE-P model. Although the steady-state Trouton ratio asymptotically reaches at high extensional rates the same limit as the FENE-P model, the transition from the Newtonian value is sharper and faster. We use the FENE-PB model in direct numerical simulations (DNS) of viscoelastic turbulent channel flow using spectral approximations. The results for various statistics of the flow and the polymer conformation, when compared against those obtained with the original FENE-P model and the same rheological parameters, show an enhanced polymer-induced drag reduction effect and enhanced deformation of the polymer molecules. This indicates that it is not only the asymptotic but also details from the extensional rheological behavior that matter in quantitatively specifying turbulent viscoelastic flow behavior.     

Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow

The present study reports detailed statistics for velocity and transfer rates of heavy particles dispersed in turbulent boundary layers. Statistics have been extracted from a homogeneous source of data covering a large target parameter space and are used here to analyze the effects of gravity and lift on particle dispersion and deposition in a systematic way. Datasets were obtained performing Direct Numerical Simulation (DNS) of particle-laden turbulent upward/downward flow in a vertical channel. Six values for the particle timescale (the particle Stokes number, St) ranging three orders of magnitude were considered to analyze the deposition process as the controlling mechanism was shifting from turbulent diffusion to inertia-moderated crossing trajectories. For the particle timescales examined, gravity and lift do not influence the qualitative behavior of particles even though velocity profiles and deposition coefficients are modified in a non-monotonic fashion, reaching an optimum for St ? 15. Physical mechanisms for the different behavior are discussed. Raw data and statistics obtained from the present DNS are made available at http://cfd.cineca.it (mirror site: http://158.110.32.35/download/database) and can be used to test simple models and closure equations for multiphase RANS and Large Eddy simulations.     

Intense Reynolds-stress events in turbulent ducts

The aim of the present work is to investigate the role of intense Reynolds shear-stress events in the generation of the secondary flow in turbulent ducts. We consider the connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensional quadrant analysis proposed by Lozano-Durán et al. (J. Fluid Mech., vol. 694, 2012, pp. 100–130). We examine both the geometrical properties of these structures and their contribution to the mean in-plane velocity components, and we perfom a comparison with turbulent channel flow at similar Reynolds number. The contribution to a certain mean quantity is defined as the ensemble average over the detected coherent structures, weighted with their own occupied volume fraction. In the core region of the duct, the contribution of intense events to the wall-normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three-dimensional objects do not differ significantly in both flows. In the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity. However, such contribution is less relevant than that of the complementary portion of the flow not included in such objects. Our results show that strong Reynolds shear-stress events are affected by the presence of a corner but, despite the important role of these structures in the dynamics of wall-bounded turbulent flows, their contribution to the secondary flow is relatively low, both in the core and in the corner.     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

Direct numerical simulation of turbulent wake behind a surface-mounted square cylinder

In this research, direct numerical simulation has been performed to study the turbulent wake behind a wall-mounted square cylinder with aspect ratio 4 and Reynolds number 12 000 (based on the free-stream velocity and obstacle side length) in a developing boundary layer. Owing to the relatively high Reynolds number and high aspect ratio of the cylinder tested, the wake is wide spread behind the cylinder and exhibits complex and energetic vortex motions. The lateral and tip vortex shedding patterns at different frequencies, coherent structures downstream of the obstacle, the production rate and distribution of turbulent kinetic energy, and the instantaneous pressure distribution in the wake region have been thoroughly investigated. In order to validate the numerical results, the first- and second-order flow statistics obtained from the simulations have been carefully compared against available wind-tunnel measurement data.     

Physics and modelling of turbulent particle deposition and entrainment: Review of a systematic study

Deposition and entrainment of particles in turbulent flows are crucial in a number of technological applications and environmental processes. We present a review of recent results from our previous works, which led to physical insights on these phenomena. These results were obtained from a systematic numerical study based on the accurate resolution – Direct Numerical Simulation via a pseudo-spectral approach – of the turbulent flow field, and on Lagrangian tracking of particles under different modelling assumptions. We underline the multiscale aspect of wall turbulence, which has challenged scientists to devise simple theoretical models adequate to fit experimental data, and we show that a sound rendering of wall turbulence mechanisms is required to produce a physical understanding of particle deposition and re-entrainment. This physical understanding can be implemented in more applied simulation techniques, such as Large-Eddy Simulation. Our arguments are based also on the phenomenology of coherent structures and on the examination of flow topology in connection with particle preferential distribution. Starting from these concepts, reasons why theoretical predictions may fail are examined together with the requirements which must be fulfilled by suitable predictive models.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

The low Reynolds number turbulent flow and mixing in a confined impinging jet reactor

Turbulent mixing takes an important role in chemical engineering, especially when the chemical reaction is fast compared to the mixing time. In this context a detailed knowledge of the flow field, the distribution of turbulent kinetic energy (TKE) and its dissipation rate is important, as these quantities are used for many mixing models. For this reason we conduct a direct numerical simulation (DNS) of a confined impinging jet reactor (CIJR) at Re = 500 and Sc = 1. The data is compared with particle image velocimetry (PIV) measurements and the basic flow features match between simulation and experiment. The DNS data is analysed and it is shown that the flow is dominated by a stable vortex in the main mixing duct. High intensities of turbulent kinetic energy and dissipation are found in the impingement zone which decrease rapidly towards the exit of the CIJR. In the whole CIJR the turbulence is not in equilibrium. The strong mixing in the impingement zone leads to a rapid development of a monomodal PDF. Due to the special properties of the flow field, a bimodal PDF is generated in cross-sections downstream the impingement zone, that slowly relaxes under relaminarising conditions. The time required for meso-mixing is dominating the overall mixing performance.     

Direct simulation of turbulent flow and heat transfer in a channel. Part I: Smooth walls

Interest in the use of supercomputers for the direct numerical calculation of turbulence prompts the development of efficient numerical techniques so that calculation at higher Reynolds numbers might be made. This paper presents an efficient pseudo-spectral technique, similar to but different from others that have recently appeared, for the calculation of momentum and heat transfer to a constant-property, turbulent fluid in a two-dimensional channel with walls at different, uniform temperature. The code uses no empiricism, although periodic boundary conditions are used for fluctuating quantities in the streamwise and spanwise directions. Calculations were made for a Prandtl number of 0·72 and Reynolds number based on friction velocity and channel half-height of 180 or 2800 based on channel half-height and average velocity. Calculations of mean velocity profile, turbulence intensities, skewness, flatness, Reynolds stress and eddy diffusivity of heat near a wall compare favourably with experimental results. Representative contour plots of the temperature field near the wall and of the spanwise and streamwise two-point velocity correlations are given. Deficiencies are that the calculation requires many hours on a fast computer with a large high-speed memory and that the grid size in each direction for appropriate resolution is approximately proportional to the square of the Reynolds number and to the Prandtl number raised to some power greater than one.     

Direct numerical simulation of turbulent heat transfer in a square duct with transverse ribs mounted on one wall

Turbulent heat transfer in a ribbed square duct of three different blockage ratios are investigated using direct numerical simulation (DNS). The results of ribbed duct cases are compared with those of a heated smooth duct flow. It is observed that owing to the existence of the ribs and confinement of the duct, organized secondary flows appear as large streamwise-elongated vortices, which intensely interact with the rib elements and four sidewalls and have profound influences on the transport of momentum and thermal energy. This study also shows that the drag and heat transfer coefficients are highly sensitive to the rib height. It is observed that as the rib height increases, the impinging effect of the flow on the windward face of the rib strengthens, leading to enhanced rates of turbulent mixing and heat transfer. The influence of sidewalls and rib height on the turbulence structures associated with temperature fluctuations are analyzed based on multiple tools such as vortex swirling strengths, temporal auto-correlations, spatial two-point cross-correlations, joint probability density functions (JPDF) between the temperature and velocity fluctuations, statistical moments of different orders, and temperature spectra.     

Direct numerical simulation of a turbulent plane Couette-Poiseuille flow with zero-mean shear

Direct numerical simulations of a turbulent Couette-Poiseuille flow with zero-mean-shear at the moving wall (SL-flow) is performed to examine flow features compared to those for a turbulent pure Poiseuille flow (P-flow). Profiles of the streamwise mean velocity, indicator function and ratio of production to dissipation show that the logarithmic region is significantly elongated for the SL-flow compared to that for the P-flow at a similar Reynolds number. In addition, the magnitudes of the Reynolds stresses are found to be larger in both inner and outer layers for the SL-flow than those for the P-flow. The spanwise spectra of the production term in the turbulent kinetic energy equation are examined to provide a structural basis for explaining the statistical behaviors. In addition, because the growth of the energy-containing motions extends to the outer layer further for the SL-flow due to the presence of a positive mean shear throughout the entire wall layer, the self-similar behavior of the energy balance between the production and transport terms with respect to the self-similar wavenumber is found far from the wall. We also find the increase in the number of uniform momentum zones in the SL-flow, revealing the hierarchical distribution of the energy-containing eddies which are composed of multiple uniform momentum zones. These coherent motions lead to the elongation of the logarithmic region for the SL-flow. Finally, investigation of the turbulent energy transfer process in a spectral domain for the SL-flow demonstrates importance of outer layer very-long structures, and these structures attribute to the energy transport process in an entire flow field.     

Influence of curvature on drag reduction by opposition control in turbulent flow along a thin cylinder

Opposition controlled fully developed turbulent flow along a thin cylinder is analyzed by means of direct numerical simulations. The influence of cylinder curvature on the skin-friction drag reduction effect by the classical opposition control (i.e., the radial velocity control) is investigated. The curvature of the cylinder affects the uncontrolled flow statistics; for instance, skin-friction coefficient increases while Reynolds shear stress (RSS) and turbulent intensity decrease. However, the control effect in the case of a small curvature is similar to that in channel flow. When the curvature is large, the maximum drag reduction rate decreased. However, the optimal location of the detection plane is the same as that in a flat plate. Further, the drag reduction effect is achieved even on a high detection plane where the drag increases in the flat plate. Although a difference in the drag reduction effect can be observed with a change in the curvature, its mechanism considered in this analysis based on the transport of the Reynolds stress is similar to that of the flat plate.     

Interfacial instability in turbulent flow over a liquid film in a channel

We revisit the stability of a deformable interface that separates a fully-developed turbulent gas flow from a thin layer of laminar liquid. Although this problem has received considerable attention previously, a model that requires no fitting parameters and that uses a base-state profile that has been validated against experiments is, as yet, unavailable. Furthermore, the significance of wave-induced perturbations in turbulent stresses remains unclear. To address these outstanding issues, we investigate this problem and introduce a turbulent base-state velocity that requires specification of a flow rate or a pressure drop only; no adjustable parameters are necessary. This base state is validated extensively against available experimental data as well as the results of direct numerical simulations. In addition, the effect of perturbations in the turbulent stress distributions is investigated, and demonstrated to be small for cases wherein the liquid layer is thin. The detailed modelling of the liquid layer also elicits two unstable modes, ‘interfacial’ and ‘internal’, with the former being the more dominant of the two. We show that it is possible for interfacial roughness to reduce the growth rate of the interfacial mode in relation to that of the internal one, promoting the latter, to the status of most dangerous mode. Additionally, we introduce an approximate measure to distinguish between ‘slow’ and ‘fast’ waves, the latter being the case for ‘critical-layer’-induced instabilities; we demonstrate that for the parameter ranges studied, the large majority of the waves are ‘slow’. Finally, comparisons of our linear stability predictions are made with experimental data in terms of critical parameters for onset of wave-formation, wave speeds and wavelengths; these yield agreement within the bounds of experimental error.     

DNS of heat transfer in turbulent and transitional channel flow obstructed by rectangular prisms

Direct numerical simulation (DNS) of heat transfer in a channel flow obstructed by rectangular prisms has been performed for Re_τ = 80–20, where Re_τ is based on the friction velocity, the channel half width and the kinematic viscosity. The molecular Prandtl number is set to be 0.71. The flow remains unsteady down to Re_τ = 40 owing to the disturbance induced by the prism. For Re_τ = 30 and 20, the flow results in a steady laminar flow. In the vicinity of the prism, the three-dimensional complex vortices are generated and heat transfer is enhanced. The Reynolds number effect on the time-averaged vortex structure and the local Nusselt number are investigated. The mechanism of the heat transfer enhancement is discussed. In addition, the mean flow parameters such as the friction factor and the Nusselt number are examined in comparison with existing DNS and experimental data.