Flow visualization and scanning PIV measurement of three-dimensional structure in near field of strongly buoyant jet

Abstract  

The near-field structure of strongly buoyant jet issuing from a square nozzle at low Froude and Reynolds numbers is studied by using LIF flow visualization and time-resolved scanning PIV. These experimental techniques allow the visualization of unsteady three-dimensional flow phenomenon occurring in the near-field of strongly buoyant jet. It is found that the buoyant jet is unstable to the positive buoyancy forces, which promote the inflow motion near the nozzle exit. The surrounding low temperature fluid moves into the nozzle inside along the nozzle corner and mixes with the high temperature fluid deep into the nozzle. Then, the flow pattern inside the nozzle becomes highly complex to promote the laminar to turbulent transition of the jet. The statistical flow characteristics of the strongly buoyant jet are evaluated from the scanning PIV measurement, and the result indicates the presence of axisymmetric distributions of mean flow and velocity fluctuations in the circle of diameter equal to the square side of the nozzle.     

Investigation of the pseudo-shock wave in a two-dimensional supersonic inlet

Abstract  

This paper describes experimental and numerical investigations into the multiple shock waves/turbulent boundary layer interaction in a supersonic inlet. The test model has a rectangular shape with an asymmetric subsonic diffuser of 5°. Experiments were conducted to obtain the visualization images and static pressure data by using supersonic wind tunnel. Numerical simulation was performed by solving the RANS equations with the Menter’s SST turbulent model. The inflow condition was a free-stream Mach number of 2.5 and a unit Reynolds number of 7.6 × 10⁷/m. Numerical results showed good agreement with the experimental results. Based on this agreement, the flow characteristics which are often very difficult to obtain experimentally alone were analyzed with the aid of numerical simulation. The structures, pressure and velocity distributions, and total pressure loss of the pseudo-shock wave in the supersonic inlet were presented in detail from flow visualization images and static pressures.     

PIV measurements of the flow and turbulent characteristics of a round jet in crossflow

The instantaneous and ensemble averaged flow characteristics of a round jet issuing normally into a crossflow was studied using a flow visualization technique and Particle Image Velocimetry measurements. Experiments were performed at a jet-to-crossflow velocity ratio, 3.3 and two Reynolds numbers, 1,050 and 2,100, based on crossflow velocity and jet diameter. Instantaneous laser tomographic images of the vertical center plane of the crossflow jet show that there exists very different natures in the flow structures of the near field jet due to Reynolds number effect even though the velocity ratio is same. It is found that the shear layer becomes much thicker when the Reynolds number is 2,100 because of the strong entrainment of the inviscid fluid by turbulent interaction between the jet and crossflow. The mean and second order statistics are calculated by ensemble averaging over 1,000 realizations of instantaneous velocity fields. The detail characteristics of mean flow field, streamwise and vertical rms velocity fluctuations, and Reynolds shear stress distributions are presented. The new PIV results are compared with those from previous experimental and LES studies.     

Large eddy simulation of turbulent horizontal buoyant jets

Large eddy simulations (LESs) of turbulent horizontal buoyant jets are carried out using a high-order numerical method and Sigma subgrid-scale (SGS) eddy-viscosity model, for a number of different Reynolds (Re) and Richardson (Ri) numbers. Simulations at previous experimental flow conditions (Re = 3200, 24, 000 and Ri = 0, 0.01) are carried out first, and the results are found to be qualitatively and quantitatively similar to the experimental results, thus validating the numerical methodology. The effect of varying Ri (values 2×10^?4, 0.001, 0.005, and 0.01) and Re (3200 and 24, 000) is studied next. The presence of stable stratification on one side and unstable stratification on the other side of the jet centreline leads to an asymmetric development of horizontal buoyant jets. It is found that this asymmetry, the total radial spread and the vertical deflection are significantly affected by Ri, while Re affects only the radial asymmetry. The need for developing improved integral models, accounting for this asymmetry, is pointed out. Turbulent production and dissipation rates are investigated, and are found to be symmetric in the horizontal plane, but asymmetric in the mid-vertical plane. A previously proposed model, for correlation between the vertical component of the fluctuating scalar flux vector and the vertical cross-correlation component of the Reynolds tensor, is modified based on the current LES results. Instantaneous scalar and velocity fields are analysed to reveal the structure of horizontal buoyant jets. Similar to the developed turbulent jet, the flow close to the nozzle too is found to be markedly different in the stable and unstable stratification regions. Persistent coherent vortex rings are found in the stable stratification region, while intermittent breakdown of vortex rings into small-scale structures is observed in the unstable stratification region. Similarities and differences between the flow structures in the horizontal buoyant jet configuration and those in the jet in crossflow configuration are discussed. Finally, a dynamic mode decomposition analysis is carried out, which indicates that the flow in the unstable stratification region is more energetic and prone to instabilities, as compared to the flow in the stable stratification region.     

