.php xmlns="http://www.w3.org/1999/.php"> - k-Wave MATLAB Toolbox
k-Wave Toolbox


Compute axial pressure for focused annulus transducer using O'Neil's solution.


p_axial = focusedAnnulusONeil(radius, diameter, amplitude, phase, frequency, sound_speed, density, axial_position)


focusedAnnulusONeil calculates the axial pressure for a focused annulus transducer using O'Neil's solution (O'Neil, H. Theory of focusing radiators. J. Acoust. Soc. Am., 21(5), 516-526, 1949). The annuluar elements are uniformly driven by a continuous wave sinusoid at a given frequency and normal surface velocity.

The solution is evaluated at the positions along the beam axis given by axial_position (where 0 corresponds to the transducer surface).

Note, O'Neil's formulae are derived under the assumptions of the Rayleigh integral, which are valid when the transducer diameter is large compared to both the transducer height and the acoustic wavelength.


radius transducer radius of curvature [m]
diameter 2 x num_elements array containing pairs of inner and outer aperture diameter (diameter of opening) [m]
amplitude array containing the normal surface velocities for each element [m/s]
phase array containing the phase for each element [rad]
frequency driving frequency [Hz]
sound_speed speed of sound in the propagating medium [m/s]
density density in the propagating medium [kg/m^3]
axial_position vector of positions along the beam axis where the pressure amplitude is calculated [m]


p_axial pressure amplitude at the positions specified by axial_position [Pa]

See Also

focusedBowlONeil, mendousse