On the applicability of Kramers-Kronig relations for ultrasonic attenuation obeying a frequency power law

An atlas of the heterogeneous viscoelastic brain with local power-law attenuation synthesised using Prony-series

This review addresses the acute need to acknowledge the mechanical heterogeneity of brain matter and to accurately calibrate its local viscoelastic material properties accordingly. Specifically, it is important to compile the existing and disparate literature on attenuation power-laws and dispersion to make progress in wave physics of brain matter, a field of research that has the potential to explain the mechanisms at play in diffuse axonal injury and mild traumatic brain injury in general. Currently, viscous effects in the brain are modelled using Prony-series, i.e., a sum of decaying exponentials at different relaxation times. Here we collect and synthesise the Prony-series coefficients appearing in the literature for twelve regions: brainstem, basal ganglia, cerebellum, corona radiata, corpus callosum, cortex, dentate gyrus, hippocampus, thalamus, grey matter, white matter, homogeneous brain, and for eight different mammals: pig, rat, human, mouse, cow, sheep, monkey and dog. Using this data, we compute the fractional-exponent attenuation power-laws for different tissues of the brain, the corresponding dispersion laws resulting from causality, and the averaged Prony-series coefficients. STATEMENT OF SIGNIFICANCE: Traumatic brain injuries are considered a silent epidemic and finite element methods (FEMs) are used in modelling brain deformation, requiring access to viscoelastic properties of brain. To the best of our knowledge, this work presents 1) the first multi-frequency viscoelastic atlas of the heterogeneous brain, 2) the first review focusing on viscoelastic modelling in both FEMs and experimental works, 3) the first attempt to conglomerate the disparate existing literature on the viscoelastic modelling of the brain and 4) the largest collection of viscoelastic parameters for the brain (212 different Prony-series spanning 12 different tissues and 8 different animal surrogates). Furthermore, this work presents the first brain atlas of attenuation power-laws essential for modelling shear waves in brain.

Three-dimensional time-domain Green's functions and spatial impulse responses for the van Wijngaarden wave equation

An exact analytical three-dimensional time-domain Green's function is introduced for the van Wijngaarden wave equation when the coefficients of the two loss terms satisfy a specific relationship. This analytical Green's function, which describes frequency-squared attenuation in acoustic media such as water, enables the subsequent derivation of new expressions that describe the lossy spatial impulse response for a circular piston. Initial time-domain assessments, which compare the Green's functions for the van Wijngaarden, Stokes, and power law wave equations using the attenuation and sound speed for water, indicate that these three lossy wave equations yield nearly identical results at distances greater than or equal to 10 μm. Lossy spatial impulse responses are also evaluated with increasing distance in and near the paraxial region of a circular piston radiating in water to reveal some interesting time-domain interactions between frequency-squared attenuation and diffraction. Similar behaviors are also demonstrated for the lossy far-field spatial impulse. In addition, the convergence is demonstrated for two analytically equivalent expressions applied to numerical computations of the lossy spatial impulse response. The results show that these new expressions are ideal for describing and explaining fundamental interactions between frequency-squared attenuation and diffraction in the time-domain.

Analytical Methods for Causality Evaluation of Photonic Materials

We comprehensively review several general methods and analytical tools used for causality evaluation of photonic materials. Our objective is to call to mind and then formulate, on a mathematically rigorous basis, a set of theorems which can answer the question whether a considered material model is causal or not. For this purpose, a set of various distributional theorems presented in literature is collected as the distributional version of the Titchmarsh theorem, allowing for evaluation of causality in complicated electromagnetic systems. Furthermore, we correct the existing material models with the use of distribution theory in order to obtain their causal formulations. In addition to the well-known Kramers–Krönig (K–K) relations, we overview four further methods which can be used to assess causality of given dispersion relations, when calculations of integrals involved in the K–K relations are challenging or even impossible. Depending on the given problem, optimal approaches allowing us to prove either the causality or lack thereof are pointed out. These methodologies should be useful for scientists and engineers analyzing causality problems in electrodynamics and optics, particularly with regard to photonic materials, when the involved mathematical distributions have to be invoked.

Photoacoustic computed tomography for functional human brain imaging [Invited]

The successes of magnetic resonance imaging and modern optical imaging of human brain function have stimulated the development of complementary modalities that offer molecular specificity, fine spatiotemporal resolution, and sufficient penetration simultaneously. By virtue of its rich optical contrast, acoustic resolution, and imaging depth far beyond the optical transport mean free path (∼1 mm in biological tissues), photoacoustic computed tomography (PACT) offers a promising complementary modality. In this article, PACT for functional human brain imaging is reviewed in its hardware, reconstruction algorithms, in vivo demonstration, and potential roadmap.

Signal Propagation in Electromagnetic Media Modelled by the Two-Sided Fractional Derivative

In this paper, wave propagation is considered in a medium described by a fractional-order model, which is formulated with the use of the two-sided fractional derivative of Ortigueira and Machado. Although the relation of the derivative to causality is clearly specified in its definition, there is no obvious relation between causality of the derivative and causality of the transfer function induced by this derivative. Hence, causality of the system is investigated; its output is an electromagnetic signal propagating in media described by the time-domain two-sided fractional derivative. It is demonstrated that, for the derivative order in the range [1,+∞), the transfer function describing attenuated signal propagation is not causal for any value of the asymmetry parameter of the derivative. On the other hand, it is shown that, for derivative orders in the range (0,1), the transfer function is causal if and only if the asymmetry parameter is equal to certain specific values corresponding to the left-sided Grünwald–Letnikov derivative. The results are illustrated by numerical simulations and analyses. Some comments on the Kramers–Krönig relations for logarithm of the transfer function are presented as well.

Comprehensive Viscoelastic Characterization of Tissues and the Inter-relationship of Shear Wave (Group and Phase) Velocity, Attenuation and Dispersion

We report shear wave phase and group velocity, dispersion and attenuation in oil-in-gelatin viscoelastic phantoms and in vivo liver data. Moreover, we measured the power law coefficient from each dispersion curve and used it, together with the shear wave velocity, to calculate an approximate value for attenuation that agrees with independent attenuation measurements. Results in phantoms exhibit good agreement for all parameters with respect to independent mechanical measurements. For in vivo data, the livers of 20 patients were scanned. Results were compared with pathology scores obtained from liver biopsies. Across these cases, increases in shear wave dispersion and attenuation were related to increased steatosis score. It was found that shear wave dispersion and attenuation are experimentally linked, consistent with simple predictions based on the rheology of tissues, and can be used individually or jointly to assess tissue viscosity. Thus, this study indicates the possible utility of using shear wave dispersion and attenuation to non-invasively and quantitatively assess steatosis.