Numerical simulation and flow visualization using soap film of the self-organized vortex structure in the wake of an array of cylinders

Abstract  

With the aid of computational fluid dynamics (CFD) and simple flow visualization technique using flowing soap-film, we present here the wake structures behind an array of cylinders for Reynolds numbers corresponding to both laminar and turbulent flow regimes. The image results illustrate interesting vortex interactions past these equally spaced cylinders; for low Reynolds number flow, well-organized wake pattern persists and manifests unsteadily to different symmetry states. An increase of free stream flow velocity causes the wake transition, resulting in the formation of asymmetric flow wake with chaotic mixing at the far wake. Observations from both the numerical simulations and soap-film are in good agreement at least qualitatively.     

Observation of vortex packets in direct numerical simulation of fully turbulent channel flow

A number of experimental studies have inferred the existence of packets of inclined, hairpinlike vortices in wall turbulence on the basis of observations made in two-dimensional x−y planes using visualization and particle image velocimetry (PIV). However, there are very few observations of hairpins in existing three-dimensional studies made using direct numerical simulation (DNS), and no such study claims to have revealed packets. We demonstrate, for the first time, the existence of hairpin vortex packets in DNS of turbulent flow. The vortex packet structure found in the present study at low Reynolds number,Re _t=300, is consistent with and substantiates the observations and the results from twodimensional PIV measurements at higher Reynolds numbers in channel, pipe and boundary layer flows. Thus, the evidence supports the view that vortex packets are a universal feature of wall turbulence, independent of effects due to boundary layer trips or critical conditions in the aforementioned numerical studies. Visualization of the DNS velocity field and vortices also shows the close association of hairpin packets with long low-momentum streaks and the regions of high Reynolds shear stress.     

圆柱绕流问题的初期流动数值模拟

本文继Zheng-huan Teng[1]中介绍的解Navier-stokes方程的椭圆涡团法,研究了一种新的变形涡团法,用以模拟不可压粘性流体绕圆柱的不定常流动。圆柱在静止流体中突然起动并做匀速直线运动。对整个流动区域构造完全Navier-Stokes方程的解并不容易,近十年出现很多数值研究,本文对算法有所推进。把本文的方法称作变形涡团法是因为圆柱边界附近的流体中用椭圆涡团,远离边界时用圆形涡团。计算圆柱绕流比平板绕流在满足附着条件上更为困难,本文分析了怎样在圆柱边界上给出适当的附着条件的数值方法。在算例中雷诺数分别取200、550、3000,得到了不定常边界层分离,二次涡等复杂的物理现象,这些数值结果与近年实验结果[2]是一致的。     

Flame-acoustic coupling of combustion instability in a non-premixed backward-facing step combustor: the role of acoustic-Reynolds stress

Combustion instability in a laboratory scale backward-facing step combustor is numerically investigated by carrying out an acoustically coupled incompressible large eddy simulation of turbulent reacting flow for various Reynolds numbers with fuel injection at the step. The problem is mathematically formulated as a decomposition of the full compressible Navier–Stokes equations using multi-scale analysis by recognising the small length scale and large time scale of the flow field relative to a longitudinal mode acoustic field for low mean Mach numbers. The equations are decomposed into those for an incompressible flow with temperature-dependent density to zeroth order and linearised Euler equations for acoustics as a first order compressibility correction. Explicit coupling terms between the two equation sets are identified to be the flow dilatation as a source of acoustic energy and the acoustic Reynolds stress (ARS) as a source of flow momentum. The numerical simulations are able to capture the experimentally observed flow–acoustic lock-on that signifies the onset of combustion instability, marked by a shift in the dominant frequency from an acoustic to a hydrodynamic mode and accompanied by a nonlinear variation of pressure amplitude. Attention is devoted to flow conditions at two Reynolds numbers before and after lock-on to show that, after lock-on, the ARS causes large-scale vortical rollup resulting in the evolution of a compact flame. As compared to acoustically uncoupled simulations at these Reynolds numbers that show an elongated flame with no significant roll up and disturbance in the upstream flow field, the ARS is seen to alter the shear layer dynamics by affecting the flow field upstream of the step as well, when acoustically coupled.     

