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.     

Wake response of a stationary finite-sized particle in a turbulent channel flow

Here we consider the effect of a finite-sized stationary particle in a channel flow of modest turbulence at _Reτ=178.12${\mathit{Re}}_{\tau }=178.12$. The size of particle is varied such that the particle Reynolds number ranges from about 40 to 450. The location of the particle is chosen to be either in the buffer layer (_yp+=17.81)$(y_{p+}=17.81)$ or at the channel center. Fully resolved direct numerical simulations of the turbulent channel flow around the particles is performed. Here the ambient turbulence intensity relative to the mean velocity seen by the particle is large (I=23.16%)$(I=23.16\%)$ in the buffer region, while it is substantially lower (I=4.09%)$(I=4.09\%)$ at the channel center. We present results on turbulence modulation due to the particle in terms of wake dynamics and vortex shedding.     

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.     

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.     

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.     

Numerical study of turbulent plane Couette flow at low Reynolds number

A direct numerical simulation of fully developed turbulent plane Couette flow has been performed. Unsteady large-scale structures, which contributed to the instantaneous energy level, were observed. These evolving and drifting vortical structures vanished after time-averaging, and the resulting mean velocity and streamwise turbulence intensity compared favourably with recent laboratory data.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

A Reynolds stress model for turbulent flow of homogeneous polymer solutions

Using a priori analyses of direct numerical simulation (DNS) data, a Reynolds stress model (RSM) is developed to account for the influence of polymer additives on turbulent flow over a wide range of flow conditions. The Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model is utilized to evaluate the polymer contribution to the stress tensor. Thirteen DNS data sets are used to analyze the budgets of elastic stress–velocity gradient correlations as well as Reynolds stress and dissipation transport. Closures are developed in the framework of the RSM model for all the required unknown and non-linear terms. The polymer stresses, velocity profiles, turbulent flow statistics and the percentage of friction drag reduction predicted by the RSM model are in good agreement with present and those obtained from independent DNS data over a wide range of rheological and flow parameters.     

Large eddy simulation of a developing turbulent boundary layer at a low Reynolds number

A spectral code has been used to simulate a developing turbulent boundary layer at low Reynolds number Re_θ (based on free stream velocity and momentum thickness) between 353 and 576. The starting field was generated by allowing a step change of temperature to diffuse outwards from one wall in a fully developed channel flow. The thermal boundary layer so created was conditionally sampled to convert it into a momentum boundary layer with an irrotational free stream region, a process which is justified by appeal to experiments. This initial field was allowed to develop until the momentum boundary layer thickness δ₉₉₅ had grown to about 1·5 times its original thickness. The results of the simulation have been compared with a wide range of experimental data. The outcome of this comparison is generally very satisfactory; the main trends of the experiments are well reproduced and our simulation supplements and extends the existing sets of experimental data. The simulation also gives pressure statistics which cannot be obtained experimentally. In particular, it gives the contribution of pressure diffusion to the balance equations for the Reynolds stress and indicates the error produced by omitting this term.     

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.     

Turbulent transport dissimilarity with modulated turbulence structure in channel flow of viscoelastic fluid

This experimental study investigated the turbulent transport dissimilarity with a modulated turbulence structure in a channel flow of a viscoelastic fluid using simultaneous particle image velocimetry and planar laser-induced fluorescence measurements. An instantaneous dye concentration field with fluctuating velocity vectors showed that mass was transferred by hierarchically large-scale wavy motions with inclination. A co-spectral analysis showed that the spatial phase modulation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation corresponded to the relaxation time. The occurrence of intense dye concentration fluctuation and small streamwise velocity fluctuation in a thin boundary layer caused dissimilar turbulent transport because of the non-zero negative correlation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation only on large scales. This explains the turbulent transport dissimilarity which leads to the zero averaged Reynolds shear stress and non-zero wall-normal turbulent mass flux.     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Measurements of turbulent flow in a channel at low Reynolds numbers

Normal and streamwise components of the velocity fields of turbulent flow in a channel at low Reynolds numbers have been measured with laser-Doppler techniques. The experiments duplicate the conditions used in current direct numerical simulations of channel flow, and good, but not exact, agreement is found for single-point moments through fourth order. In order to eliminate LDV velocity bias and to measure velocity spectra, the mean time interval between LDV signals was adjusted to be much smaller than the smallest turbulence time scale. Spectra of the streamwise and normal components of velocity at locations spanning the channel are presented.     

Effects of dimensional wall temperature on Reynolds stress budgets in a supersonic turbulent channel flow with thermally perfect gas

ABSTRACT

