Lacunae Based Stabilization of PMLs

hyperbolic PDE;"Maxwells Equations";Artificial Boundary Conditions;Lacunae;Absorbing Layers;PML

Qasimov, HeydarRafiq oglu ; Pierre Gremaud, Committee Member,Shamim Rahman, Committee Member,Dmitri Zenkov, Committee Member,Ralph Smith, Committee Member,Semyon Tsynkov, Committee Chair,Qasimov, HeydarRafiq oglu ; Pierre Gremaud ; Committee Member ; Shamim Rahman ; Committee Member ; Dmitri Zenkov ; Committee Member ; Ralph Smith ; Committee Member ; Semyon Tsynkov ; Committee Chair

University:North Carolina State University

关键词: hyperbolic PDE;    "Maxwells Equations";    Artificial Boundary Conditions;    Lacunae;    Absorbing Layers;    PML;   

Others  :  https://repository.lib.ncsu.edu/bitstream/handle/1840.16/4820/etd.pdf?sequence=1&isAllowed=y

美国|英语

来源: null

【 摘 要 】

Perfectly matched layers (PML) enclose the computational domain for simulating electromagnetic phenomena defined over unbounded regions. While being overall very successful, this procedure has sometimes been reported to develop instabilities and exhibit a physically unaccountable growth of the solution inside the layer over long integrationtimes. In the thesis, we conduct a numerical as well as analytical study of the PML's response after a long time of integration. The physical and mathematical PMLsare implemented with several well-known schemes: Yee, leap-frog, Lax-Wendroff, and Runge-Kutta in time with central differences in space. Then, the eigen-structure of each discretization at quiescent state is studied to gain an insight into the nature of instability, the sources of growth of the solution and the potential contamination of the domain of interest. The results of this investigation provide useful information regarding the better and worse performers among thespecific combinations of schemes and PMLs, yet they do not precisely identify the mechanism behind the growth of the solution. Therefore, the main focus of the thesis is to build a methodology that would inhibit the instability of the PML regardless of its source. The approach is based on the concept of numerical integration that exploitsthe presence of lacunae in the solutions. It applies to hyperbolic partial differential equations and systemsthat satisfy the Huygens' principle, in particular, the Maxwell's system of equations that governs the propagation of electromagnetic waves. The methodology does not modify the equations inside the layer and hence, while eliminating the undesirable growth, it fully preserves all the advantageous properties of a given PML, such as matching at the interface and the degree of absorption.A practical algorithm is constructed in the thesis, and atheorem is proved that guarantees a temporally uniform error bound over the domain of interest. The theoretical findings are corroborated numerically, and the potential for extending the methodology is discussed