Measurement of the spectral directivity of optoacoustic and ultrasonic transducers with a laser ultrasonic source

The effect of copper grain size on laser ultrasonic backscattered signal

Grain size has an essential influence on the serviceability of metallic materials. In this paper, a noncontact laser ultrasonic testing platform is built to study the effect of copper grain size on the laser ultrasonic backscattered signal. According to the correlation between grain size and ultrasonic wavelength, the ultrasonic scattering by copper grains in the experiment contains not only Rayleigh scattering but also the transition region from Rayleigh scattering to stochastic scattering. Using time-frequency analysis, the influence of copper grain size on the characteristic parameters of backscattering was explored, and a prediction model of grain size was established, which was compared with the prediction model based on the attenuation method to verify the accuracy of the backscattering model. The results show that the backscattered signal can adequately characterize the grain size information and laser ultrasonics is a method that can realize on-line detection of grain size.

Directivity and Frequency-Dependent Effective Sensitive Element Size of Membrane Hydrophones: Theory Versus Experiment

It is important to know hydrophone frequency-dependent effective sensitive element size in order to account for spatial averaging artifacts in acoustic output measurements. Frequency-dependent effective sensitive element size may be obtained from hydrophone directivity measurements. Directivity was measured at 1, 2, 3, 4, 6, 8, and 10 MHz from ±60° in two orthogonal planes for eight membrane hydrophones with nominal geometrical sensitive element radii ( a_g ) ranging from 100 to [Formula: see text]. The mean precision of directivity measurements (obtained from four repeated measurements at each frequency and angle) averaged over all frequencies, angles, and hydrophones was 5.8%. Frequency-dependent effective hydrophone sensitive element radii a_eff(f) were estimated by fitting the theoretical directional response for a disk receiver to directivity measurements using the sensitive element radius ( a ) as an adjustable parameter. For the eight hydrophones in aggregate, the relative difference between effective and geometrical sensitive element radii, ( a_eff - a_g)/a_g , was fit to C /( ka_g)ⁿ , where k = 2π/λ and λ = wavelength. The functional fit yielded C = 1.89 and n = 1.36 . The root mean square difference between the data and the model was 34%. It was shown that for a given value for a_g , [Formula: see text] for membrane hydrophones far exceeds that for needle hydrophones at low frequencies (e.g., < 4 MHz when [Formula: see text]). This empirical model for [Formula: see text] provides information required for the compensation of spatial averaging artifacts in acoustic output measurements and is useful for choosing an appropriate sensitive element size for a given experiment.

Considerations for Choosing Sensitive Element Size for Needle and Fiber-Optic Hydrophones-Part I: Spatiotemporal Transfer Function and Graphical Guide

The spatiotemporal transfer function for a needle or reflectance-based fiber-optic hydrophone is modeled as separable into the product of two filters corresponding to frequency-dependent sensitivity and spatial averaging. The separable hydrophone transfer function model is verified numerically by comparison to a more general rigid piston spatiotemporal response model that does not assume separability. Spatial averaging effects are characterized by frequency-dependent "effective" sensitive element diameter, which can be more than double the geometrical sensitive element diameter. The transfer function is tested in simulation using a nonlinear focused pressure wave model based on Gaussian harmonic radial pressure distributions. The pressure wave model is validated by comparing to experimental hydrophone scans of nonlinear beams produced by three source transducers. An analytic form for the spatial averaging filter, applicable to Gaussian harmonic beams, is derived. A second analytic form for the spatial averaging filter, applicable to quadratic harmonic beams, is derived by extending the spatial averaging correction recommended by IEC 62127-1 Annex E to nonlinear signals with multiple harmonics. Both forms are applicable to all hydrophones (not just needle and fiber-optic hydrophones). Simulation analysis performed for a wide variety of transducer geometries indicates that the Gaussian spatial averaging filter formula is more accurate than the quadratic formula over a wider range of harmonics. Additional experimental validation is provided in Part II. Readers who are uninterested in hydrophone theory may skip the theoretical and experimental sections of this paper and proceed to the graphical guide for practical information to inform and support selection of hydrophone sensitive element size (but might be well advised to read the Introduction).

