Response of the Sphere Wake to Freestream Fluctuations

Direct numerical simulations have been used to investigate the response of the wake of a sphere to freestream fluctuations. This study has been motivated by the need to understand particle-induced turbulence enhancement in particulate flows. A sequence of simulations of flow past a sphere have been carried out where the frequency and amplitude of the freestream fluctuations and the flow Reynolds number has been varied systematically. It has been suggested that turbulence enhancement is primarily caused by vortex shedding from particles (Gore and Crowe, 1989; Hetsroni, 1989). Our simulations of the forced wake indicate that turbulence enhancement may be attributed to natural vortex shedding only when the freestream fluctuation level is low and the Reynolds number is greater than about 300. In addition to natural vortex shedding, the current simulations also suggest another mechanism for turbulence enhancement. It is found that in the presence of freestream fluctuations, the wake behaves like an oscillator and returns large amounts of kinetic energy to the surrounding fluid at resonance. This mechanism is not associated with natural vortex shedding and is therefore capable of enhancing freestream turbulence even at Reynolds numbers less than 300. Simulations also indicate that when the turbulence intensity of the carrier fluid is high, this resonance mechanism might be solely responsible for turbulence enhancement. Finally, our simulations also suggest a possible explanation for the correlation between turbulence enhancement and the ratio of the particle size to the size of energy containing eddies of turbulence found by Gore and Crowe (1989). Received 5 October 1999 and accepted 14 October 1999     

Simulation of the Effect of the Freestream Turbulence Parameters on Heat Transfer in an Unsteady Boundary Layer

A variant of the two-parameter turbulence model which makes it possible continuously to calculate a flow region with laminar, transition and turbulent regimes is proposed for investigating the flow under conditions of high freestream turbulence intensity. It is shown that the properties of the thermal transition can be theoretically described using the quasi-steady turbulence model in the case of periodic freestream velocity distribution. The numerical results are compared with theoretical and experimental data. The approach proposed is developed for determining the combined effect of the parameters of harmonic fluctuations of the external velocity and freestream turbulence on the heat transfer characteristics on a flat plate with different boundary conditions for the enthalpy.     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

亚临界雷诺数下圆柱绕流的大涡模拟

本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟．求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算．计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的     

Simulations for gas flows in microgeometries using the direct simulation Monte Carlo method

Micro gas flows are often encountered in MEMS devices and classical CFD could not accurately predict the flow and thermal behavior due to the high Knudsen number. Therefore, the gas flow in microgeometries was investigated using the direct simulation Monte Carlo (DSMC) method. New treatments for boundary conditions are verified by simulations of micro-Poiseuille flow, compared with the previous boundary treatments, and slip analytical solutions of the continuum theory. The orifice flow and the corner flow in microchannels are simulated using the modified DSMC codes. The predictions were compared with existing experimental phenomena as well as predictions using continuum theory. The results showed that the channel geometry significantly affects the microgas flow. In the orifice flow, the flow separation occurred at very small Reynolds numbers. In the corner flow, no flow separation occurred even with a high driving pressure. The DSMC results agreed well with existing experimental information.     

水中开孔腔流激振荡控制实验研究

水中开孔腔流激振荡是水下航行器的一类突出噪声源. 为探究有效抑制水中开孔腔流激振荡的控制方法和作用特性, 首先以水下航行器的表面开孔结构为对象设计了开孔腔模型, 并提出一种基于来流边界层分流原理的流动控制装置——前缘分流体, 借助循环水洞装置对水中开孔腔的流激振荡特性及其控制进行了实验研究. 通过沿流向和展向安装于腔底的动态压力传感器测量腔内脉动压力, 分别从腔内脉动压力的频谱特性和空间分布特性两方面, 探讨了水中开孔腔在不同流速下的流激振荡特性和前缘分流体对水中开孔腔流激振荡的控制效果, 并对前缘分流体的主要作用机理进行了分析. 研究结果表明: 水中开孔腔流激振荡形式以剪切层自持振荡为主, 在流速较低时, 如2.4 m/s, 就会产生稳定的自持振荡, 且具有随流速升高而急剧增大的趋势; 前缘分流体对水中开孔腔绕流自持振荡具有良好的抑制效果, 且抑制效果随流速增加而显著提升, 对腔内脉动压力频谱峰值和总级的最大抑制量分别达到25.3 dB和15.6 dB; 此外, 前缘分流体对开孔腔流激振荡具有低频频移作用, 有益于避免发生流激空腔共振; 脉动压力空间分布特性表明, 前缘分流体对水中开孔腔流激振荡抑制机理主要在于破坏了腔内流场受到的周期调制作用.      

Controlling the near vortex wake of a swept half-wing by means of mini-flaps

