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

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

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

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

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

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

表面粗糙度的分形特征及其对微通道内层流流动的影响

基于分形几何学,研究了表面粗糙度的分形特征.采用Weierstrass- Mandelbrot函数对多尺度自仿射的表面粗糙度进行了描述;建立了微通道内层流流动的三维模型并对表面粗糙度的影响进行了数值模拟,分析了雷诺数、相对粗糙度和分形维数对流动阻力特性的影响.研究结果表明,与常规尺度通道不同,粗糙微通道的Poiseuille数不再是常数,而是随雷诺数近似线性增加;相对粗糙度越大,流动产生的回流和分离所导致的流动压降越明显.在相同的相对粗糙度下,粗糙表面的分形维数越大,表面轮廓变化就越频繁,这也将导致流动阻 关键词： 粗糙度 层流阻力系数 微通道 分形     

Effects of upstream roughness and Reynolds number on separated and reattached turbulent flow

This paper presents experimental investigation of upstream roughness and Reynolds number effects on the recirculation region over a smooth forward facing step. The upstream rough wall was produced from 1.5 mm sand grains and the Reynolds number based on step height, Re_h, was varied from 2040 to 9130 for both the upstream smooth and rough walls. For the smooth wall, the reattachment length increased monotonically with Re_h to an asymptotic value of 2.2 step heights for Re_h ≥ 6380. Upstream roughness reduced the reattachment length by 44% because of larger momentum deficit and higher turbulence level in the rough wall boundary layer. The mean velocities and Reynolds stresses were also reduced by roughness. The Reynolds shear stress and production of turbulent kinetic energy showed high negative values at the leading edge of the step indicating counter-gradient diffusion. The implications of these results for standard eddy viscosity models are discussed.     

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

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

The fractal measurement of experimental images of supersonic turbulent mixing layer

Flow visualization of supersonic mixing layer has been studied based on the high spatiotemporal resolution Nano-based Planar Laser Scattering(NPLS) method in SML-1 wind tunnel. The corresponding images distinctly reproduced the flow structure of laminar,transitional and turbulent region,with which the fractal measurement can be implemented. Two methods of measuring fractal dimension were introduced and compared. The fractal dimension of the transitional region and the fully developing turbulence region of supersonic mixing layer were measured based on the box-counting method. In the transitional region,the fractal dimension will increase with turbulent intensity. In the fully developing turbulent region,the fractal dimension will not vary apparently for different flow structures,which em-bodies the self-similarity of supersonic turbulence.     

基于三角波瓣混合器的超声速流场精细结构和掺混特性

在超声速吸气式混合层风洞中,采用基于纳米粒子的平面激光散射(NPLS)技术对平板混合层和三角波瓣混合器诱导的混合层流场精细结构进行了对比实验研究.上下两层来流的实测马赫数分别为1.98和2.84,对流马赫数为0.2.NPLS图像清晰地展示了Kelvin-Helmholtz涡、流向涡、波系结构以及大尺度涡结构的配对合并过程.通过对比分析时间相关的NPLS流场图像,发现了大尺度拟序结构随时间发展演化的非定常特性.基于流动显示结果,采用分形维数和间歇因子指标对流场结构和混合特性进行了定量分析.实验研究表明,三角波瓣混合器诱导的流向涡结构显著提高了上下两层来流的掺混效率,其流动远场的分形维数突破了平板混合层中完全湍流区的分形维数值,达到了1.88,流场结构表现出明显的破碎性,有利于流动在标量层面的扩散和掺混.流动间歇性分析表明,流向涡与展向涡的相互剪切作用主导着混合层的掺混特性,同时由于流向涡的卷吸作用,三角波瓣混合器诱导的混合层混合区域更大,更多的流质被卷入混合区完成混合.     

基于Ramp-VG阵列的超声速混合层流动控制实验研究

利用基于纳米示踪的平面激光散射(nano-tracer-based planar laser scattering,NPLS)技术研究了Ramp-VG阵列对超声速混合层流场的控制效果.对流Mach数Mac=0.17.通过比较无控和控制状态下的混合层NPLS图像,发现控制状态下混合层流动速度提高了5％～15％,K-H不稳...     

A flow control study of a supersonic mixing layer via NPLS

The flow control of a supersonic mixing layer was studied in a supersonic mixing layer wind tunnel with convective Mach number (Mc) at 0.5. The passive control of the mixing layer was achieved by perturbation tapes on the trailing edge of the splitter plate. The control effects of 2D and 3D perturbation tapes with different sizes were compared. The mixing layer was visualized via NPLS, and the transient fine structures were identifiable in NPLS images, which were used to analyze the effects of flow control. The results show that the 2D tapes can enhance the 2D characteristic of the mixing layer, delaying mixing layer transition; and the 3D tapes can enhance the 3D characteristic of the mixing layer, advancing mixing layer transition. 3D structures of the mixing layer were visualized, and the H-type Λ vortexes were found with 3D tapes control.     

