Spatio-temporal invariants of the time reversal operator

The decomposition of the time reversal operator, known by the French acronym DORT, is a technique to extract point scatterers' monochromatic Green's functions from a medium. It is used to detect, locate, and focus on scatterers in various domains such as underwater acoustics, medical ultrasound, and nondestructive evaluation. A limitation of the method arises from its single-frequency nature, when the signals used in acoustics are often broadband. Reconstruction of the broadband Green's functions from the single-frequency Green's functions can be very difficult when numerous scatterers are present in the medium. Moreover, the method does not take advantage of the axial resolution associated with broadband signals. Time domain methods are investigated here as an answer to these problems. It is shown that the time reversal operator in the time domain takes the form of a tensor. The properties of the invariants are discussed. It is shown they do not have all the expected properties. Another method is proposed that requires a priori information on the medium.     

Revisiting iterative time reversal processing: application to detection of multiple targets

The iterative time reversal processing represents a high speed and easy way to self-focus on the strongest scatterer in a multitarget medium. However, finding weaker scatterers is a more difficult task that can be solved by computing the eigenvalue and eigenvector decomposition of the time reversal operator, the so-called DORT method. Nevertheless, as it requires the measurement of the complete interelements response matrix and time-consuming computation, the separation of multiple targets may not be achieved in real time. In this study, a new real time technique is proposed for multitarget selective focusing that does not require the experimental acquisition of the time reversal operator. This technique achieves the operator decomposition using a particular sequence of filtered waves propagation instead of computational power. Due to its simplicity of implementation, this iterative process can be achieved in real time. This high speed selective focusing is experimentally demonstrated by detecting targets through a heterogeneous medium and in a speckle environment. A theoretical analysis compares this technique to the DORT formalism.     

Time reversal of speckle noise

Focusing a wave in an unknown inhomogeneous medium is an open problem in wave physics. This work presents an iterative method able to focus in pulse-echo mode in an inhomogeneous medium containing a random distribution of scatterers. By performing a coherent summation of the random echoes backscattered from a set of points surrounding the desired focus, a virtual bright pointlike reflector is generated. A time-reversal method enables an iterative convergence towards the optimal wave field focusing at the location of this virtual scatterer. Thanks to this iterative time-reversal process, it is possible to focus at any arbitrary point in the heterogeneous medium even in the absence of pointlike source. An experimental demonstration is given for the correction of strongly distorted images in the field of medical ultrasound imaging. This concept enables envisioning many other applications in wave physics.     

Optimal adaptive focusing through heterogeneous media with the minimally invasive inverse filter

The inverse filter is a technique used to adaptively focus waves through heterogeneous media. It is based on the inversion of the Green's functions matrix between the M transducers of a focusing array and N control points in the focal area. The inverse filter minimizes the pressure deposited around the focal point. However it is highly invasive, requiring the presence of N transducers or hydrophones in the focal area at the control points' locations to measure the Green's functions. This paper presents a way of reaching the inverse filter's focusing quality with a minimally invasive setup: only one transducer (at the desired focal location) is needed. This minimally invasive inverse filter takes advantage of the fact all the information about the propagation medium can be retrieved from the signals backscattered by the medium towards the focusing array, if the propagation medium is lossless. A numerical simulation is performed to test this minimally invasive inverse filter through a scattering, lossless medium. The focusing quality equals the conventional, highly invasive inverse filter's. The average spatial and temporal contrast is increased by up to 10 dB compared to the time reversal focusing.     

Experimental detection and focusing in shallow water by decomposition of the time reversal operator

A rigid 24-element source-receiver array in the 10-15 kHz frequency band, connected to a programmable electronic system, was deployed in the Bay of Brest during spring 2005. In this 10- to 18-m-deep environment, backscattered data from submerged targets were recorded. Successful detection and focusing experiments in very shallow water using the decomposition of the time reversal operator (DORT method) are shown. The ability of the DORT method to separate the echo of a target from reverberation as well as the echo from two different targets at 250 m is shown. An example of active focusing within the waveguide using the first invariant of the time reversal operator is presented, showing the enhanced focusing capability. Furthermore, the localization of the scatterers in the water column is obtained using a range-dependent acoustic model.     