The effect of mini-flaps on the flow pattern in the near vortex wake behind a model swept half-wing is investigated. The distributions of the time-average flow velocity were measured in a subsonic wind tunnel, in a section normal to the freestream velocity vector located at a distance of 3.8 wing half-spans from its trailing edge. When mini-flaps are mounted on both upper and lower wing surfaces, two vortices (tip and auxiliary) of the same sign are observable in the above-mentioned flow section; they are separated by an extended region of vorticity of the opposite sign. The model angle-of-attack effect on the intensities of the tip and auxiliary vortices is estimated.     

Large Eddy Simulation of boundary layer transition over an isolated ramp-type micro roughness element

Boundary layer transition over an isolated surface roughness element is investigated by means of numerical simulation. Large Eddy Simulation (LES) flow-modeling approach is employed to study flow characteristics and transition phenomenon past a roughness element immersed within an incoming developing boundary layer, at a height-based Reynolds number of 1170. LES numerical results are compared to experimental data from literature showing the time-averaged velocity distribution, the velocity fluctuation statistics and the instantaneous flow topology.Despite slight difference in the intensity of streamwise velocity fluctuations, the present LES results and experimental data show very good agreement. The mean flow visualization shows streamwise counter-rotating vortices pairs formation downstream of the obstacle. The primary pair induces an upwash motion and a momentum deficit that creates a Kelvin-Helmholtz type flow instability. The instantaneous flow topology reveals the formation of coherent K-H vortices downstream that produce turbulent fluctuations in the wake of the roughness element. These vortices are streched and lifted up when moving downstream. The velocity fluctuations results show that the onset of the turbulence is dominated by the energy transfer of large-scale vortices.     

Stability of a supersonic flat-plate wake (Comparison of numerical and experimental results)

The stability of the wake downstream of a plat plate with a blunt trailing edge is numerically investigated at a supersonic freestream velocity. The calculated stability characteristics are compared with the experimental data.     

Unsteady force measurements in sphere flow from subcritical to supercritical Reynolds numbers

The flow over a smooth sphere is examined in the Reynolds number range of 5.0 × 10⁴ < Re < 5.0 × 10⁵ via measurements of the fluctuating forces and particle image velocimetry measurements in a planar cut of the velocity field. Comprehensive studies of the statistics and spectra of the forces are presented for a range of subcritical and supercritical Reynolds numbers. While the subcritical lateral force spectra are dominated by activity corresponding to the large-scale vortex shedding frequency at a Strouhal number of approximately 0.18, there is no such peak apparent in the supercritical spectra, although resolution effects may become important in this region. Nor does the large-scale vortex shedding appear to have a significant effect on the drag force fluctuations at either sub- or super-critical Reynolds numbers. A simple double spring model is shown to capture the main features of the lateral force spectra. The low-frequency force fluctuations observed in earlier computational studies are shown to have important implications for statistical convergence, and in particular, the apparent mean side force observed in earlier studies. At least one thousand dimensionless time units are required for reasonable estimates of the second and higher moments below the critical Reynolds number and even more for supercritical flow, stringent conditions for computational studies. Lastly, investigation of the relationship between the motion of the instantaneous wake shape, defined via the local position where the streamwise velocity is equal to half the freestream value, and the in-plane lateral force for subcritical flow reveals a significant negative correlation throughout the near wake, which is shown to be related to a structure inferred to arise from the large-scale vortex shedding convecting downstream at 61% of the freestream velocity. In addition to its utility in understanding basic sphere flow, the apparatus is also a testbed that will be used in future studies, examining the effect of both static and dynamic changes to the surface morphology.     

Experimental study of flow past a square cylinder at high Reynolds numbers

 Measurements of two-components of velocity in the wake of a square cylinder using a hot-wire anemometer are reported. Two Reynolds numbers, namely 8700 and 17,625, have been considered. The measurements were carried out in a low-speed, low-turbulence wind tunnel. Benchmark experiments at much lower Reynolds numbers show good agreement between the present experiments and published results. At higher Reynolds numbers, the experimental data reveal anticipated trends in terms of wake recovery and turbulence decay. Both velocity and velocity fluctuations show symmetry about the wake axis. The experimental data have been compared with the large eddy simulation (LES) calculation reported by Wang et al. [University of Illinois at Urbana – Champaign (1996) Report CFD 96-03] and LDV measurements of Lyn et al. [J Fluid Mech (1995) 304: 285–319]. The agreement among the three sets is generally acceptable in terms of the time-averaged velocity components, but not the velocity fluctuations. The turbulence fluctuations in the present experiments are seen to be lower than in the referred work. The differences have been traced to factors such as the aspect ratio, blockage ratio and upstream turbulence. Experiments with increased upstream turbulence did show a reduction in the discrepancy between the present experiments and the published data. An assessment of the experimental data in terms of physical mechanisms revealed that (a) streamwise normal stresses were correlated with the vortex centers, and (b) the turbulence kinetic energy profiles are similar to the turbulence shear stress. Spectral analysis of the velocity signals was carried out in the present work. Energy transfer from the mean flow to the streamwise velocity fluctuation was confirmed in the near wake. A redistribution of the kinetic energy between the streamwise and transverse components of velocity over a longer distance downstream was subsequently observed. Received: 17 May 1999/Accepted: 29 December 1999     

