Non-Gaussian statistics and temporal variations of the ultrasound signal backscattered by blood at frequencies between 10 and 58 MHz

Very little is known about the blood backscattering behavior and signal statistics following flow stoppage at frequencies higher than 10 MHz. Measurements of the radio frequency (rf) signals backscattered by normal human blood (hematocrit = 40%, temperature = 37 degrees C) were performed in a tube flow model at mean frequencies varying between 10 and 58 MHz. The range of increase of the backscattered power during red blood cell (RBC) rouleau formation was close to 15 dB at 10 and 36 MHz, and dropped, for the same blood samples, below 8 dB at 58 MHz. Increasing the frequency from 10 to 58 MHz raised the slope of the power changes at the beginning of the kinetics of aggregation, and could emphasize the non-Gaussian behavior of the rf signals interpreted in terms of the K and Nakagami statistical models. At 36 and 58 MHz, significant increases of the kurtosis coefficient, and significant reductions of the Nakagami parameter were noted during the first 30 s of flow stoppage. In conclusion, increasing the transducer frequency reduced the magnitude of the backscattered power changes attributed to the phenomenon of RBC aggregation, but improved the detection of rapid growth in aggregate sizes and non-Gaussian statistical behavior.     

Experimental ultrasound characterization of red blood cell aggregation using the structure factor size estimator

The frequency dependence of the ultrasonic backscattering coefficient (BSC) was studied to assess the level of red blood cell (RBC) aggregation. Three monoelement focused wideband transducers were used to insonify porcine blood sheared in a Couette flow from 9 to 30 MHz. A high shear rate was first applied to promote disaggregation. Different residual shear rates were then used to promote formation of RBC aggregates. The structure factor size estimator (SFSE), a second-order data reduction model based on the structure factor, was applied to the frequency-dependent BSC. Two parameters were extracted from the model to describe the level of aggregation at 6% and 40% hematocrits: W, the packing factor, and D the aggregate diameter, expressed in number of RBCs. Both parameters closely matched theoretical values for nonaggregated RBCs. W and D increased during aggregation with stabilized values modulated by the applied residual shear rate. Furthermore, parameter D during the kinetics of aggregation at 6% hematocrit under static conditions correlated with an optical RBC aggregate size estimation from microscopic images (r(2)=0.76). To conclude, the SFSE presents an interesting framework for tissue characterization of partially correlated dense tissues such as aggregated RBCs.     

Effects of red blood cell aggregates dissociation on the estimation of ultrasound speckle image velocimetry

Ultrasound speckle image of blood is mainly attributed by red blood cells (RBCs) which tend to form RBC aggregates. RBC aggregates are separated into individual cells when the shear force is over a certain value. The dissociation of RBC aggregates has an influence on the performance of ultrasound speckle image velocimetry (SIV) technique in which a cross-correlation algorithm is applied to the speckle images to get the velocity field information. The present study aims to investigate the effect of the dissociation of RBC aggregates on the estimation quality of SIV technique. Ultrasound B-mode images were captured from the porcine blood circulating in a mock-up flow loop with varying flow rate. To verify the measurement performance of SIV technique, the centerline velocity measured by the SIV technique was compared with that measured by Doppler spectrograms. The dissociation of RBC aggregates was estimated by using decorrelation of speckle patterns in which the subsequent window was shifted as much as the speckle displacement to compensate decorrelation caused by in-plane loss of speckle patterns. The decorrelation of speckles is considerably increased according to shear rate. Its variations are different along the radial direction. Because the dissociation of RBC aggregates changes ultrasound speckles, the estimation quality of SIV technique is significantly correlated with the decorrelation of speckles. This degradation of measurement quality may be improved by increasing the data acquisition rate. This study would be useful for simultaneous measurement of hemodynamic and hemorheological information of blood flows using only speckle images.     

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.     

Acoustic characterization of echogenic liposomes: frequency-dependent attenuation and backscatter

