Comparison of WKB and finite difference calculations for a two-dimensional cochlear model

Abstract

There are many points of uncertainty in the subject of cochlear models. In this paper only the question of efficient computing methods is addressed. For the cochlear model with a one-dimensional approximation for the fluid motion, Zweig, Lipes, and Pierce [J. Acoust. Soc. Am. 59, 975-982 (1976)] have shown that the WKB method agrees well with a direct numerical integration. For the two-dimensional fluid model, Neely [E.D. thesis, California Institute of Technology, Pasadena, CA (1977)] has shown that a direct finite difference solution is an order of magnitude faster than the integral equation approach used by Allen [J. Acoust. Soc. Am 61, 110-119 (1977)]. In the present work, a formal WKB solution is derived following Whitham [Linear and Nonlinear Waves (Wiley, New York, 1974)]. The advantage of this formulation is simplicity, but the disadvantage is that no error estimate is available. We find that the numerical results from the WKB solution agree well with those of Neely (1977), while the computer time is reduced by another order of magnitude. Thus, the WKB method seems to offer the satisfactory accuracy, efficiency, and flexibility for treating the more realistic cochlear models.