Investigation on a coupled Navier–Stokes–Direct Simulation Monte Carlo method for the simulation of plume flowfield of a conical nozzle

To assess the plume effects of space thrusters, the accurate plume flowfield is indispensable. The plume flow of thrusters involves both continuum and rarefied flow regimes. Coupled Navier–Stokes–Direct Simulation Monte Carlo (NS–DSMC) method is a major approach to the simulation of continuum‐rarefied flows. An axisymmetric coupled NS–DSMC solver, possessing adaptive‐interface and two‐way coupling features, is investigated in this paper for the simulation of the nozzle and plume flows of thrusters. The state‐based coupling scheme is adopted, and the gradient local Knudsen number is used to indicate the breakdown of continuum solver. The nitrogen flows in a conical nozzle and its plume are chosen as the reference case to test the coupled solver. The threshold value of the continuum breakdown parameter is studied based on both theoretical kinetic velocity sampling and coupled numerical tests. Succeeding comparisons between coupled and full DSMC results demonstrate their conformities, meanwhile, the former saves 58.8% computational time. The pitot pressure evaluated from the coupled simulation result is compared with the experimental data proposed in literatures, revealing that the coupled method makes precise predictions on the experimental pitot pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

Particle resolved direct numerical simulation of a liquid–solid fluidized bed: Comparison with experimental data

Particle-resolved direct numerical simulations of a 3-D liquid–solid fluidized bed experimentally investigated by Aguilar-Corona (2008) have been performed at different fluidization velocities (corresponding to a range of bed solid volume fraction between 0.1 and 0.4), using Implicit Tensorial Penalty Method. Particle Reynolds number and Stokes number are O(100) and O(10), respectively. In this paper, we compare the statistical quantities computed from numerical results with the experimental data obtained with 3-D trajectography and High Frequency PIV. Fluidization law predicted by the numerical simulations is in very good agreement with the experimental curve and the main features of trajectories and Lagrangian velocity signal of the particles are well reproduced by the simulations. The evolution of particle and flow velocity variances as a function of bed solid volume fraction is also well captured by the simulations. In particular, the numerical simulations predict the right level of anisotropy of the dispersed phase fluctuations and its independence of bed solid volume fraction. They also confirm the high value of the ratio between the fluid and the particle phase fluctuating kinetic energy. A quick analysis suggests that the fluid velocity fluctuations are mainly driven by fluid–particle wake interactions (pseudo-turbulence) whereas the particle velocity fluctuations derive essentially from the large scale flow motion (recirculation). Lagrangian autocorrelation function of particle fluctuating velocity exhibits large-scale oscillations, which are not observed in the corresponding experimental curves, a difference probably due to a statistical averaging effect. Evolution as a function of the bed solid volume fraction and the collision frequency based upon transverse component of particle kinetic energy correctly matches the experimental trend and is well fitted by a theoretical expression derived from Kinetic Theory of Granular Flows.     

Simulation of the flow past a circular cylinder in the supercritical regime by blending RANS and variational-multiscale LES models

A strategy which blends a variational multiscale large eddy simulation (VMS-LES) model and a RANS model in a hybrid approach is investigated. A smooth blending function, which is based on the value of a blending parameter, is used for switching from VMS-LES to RANS. Different definitions of the blending parameter are investigated. The capabilities of the novel hybrid approach are appraised in the simulation of the flow around a circular cylinder at a Reynolds number 1.4×10⁵, based on the freestream velocity and on the cylinder diameter, in the presence of turbulent boundary-layer due to turbulent inflow conditions. A second study at Reynolds numbers from Re=6.7×10⁵ to 1.25×10⁶ is also presented. The effect of using the VMS-LES approach in the hybrid model is evaluated. Results are compared to those of other RANS, LES and hybrid simulations in the literature and with experimental data     

Large-eddy simulation of a turbulent boundary layer with blowing

The modifications of a turbulent boundary layer induced by blowing through a porous plate were investigated using large-eddy simulation. The Reynolds number (based on the length of the plate) of the main flow was about 850000. Large-eddy simulations of such a boundary layer needs a turbulent inflow condition. After a review of available turbulent inflow, we describe in details the condition we developed, which consisted of recycling the velocity fluctuations. Then we show the necessity for this inflow to be non-stationary and to be three dimensional with respect to the mass conservation equation. If these properties are not achieved, we found that the velocity fluctuations do not grow as expected along the domain. Finally, the results of simulations of the boundary layer submitted to blowing are compared with experimental measurements. The good agreement obtained validate our turbulent inflow conditions and also the blowing model used. PACS 47.27.Eq, 47.27.Te, 44.20.+b     

