Perfectly Matched Layers in a Divergence Preserving ADI Scheme for Electromagnetics

For numerical simulations of highly relativistic and transversely accelerated charged particles including radiation fast algorithms are needed. While the radiation in particle accelerators has wavelengths in the order of 100 um the computational domain has dimensions roughly 5 orders of magnitude larger resulting in very large mesh sizes. The particles are confined to a small area of this domain only. To resolve the smallest scales close to the particles subgrids are envisioned. For reasons of stability the alternating direction implicit (ADI) scheme by D. N. Smithe et al. (J. Comput. Phys. 228 (2009) pp.7289-7299) for Maxwell equations has been adopted. At the boundary of the domain absorbing boundary conditions have to be employed to prevent reflection of the radiation. In this paper we show how the divergence preserving ADI scheme has to be formulated in perfectly matched layers (PML) and compare the performance in several scenarios.Comment: 8 pages, 6 figure

Method of lines transpose: an efficient A-stable solver for wave propagation

Building upon recent results obtained in [7,8,9], we describe an efficient second order, A-stable scheme for solving the wave equation, based on the method of lines transpose (MOLT), and the resulting semi-discrete (i.e. continuous in space) boundary value problem. In [7], A-stable schemes of high order were derived, and in [9] a high order, fast O(N) spatial solver was derived, which is matrix-free and is based on dimensional-splitting. In this work, are interested in building a wave solver, and our main concern is the development of boundary conditions. We demonstrate all desired boundary conditions for a wave solver, including outflow boundary conditions, in 1D and 2D. The scheme works in a logically Cartesian fashion, and the boundary points are embedded into the regular mesh, without incurring stability restrictions, so that boundary conditions are imposed without any reduction in the order of accuracy. We demonstrate how the embedded boundary approach works in the cases of Dirichlet and Neumann boundary conditions. Further, we develop outflow and periodic boundary conditions for the MOLT formulation. Our solver is designed to couple with particle codes, and so special attention is also paid to the implementation of point sources, and soft sources which can be used to launch waves into waveguides

System- and Data-Driven Methods and Algorithms

An increasing complexity of models used to predict real-world systems leads to the need for algorithms to replace complex models with far simpler ones, while preserving the accuracy of the predictions. This two-volume handbook covers methods as well as applications. This first volume focuses on real-time control theory, data assimilation, real-time visualization, high-dimensional state spaces and interaction of different reduction techniques

SCEE 2008 book of abstracts : the 7th International Conference on Scientific Computing in Electrical Engineering (SCEE 2008), September 28 – October 3, 2008, Helsinki University of Technology, Espoo, Finland

This report contains abstracts of presentations given at the SCEE 2008 conference.reviewe

Simulation de la propagation d'ondes électromagnétiques en nano-optique par une méthode Galerkine discontinue d'ordre élevé

The goal of this thesis is to develop a discontinuous Galerkin time-domain method to be able to handle realistic nanophotonics computations. During the last decades, the evolution of lithography techniques allowed the creation of geometrical structures at the nanometer scale, thus unveiling a variety of new phenomena arising from light-matter interactions at such levels. These effects usually occur when the device is of comparable size or (much) smaller than the wavelength of the incident field. This work relies on the development and implementation of appropriate models for dispersive materials (mostly metals), as well as on a large panel of classical computational techniques. Two major methodological developments are presented and studied in details: (i) curvilinear elements, and (ii) local order of approximation. This work is complemented with several physical studies of real-life nanophotonics applications.L’objectif de cette thèse est de développer une méthode Galerkine discontinue d’ordre élevé capable de prendre en considération des simulations réalistes liées à la nanophotonique. Au cours des dernières décennies, l’évolution des techniques de lithographie a permis la création de structure géométriques de tailles nanométriques, révélant ainsi une large gamme de phénomènes nouveaux nés de l’interaction lumière-matière à ces échelles. Ces effets apparaissent généralement pour des objets de taille égale ou (très) inférieure à la longueur d’onde du champ incident. Ce travail repose sur le développement et l’implémentation de modèles de dispersion appropriés (principalement pour les métaux), ainsi que sur un large éventail de méthodes computationnelles classiques. Deux développements méthodologiques majeurs sont présentés et étudiés en détails: (i) les éléments courbes, et (ii) l’ordre d’approximation local. Ces études sont accompagnées de plusieurs cas-tests réalistes tirés de la nanophotonique

Fully coupled CEM/CFD modelling of microwave heating in a porous medium

Computational results for the microwave heating of a porous material are presented in this paper. Coupled finite difference time domain and finite volume methods are used to solve equations that describe the electromagnetic field and heat and mass transfer in porous media. These equations are nonlinearly coupled through the dielectric properties which depend both on temperature and moisture content. By investigating the resonant behaviour in two-dimensional microwave cavities, the FD-TD scheme is validated. Validation of the microwave power distribution in 3-D microwave enclosures is compared with other numerical results available. 3-D temperature distribution in a biomaterial is validated against experimental results. Results using the proposed fully coupled approach are discussed and analyzed. The model is able to reflect the evolution of both temperature and moisture fields as well as energy penetration as the moisture in the porous medium evaporates. Moisture movement results from internal pressure gradients produced by the internal heating and phase change. The model is validated by comparison to some published results for simpler problems

Suffolk University Academic Catalog, College of Arts and Sciences and Sawyer Business School, 2012-2013

This catalog contains information for both the undergraduate and graduate programs. The catalog is a pdf version of the Suffolk website, and the pages are not in numerical order. In order to navigate the program descriptions and degree requirements, please use the menu links on the left side of each page. To view the course descriptions, refer to the A-Z list of courses starting on page 1,220 (these are also included here as separate pdf files with lists of CAS and SBS courses). You can also search for a particular course or program by clicking ctrl+f and typing in the course number or name. Please contact the Archives if you need assistance navigating this catalog or finding information on degree requirements or course descriptions.https://dc.suffolk.edu/cassbs-catalogs/1165/thumbnail.jp