Experimental investigation of supersonic flow over elliptic surface

The coherent structures of flow over a compression elliptic surface are experimentally investigated in a supersonic low-noise wind tunnel at Mach Number 3 using nano-tracer planar laser scattering (NPLS) and particle image velocimetry (PIV) techniques. High spacial resolution images and the average velocity profiles of both laminar inflow and turbulent inflow over the testing model were captured. From statistically significant ensembles, spatial correlation analysis of both cases is performed to quantify the mean size and orientation of large structures. The results indicate that the mean structure is elliptical in shape and structure angles in separated region of laminar inflow are slightly smaller than that of turbulent inflow. Moreover, the structure angle of both cases increases with its distance away from from the wall. POD analysis of velocity and vorticity fields is performed for both cases. The energy portion of the first mode for the velocity data is much larger than that for the vorticity field. For vorticity decompositions, the contribution from the first mode for the laminar inflow is slightly larger than that for the turbulent inflow and the cumulative contributions for laminar inflow converges slightly faster than that for turbulent inflow     

Flow control of micro-ramps on supersonic forward-facing step flow

The effects of the micro-ramps on supersonic turbulent flow over a forward-facing step(FFS) was experimentally investigated in a supersonic low-noise wind tunnel at Mach number 3 using nano-tracer planar laser scattering(NPLS)and particle image velocimetry(PIV) techniques. High spatiotemporal resolution images and velocity fields of supersonic flow over the testing model were captured. The fine structures and their spatial evolutionary characteristics without and with the micro-ramps were revealed and compared. The large-scale structures generated by the micro-ramps can survive the downstream FFS flowfield. The micro-ramps control on the flow separation and the separation shock unsteadiness was investigated by PIV results. With the micro-ramps, the reduction in the range of the reversal flow zone in streamwise direction is 50% and the turbulence intensity is also reduced. Moreover, the reduction in the average separated region and in separation shock unsteadiness are 47% and 26%, respectively. The results indicate that the micro-ramps are effective in reducing the flow separation and the separation shock unsteadiness.     

激波与层流/湍流边界层相互作用实验研究

在超声速风洞中,分别对层流和湍流来流条件下的边界层和斜激波(激波强度足以引起流动分离)相互干扰进行了实验研究,利用纳米粒子示踪平面激光散射(NPLS)技术获得了两种条件下流场的精细结构图像;利用粒子图像测速(PIV)技术获得了两种条件下流场的速度场和涡量场;综合运用NPLS结果和PIV结果对比分析了两种流动的瞬时流动结构和时间相关性,实验结果表明:层流边界层内的分离区呈现出狭长的条状,而湍流边界层内分离区呈现出较规则的椭圆;在入射激波上游距入射点较远的位置,层流边界层外围拟序结构会诱导出一系列压缩波系,进而汇聚成空间位置不稳定的诱导激波,而湍流边界层则是在入射激波上游较近的地方直接形成较强且稳定的诱导激波;在入射激波下游,层流边界层内的膨胀区域较小且急促,膨胀后产生的再附激波很弱,而湍流边界层内的膨胀区域较大,膨胀后产生的激波较强。     

超声速层流/湍流压缩拐角流动结构的实验研究

在Ma=3.0的超声速风洞中, 分别对上游边界层为超声速层流和湍流, 压缩角度为25°和28°的压缩拐角流动进行了实验研究. 采用纳米粒子示踪平面激光散射(NPLS)技术获得了流场整体和局部区域的精细结构, 边界层、剪切层、分离激波、回流区和再附激波等典型结构清晰可见, 测量了超声速层流压缩拐角壁面的压力系数. 从时间平均的流场结构中测量出分离激波、再附激波的角度和再附后重新发展的边界层的增长情况, 通过分析时间相关的流场NPLS图像, 可以发现流场结构随时间的演化特性. 实验结果表明: 在25°的压缩角度下, 超声速层流压缩拐角流动发生了典型的分离, 边界层迅速增长失稳转捩, 并引起一道诱导激波, 流场中出现了K-H涡、剪切层和微弱压缩波结构, 而超声速湍流压缩拐角流动没有出现分离, 湍流边界层始终表现为附着状态; 在28° 的压缩角度下, 超声速层流压缩拐角流动进一步分离, 回流区范围明显扩大, 诱导激波、分离激波向上游移动, 再附激波向下游移动, 分离区流动结构复杂, 相比之下, 超声速湍流压缩拐角流动的回流区范围明显较小, 边界层增长缓慢, 流场中没有出现诱导激波、K-H涡和压缩波, 流动分离区域的结构也相对简单, 但分离激波的强度则明显更强. 关键词： 压缩拐角 层流 湍流 流动结构     

