Search NASASearch

SEARCH · Search NASA

Results for “Climate and Earth system modelling”

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.

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

1,694 records · Page 38

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE

Integrated optical frequency division for microwave and mmWave generation

Abstract The generation of ultra-low-noise microwave and mmWave in miniaturized, chip-based platforms can transform communication, radar and sensing systems 1–3 . Optical frequency division that leverages optical references and optical frequency combs has emerged as a powerful technique to generate microwaves with superior spectral purity than any other approaches 4–7 . Here we demonstrate a miniaturized optical frequency division system that can potentially transfer the approach to a complementary metal-oxide-semiconductor-compatible integrated photonic platform. Phase stability is provided by a large mode volume, planar-waveguide-based optical reference coil cavity 8,9 and is divided down from optical to mmWave frequency by using soliton microcombs generated in a waveguide-coupled microresonator 10–12 . Besides achieving record-low phase noise for integrated photonic mmWave oscillators, these devices can be heterogeneously integrated with semiconductor lasers, amplifiers and photodiodes, holding the potential of large-volume, low-cost manufacturing for fundamental and mass-market applications 13 .

Science & Technology - Other Topics

Shear and combustion characterization of printable ammonium perchlorate composite propellant formulations

Seven formulations of ammonium perchlorate composite propellant (APCP) are developed and their properties relevant for successful additive manufacture are characterized. Extrusion in a custom-built 3D printing system and spindle viscometry are used to collect viscosity measurements of curing primary and non-curing secondary versions of the formulations. The formulations that behave similarly to a Bingham plastic, with apparent viscosities between 4 and 8 thousand Pa*s at 30 minutes post-mix, are determined to be most suitable for printing applications. Measurements of one material show a yield stress of 1 kPa. Ambient pressure burn rates of 2.0–2.9 mm/s were measured for the compositions tested. All measured burn rates were comparable to the lower end of typical burn rates for APCP. In conclusion, the results of the characterization demonstrate a propellant suitable for use in a 3D printing system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Asymmetric pathways for lithium extraction and recovery based on the two-phase equilibrium of layered oxides

Electrochemical intercalation offers a promising platform for Li + extraction. However, only limited types of electrode materials have been investigated. The challenge to broaden and tailor materials for electrochemical intercalation-based Li + extraction lies in the lack of understanding of material’s response upon co-intercalation of multiple ions, therefore, paired process design to enable reversible Li + extraction and recovery. Here, we showcase the design of asymmetric ion pathways for Li + extraction and recovery for host material with complex Li + and Na + interaction using layered cobalt oxide as a model material. The two-phase equilibrium of Na 0.48 CoO 2 and Li 0.94 CoO 2 governs Li + selectivity when a high depth of intercalation is achieved (low vacancy level). We show that the relative rate between ion exchange and intercalation is critical to determine the ion pathways. The relationship can be quantitatively compared using the average pseudo ion exchange rate (C pseudoIX ) and the intercalation rate (C inter ). The ion pathways at the three regimes with C pseudoIX > C inter , C pseudoIX ~ C inter , and C pseudoIX < C inter are constructed. By selecting the optimized ion pathway and particle size, we demonstrate 9.7×10 4 Li + selectivity with 99% purity Li + recovery from an initial 1:1000 Li: Na molar ratio solution using 115 mAh/g specific capacity.

electrochemistry

Optoelectronic polymer memristors with dynamic control for power-efficient in-sensor edge computing

Abstract As the demand for edge platforms in artificial intelligence increases, including mobile devices and security applications, the surge in data influx into edge devices often triggers interference and suboptimal decision-making. There is a pressing need for solutions emphasizing low power consumption and cost-effectiveness. In-sensor computing systems employing memristors face challenges in optimizing energy efficiency and streamlining manufacturing due to the necessity for multiple physical processing components. Here, we introduce low-power organic optoelectronic memristors with synergistic optical and mV-level electrical tunable operation for a dynamic “control-on-demand” architecture. Integrating signal sensing, featuring, and processing within the same memristors enables the realization of each in-sensor analogue reservoir computing module, and minimizes circuit integration complexity. The system achieves 97.15% fingerprint recognition accuracy while maintaining a minimal reservoir size and ultra-low energy consumption. Furthermore, we leverage wafer-scale solution techniques and flexible substrates for optimal memristor fabrication. By centralizing core functionalities on the same in-sensor platform, we propose a resilient and adaptable framework for energy-efficient and economical edge computing.

Optics

Conceptual design study of neutron detectors for safeguards measurement of an irradiated pebble