Ultrasound contrast agents (UCAs) are used clinically to aid detection and diagnosis of abnormal blood flow or perfusion. Characterization of UCAs can aid in the optimization of ultrasound parameters for enhanced image contrast. In this study echogenic liposomes (ELIPs) were characterized acoustically by measuring the frequency-dependent attenuation and backscatter coefficients at frequencies between 3 and 30 MHz using a broadband pulse-echo technique. The experimental methods were initially validated by comparing the attenuation and backscatter coefficients measured from 50-μm and 100-μm polystyrene microspheres with theoretical values. The size distribution of the ELIPs was measured and found to be polydisperse, ranging in size from 40 nm to 6 μm in diameter, with the highest number observed at 65 nm. The ELIP attenuation coefficients ranged from 3.7 ± 1.0 to 8.0 ± 3.3 dB/cm between 3 and 25 MHz. The backscatter coefficients were 0.011 ± 0.006 (cm str)(-1) between 6 and 9 MHz and 0.023?±?0.006 (cm str)(-1) between 13 and 30 MHz. The measured scattering-to-attenuation ratio ranged from 8% to 22% between 6 and 25 MHz. Thus ELIPs can provide enhanced contrast over a broad range of frequencies and the scattering properties are suitable for various ultrasound imaging applications including diagnostic and intravascular ultrasound.     

High frequency attenuation measurements of lipid encapsulated contrast agents

A number of recent studies have indicated the potential of ultrasound contrast agent imaging at high ultrasound frequencies. However, the acoustic properties of microbubbles at frequencies above 10 MHz remain poorly understood at present. In this study we characterize the high frequency attenuation properties of (1) BR14, (2) BR14 that has been mechanically filtered (1 and 2 microm pore sizes) to exclude larger bubbles, and (3) the micron to submicron agent BG2423. A narrowband pulse-echo substitution method is employed with a series of four transducers covering the frequency range from 2 to 50 MHz. For BR14, attenuation decreases rapidly from 2 to 10 MHz and then more gradually from 10 to 50 MHz. For 2 microm filtration, the attenuation peaks between 10 and 15 MHz. For 1 microm filtration, attenuation continues to rise until 50 MHz. The agent BG2423 exhibits a diffuse attenuation peak in the range of 15-25 MHz and remains high until 50 MHz. These results demonstrate a strong influence of bubble size on high frequency attenuation curves, with bubble diameters of 1-2 microm and below having more pronounced acoustic activity at frequencies above 10 MHz.     

Influence of cell packing by centrifugation on 40-MHz ultrasound backscatter

High-frequency ultrasound (HFUS) signals backscattered from RBL-2H3 cell pellets prepared under different centrifugal forces were analyzed to investigate the packing effect of cell aggregates. The measurements were performed in a pulse-echo setup with a 40-MHz transducer. The changes of ultrasound signals from cell pellet in backscattered power, statistical parameter, and pellet thickness were monitored after centrifugation at between 100g and 1600g. Experimental results showed that the HFUS backscattered power from cell pellets was inversely proportional to centrifugal force and increased to a plateau within 1-2 h after centrifugation. The initial thickness of cell pellets decreased with higher centrifugal force, but the changes in thickness and time that took to reach a plateau increased at higher centrifugal force. The envelope statistics of backscattered signals with Nakagami distribution indicates that the centrifugal force and elapsed time after centrifugation affected the backscattering characteristics. The present study suggests that centrifugal force and data acquisition time after cell pellet formation should be considered in in vitro cell packing method with centrifugation to emulate the tissue in vivo.     

Anisotropy of the backscatter coefficient of formalin-fixed ovine myocardium

The objective of this study was to measure the backscatter coefficient of formalin-fixed myocardial tissue as a function of angle of insonification relative to the myocardial fiber direction. Backscatter measurements were performed on eight cylindrical formalin-fixed lamb myocardial specimens and compensated for attenuation and diffraction effects to determine the backscatter coefficient. The backscatter coefficient at 5 MHz was found to be maximum for insonification perpendicular to the predominant myofiber orientation and minimum for parallel insonification, with values of (17+/-14) and (1.2+/-0.7) x 10(-4) cm(-1) sr(-1) (mean+/-standard deviation), respectively.     

Frequency dependence of ultrasonic backscatter from human trabecular bone: theory and experiment

A model describing the frequency dependence of backscatter coefficient from trabecular bone is presented. Scattering is assumed to originate from the surfaces of trabeculae, which are modeled as long thin cylinders with radii small compared with the ultrasonic wavelength. Experimental ultrasonic measurements at 500 kHz, 1 MHz, and 2.25 MHz from a wire target and from trabecular bone samples from human calcaneus in vitro are reported. In both cases, measurements are in good agreement with theory. For mediolateral insonification of calcaneus at low frequencies, including the typical diagnostic range (near 500 kHz), backscatter coefficient is proportional to frequency cubed. At higher frequencies, the frequency response flattens out. The data also suggest that at diagnostic frequencies, multiple scattering effects on the average are relatively small for the samples investigated. Finally, at diagnostic frequencies, the data suggest that absorption is likely to be a larger component of attenuation than scattering.     

