Velocity and density field measurements by digital speckle method

Velocity and density field measurements based on image processing of laser speckle or pseudo-speckle pattern have been developed. Laser speckle velocimetry (LSV) or white-light speckle velocimetry (WSV), which corresponds to a high-image-density PIV, gives a local velocity vector map or whole field velocity contour map of a two-dimensional flow field seeded densely with fine particles. This technique has an advantage in high-speed flow measurement without limit of frame rate except for directional ambiguity. New techniques of laser speckle photography and laser speckle interferometry by means of digital image processing have been developed recently for density field measurements. In laser speckle photography, a local density gradient vector map is reconstructed by cross-correlation evaluation between the reference and the object speckle patterns. In laser speckle interferometry, an equi-density contour map is reconstructed by image subtraction between the reference and the object interferometric speckle patterns.     

电子散斑归一化互相关法测位移

把粗糙物体表面受到相干光照明时产生的散斑场看作是一个稳态随机信号,在此基础上,本文提出了一种能做到大量程,高精度位移测量的新方法——电子散斑归一化互相关法测位移或称电子散斑照相法测位移.该法实质上是借助数字图像处理技术,通过物体位移前后二个散斑场之间的相关性来进行位移测量方法的发展.文中给出了理论分析和实验系统.实验结果表明,该法有诸如准实时,大量程(毫米级),高精度和线性可靠等优点.     

气相爆轰波在分叉管中传播现象的数值研究

数值研究气相爆轰波在分叉管中的传播现象.用二阶附加半隐龙格-库塔法和5阶WENO格式求解二维欧拉方程,用基元反应描述爆轰化学反应过程,得到了密度、压力、温度、典型组元质量分数场及数值胞格结构和爆轰波平均速度.结果表明:气相爆轰波在分叉管中传播,分叉口左尖点的稀疏波导致诱导激波后压力、温度急剧下降,诱导激波和化学反应区分离,爆轰波衰减为爆燃波(即爆轰熄灭).分离后的诱导激波在垂直支管右壁面反射,并导致二次起爆.畸变的诱导激波在水平和垂直支管中均发生马赫反射.分叉口上游均匀胞格区和分叉口附近大胞格区的边界不是直线,其起点通常位于分叉口左尖点上游或恰在左尖点.水平支管中马赫反射三波点迹线始于右尖点下游.分叉口左尖点附近的流场中出现了复杂的旋涡结构、未反应区及激波与旋涡作用.旋涡加速了未反应区的化学反应速率.反射激波与旋涡作用并使旋涡破碎.反射激波与未反应区作用,加速其反应消耗,并形成一个内嵌的射流.数值计算得到的波系演变和胞格结构与实验定性一致.     

Whole field spatial frequency analysis of double-or multiple-exposed PIV images

Speckle Velocimetry or high-image density PIV gives a velocity vector map of a twodimensional flow field by point-by-point spatial frequency analysis of local pattern at an interrogation spot in a double- or a multiple-exposed image of laser speckle or pseudo-speckle pattern generated by pulsed laser-light-sheet illuminations of a plane in the flow densely seeded with fine particles. The whole field spatial frequency analysis of the double- or multiple-exposed PIV image gives more useful information of the flow field. Optical Fourier transform is a conventional technique not only for the local spatial frequency analysis but for the whole field analysis. Filtering of spatial frequency is one of the typical techniques for the latter which can reconstruct a velocity contour or a component velocity contour map of the flow. Fundamentally this technique is a simple and efficient analogue method to get more information in the velocity field analysis of the flow than a digital image processing, but in practical applications the measurement is restricted to simple flows because of optical noise and low spatial resolution. In order to improve the technique, the fundamental characteristics of the filtering and the noise yielded in the filtering process must be investigated. Meanwhile, wavelet analysis can also be applied to the whole field spatial frequency analysis of PIV image. In this paper the filtering technique is examined by numerical convolution integral, and the results obtained are compared with ones obtained by the wavelet analysis.     

Modulated speckle simulations based on the random-walk model

The random-walk model is employed to simulate modulated speckle patterns. We demonstrate that the geometrical image approximation fails to describe the modulated speckle pattern. A new approach to analyzing this phenomenon is proposed. The validity of the approximations employed is verified by comparison of the simulation with the experimental results. Speckle metrological applications and phase measurement techniques could be improved by taking advantage of this model.     

强梯度复杂流场中的粒子动力学响应试验研究

