Tl. Szabo, TIME-DOMAIN WAVE-EQUATIONS FOR LOSSY MEDIA OBEYING A FREQUENCY POWER-LAW, The Journal of the Acoustical Society of America, 96(1), 1994, pp. 491-500
For attenuation described by a slowly varying power law function of fr
equency, alpha=alpha(0)\omega\(y), classical lossy time domain wave eq
uations exist only for the restricted cases where y=0 or y=2. For the
frequently occurring practical situation in which attenuation is much
smaller than the wave number, a lossy dispersion characteristic is der
ived that has the desired attenuation general power law dependence. In
order to obtain the corresponding time domain lossy wave equation, ti
me domain loss operators similar in function to existing derivative op
erators are developed through the use of generalized functions. Three
forms of lossy wave equations are found, depending on whether y is an
even or odd integer or a noninteger. A time domain expression of causa
lity analogous in function to the Kramers-Kronig relations in the freq
uency domain is used to derive the causal wave equations. Final causal
versions of the time domain wave equations are obtained even for the
cases where y greater than or equal to 1, which, according to the Pale
y-Wiener theorem, are unobtainable from the Kramers-Kronig relations.
Different forms of the wave equation are derived including normal time
, retarded time, and parabolic (one and three dimensional). These equa
tions compare favorably with those from the literature corresponding t
o y=0, 0.5, 1, and 2.