Hydraulic shock (colloquial: water hammer; fluid hammer) is a pressure surge or wave caused when a fluid in motion is forced to stop or change direction suddenly: a momentum change. It is usually observed in a liquid but gases can also be affected. This phenomenon commonly occurs when a valve closes suddenly at an end of a pipeline system and a pressure wave propagates in the pipe.

Effect of a pressure surge on a float gauge

This pressure wave can cause major problems, from noise and vibration to pipe rupture or collapse. It is possible to reduce the effects of the water hammer pulses with accumulators, expansion tanks, surge tanks, blowoff valves, and other features. The effects can be avoided by ensuring that no valves will close too quickly with significant flow, but there are many situations that can cause the effect.

Rough calculations can be made using the Zhukovsky (Joukowsky) equation,[1] or more accurate ones using the method of characteristics.

History

edit

In the 1st century B.C., Marcus Vitruvius Pollio described the effect of water hammer in lead pipes and stone tubes of the Roman public water supply.[2][3] Water hammer was exploited before there was even a word for it.

The Alhambra, built by Nasrid Sultan Ibn al-Ahmar of Granada beginning 1238, used a hydram to raise water. Through a first reservoir, filled by a channel from the Darro River, water emptied via a large vertical channel into a second reservoir beneath, creating a whirlpool that in turn propelled water through a much smaller pipe up six metres whilst most water drained into a second, slightly larger pipe.[4]

In 1772, Englishman John Whitehurst built a hydraulic ram for a home in Cheshire, England.[5] In 1796, French inventor Joseph Michel Montgolfier (1740–1810) built a hydraulic ram for his paper mill in Voiron.[6] In French and Italian, the terms for "water hammer" come from the hydraulic ram: coup de bélier (French) and colpo d'ariete (Italian) both mean "blow of the ram".[7] As the 19th century witnessed the installation of municipal water supplies, water hammer became a concern to civil engineers.[8][9][10] Water hammer also interested physiologists who were studying the circulatory system.[11]

Although it was prefigured in work by Thomas Young,[12][11] the theory of water hammer is generally considered to have begun in 1883 with the work of German physiologist Johannes von Kries (1853–1928), who was investigating the pulse in blood vessels.[13][14] However, his findings went unnoticed by civil engineers.[15][16] Kries's findings were subsequently derived independently in 1898 by the Russian fluid dynamicist Nikolay Yegorovich Zhukovsky (1847–1921),[1][17] in 1898 by the American civil engineer Joseph Palmer Frizell (1832–1910),[18][19] and in 1902 by the Italian engineer Lorenzo Allievi (1856–1941).[20]

Cause and effect

edit

Water flowing through a pipe has momentum. If the moving water is suddenly stopped, such as by closing a valve downstream of the flowing water, the pressure can rise suddenly with a resulting shock wave. In domestic plumbing this shock wave is experienced as a loud banging resembling a hammering noise. Water hammer can cause pipelines to break if the pressure is sufficiently high. Air traps or stand pipes (open at the top) are sometimes added as dampers to water systems to absorb the potentially damaging forces caused by the moving water.

For example, the water traveling along a tunnel or pipeline to a turbine in a hydroelectric generating station may be slowed suddenly if a valve in the path is closed too quickly. If there is 14 km (8.7 mi) of tunnel of 7.7 m (25 ft) diameter full of water travelling at 3.75 m/s (8.4 mph),[21] that represents approximately 8,000 megajoules (2,200 kWh) of kinetic energy. This energy can be dissipated by a vertical surge shaft into which the water flows[22] which is open at the top. As the water rises up the shaft its kinetic energy is converted into potential energy, avoiding sudden high pressure. At some hydroelectric power stations, such as the Saxon Falls Hydro Power Plant In Michigan, what looks like a water tower is in fact a surge drum.[23]

In residential plumbing systems, water hammer may occur when a dishwasher, washing machine or toilet suddenly shuts off water flow. The result may be heard as a loud bang, repetitive banging (as the shock wave travels back and forth in the plumbing system), or as some shuddering.

Other potential causes of water hammer:

  • A pump stopping
  • A check valve which closes quickly (i.e., "check valve slam") due to the flow in a pipe reversing direction on loss of motive power, such as a pump stopping. "Non-slam" check valves can be used to reduce the pressure surge.
  • Filling an empty pipe that has a restriction such as a partially open valve or an orifice that allows air to pass easily as the pipe rapidly fills, but with the pressure increasing once full the water encounters the restriction.