示踪粒子在(高)超声速流场中的动力学响应是粒子成像测速等粒子示踪测量技术的关键问题之一.现有文献对粒子动力学响应的试验测量往往是通过单个斜激波响应的测量方法.然而,当示踪粒子用于测量高速飞行器发动机内部复杂的激波串流场时,粒子将经历由多道激波导致的速度、压力、黏性等剧烈变化.本文结合目前(高)超声速飞行器的研究热潮,重点关注示踪粒子在应用于发动机内部具有连续激波的复杂流场测量中存在的跟随性评估方面,开展了一系列的相关试验研究.包括测量超声速风洞的喷管出口速度分布以验证测试系统的性能,在马赫4.2和3.0流场中测量了粒子对二维10°和15°单斜劈绕流中的斜激波动力响应,并测量了模拟发动机内部连续梯度的双斜劈粒子斜激波动力响应.结合粒子动力学的理论模型,得到了各状态的粒子弛豫时间、弛豫距离、Stokes数.基于图像方法、统计学规律分析了激波非定常抖动对测量结果的影响,并对测量结果进行了修正.结果显示,相同斜劈角度下,马赫数越高,粒子的弛豫时间、弛豫距离就越大.但是在相同的来流马赫数下,斜劈角度越大,粒子的弛豫时间、弛豫距离反而减小.在强梯度之后由于流场的雷诺数和黏性系数变化剧烈,粒子的跟随性降低了大约5.7%,stokes数增加了约1%.虽然在本文条件下Stokes数仍满足超声速流场对粒子跟随性的要求,但粒子响应的降低无疑是值得关注的,尤其是当其被应用于具有更多连续梯度的复杂流场测量中.     

Numerical study of shock wave interaction in steady flows of a viscous heat-conducting gas with a low ratio of specific heats

Specific features of shock wave interaction in a viscous heat-conducting gas with a low ratio of specific heats are numerically studied. The case of the Mach reflection of shock waves with a negative angle of the reflected wave with respect to the free-stream velocity vector is considered, and the influence of viscosity on the flow structure is analyzed. Various issues of nonuniqueness of the shock wave configuration for different Reynolds numbers are discussed. Depending on the initial conditions and Reynolds numbers, two different shock wave configurations may exist: regular configuration interacting with an expansion fan and Mach configuration. In the dual solution domain, a possibility of the transition from regular to the Mach reflection of shock waves is considered.     

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

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

Improvement of ultrasound speckle image velocimetry using image enhancement techniques

Ultrasound-based techniques have been developed and widely used in noninvasive measurement of blood velocity. Speckle image velocimetry (SIV), which applies a cross-correlation algorithm to consecutive B-mode images of blood flow has often been employed owing to its better spatial resolution compared with conventional Doppler-based measurement techniques. The SIV technique utilizes speckles backscattered from red blood cell (RBC) aggregates as flow tracers. Hence, the intensity and size of such speckles are highly dependent on hemodynamic conditions. The grayscale intensity of speckle images varies along the radial direction of blood vessels because of the shear rate dependence of RBC aggregation. This inhomogeneous distribution of echo speckles decreases the signal-to-noise ratio (SNR) of a cross-correlation analysis and produces spurious results. In the present study, image-enhancement techniques such as contrast-limited adaptive histogram equalization (CLAHE), min/max technique, and subtraction of background image (SB) method were applied to speckle images to achieve a more accurate SIV measurement. A mechanical sector ultrasound scanner was used to obtain ultrasound speckle images from rat blood under steady and pulsatile flows. The effects of the image-enhancement techniques on SIV analysis were evaluated by comparing image intensities, velocities, and cross-correlation maps. The velocity profiles and wall shear rate (WSR) obtained from RBC suspension images were compared with the analytical solution for validation. In addition, the image-enhancement techniques were applied to in vivo measurement of blood flow in human vein. The experimental results of both in vitro and in vivo SIV measurements show that the intensity gradient in heterogeneous speckles has substantial influence on the cross-correlation analysis. The image-enhancement techniques used in this study can minimize errors encountered in ultrasound SIV measurement in which RBCs are used as flow tracers instead of exogenous contrast agents.     

Application of improved Speckle contouring technique to surface roughness measurements