Characterization of laser ultrasound source signals in biological tissues for imaging applications

Short optical pulses emitted from a tunable Q-switched laser (800 to 2000 nm) generate laser ultrasound (LUS) signals at the surface of biological tissue. The LUS signal's acoustic frequency content, dependence on sample type, and optical wavelength are observed in the far field. The experiments yield a reference dataset for the design of noncontact LUS imaging systems. Measurements show that the majority of LUS signal energy in biological tissues is within the 0.5 and 3 MHz frequency bands and the total acoustic energy generated increases with the optical absorption coefficient of water, which governs tissue optical absorption in the infrared range. The experimental results also link tissue surface roughness and acoustic attenuation with limited LUS signal bandwidth in biological tissue. Images constructed using 810-, 1064-, 1550-, and 2000-nm generation laser wavelengths and a contact piezoelectric receiver demonstrates the impact of the generation laser wavelength on image quality. A noncontact LUS-based medical imaging system has the potential to be an effective medical imaging device. Such a system may mitigate interoperator variability associated with current medical ultrasound imaging techniques and expand the scope of imaging applications for ultrasound.

Directivity and Frequency-Dependent Effective Sensitive Element Size of a Reflectance-Based Fiber-Optic Hydrophone: Predictions From Theoretical Models Compared With Measurements

The goal of this work was to measure the directivity of a reflectance-based fiber-optic hydrophone at multiple frequencies and to compare it to four theoretical models: rigid baffle (RB), rigid piston (RP), unbaffled (UB), and soft baffle (SB). The fiber had a nominal 105- [Formula: see text] diameter core and a 125- [Formula: see text] overall diameter (core + cladding). Directivity measurements were performed at 2.25, 3.5, 5, 7.5, 10, and 15 MHz from ±90° in two orthogonal planes. Effective hydrophone sensitive element radius was estimated by least-squares fitting the four models to the directivity measurements using the sensitive element radius as an adjustable parameter. Over the range from 2.25 to 15 MHz, the average magnitudes of differences between the effective and nominal sensitive element radii were 59% ± 49% (RB), 10% ± 5% (RP), 46% ± 38% (UB), and 71% ± 19% (SB). Therefore, the directivity of a reflectance-based fiber-optic hydrophone may be best estimated by the RP model.

Directivity and Frequency-Dependent Effective Sensitive Element Size of Needle Hydrophones: Predictions From Four Theoretical Forms Compared With Measurements

Directivity is a hydrophone specification that describes response as a function of angle of incidence. The goal of this study was to compare, in the context of needle hydrophones, the commonly used rigid baffle model for hydrophone directivity to three alternative models: soft baffle, unbaffled (UB), and rigid piston (RP). Directivity measurements were performed at 1, 2, 3, 4, 6, 8, and 10 MHz from ±7° in two orthogonal planes for two ceramic and two polymer needle hydrophones with nominal geometrical sensitive element diameters of 200, 400, 600, and 1000 . Effective hydrophone sensitive element radius was estimated by least-squares fitting the four models to directivity measurement data using the sensitive element radius (a) as an adjustable parameter. For > 4 (where and = wavelength), the RP model outperformed the other three models. For , the average error in estimated sensitive element radius was 7% [95% confidence interval (CI): 3%-12%] for the RP model while the lowest average error by the other three models was 46% (95% CI: 38%-54%) for the UB model.

2D Analytical Model for the Directivity Prediction of Ultrasonic Contact Type Transducers in the Generation of Guided Waves

