@article{WagnerPinskyMalhotra, author = {Wagner, Marcus and Pinsky, Peter M. and Malhotra, Manish}, title = {A multiple-frequency partial-field method for exterior acoustics based on Pad{\´e} via lanczos approximants}, series = {Proceedings of the ASME Noise Control and Acoustics Division - 2001 : presented at the 2001 ASME International Mechanical Engineering Congress and Exposition, November 11 - 16, 2001, New York, New York}, journal = {Proceedings of the ASME Noise Control and Acoustics Division - 2001 : presented at the 2001 ASME International Mechanical Engineering Congress and Exposition, November 11 - 16, 2001, New York, New York}, number = {28}, publisher = {ASME}, address = {New York, NY}, isbn = {0-7918-3551-0}, doi = {10.1115/IMECE2001/NCA-23524}, pages = {229 -- 234}, abstract = {A solution methodology is introduced for the efficient computation of the acoustic field over restricted domains and for a frequency window. Typically, such partial field solutions include, for example, surfaces enclosing the radiating structure or even single points in the computational domain. The multiple-frequency partial-field (MFPF) method starts out by reformulating the finite element matrix system into a suitable shifted form. The DtN map is used as a radiation boundary condition and is interpreted as a low rank update of the matrix problem. The shifted standard form is then approximated by a rational matrix-valued Pad{\´e} approximant and solved simultaneously over a frequency range. To obtain the Pad{\´e} approximation, a banded unsymmetric Lanczos process is applied on the standard shifted form exploiting the matrix Pad{\´e}-via-Lanczos connection. Numerical examples show the feasibility of the outlined procedure.}, language = {en} } @article{WagnerPinskyOberaietal., author = {Wagner, Marcus and Pinsky, Peter M. and Oberai, Assad and Malhotra, Manish}, title = {A Krylov subspace projection method for simultaneous solution of Helmholtz problems at multiple frequencies}, series = {Computer Methods in Applied Mechanics and Engineering}, volume = {192}, journal = {Computer Methods in Applied Mechanics and Engineering}, number = {41-42}, publisher = {Elsevier}, doi = {10.1016/S0045-7825(03)00429-8}, pages = {4609 -- 4640}, abstract = {A Krylov subspace projection method which provides simultaneous solutions of the Helmholtz equation at multiple frequencies in one solution step is presented. The projector is obtained with an unsymmetric block Lanczos algorithm applied to a transfer function derived from a finite element discretization. This approach is equivalent to a matrix-valued Pad{\´e} approximation of the transfer function. The proposed method is an extension of the formulation presented in [J. Comput. Acoust. 8 (2000) 223] to unsymmetric systems and allows the treatment of a much wider range of practical problems, including near-field and fluid-structure interaction computations}, language = {en} } @article{WagnerPinskyMalhotra, author = {Wagner, Marcus and Pinsky, Peter M. and Malhotra, Manish}, title = {Application of Pad{\´e} via Lanczos approximations for efficient multifrequency solution of Helmholtz problems}, series = {Journal of the Acoustical Society of America}, volume = {113}, journal = {Journal of the Acoustical Society of America}, number = {1}, publisher = {American Inst. of Physics}, issn = {1520-8524}, doi = {10.1121/1.1514932}, pages = {313 -- 319}, abstract = {This paper addresses the efficient solution of acoustic problems in which the primary interest is obtaining the solution only on restricted portions of the domain but over a wide range of frequencies. The exterior acoustics boundary value problem is approximated using the finite element method in combination with the Dirichlet-to-Neumann (DtN) map. The restriction domain problem is formally posed in transfer function form based on the finite element solution. In order to obtain the solution over a range of frequencies, a matrix-valued Pad{\´e} approximation of the transfer function is employed, using a two-sided block Lanczos algorithm. This approach provides a stable and efficient representation of the Pad{\´e} approximation. In order to apply the algorithm, it is necessary to reformulate the transfer function due to the frequency dependency in the nonreflecting boundary condition. This is illustrated for the case of the DtN boundary condition, but there is no restriction on the approach which can also be applied to other radiation boundary conditions. Numerical tests confirm that the approach offers significant computational speed-up.}, language = {en} }