A Comprehensive Report on Ultrasonic Attenuation of Engineering Materials, Including Metals, Ceramics, Polymers, Fiber-Reinforced Composites, Wood, and Rocks

In this paper, ultrasonic attenuation of engineering materials is evaluated comprehensively, covering metals, ceramics, polymers, fiber-reinforced composites, wood, and rocks. After verifying two reliable experimental methods, 336 measurements are conducted and their results are tabulated. Attenuation behavior is determined over broadband spectra, extending up to 15 MHz in low attenuating materials. The attenuation spectra are characterized in combination with four power law terms, with many showing linear frequency dependence, with or without Rayleigh scattering. Dislocation damping effects are re-evaluated and a new mechanism is proposed to explain some of the linear frequency dependencies. Additionally, quadratic and cubic dependencies due to Datta–Kinra scattering and Biwa scattering, respectively, are used for some materials to construct model relations. From many test results, some previously hidden behaviors emerged upon data evaluation. Effects of cold working, tempering, and annealing are complex and sometimes contradictory. Comparison to available literature was attempted for some, but most often prior data were unavailable. This collection of new attenuation data will be of value in materials selection and in designing structural health monitoring and non-destructive inspection protocols.

Measurement of Viscoelastic Material Model Parameters Using Fractional Derivative Group Shear Wave Speeds in Simulation and Phantom Data

While ultrasound shear wave elastography originally focused on tissue stiffness under the assumption of elasticity, recent work has investigated the higher order, viscoelastic properties of the tissue. This article presents a method to use group shear wave speeds (gSWSs) at a series of derivative orders to characterize viscoelastic materials. This method uses a least squares fitting algorithm to match the experimental data to the calculated gSWS data, using an assumed material model and excitation geometry matched to the experimental imaging configuration. Building on a previous study that used particle displacement, velocity, and acceleration signals, this study extends the analysis to a continuous range of fractional derivative orders between 0 and 2. The method can be applied to any material model. Herein, material characterization was performed for three different two-parameter models and three different three-parameter models. This group speed-based method was applied to both shear wave simulations with ultrasonic tracking and experimental acquisitions in viscoelastic phantoms [similar to the Phase II Quantitative Imaging Biomarkers Alliance (QIBA) phantoms]. In both the cases, the group speed method produced more repeatable characterization overall than fitting the phase velocity results from the peak of the 2-D Fourier transform. Results suggest that the linear attenuation model is a better fit than the Voigt model for the viscoelastic QIBA phantoms.

Tractable calculation of the Green's tensor for shear wave propagation in an incompressible, transversely isotropic material

Assessing material properties from observations of shear wave propagation following an acoustic radiation force impulse (ARFI) excitation is difficult in anisotropic materials because of the complex relations among the propagation direction, shear wave polarizations, and material symmetries. In this paper, we describe a method to calculate shear wave signals using Green’s tensor methods in an incompressible, transversely isotropic (TI) material characterized by three material parameters. The Green’s tensor is written as the sum of an analytic expression for the SH propagation mode, and an integral expression for the SV propagation mode that can be evaluated by interpolation within precomputed integral functions with an efficiency comparable to the evaluation of a closed-form expression. By using parametrized integral functions, the number of requried numerical integrations is reduced by a factor of 10² − 10⁹ depending on the specific problem under consideration. Results are presented for the case of a point source positioned at the origin and a tall Gaussian source similar to an ARFI excitation. For an experimental configuration with a tilted material symmetry axis, results show that shear wave signals exhibit structures that are sufficiently complex to allow measurement of all three material parameters that characterize an incompressible, TI material.

Towards a consensus on rheological models for elastography in soft tissues

A rising wave of technologies and instruments are enabling more labs and clinics to make a variety of measurements related to tissue viscoelastic properties. These instruments include elastography imaging scanners, rheological shear viscometers, and a variety of calibrated stress-strain analyzers. From these many sources of disparate data, a common step in analyzing results is to fit the measurements of tissue response to some viscoelastic model. In the best scenario, this places the measurements within a theoretical framework and enables meaningful comparisons of the parameters against other types of tissues. However, there is a large set of established rheological models, even within the class of linear, causal, viscoelastic solid models, so which of these should be chosen? Is it simply a matter of best fit to a minimum mean squared error of the model to several data points? We argue that the long history of biomechanics, including the concept of the extended relaxation spectrum, along with data collected from viscoelastic soft tissues over an extended range of times and frequencies, and the theoretical framework of multiple relaxation models which model the multi-scale nature of physical tissues, all lead to the conclusion that fractional derivative models represent the most succinct and meaningful models of soft tissue viscoelastic behavior. These arguments are presented with the goal of clarifying some distinctions between, and consequences of, some of the most commonly used models, and with the longer term goal of reaching a consensus among different sub-fields in acoustics, biomechanics, and elastography that have common interests in comparing tissue measurements.

Exact and approximate analytical time-domain Green's functions for space-fractional wave equations

The Chen-Holm and Treeby-Cox wave equations are space-fractional partial differential equations that describe power law attenuation of the form α(ω)≈α₀|ω|^y. Both of these space-fractional wave equations are causal, but the phase velocities differ, which impacts the shapes of the time-domain Green's functions. Exact and approximate closed-form time-domain Green's functions are derived for these space-fractional wave equations, and the resulting expressions contain symmetric and maximally skewed stable probability distribution functions. Numerical results are evaluated with ultrasound parameters for breast and liver at different times as a function of space and at different distances as a function of time, where the reference calculations are computed with the Pantis method. The results show that the exact and approximate time-domain Green's functions contain both outbound and inbound propagating terms and that the inbound component is negligible a short distance from the origin. Exact and approximate analytical time-domain Green's functions are also evaluated for the Chen-Holm wave equation with power law exponent y = 1. These comparisons demonstrate that single term analytical expressions containing stable probability densities provide excellent approximations to the time-domain Green's functions for the Chen-Holm and Treeby-Cox wave equations.

Piecewise parabolic method for propagation of shear shock waves in relaxing soft solids: One-dimensional case