Instantaneous and time-averaged flow structures around a blunt double-cone with or without supersonic film cooling visualized via nano-tracer planar laser scattering

In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering (NPLS), which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 μs revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H (Kelvin-Helmholtz) vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward-facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.     

带喷流超声速后台阶流场精细结构及其运动特性研究

采用基于纳米示踪的平面激光散射技术(NPLS)对带超声速喷流的后台阶流动精细结构进行了研究. 来流马赫数为3.4, 喷流实测马赫数为2.45, 而名义马赫数为2.5. 结果清晰地揭示了激波、剪切层、混合层、Kelvin-Helmholtz涡、羊角涡及湍流大尺度结构等大量典型流场结构. 基于大量流场精细结构图像, 对典型位置处的结构进行了空间两点相关性分析, 在喷流混合层前端涡结构小于湍流充分发展的尾端, 结构角相对也小. 喷流工作时, 模型台阶下游表面由一薄层气膜覆盖. 获得了模型流向和不同高度展向平面内的流场结构, 对照纹影试验结果, 分析了流动特点及时间演化规律. 采用微型压力扫描系统测试了模型表面的压力系数分布, 靠近喷流下游处压力系数区域0.0146. 针对NPLS图像做了流动的分形维数的分析, 发现在流动初始阶段分形维数接近于1, 越靠下游分形维数越高.     

Schlieren “PIV” for turbulent flows

The possibility of using commercial PIV equipment combined with schlieren optics to measure the velocity fields of turbulent flows is explored. Given a sufficiently high Reynolds number and adequate refractive flow differences, turbulent eddies can serve as the PIV “particles” in a schlieren image or shadowgram. The PIV software analyzes motion between consecutive schlieren or shadowgraph frames to obtain velocity fields. Velocimetry examples of an axisymmetric sonic helium jet in air and a 2D turbulent boundary layer at Mach 3 are shown. Due to optical path integration, axisymmetric flows require the inverse Abel transform to extract center-plane velocity data. Conditions for optimum schlieren sensitivity are examined. In its present embodiment, “schlieren PIV” is not useful for laminar flows nor for fully 3D flows. Otherwise it functions much like standard PIV under conditions where individual particles are not resolved and velocimetry is instead based on correlation of the motion of turbulent structures. “Schlieren PIV” shows significant promise for general refractive turbulent flow velocimetry if its integrative nature can be overcome through sharp-focusing optics.     

Simultaneous density and velocity measurements in a supersonic turbulent boundary layer

A novel technique for simultaneous measurements of instantaneous whole-field density and velocity fields of supersonic flows has been developed.The density measurement is performed based on the nano-tracer planar laser scattering(NPLS) technique,while the velocity measurement is carried out using particle image velocimetry(PIV).The present experimental technique has been applied to a flat-plate turbulent boundary layer at Mach 3,and the measurement accuracy of the density and velocity are discussed.Based on this new technique,the Reynolds stress distributions were also obtained,demonstrating that this is an effective means for measuring Reynolds stresses under compressible conditions.     

Dynamic analysis of bubble-driven liquid flows using time-resolved particle image velocimetry and proper orthogonal decomposition techniques

Abstract  

An experimental study to evaluate dynamic structures of flow motion and turbulence characteristics in bubble-driven water flow in a rectangular tank with a varying flow rate of compressed air is conducted. Liquid flow fields are measured by time-resolved particle image velocimetry (PIV) with fluorescent tracer particles to eliminate diffused reflections, and by an image intensifier to acquire enhanced clean particle images. By proper orthogonal decomposition (POD) analysis, the energy distributions of spatial and temporal modes are acquired. Time-averaged velocity and turbulent kinetic energy distributions are varied with the air flow rates. With increasing Reynolds number, bubble-induced turbulent motion becomes dominant rather than the recirculating flow near the side wall. Detailed spatial structures and the unsteady behavior of dominant dynamic modes associated with turbulent kinetic energy distributions are addressed.     

Direct numerical simulation of a fully developed compressible wall turbulence over a wavy wall