High-frequency ultrasound backscattering by blood: analytical and semianalytical models of the erythrocyte cross section

This paper proposes analytical and semianalytical models of the ultrasonic backscattering cross section (BCS) of various geometrical shapes mimicking a red blood cell (RBC) for frequencies varying from 0 to 90 MHz. By assuming the first-order Born approximation and by modeling the shape of a RBC by a realistic biconcave volume, different scattering behaviors were identified for increasing frequencies. For frequencies below 18 MHz, a RBC can be considered a Rayleigh scatterer. For frequencies less than 39 MHz, the general concept of acoustic inertia tensor is introduced to describe the variation of the BCS with the frequency and the incidence direction. For frequencies below 90 MHz, ultrasound backscattering by a RBC is equivalent to backscattering by a cylinder of height 2 microm and diameter 7.8 microm. These results lay the basis of ultrasonic characterization of RBC aggregation by proposing a method that distinguishes the contribution of the individual RBC acoustical characteristics from collective effects, on the global blood backscattering coefficient. A new method of data reduction that models the frequency dependence of the ultrasonic BCS of micron-sized weak scatterers is also proposed. Applications of this method are in tissue characterization as well as in hematology.     

Assessment of accuracy of the structure-factor-size-estimator method in determining red blood cell aggregate size from ultrasound spectral backscatter coefficient

A computer simulation study to produce ultrasonic backscatter coefficients (BSCs) from red blood cell (RBC) clusters is discussed. The simulation algorithm is suitable for generating non-overlapping, isotropic, and fairly identical RBC clusters. RBCs were stacked following the hexagonal close packing (HCP) structure to form a compact spherical aggregate. Such an aggregate was repeated and placed randomly under non-overlapping condition in the three-dimensional space to mimic an aggregated blood sample. BSCs were computed between 750 KHz and 200 MHz for samples of various cluster sizes at different hematocrits. Magnitudes of BSCs increased with mean aggregate sizes at low frequencies (<20 MHz). The accuracy of the structure-factor-size-estimator (SFSE) method in determining mean aggregate size and packing factor was also examined. A good correlation (R(2) ≥ 0.94) between the mean size of aggregates predicted by the SFSE and true size was found for each hematocrit. This study shows that for spherical aggregates there exists a region for each hematocrit where SFSE works most accurately. Typically, error of SFSE in estimating mean cluster size was <20% for dimensions between 14 and 17 μm at 40% hematocrit. This study suggests that the theoretical framework of SFSE is valid under the assumption of isotropic aggregates.     

Ultrasound attenuation estimation using backscattered echoes from multiple sources

The objective of this study was to devise an algorithm that can accurately estimate the attenuation along the propagation path (i.e., the total attenuation) from backscattered echoes. It was shown that the downshift in the center frequency of the backscattered ultrasound echoes compared to echoes obtained in a water bath was calculated to have the form Deltaf=mf(o)+b after normalizing with respect to the source bandwidth where m depends on the correlation length, b depends on the total attenuation, and f(o) is the center frequency of the source as measured from a reference echo. Therefore, the total attenuation can be determined independent of the scatterer correlation length by measuring the downshift in center frequency from multiple sources (i.e., different f(o)) and fitting a line to the measured shifts versus f(o). The intercept of the line gives the total attenuation along the propagation path. The calculations were verified using computer simulations of five spherically focused sources with 50% bandwidths and center frequencies of 6, 8, 10, 12, and 14 MHz. The simulated tissue had Gaussian scattering structures with effective radii of 25 mum placed at a density of 250 mm(3). The attenuation of the tissue was varied from 0.1 to 0.9 dB / cm-MHz. The error in the attenuation along the propagation path ranged from -3.5+/-14.7% for a tissue attenuation of 0.1 dB / cm-MHz to -7.0+/-3.1% for a tissue attenuation of 0.9 dB / cm-MHz demonstrating that the attenuation along the propagation path could be accurately determined using backscattered echoes from multiple sources using the derived algorithm.     

Ultrasound generation with high power and coil only EMAT concepts