Shear shock waves can be generated spontaneously deep within the brain during a traumatic injury. This recently observed behavior could be a primary mechanism for the generation of traumatic brain injuries. However, shear shock wave physics and its numerical modeling are relatively unstudied. Existing numerical solvers used in biomechanics are not designed for the extremely large Mach numbers (greater than 1) observed in the brain. Furthermore, soft solids, such as the brain, have a complex nonclassical viscoleastic response, which must be accurately modeled to capture the nonlinear wave behavior. Here, we develop a 1D inviscid velocity-stress-like system to model the propagation of shear shock waves in a homogeneous medium. Then a generalized Maxwell body is used to model a relaxing medium that can describe experimentally determined attenuation laws. Finally, the resulting system is solved numerically with the piecewise parabolic method, a high-order finite volume method. The nonlinear and the relaxing components of this method are validated with theoretical predictions. Comparisons between numerical solutions obtained for the proposed model and the experiments of plane shear shock wave propagation based on high frame-rate ultrasound imaging and tracking are shown to be in excellent agreement.

Verification of the Kramers-Kronig relations between ultrasonic attenuation and phase velocity in a finite spectral range for CFRP composites

The aim of this work is to verify applicability of the Kramers-Kronig relations between the attenuation coefficient and phase velocity of longitudinal acoustic waves in carbon fiber reinforced plastic (CFRP) composites in a finite megahertz frequency range. To measure these characteristics, the method of broadband acoustic spectroscopy with a laser source of ultrasound is used. We have experimentally shown that absolute attenuation in CFRPs is determined by both absorption in a polymer matrix and scattering on carbon fibers and gas pores of several tenths of microns in size. We have also found that the increasing composite porosity leads to the decrease in the absolute value of the ultrasonic velocity and to the increase in its relative dispersion in the entire studied frequency range of 1-10 MHz. The experimental results have shown that for the studied CFRPs in this range, features of the mechanism causing attenuation and dispersion of longitudinal acoustic waves by propagation in composites do not influence on applicability of the local Kramers-Kronig relations.

Time domain model and experimental validation of non-contact surface wave scanner

This paper presents a time-domain model for the prediction of acoustic field in an air-coupled, non-contact, ultrasonic surface wave scanner, which includes an air-coupled Emitter, the Propagation space, and an air-coupled Receiver (EPR). The computation is divided into three steps, with each step being modeled in the time domain by its spatio-temporal transfer function. The latter are then used in turn, to find the pulse response of the overall system. The model takes the finite size of the aperture receiver, the attenuation in both air and the tested solid sample, as well as the electric response of the emitter-receiver set h_e into account. The attenuation is characterized by a causal time-domain Green's function, allowing wideband attenuation of a lossy medium, obeying the power law αω=α₀ω^η,0⩽η⩽2, to be used. The model is implemented numerically using a Discrete Representation approach. It is then validated quantitatively by comparing the predicted acoustic field with experiment. The prediction error for three typical field features, the system's impulse response, the on-axis field distribution, and the directivity pattern, is globally smaller than 3%. In order to obtain this high level of accuracy in the model, the parameters characterizing the solid sample used during the experiment were measured experimentally, with a specifically developed experimental setup. Overall, the proposed model is approximately 100 times faster than 3D FEM with an equivalent spatio-temporal resolution. In parallel, a simplified model is proposed, which neglects the attenuation in air and assumes the emitter inclination angle to be perfectly adjusted. This approach makes it possible to further shorten the computational time by a factor of about ten, whilst maintaining good accuracy. Thanks to its computational efficiency, the proposed model can be used to formulate various recommendations concerning the scanner settings, in particular the inclination angles of the emitter and receiver, and their distance from the sample.

Skull's Photoacoustic Attenuation and Dispersion Modeling with Deterministic Ray-Tracing: Towards Real-Time Aberration Correction

Although transcranial photoacoustic imaging has been previously investigated by several groups, there are many unknowns about the distorting effects of the skull due to the impedance mismatch between the skull and underlying layers. The current computational methods based on finite-element modeling are slow, especially in the cases where fine grids are defined for a large 3-D volume. We develop a very fast modeling/simulation framework based on deterministic ray-tracing. The framework considers a multilayer model of the medium, taking into account the frequency-dependent attenuation and dispersion effects that occur in wave reflection, refraction, and mode conversion at the skull surface. The speed of the proposed framework is evaluated. We validate the accuracy of the framework using numerical phantoms and compare its results to k-Wave simulation results. Analytical validation is also performed based on the longitudinal and shear wave transmission coefficients. We then simulated, using our method, the major skull-distorting effects including amplitude attenuation, time-domain signal broadening, and time shift, and confirmed the findings by comparing them to several ex vivo experimental results. It is expected that the proposed method speeds up modeling and quantification of skull tissue and allows the development of transcranial photoacoustic brain imaging.

Versatile Ultrasonic Spectrometer for Liquids With Practical Sample Handling by Using Standard Cuvettes

This work describes the development and tests of a new ultrasonic spectrometer for liquids based on the use of commercial cuvettes for sample handling. The spectrometer operates in the frequency range from 20 to 80 MHz and gathers some important characteristics, which are its high thermal stability (better than 0.01 °C), by the use of Peltier cells, and practical sample handling with small volume (≤3 ml) samples placed inside cuvettes which can be easily removed from the spectrometer, cleaned/sterilized, or simply discarded. Through-transmission operation is used to measure propagation velocity and attenuation coefficient, and the spectrometer was tested with mixtures of water and NaCl, which have attenuations smaller than that of distilled water, and higher attenuation samples of silicone and castor oil. Backscattering studies of polystyrene particles of 10- and 15- [Formula: see text] diameters were also conducted, showing the versatility of the instrument.

Effects of Nonlinear Propagation of Focused Ultrasound on the Stable Cavitation of a Single Bubble

Many biomedical applications such as ultrasonic targeted drug delivery, gene therapy, and molecular imaging entail the problems of manipulating microbubbles by means of a high-intensity focused ultrasound (HIFU) pressure field; namely stable cavitation. In high-intensity acoustic field, bubbles demonstrate translational instability, the well-known erratic dancing motion, which is caused by shape oscillations of the bubbles that are excited by their volume oscillations. The literature of bubble dynamics in the HIFU field is mainly centered on experiments, lacking a systematic study to determine the threshold for shape oscillations and translational motion. In this work, we extend the existing multiphysics mathematical modeling platform on bubble dynamics for taking account of (1) the liquid compressibility which allows us to apply a high-intensity acoustic field; (2) the mutual interactions of volume pulsation, shape modes, and translational motion; as well as (3) the effects of nonlinearity, diffraction, and absorption of HIFU to incorporate the acoustic nonlinearity due to wave kinematics or medium—all in one model. The effects of acoustic nonlinearity on the radial pulsations, axisymmetric modes of shape oscillations, and translational motion of a bubble, subjected to resonance and off-resonance excitation and various acoustic pressure, are examined. The results reveal the importance of considering all the involved harmonics and wave distortion in the bubble dynamics, to accurately predict the oscillations, translational trajectories, and the threshold for inertial (unstable) cavitation. This result is of interest for understanding the bubble dynamical behaviors observed experimentally in the HIFU field.

