Nonlinear full wave time domain solutions using FDTD_SPICE for high speed digital

Abstract With high speed digital designs extending well into the GHz range, quasi-static and TEM solutions are no longer valid. Digital designers like to see time domain waveforms as provided by SPICE. However, full wave effects are often accounted for only by frequency domain electromagnetic solvers. A simulator that combines FDTD and SPICE is used in this paper as a solution to this problem. The nonlinear circuit elements are simulated in SPICE and the structural elements in FDTD. The simulators are tightly integrated and talk to each other at each time-step. The approach is used to analyze difficult problems such as meanders, ground bounce and placement of decoupling capacitors in a PCB, and dispersion in microstrips. Authors/Speakers Neven Orhanovic Current Activities He is currently with Applied Simulation Technology working on time domain full wave methods. Backgroun

Recent development of surface integral equation solvers for multiscale interconnects and circuits

This paper presents a brief review and recent development of surface integral equation solvers for multiscale interconnects and circuits modeling. As the future production processes down to 5 nm and the operating frequency increases, both multi-scale and large-scale natures should be taken into account in the electromagnetic simulations. Fast, efficient, stable, and broadband integral equation based solvers become indispensable when millions or ten s of millions of unknowns might be involved in the simulation of the integrated circuit. Recent progress and our latest researches in the development of broadband fast electromagnetic solvers will be demonstrated.published_or_final_versio

Parallel 3-D marine controlled-source electromagnetic modelling using high-order tetrahedral Nédélec elements

We present a parallel and high-order Nédélec finite element solution for the marine controlled-source electromagnetic (CSEM) forward problem in 3-D media with isotropic conductivity. Our parallel Python code is implemented on unstructured tetrahedral meshes, which support multiple-scale structures and bathymetry for general marine 3-D CSEM modelling applications. Based on a primary/secondary field approach, we solve the diffusive form of Maxwell’s equations in the low-frequency domain. We investigate the accuracy and performance advantages of our new high-order algorithm against a low-order implementation proposed in our previous work. The numerical precision of our high-order method has been successfully verified by comparisons against previously published results that are relevant in terms of scale and geological properties. A convergence study confirms that high-order polynomials offer a better trade-off between accuracy and computation time. However, the optimum choice of the polynomial order depends on both the input model and the required accuracy as revealed by our tests. Also, we extend our adaptive-meshing strategy to high-order tetrahedral elements. Using adapted meshes to both physical parameters and high-order schemes, we are able to achieve a significant reduction in computational cost without sacrificing accuracy in the modelling. Furthermore, we demonstrate the excellent performance and quasi-linear scaling of our implementation in a state-of-the-art high-performance computing architecture.This project has received funding from the European Union's Horizon 2020 programme under the Marie Sklodowska-Curie grant agreement No. 777778. Furthermore, the research leading to these results has received funding from the European Union's Horizon 2020 programme under the ChEESE Project (https://cheese-coe.eu/ ), grant agreement No. 823844. In addition, the authors would also like to thank the support of the Ministerio de Educación y Ciencia (Spain) under Projects TEC2016-80386-P and TIN2016-80957-P. The authors would like to thank the Editors-in-Chief and to both reviewers, Dr. Martin Cuma and Dr. Raphael Rochlitz, for their valuable comments and suggestions which helped to improve the quality of the manuscript. This work benefited from the valuable suggestions, comments, and proofreading of Dr. Otilio Rojas (BSC). Last but not least, Octavio Castillo-Reyes thanks Natalia Gutierrez (BSC) for her support in CSEM modeling with BSIT.Peer ReviewedPostprint (author's final draft

Advanced electrode models and numerical modelling for high frequency Electrical Impedance Tomography systems

The thesis discusses various electrode models and finite element analysis methods for Electrical Impedance Tomography (EIT) systems. EIT is a technique for determining the distribution of the conductivity or admittivity in a volume by injecting electrical currents into the volume and measuring the corresponding potentials on the surface of the volume. Various electrode models were investigated for operating EIT systems at higher frequencies in the beta-dispersion band. Research has shown that EIT is potentially capable to distinguish malignant and benign tumours in this frequency band. My study concludes that instrumental effects of the electrodes and full Maxwell effects of EIT systems are the major issues, and they have to be addressed when the operating frequency increases. In the thesis, I proposed 1) an Instrumental Electrode Model (IEM) for the quasi-static EIT formula, based on the analysis of the hardware structures attached to electrodes; 2) a Complete Electrode Model based on Impedance Boundary Conditions (CEM-IBC) that introduces the contact impedances into the full Maxwell EIT formula; 3) a Transmission line Port Model (TPM) for electrode pairs with the instrumental effects, the contact impedance, and the full Maxwell effects considered for EIT systems. Circuit analysis, Partial Differential Equations (PDE) analysis, numerical analysis and finite element methods were used to develop the models. The results obtained by the proposed models are compared with widely used Commercial PDE solvers. This thesis addresses the two major problems (instrumental effects of the electrodes and full Maxwell effects of EIT systems) with the proposed advanced electrode models. Numerical experiments show that the proposed models are more accurate in the high frequency range of EIT systems. The proposed electrode models can be also applicable to inverse problems, and the results show promising. Simple hardware circuits for verifying the results experimentally have been also designed