Compressible turbulent channel flow over a wavy surface is investigated by direct numerical simulations using high-resolution finite difference schemes. The Reynolds number considered in the present paper is 3380 based on the bulk velocity, the channel half-width and the kinetic viscosity at the wall. Four test cases are simulated and analysed at Ma_m = 0.33, 0.8, 1.2, 1.5 based on the bulk velocity and the speed of sound at the wall. We mainly focus on the curvature and the Mach number effects on the compressible turbulent flows. Numerical results show that although the wavy wall has effects on the mean and fluctuation quantities, log law still exists in the distribution of the wave-averaged streamwise velocity if the roughness effects are taken into consideration in the scaling of it. Near-wall streaks are broken by the wavy surface and near-wall quasi-streamwise vortices mostly begin at the upslope of the wave and pass over the crest of it. The wavy wall makes the turbulence more active and the flow easier to be blended. From the viewpoint of turbulent kinetic budgets, curvature effects strengthen both the diffusion terms and the dissipation terms. At the same time, they change the properties of the compressibility-related terms and promote more inner energy transferring into turbulent kinetic energy. As the Mach number increases, the reattachment of the mean flow is delayed, which indicates the mean separation bubble becomes larger. Concerning the near-wall coherent structures, the vortices are more sparsely distributed with the increasing of the Mach number. For the supersonic cases, shock waves appear. Though they have little effects on the mean turbulent quantities, they change the structures of the flow fields and induce local separations at the upper wall of the channel.     

Stereoscopic PIV measurements of flow around a marine propeller

A stereoscopic PIV (Particle Image Velocimetry) technique has been employed to measure the 3 dimensional flow structure of turbulent wake behind a marine propeller with 5 blades. The out-of-plane velocity component was measured using particle images captured simultaneously by two CCD cameras installed in the angular displacement configuration. 400 instantaneous velocity fields were acquired for each of four different blade phases of 0°, 18°, 36° and 54°. They were ensemble averaged to investigate the spatial evolution of propeller wake in the near wake region up to one propeller diameter (D) downstream. The phase-averaged velocity fields show clearly the viscous wake formed by the boundary layers developed along both surfaces of the blade. Tip vortices were generated periodically and the slipstream contraction occurs in the near-wake region. The out-of-plane velocity component has large values at the locations of tip and trailing vortices. With going downstream, the axial turbulence intensity and the strength of tip vortices were decreased due to the viscous dissipation, turbulence diffusion and blade-to-blade interaction. The difference in the mean velocity fields measured by SPIV and 2-D PIV methods was about 5% ≈ 10%. However, the 2-D PIV results also give sufficient information on propeller wake beyond the region of X/D=0.2.     

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.     

Convection and correlation of coherent structure in turbulent boundary layer using tomographic particle image velocimetry

The present experimental work focuses on a new model for space–time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wall. A turbulent boundary layer flow at Reθ= 2460 is measured by tomographic particle image velocimetry(tomographic PIV). It is demonstrated that arch, cane,and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space–time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space–time correlation instead of Taylor hypothesis. The convection velocities derived from the space–time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition(FD) shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition(POD) has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.     

A study on large coherent structures and noise emission in a turbulent round jet

A moderate Reynolds number,and high subsonic turbulent round jet is investigated by large eddy simulation.The detailed results(e.g.mean flow properties,turbulence intensities,etc.)are validated against the experimental data,and special attention is paid to study motions of coherent structures and their contributions to far-field noise.Eulerian methods(e.g.Q-criteria andλ2criteria)are utilized for visualizing coherent structures directly for instantaneous flow fields,and Lagrangian coherent structures accounting for integral effect are shown via calculating fields of finite time Lyapunov exponents based on bidimensional velocity fields.All visualizations demonstrate that intrusion of three-dimensional vortical structures into jet core occurs intermittently at the end of the potential core,resulting from the breakdown of helical vortex rings in the shear layer.Intermittencies in the shear layer and on the centerline are studied quantitatively,and distinctively different distributions of probability density function are observed.Moreover,the physical sound sources are obtained through a filtering operation of defined sources in Lighthill’s analogy,and their distributions verify that intrusion of vortical structures into the core region serves as important sound sources,in particular for noise at aft angles.The facts that intermittent behaviors are caused by motions of coherent structures and correlated with noise generation imply that to establish reasonable sound sources in active noise production region based on intermittent coherent structures is one of the key issues for far-field noise prediction.     

Three-dimensional turbulent flow over cube-obstacles