Development of a Tissue-Mimicking Phantom of the Brain for Ultrasonic Studies

Constructing tissue-mimicking phantoms of the brain for ultrasonic studies is complicated by the low backscatter coefficient of brain tissue, causing difficulties in simultaneously matching the backscatter and attenuation properties. In this work, we report on the development of a polyvinyl alcohol-based tissue-mimicking phantom with properties approaching those of human brain tissue. Polyvinyl alcohol was selected as the base material for the phantom as its properties can be varied by freeze-thaw cycling, variations in concentration and the addition of scattering inclusions, allowing some independent control of backscatter and attenuation. The ultrasonic properties (including speed of sound, attenuation and backscatter) were optimized using these methods with talc powder as an additive. It was determined that the ultrasonic properties of the phantom produced in this study are best matched to brain tissue in the frequency range 1-3 MHz, indicating its utility for laboratory ultrasonic studies in this frequency range.

The dispersion formula and the Green's function associated with an attenuation obeying a frequency power law

An attenuation obeying a frequency power law scales as |ω|β , where ω is angular frequency and β is a real constant. A recently developed dispersion formula predicts that the exponent β can take only certain values in well defined, disjoint intervals. It is shown here that these admissible values of β are consistent with the physical requirement, stemming from the second law of thermodynamics, that the work done during the passage of a wave must always be positive. Since the dispersion formula, which is derived from the strain-hardening wave equation, is a causal transform, it is expected that the associated Green's function should also satisfy causality for all the permitted values of β. Such is not the case, however: the Green's function is maximally flat at the time of source activation, and hence is causal, but only for values of β in the interval (0.5, 1). This restriction supersedes the weaker constraints on β derived from the dispersion formula alone. For the previously admissible values of β outside the interval (0.5, 1), although the dispersion formula satisfies causality, the Green's function is non-causal. Evidently, causality may be satisfied by the dispersion formula but violated by the Green's function.

Time-domain analysis of power law attenuation in space-fractional wave equations

Ultrasound attenuation in soft tissue follows a power law as a function of the ultrasound frequency, and in medical ultrasound, power law attenuation is often described by fractional calculus models that contain one or more time- or space-fractional derivatives. For certain time-fractional models, exact and approximate time-domain Green's functions are known, but similar expressions are not available for the space-fractional models that describe power law attenuation. To address this deficiency, a numerical approach for calculating time-domain Green's functions for the Chen-Holm space-fractional wave equation and Treeby-Cox space-fractional wave equation is introduced, where challenges associated with the numerical evaluation of a highly oscillatory improper integral are addressed with the Filon integration formula combined with the Pantis method. Numerical results are computed for both of these space-fractional wave equations at different distances in breast and liver with power law exponents of 1.5 and 1.139, respectively. The results show that these two space-fractional wave equations are causal and that away from the origin, the time-domain Green's function for the Treeby-Cox space-fractional wave equation is very similar to the time-domain Green's function for the time-fractional power law wave equation.

Characterization of Viscoelastic Materials Using Group Shear Wave Speeds

Recent investigations of viscoelastic properties of materials have been performed by observing shear wave propagation following localized, impulsive excitations, and Fourier decomposing the shear wave signal to parameterize the frequency-dependent phase velocity using a material model. This paper describes a new method to characterize viscoelastic materials using group shear wave speeds , , and determined from the shear wave displacement, velocity, and acceleration signals, respectively. Materials are modeled using a two-parameter linear attenuation model with phase velocity and dispersion slope at a reference frequency of 200 Hz. Analytically calculated lookup tables are used to determine the two material parameters from pairs of measured group shear wave speeds. Green's function calculations are used to validate the analytic model. Results are reported for measurements in viscoelastic and approximately elastic phantoms and demonstrate good agreement with phase velocities measured using Fourier analysis of the measured shear wave signals. The calculated lookup tables are relatively insensitive to the excitation configuration. While many commercial shear wave elasticity imaging systems report group shear wave speeds as the measures of material stiffness, this paper demonstrates that differences , , and of group speeds are first-order measures of the viscous properties of materials.

Focused ultrasound actuation of shape memory polymers; acoustic-thermoelastic modeling and testing

Controlled drug delivery (CDD) technologies have received extensive attention recently.

Fractional biharmonic operator equation model for arbitrary frequency-dependent scattering attenuation in acoustic wave propagation

This paper proposes a fractional biharmonic operator equation model in the time-space domain to describe scattering attenuation of acoustic waves in heterogeneous media. Compared with the existing models, the proposed fractional model is able to describe arbitrary frequency-dependent scattering attenuation, which typically obeys an empirical power law with an exponent ranging from 0 to 4. In stark contrast to an extensive and rapidly increasing application of the fractional derivative models for wave absorption attenuation in the literature, little has been reported on frequency-dependent scattering attenuation. This is largely because the order of the fractional Laplacian is from 0 to 2 and is infeasible for scattering attenuation. In this study, the definition of the fractional biharmonic operator in space with an order varying from 0 to 4 is proposed, as well as a fractional biharmonic operator equation model of scattering attenuation which is consistent with arbitrary frequency power-law dependency and obeys the causal relation under the smallness approximation. Finally, the correlation between the fractional order and the ratio of wavelength to the diameter of the scattering heterogeneity is investigated and an expression on exponential form is also provided.

Approximate analytical time-domain Green's functions for the Caputo fractional wave equation

The Caputo fractional wave equation [Geophys. J. R. Astron. Soc. 13, 529-539 (1967)] models power-law attenuation and dispersion for both viscoelastic and ultrasound wave propagation. The Caputo model can be derived from an underlying fractional constitutive equation and is causal. In this study, an approximate analytical time-domain Green's function is derived for the Caputo equation in three dimensions (3D) for power law exponents greater than one. The Green's function consists of a shifted and scaled maximally skewed stable distribution multiplied by a spherical spreading factor 1/(4πR). The approximate one dimensional (1D) and two dimensional (2D) Green's functions are also computed in terms of stable distributions. Finally, this Green's function is decomposed into a loss component and a diffraction component, revealing that the Caputo wave equation may be approximated by a coupled lossless wave equation and a fractional diffusion equation.