In this paper, a novel 2D analytical model based on the Huygens’s principle of wave propagation is proposed in order to predict the directivity patterns of contact type ultrasonic transducers in the generation of guided waves (GWs). The developed model is able to estimate the directivity patterns at any distance, at any excitation frequency and for any configuration and shape of the transducers with prior information of phase dispersive characteristics of the guided wave modes and the behavior of transducer. This, in turn, facilitates to choose the appropriate transducer or arrays of transducers, suitable guided wave modes and excitation frequency for the nondestructive testing (NDT) and structural health monitoring (SHM) applications. The model is demonstrated for P1-type macro-fiber composite (MFC) transducer glued on a 2 mm thick aluminum (Al) alloy plate. The directivity patterns of MFC transducer in the generation of fundamental guided Lamb modes (the S0 and A0) and shear horizontal mode (the SH0) are successfully obtained at 80 kHz, 5-period excitation signal. The results are verified using 3D finite element (FE) modelling and experimental investigation. The results obtained using the proposed model shows the good agreement with those obtained using numerical simulations and experimental analysis. The calculation time using the analytical model was significantly shorter as compared to the time spent in experimental analysis and FE numerical modelling.

A Method for Measuring the Directional Response of Ultrasound Receivers in the Range 0.3-80 MHz Using a Laser-Generated Ultrasound Source

A simple method for measuring the directivity of an ultrasound receiver is described. The method makes use of a custom-designed laser ultrasound source which generates a large diameter (>1 cm) broadband monopolar plane wave with a continuous frequency content extending from to . The plane wave is highly uniform in amplitude (±5% over >8 mm) and phase (equivalent to at 80 MHz over ). To measure directivity, the source is rotated around the receiver under test in a compact centimeter-scale setup. To demonstrate the method, it was used to measure the directivity of two broadband small aperture Fabry-Perot ultrasound sensors over an angular range of ±50° at frequencies up to 80 MHz. Measurements were found to be highly repeatable with an estimated typical repeatability <4% in the range of 0.5-25 MHz. Due to the broad bandwidth, large size, and adjustable nature of the source, the method is widely applicable and could aid the characterization of receivers used in medical ultrasound, ultrasonic nondestructive testing. and ultrasound metrology.

Nonconfocal all-optical laser-ultrasound and photoacoustic imaging system for angle-dependent deep tissue imaging

Biomedical imaging systems incorporating both photoacoustic (PA) and ultrasound capabilities are of interest for obtaining optical and acoustic properties deep in tissue. While most dual-modality systems utilize piezoelectric transducers, all-optical systems can obtain broadband high-resolution data with hands-free operation. Previously described reflection-mode all-optical laser-ultrasound (LUS) systems use a confocal source and detector; however, angle-dependent raypaths are lost in this configuration. As a result, the overall imaging aperture is reduced, which becomes increasingly problematic with depth. We present a reflection-mode nonconfocal LUS and PA imaging system that uses signals recorded on all-optical hardware to create angle-dependent images. We use reverse-time migration and time reversal to reconstruct the LUS and PA images. We demonstrate this methodology with both a numerical model and tissue phantom experiment to image a steep-curvature vessel with a limited aperture 2-cm beneath the surface. Nonconfocal imaging demonstrates improved focusing by 30% and 15% compared to images acquired with a single LUS source in the numerical and experimental LUS images, respectively. The appearance of artifacts is also reduced. Complementary PA images are straightforward to acquire with the nonconfocal system by tuning the source wavelength and can be further developed for quantitative multiview PA imaging.

Joint Reconstruction of Absorbed Optical Energy Density and Sound Speed Distributions in Photoacoustic Computed Tomography: A Numerical Investigation

