Thesis (Ph. D.)--University of Wisconsin--Madison, 1973. / Vita. Typescript. eContent provider-neutral record in process. Description based on print version record. Includes bibliographical references.

<p>In this thesis, we propose and analyze a fast method for computing the solution of the Helmholtz equation in a bounded domain with a variable wave speed function. The method is based on wave splitting. The Helmholtz equation is first split into one--way wave equations which are then solved iteratively for a given tolerance. The source functions depend on the wave speed function and on the solutions of the one--way wave equations from the previous iteration. The solution of the Helmholtz equation is then approximated by the sum of the one--way solutions at every iteration. To improve the computational cost, the source functions are thresholded and in the domain where they are equal to zero, the one--way wave equations are solved with GO with a computational cost independent of the frequency. Elsewhere, the equations are fully resolved with a Runge-Kutta method. We have been able to show rigorously in one dimension that the algorithm is convergent and that for fixed accuracy, the computational cost is just O(ω^1/p) for a p-th order Runge-Kutta method. This has been confirmed by numerical experiments.</p>

This dissertation presents the Generalized Finite Element Method (GFEM) for the scalar Helmholtz equation, which describes the time harmonic acoustic wave propagation problem. We introduce several handbook functions for the Helmholtz equation, namely the planewave, wave-band, and Vekua functions, and we use these handbook functions to enrich the Finite Element space via the Partition of Unity Method to create the GFEM space. The enrichment of the approximation space by these handbook functions reduces the pollution effect due to wave number and we are able to obtain a highly accurate solution with a much smaller number of degrees-of-freedom compared with the classical Finite Element Method. The q-convergence of the handbook functions is investigated, where q is the order of the handbook function, and it is shown that asymptotically the handbook functions exhibit the same rate of exponential convergence. Hence we can conclude that the selection of the handbook functions from an admissible set should be dictated only by the ease of implementation and computational costs. Another issue addressed in this dissertation is the error coming from the artificial truncation boundary condition, which is necessary to model the Helmholtz problem set in the unbounded domain. We observe that for high q, the most significant component of the error is the one due to the artificial truncation boundary condition. Here we propose a method to assess this error by performing an additional computation on the extended domain using GFEM with high q.

No description available.

Thesis (Ph. D.)--University of Hong Kong, 2010. / Includes bibliographical references (leaves 225-235). Also available in print.

In this thesis, we propose and analyze a fast method for computing the solution of the Helmholtz equation in a bounded domain with a variable wave speed function. The method is based on wave splitting. The Helmholtz equation is first split into one--way wave equations which are then solved iteratively for a given tolerance. The source functions depend on the wave speed function and on the solutions of the one--way wave equations from the previous iteration. The solution of the Helmholtz equation is then approximated by the sum of the one--way solutions at every iteration. To improve the computational cost, the source functions are thresholded and in the domain where they are equal to zero, the one--way wave equations are solved with GO with a computational cost independent of the frequency. Elsewhere, the equations are fully resolved with a Runge-Kutta method. We have been able to show rigorously in one dimension that the algorithm is convergent and that for fixed accuracy, the computational cost is just O(ω1/p) for a p-th order Runge-Kutta method. This has been confirmed by numerical experiments.

published_or_final_version / Mechanical Engineering / Doctoral / Doctor of Philosophy

No description available.

"In this dissertation we study the so-called “city effect” problem. This effect occurs when earthquakes strike large cities. In earlier studies, seismic wave propagation was evaluated in a separate step and then impacts on man made structures above ground were calculated. The 1985 Michoacan earthquake in Mexico City led Wirgin and Bard (1996) to hypothesize that city buildings may collectively affect the ground motion during an earthquake. Ghergu and Ionescu (2009) proposed a model of this phenomenon and a solution algorithm. Our contribution is to extend their work and to provide a mathematical analysis for proving the existence of preferred frequencies coupling vibrations of buildings to underground seismic waves. Given the geometry and the specific physical constants of an idealized two dimensional city, Ghergu and Ionescu computed a frequency that will couple vibrating buildings to underground seismic waves. This frequency was obtained by increasing the number of buildings at the expense of solving larger and larger systems. Our idea is to use a periodic Green's function and perform computations on a single period. That allows for much faster computations, and makes it possible to consider more complex geometries within a single period. We provide a rather in depth and proof based account of different formulations for the periodic Green's function that we need. We show that they are indeed fundamental solutions to the Helmholtz operator and we analyze their convergence rate. Finally, we give a mathematical proof of existence of preferred frequencies coupling vibrations of buildings to underground seismic waves."

Thesis (M.S. in Engineering Acoustics)--Naval Postgraduate School, December 2003. / Thesis advisor(s): S.E. Forsythe, S.R. Baker. Includes bibliographical references (p. 167). Also available online.