Multi-scale wave propagation problems are computationally costly to solve by
traditional techniques because the smallest scales must be represented over a
domain determined by the largest scales of the problem. We have developed and
analyzed new numerical methods for multi-scale wave propagation in the
framework of heterogeneous multi-scale method. The numerical methods couples
simulations on macro- and micro-scales for problems with rapidly oscillating
coefficients.
We expose in full detail a constructive procedure to invert the so--called
"finite Markov moment problem". The proofs rely on the general theory of
Toeplitz matrices together with the classical Newton's relations.
We consider a class of finite Markov moment problems with arbitrary number of
positive and negative branches. We show criteria for the existence and
uniqueness of solutions, and we characterize in detail the non-unique solution
families. Moreover, we present a constructive algorithm to solve the moment
problems numerically and prove that the algorithm computes the right solution.
The Gaussian beam superposition method is an asymptotic method for computing
high frequency wave fields in smoothly varying inhomogeneous media. In this
paper we study the accuracy of the Gaussian beam superposition method and
derive error estimates related to the discretization of the superposition
integral and the Taylor expansion of the phase and amplitude off the center of
the beam. We show that in the case of odd order beams, the error is smaller
than a simple analysis would indicate because of error cancellation effects
between the beams.