SEARCH · Search NASA

Results for “in situ XAFS spectroscopy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

2,765 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Investigating the Combustion Performance of Dual Fuel Combustion with Diesel and Port Injected Hydrogen in a Large Bore Locomotive Engine

The heavy-duty transportation sector has primarily relied on conventional diesel combustion engines given their reliability and high thermal efficiency relative to spark ignition engines, but increased focus on reducing greenhouse gas emissions has led to investigation into alternative fuels. Gaseous hydrogen fuel has garnered a great deal of recent interest in the engine community given it has zero carbon, but hydrogen is not available at the scale and cost that petroleum fuels are currently available, and this is a barrier to adoption for industries that are looking to decarbonize their operations. Because of the fuel flexibility provided, dual fuel technology offers a pathway for some industries to adopt hydrogen as a fuel source while maintaining sufficient flexibility in times and locations where the new fuel is not yet available. This computational study investigates dual fuel combustion in a large bore locomotive engine architecture using direct injected diesel and port injected gaseous hydrogen fuel. With an optimal port fuel injection configuration from previous work, simulations of varying substitution ratio, compression ratio, manifold air temperature, diesel injection timing, and diesel injection pressure were performed to understand their effect on combustion performance. Results indicated that both increased substitution ratio and higher intake air temperature accelerates hydrogen flame propagation and can result in high peak cylinder pressures. Additionally, diesel injection timing and injection pressure were demonstrated as effective methods for controlling dual fuel combustion heat release rates.

ODonnell, Patrick Christopher↗

DOE OSTI2026

Final Technical Report: Transport of Complex Mixtures in Ion-Containing Polymer Membranes

Permselective ion-containing membranes are an integral component for many applications from water treatment, fuel cells, and solar fuels devices where the selective transport of molecules and ions is desired. In solar fuels devices, ion-containing polymer membranes are responsible for permitting selective transport of ions between electrodes to maintain overall charge neutrality yet limit transport of reaction products produced at the electrodes. While the transport of single solutes through such membranes has been fairly well described, binary and multicomponent transport is poorly understood due to the myriad of interactions that occur in these systems (i.e. between co-permeants and between permeants and the membrane). Solar fuels devices are just one example of an application where understanding the transport of multiple simultaneous species is critically important to improving device performance as product crossover leads to reductions in overall device performance. The objectives of this research was to improve our understanding of the complex array of factors that influence transport behavior of multiple solutes within ion-containing polymer membranes. This experimental project addressed the lack of fundamental understanding of multicomponent transport behavior by synthesizing ion exchange membranes with varied incorporation of comonomers (ionic and neutral moieties) to investigate fundamental relationships between membrane structure, membrane physiochemical properties, and transport behavior of solutes and complex solute mixtures through dense, hydrated membranes.

25 ENERGY STORAGE↗

DOE OSTI2026

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2026

Joint scheduling of energy, fast and primary frequency response reserves in integrated transmission–distribution networks

Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.

Noh, Seung-Gil↗

DOE OSTI2026

Universal Relationship between Limiting Current and Electrochemical Transport Properties in Malonate-Based Polymer Electrolytes

There is considerable interest in developing high-performance electrolytes for rechargeable lithium batteries. For practical applications, the electrolyte must support large dc currents. However, the parameters most often reported in the literature, conductivity, κ, and current fraction, ρ + , reflect ion transport in the limit of infinitesimal currents. In this limit, the efficacy of an electrolyte is given by the product κρ + . The limiting current density, i lim , is the maximum current density that can be applied across an electrolyte; the cell voltage diverges if the applied current density exceeds ilim. This parameter reflects ion transport in the limit of large dc currents and is therefore of practical interest. It would therefore be convenient if i lim could be predicted from measurements of κρ + . In order to explore this possibility, we studied six malonate-based polymers and PEO at a fixed salt concentration (r = 0.08) and temperature (90°C) using symmetric cells with planar electrodes. Unfortunately, there is no correlation between ilim and κρ + . When the applied current density, i, is less than ilim, the cell voltage approaches a stable plateau, ϕ plateau . Here, we found a linear dependence between i and thickness-normalized plateau potential, ϕ plateau L –1 , irrespective of the magnitude of the applied current. In all seven polymer electrolytes, we found a linear correlation between ilim and the slopes of these lines, σ. In other words, measurements of σ can be used to predict the limiting current.