Nuclear material control and accounting (MC&A) of pebble-bed reactors (PBRs) is challenging because a PBR utilizes hundreds of thousands of identical, unmarked pebbles that are continuously recirculated through the core. To develop tools that enable the implementation of international safeguards, especially in the context of MC&A of spent pebbles, we designed and simulated three neutron detection concepts to determine fissile content in individual pebbles: a differential die-away (DDA) detector, a californium interrogation prompt neutron (CIPN) detector, and a passive neutron albedo reactivity (PNAR) detector using Monte Carlo calculations. Burnup calculations were performed on the spent pebbles from the PBMR-400 classic PBR. The varying neutron and gamma source terms, and isotopic compositions in the spent pebbles calculated at various burnup levels were used in the neutron detector models. DDA was found to be sensitive to the number of passes a pebble has had through the core and to the fissile content contained in a spent pebble. Optimization in the DDA design further increased the neutron count rates and thus reduced counting uncertainty. Meanwhile, passive neutron counting using the same detector body could distinguish pebbles with different numbers of passes, but its response was dominated by neutron-emitting actinides and was not sensitive to fissile content. On the other hand, the PNAR technique was not viable for a single pebble but performed reasonably for a 27-pebble array, which suggested potential use for verification measurements of containers filled with 27 or more spent pebbles.

CIPN

Excitons in van der Waals magnetic materials

Two-dimensional magnetic semiconductors provide a unique platform where long-range magnetic order coexists with strongly bound excitons. Because excitonic states and magnetic moments originate from the same electronic orbitals and couple via intrinsic exchange interactions, optical excitations in these systems exhibit pronounced sensitivity to magnetic order. Recent experiments show unusually strong magneto-optical responses and direct exciton–magnon coupling, establishing new routes for controlling light–matter interactions with spin degrees of freedom. Here, this Review surveys key developments, focusing on representative material systems, experimental signatures, and theoretical frameworks used to describe these phenomena. We conclude with perspectives on how this rapidly evolving field could enable next-generation optoelectronic and quantum technologies leveraging the coupled dynamics of light, charge and spin.

36 MATERIALS SCIENCE

Nuclear Safety [Vol. 34, No. 1, January-March 1993]

Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 1 An Integrated Regional Approach To Risk Management of Industrial Systems, S. Chakraborty and R. Stratton; 9 Annual Technical Meeting of the NRC Cooperative Severe Accident Research Program, E. G. Silver; ACCIDENT ANALYSIS: 13 Deterministic Severe Accident Criteria as Severe Accident Design Criteria and Policy for the New Production Reactor-Heavy Water Reactor, P. T. Rhoads; 20 Proposed Deterministic Severe Accident Criteria for the Heavy Water Reactor-New Production Reactor Containment, K. D. Bergeron, S. E. Slezak, and C. E. Leach; 33 Quasi-Static Core Liquid Level Depression and Long-Term Core Uncovery During a PWR LOCA, Y. Kukita, R. R. Schultz, H. Nakamura, and J. Katayama; 49 LIRA: An Advanced Containment System to Minimize the Accidental Radioactivity Releases, A. Turricchia; 63 Errata to “A Review of Hydrogen Production During Melt/Water Interaction in LWRs,” Vol. 33, No. 4; CONTROL AND INSTRUMENTATION: 64 The Nuplex 80+™ Advanced Control Complex from ABB Combustion Engineering, F. Ridolfo, D. Harmon, and K. Scarola; DESIGN FEATURES: 76 A Scheme for Passive Isolation of the Containment of a Reactor, A. K. Ghosh, V. V. Raj, and K. Kakodkar; 84 Utility Requirements for Safety in the Passive Advanced Light-Water Reactor, T. U. Marston, W. H. Layman, and G. Bockhold, Jr.; WASTE AND SPENT FUEL MANAGEMENT: 91 Activities Related to Waste and Spent Fuel Management, Compiled by M. D. Muhlheim and E. G. Silver; OPERATING EXPERIENCES: 103 Testing Deficiencies in Auxiliary Systems Feedwater Systems, J. D. Kueck; 110 Reactor Shutdown Experience, Compiled by J. W. Cletcher; 113 Selected Safety-Related Events, Compiled by G. A. Murphy; 115 Operating U.S. Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver; RECENT DEVELOPMENTS: 133 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 145 Reports, Standards, and Safety Guides, D. S. Queener; 150 Proposed Rule Changes as of Sept. 30, 1992; ANNOUNCEMENTS: 75 Harvard Short Course on Risk: Science, Assessment, and Management; 102 International Conference on Reactor Physics and Reactor Computations; 102 The Fourth International Conference on Nuclear Waste Reprocessing and Waste Management [RECOD ’94]; 132 International Symposium on the Radiological Exposure of the Population of the European Community from Radioactivity in the Mediterranean Sea [MARINA-MED]; 162 International High-Level Radioactive Waste Conference ‘Technology Integration—Above and Below”; 154 The Authors; 158 Indexes to Nuclear Safety, Volume 33.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Unlocking superplasticity in medium and high-entropy alloys