Electro-magnetic acoustic transducers (EMATs) are intended as non-contact and non-destructive ultrasound transducers for metallic material. The transmitted intensities from EMATS are modest, particularly at notable lift off distances. Some time ago a concept for a “coil only EMAT” was presented, without static magnetic field. In this contribution, such compact “coil only EMATs” with effective areas of 1–5 cm² were driven to excessive power levels at MHz frequencies, using pulsed power technologies. RF induction currents of 10 kA and tens of Megawatts are applied. With increasing power the electroacoustic conversion efficiency also increases. The total effect is of second order or quadratic, therefore non-linear and progressive, and yields strong ultrasound signals up to kW/cm² at MHz frequencies in the metal. Even at considerable lift off distances (cm) the ultrasound can be readily detected. Test materials are aluminum, ferromagnetic steel and stainless steel (non-ferromagnetic). Thereby, most metal types are represented. The technique is compared experimentally with other non-contact methods: laser pulse induced ultrasound and spark induced ultrasound, both damaging to the test object’s surface. At small lift off distances, the intensity from this EMAT concept clearly outperforms the laser pulses or heavy spark impacts.     

Ultrasonic visualization of dynamic behavior of red blood cells in flowing blood

It is well known that the scatter of ultrasound by blood is mainly attributed to red blood cells (RBCs) and RBC aggregation. In the present review, researches of hemodynamic influence on RBC aggregation and ultrasound backscatter from blood were overviewed. A mock flow loop and a cylindrical chamber were employed to produce various blood flows, such as pulsatile, oscillatory, and rotational flow. The “black hole” (BLH), a dark hole at the tube center surrounded by bright zone in the cross sectional B-mode image and “bright collapsing ring” (BRCR) phenomena, appearance of bright ring at the periphery and collapse of it at the center during a pulsatile cycle, were observed under pulsatile flow. The combined effects of shear rate and flow acceleration on RBC aggregation were suggested as a possible mechanism for these phenomena. The stroke volume-dependence of the “bright ring” phenomenon under oscillatory flow could also be explained by flow acceleration. The enveloped echo images from rotational flow in a compact blood chamber showed the spatial and temporal variations of RBC aggregation, which varied with the mammalian species. In the stenotic model, it was found that the echogenic variation increased locally at a distance of three tube diameters downstream from the stenosis during decelerating period, which was proposed to be mainly due to flow turbulence. The similar ‘bright ring’ was also observed fromin vivo human carotid artery in harmonic imaging.     

Accurate coarse-grained modeling of red blood cells

We develop a systematic coarse-graining procedure for modeling red blood cells (RBCs) using arguments based on mean-field theory. The three-dimensional RBC membrane model takes into account the bending energy, in-plane shear energy, and constraints of fixed surface area and fixed enclosed volume. The coarse-graining procedure is general, it can be used for arbitrary level of coarse-graining and does not employ any fitting parameters. The sensitivity of the coarse-grained model is investigated and its behavior is validated against available experimental data and in dissipative particle dynamics (DPD) simulations of RBCs in capillary and shear flows.     

Numerical analysis of ultrasonic transmission and absorption of oblique plane waves through the human skull.

Ultrasonic transmission and absorption of oblique plane waves through the human skull are analyzed numerically for frequencies ranging from 1/2 to 1 MHz. These frequencies are optimum for noninvasive ultrasound therapy of brain disorders where numerical predictions of skull transmission are used to set the phase and amplitude of source elements in the phased array focusing system. The idealized model of the skull is a three-layer solid with ivory outer and inner layers and a middle marrow layer. Each layer is modeled as a flat, homogeneous, isotropic, linear solid with effective complex wave speeds to account for focused energy losses due to material damping and scattering. The model is used to predict the amplitude and phase of the transmitted wave and volumetric absorption. Results are reported for three different skull thicknesses: 3 mm, 6 mm, and 9 mm. Thickness resonances are observed in the transmitted wave for 3 mm skulls at all frequencies and for the 6 mm skulls below 0.75 MHz. Otherwise, the transmission is dominated by the direct wave. Skull phase errors due to shear waves are shown to minimally degrade the power at the focus for angles of incidence up to 20 degrees from normal even for low material damping. The location of the peak volumetric absorption occurs either in the outer ivory or middle marrow layer and shown to vary due to wave interference.     

Numerical simulation of the red blood cell aggregation and deformation behaviors in ultrasonic field