edit
 
Expansion joints on a steam line that have been destroyed by steam hammer

Steam hammer can occur in steam systems when some of the steam condenses into water in a horizontal section of the piping. The steam forcing the liquid water along the pipe forms a "slug" which impacts a valve of pipe fitting, creating a loud hammering noise and high pressure. Vacuum caused by condensation from thermal shock can also cause a steam hammer. Steam hammer or steam condensation induced water hammer (CIWH) was exhaustively investigated both experimentally and theoretically more than a decade ago because it can have radical negative effects in nuclear power plants.[24] It is possible to theoretically explain the 2 millisecond duration 130 bar overpressure peaks with a special 6 equation multiphase thermohydraulic model,[25] similar to RELAP.

Steam hammer can be minimized by using sloped pipes and installing steam traps.

On turbocharged internal combustion engines, a "gas hammer" can take place when the throttle is closed while the turbocharger is forcing air into the engine. There is no shockwave but the pressure can still rapidly increase to damaging levels or cause compressor surge. A pressure relief valve placed before the throttle prevents the air from surging against the throttle body by diverting it elsewhere, thus protecting the turbocharger from pressure damage. This valve can either recirculate the air into the turbocharger's intake (recirculation valve), or it can blow the air into the atmosphere and produce the distinctive hiss-flutter of an aftermarket turbocharger (blowoff valve).

Mitigation measures

edit

Water hammers have caused accidents and fatalities, but usually damage is limited to breakage of pipes or appendages. An engineer should always assess the risk of a pipeline burst. Pipelines transporting hazardous liquids or gases warrant special care in design, construction, and operation. Hydroelectric power plants especially must be carefully designed and maintained because the water hammer can cause water pipes to fail catastrophically.

The following characteristics may reduce or eliminate water hammer:

  • Reduce the pressure of the water supply to the building by fitting a regulator.
  • Lower fluid velocities. To keep water hammer low, pipe-sizing charts for some applications recommend flow velocity at or below 1.5 m/s (4.9 ft/s).
  • Fit slowly closing valves. Toilet fill valves are available in a quiet fill type that closes quietly.
  • Non-slam check valves do not rely on fluid flow to close and will do so before the water flow reaches significant velocity.
  • High pipeline pressure rating (does not reduce the effect but protects against damage).
  • Good pipeline control (start-up and shut-down procedures).
  • Water towers (used in many drinking water systems) or surge tanks help maintain steady flow rates and trap large pressure fluctuations.
  • Air vessels such as expansion tanks and some types of hydraulic accumulators work in much the same way as water towers, but are pressurized. They typically have an air cushion above the fluid level in the vessel, which may be regulated or separated by a bladder. Sizes of air vessels may be up to hundreds of cubic meters on large pipelines. They come in many shapes, sizes and configurations. Such vessels often are called accumulators or expansion tanks.
  • A hydropneumatic device similar in principle to a shock absorber called a 'Water Hammer Arrestor' can be installed between the water pipe and the machine, to absorb the shock and stop the banging.
  • Air valves often remediate low pressures at high points in the pipeline. Though effective, sometimes large numbers of air valves need be installed. These valves also allow air into the system, which is often unwanted. Blowoff valves may be used as an alternative.
  • Shorter branch pipe lengths.
  • Shorter lengths of straight pipe, i.e. add elbows, expansion loops. Water hammer is related to the speed of sound in the fluid, and elbows reduce the influences of pressure waves.
  • Arranging the larger piping in loops that supply shorter smaller run-out pipe branches. With looped piping, lower velocity flows from both sides of a loop can serve a branch.
  • Flywheel on a pump.
  • Pumping station bypass.

Magnitude of the pulse

edit
 
Typical pressure wave caused by closing a valve in a pipeline

One of the first to successfully investigate the water hammer problem was the Italian engineer Lorenzo Allievi.

Water hammer can be analyzed by two different approaches—rigid column theory, which ignores compressibility of the fluid and elasticity of the walls of the pipe, or by a full analysis that includes elasticity. When the time it takes a valve to close is long compared to the propagation time for a pressure wave to travel the length of the pipe, then rigid column theory is appropriate; otherwise considering elasticity may be necessary.[26] Below are two approximations for the peak pressure, one that considers elasticity, but assumes the valve closes instantaneously, and a second that neglects elasticity but includes a finite time for the valve to close.