低速离心压气机的全工况数值模拟和校核

本文选用不同的计算模型和网格数目在全流量域对NASA低速离心压气机进行了数值模拟,通过与实验数据的对比检验了计算结果的精度,为具有无叶扩压器结构的低压离心压气机确定了合理的计算方案。同时通过计算结果的观察和对比,分析了小流量下叶轮内部流场的部分变化特征。     

Computational study of the laminar to turbulent transition over the SD7003 airfoil in ground effect

In this work, we present a numerical study of the laminar-turbulence transition flow around a symmetrical air-foil at a low Reynolds number in free flow and near the ground surface at different angles of attack. Finite volume method is employed to solve the unsteady Reynolds-averaged Navier–Stokes (RANS) equation. In this way, the Transition SST turbulence model is used for modeling the flow turbulence. Flow around the symmetrical airfoil SD7003 is numerically simulated in free stream and near the ground surface. Our numerical method can detect different aspects of flow such as adverse pressure gradient, laminar separation bubble and laminar to turbulent transition onset and the numerical results are in good agreement with the experimental data.     

DNS study of swirling intensity effect on flow pattern of a circular jet

Abstract  

The figures show the 3D flow pattern of a circular jet with different swirling intensity. Reynolds number is approximately 4300 computed based on the nozzle diameter (d), jet velocity (U), and air fluid property at 1 atm and 300 K. The overall computational domain is set to be 4 × 4 × 12 d in spanwise, height, and streamwise direction. The governing equations are the fully compressible Navier–Stokes equations, firstly differenced by eighth-order explicit scheme and then advanced temporarily by using the fourth-order explicit Runge–Kutta method. 3D characteristics non-reflecting boundary condition including transverse source contribution is imposed on all other boundaries except the inflow boundary handled by assigning fixed profiles of temperature and velocity. To ensure the simulation resolution, here over 16 million grids are employed in sum, combined with a handful of grids located at buffer zones of outflow boundaries. To correctly represent the vortex in the flow field, velocity gradient tensor invariant Q is used here. And ψ refers to the swirling intensity defined as the ratio of tangential momentum to axis momentum. As shown in velocity profile, the flow pattern of the jet changes from a close mode to a totally open mode as ψ increases from 0.4 to 1.5. Accordingly, the recirculation zone gradually moves upstream and backflow velocity is enlarged as well. It is inteseting to found that the obvious drops of the momentums in two shown directions always occur at the same position downstream, no matter how large the ψ value is. Therefore, a momentum compensatory mechanism is expected to exist in the vortex-abundant zone. With the increase of ψ value, the increased strain rate in tangential direction can induce vortex more quickly, intensifying the entrainment and velocity-attenuation, which can be observed in Q value profile.     

Experimental Investigation of Opposing Mixed Convection in a Channel with an open Cavity Below

Abstract

In this article, mixed convection in an open cavity with a heated wall bounded by a horizontal unheated plate is investigated experimentally. The heated wall is on the opposite side of the forced inflow. The results are reported in terms of wall temperature profiles of the heated wall and flow visualization. The range of pertinent parameters used in this experiment are Reynolds numbers (Re) from 100 to 2,000 and Richardson numbers (Ri) from 4.3 to 6,400. Also, the ratio between the length and the height of cavity (L/D) ranges from 0.5–2.0, and the ratio between the channel and cavity height (H/D) is equal to 1.0. The lack of experimental results on mixed convection in a channel with an open cavity below was an impetus for investigating this configuration when one cavity vertical wall is heated at uniform heat flux. The present results show that at the lowest investigated Reynolds number, the surface temperatures are lower than the corresponding surface temperatures for Re = 2,000 at the same ohmic heat flux. The flow visualization shows that for Re = 1,000, there are two nearly distinct fluid motions: a parallel forced flow in the channel and a recirculation flow inside the cavity. For Re = 100, the effect of a stronger buoyancy determines a penetration of thermal plumes from the heated plate wall into the upper channel. Moreover, the flow visualization shows that for lower Reynolds numbers, the forced motion penetrates inside the cavity, and a vortex structure is adjacent to the unheated vertical plate. At higher Reynolds numbers, the vortex structure has a larger extension while L/D is held constant.     