Jana, Rounak [Lawrence Berkeley National Laborator↗

DOE OSTI2025

Interaction of magnetic fields with spinons in a fractionalized state

The 4d-electron trimer lattice Ba₄Nb₁₋ₓRu₃₊ₓO₁₂ exhibits either a quantum spin liquid (QSL) or a heavy-fermion strange metal (HFSM) phase, depending on Nb content. In the QSL state, itinerant spinons act as effective heat carriers, enhancing thermal conductivity. Strikingly, applying a magnetic field up to 14 T causes an abrupt, up-to-5000% increase in heat capacity below 150 mK, disrupting the linear temperature dependence typical of both phases. Meanwhile, AC susceptibility and electrical resistivity remain nearly unchanged, while thermal conductivity drops by up to 40% below 4 K. These results suggest spinons, despite being charge-neutral, are highly sensitive to magnetic fields at low temperatures. We propose that the magnetic field could induce Anderson localization of spinons, creating emergent non-magnetic two-level systems responsible for the Schottky-like anomaly in heat capacity. These findings point to a previously unexplored regime of spinon dynamics, potentially governed by field-induced localization and distinct from conventional magnetic or transport signatures.

Condensed-matter physics↗

DOE OSTI2024

O Corona, where art thou? eROSITA’s view of UV-optical-IR variability-selected massive black holes in low-mass galaxies

Finding massive black holes (MBHs,M BH ≈ 10 4 –10 7 M ⊙ ) in the nuclei of low-mass galaxies $\left( {{M_*}\mathop {\mathop < \limits_ }\limits_ {{10}^{10}}{M_ \odot }} \right)$ is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of ~200 UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available SRG/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. This result reveals the potential for severe biases in occupation fractions derived from data from only one waveband combined with SEDs and scaling relations of more massive black holes and galaxies.

Astronomy & Astrophysics↗

DOE OSTI2026

Generalizable Web User Interface for Scalable and Streamlined Deployment of Building Energy Management Systems in Small and Medium-Sized Commercial Buildings

Small and medium-sized commercial buildings (SMCBs) comprise 94% of US commercial buildings yet face significant barriers to implementing building energy management systems despite advances in smart device technology. Existing solutions present critical limitations: cloud-based API solutions simplify deployment but create vendor lock-in constraints; commercial integrated software solutions ensure compatibility via standardized protocols but require substantial cost and technical expertise; open-source IoT platforms offer cost-effective vendor independence but provide insufficient standardized protocol support for commercial building automation. This research presents a generalizable web user interface framework that bridges the gap between evolving smart device capabilities and lagging software infrastructure for SMCBs. The proposed system integrates VOLTTRON open-source middleware with an automated configuration converter that transforms unified specifications written in YAML, a human-readable data-serialization format, into system-specific files, streamlining manual setup processes. The vendor-agnostic architecture supports industry-standard protocols (BACnet and Modbus) and semantic building models while providing adaptive web interfaces that dynamically adjust to various building configurations. Demonstrations through simulation-based testing and a field deployment show automatic interface adaptation across heterogeneous HVAC systems and multizone monitoring. The automated configuration converter also substantially reduces labor-intensive setup.

Chung, Jihoon [ORNL] (ORCID:0000000184880815)↗

DOE OSTI2026

Giant Exfoliation Induced Magnetic Coercivity in Fe 3 GaTe 2

Permanent magnets with strong anisotropy and high coercivity underpin modern information and energy technologies, yet rare-earth-free alternatives remain limited. Here, we show that thickness engineering via mechanical exfoliation induces hard magnetic behavior in the van der Waals ferromagnet Fe 3 GaTe 2 . Bulk crystals exhibit Curie temperatures above 350 K but negligible room-temperature coercivity. When thinned below ∼100 nm, the coercive field is dramatically enhanced, reaching nearly 1 T at room temperature for in-plane fields—comparable to conventional hard magnets. Micromagnetic analysis reveals a crossover in magnetization reversal from domain-mediated processes in bulk samples to quasi-coherent rotation in thin flakes, driven by increased effective anisotropy and suppressed domain formation. This thickness-dependent transition enables tuning of magnetic hardness without chemical modification. Combined with high saturation magnetization and robust room-temperature performance, Fe 3 GaTe 2 emerges as a promising rare-earth-free material for spintronic applications. Its layered structure further allows integration into van der Waals heterostructures, where large in-plane coercivity can stabilize magnetic states against perturbations and interlayer coupling, offering potential for high-density nonvolatile memory and domain-wall-based devices.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Atomic faulting drives exceptional toughness in low thermal expansion chromium alloys

Endowing functional properties with mechanical responses in traditional metals has been a frontier topic, akin to transforming base metal into gold. Chromium and its alloys, with their functional deficiencies and limited ductility, serve as typical examples. Herein, we report a Cr96Fe4Ge1.3B1 alloy that unifies low thermal expansion (LTE, αl = 1.79 × 10-6 K-1, 200 − 315 K) with exceptional toughness (240.2 J·cm-3). The enhancement in mechanical responses is primarily attributed to layered Cr2B intermetallic precipitates, which ameliorate interfacial cohesion and simultaneously refine the grain structure. The weakened interlayer interactions within the Cr-B layers facilitate the nucleation and movement of numerous tiny stacking faults in precipitates, efficiently alleviating strain energy and resulting in marked work-hardening ability. Additionally, antiferromagnetic fluctuations in the BCC matrix contribute to the unique LTE behavior. This paves the way for the design of high-performance alloys featuring layered-symmetry precipitates.

Yu, Chengyi [University of Science and Technology ↗

DOE OSTI2026

Electron/hadron separation with the TRD based on GaAs pixel detectors

Transition Radiation Detectors (TRDs) are widely used for particle identification in both high-energy physics and astroparticle physics. They are typically equipped with gaseous detectors. The main limitation of these types of detectors is that TR photons and ionization losses cannot be decoupled, which significantly reduces the particle separation power. Recent advancements in the development of pixel detectors based on GaAs sensors offer a unique opportunity to effectively detect TR photons and separate them from ionization losses. Such detectors represent novel devices that combine precise tracking capabilities with particle identification (PID) properties. The present work is dedicated to an experimental study of particle identification properties of the TRD prototype based on 500 μ m-thick GaAs sensor bonded to a Timepix3 chip. Studies were performed at the CERN SPS and it was shown that at a particle momentum of 20 GeV/c, the probability of misidentifying a hadron as an electron is below 10 −2 for an electron detection efficiency of 98%–99%, and it is 1 . 6⋅10 −4 for electron detection efficiency of 90%. This performance surpasses the electron/hadron rejection power of all known TRDs by an order of magnitude for the same detector length.

GaAs detectors↗

NASA NTRS1998

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff↗

DOE OSTI2024

Massive red spiral galaxies in SDSS-IV MaNGA survey

ABSTRACT Massive red spiral galaxies (MRSGs) are supposed to be the possible progenitors of lenticular galaxies (S0s). We select a large sample of MRSGs ($M_*\gt 10^{10.5}\rm {\rm M}_{\odot }$) from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) DR17 using the g − r colour versus stellar mass diagram, along with control samples of blue spirals and S0s. Our main results are as follows: (1) After comparing the Sérsic index, concentration parameter, asymmetry parameter distribution, size–mass relation, and Σ1 (stellar mass surface density within the central 1 kpc)−mass relation, we find MRSGs are similar to S0s and have more compact and symmetric structures than blue spirals. MRSGs also resemble S0s in Dn4000, metallicity, Mgb/$\rm \left\langle Fe \right\rangle$, and V/σ radial profile. (2) By using MaNGA 2D spectra data, we separate the spatial regions into inner (R < 0.8Re) and outer (0.8 < R < 1.5Re) regions, and detect residual star formation in the outer regions of MRSGs. (3) When we select a sub-sample of MRSGs with NUV − r > 5, we find that they are completely star formation quenched in both inner and outer regions. Compared to optically selected MRSGs, NUV − r selected MRSGs appear to be more concentrated and have more massive dark matter haloes. The similarities between S0s and MRSGs suggest the possible evolutionary trend between MRSGs and S0s.

Astronomy & Astrophysics↗

DOE OSTI2024

Trust-Based Detection and Mitigation of Cyber Attacks in Distributed Cooperative Control of Islanded AC Microgrids

In this study, we address the challenge of detecting and mitigating cyber attacks in the distributed cooperative control of islanded AC microgrids, with a particular focus on detecting False Data Injection Attacks (FDIAs), a significant threat to the Smart Grid (SG). The SG integrates traditional power systems with communication networks, creating a complex system with numerous vulnerable links, making it a prime target for cyber attacks. These attacks can lead to the disclosure of private data, control network failures, and even blackouts. Unlike machine learning-based approaches that require extensive datasets and mathematical models dependent on accurate system modeling, our method is free from such dependencies. To enhance the microgrid’s resilience against these threats, we propose a resilient control algorithm by introducing a novel trustworthiness parameter into the traditional cooperative control algorithm. Our method evaluates the trustworthiness of distributed energy resources (DERs) based on their voltage measurements and exchanged information, using Kullback-Leibler (KL) divergence to dynamically adjust control actions. We validated our approach through simulations on both the IEEE-34 bus feeder system with eight DERs and a larger microgrid with twenty-two DERs. The results demonstrated a detection accuracy of around 100%, with millisecond range mitigation time, ensuring rapid system recovery. Additionally, our method improved system stability by up to almost 100% under attack scenarios, showcasing its effectiveness in promptly detecting attacks and maintaining system resilience. These findings highlight the potential of our approach to enhance the security and stability of microgrid systems in the face of cyber threats.

Computer Science↗

NASA NTRS1964

Performance of Solar Reflectors When Applied to Storage of Cryogenic Propellants in Space

An analysis was performed to determine the effectiveness of solar reflectors in reducing heat transfer through the insulation of cryogenic fuel tanks located in deep space. A solar reflector, as used herein, is defined as a surface having a ratio of absorptance with respect to solar radiation to hemispherical emittance (α/ϵ) less than 1.0. This analysis considers the effect of surface coatings with α/ϵ values varying from 0.1 to 1.0 when applied to the external surfaces of bulk insulation or multifoil insulation. Results were obtained in terms of heat leak to two cryogenic propellants, liquid hydrogen and liquid oxygen, when stored at Mars, Earth, and Mercury orbital distances from the sun. The results show that, with sufficiently good insulation and low heat leaks to the propellant, α/ϵ is a reasonable criterion of merit for these coatings, but that, at larger heat leaks, consideration of ϵ (or α) is necessary. It was observed that selective finishes of α/ϵ = 0.1 on a given amount of bulk insulation material can reduce the heat leak to 40 to 55 percent of the heat leak obtained when gray finishes are considered. Or, for equivalent heat leaks, a reduction in bulk insulation material thickness of about one-half is possible with high-quality selective finishes (α/ϵ = 0.1) compared with flat reflectors. Larger benefits may be achieved from selective finishes painted on foils; the analysis shows that heat leak may be reduced tenfold for surfaces with an α/ϵ of 0.1 compared with gray surfaces. Again, consideration of equivalent heat leaks indicates that the use of approximately 10 times fewer foils of emissivity 0.05 is possible when spectrally selective surfaces are employed. The present state of the art in achieving such surfaces is also briefly reviewed.

Propellant Storage↗

DOE OSTI2025

Fast lithium ion diffusion in brownmillerite Li x Sr 2 Co 2 O 5

Transition metal oxides not only exhibits novel magnetic properties but also provides outstanding ionic transports. Ionic conductors have great potential for interesting tunable physical properties via ionic liquid gating and novel energy storage applications such as all-solid-state lithium batteries. In particular, low migration barriers and high hopping attempt frequency are the keys to achieve fast ion diffusion in solids. Taking advantage of the oxygen-vacancy channel in Li x Sr 2 Co 2 O 5 , we show that migration barriers of lithium ion are as small as 0.28–0.17 eV depending on the lithium concentration rates. Our first-principles calculation also investigated hopping attempt frequency and concluded the room temperature ionic diffusivity and ion conductivity are high as 10 −7 –10 −6 cm 2 s −1 and 10 −3 –10 −2 Scm −1 , respectively, which outperform most of perovskite-type, garnet-type, and sulfide Li-ion solid-state electrolytes. This work proves Li x Sr 2 Co 2 O 5 as a promising super-ionic conductor.

Crystallographic defects↗

NASA NTRS2012

Returning an Entire Near-Earth Asteroid in Support of Human Exploration Beyond Low-Earth Orbit

This paper describes the results of a study into the feasibility of identifying, robotically capturing, and returning an entire Near-Earth Asteroid (NEA) to the vicinity of the Earth by the middle of the next decade. The feasibility of such an asteroid retrieval mission hinges on finding an overlap between the smallest NEAs that could be reasonably discovered and characterized and the largest NEAs that could be captured and transported in a reasonable flight time. This overlap appears to be centered on NEAs roughly 7 m in diameter corresponding to masses in the range of 250,000 kg to 1,000,000 kg. The study concluded that it would be possible to return a approx.500,000-kg NEA to high lunar orbit by around 2025. The feasibility is enabled by three key developments: the ability to discover and characterize an adequate number of sufficiently small near-Earth asteroids for capture and return; the ability to implement sufficiently powerful solar electric propulsion systems to enable transportation of the captured NEA; and the proposed human presence in cislunar space in the 2020s enabling exploration and exploitation of the returned NEA. Placing a 500-t asteroid in high lunar orbit would provide a unique, meaningful, and affordable destination for astronaut crews in the next decade. This disruptive capability would have a positive impact on a wide range of the nation's human space exploration interests. It would provide a high-value target in cislunar space that would require a human presence to take full advantage of this new resource. It would offer an affordable path to providing operational experience with astronauts working around and with a NEA that could feed forward to much longer duration human missions to larger NEAs in deep space. It represents a new synergy between robotic and human missions in which robotic spacecraft would retrieve significant quantities of valuable resources for exploitation by astronaut crews to enable human exploration farther out into the solar system. The capture, transportation, examination, and dissection of an entire NEA would provide valuable information for planetary defense activities that may someday have to deflect a much larger near-Earth object. Transportation of the NEA to lunar orbit with a total flight time of 6 to 10 years would be enabled by a ~40-kW solar electric propulsion system with a specific impulse of 3,000 s. The flight system could be launched to low-Earth orbit (LEO) on a single Atlas V-class launch vehicle, and return to lunar orbit a NEA with at least 28 times the mass launched to LEO. Longer flight times, higher power SEP systems, or a target asteroid in a particularly favorable orbit could increase the mass amplification factor from 28-to-1 to 70-to-1 or greater. The NASA GRC COMPASS team estimated the full life-cycle cost of an asteroid capture and return mission at approx.$2.6B.

asteroid mining↗

DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