Instant valve closure; compressible fluid

edit

The pressure profile of the water hammer pulse can be calculated from the Joukowsky equation[27]

 

So for a valve closing instantaneously, the maximal magnitude of the water hammer pulse is

 

where ΔP is the magnitude of the pressure wave (Pa), ρ is the density of the fluid (kg/m3), a0 is the speed of sound in the fluid (m/s), and Δv is the change in the fluid's velocity (m/s). The pulse comes about due to Newton's laws of motion and the continuity equation applied to the deceleration of a fluid element.[28]

Equation for wave speed

edit

As the speed of sound in a fluid is  , the peak pressure depends on the fluid compressibility if the valve is closed abruptly.

 

where

a = wave speed,
B = equivalent bulk modulus of elasticity of the system fluid–pipe,
ρ = density of the fluid,
K = bulk modulus of elasticity of the fluid,
E = elastic modulus of the pipe,
D = internal pipe diameter,
t = pipe wall thickness,
c = dimensionless parameter due to system pipe-constraint condition[clarify] on wave speed.[28][page needed]

Slow valve closure; incompressible fluid

edit

When the valve is closed slowly compared to the transit time for a pressure wave to travel the length of the pipe, the elasticity can be neglected, and the phenomenon can be described in terms of inertance or rigid column theory:

 

Assuming constant deceleration of the water column (dv/dt = v/t), this gives

 

where:

F = force [N],
m = mass of the fluid column [kg],
a = acceleration [m/s2],
P = pressure [Pa],
A = pipe cross-section [m2],
ρ = fluid density [kg/m3],
L = pipe length [m],
v = flow velocity [m/s],
t = valve closure time [s].

The above formula becomes, for water and with imperial unit,

 

For practical application, a safety factor of about 5 is recommended:

 

where P1 is the inlet pressure in psi, V is the flow velocity in ft/s, t is the valve closing time in seconds, and L is the upstream pipe length in feet.[29]

Hence, we can say that the magnitude of the water hammer largely depends upon the time of closure, elastic components of pipe & fluid properties.[30]

Expression for the excess pressure due to water hammer

edit

When a valve with a volumetric flow rate Q is closed, an excess pressure ΔP is created upstream of the valve, whose value is given by the Joukowsky equation:

 

In this expression:[31]

ΔP is the overpressurization in Pa;
Q is the volumetric flow in m3/s;
Z is the hydraulic impedance, expressed in kg/m4/s.

The hydraulic impedance Z of the pipeline determines the magnitude of the water hammer pulse. It is itself defined by

 

where

ρ the density of the liquid, expressed in kg/m3;
A cross sectional area of the pipe, m2;
B equivalent modulus of compressibility of the liquid in the pipe, expressed in Pa.

The latter follows from a series of hydraulic concepts:

  • compressibility of the liquid, defined by its adiabatic compressibility modulus Bl, resulting from the equation of state of the liquid generally available from thermodynamic tables;
  • the elasticity of the walls of the pipe, which defines an equivalent bulk modulus of compressibility for the solid Bs. In the case of a pipe of circular cross-section whose thickness t is small compared to the diameter D, the equivalent modulus of compressibility is given by the formula  , in which E is the Young's modulus (in Pa) of the material of the pipe;
  • possibly compressibility Bg of gas dissolved in the liquid, defined by  
    γ being the specific heat ratio of the gas,
    α the rate of ventilation (the volume fraction of undissolved gas),
    and P the pressure (in Pa).

Thus, the equivalent elasticity is the sum of the original elasticities:

 

As a result, we see that we can reduce the water hammer by:

  • increasing the pipe diameter at constant flow, which reduces the flow velocity and hence the deceleration of the liquid column;
  • employing the solid material as tight as possible with respect to the internal fluid bulk (solid Young modulus low with respect to fluid bulk modulus);
  • introducing a device that increases the flexibility of the entire hydraulic system, such as a hydraulic accumulator;
  • where possible, increasing the fraction of undissolved gases in the liquid.

Dynamic equations

edit

The water hammer effect can be simulated by solving the following partial differential equations.

 
 