基于分形多孔介质三维网络模型的非混溶两相流驱替数值模拟

提出了一个描述多孔介质孔隙尺寸分布的三维分形网络模型,利用该模型对多孔介质中的非混溶两相流驱替进行了数值模拟,研究了孔隙尺寸分布分维D_f和两相流黏滞比M对驱替前沿指进型的影响,结果表明指进型容量维数D_h随着孔隙尺寸分布分维D_f以及黏滞比M的增大而减少,并通过曲线拟合得到了它们之间的定量关系. 关键词： 多孔介质 三维网络 黏滞指进 非混溶两相流     

Turbulent plane wall jets over smooth and rough surfaces

Large-eddy simulations were carried out to study the effects of surface roughness on a plane wall-jet using the Lagrangian dynamic eddy-viscosity subgrid-scale model, at Re = 7500 (based on the jet bulk velocity and height). Results over both smooth and rough surfaces were validated by experimental data at the same Reynolds number. As the jet is injected into the still environment, large-scale rollers are generated in the shear layer between the high-momentum fluid of the jet and the surrounding and are convected downstream with the flow. To understand the extent to which the outer-layer structures modify the flow in the inner layer and the extent to which the effect of roughness spreads away from the wall, both instantaneous and mean flow fields were investigated. The results revealed that, for the Reynolds number and roughness height considered in this study, the effect of roughness is mostly confined to the near-wall region of the wall jet. There is no structural difference between the outer layer of the wall jet over the smooth and rough surfaces. Roughness does not affect the size of the outer-layer structures or the scaling of the profiles of Reynolds stresses in the outer layer. However, in the inner layer, roughness redistributes stresses from streamwise to wall-normal and spanwise directions toward isotropy. Contours of joint probability-density function of the streamwise and wall-normal velocity fluctuations at the bottom of the logarithmic region match those of the turbulent boundary layer at the same height; while the traces of the outer-layer structure were detected at the top of the logarithmic region, indicating that they do not affect the flow very close to the wall, but still modify a major portion of the inner layer. This modification must be taken into consideration when the inner layer of a wall jet is compared with the conventional turbulent boundary layer.     

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.     

高超声速平板边界层流动显示的试验研究

本文在高超声速脉冲式风洞内对基于纳米示踪的平面激光散射技术(nano-based planar laser scattering, NPLS)的应用进行了探索, 并在此基础上对平板边界层流动结构的精细测量进行了研究. 试验来流Ma=7.3, 总压4.8 MPa, 总温680 K. 通过时序的分析和调试, 对各分系统实现了高精度的同步控制; 定量的粒子注入及混合, 实现了粒子的均匀撒播, 对主流获得了均匀的显示效果; 对于边界层流动, 获得了精细的瞬态流动结构图像, 显示了层流到湍流的转捩过程, 并分析了其时空演化特性.     

Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

The technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale l_ν = ν/u_τ (where ν is the kinematic viscosity and u_τ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale l_d, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale l_d is used for rescaling in the inner layer, and the boundary layer thickness δ is used in the outer region. Both parts are then combined for recycling using a blending function. Unlike the blending function proposed by Lund et al. which transitions from the inner layer to the outer layer at approximately 0.2δ, here the location of blending is shifted upwards to enable simulations of very rough surfaces in which the roughness length may exceed the height of 0.2δ assumed in the traditional method. The extended rescaling–recycling method is tested in large eddy simulation of flow over surfaces with various types of roughness element shapes.     

超声速平板圆台突起物绕流实验和数值模拟研究

高速飞行器表面不可避免的存在突起物并形成复杂流场, 从而引起飞行器气动特性和热载荷的变化; 同时, 突起物是流动控制的重要方法之一, 合适的突起物形状及安装位置对于改善冲压发动机进气道性能有重要意义. 本文采用基于纳米粒子的平面激光散射技术(NPLS)研究了马赫3.0来流边界层为层流的平板上三个不同高度圆台突起物绕流流场, 主要关注了突起物后方的尾迹边界层, 并采用高精度的显式五阶精度加权紧致非线性格式(WCNS-E-5)离散求解Navier-Stokes方程模拟了该流场. 获得了超声速圆台绕流精细流场结构, 观察到突起物后方尾迹区域边界层发展的过程. 结合实验和数值模拟结果可以发现, 当圆台高度接近或者小于当地边界层厚度时, 突起物对边界层的扰动非常弱, 圆台后方尾迹边界层能够维持较长距离的层流状态, 在边界层转捩阶段也有清晰的发卡涡结构出现; 反之, 边界层受到的扰动明显增大, 在突起物后方很快发展为湍流; 风洞噪声对本文研究圆台引起的边界层扰动有一定影响, 实验获得的边界层转捩位置要比数值结果靠前. 基于NPLS流场图像, 采用间歇性方法分析了圆台突起物后方边界层的特性, 对于高度大于边界层厚度的圆台其间歇性曲线较为接近并且更加饱满, 边界层的脉动也更为强烈.     

Measurement of Aggregate Fractal Dimensions Using Static Light Scattering

Aggregates formed from colloidal particles will vary in shape according to the aggregation regime prevalent. Compact structures are formed when the aggregation is slow, whilst loose tenuous structures are formed when rapid (or diffusion limited) aggregation prevails. These structures can be fractal in nature, that is, there is a relationship between porosity and the number of primary particles making up the aggregate, and is described by the fractal dimension, d_F. Fractal dimensions of hematite aggregates have been measured experimentally using the static light scattering technique. Fractal dimensions varied with aggregation regimes; for the rapid aggregation regime, d_F was found to be 2.8, whilst for conditions in which aggregation was slow (retardation forces prevail), d_F's of 2.3 were measured. For conditions which lead to aggregation in which both diffusion and retardation forces play a part, structures with fractal dimensions such that 2.3 < d_F < 2.8 were found. The effects of adsorbed fulvic acid, a naturally occuring organic acid, on the kinetics of hematite aggregation and on the resulting structure of hematite aggregates were also investigated. The study of aggregate structure shows that the fractal dimensions of hematite aggregates which are partially coated with fulvic acid molecules are higher than those obtained with no adsorbed fulvic acid. The scattering exponents obtained from static light scattering experiments of these aggregates range from 2.83 ± 0.08 to 3.42 ± 0.1. The scattering exponents of greater than 3 indicate that the scattering is the result of objects that contains pores which are bounded by surfaces with a fractal structure, and can be related only to surface fractal dimension. The high fractal dimensions are due to restructuring within the aggregates, which only occured at low coverage by the organic acid.     

Frequency and spatial characteristics of the electrophysical parameters of a living tree trunk

The method of vertical electrical sounding is applied to derive the spatial characteristic of the dc resistivity ρ of a living pine trunk, ρ = L-0.85, and the frequency characteristic of the ac resistance R of the same pine trunk, R = f-0.053. The results are simulated by the methods of fractal geometry, according to which ρ = L ^{−h + 2} and $ R \sim f^{ - \tfrac{3} {h} + 1} $ R \sim f^{ - \tfrac{3} {h} + 1} , where h is the walk dimension of the electric current in the trunk. From comparison between the experimental and theoretical results, it follows that h = 2.85 and fractal dimension D = 1/h = 0.35. Since D = 1, conducting layers form a Cantor set.     

高超声速圆锥边界层转捩实验研究

在高超声速静音风洞内, 通过基于纳米粒子示踪的平面激光散射(nano-tracer-based planar laser scattering, NPLS)技术、高频压力传感器和温敏漆(temperature sensitive paints, TSP)技术开展了0°攻角条件下7°直圆锥高超声速边界层转捩相关实验研究, 得到了圆锥边界层由层流发展至湍流完整过程的NPLS图像, 清晰地展示了第2模态波的"绳状"结构, 尖锥与钝锥边界层的NPLS结果表明尖锥边界层转捩中第2模态波占主导, 而钝锥边界层在转捩前出现波长约为第2模态波波长5倍(甚至更长)、特征频率不高于31 kHz的狭长涡结构; 采用功率谱密度(power spectrum density, PSD)分析、互相关和N值计算对高频脉动压力数据进行分析, 得到了边界层内扰动波的发展规律, 在尖锥和钝锥中均观察到了沿流向第2模态波幅值先增大后减小、特征频率逐渐降低, 低频成分逐渐增加, 表明边界层发展过程中第2模态率先发展达到饱和, 而后逐渐衰减, 而低频模态则逐渐发展; 通过TSP技术得到了不同单位Reynolds数下的圆锥表面温升分布, 结果表明, 随单位Reynolds数增大, 边界层转捩阵面前移.      

Study on the dynamic characteristics of flow over building cluster at high Reynolds number by large eddy simulation

In this paper,the dynamic characteristics of building clusters are simulated by large eddy simulation at high Reynolds number for both homogeneous and heterogeneous building clusters.To save the computational cost a channel-like flow model is applied to the urban canopy with free slip condition at the upper boundary.The results show that the domain height is an important parameter for correct evaluation of the dynamic characteristics.The domain height must be greater than 8h(h is the average building height)in order to obtain correct roughness height while displacement height and roughness sublayer are less sensitive to the domain height.The Reynolds number effects on the dynamic characteristics and flow patterns are investigated.The turbulence intensity is stronger inside building cluster at high Reynolds number while turbulence intensity is almost unchanged with Reynolds number above the building cluster.Roughness height increases monotonously with Reynolds number by 20%from Re*=103 to Re*=105 but displacement height is almost unchanged.Within the canopy layer of heterogeneous building clusters,flow structures vary between buildings and turbulence is more active at high Reynolds number.