Direct numerical simulation of the flow around a wall-mounted square cylinder under various inflow conditions

The flow around a wall-mounted square cylinder of side d is investigated by means of direct numerical simulation (DNS). The effect of inflow conditions is assessed by considering two different cases with matching momentum-thickness Reynolds numbers Re_θ ? 1000 at the obstacle: the first case is a fullyturbulent zero pressure gradient boundary layer, and the second one is a laminar boundary layer with prescribed Blasius inflow profile further upstream. An auxiliary simulation carried out with the pseudo-spectral Fourier–Chebyshev code SIMSON is used to obtain the turbulent time-dependent inflow conditions which are then fed into the main simulation where the actual flow around the cylinder is computed. This main simulation is performed, for both laminar and turbulent-inflows, with the spectral-element method code Nek5000. In both cases the wake is completely turbulent, and we find the same Strouhal number St ? 0.1, although the two wakes exhibit structural differences for x > 3d downstream of the cylinder. Transition to turbulence is observed in the laminar-inflow case, induced by the recirculation bubble produced upstream of the obstacle, and in the turbulent-inflow simulation the streamwise fluctuations modulate the horseshoe vortex. The wake obtained in our laminar-inflow case is in closer agreement with reference particle image velocimetry measurements of the same geometry, revealing that the experimental boundary layer was not fully turbulent in that dataset, and highlighting the usefulness of DNS to assess the quality of experimental inflow conditions.     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

Direct numerical simulation of transitional flow in a finite length curved pipe

Transition to turbulent flow in a curved pipe has been well studied through experiments and numerical simulations. Numerical simulations often use a helical pipe with an infinite length such that the inlet and outlet boundary conditions can be modelled as periodic which greatly reduces computational time. In this study, we examined a finite length curved pipe with Poiseuille flow imposed at the inlet and a stress-free boundary condition at the outlet. Direct numerical simulation of the Navier-Stokes equations for rigid walls and a Newtonian fluid was performed using nek5000. Straight extensions were added to the inlet and outlet such to diminish the impact of boundary conditions on the flow field in the region with curvature. The examined model has a pipe radius of curvature that is three times the pipe radius. The model has ～355 million nodes and required an order of magnitude greater computational time when compared with an infinite length curved pipe. Results show that the critical Reynolds number, the lowest value with instabilities present in the flow, is much greater than that of a straight pipe and occurs near Re=5000–5200. This is larger than the critical Reynolds number typically reported for an infinite length curved pipe (Re=4200–4300).     

Application of RANS-based method to predict acoustic noise of chevron nozzles

Passive noise control devices for jet flows, such as chevron nozzles, have been studied for a long time due to their large application in turbofan engines. The main purpose of their usage is the reduction of low-frequency noise generation and thus decreasing the noise perceived at the far field. This work is a numerical study of acoustic noise generated by jet flow operating at Mach number 0.9 and Reynolds number 1.38 × 10⁶, considering two chevron nozzle geometries that differ from each other by the penetration angle into the flow. The main aim was to demonstrate that Reynolds averaged Navier Stokes (RANS)-based methods are reliable means to obtain acoustical noise predictions for the industry with a considerably low computational cost. In order to achieve this objective, computational fluid dynamics (CFD) RANS simulations were performed with a cubic k-ɛ model and the acoustic noise spectrum for different angles of radiation was obtained via the Lighthill ray-tracing (LRT) method. The numerical results for the acoustic and flow fields were seen to be in reasonable agreement with the experimental data, suggesting that this methodology can be used as a fast and useful option to predict acoustic noise of jet flows from chevron nozzles.     

A numerical study of rotation effects on jets effusing from inclined holes into a cross-flow