In order to investigate the influence of surface roughness on turbulent flow and examine the wall-similarity hypothesis of Townsend, three-dimensional numerical study of turbulent channel flow over smooth and cube-rough walls with different roughness height has been carried out by using large eddy simulation(LES) coupled with immersed boundary method(IBM). The effects of surface roughness array on mean and fluctuating velocity profiles, Reynolds shear stress, and typical coherent structures such as quasi-streamwise vortices(QSV) in turbulent channel flow are obtained. The significant influences on turbulent fluctuations and structures are observed in roughness sub-layer(five times of roughness height).However, no dramatic modification of the log-law of the mean flow velocity and turbulence fluctuations can be found by surface cube roughness in the outer layer. Therefore, the results support the wall-similarity hypothesis. Moreover, the von Karman constant decreases with the increase of roughness height in the present simulation results. Besides, the larger size of QSV and more intense ejections are induced by the roughness elements, which is crucial for heat and mass transfer enhancement.     

Visualization of the inflow ahead of a rotating propeller attached to a container ship model

The inflow ahead of a rotating propeller attached to a container ship model was visualized using a two-frame particle image velocimetry (PIV) technique. For illuminating the inflow region, a transparent window was installed at the stern of the ship model. Ensemble-averaged mean velocity fields were measured at four different blade phases under the design loading condition. The characteristics of the inflow in the upper plane above the propeller axis are quite different from those below the propeller axis. In the far upstream region above the propeller axis, most of the inflow comes from the hull wake and the axial velocity is very small. As the inflow moves toward the propeller plane, its axial velocity component increases rapidly. In addition, the variation of the inflow characteristics with respect to phase angle becomes apparent. The thick hull boundary layer and out-of-plane motion resulting from the propeller rotation produce a large turbulent kinetic energy around the tip of the propeller blade in the upper inflow region. The axial velocity distribution of the propeller inflow is asymmetric with respect to the vertical center axis, exhibiting different axial velocities on the port and starboard sides.     

Visualisation on supersonic flow over backward-facing step with or without roughness

At Mach number 3.4, visualisation experiments of flow over backward-facing step (BFS) with or without roughness band attached on upstream wall are carried out via traditional schlieren and newly developed nano-tracer-based planar laser scattering (NPLS). The time-averaged flow characteristic of the reattachment region and the instantaneous rich structures of the redeveloping boundary layer in the steamwise-normal plane are both revealed. Additionally, top views in the different planes (y/h = 0.67, 1.00, 1.33, 1.67, 2.00) are imaged with a resolution of 0.064 mm/pixel. By contrasting the NPLS images at different times, the unsteady evolution characteristic of the coherent vortices in the redeveloping boundary layer was discussed. Static wall pressure is measured by a micro-pressure scanning system. The incipient formation positions are pointed out statistically. Without roughness, the longitudinal structures with scales of 1.0h and 1.2h form later and distribute in a longer region compared to that with roughness. Fractal analysis is applied and the averaged fractal dimensions of the overall and sectional flow structures are calculated. If roughness is adopted, the fractal dimension will be larger and the turning point in the sectional fractal dimensions is earlier. However, the dimension tends to be one coincided value in the farther downstream.     

基于流动显示的压缩拐角流动结构定量研究

在Ma = 3.0的超声速风洞中, 采用NPLS技术对上游边界层为层流的25° 压缩拐角进行了流动显示实验, 获得了压缩拐角的精细流动结构, 边界层、剪切层和激波等结构清晰可见. 基于流动显示数据, 采用间歇性、空间相关性和分形分析对流动结构进行了定量研究, 计算了边界层和分离区的间歇因子分布, 获取了边界层中拟序结构和结构角的大小, 给出了边界层分形维数的分布, 并与Ringuette和Bookey等的实验结果进行比较, 阐述了压缩拐角流动结构的定量特征.     

湍流边界层拟序结构的大涡模拟研究

采用动力亚格子模型,利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。从瞬时速度和脉动 速度场、脉动速度相关、均方根脉动涡量分布、以及瞬时涡量场等多个方面,对湍流边界层流动的拟序结构进行了分析, 包括近壁区小尺度湍流结构和瞬态过程,如条纹结构、喷射和扫掠过程、以及近壁旋涡结构等。     

Structure of the channel flow behind a surface-mounted rib under conditions of laminar-turbulent transition

The PIV technique was used to measure the instantaneous vector fields of flow velocity and vorticity behind a thin cross-flow rib installed in a channel with laminar, transient, or turbulent flow. The data were treated statistically to determine the fields of mean longitudinal flow velocity and the correlations of pulsating velocity components 〈u′ν′〉 and 〈u′u′〉. Some features of the flow structure developing under conditions of laminar-turbulent transition behind the rib have been revealed.