Photoacoustic computed tomography (PACT) is a rapidly emerging bioimaging modality that seeks to reconstruct an estimate of the absorbed optical energy density within an object. Conventional PACT image reconstruction methods assume a constant speed-of-sound (SOS), which can result in image artifacts when acoustic aberrations are significant. It has been demonstrated that incorporating knowledge of an object's SOS distribution into a PACT image reconstruction method can improve image quality. However, in many cases, the SOS distribution cannot be accurately and/or conveniently estimated prior to the PACT experiment. Because variations in the SOS distribution induce aberrations in the measured photoacoustic wavefields, certain information regarding an object's SOS distribution is encoded in the PACT measurement data. Based on this observation, a joint reconstruction (JR) problem has been proposed in which the SOS distribution is concurrently estimated along with the sought-after absorbed optical energy density from the photoacoustic measurement data. A broad understanding of the extent to which the JR problem can be accurately and reliably solved has not been reported. In this work, a series of numerical experiments is described that elucidate some important properties of the JR problem that pertain to its practical feasibility. To accomplish this, an optimization-based formulation of the JR problem is developed that yields a non-linear iterative algorithm that alternatively updates the two image estimates. Heuristic analytic insights into the reconstruction problem are also provided. These results confirm the ill-conditioned nature of the joint reconstruction problem that will present significant challenges for practical applications.

Three-Dimensional Optoacoustic and Laser-Induced Ultrasound Tomography System for Preclinical Research in Mice: Design and Phantom Validation

In this work, we introduce a novel three-dimensional imaging system for in vivo high-resolution anatomical and functional whole-body visualization of small animal models developed for preclinical and other type of biomedical research. The system (LOUIS-3DM) combines a multiwavelength optoacoustic tomography (OAT) and laser-induced ultrasound tomography (LUT) to obtain coregistered maps of tissue optical absorption and speed of sound, displayed within the skin outline of the studied animal. The most promising applications of the LOUIS-3DM include 3D angiography, cancer research, and longitudinal studies of biological distributions of optoacoustic contrast agents.

A Constrained Variable Projection Reconstruction Method for Photoacoustic Computed Tomography Without Accurate Knowledge of Transducer Responses

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the absorbed optical energy density within tissue. When the imaging system employs conventional piezoelectric ultrasonic transducers, the ideal photoacoustic (PA) signals are degraded by the transducers' acousto-electric impulse responses (EIRs) during the measurement process. If unaccounted for, this can degrade the accuracy of the reconstructed image. In principle, the effect of the EIRs on the measured PA signals can be ameliorated via deconvolution; images can be reconstructed subsequently by application of a reconstruction method that assumes an idealized EIR. Alternatively, the effect of the EIR can be incorporated into an imaging model and implicitly compensated for during reconstruction. In either case, the efficacy of the correction can be limited by errors in the assumed EIRs. In this work, a joint optimization approach to PACT image reconstruction is proposed for mitigating errors in reconstructed images that are caused by use of an inaccurate EIR. The method exploits the bi-linear nature of the imaging model and seeks to refine the measured EIR during the process of reconstructing the sought-after absorbed optical energy density. Computer-simulation and experimental studies are conducted to investigate the numerical properties of the method and demonstrate its value for mitigating image distortions and enhancing the visibility of fine structures.

Discrete imaging models for three-dimensional optoacoustic tomography using radially symmetric expansion functions

Optoacoustic tomography (OAT), also known as photoacoustic tomography, is an emerging computed biomedical imaging modality that exploits optical contrast and ultrasonic detection principles. Iterative image reconstruction algorithms that are based on discrete imaging models are actively being developed for OAT due to their ability to improve image quality by incorporating accurate models of the imaging physics, instrument response, and measurement noise. In this work, we investigate the use of discrete imaging models based on Kaiser-Bessel window functions for iterative image reconstruction in OAT. A closed-form expression for the pressure produced by a Kaiser-Bessel function is calculated, which facilitates accurate computation of the system matrix. Computer-simulation and experimental studies are employed to demonstrate the potential advantages of Kaiser-Bessel function-based iterative image reconstruction in OAT.

Fast spatiotemporal image reconstruction based on low-rank matrix estimation for dynamic photoacoustic computed tomography

