<p>Domain decomposition method provides a solution for the very large electromagnetic</p><p>system which are impossible for single domain methods. Discontinuous Galerkin</p><p>(DG) method can be viewed as an extreme version of the domain decomposition,</p><p>i.e., each element is regarded as one subdomain. The whole system is solved element</p><p>by element, thus the inversion of the large global system matrix is no longer necessary,</p><p>and much larger system can be solved with the DG method compared to the</p><p>continuous Galerkin (CG) method.</p><p>In this work, the DG method is implemented on a subdomain level, that is, each subdomain contains multiple elements. The numerical flux only applies on the</p><p>interfaces between adjacent subdomains. The subodmain level DG method divides</p><p>the original large global system into a few smaller ones, which are easier to solve,</p><p>and it also provides the possibility of parallelization. Compared to the conventional</p><p>element level DG method, the subdomain level DG has the advantage of less total</p><p>DoFs and fexibility in interface choice. In addition, the implicit time stepping is </p><p>relatively much easier for the subdomain level DG, and the total CPU time can be</p><p>much less for the electrically small or multiscale problems.</p><p>The hybrid of elements are employed to reduce the total DoF of the system.</p><p>Low-order tetrahedrons are used to catch the geometry ne parts and high-order</p><p>hexahedrons are used to discretize the homogeneous and/or geometry coarse parts.</p><p>In addition, the non-conformal mesh not only allow dierent kinds of elements but</p><p>also sharp change of the element size, therefore the DoF can be further decreased.</p><p>The DGTD method in this research is based on the EB scheme to replace the</p><p>previous EH scheme. Dierent from the requirement of mixed order basis functions</p><p>for the led variables E and H in the EH scheme, the EB scheme can suppress the</p><p>spurious modes with same order of basis functions for E and B. One order lower in</p><p>the basis functions in B brings great benets because the DoFs can be signicantly</p><p>reduced, especially for the tetrahedrons parts.</p><p>With the basis functions for both E and B, the EB scheme upwind </p><p>ux and</p><p>EB scheme Maxwellian PML, the eigen-analysis and numerical results shows the</p><p>eectiveness of the proposed DGTD method, and multiscale problems are solved</p><p>eciently combined with the implicit-explicit hybrid time stepping scheme and multiple</p><p>kinds of elements.</p><p>The EB scheme DGTD method is further developed to allow arbitrary anisotropic</p><p>media via new anisotropic EB scheme upwind </p><p>ux and anisotropic EB scheme</p><p>Maxwellian PML. The anisotropic M-PML is long time stable and absorb the outgoing</p><p>wave eectively. A new TF/SF boundary condition is brought forward to</p><p>simulate the half space case. The negative refraction in YVO4 bicrystal is simulated</p><p>with the anisotropic DGTD and half space TF/SF condition for the rst time with</p><p>numerical methods.</p> / Dissertation
Identifer | oai:union.ndltd.org:DUKE/oai:dukespace.lib.duke.edu:10161/11357 |
Date | January 2015 |
Creators | Ren, Qiang |
Contributors | Liu, Qing H |
Source Sets | Duke University |
Detected Language | English |
Type | Dissertation |
Page generated in 0.0118 seconds