In this work, the complexity of the flow field arising from the impact of the interaction of coolant jets with a hot cross-flow under rotation conditions was numerically simulated using large eddy simulation with artificial inflow boundary condition. The finite-volume method and the unsteady PISO (Pressure Implicit with Splitting of Operators) algorithm were applied on a non-uniform staggered grid. The simulations were performed for four different values of rotation number (Ro) of 0.0, 0.03021, 0.06042, and 0.12084, a jet Reynolds number of 4700, based on the hole width and the jet exit velocity. The air jet was injected at 30° and 90° in the streamwise direction with a density ratio of 1.04 and a velocity ratio of 0.5. The flow fields of the present study were compared with experimental data in order to validate the reliability of the LES technique. It was shown that the rotation has a strong impact on the jet trajectory behaviour and the film cooling effectiveness. The film trajectory always inclines centrifugally. Under rotating conditions, the film trajectory departs from the centreline to the left boundary. The deflection becomes greater as Ro increases. Furthermore, it was also found that the injection angle has a strong impact on separation and reattachment behaviour as well as the strength of the penetration into the cross-flow. As it increases, the distribution of the film cooling downstream the jet exit is more non-uniform and the film cooling effectiveness level slightly decreases.     

Arterial stenosis murmurs: an analysis of flow and pressure fields

The flow field distal to an arterial stenosis is simulated by a confined turbulent jet with moderate Reynolds numbers. The wall pressure fluctuations are related to the momentum fluctuations of the jet by the Poisson equation. A Green's function was derived to satisfy the boundary conditions on a cylindrical surface. This allows the solution of the Poisson's equation to include only a volume integral of the fluctuating momentum, weighed by the relative distance between the source and the sensor. The velocity fluctuations on the jet centerline and at the middle of the shear layer were measured using a laser Doppler anemometer. The wall pressure fluctuations were detected simultaneously by an array of nine wall-mounted pressure transducers along the axial direction. Cross correlation performed between the velocity and pressure fluctuations reveals that the pressure fluctuations were mostly imposed by the passage of turbulent eddies with a convective velocity that is a function of the jet exit velocity. The Strouhal number, defined by the frequency of the passing large-scale structure, is a function of the initial conditions only very close to the jet exit. Further downstream, where the effect of the initial conditions is lost, the Strouhal number approaches a constant irrespect of the jet Reynolds number. The contribution of a source near the jet exit to wall pressure fluctuation near the reattachment is rather weak due to the rapidly decaying weighting function in the axial direction. However, for sources located within one nozzle diameter from the sensor, the cross-spectral density function has a high magnitude with maximum coherence where the pressure spectral changes its slope.     

出口雷诺数对平面射流自保持性的影响

通过实验研究出口雷诺数对平面湍流射流自保持性的影响. 测量的射流来自相同的喷嘴但不同的雷诺数Re（≡U_jh/ν,其中U_j是出口平均速度、h是窄缝出口的厚度和ν是黏性系数）,其变化范围是Re=4582—57735.所得的数据包括沿轴线的平均速度、湍流强度、积分尺度、高阶矩和能谱. 实验发现,随着Re的增大,平面射流发展减慢,平均速度和湍流强 关键词： 平面射流 雷诺数 自保持性     

Experimental investigation on laminar separation control for flow over a two-dimensional bump

The laminar boundary layer separation flow over a two-dimensional bump controlled by synthetic jets is experimentally investigated in a water channel with hydrogen-bubble visualisation and particle image velocimetry (PIV) techniques. The two-dimensional synthetic jet is applied near the separation point. Two Reynolds numbers (Re = 700 and 1120) based on the bump height and free-stream velocity are adopted in this experiment, and seven different excitation frequencies at each Reynolds number are considered, focusing on the separation control as well as the vortex dynamics. The experimental results show that the optimal control can only be achieved within some excitation frequencies at both Reynolds numbers. However, beyond this range, further increasing the excitation frequency leads to an increase in the separation region. The proper orthogonal decomposition (POD) technique and vortex identification by swirling strength (Λ_ci) are applied for the deeper analysis of the separated flow. The reconstructed Λ_ci field by the first four POD modes is used and vortex lock-on phenomenon is observed. It is found that the negative synthetic jet vortex with clockwise rotation draws the separated wake shear layer as it is convected downstream, and then they syncretise together. Thus, the new vortex is induced and shedding downstream periodically.