Time-domain comparisons of power law attenuation in causal and noncausal time-fractional wave equations

The attenuation of ultrasound propagating in human tissue follows a power law with respect to frequency that is modeled by several different causal and noncausal fractional partial differential equations. To demonstrate some of the similarities and differences that are observed in three related time-fractional partial differential equations, time-domain Green's functions are calculated numerically for the power law wave equation, the Szabo wave equation, and for the Caputo wave equation. These Green's functions are evaluated for water with a power law exponent of y = 2, breast with a power law exponent of y = 1.5, and liver with a power law exponent of y = 1.139. Simulation results show that the noncausal features of the numerically calculated time-domain response are only evident very close to the source and that these causal and noncausal time-domain Green's functions converge to the same result away from the source. When noncausal time-domain Green's functions are convolved with a short pulse, no evidence of noncausal behavior remains in the time-domain, which suggests that these causal and noncausal time-fractional models are equally effective for these numerical calculations.

Nonlinear acoustic pulse propagation in dispersive sediments using fractional loss operators

The nonlinear progressive wave equation (NPE) is a time-domain formulation of the Euler fluid equations designed to model low-angle wave propagation using a wave-following computational domain. The wave-following frame of reference permits the simulation of long-range propagation and is useful in modeling blast wave effects in the ocean waveguide. Existing models do not take into account frequency-dependent sediment attenuation, a feature necessary for accurately describing sound propagation over, into, and out of the ocean sediment. Sediment attenuation is addressed in this work by applying lossy operators to the governing equation that are based on a fractional Laplacian. These operators accurately describe frequency-dependent attenuation and dispersion in typical ocean sediments. However, dispersion within the sediment is found to be a secondary process to absorption and effectively negligible for ranges of interest. The resulting fractional NPE is benchmarked against a Fourier-transformed parabolic equation solution for a linear case, and against the analytical Mendousse solution to Burgers' equation for the nonlinear case. The fractional NPE is then used to investigate the effects of attenuation on shock wave propagation.

Wave-speed dispersion associated with an attenuation obeying a frequency power law

An attenuation scaling as a power of frequency, |ω|(β), over an infinite bandwidth is neither analytic nor square-integrable, thus calling into question the application of the Kramers-Krönig dispersion relations for determining the frequency dependence of the associated phase speed. In this paper, three different approaches are developed, all of which return the dispersion formula for the wavenumber, K(ω). The first analysis relies on the properties of generalized functions and the causality requirement that the impulse response, k(t), the inverse Fourier transform of -iK(ω), must vanish for t < 0. Second, a wave equation is introduced that yields the phase-speed dispersion associated with a frequency-power-law attenuation. Finally, it is shown that, with minor modification, the Kramers-Krönig dispersion relations with no subtractions (the Plemelj formulas) do in fact hold for an attenuation scaling as |ω|(β), yielding the same dispersion formula as the other two derivations. From this dispersion formula, admissible values of the exponent β are established. Physically, the inadmissible values of β, which include all the integers, correspond to attenuation-dispersion pairs whose Fourier components cannot combine in such a way as to make the impulse response, k(t), vanish for t < 0. There is no upper or lower limit on the value that β may take.

An improved wave-vector frequency-domain method for nonlinear wave modeling

In this paper, a recently developed wave-vector frequency-domain method for nonlinear wave modeling is improved and verified by numerical simulations and underwater experiments. Higher order numeric schemes are proposed that significantly increase the modeling accuracy, thereby allowing for a larger step size and shorter computation time. The improved algorithms replace the left-point Riemann sum in the original algorithm by the trapezoidal or Simpson's integration. Plane waves and a phased array were first studied to numerically validate the model. It is shown that the left-point Riemann sum, trapezoidal, and Simpson's integration have first-, second-, and third-order global accuracy, respectively. A highly focused therapeutic transducer was then used for experimental verifications. Short high-intensity pulses were generated. 2-D scans were conducted at a prefocal plane, which were later used as the input to the numerical model to predict the acoustic field at other planes. Good agreement is observed between simulations and experiments.

Estimation of fast and slow wave properties in cancellous bone using Prony's method and curve fitting

The presence of two longitudinal waves in poroelastic media is predicted by Biot's theory and has been confirmed experimentally in through-transmission measurements in cancellous bone. Estimation of attenuation coefficients and velocities of the two waves is challenging when the two waves overlap in time. The modified least squares Prony's (MLSP) method in conjuction with curve-fitting (MLSP + CF) is tested using simulations based on published values for fast and slow wave attenuation coefficients and velocities in cancellous bone from several studies in bovine femur, human femur, and human calcaneus. The search algorithm is accelerated by exploiting correlations among search parameters. The performance of the algorithm is evaluated as a function of signal-to-noise ratio (SNR). For a typical experimental SNR (40 dB), the root-mean-square errors (RMSEs) for one example (human femur) with fast and slow waves separated by approximately half of a pulse duration were 1 m/s (slow wave velocity), 4 m/s (fast wave velocity), 0.4 dB/cm MHz (slow wave attenuation slope), and 1.7 dB/cm MHz (fast wave attenuation slope). The MLSP + CF method is fast (requiring less than 2 s at SNR = 40 dB on a consumer-grade notebook computer) and is flexible with respect to the functional form of the parametric model for the transmission coefficient. The MLSP + CF method provides sufficient accuracy and precision for many applications such that experimental error is a greater limiting factor than estimation error.

Cancellous bone fast and slow waves obtained with Bayesian probability theory correlate with porosity from computed tomography

A Bayesian probability theory approach for separating overlapping ultrasonic fast and slow waves in cancellous bone has been previously introduced. The goals of this study were to investigate whether the fast and slow waves obtained from Bayesian separation of an apparently single mode signal individually correlate with porosity and to isolate the fast and slow waves from medial-lateral insonification of the calcaneus. The Bayesian technique was applied to trabecular bone data from eight human calcanei insonified in the medial-lateral direction. The phase velocity, slope of attenuation (nBUA), and amplitude were determined for both the fast and slow waves. The porosity was assessed by micro-computed tomography (microCT) and ranged from 78.7% to 94.1%. The method successfully separated the fast and slow waves from medial-lateral insonification of the calcaneus. The phase velocity for both the fast and slow wave modes showed an inverse correlation with porosity (R(2) = 0.73 and R(2) = 0.86, respectively). The slope of attenuation for both wave modes also had a negative correlation with porosity (fast wave: R(2) = 0.73, slow wave: R(2) = 0.53). The fast wave amplitude decreased with increasing porosity (R(2) = 0.66). Conversely, the slow wave amplitude modestly increased with increasing porosity (R(2) = 0.39).

