Browsing by Subject "Wave scattering"
Now showing 1 - 3 of 3
Results Per Page
Sort Options
Item Fast numerical methods for high frequency wave scattering(2012-05) Tran, Khoa Dang; Engquist, Björn, 1945-; Ling, Hao; Ghattas, Omar; Tsai, Richard; Ying, LexingComputer simulation of wave propagation is an active research area as wave phenomena are prevalent in many applications. Examples include wireless communication, radar cross section, underwater acoustics, and seismology. For high frequency waves, this is a challenging multiscale problem, where the small scale is given by the wavelength while the large scale corresponds to the overall size of the computational domain. Research into wave equation modeling can be divided into two regimes: time domain and frequency domain. In each regime, there are two further popular research directions for the numerical simulation of the scattered wave. One relies on direct discretization of the wave equation as a hyperbolic partial differential equation in the full physical domain. The other direction aims at solving an equivalent integral equation on the surface of the scatterer. In this dissertation, we present three new techniques for the frequency domain, boundary integral equations.Item Multiple-grid adaptive integral method for general multi-region problems(2011-08) Wu, Mingfeng; Yilmaz, Ali E.; Ling, Hao; Pearce, John; Alu, Andrea; Ying, LexingEfficient electromagnetic solvers based on surface integral equations (SIEs) are developed for the analysis of scattering from large-scale and complex composite structures that consist of piecewise homogeneous magnetodielectric and perfect electrically/magnetically conducting (PEC/PMC) regions. First, a multiple-grid extension of the adaptive integral method (AIM) is presented for multi-region problems. The proposed method accelerates the iterative method-of-moments solution of the pertinent SIEs by employing multiple auxiliary Cartesian grids: If the structure of interest is composed of K homogeneous regions, it introduces K different auxiliary grids. It uses the k^{th} auxiliary grid first to determine near-zones for the basis functions and then to execute AIM projection/anterpolation, propagation, interpolation, and near-zone pre-correction stages in the k^{th} region. Thus, the AIM stages are executed a total of K times using different grids and different groups of basis functions. The proposed multiple-grid AIM scheme requires a total of O(N^{nz,near}+sum({N_k}^Clog{N_k}^C)) operations per iteration, where N^{nz,near} denotes the total number of near-zone interactions in all regions and {N_k}^C denotes the number of nodes of the k^{th} Cartesian grid. Numerical results validate the method’s accuracy and reduced complexity for large-scale canonical structures with large numbers of regions (up to 10^6 degrees of freedom and 10^3 regions). Then, a Green function modification approach and a scheme of Hankel- to Teoplitz-matrix conversions are efficiently incorporated to the multiple-grid AIM method to account for a PEC/PMC plane. Theoretical analysis and numerical examples show that, compared to a brute-force imaging scheme, the Green function modification approach reduces the simulation time and memory requirement by a factor of (almost) two or larger if the structure of interest is terminated on or resides above the plane, respectively. In addition, the SIEs are extended to cover structures composed of metamaterial regions, PEC regions, and PEC-material junctions. Moreover, recently introduced well-conditioned SIEs are adopted to achieve faster iterative solver convergence. Comprehensive numerical tests are performed to evaluate the accuracy, computational complexity, and convergence of the novel formulation which is shown to significantly reduce the number of iterations and the overall computational work. Lastly, the efficiency and capabilities of the proposed solvers are demonstrated by solving complex scattering problems, specifically those pertinent to analysis of wave propagation in natural forested environments, the design of metamaterials, and the application of metamaterials to radar cross section reduction.Item Wave Interactions with Arrays of Bottom-Mounted Circular Cylinders: Investigation of Optical and Acoustical Analogies(2010-10-12) Baquet, AldricWave scattering by arrays of cylinders has received special attention by many authors and analytical solutions have been derived. The investigation of optical and acoustical analogies to the problem of interaction of water waves with rigid and flexible cylinder arrays is the main focus of this thesis. In acoustics, a sound may be attenuated while it propagates through a layer of bubbly liquid. In fact, if the natural frequency of the bubbles is in the range of the wave periods, the attenuation becomes more evident. The ultimate objective of the research described herein is to determine if this phenomenon may also be found in the interaction between water waves and arrays of flexible cylinders. In a first approach, arrays of rigid cylinders are studied in shallow water. The array is treated as an effective medium, which allows for the definition of reflection and transmission coefficients for the array, and theories from Hu and Chan (2005) associated with the Fabry-Perot interferometer are compared against direct computations of wave scattering using the commercial code WAMIT. Reflection and transmission coefficients from WAMIT are evaluated by applying a Maximum Likelihood Method. The results from WAMIT were found to be in good agreement with those obtained from the effective medium theory. Due to observed inconsistencies for short wave periods and small incident angles, the effective width of the medium is defined and corrected. For the case of a flexible cylinder, generalized modes corresponding to deformations of the cylinder's surface are formulated and added to WAMIT's subroutine. Equations of motion are derived from the theory of vibration for thin shells and mass and stiffness matrices are defined. The objective is to maximize wave attenuation from the array of flexible cylinders. Therefore, the natural periods of the "breathing" mode for these cylinders is set in the range of the studied wave periods. Then, material properties, as well as mass and stiffness matrices, are chosen to achieve this effect.