where V is the fluid velocity inside pipe,   is the fluid density, B is the equivalent bulk modulus, and f is the Darcy–Weisbach friction factor.[32]

Column separation

edit

Column separation is a phenomenon that can occur during a water-hammer event. If the pressure in a pipeline drops below the vapor pressure of the liquid, cavitation will occur (some of the liquid vaporizes, forming a bubble in the pipeline, keeping the pressure close to the vapor pressure). This is most likely to occur at specific locations such as closed ends, high points or knees (changes in pipe slope). When subcooled liquid flows into the space previously occupied by vapor the area of contact between the vapor and the liquid increases. This causes the vapor to condense into the liquid reducing the pressure in the vapor space. The liquid on either side of the vapor space is then accelerated into this space by the pressure difference. The collision of the two columns of liquid (or of one liquid column if at a closed end) causes a large and nearly instantaneous rise in pressure. This pressure rise can damage hydraulic machinery, individual pipes and supporting structures. Many repetitions of cavity formation and collapse may occur in a single water-hammer event.[33]

Simulation software

edit

Most water hammer software packages use the method of characteristics[28] to solve the differential equations involved. This method works well if the wave speed does not vary in time due to either air or gas entrainment in a pipeline. The wave method (WM) is also used in various software packages. WM lets operators analyze large networks efficiently. Many commercial and non-commercial packages are available.

Software packages vary in complexity, dependent on the processes modeled. The more sophisticated packages may have any of the following features:

  • Multiphase flow capabilities.
  • An algorithm for cavitation growth and collapse.
  • Unsteady friction: the pressure waves dampen as turbulence is generated and due to variations in the flow velocity distribution.
  • Varying bulk modulus for higher pressures (water becomes less compressible).
  • Fluid structure interaction: the pipeline reacts on the varying pressures and causes pressure waves itself.

Applications

edit
  • The water hammer principle can be used to create a simple water pump called a hydraulic ram.
  • Leaks can sometimes be detected using water hammer.
  • Enclosed air pockets can be detected in pipelines.
  • The water hammer from a liquid jet created by a collapsing microcavity is studied for potential applications noninvasive transdermal drug delivery.[34]

See also

edit

References