The objective of this paper is to propose an immersed boundary lattice Boltzmann method (IB-LBM) considering the ultrasonic effect to simulate red blood cell (RBC) aggregation and deformation in ultrasonic field. Numerical examples involving the typical streamline, normalized out-of-plane vorticity contours and vector fields in pure plasma under three different ultrasound intensities are presented. Meanwhile, the corresponding transient aggregation behavior of RBCs, with special emphasis on the detailed process of RBC deformation, is shown. The numerical results reveal that the ultrasound wave acted on the pure plasma can lead to recirculation flow, which contributes to the RBCs aggregation and deformation in microvessel. Furthermore, increasing the intensity of the ultrasound wave can significantly enhance the aggregation and deformation of the RBCs. And the formation of the RBCs aggregation leads to the fluctuated and dropped vorticity value of plasma in return.     

H-22细胞声孔效应的应力阈值

理论和实验研究了超声空化场中的H-22型肝癌细胞产生可逆声孔效应的剪应力阈值.本文用1.37 MHz的聚焦声场,当超顺磁性纳米氧化铁在细胞悬液中的终浓度为410 μg/mL,换能器负载电功率为2 W,超声辐照60 s,细胞存活率90%以上时,有45.9±13.5%的细胞显示普鲁士蓝染阳性,暗示超声作用下,这些细胞表面曾出现可逆性微孔而使磁性微粒由此进入细胞内.利用无界自由空间微泡运动方程的球对称稳态解对实验条件下细胞膜表面的切变应力进行数值估算,结果表明,使H-22细胞产生可逆性声孔效应的微流剪应力阈值为697 Pa. 关键词： 声孔效应 磁性标记 微流 剪应力     

Effect of ultrasound and temperature on tomato peroxidase

Tomato peroxidase (POD) was inactivated by heat and ultrasound. Thermal inactivation of tomato POD was performed at temperatures of 63, 64, 65, 66 and 67°C. Thermal inactivation of POD in tomato showed apparent first-order kinetics. E(a) was calculated from the slopes of Arrhenius plot and found as 14.9×10(4) J/mol. The effect of ultrasound on tomato peroxidase inactivation was investigated at 15%, 25%, 40%, 50% and 75% ultrasonic powers for 20-150 s. It was observed that as the ultrasonic power increased, inactivation rate increased. 100% POD inactivation was observed at 50% power for 150 s and at 75% power for 90 s of ultrasonication. Regeneration of POD activity was investigated for the samples exposed to ultrasound at different ultrasonic powers. It was observed that at 15% and 25% ultrasonic powers, residual enzyme activity increased linearly, whereas at powers of 40%, 50% and 75% residual enzyme activity increased nonlinearly. There was no regeneration in the samples which 100% enzyme inactivation obtained by ultrasound. A significant decrease in vitamin C content of tomato extract was observed as a function of temperature and treatment time, whereas ultrasonic treatment had no significant effect on vitamin C content of tomato extract.     

Determination of object resonances by vibro-acoustography and their associated modes

Vibro-acoustography technique known by its noncontact excitation was used to detect resonance frequencies of objects in water. Two intersecting ultrasound beams generated by a 40 mm-diameter annular array transducer, focused at 35 mm and driven at f1=2.2 MHz and f2=2.22 MHz respectively, were targeted inside the object under test to produce a radiation force beating at the difference frequency f2-f1. This low frequency radiation force was used to excite the resonance vibration modes of the object by sweeping the frequency f2 between 2.22 and 2.275 MHz. The amplitude of the acoustic emission produced by the vibrations of the object was detected by a low frequency hydrophone (BW=60 kHz). By this approach, it was possible to detect resonance frequencies through amplitude variations of the measured acoustic emission. Experiments were conducted in a water tank for objects of different shapes and sizes. With a chalk sphere (15 mm-diameter) two resonance frequencies were detected at 45.75 and 68.75 kHz, and with a cylinder (10.38 mm-diameter and 32.20 mm-length) four principal resonance frequencies were identified in the 60 kHz-bandwidth of the hydrophone. It was shown with finite element calculations performed with Ansys, in which both solid and fluid parts were modelled, that the measured resonance frequencies corresponded to compressional or dilatation vibration modes of the object. It was verified that shear waves generated by torsional vibration modes were not propagated in water, as it is well known. The use of this technique to characterize heterogeneities in different media seems to be relatively more advantageous to other ultrasonic methods.