Browsing by Subject "Perfectly matched layers"
Now showing 1 - 2 of 2
Results Per Page
Sort Options
Item Fast algorithms for frequency domain wave propagation(2012-12) Tsuji, Paul Hikaru; Ying, Lexing; Ghattas, Omar N.; Engquist, Bjorn; Fomel, Sergey; Ren, KuiHigh-frequency wave phenomena is observed in many physical settings, most notably in acoustics, electromagnetics, and elasticity. In all of these fields, numerical simulation and modeling of the forward propagation problem is important to the design and analysis of many systems; a few examples which rely on these computations are the development of metamaterial technologies and geophysical prospecting for natural resources. There are two modes of modeling the forward problem: the frequency domain and the time domain. As the title states, this work is concerned with the former regime. The difficulties of solving the high-frequency wave propagation problem accurately lies in the large number of degrees of freedom required. Conventional wisdom in the computational electromagnetics commmunity suggests that about 10 degrees of freedom per wavelength be used in each coordinate direction to resolve each oscillation. If K is the width of the domain in wavelengths, the number of unknowns N grows at least by O(K^2) for surface discretizations and O(K^3) for volume discretizations in 3D. The memory requirements and asymptotic complexity estimates of direct algorithms such as the multifrontal method are too costly for such problems. Thus, iterative solvers must be used. In this dissertation, I will present fast algorithms which, in conjunction with GMRES, allow the solution of the forward problem in O(N) or O(N log N) time.Item The inverse medium problem in PML-truncated elastic media(2010-12) Kucukcoban, Sezgin; Kallivokas, Loukas F.; Demkowicz, Leszek F.; Ghattas, Omar; Kinnas, Sypros A.; Rodin, Gregory J.; Torres-Verdin, CarlosWe introduce a mathematical framework for the inverse medium problem arising commonly in geotechnical site characterization and geophysical probing applications, when stress waves are used to probe the material composition of the interrogated medium. Specifically, we attempt to recover the spatial distribution of Lame's parameters ( and μ) of an elastic semi-infinite arbitrarily heterogeneous medium, using surface measurements of the medium's response to prescribed dynamic excitations. The focus is on characterizing near-surface deposits, and to this end, we develop a method that is implemented directly in the time-domain, is driven by the full waveform response collected at receivers on the surface, while the domain of interest is truncated using Perfectly-Matched-Layers (PMLs) to limit the originally semi-infinite extent of the physical domain. There are two key issues associated with the problem at hand: (a) the forward problem, namely the numerical simulation of the wave motion in the domain of interest; and (b) the framework and strategies for tackling the inverse problem. To address the forward problem, it is necessary that the domain of interest be truncated, and the resulting finite domain be forced to mimic the physics of the original problem: to this end, we introduce unsplit-field PMLs, and develop and implement two new formulations, one fully-mixed and one hybrid (mixed coupled with a non-mixed approach) that model wave motion within the, now PML-truncated, domain. To address the inverse problem, we adopt a partial-differential-equation-constrained optimization framework that results in the usual triplet of an initial-and-boundary-value forward problem, a final-and-boundary-value adjoint problem, and a time-independent boundary-value control problem. This triplet of boundary-value-problems is used to guide the optimizer to the target profile of the spatially distributed Lame parameters. Given the multiplicity of solutions, we assist the optimizer, by deploying regularization schemes, continuation schemes (regularization factor and source-frequency content), as well as a physics-driven simple procedure to bias the search directions. We report numerical examples attesting to the quality, stability, and efficiency of the forward wave modeling. We also report moderate success with numerical experiments targeting inversion of both smooth and sharp profiles in two dimensions.