Ultrasound speckle and equivalent scatterers

B-mode ultrasound images are characterised by the speckle artefact, which introduces fine-false structures whose apparent resolution is beyond the imaging system capabilities. Speckle presence is due to interference effects between overlapping echoes and its occurrence is related to a great number of randomly distributed structure scatterers within a resolution cell. Basing our analysis on linear system theory, we show that a dense random set of scatterers can be substituted by an equivalent one with a much smaller number of periodic scatterers. This new structure with regularly distributed scatterers is able to give rise to the same B-mode image and the same speckle pattern, for a given ultrasound pulse. This new approach helps the understanding of the deterministic nature of speckle and may reduce drastically the computing time in numerical simulations. Additionally, it can contribute to periodicity analysis used in tissue characterisation.     

Real time inverse filter focusing through iterative time reversal

In order to achieve an optimal focusing through heterogeneous media we need to build the inverse filter of the propagation operator. Time reversal is an easy and robust way to achieve such an inverse filter in nondissipative media. However, as soon as losses appear in the medium, time reversal is not equivalent to the inverse filter anymore. Consequently, it does not produce the optimal focusing and beam degradations may appear. In such cases, we showed in previous works that the optimal focusing can be recovered by using the so-called spatiotemporal inverse filter technique. This process requires the presence of a complete set of receivers inside the medium. It allows one to reach the optimal focusing even in extreme situations such as ultrasonic focusing through human skull or audible sound focusing in strongly reverberant rooms. But, this technique is time consuming and implied fastidious numerical calculations. In this paper we propose a new way to process this inverse filter focusing technique in real time and without any calculation. The new process is based on iterative time reversal process. Contrary to the classical inverse filter technique, this iteration does not require any computation and achieves the inverse filter in an experimental way using wave propagation instead of computational power. The convergence from time reversal to inverse filter during the iterative process is theoretically explained. Finally, the feasibility of this iterative technique is experimentally demonstrated for ultrasound applications.     

Detection of cracks in a thin air-filled hollow cylinder by application of the DORT method to elastic components of the echo

