Physical Effects of Aleph's Semi-implicit Particle-in-cell Algorithm
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Proceedings - International Symposium on Discharges and Electrical Insulation in Vacuum, ISDEIV
This presentation describes a new effort to better understand insulator flashover in high current, high voltage pulsed power systems. Both experimental and modeling investigations are described. Particular emphasis is put upon understand flashover that initiate in the anode triple junction (anode-vacuum-dielectric).
High-enthalpy hypersonic flight represents an application space of significant concern within the current national-security landscape. The hypersonic environment is characterized by high-speed compressible fluid mechanics and complex reacting flow physics, which may present both thermal and chemical nonequilibrium effects. We report on the results of a three-year LDRD effort, funded by the Engineering Sciences Research Foundation (ESRF) investment area, which has been focused on the development and deployment of new high-speed thermochemical diagnostics capabilities for measurements in the high-enthalpy hypersonic environment posed by Sandia's free-piston shock tunnel. The project has additionally sponsored model development efforts, which have added thermal nonequilibrium modeling capabilities to Sandia codes for subsequent design of many of our shock-tunnel experiments. We have cultivated high-speed, chemically specific, laser-diagnostic approaches that are uniquely co-located with Sandia's high-enthalpy hypersonic test facilities. These tools include picosecond and nanosecond coherent anti-Stokes Raman scattering at 100-kHz rates for time-resolved thermometry, including thermal nonequilibrium conditions, and 100-kHz planar laser-induced fluorescence of nitric oxide for chemically specific imaging and velocimetry. Key results from this LDRD project have been documented in a number of journal submissions and conference proceedings, which are cited here. The body of this report is, therefore, concise and summarizes the key results of the project. The reader is directed toward these reference materials and appendices for more detailed discussions of the project results and findings.
Abstract not provided.
Abstract not provided.
Physics of Fluids
For decades, it has been observed that the commonly used Borgnakke-Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke-Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke-Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke-Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). The consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.
Journal of Physics D: Applied Physics
Kinetic simulations of plasma phenomena during and after formation of the conductive plasma channel of a nanosecond pulse discharge are analyzed and compared to existing experimental measurements. Particle-in-cell with direct simulation Monte Carlo collisions (PIC-DSMC) modeling is used to analyze a discharge in helium at 200 Torr and 300 K over a 1 cm gap. The analysis focuses on physics that would not be reproduced by fluid models commonly used at this high number density and collisionality, specifically non-local and stochastic phenomena. Similar analysis could be used to improve the predictive capability of lower fidelity or reduced order models. First, the modeling results compare favorably with experimental measurements of electron number density, temperature, and 1D electron energy distribution function at the same conditions. Second, it is shown that the ionization wave propagates in a stochastic, stepwise manner, dependent on rare, random ionization events ahead of the ionization wave when the ionization fraction in front of the ionization wave is very low, analagous to the stochastic branching of streamers in 3D. Third, analysis shows high-energy runaway electrons accelerated in the cathode layer produce electron densities in the negative glow region over an order of magnitude above those in the positive column. Future work to develop reduced order models of these two phenomena would improve the accuracy of fluid plasma models without the cost of PIC-DSMC simulations.
Abstract not provided.
Abstract not provided.
IEEE International Conference on Plasma Science
Nanosecond pulsed discharges provide versatile experimental and computational testbeds for the exploration of fundamental plasma physics. In particular, the fast rise time and short duration produce plasmas which are both spatially diffuse and uniform enough to probe experimentally and confine the kinetics of interest to sufficiently short time scales to be computationally tractable. This work will focus on validation of particle-in-cell with Monte Carlo collisions (PIC-MCC) modeling and analysis of plasma phenomenon during and after formation of the conductive plasma channel of a nanosecond pulse discharge in helium at 200 Torr and 300 K over a 1 cm gap. The validation will compare results of the simulation to measurements of electron number density, temperature, 1D electron energy distribution function, and Townsend ionization coefficient, as well as ion mobility. Analysis of the stochastic nature of the electron avalanche ahead of the ionization wave front and of significant ionization overshoot in the presheath region is also performed.