Superplasticity, the capacity of materials to sustain extraordinary tensile elongations at elevated temperatures, underpins a range of advanced metal-forming technologies. Conventionally, it is achieved in fine-grained alloys where deformation is dominated by grain-boundary sliding, accommodated by diffusion and dislocation activity. The advent of medium- and high-entropy alloys (M/HEAs), with their high chemical complexity and unconventional phase stability, offers new pathways to superplastic behavior beyond traditional alloy systems. Although investigated only recently, several M/HEAs already exhibit elongations that rival or exceed those of classical superplastic materials, particularly when ultrafine or metastable microstructures are engineered. Here, we review progress in understanding superplastic deformation in M/HEAs, emphasizing the interplay among composition, initial microstructure, thermomechanical processing, and microstructural evolution during high-temperature deformation. We discuss approaches to generating the fine-grained structures necessary for grain-boundary sliding, including severe plastic deformation and tailored heat treatments. We further highlight dynamic phenomena such as phase transformations, evolving grain-boundary chemistry, and deformation-induced grain refinement that can enhance plasticity in these systems. These mechanisms often shift the balance of deformation processes, enabling large elongations even outside classical criteria. Finally, we outline key challenges for application, including cost, scalability, recyclability, and microstructural stability.

klenam, Desmond [University of the Witwatersrand,

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (

Design Strategies Based on Electronic Interactions for Effective Catalysts in Lithium–Sulfur Batteries

Abstract Lithium–sulfur batteries (LSBs) are considered promising next‐generation batteries due to their high energy density (>500 W h kg −1 ). However, LSBs exhibit an unsatisfactory energy density (<400 W h kg −1 ) and cycle life (<300 cycles) because of the shuttle effect caused by soluble lithium polysulfide (LiPS) intermediates and the sluggish conversion reaction kinetics caused by insulating sulfur (S 8 ) and lithium sulfide (Li 2 S). Although various types of catalysts, including metal‐based compounds to single‐atom catalysts, have been reported to address these issues, most catalysts exhibited limited catalytic activity under practical lean electrolyte conditions (<5 µL mg −1 ). A comprehensive understanding of the synthetic strategy and catalytic mechanism of catalysts is essential for their design, but understanding the electronic effects of the catalysts and LiPS is more important. Furthermore, the electronic design of these catalysts is not well understood. In this review, we introduce the catalytic mechanisms in LSBs and discuss catalyst design strategies in terms of electronic effects on the interactions between reactants and catalysts, with a primary focus on heterogeneous catalytic systems. We additionally consider how the electronic property of homogeneous systems, particularly redox mediators, affects catalytic behavior under lean electrolyte conditions and propose future research directions for catalyst development in LSBs.

Chemistry

Queued Up: 2026 Edition, Characteristics of Power Plants Seeking Transmission Interconnection As of the End of 2025 [Slides]

Electric transmission system operators (ISOs, RTOs, or utilities) require proposed power plants seeking to connect to the transmission grid to undergo a series of impact studies before they can be built. This process establishes what new transmission equipment or upgrades may be needed before a project can connect to the system and assigns the costs of that equipment. The lists of projects in this process are known as “interconnection queues”. In collaboration with https://www.interconnection.fyi. Berkeley Lab compiled, aggregated, and cleaned interconnection queue data from >50 transmission grid operators (7 ISO/RTOs and 50 non-ISO balancing areas), which collectively represent ~98% of currently installed U.S. electric generating capacity.

24 POWER TRANSMISSION AND DISTRIBUTION

Tuning Anisotropic Optical Properties of Inorganic and Hybrid Organic–Inorganic MXenes via Topochemical Surface Modification

Surface groups are central to the properties of MXenes, yet their role in optical anisotropy remains largely unexplored. Here, we use a topochemical route to synthesize single crystals of stacked Ti 3 C 2 Cl 2 and hybrid organic–inorganic MXenes (h-MXenes) with lateral sizes of 38–75 μm, rotational registry, and tunable interlayer spacing. Solid-state NMR spectroscopy shows that topochemical substitution generates mixed amido, imido, and hydride surface motifs, which modify the electronic structure of the Ti 3 C 2 inorganic core. Imaging spectroscopic ellipsometry with micron-scale spatial resolution enables reconstruction of the complex dielectric tensor of individual multilayer crystals. Ti 3 C 2 Cl 2 exhibits a type-II hyperbolicity above 930 nm, whereas h-MXenes do not display hyperbolicity within the measured 300–1700 nm window, instead showing reduced in-plane conductivity, suppressed out-of-plane light absorption, and a chain-length-dependent blue shift of a near-infrared absorption feature. These results demonstrate topochemical surface modification as a direct handle for engineering MXenes as surface-programmable optical media.

Hybrid materials

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato

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