Direct numerical simulations of temporally evolving supersonic turbulent channel flow of thermally perfect gas are conducted at Mach number 3.0 and Reynolds number 4800, combined with constant dimensional wall temperatures from 149.075 to 1788.90?K to study the influence of dimensional wall temperature on the characteristics of Reynolds stress budgets. It is found that, as the dimensional wall temperature increases, the production, diffusion, pressure–velocity gradient correlation and dissipation terms increase, whereas the compressibility-related term decreases. This is mainly due to variations in mean flow properties. The mechanism of inter-component transfer (ICT) is insensitive to the dimensional wall temperature. The ICT relating to the pressure–velocity gradient correlation term can be divided into inner and outer regions, and the critical position separating these regions is at the semi-local scaling of approximately 16 irrespective of the different dimensional wall temperature.     

Sensitivity analysis of low Reynolds number channel flow using the finite volume method

An unsteady finite volume‐based fractional step algorithm solved on a staggered grid has been developed for computing design sensitivity parameters in two‐dimensional flows. Verification of the numerical code is performed for the case of low Reynolds number, pressure‐driven flow through a straight channel, which has an exact steady‐state solution to the Navier–Stokes equations. Sensitivity of the flow to the channel height, fluid viscosity, and imposed pressure gradient is considered. Three different numerical techniques for computing the design sensitivity parameters: finite difference, complex‐step differentiation, and sensitivity equation method (SEM), are compared in terms of numerical error (relative to the exact solution), computational expense, and ease of implementation. Results indicate that, of all the three methods, complex step is the most accurate and requires the least computational time. In addition, treatment of the boundary conditions in SEM is addressed, within the framework of the present finite volume approach, with special attention given to parameter dependence in the boundary conditions. Error estimation based on the Grid Convergence Index provides a good indication of the exact error in the SEM solutions. An example of application of the use of sensitivity parameters to estimate the propagation of input uncertainty through the numerical simulation is also provided. Copyright © 2007 John Wiley & Sons, Ltd.     

Modeling of drag reduction in turbulent channel flow with hydrophobic walls by FVM method and weakly-compressible flow equations

In this paper the effects of hydrophobic wall on skin-friction drag in the channel flow are investigated through large eddy simulation on the basis of weaklycompressible flow equations with the MacCormack's scheme on collocated mesh in the FVM framework. The slip length model is adopted to describe the behavior of the slip velocities in the streamwise and spanwise directions at the interface between the hydrophobic wall and turbulent channel flow. Simulation results are presented by analyzing flow behaviors over hydrophobic wall with the Smagorinky subgrid-scale model and a dynamic model on computational meshes of different resolutions. Comparison and analysis are made on the distributions of timeaveraged velocity, velocity fluctuations, Reynolds stress as well as the skin-friction drag. Excellent agreement between the present study and previous results demonstrates the accuracy of the simple classical second-order scheme in representing turbulent vertox near hydrophobic wall. In addition, the relation of drag reduction efficiency versus time-averaged slip velocity is established. It is also foundthat the decrease of velocity gradient in the close wall region is responsible for the drag reduction. Considering its advantages of high calculation precision and efficiency, the present method has good prospect in its application to practical projects.     

Acoustic upwinding for sub‐ and super‐sonic turbulent channel flow at low Reynolds number

A recently developed asymmetric implicit fifth‐order scheme with acoustic upwinding for the spatial discretization for the characteristic waves is applied to the fully compressible, viscous and non‐stationary Navier–Stokes equations for sub‐ and super‐sonic, mildly turbulent, channel flow (Re_τ=360). For a Mach number of 0.1, results are presented for uniform (32³, 64³ and 128³) and non‐uniform (expanding wall‐normal, 32³ and 64³) grids and compared to the (incompressible) reference solution found in (J. Fluid. Mech. 1987; 177 :133–166). The results for uniform grids on 128³ and 64³ nodes show high resemblance with the reference solution. Expanding grids are applied on 64³‐ and 32³‐node grids. The capability of the proposed technique to solve compressible flow is first demonstrated by increasing the Mach number to 0.3, 0.6 and 0.9 for isentropic flow on the uniform 64³‐grid. Next, the flow speed is increased to Ma=2. The results for the isothermal‐wall supersonic flows give very good agreement with known literature results. The velocity field, the temperature and their fluctuations are well resolved. This means that in all presented (sub‐ and super‐sonic) cases, the combination of acoustic upwinding and the asymmetric high‐order scheme provides sufficient high wave‐number damping and low wave‐number accuracy to give numerically stable and accurate results. Copyright © 2007 John Wiley & Sons, Ltd.