On the Convergence Behavior of IDR(s) and Related Methods

An explanation is given of the convergence behavior of IDR(s) methods. The convergence mechanism of these algorithms has two components. The first consists of damping properties of certain factors in the residual polynomials, which becomes less important for large values of s. The second component depends on the behavior of Lanczos polynomials...

Comparison of fundamental, second harmonic, and superharmonic imaging: A simulation study

In medical ultrasound, fundamental imaging (FI) uses the reflected echoes from the same spectral band as that of the emitted pulse. The transmission frequency determines the trade-off between penetration depth and spatial resolution. Tissue harmonic imaging (THI) employs the second harmonic of the emitted frequency band to construct images....

An iterative method for the computation of nonlinear, wide-angle, pulsed acoustic fields of medical diagnostic transducers

The development and optimization of medical ultrasound transducers and imaging modalities require a computational method that accurately predicts the nonlinear acoustic pressure field. A prospective method should provide the wide-angle, pulsed field emitted by an arbitrary planar source distribution and propagating in a three-dimensional, large...

A filtered convolution method for the computation of acoustic wave fields in very large spatiotemporal domains

The full-wave computation of transient acoustic fields with sizes in the order of 100x100x100 wavelengths by 100 periods requires a numerical method that is extremely efficient in terms of storage and computation. Iterative integral equation methods offer a good performance on these points, provided that the recurring spatiotemporal convolutions...

A new iterative solver for the time-harmonic wave equation

The time-harmonic wave equation, also known as the Helmholtz equation, is obtained if the constant-density acoustic wave equation is transformed from the time domain to the frequency domain. Its discretization results in a large, sparse, linear system of equations. In two dimensions, this system can be solved efficiently by a direct method. In...