k-Wave Toolbox


Calculate ultrasound absorption in distilled water.


abs = waterAbsorption(f, T)


waterAbsorption calculates the ultrasonic absorption in distilled water at a given temperature and frequency using a 7th order polynomial fitted to the data given by Pinkerton (1949, The Absorption of Ultrasonic Waves in Liquids and its Relation to Molecular Constitution, Proceedings of the Physical Society. Section B, 2, 129-141). A plot of the absorption at 20degC returned by waterAbsorption along with the absorption calculated by the corresponding power law given by Szabo (2004, Diagnostic Ultrasound Imaging) is given below.

% calculate absorption between 0 and 50 MHz
f = 0:2:50;                     % [MHz]
abs = waterAbsorption(f, 20);  % [dB/cm]

% plot absorption and power law
plot(f, abs, 'k-', f, 2.17e-3 * f.^2, 'bx');
xlabel('Frequency [MHz]');
ylabel('Absorption Coefficient [dB/cm]');


f array of frequency values [MHz]
T array of water temperature values [degC]


abs absorption [dB/cm]

See Also

waterDensity, waterNonlinearity, waterSoundSpeed