Some Key Developments in Computational Electromagnetics and their Attribution

Key developments in computational electromagnetics are proposed. Historical highlights are summarized concentrating on the two main approaches of differential and integral methods. This is seen as timely as a retrospective analysis is needed to minimize duplication and to help settle questions of attribution

Hypercube matrix computation task

A major objective of the Hypercube Matrix Computation effort at the Jet Propulsion Laboratory (JPL) is to investigate the applicability of a parallel computing architecture to the solution of large-scale electromagnetic scattering problems. Three scattering analysis codes are being implemented and assessed on a JPL/California Institute of Technology (Caltech) Mark 3 Hypercube. The codes, which utilize different underlying algorithms, give a means of evaluating the general applicability of this parallel architecture. The three analysis codes being implemented are a frequency domain method of moments code, a time domain finite difference code, and a frequency domain finite elements code. These analysis capabilities are being integrated into an electromagnetics interactive analysis workstation which can serve as a design tool for the construction of antennas and other radiating or scattering structures. The first two years of work on the Hypercube Matrix Computation effort is summarized. It includes both new developments and results as well as work previously reported in the Hypercube Matrix Computation Task: Final Report for 1986 to 1987 (JPL Publication 87-18)

A simulation method for determining the optical response of highly complex photonic structures of biological origin

We present a method based on a time domain simulation of wave propagation that allows studying the optical response of a broad range of dielectric photonic structures. This method is particularly suitable for dealing with complex biological structures. One of the main features of the proposed approach is the simple and intuitive way of defining the setup and the photonic structure to be simulated, which can be done by feeding the simulation with a digital image of the structure. We also develop a set of techniques to process the behavior of the evolving waves within the simulation. These techniques include a direction filter, that permits decoupling of waves travelling simultaneously in different directions, a dynamic differential absorber, to cancel the waves reflected at the edges of the simulation space, a multi-frequency excitation scheme based on a filter that allows decoupling waves of different wavelengths travelling simultaneously, and a near-to-far-field approach to evaluate the resulting wavefield outside the simulation domain. We validate the code and, as an example, apply it to the complex structure found in a microorganism called Diachea leucopoda, which exhibits a multicolor iridescent appearance.Comment: 43 pages, 19 figure

Vademecum-based GFEM (V-GFEM): optimal enrichment for transient problems

This is the accepted version of the following article: [Canales, D., Leygue, A., Chinesta, F., González, D., Cueto, E., Feulvarch, E., Bergheau, J. -M., and Huerta, A. (2016) Vademecum-based GFEM (V-GFEM): optimal enrichment for transient problems. Int. J. Numer. Meth. Engng, 108: 971–989. doi: 10.1002/nme.5240.], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5240/fullThis paper proposes a generalized finite element method based on the use of parametric solutions as enrichment functions. These parametric solutions are precomputed off-line and stored in memory in the form of a computational vademecum so that they can be used on-line with negligible cost. This renders a more efficient computational method than traditional finite element methods at performing simulations of processes. One key issue of the proposed method is the efficient computation of the parametric enrichments. These are computed and efficiently stored in memory by employing proper generalized decompositions. Although the presented method can be broadly applied, it is particularly well suited in manufacturing processes involving localized physics that depend on many parameters, such as welding. After introducing the vademecum-generalized finite element method formulation, we present some numerical examples related to the simulation of thermal models encountered in welding processes.Peer ReviewedPostprint (author's final draft

Electromagnetic modeling and simulations of interest for low-level laser therapies or for problems involving moving objects

This work presents results on the computed electromagnetic field within isolated mitochondria when exposed to near-infrared illuminations concerning photobiomodulation experiments. The accuracy of the electromagnetic models implemented for dosimetry is important. The mechanism of interaction of light with these organelles is still unclear, so it is important to improve our knowledge with reliable simulations and experiments. To obtain such results, we present several models. Although they refer to a well-defined experimental setup, the different models must consider the different possible arrangements of the mitochondria and the differences in their dimensions and constitutive parameters. Different wavelengths and polarizations are considered too. The effects of all the parameters on the electromagnetic field inside the mitochondria and the internal morphology are studied. Computational techniques are also exploited to find reliable approximations of the solutions of problems involving moving objects. The effect of the motion of objects on the electromagnetic field is studied as well; for this type of problem, it is necessary to understand that the constitutive relations are modified because the movement makes any material to appear as a bianisotropic medium. In particular, the rotation movement of symmetrical bodies is studied, the electromagnetic field is calculated and an inversion algorithm for the estimate of the rotation speed is proposed

Fourth-order finite volume algorithm with adaptive mesh refinement in space and time for multi-fluid plasma modeling, A

2022 Spring.Includes bibliographical references.Improving our fundamental understanding of plasma physics using numerical methods is pivotal to the advancement of science and the continual development of cutting-edge technologies such as nuclear fusion reactions for energy production or the manufacturing of microelectronic devices. An elaborate and accurate approach to modeling plasmas using computational fluid dynamics (CFD) is the multi-fluid method, where the full set of fluid mechanics equations are solved for each species in the plasma simultaneously with Maxwell's equations in a coupled fashion. Nevertheless, multi-fluid plasma modeling is inherently multiscale and multiphysics, presenting significant numerical and mathematical stiffness. This research aims to develop an efficient and accurate multi-fluid plasma model using higher-order, finite-volume, solution-adaptive numerical methods. The algorithm developed herein is verified to be fourth-order accurate for electromagnetic simulations as well as those involving fully-coupled, multi-fluid plasma physics. The solutions to common plasma test problems obtained by the algorithm are validated against exact solutions and results from literature. The algorithm is shown to be robust and stable in the presence of complex solution topology and discontinuities, such as shocks and steep gradients. The optimizations in spatial discretization provided by the fourth-order algorithm and adaptive mesh refinement are demonstrated to improve the solution time by a factor of 10 compared to lower-order methods on fixed-grid meshes. This research produces an advanced, multi-fluid plasma modeling framework which allows for studying complex, realistic plasmas involving collisions and practical geometries