One of the well known drawbacks in Electronic Speckle Pattern Interferometry (ESPI) is the poor fringe density due to electronic and speckle noise and decreasing fringe visibility with increasing speckle decorrelation. As recently reported we suggested improving the fringe density and enlarging the measurement range by the incremental addition of phase images. The technique has already been successfully applied to in-plane and out-of-plane displacement field measurements of carbon-fibre reinforced composites. This paper demonstrates how the technique can be used to improve the sensitivity and to extend the measurement range for speckle contouring applications. The angular correlation of the speckle phase images is used to get the desired 3D information of the investigated object surface. It is demonstrated how the angular speckle contouring techinque can easily be automated with respect, for example, to possible industrial applications. Roughness measurements are carried out on milled Rugo test surfaces. All measurements are performed without the usual vibration isolation, and of course, without the treatment of the surface with fine, white powder for contrast enhancement, which would falsify the results.     

凝聚炸药爆轰波在高声速材料界面上的折射现象分析

凝聚炸药爆轰在边界高声速材料约束下传播时,爆轰波会在约束材料界面上产生复杂的折射现象.本文针对凝聚炸药爆轰波在高声速材料界面上的折射现象展开理论和数值模拟分析.首先通过建立在爆轰ZND模型上的改进爆轰波极曲线理论给出爆轰波折射类型,然后发展一种求解爆轰反应流动方程的基于特征理论的二阶单元中心型Lagrange计算方法来数值模拟典型的爆轰波折射过程.从改进爆轰波极曲线理论和二阶Lagrange方法数值模拟给出的结果看出,凝聚炸药爆轰波在高声速材料界面上的折射类型有四种:反射冲击波的正规折射、带束缚前驱波的非正规折射、带双Mach反射的非正规折射、带λ波结构的非正规折射.     

Temporal heterodyne shearing speckle pattern interferometry

Shearing speckle pattern interferometry is a full-field speckle interferometric technique used to determine surface displacement derivatives. In this paper, a new measurement system of real-time heterodyne shearography interferometry is presented. This system combined with heterodyne measurement, shearography interferometry and time domain signal processing technology can dynamically detect the out-of-plane displacement gradient. The principles and system arrangement are described. Using the Jones matrix, the mathematical expression of light intensity distribution passing through this system is deduced. A preliminary experiment was performed to demonstrate the performance of this new device, and simulations were conducted using the finite element method. Comparison of results shows that quantitative measurement of the displacement derivative has been achieved.     

激波聚焦反射的实验和数值研究

 采用高速摄影技术和数值模拟方法，对入射激波在两种不同形状的抛物形反射器表面聚焦和反射的过程进行了研究，得到了激波聚焦反射过程的波系结构的实验阴影照片和数值计算结果，两者符合得很好。对激波聚焦形成的气体动力学焦点的特性进行了分析，结果表明，入射激波在两种反射器反射后聚焦所形成的气体动力学焦点均是由三波点在轴心处的会聚所导致的，气体动力学焦点位于相应的反射器壁面的几何焦点附近。不同的反射器中，激波聚焦前后的波系结构也不同，对较浅的反射器，入射激波反射前在反射器壁面形成了弓形激波，反射之后需要相对较长的时间完成聚焦，形成聚焦反射激波之后弓形激波仍未相交；对较深的反射器，入射激波反射后在更短的时间内聚焦，聚焦时弓形激波已经相交，聚焦反射激波之后的流场波系结构更加复杂。     

Generalized hydrodynamic theory of shock waves: Mach-number dependence of inverse shock width for nitrogen gas

A generalized hydrodynamic theory of shock wave structures is presented for rigid diatomic gases. The theory yields shock solutions for all Mach numbers. The accuracy of the theory is examined for nitrogen gas as an example. Excellent agreement with experiment is obtained for the inverse shock widths for the entire range of Mach number studied by experiment. The comparison made shows that we are now in possession of a successful continuum theory of shock waves in molecular gases for all Mach numbers.     

Flow visualization by means of single-exposure speckle photography

A photograph of a flow field is taken under laser illumination. If the exposure is short enough, the velocity distribution in the field will be mapped on the photograph as variations in speckle contrast. These contrast variations can be converted to intensity variations by means of a simple spatial filtering technique, to give a direct picture of the velocity distribution in the flow field. A potential application of this technique to the mapping of retinal blood flow is described.     

EPSI validation for cylindrical gears

Shape deviations of spur gears obtained by using Electronic Speckle Pattern Interferometry are compared with the results of numerical simulation. The dual illumination method produces a speckle pattern by simultaneous illumination of a sample by two laser waves symmetrically in the direction of observation. Two images obtained, respectively, before and after deformation are subtracted to obtain the field of displacement. Furthermore, a study of meshing stiffness was performed for different types of setting for an internal gear.     