Numerical investigation of the effects of shear waves in transcranial photoacoustic tomography with a planar geometry

Using a recently developed reconstruction method for photoacoustic tomography (PAT) valid for a planar measurement geometry parallel to a layered medium, we investigate the effects of shear wave propagation in the solid layer upon the ability to estimate Fourier components of the object. We examine this ability as a function of the thickness of the layer supporting shear waves as well as of the incidence angle of the field in the planewave representation. Examples are used to demonstrate the importance of accounting for shear waves in transcranial PAT. Error measures are introduced to quantify the error found when omitting shear waves from the forward model in PAT.

Determining attenuation properties of interfering fast and slow ultrasonic waves in cancellous bone

Previous studies have shown that interference between fast waves and slow waves can lead to observed negative dispersion in cancellous bone. In this study, the effects of overlapping fast and slow waves on measurements of the apparent attenuation as a function of propagation distance are investigated along with methods of analysis used to determine the attenuation properties. Two methods are applied to simulated data that were generated based on experimentally acquired signals taken from a bovine specimen. The first method uses a time-domain approach that was dictated by constraints imposed by the partial overlap of fast and slow waves. The second method uses a frequency-domain log-spectral subtraction technique on the separated fast and slow waves. Applying the time-domain analysis to the broadband data yields apparent attenuation behavior that is larger in the early stages of propagation and decreases as the wave travels deeper. In contrast, performing frequency-domain analysis on the separated fast waves and slow waves results in attenuation coefficients that are independent of propagation distance. Results suggest that features arising from the analysis of overlapping two-mode data may represent an alternate explanation for the previously reported apparent dependence on propagation distance of the attenuation coefficient of cancellous bone.

Poly(vinyl alcohol) gels as photoacoustic breast phantoms revisited

A popular phantom in photoacoustic imaging is poly(vinyl alcohol) (PVA) hydrogel fabricated by freezing and thawing (F-T) aqueous solutions of PVA. The material possesses acoustic and optical properties similar to those of tissue. Earlier work characterized PVA gels in small test specimens where temperature distributions during F-T are relatively homogeneous. In this work, in breast-sized samples we observed substantial temperature differences between the shallow regions and the interior during the F-T procedure. We investigated whether spatial variations were also present in the acoustic and optical properties. The speed of sound, acoustic attenuation, and optical reduced scattering coefficients were measured on specimens sampled at various locations in a large phantom. In general, the properties matched values quoted for breast tissue. But while acoustic properties were relatively homogeneous, the reduced scattering was substantially different at the surface compared with the interior. We correlated these variations with gel microstructure inspected using scanning electron microscopy. Interestingly, the phantom's reduced scattering spatial distribution matches the optical properties of the standard two-layer breast model used in x ray dosimetry. We conclude that large PVA samples prepared using the standard recipe make excellent breast tissue phantoms.

A k-space Green's function solution for acoustic initial value problems in homogeneous media with power law absorption

An efficient Green's function solution for acoustic initial value problems in homogeneous media with power law absorption is derived. The solution is based on the homogeneous wave equation for lossless media with two additional terms. These terms are dependent on the fractional Laplacian and separately account for power law absorption and dispersion. Given initial conditions for the pressure and its temporal derivative, the solution allows the pressure field for any time t>0 to be calculated in a single step using the Fourier transform and an exact k-space time propagator. For regularly spaced Cartesian grids, the former can be computed efficiently using the fast Fourier transform. Because no time stepping is required, the solution facilitates the efficient computation of the pressure field in one, two, or three dimensions without stability constraints. Several computational aspects of the solution are discussed, including the effect of using a truncated Fourier series to represent discrete initial conditions, the use of smoothing, and the properties of the encapsulated absorption and dispersion.

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.

Dynamic acoustoelastic testing of weakly pre-loaded unconsolidated water-saturated glass beads

Dynamic acoustoelastic testing is applied to weakly pre-loaded unconsolidated water-saturated glass beads. The gravitational acceleration produces, on the probed beads, a static stress of order 130 Pa, thus the granular medium is close to the jamming transition. A low-frequency (LF) acoustic wave gently disturbs the medium, inducing successively slight expansion and compaction of the granular packing expected to modulate the number of contacts between beads. Ultrasound (US) pulses are emitted simultaneously to dynamically detect the induced modification of the granular skeleton. US propagation velocity and attenuation both increase when the LF pressure increases. The quadratic nonlinear elastic parameter β, related to the pressure dependence of US propagation velocity, was measured in the range 60-530 if water-saturated glass beads are considered as an effective medium. A dynamic modification of US scattering induced by beads is proposed to modulate US attenuation. Complex hysteretic behaviors and tension-compression asymmetry are also observed and analyzed by time-domain and spectral analyses. Furthermore acoustic nonlinearities are measured in cases of quasi-static and dynamic variations of the LF wave amplitude, providing quantitatively similar acoustic nonlinearities but qualitatively different.

Inverse problems in cancellous bone: estimation of the ultrasonic properties of fast and slow waves using Bayesian probability theory

Quantitative ultrasonic characterization of cancellous bone can be complicated by artifacts introduced by analyzing acquired data consisting of two propagating waves (a fast wave and a slow wave) as if only one wave were present. Recovering the ultrasonic properties of overlapping fast and slow waves could therefore lead to enhancement of bone quality assessment. The current study uses Bayesian probability theory to estimate phase velocity and normalized broadband ultrasonic attenuation (nBUA) parameters in a model of fast and slow wave propagation. Calculations are carried out using Markov chain Monte Carlo with simulated annealing to approximate the marginal posterior probability densities for parameters in the model. The technique is applied to simulated data, to data acquired on two phantoms capable of generating two waves in acquired signals, and to data acquired on a human femur condyle specimen. The models are in good agreement with both the simulated and experimental data, and the values of the estimated ultrasonic parameters fall within expected ranges.

Spatial backward planar projection in absorbing media possessing an arbitrary dispersion relation

Planar projection methods have been shown to rapidly relate fields between two planes. Such an approach is particularly useful for characterizing transducers, since only a single plane needs to be measured in order to characterize an entire field. The present work considers the same approach in the presence of an arbitrary dispersion relation. Unlike traditional methods that use Fourier solutions of the time-domain wave equation, the approach starts from a frequency-domain Helmholtz equation for waves in a dispersive medium. It is shown that a transfer function similar to that derived from time domain equations can be utilized. Both the forward- and backward-projection behaviors are examined and it is demonstrated that the approach is invariant to propagation direction.