edit
  1. ^ a b Joukowsky, Nikolay (1900), "Über den hydraulischen Stoss in Wasserleitungsröhren" [On hydraulic shock in water pipes], Mémoires de l'Académie Impériale des Sciences de St.-Pétersbourg, 8th series (in German), 9 (5): 1–71
  2. ^ Vitruvius Pollio with Morris Hicky Morgan, trans. The Ten Books on Architecture (Cambridge, Massachusetts: Harvard University Press, 1914) ; Book 8, Chapter 6, sections 5-8 , pp. 245-246. Archived 2012-07-11 at the Wayback Machine Vitruvius states that when a water pipe crosses a wide valley, it must sometimes be constructed as an inverted siphon. He states that cavities ("venters") must be constructed periodically along the pipe "and in the venter, water cushions must be constructed to relieve the pressure of the air." "But if there is no such venter made in the valleys, nor any substructure built on a level, but merely an elbow, the water will break out, and burst the joints of the pipes." Swiss engineer Martin Schwarz — Martin Schwarz, "Neue Forschungsergebnisse zu Vitruvs colliviaria" [New research results on Vitruvius' colliviaria], pp. 353-357, in: Christoph Ohlig, ed., Cura Aquarum in Jordanien (Siegburg, Germany: Deutschen Wasserhistorischen Gesellschaft, 2008) — argues that Vitruvius' phrase vis spiritus referred not to air pressure, but to pressure transients (water hammer) in the water pipes. He found stone plugs (colliviaria) in Roman water pipes, which could be expelled by water hammer, allowing water in the pipe to flood the air chamber above the pipe, instead of rupturing the pipe.
  3. ^ Ismaier, Andreas (2011), Untersuchung der fluiddynamischen Wechselwirkung zwischen Druckstößen und Anlagenkomponenten in Kreiselpumpensystemen [Investigation of the fluid dynamic interaction between pressure surges and system components in centrifugal pumping systems], Schriftenreihe des Lehrstuhls für Prozessmaschinen und Anlagentechnik, Universität Erlangen; Nürnberg Lehrstuhl für Prozessmaschinen und Anlagentechnik (in German), vol. 11, Shaker, ISBN 978-3-8322-9779-4
  4. ^ The hidden world beneath the ancient Alhambra fortress. BBC 2020. Film Grenada, BBC and youtube
  5. ^ Whitehurst, John (1775), "Account of a machine for raising water, executed at Oulton, in Cheshire, in 1772", Philosophical Transactions of the Royal Society of London, 65: 277–279, doi:10.1098/rstl.1775.0026, archived from the original on 2017-03-28 See also plate preceding page 277.
  6. ^ Montgolfier, J. M. de (1803), "Note sur le bélier hydraulique, et sur la manière d'en calculer les effets" [Note on the hydraulic ram, and on the method of calculating its effects] (PDF), Journal des Mines (in French), 13 (73): 42–51, archived (PDF) from the original on 2013-10-18
  7. ^ Tijsseling, A. S.; Anderson, A. (2008), "Thomas Young's research on fluid transients: 200 years on" (PDF), Proceedings of the 10th International Conference on Pressure Surges, Edinburgh, UK: 21–33, archived (PDF) from the original on 2013-10-24 see page 22.
  8. ^ Ménabréa*, L. F. (1858), "Note sur les effects de choc de l'eau dans les conduites" [Note on the effects of water shocks in pipes], Comptes rendus (in French), 47: 221–224, archived from the original on 2017-03-28 *Luigi Federico Menabrea (1809–1896), Italian general, statesman and mathematician.
  9. ^ Michaud*, J. (1878), "Coups de bélier dans les conduites. Étude des moyens employés pour en atténeur les effects" [Water hammer in pipes. Study of means used to mitigate its effects], Bulletin de la Société Vaudoise des Ingénieurs et des Architectes (in French), 4 (3, 4): 56–64, 65–77 Available at: E.T.H. (Eidgenössische Technische Hochschule, Federal Institute of Technology) (Zürich, Switzerland). *Jules Michaud (1848–1920), Swiss engineer.
  10. ^ Castigliano, Alberto (1874). "Intorno alla resistenza dei tubi alle pressioni continue e ai colpi d'ariete" [Regarding the resistance of pipes to continuous pressures and water hammer]. Atti della Reale accademia della scienze di Torino [Proceedings of the Royal Academy of Sciences of Turin] (in Italian). 9: 222–252. *Carlo Alberto Castigliano (1847–1884), Italian mathematician and physicist.
  11. ^ a b Tijsseling, A. S.; Anderson, A. (2008). Hunt, S. (ed.). "Thomas Young's research on fluid transients: 200 years on". Proc. of the 10th Int. Conf. On Pressure Surges. Edinburgh, United Kingdom: BHR Group: 21–33. ISBN 978-1-85598-095-2.
  12. ^ Young, Thomas (1808). "Hydraulic investigations, subservient to an intended Croonian lecture on the motion of the blood". Philosophical Transactions of the Royal Society of London. 98: 164–186.
  13. ^ von Kries, J. (1883), "Ueber die Beziehungen zwischen Druck und Geschwindigkeit, welche bei der Wellenbewegung in elastischen Schläuchen bestehen" [On the relationships between pressure and velocity, which exist in connection with wave motion in elastic tubing], Festschrift der 56. Versammlung Deutscher Naturforscher und Ärzte (Festschrift of the 56th Convention of German Scientists and Physicians) (in German), Tübingen, Germany: Akademische Verlagsbuchhandlung: 67–88, archived from the original on 2017-03-28
  14. ^ von Kries, J. (1892), Studien zur Pulslehre [Studies in Pulse Science] (in German), Tübingen, Germany: Akademische Verlagsbuchhandlung, archived from the original on 2017-03-28
  15. ^ Tijsseling, Arris S.; Anderson, Alexander (2004), "A precursor in waterhammer analysis – rediscovering Johannes von Kries" (PDF), Proceedings of the 9th International Conference on Pressure Surges, Chester, UK: 739–751, archived (PDF) from the original on 2016-03-04
  16. ^ Tijsseling, Arris S.; Anderson, Alexander (2007), "Johannes von Kries and the history of water hammer", Journal of Hydraulic Engineering, 133 (1): 1–8, doi:10.1061/(ASCE)0733-9429(2007)133:1(1)
  17. ^ Tijsseling, Arris S.; Anderson, Alexander (2006), The Joukowsky equation for fluids and solids (PDF), archived (PDF) from the original on 2012-09-12
  18. ^ Frizell, J. P. (1898), "Pressures resulting from changes of velocity of water in pipes", Transactions of the American Society of Civil Engineers, 39: 1–18, doi:10.1061/TACEAT.0001315, archived from the original on 2017-03-28
  19. ^ Hale, R. A. (September 1911), "Obituary: Joseph Palmer Frizell, M. Am. Soc. C. E.", Transactions of the American Society of Civil Engineers, 73: 501–503, archived from the original on 2017-03-29
  20. ^ See:
    • Allievi, L. (1902), "Teoria generale del moto perturbato dell'acqua nei tubi in pressione (colpo d'ariete)" [General theory of the perturbed motion of water in pipes under pressure (water hammer)], Annali della Società degli Ingegneri ed Architetti Italiani (Annals of the Society of Italian Engineers and Architects) (in Italian), 17 (5): 285–325
    • Reprinted: Allievi, L. (1903). "Teoria generale del moto perturbato dell'acqua nei tubi in pressione (colpo d'ariete)". Atti dell'Associazione elettrotecnica italiana [Proceedings of the Italian Electrotechnical Association] (in Italian). 7 (2–3): 140–196.
  21. ^ "Surge shaft in 14km water conductor system in hyro project - Hydraulics and Hydrology Forum - Hydraulics and Hydrology - be Communities by Bentley". communities.bentley.com. Archived from the original on 18 January 2013. Retrieved 3 February 2022.
  22. ^ "CR4 - Thread: Pressure Shaft and Surge Shaft". Archived from the original on 2011-12-20. Retrieved 2012-07-16.
  23. ^ "Saxon Falls Hydro Generating Station | Xcel Energy". www.xcelenergy.com. Archived from the original on 2017-08-16. Retrieved 2017-08-16.
  24. ^ Barna, I.F.; Imre, A.R.; Baranyai, G.; Ézsöl, Gy. (2010). "Experimental and theoretical study of steam condensation induced water hammer phenomena". Nuclear Engineering and Design. 240 (1): 146–150. Bibcode:2010NuEnD.240..146B. doi:10.1016/j.nucengdes.2009.09.027 – via Elsevier.
  25. ^ Tiselj, I.; Petelin, S. (1997). "Modelling of Two-Phase Flow with Second-Order Accurate Scheme". Journal of Computational Physics. 136 (2): 503–521. Bibcode:1997JCoPh.136..503T. doi:10.1006/jcph.1997.5778 – via Elsevier.
  26. ^ Bruce, S.; Larock, E.; Jeppson, R. W.; Watters, G. Z. (2000), Hydraulics of Pipeline Systems, CRC Press, ISBN 0-8493-1806-8
  27. ^ Thorley, A. R. D. (2004), Fluid Transients in Pipelines (2nd ed.), Professional Engineering Publishing, ISBN 0-79180210-8[page needed]
  28. ^ a b c Streeter, V. L.; Wylie, E. B. (1998), Fluid Mechanics (International 9th Revised ed.), McGraw-Hill Higher Education, ISBN 0070625379[page needed]
  29. ^ "Water Hammer & Pulsation". Archived 2008-07-01 at the Wayback Machine
  30. ^ "What is Water Hammer/Steam Hammer ?". www.forbesmarshall.com. Retrieved 2019-12-26.
  31. ^ Faisandier, J., Hydraulic and Pneumatic Mechanisms, 8th edition, Dunod, Paris, 1999, ISBN 2100499483.
  32. ^ Chaudhry, Hanif (1979). Applied Hydraulic Transients. New York: Van Nostrand Reinhold.
  33. ^ Bergeron, L., 1950. Du Coup de Bélier en Hydraulique - Au Coup de Foudre en Electricité. (Waterhammer in hydraulics and wave surges in electricity.) Paris: Dunod (in French). (English translation by ASME Committee, New York: John Wiley & Sons, 1961.)
  34. ^ Postema M, van Wamel A, Lancée CT, de Jong N (2004). "Ultrasound-induced encapsulated microbubble phenomena". Ultrasound in Medicine & Biology. 30 (6): 827–840. doi:10.1016/j.ultrasmedbio.2004.02.010. PMID 15219962. S2CID 33442395.
edit