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.