Quasi-linear analysis of the extraordinary electron wave destabilized by runaway electrons

Runaway electrons with strongly anisotropic distributions present in post-disruption tokamak plasmas can destabilize the extraordinary electron (EXEL) wave. The present work investigates the dynamics of the quasi-linear evolution of the EXEL instability for a range of different plasma parameters using a model runaway distribution function valid for highly relativistic runaway electron beams produced primarily by the avalanche process. Simulations show a rapid pitch-angle scattering of the runaway electrons in the high energy tail on the $100-1000\;\rm \mu s$ time scale. Due to the wave-particle interaction, a modification to the synchrotron radiation spectrum emitted by the runaway electron population is foreseen, exposing a possible experimental detection method for such an interaction

GRMHD in axisymmetric dynamical spacetimes: the X-ECHO code

We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an important role, such as for the evolution of magnetized neutron stars or for the gravitational collapse of the magnetized rotating cores of massive stars, which is the astrophysical scenario believed to eventually lead to (long) GRB events. The code is based on the extension of the Eulerian conservative high-order (ECHO) scheme [Del Zanna et al., A&A 473, 11 (2007)] for GRMHD, here coupled to a novel solver for the Einstein equations in the extended conformally flat condition (XCFC). We fully exploit the 3+1 Eulerian formalism, so that all the equations are written in terms of familiar 3D vectors and tensors alone, we adopt spherical coordinates for the conformal background metric, and we consider axisymmetric spacetimes and fluid configurations. The GRMHD conservation laws are solved by means of shock-capturing methods within a finite-difference discretization, whereas, on the same numerical grid, the Einstein elliptic equations are treated by resorting to spherical harmonics decomposition and solved, for each harmonic, by inverting band diagonal matrices. As a side product, we build and make available to the community a code to produce GRMHD axisymmetric equilibria for polytropic relativistic stars in the presence of differential rotation and a purely toroidal magnetic field. This uses the same XCFC metric solver of the main code and has been named XNS. Both XNS and the full X-ECHO codes are validated through several tests of astrophysical interest.Comment: 18 pages, 9 figures, accepted for publication in A&

Plasma Edge Kinetic-MHD Modeling in Tokamaks Using Kepler Workflow for Code Coupling, Data Management and Visualization

A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localized modes (ELMs) is presented in this report. This tool brings together, in a coordinated and effective manner, several first-principles physics simulation codes, stability analysis packages, and data processing and visualization tools. A Kepler workflow is used in order to carry out an edge plasma simulation that loosely couples the kinetic code, XGC0, with an ideal MHD linear stability analysis code, ELITE, and an extended MHD initial value code such as M3D or NIMROD. XGC0 includes the neoclassical ion-electron-neutral dynamics needed to simulate pedestal growth near the separatrix. The Kepler workflow processes the XGC0 simulation results into simple images that can be selected and displayed via the Dashboard, a monitoring tool implemented in AJAX allowing the scientist to track computational resources, examine running and archived jobs, and view key physics data, all within a standard Web browser. The XGC0 simulation is monitored for the conditions needed to trigger an ELM crash by periodically assessing the edge plasma pressure and current density profiles using the ELITE code. If an ELM crash is triggered, the Kepler workflow launches the M3D code on a moderate-size Opteron cluster to simulate the nonlinear ELM crash and to compute the relaxation of plasma profiles after the crash. This process is monitored through periodic outputs of plasma fluid quantities that are automatically visualized with AVS/Express and may be displayed on the Dashboard. Finally, the Kepler workflow archives all data outputs and processed images using HPSS, as well as provenance information about the software and hardware used to create the simulation. The complete process of preparing, executing and monitoring a coupled-code simulation of the edge pressure pedestal buildup and the ELM cycle using the Kepler scientific workflow system is described in this paper