The DORT method (Decomposition de l'Opérateur de Retournement Temporel in French) is a scattering analysis technique which uses arrays of transducers. This method is efficient for detection of selective focusing on point-like scatterers. It has been also applied to analyze the scattering by an air-filled cylindrical steel shell immersed in water. It was shown that the diagonalization of the time reversal operator allows us to separate the different elastic components of the scattered field. Here, we apply the method to detect flaws in hollow cylinders. In this case, the dominant components are the three circumferential waves (A0, A1 and S0 Lamb modes). Each Lamb mode corresponds to an invariant of the time reversal operator. The dispersion curves of these waves are calculated from these invariants. Resonance frequencies of the shell are deduced from the frequency dependence of the eigenvalues of the time reversal operator. It is shown that the presence of a crack (0.2 mm in depth) affects significantly the eigenvalue distribution of the time reversal operator. Thus, the DORT method offers a new means for detecting defects in a shell.     

Time reversal and the inverse filter

To focus ultrasonic waves in an unknown inhomogeneous medium using a phased array, one has to calculate the optimal set of signals to be applied on the transducers of the array. In the case of time-reversal mirrors, one assumes that a source is available at the focus, providing the Green's function of this point. In this paper, the robustness of this time-reversal method is investigated when loss of information breaks the time-reversal invariance. It arises in dissipative media or when the field radiated by the source is not entirely measured by the limited aperture of a time-reversal mirror. However, in both cases, linearity and reciprocity relations ensure time reversal to achieve a spatiotemporal matched filtering. Nevertheless, though it provides robustness to this method, no constraints are imposed on the field out of the focus and sidelobes may appear. Another approach consists of measuring the Green's functions associated to the focus but also to neighboring points. Thus, the whole information characterizing the medium is known and the inverse source problem can be solved. A matrix formalism of the propagation operator is introduced to compare the time-reversal and inverse filter techniques. Moreover, experiments investigated in various media are presented to illustrate this comparison.     

Time reversed reverberation focusing in a waveguide

Time reversal mirrors have been applied to focus energy at probe source locations and point scatterers in inhomogeneous media. In this paper, we investigate the application of a time reversal mirror to rough interface reverberation processing in a waveguide. The method is based on the decomposition of the time reversal operator which is computed from the transfer matrix measured on a source-receiver array [Prada et al., J. Acoust. Soc. Am. 99, 2067-2076 (1996)]. In a similar manner, reverberation data collected on a source-receiver array can be filtered through an appropriate temporal window to form a time reversal operator. The most energetic eigenvector of the time reversal operator focuses along the interface at the range corresponding to the filter delay. It is also shown that improved signal-to-noise ratio measurement of the time reversal operator can be obtained by ensonifying the water column with a set of orthogonal array beams. Since these methods do not depend upon a priori environmental information, they are applicable to complex shallow water environments. Numerical simulations with a Pekeris waveguide demonstrate this method.     

亚波长间距理想导体球阵列近区时间反演电磁场的快速求解

为快速求解亚波长间距分布的理想导体球阵列近区的时间反演电磁场,提出一种基于等效偶极子模型的解析分析方法.首先,通过分析球面波照射理想导体小球的散射场解析解发现,散射场可以近似等效为电磁偶极子辐射场的叠加.等效偶极子的强度与初始激励源的幅度成正比关系.通过建立不同小球等效偶极子矢量间的耦合方程组可以直接求解得到相应矢量的大小.然后,结合时间反演腔理论得到相应的时间反演并矢格林函数,继而得到小球阵列近区的时间反演场分布.最后,通过与数值仿真软件的计算结果进行对比,验证了方法的正确性及高效性.研究表明,时间反演技术结合近场亚波长间距小散射体加载能够实现超分辨率的场聚焦.     

Specific features of the manifestation of the Doppler effect in diffraction of focused coherent beams in a scattering flow

Intensity fluctuations of a dynamic speckle field formed by diffraction of sharply focused coherent beams in a random flow are studied. Dependences of the form of the spectrum of scattered light intensity fluctuations on the speckle observation angle, flow velocity, and concentration of scatterers are found. The influence of the unscattered component on the formation of the statistically inhomogeneous dynamic speckle field is determined.     

超声多波聚焦及声偏振方向控制方法*

在超声多波聚焦思想的基础上,通过数值模拟计算,分析了时间反转法的多波聚焦特性以及对声场偏振方向进行控制的可行性。结果表明,在待测目标的不同位置处,时间反转法都能够实现多波聚焦的效果,使具有不同传播速度、不同偏振特性的多种声波自适应聚焦。但是,在介质的近表面处,由于受到表面波的影响,多波聚焦声场仍然具有椭圆偏振特性,无法实现声场偏振方向精确控制的目的;而位于介质内部的多波聚焦点受到表面波影响很小,数值计算结果表明此时多波聚焦声场具有线偏振特性,通过改变声源前后两个脉冲的激发幅度和相位,可以控制声场的偏振方向,达到偏振方向扫描的目的。该文的研究为精确检测裂纹方向或界面性质提供一种可能的途径。     

Imaging in the presence of grain noise using the decomposition of the time reversal operator

In this paper, we are interested in detecting and imaging defects in samples of cylindrical geometry with large speckle noise due to the microstructure. The time reversal process is an appropriate technique for detecting flaws in such heterogeneous media as titanium billets. Furthermore, time reversal can be iterated to select the defect with the strongest reflectivity and to reduce the contribution of speckle noise. The DORT (the French acronym for Decomposition of the Time Reversal Operator) method derives from the mathematical analysis of the time reversal process. This detection technique allows the determination of a set of signals to be applied to the transducers in order to focus on each defect separately. In this paper, we compare three immersion techniques on a titanium sample, standard transmit/receive focusing, the time reversal mirror (TRM), and the DORT method. We compare the sensitivity of these three techniques, especially the sensitivity to a poor alignment of the array with the front face of the sample. Then we show how images of the sample can be obtained with the TRM and the DORT method using backpropagation algorithm.     

时间反转方法用于高强度聚焦超声的焦点偏移补偿

传统的高强度聚焦超声(HIFU)治疗中实际焦点和预设焦点容易出现偏移,为考察时间反转方法对HIFU治疗中焦点偏移的补偿效果,采用时域有限差分方法求解Westervelt方程,建立高强度聚焦声场数值模型。数值计算得到在人体软组织中进行HIFU治疗时,采用时间反转方法后焦点偏移距离最大仅为1.6 mm。脂肪层厚度及声源强度改变对时间反转聚焦精度影响不大,F数(焦点距离同换能器孔径的比值)降低时,焦点偏移减小。研究表明在人体软组织吸收系数和非线性系数范围内,时间反转方法可有效补偿焦点偏移,达到更好的聚焦效果。      

Time reversal of flexural waves in a beam at audible frequency

There has been very limited work on the application of time reversal to the propagation of audible frequency waves in mechanical structures. The present work concentrates on the application of time reversal to the focusing of audible range, flexural waves in an infinite beam, and to the detection of local heterogeneity in such a beam. Practical applications of time reversal of flexural waves in structures include vibration energy focusing, detection of vibratory or acoustic sources, and detection of defects in mechanical structures. An analytical model of flexural wave propagation in the beam as well as sensing and emission using piezoelectric transducers is presented. Time reversal experiments are conducted and compared to the model results in either a homogeneous beam or a beam with point mass heterogeneities. In the various situations tested, it is shown that time reversal effectively compensates the spreading in time of the impulse due to the dispersive propagation of flexural waves. One interesting aspect of this property is the generation of large amplitude impulsive responses in the beam using remote actuators. Finally, the "Decomposition de l'Operateur de Retournement Temporel" approach is examined to detect and localize point mass scatterers in the beam.     

Optimal focusing by spatio-temporal inverse filter. II. Experiments. Application to focusing through absorbing and reverberating media

To focus ultrasonic waves in an unknown heterogeneous medium using a phased array, one has to calculate the optimal set of signals to be applied on the transducers of the array. (In most applications of ultrasound, medical imaging, medical therapy, nondestructive testing, the first step consists of focusing a broadband ultrasound beam deeply inside the medium to be investigated.) Focusing in a homogeneous medium simply requires to compensate for the varying focus-array elements geometrical distances. Nevertheless, heterogeneities in the medium, in terms of speed of sound, density, or absorption, may strongly degrade the focusing. Different techniques have been developed in order to correct such aberrations induced by heterogeneous media (time reversal, speckle brightness, for example). In the companion to this paper, a new broadband focusing technique was investigated: the spatio-temporal inverse filter. Experimental results obtained in various media, such as reverberating and absorbing media, are presented here. In particular, intraplate echoes suppression and high-quality focusing through a human skull, as well as hyper-resolution in a reverberating medium, will be shown. It is important to notice that all these experiments were performed with fully programmable multichannel electronics whose use is required to fully exploit the spatio-temporal technique.     

Fourier magnitude of the field incident on a random scattering medium from spatial speckle intensity correlations

Spatial speckle intensity correlations are used to determine the spatial Fourier magnitude of a field incident on a random scattering medium. The patterned beam is scanned across the scattering medium, and the speckle pattern on the opposite side is imaged at each beam position. A theory based on a Green's function representation is used to reconstruct the spatial Fourier magnitude of the patterned incident field.