HYPERS simulations of solar wind interactions with the Earth's magnetosphere and the Moon

The hybrid simulations, where the ions are treated kinetically and the electrons as a fluid, seek to describe ion microphysics with maximum physical fidelity. The hybrid approach addresses the fundamental need for space plasma models to incorporate physics beyond magnetohydrodynamics. Global hybrid simulations must account for a wide range of both kinetic ion and whistler/Alfvén wave spatio-temporal scales in strongly inhomogeneous plasmas. We present results from two three-dimensional hybrid simulations performed with a novel asynchronous code, HYPERS designed to overcome computational bottlenecks that typically arise in such multiscale simulations. First, we demonstrate an excellent match between simulated lunar wake profiles and observations. We also compare our simulations with two other simulations performed with conventional (time-stepped) hybrid codes. Second, we investigate the interaction of the solar wind with the Earth's dayside magnetosphere under conditions when the orientation of the interplanetary magnetic field is quasi-radial. In this high-resolution simulation we highlight three-dimensional properties of foreshock perturbations formed by the backstreaming ions

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

The physics of streamer discharge phenomena

In this review we describe a transient type of gas discharge which is commonly called a streamer discharge, as well as a few related phenomena in pulsed discharges. Streamers are propagating ionization fronts with self-organized field enhancement at their tips that can appear in gases at (or close to) atmospheric pressure. They are the precursors of other discharges like sparks and lightning, but they also occur in for example corona reactors or plasma jets which are used for a variety of plasma chemical purposes. When enough space is available, streamers can also form at much lower pressures, like in the case of sprite discharges high up in the atmosphere. We explain the structure and basic underlying physics of streamer discharges, and how they scale with gas density. We discuss the chemistry and applications of streamers, and describe their two main stages in detail: inception and propagation. We also look at some other topics, like interaction with flow and heat, related pulsed discharges, and electron runaway and high energy radiation. Finally, we discuss streamer simulations and diagnostics in quite some detail. This review is written with two purposes in mind: First, we describe recent results on the physics of streamer discharges, with a focus on the work performed in our groups. We also describe recent developments in diagnostics and simulations of streamers. Second, we provide background information on the above-mentioned aspects of streamers. This review can therefore be used as a tutorial by researchers starting to work in the field of streamer physics.Comment: 89 pages, 29 figure

Experimente zur Strukturbildung in Plasmen

Since the early days of plasma physics, it is known that plasma is more than a hotgaseous state of matter. The electric charge of the plasma constituents (usually electrons and ions) allows for collective behavior and the time scales of electron and iondynamics are well separated. Thus, plasmas exhibit a rich variety of phenomena whichare unknown to ordinary gases and contribute to the compelling properties of thisstate of matter. The wealth of phenomena, their relevance for astrophysical and technological applications, and their impact on other fields have driven plasma researchfor almost a century now and have provided much insight into the formation of staticstructures, oscillatory instabilities, waves, and turbulence

Spin polarized proton beam generation from gas-jet targets by intense laser pulses

A method of generating spin polarized proton beams from a gas jet by using a multi-petawatt laser is put forward. With currently available techniques of producing pre-polarized monatomic gases from photodissociated hydrogen halide molecules and petawatt lasers, proton beams with energy ~ 50 MeV and ~ 80 % polarization are proved to be obtained. Two-stage acceleration and spin dynamics of protons are investigated theoretically and by means of fully self-consistent three dimensional particle-in-cell simulations. Our results predict the dependence of the beam polarization on the intensity of the driving laser pulse. Generation of bright energetic polarized proton beams would open a domain of polarization studies with laser driven accelerators, and have potential application to enable effective detection in explorations of quantum chromodynamics

Computational Electromagnetic Methods for Transcranial Magnetic Stimulation.

Transcranial magnetic stimulation (TMS) is a noninvasive technique used both as a research tool for cognitive neuroscience and as a FDA approved treatment for depression. During TMS, coils positioned near the scalp generate electric fields and activate targeted brain regions. In this thesis, several computational electromagnetics methods that improve the analysis, design, and uncertainty quantification of TMS systems were developed. Analysis: A new fast direct technique for solving the large and sparse linear system of equations (LSEs) arising from the finite difference (FD) discretization of Maxwell’s quasi-static equations was developed. Following a factorization step, the solver permits computation of TMS fields inside realistic brain models in seconds, allowing for patient-specific real-time usage during TMS. The solver is an alternative to iterative methods for solving FD LSEs, often requiring run-times of minutes. A new integral equation (IE) method for analyzing TMS fields was developed. The human head is highly-heterogeneous and characterized by high-relative permittivities (10^7). IE techniques for analyzing electromagnetic interactions with such media suffer from high-contrast and low-frequency breakdowns. The novel high-permittivity and low-frequency stable internally combined volume-surface IE method developed. The method not only applies to the analysis of high-permittivity objects, but it is also the first IE tool that is stable when analyzing highly-inhomogeneous negative permittivity plasmas. Design: TMS applications call for electric fields to be sharply focused on regions that lie deep inside the brain. Unfortunately, fields generated by present-day Figure-8 coils stimulate relatively large regions near the brain surface. An optimization method for designing single feed TMS coil-arrays capable of producing more localized and deeper stimulation was developed. Results show that the coil-arrays stimulate 2.4 cm into the head while stimulating 3.0 times less volume than Figure-8 coils. Uncertainty quantification (UQ): The location/volume/depth of the stimulated region during TMS is often strongly affected by variability in the position and orientation of TMS coils, as well as anatomical differences between patients. A surrogate model-assisted UQ framework was developed and used to statistically characterize TMS depression therapy. The framework identifies key parameters that strongly affect TMS fields, and partially explains variations in TMS treatment responses.PhDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/111459/1/luisgo_1.pd