Cancellous bone analysis with modified least squares Prony's method and chirp filter: phantom experiments and simulation

The presence of two longitudinal waves in porous media is predicted by Biot's theory and has been confirmed experimentally in cancellous bone. When cancellous bone samples are interrogated in through-transmission, these two waves can overlap in time. Previously, the Modified Least-Squares Prony's (MLSP) method was validated for estimation of amplitudes, attenuation coefficients, and phase velocities of fast and slow waves, but tended to overestimate phase velocities by up to about 5%. In the present paper, a pre-processing chirp filter to mitigate the phase velocity bias is derived. The MLSP/chirp filter (MLSPCF) method was tested for decomposition of a 500 kHz-center-frequency signal containing two overlapping components: one passing through a low-density-polyethylene plate (fast wave) and another passing through a cancellous-bone-mimicking phantom material (slow wave). The chirp filter reduced phase velocity bias from 100 m/s (5.1%) to 69 m/s (3.5%) (fast wave) and from 29 m/s (1.9%) to 10 m/s (0.7%) (slow wave). Similar improvements were found for 1) measurements in polycarbonate (fast wave) and a cancellous-bone-mimicking phantom (slow wave), and 2) a simulation based on parameters mimicking bovine cancellous bone. The MLSPCF method did not offer consistent improvement in estimates of attenuation coefficient or amplitude.

Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian

The efficient simulation of wave propagation through lossy media in which the absorption follows a frequency power law has many important applications in biomedical ultrasonics. Previous wave equations which use time-domain fractional operators require the storage of the complete pressure field at previous time steps (such operators are convolution based). This makes them unsuitable for many three-dimensional problems of interest. Here, a wave equation that utilizes two lossy derivative operators based on the fractional Laplacian is derived. These operators account separately for the required power law absorption and dispersion and can be efficiently incorporated into Fourier based pseudospectral and k-space methods without the increase in memory required by their time-domain fractional counterparts. A framework for encoding the developed wave equation using three coupled first-order constitutive equations is discussed, and the model is demonstrated through several one-, two-, and three-dimensional simulations.

Simplified expressions of the subtracted Kramers-Kronig relations using the expanded forms applied to ultrasonic power-law systems

The Kramers-Kronig (KK) relations are a large class of integral transformations that exploit the broad principle of simple causality in order to link the physical properties of matter and materials. In applications to the complex-valued wavenumber for acoustic propagation, the method of subtractions is used to form convergent integral relations between the phase velocity and the attenuation coefficient. When the method of subtractions is applied in the usual manner, the integrands in the relations become unnecessarily complicated. In this work, an expanded form of the subtracted relations is presented, which is essentially a truncated Taylor series expansion of the Hilbert transforms. The implementation of the relations only requires the explicit evaluation of two simply expressed integrals involving the Hilbert transform kernel. These two integrals determine the values of the other terms in the subtracted relations, demonstrating the computational efficiency of the technique. The method is illustrated analytically through its application to power-law attenuation coefficients and its associated dispersion, which are observed in a wide variety of materials. This approach explicitly shows the central role of the Hilbert transform kernel in the KK relations, which can become obscured in other formulations.

Fractal ladder models and power law wave equations

The ultrasonic attenuation coefficient in mammalian tissue is approximated by a frequency-dependent power law for frequencies less than 100 MHz. To describe this power law behavior in soft tissue, a hierarchical fractal network model is proposed. The viscoelastic and self-similar properties of tissue are captured by a constitutive equation based on a lumped parameter infinite-ladder topology involving alternating springs and dashpots. In the low-frequency limit, this ladder network yields a stress-strain constitutive equation with a time-fractional derivative. By combining this constitutive equation with linearized conservation principles and an adiabatic equation of state, a fractional partial differential equation that describes power law attenuation is derived. The resulting attenuation coefficient is a power law with exponent ranging between 1 and 2, while the phase velocity is in agreement with the Kramers-Kronig relations. The fractal ladder model is compared to published attenuation coefficient data, thus providing equivalent lumped parameters.

Measurement of broadband temperature-dependent ultrasonic attenuation and dispersion using photoacoustics

The broadband ultrasonic characterization of biological fluids and tissues is important for the continued development and application of high-resolution ultrasound imaging modalities. Here, a photoacoustic technique for the transmission measurement of temperature-dependent ultrasonic attenuation and dispersion is described. The system uses a photoacoustic plane wave source constructed from a polymethylmethacrylate substrate with a thin optically absorbent layer. Broadband ultrasonic waves are generated by illuminating the absorbent layer with nanosecond pulses of laser light. The transmitted ultrasound waves are detected by a planar 7-microm high-finesse Fabry-Perot interferometer. Temperature-induced thickness changes in the Fabry-Perot interferometer are tracked to monitor the sample temperature and maintain the sensor sensitivity. The measured -6 dB bandwidth for the combined source and sensor is 1 to 35 MHz, with an attenuation corrected signal level at 100 MHz of -10 dB. The system is demonstrated through temperature-dependent ultrasound measurements in castor oil and olive oil. Power law attenuation parameters are extracted by fitting the experimental attenuation data to a frequency power law while simultaneously fitting the dispersion data to the corresponding Kramers-Krönig relation. The extracted parameters are compared with other calibration measurements previously reported in the literature.

Determination of power-law attenuation coefficient and dispersion spectra in multi-wall carbon nanotube composites using Kramers-Kronig relations

Using a broadband through-transmission technique, the attenuation coefficient and phase velocity spectra have been measured for a set of multi-wall carbon nanotube (MWCNT)-nylon composites (from pure nylon to 20% MWCNT by weight) in the ultrasonic frequency band from 4 to 14 MHz. The samples were found to be effectively homogeneous on spatial scales from the low end of ultrasonic wavelengths investigated and up (>0.2 mm). Using Kramers-Kronig relations, the attenuation and dispersion data were found to be consistent with a power-law attenuation model with a range of exponents from y=1.12 to y=1.19 over the measurement bandwidth. The attenuation coefficients of the respective samples are found to decrease with increasing MWCNT content and a similar trend holds also for the dispersion. In contrast, the mean phase velocities for the samples rise with increasing MWCNT content indicating an increase in the mechanical moduli.

Analytical time-domain Green's functions for power-law media

