No description available.

Thesis (Ph.D.)--Boston University / PLEASE NOTE: Boston University Libraries did not receive an Authorization To Manage form for this thesis or dissertation. It is therefore not openly accessible, though it may be available by request. If you are the author or principal advisor of this work and would like to request open access for it, please contact us at open-help@bu.edu. Thank you. / This thesis studies the question of control of Maxwell's equations in a medium with positive conductivity by means of boundary surface currents. Two types of domains and media are considered in connection with this question. First is a bounded simply connected star-shaped domain in R^3 which is made up of a heterogeneous medium with small conductivity, with controls being applied over the entire boundary. Using the Hilbert Uniqueness Method of Lions, the exact boundary controllability over a sufficiently long time period is established for this case, provided the conductivity is small enough to satisfy a certain technical inequality. It is also found that the requirement for the conductivity term to be very small remains in place even if the medium considered is homogenous. In order to remove this constraint, a special domain type is considered next - a cube - made up of a homogenous medium where the conductivity is allowed to take on any non-negative value. An additional restriction imposed here in order to make this problem more suitable for practical implementations is that the controls are applied over only one face of the cube. Employing the Method of Moments the spectral controllability is established for this case. It is also established that the exact controllability fails for this geometry regardless of the size of the conductivity term. This thesis will also consider the question of reconstructing the source of electromagnetic radiation, which is related to the controllability problem. / 2031-01-02

Application and implication of using complex vectors and complex transformations in solutions of Maxwell’s equations is investigated. Complex vectors are used in complex plane waves and help to represent this type of waves geometrically. It is shown that they are also useful in representing inhomogeneous plane waves in plasma, single-negative and double-negative metamaterials. In specific I will investigate the Otto configuration and Kretschmann configuration and I will show that in order to observe the minimum in reflection coefficient it is necessary for the metal to be lossy. We will compare this to the case of plasmon-like resonance when a PEC periodic structure is illuminated by a plane wave. Complex transformations are crucial in deriving Gaussian beam solutions of paraxial Helmholtz equation from spherical wave solution of Helmholtz equation. Vector Gaussian beams also will be discussed shortly.

Application and implication of using complex vectors and complex transformations in solutions of Maxwell’s equations is investigated. Complex vectors are used in complex plane waves and help to represent this type of waves geometrically. It is shown that they are also useful in representing inhomogeneous plane waves in plasma, single-negative and double-negative metamaterials. In specific I will investigate the Otto configuration and Kretschmann configuration and I will show that in order to observe the minimum in reflection coefficient it is necessary for the metal to be lossy. We will compare this to the case of plasmon-like resonance when a PEC periodic structure is illuminated by a plane wave. Complex transformations are crucial in deriving Gaussian beam solutions of paraxial Helmholtz equation from spherical wave solution of Helmholtz equation. Vector Gaussian beams also will be discussed shortly.

No description available.

<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

In the Large Hadron Collider (LHC), particle physics events are created by colliding high energy proton beams at a number of interaction points around the ring. One of the main performance indicating parameters of the LHC is the luminosity. The luminosity is limited by, amongst other things, the strength of the beam-beam interaction. In this thesis, the effect of the beam-beam interaction on the luminosity performance of the LHC and the proposed High Luminosity Large Hadron Collider (HL-LHC) is investigated. Results from a number of dedicated, long-range beam-beam machine studies are presented and analysed. In these studies, the minimum beam-beam separation for two different beta star— optics are identified. This separation defines the minimum operational crossing angle in the LHC. The data from these studies are then compared to simulation of the dynamic aperture and the results are discussed. In addition to studies of the LHC, an analytical approach is derived in order to describe the hourglass effect, which may become a contributing factor in limiting the luminosity performance of the HL-LHC.

Herein a fully parallel, upwind and flux-split Finite-Volume Time-Domain (FVTD) numerical engine for solving Maxwell's Equations on unstructured grids is developed. The required background theory for solving Maxwell's Equations using FVTD is given in sufficient detail, including a description of both the temporal and spatial approximations used. The details of the local-time stepping strategy of Fumeaux et al. is included. A global mesh-truncation scheme using field integration over a Huygens' surface is also presented. The capabilities of the FVTD algorithm are augmented with thin-wire and subcell circuit models that permit very flexible and accurate simulations of circuit-driven wire structures. Numerical and experimental validation shows that the proposed models have a wide-range of applications. Specifically, it appears that the thin-wire and subcell circuit models may be very well suited to the simulation of radio-frequency coils used in magnetic resonance imaging systems. A parallelization scheme for the volumetric field solver, combined with the local-time stepping, global mesh-truncation and subcell models is developed that theoretically provides both linear time- and memory scaling in a distributed parallel environment. Finally, the FVTD code is converted to the frequency domain and the possibility of using different flux-reconstruction schemes to improve the iterative convergence of the Finite-Volume Frequency-Domain algorithm is investigated.

Herein a fully parallel, upwind and flux-split Finite-Volume Time-Domain (FVTD) numerical engine for solving Maxwell's Equations on unstructured grids is developed. The required background theory for solving Maxwell's Equations using FVTD is given in sufficient detail, including a description of both the temporal and spatial approximations used. The details of the local-time stepping strategy of Fumeaux et al. is included. A global mesh-truncation scheme using field integration over a Huygens' surface is also presented. The capabilities of the FVTD algorithm are augmented with thin-wire and subcell circuit models that permit very flexible and accurate simulations of circuit-driven wire structures. Numerical and experimental validation shows that the proposed models have a wide-range of applications. Specifically, it appears that the thin-wire and subcell circuit models may be very well suited to the simulation of radio-frequency coils used in magnetic resonance imaging systems. A parallelization scheme for the volumetric field solver, combined with the local-time stepping, global mesh-truncation and subcell models is developed that theoretically provides both linear time- and memory scaling in a distributed parallel environment. Finally, the FVTD code is converted to the frequency domain and the possibility of using different flux-reconstruction schemes to improve the iterative convergence of the Finite-Volume Frequency-Domain algorithm is investigated.

Cette thèse présente une méthode parallèle de résolution de systèmes linéaires creux basée sur un algorithme multigrille géométrique. Les estimations de la solution sont calculées par méthode directe sur le niveau grossier ou par méthode itérative de type splitting sur les maillages raffinés; des opérateurs inter-grilles sont définis pour interpoler les solutions approximatives entre les différents niveaux de raffinements. Ce solveur est utilisé dans le cadre de simulations électromagnétiques en 3D (équations de Maxwell en régime harmonique discrétisées par éléments finis de Nédélec de premier ordre) en tant que méthode stationnaire ou comme préconditionneur d’une méthode de Krylov (GMRES). / Multigrid algorithm. The system is solved thanks to a direct method on the coarse mesh anditerative splitting method on refined meshes; inter-grid operators are defined to interpolate theapproximate solutions on the different refinement levels. Applied to 3D electromagnetic simulations(Nédélec first order finite element approximation of time harmonic Maxwell equations) thissolver is used either as a stationary method or as a preconditioner for a Krylov subspace method(GMRES).