Navier-Stokes方程二阶速度滑移边界条件的检验

对微尺度气体流动,Navier-Stokes方程和一阶速度滑移边界条件的结果与实验数据相比,在滑移区相互符合,在过渡领域则显著偏离.为改善Navier-Stokes方程在过渡领域的表现,有些研究者尝试引入二阶速度滑移边界条件,如Cercignani模型,Deissler模型和Beskok-Karniadakis模型.以微槽道气体流动为例,将Navier-Stokes方程在不同的二阶速度滑移模型下的结果与动理论的直接模拟Monte Carlo(DSMC)方法和信息保存(IP)方法以及实验数据进行比较.在所考察的3种具有代表性的二阶速度滑移模型中,Cercignani模型表现最好,其所给出的质量流率在Knudsen数为0.4时仍与DSMC和IP结果相符;然而,细致比较表明,Cercignani模型给出的物面滑移速度及其附近的速度分布在滑流区和过渡领域的分界处(Kn=0.1)已明显偏离DSMC和IP的结果.     

A computational investigation of unsteady turbulent wake–boundary-layer interaction

The relative motion of rotor and stator blade rows in a turbomachine generates periodically unsteady flow on the blades due to travelling wake perturbations. To better understand the attendant wake–boundary-layer interaction a calculation procedure was developed to model the behaviour of this complex unsteady flow. Due to nonlinear interactions with the boundary layer, the travelling discrete frequency wakes were found to decrease the velocity profile shape factor. For the range of reduced frequencies examined (=0.33–9.33) the skin-friction coefficient was found to be frequency dependent. The calculated results for both steady and unsteady velocity profiles, and for skin friction compared well with experimental data. Although the agreement between measured and calculated velocity phase shift was poor, in both experimental and model results the negative phase shift throughout the boundary layer due to the travelling-wave fluctuations has been captured.     

Modelling of particle-wall collisions in confined gas-particle flows

This paper demonstrates that numerical simulations of confined particulate two-phase flows require a detailed modelling of particle—wall collisions which includes the wall surface structure and the particle shape. These effects are taken into account by “irregular bouncing” models which are based on the statistical treatment of the collision process. In this study, results obtained using various “irregular bouncing” models based on the impulse equations for a particle—wall collision are considered and compared with experimental observations. The wall roughness is simulated by assuming that the particle collides with a virtual wall which has a randomly distributed inclination with respect to the plane, smooth wall. A Gaussian distribution for this random inclination showed the best agreement with experimental results. Numerical predictions of a turbulent two—phase flow in a vertical channel, where the particle phase is treated using a Lagrangian approach, showed that the different models applied for a particle-wall collision have a strong effect on the particle velocity fluctuations and the mass flux profiles in the region of fully developed flow. The numerical simulations using the irregular bouncing models yielded considerably higher values for the particle velocity fluctuations, which also agreed better with the experimental values. This effect was most pronounced for large particles, where the distance they need to respond to the fluid flow is larger than the characteristic dimension of the confinement. On the other hand, the motion of small particles is less affected by the choice of the wall-collision model. These effects of the wall roughness on the velocity fluctuations of the dispersed phase have not been considered in previous studies using irregular bouncing models.     

URANS prediction of flow fluctuations in rod bundle with split-type spacer grid

Turbulent flow in a rod bundle with split-type spacer grid has been studied using Unsteady Reynolds-Averaged Navier–Stokes (URANS) approach. In the previous studies of turbulent flow in rod bundles URANS (as well as steady-state RANS) simulations predicted mean velocity profiles fairly well. However, they severely underpredicted velocity fluctuations, which is investigated in the present study. Our simulations were performed with the Shear Stress Transport (SST) turbulence model and automatic wall-treatment using OpenFOAM, an open-source CFD code. Results of URANS simulations are compared with the measurements of the MATiS-H experiment, which was performed at Korean Atomic Energy Research Institute (KAERI) in 2011–2012.The URANS predictions of velocity fluctuations have been improved by appropriately summing up fluctuations resolved by the basic URANS model and non-resolved fluctuations, which were modelled with the turbulence model. This treatment of turbulent fluctuations, which are directly measured in high-quality experiments, allows more detailed evaluation of various URANS turbulence models. It was found out that the best agreement is achieved when resolved and modelled fluctuations are assumed to be uncorrelated, which indicates that the large-scale structures in this particular flow are distinct in the spectral space from the rest of turbulence. Turbulent flow in the MATiS-H experiment was reproduced by numerous authors using different approaches and our results are among the most accurate.