Influence of the adiabatic index on switching between different types of shock wave reflection in a steady supersonic gas flow

A comprehensive pattern of different types of shock wave reflection in a steady supersonic gas flow is analytically constructed with regard to a new wave configuration found by the authors-negative-angle irregular reflection. A double Mach reflection with a negative reflection angle in a steady supersonic gas flow is numerically obtained for the first time.     

Detonation re-initiation mechanism following the Mach reflection of a quenched detonation

This experimental study addresses the re-initiation mechanism of detonation waves following the Mach reflection of a shock–flame complex. The detonation diffraction around a cylinder is used to reproducibly generate the shock–flame complex of interest. The experiments are performed in methane–oxygen. We use a novel experimental technique of coupling a two-in-line-spark flash system with a double-frame camera in order to obtain microsecond time resolution permitting accurate schlieren velocimetry. The first series of experiments compares the non-reactive sequence of shock reflections with the reflection over a rough wall under identical conditions. It was found that the hot reaction products generated along the rough wall are entrained by the wall jet into a large vortex structure behind the Mach stem. The second series of experiments performed in more sensitive mixtures addressed the sequence of events leading to the detonation establishment along the Mach and transverse waves. Following ignition and jet entrainment, a detonation first appears along the Mach stem while the transverse wave remains non-reactive. The structure of the unburned tongue however indicates local instabilities and hot spot formation, leading to the rapid reaction of this gas. Numerical simulations are also reported, confirming the sequence of ignition events obtained experimentally.     

Stationary collisionless shock waves in an initial plasma with high ion temperature

Stationary collisonless shock waves propagating perpendicularly to an initial magnetic field are produced by the fast-rising magnetic field \((\dot B = 7 \cdot 10^{10} G/sec)\) of a theta pinch (coil diameter 16 cm, coil length 60 cm). The initial plasma is produced by a fast theta pinch discharge (810 kHz). At filling pressures between 5 and 15 mtorr H₂ or D₂ the degree of ionization is about 50%. By choosing the filling pressure properly it is possible to trap a homogeneous magnetic field. The ions of this plasma have a temperature of a few 10 eV. This value is much higher than the electron temperature and results in a local plasmaβ between 0.3 and 5. In this initial plasma stationary collisionless shock waves with Mach numbers between 1.5 and 5 are observed. The snow-plough model is used to derive conditions for the stationary state, attainable Mach number, and velocity of the front which relate the external magnetic field and the parameters of the initial plasma. Strong collisionless dissipation can be demonstrated by measuring the profiles of magnetic field, density, and electron temperature of the shock waves. For the electrons this dissipation mechanism can be described by an effective collison frequency. This phenomenologically introduced frequency determines the width of the shock front at least for subcritical shock waves. It exceeds the classical electron-ion collision frequency by 1–2 orders of magnitude and is roughly equal to one-third of the ion plasma frequency. The ion temperature can be estimated from the steady state conservation relations. The ions are heated in the two degrees of freedom perpendicular to the magnetic field. For shock waves with Mach numbers below the critical one the ions seem to be heated merely adiabatically. In strong shock waves this heating is considerably exceeded, and for high Mach numbers it yields ion temperatures up to about 500 eV. Finally, semi-empirical formulas are derived to estimate the possible temperatures of electrons and ions behind the shock front.     

Detached shock waves around cylinders flying at mach number ranging from 1 to 2

Using a ballistic range with a small single-stage powder gun, characteristics of detached shock waves around cylindrical projectiles were investigated. The detaches shock waves were analyzed quantitatively by visualization method. The projectiles were made of steel, with a diameter of 5.56 mm and with a mass of 1.77 g. The projectiles were ejected into the test section of the ballistic range at supersonic speeds. Direction-indicating color schlieren method (DInCS) was used to visualize the flow field. The stand-off distance of the detached shock waves and the shape were obtained as the function of Mach number. The functions were compared with the results obtained by wind tunnel experiments. Analyzed images provided the quantitative information on shock shapes. The results were as follows: (1) the stand-off distances of the detached shock waves have been shown as exponential functions of Mach number; (2) the shapes of the shock waves have been found as quadratic functions of Mach number. The differences in the functions between this paper and a wind tunnel test were shown.