Frequency-dependent loss and dispersion are typically modeled with a power-law attenuation coefficient, where the power-law exponent ranges from 0 to 2. To facilitate analytical solution, a fractional partial differential equation is derived that exactly describes power-law attenuation and the Szabo wave equation ["Time domain wave-equations for lossy media obeying a frequency power-law," J. Acoust. Soc. Am. 96, 491-500 (1994)] is an approximation to this equation. This paper derives analytical time-domain Green's functions in power-law media for exponents in this range. To construct solutions, stable law probability distributions are utilized. For exponents equal to 0, 1/3, 1/2, 2/3, 3/2, and 2, the Green's function is expressed in terms of Dirac delta, exponential, Airy, hypergeometric, and Gaussian functions. For exponents strictly less than 1, the Green's functions are expressed as Fox functions and are causal. For exponents greater than or equal than 1, the Green's functions are expressed as Fox and Wright functions and are noncausal. However, numerical computations demonstrate that for observation points only one wavelength from the radiating source, the Green's function is effectively causal for power-law exponents greater than or equal to 1. The analytical time-domain Green's function is numerically verified against the material impulse response function, and the results demonstrate excellent agreement.

Acoustic characterization of high intensity focused ultrasound fields: a combined measurement and modeling approach

Acoustic characterization of high intensity focused ultrasound (HIFU) fields is important both for the accurate prediction of ultrasound induced bioeffects in tissues and for the development of regulatory standards for clinical HIFU devices. In this paper, a method to determine HIFU field parameters at and around the focus is proposed. Nonlinear pressure waveforms were measured and modeled in water and in a tissue-mimicking gel phantom for a 2 MHz transducer with an aperture and focal length of 4.4 cm. Measurements were performed with a fiber optic probe hydrophone at intensity levels up to 24,000 W/cm(2). The inputs to a Khokhlov-Zabolotskaya-Kuznetsov-type numerical model were determined based on experimental low amplitude beam plots. Strongly asymmetric waveforms with peak positive pressures up to 80 MPa and peak negative pressures up to 15 MPa were obtained both numerically and experimentally. Numerical simulations and experimental measurements agreed well; however, when steep shocks were present in the waveform at focal intensity levels higher than 6000 W/cm(2), lower values of the peak positive pressure were observed in the measured waveforms. This underrepresentation was attributed mainly to the limited hydrophone bandwidth of 100 MHz. It is shown that a combination of measurements and modeling is necessary to enable accurate characterization of HIFU fields.

Interference between wave modes may contribute to the apparent negative dispersion observed in cancellous bone

Previous work has shown that ultrasonic waves propagating through cancellous bone often exhibit a linear-with-frequency attenuation coefficient, but a decrease in phase velocity with frequency (negative dispersion) that is inconsistent with the causality-imposed Kramers-Kronig relations. In the current study, interfering wave modes similar to those observed in bone are shown to potentially contribute to the observed negative dispersion. Biot theory, the modified Biot-Attenborogh model, and experimental results are used to aid in simulating multiple-mode wave propagation through cancellous bone. Simulations entail constructing individual wave modes exhibiting a positive dispersion using plausible velocities and amplitudes, and then summing the individual modes to create mixed-mode output wave forms. Results of the simulations indicate that mixed-mode wave forms can exhibit negative dispersion when analyzed conventionally under the assumption that only one wave is present, even when the individual interfering waves exhibit positive dispersions in accordance with the Kramers-Kronig relations. Furthermore, negative dispersion is observed when little or no visual evidence of interference exists in the time-domain data. Understanding the mechanisms responsible for the observed negative dispersion could aid in determining the true material properties of cancellous bone, as opposed to the apparent properties measured using conventional data analysis techniques.

Ultrasonic velocity dispersion in bovine cortical bone: an experimental study

Cortical bone quality is determinant in bone fragility and its ultrasonic evaluation has become possible in clinical practice. However, the interaction between a broadband ultrasonic pulse and this complex multiscale medium remains poorly understood. The frequency dependence of phase velocity, which may impact clinical measurements, has been sparsely investigated. Our objective is to evaluate the determinants of the frequency dependence of phase velocity in bovine femoral cortical bone samples using an in vitro ultrasonic transmission device. The apparent phase velocity varies quasilinearly on a 1 MHz restricted bandwidth around 4 MHz. After compensating for diffraction effects, significant differences in velocity dispersion are obtained according to the anatomical location. The microstructure of each sample is determined using an optical microscope, which allows assessing the dependence of dispersion on the type of bone microstructure. Mostly positive but also negative values of dispersion are measured. Negative dispersion is mainly obtained in samples constituted of mixed microstructure, which may be explained by phase cancellation effects due to the presence of different microstructures within the same sample. Dispersion is shown to be related to broadband ultrasonic attenuation values, especially in the radial direction. Results are compared with the local Kramers-Kronig relationships.

Negative dispersion in bone: the role of interference in measurements of the apparent phase velocity of two temporally overlapping signals

In this study the attenuation coefficient and dispersion (frequency dependence of phase velocity) are measured using a phase sensitive (piezoelectric) receiver in a phantom in which two temporally overlapping signals are detected, analogous to the fast and slow waves typically found in measurements of cancellous bone. The phantom consisted of a flat and parallel Plexiglas plate into which a step discontinuity was milled. The phase velocity and attenuation coefficient of the plate were measured using both broadband and narrowband data and were calculated using standard magnitude and phase spectroscopy techniques. The observed frequency dependence of the phase velocity and attenuation coefficient exhibit significant changes in their frequency dependences as the interrogating ultrasonic field is translated across the step discontinuity of the plate. Negative dispersion is observed at specific spatial locations of the plate at which the attenuation coefficient rises linearly with frequency, a behavior analogous to that of bone measurements reported in the literature. For all sites investigated, broadband and narrowband data (3-7 MHz) demonstrate excellent consistency. Evidence suggests that the interference between the two signals simultaneously reaching the phase sensitive piezoelectric receiver is responsible for this negative dispersion.

Prediction of negative dispersion by a nonlocal poroelastic theory

The objective of this work is to show that the negative dispersion of ultrasonic waves propagating in cancellous bone can be explained by a nonlocal version of Biot's theory of poroelasticity. The nonlocal poroelastic formulation is presented in this work and the exact solutions for one- and two-dimensional systems are obtained by the method of Fourier transform. The nonlocal phase speeds for solid- and fluid-borne waves show the desired negative dispersion where the magnitude of dispersion is strongly dependent on the nonlocal parameters and porosity. Dependence of the phase speed and attenuation is studied for both porosity and frequency variation. It is shown that the nonlocal parameter can be easily estimated by comparing the theoretical dispersion rate with experimental observations. It is also shown that the modes of Lamb waves show similar negative dispersion when predicted by the nonlocal poroelastic theory.