In order to monitor dynamic physiological events in near-real time, a variety of photoacoustic computed tomography (PACT) systems have been developed that can rapidly acquire data. Previously reported studies of dynamic PACT have employed conventional static methods to reconstruct a temporally ordered sequence of images on a frame-by-frame basis. Frame-by-frame image reconstruction (FBFIR) methods fail to exploit correlations between data frames and are known to be statistically and computationally suboptimal. In this study, a low-rank matrix estimation-based spatiotemporal image reconstruction (LRME-STIR) method is investigated for dynamic PACT applications. The LRME-STIR method is based on the observation that, in many PACT applications, the number of frames is much greater than the rank of the ideal noiseless data matrix. Using both computer-simulated and experimentally measured photoacoustic data, the performance of the LRME-STIR method is compared with that of conventional FBFIR method followed by image-domain filtering. The results demonstrate that the LRME-STIR method is not only computationally more efficient but also produces more accurate dynamic PACT images than a conventional FBFIR method followed by image-domain filtering.

Using optoacoustic imaging for measuring the temperature dependence of Grüneisen parameter in optically absorbing solutions

Grüneisen parameter is a key temperature-dependent physical characteristic responsible for thermoelastic efficiency of materials. We propose a new methodology for accurate measurements of temperature dependence of Grüneisen parameter in optically absorbing solutions. We use two-dimensional optoacoustic (OA) imaging to improve accuracy of measurements. Our approach eliminates contribution of local optical fluence and absorbance. To validate the proposed methodology, we studied temperature dependence of aqueous cupric sulfate solutions in the range from 22 to 4 °C. Our results for the most diluted salt perfectly matched known temperature dependence for the Grüneisen parameter of water. We also found that Grüneisen-temperature relationship for cupric sulfate exhibits linear trend with respect to the concentration. In addition to accurate measurements of Grüneisen changes with temperature, the developed technique provides a basis for future high precision OA temperature monitoring in live tissues.

Optoacoustic multispectral imaging of radiolucent foreign bodies in tissue

Optoacoustic imaging is an emerging medical technology that uniquely combines the absorption contrast of optical imaging and the penetration depth of ultrasound. While it is not currently employed as a clinical imaging modality, the results of current research strongly support the use of optoacoustic-based methods in medical imaging. One such application is the diagnosis of the presence of soft tissue foreign bodies. Because many radiolucent foreign bodies have sufficient contrast for imaging in the optical domain, laser-induced optoacoustic imaging could be advantageous for the detection of such objects. Common foreign bodies have been scanned over a range of visible and near infrared wavelengths by using an optoacoustic method to obtain the spectroscopic properties of the materials commonly associated with these foreign bodies. The derived optical absorption spectra compared quite closely to the absorption spectra generated when using a conventional spectrophotometer. By using the probe-beam deflection technique, a novel, pressure-wave detection method, we successfully generated optoacoustic spectroscopic plots of a wooden foreign body embedded in a tissue phantom, which closely resembled the spectrum of the same object obtained in isolation. A practical application of such spectra is to assemble a library of spectroscopic data for radiolucent materials, from which specific characteristic wavelengths can be selected for use in optimizing imaging instrumentation and provide a basis for the identification of the material properties of particular foreign bodies.

An imaging model incorporating ultrasonic transducer properties for three-dimensional optoacoustic tomography

Optoacoustic tomography (OAT) is a hybrid imaging modality that combines the advantages of optical and ultrasound imaging. Most existing reconstruction algorithms for OAT assume that the ultrasound transducers employed to record the measurement data are point-like. When transducers with large detecting areas and/or compact measurement geometries are utilized, this assumption can result in conspicuous image blurring and distortions in the reconstructed images. In this work, a new OAT imaging model that incorporates the spatial and temporal responses of an ultrasound transducer is introduced. A discrete form of the imaging model is implemented and its numerical properties are investigated. We demonstrate that use of the imaging model in an iterative reconstruction method can improve the spatial resolution of the optoacoustic images as compared to those reconstructed assuming point-like ultrasound transducers.