Search NASASearch

SEARCH · Search NASA

Results for “high mass loading”

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.

773 records · Page 39

Influence of functional additives, fillers, and pigments on thermal and catalytic pyrolysis of polyethylene for waste plastic upcycling

Pyrolysis offers a relatively green and economical method to convert waste plastics into valuable chemicals and fuels without the need for harmful solvents, toxic chemicals, or costly high-pressure reactors. Despite its popularity among chemical upcycling technologies, industrial adoption suffers from feedstock heterogeneity, low-quality products, and catalyst deactivation. Most plastics in our daily lives are formulated with functional additives, fillers, and colorants. These additives remaining in end-of-life waste streams increase feedstock heterogeneity, creating a challenging issue in recycling plastics. Still, the potential impacts of additives on the chemical upcycling of plastics have been poorly understood. In this study, polyethylene compounded with a range of widely used additives (antioxidants, stabilizers, pigments, fillers, slip agents, and flame retardants) was subjected to both thermal pyrolysis and catalytic pyrolysis in different catalyst-to-feedstock contact modes. It showed that many inorganic additives, such as talc, kaolin, CaCO 3 , TiO 2 , carbon black, and zinc stearate, facilitated polymer decomposition during pyrolysis, increasing light hydrocarbons while also promoting aromatic and carbon residue formation. Conversely, antioxidants and stabilizers inhibited depolymerization, favoring heavier hydrocarbons. During catalytic pyrolysis with HZSM-5 zeolite, additives strongly enhanced aromatic and catalytic coke formation, especially when there was direct contact between plastics and catalysts. Although certain additives seem beneficial in the short term by promoting polymer cracking and improving the selectivity of aromatics, the transport of the additives and their degradation products and increased carbon coking can contaminate products, deactivate or modify catalysts, and foul reactors. These findings address a critical knowledge gap in effectively converting waste plastics via a greener route.

42 ENGINEERING

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

Decorated Honeycomb Lattice and Successive Magnetic Transitions of the Delafossite Derivative K 4 Ni 5 Te 3 O 16

Here, we report the synthesis and characterization of the delafossite derivative K 4 Ni 5 Te 3 O 16 , which hosts a decorated honeycomb Ni 2+ lattice. Single-crystal X-ray diffraction, electron microscopy, and neutron diffraction establish a noncentrosymmetric Pmm2 crystal structure that combines regular honeycomb units with elongated motifs formed by Ni 2+ trimers along extended edges. Magnetic susceptibility reveals strong antiferromagnetic interactions, while heat capacity measurements identify two successive magnetic transitions at 30 and 12 K. Neutron diffraction shows that these transitions correspond to a progression from partial to full long-range magnetic order on distinct Ni sites, reflecting competition among the various exchange interactions. The resulting magnetic point group, mm2.1′, inherits the lattice polarity and permits high-order magnetoelectric coupling. Consistent with this symmetry, a quadratic magnetic-field-induced polarization is experimentally observed below the magnetic ordering temperature, likely arising from exchange-driven magnetostriction coupled to the polar lattice. These results establish K 4 Ni 5 Te 3 O 16 as an interesting platform for engineering B-site ordering in delafossite derivatives to realize complex quantum magnetism and functional responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang

3D-Printed mmWave Quasi-Holographic Antenna for 2-D Beamforming

This paper presents the design and manufacturing of a novel 2D-scanning antenna that integrates a 3D-printed Rotman lens with a quasi-holographic leaky-wave antenna (HLWA). The proposed design achieved beam-scanning capabilities by leveraging the beamforming of the Rotman lens and the high-gain directive properties of the quasi-HLWA. The Rotman lens (RL) enables beam steering in the elevation plane by switching between input ports. The quasi-HLWA, designed using holographic principles, achieves frequency-controlled beam scanning in the azimuth plane. The entire antenna structure was fabricated using additive manufacturing with an Ink1092 substrate and silver ink for the conductive traces. This approach provides greater control over material placement and design freedom compared to traditional methods. A 25° transmission linear substrate taper was used to ensure good impedance matching between the Rotman lens and the quasi-HLWA, allowing greater gain while maintaining a good scanning range. The experimental results validate the 2D scanning capability of the proposed antenna. The antenna system provides coverage from −54° to 54° in the elevation θ plane and −28° to 28° in the azimuth plane ϕ , with a maximum measured gain of 21.3 dBi at 28 GHz with an average radiation efficiency η=60 %. The fabricated prototype was tested, and the performance was in good agreement with the simulated performance.

Rotman lens

Well-Defined Iron Sites in Crystalline Carbon Nitride

Carbon nitride materials can be hosts for transition metal sites, but Mössbauer studies on iron complexes in carbon nitrides have always shown a mixture of environments and oxidation states. Here we describe the synthesis and characterization of a crystalline carbon nitride with stoichiometric iron sites that all have the same environment. The material (formula C 6 N 9 H 2 Fe 0.4 Li 1.2 Cl, abbreviated PTI/FeCl 2 ) is derived from reacting poly(triazine imide)·LiCl (PTI/LiCl) with a low-melting FeCl 2 /KCl flux, followed by anaerobic rinsing with methanol. X-ray diffraction, X-ray absorption and Mössbauer spectroscopies, and SQUID magnetometry indicate that there are tetrahedral high-spin iron(II) sites throughout the material, all having the same geometry. As a result, the material is active for electrocatalytic nitrate reduction to ammonia, with a production rate of ca. 0.1 mmol cm –2 h –1 and Faradaic efficiency of ca. 80% at −0.80 V vs RHE.

Anions

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Mesoscale Magnetostructural Phase Separation in Fe‐deficient Fe 5 GeTe 2

Two-dimensional van der Waals ferromagnet Fe 5-x GeTe 2 (F5GT) is promising for spintronic applications due to its high Curie temperature, layered structure, and ability to host complex magnetic textures. However, the origin of its sample-dependent magnetic anisotropy remains unclear, hindering control of its magnetic behavior. Here, we use spatially resolved cryogenic scanning transmission electron microscopy (STEM) to correlatively map magnetism, lattice structure, and chemistry across atomic-to-micron scales. We reveal that only mesoscale, not nanoscale, inclusions of a Fe-deficient secondary phase significantly modify magnetic behavior, establishing a previously unrecognized critical length scale. This phase separation, induced by quenching, leads to in-plane magnetic anisotropy, while slow cooling confines separation to a few nanometers and preserves out-of-plane anisotropy. These findings reconcile prior inconsistencies and establish a predictive framework for tuning magnetism in F5GT through thermal processing, with broader implications for controlling anisotropy in other two-dimensional magnetic materials.

2D ferromagnets

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

Biodegradable Biocomposite Development for Large Scale Additive Manufacturing

Oak Ridge National Laboratory (ORNL) worked with Mitsubishi Chemical America, Inc. to develop a highly biodegradable ~100% biobased feedstock for large format additive manufacturing (AM). Together, they developed an AM feedstock based on the biobased polyester, polybutylene succinate (PBS), that is typically used in agriculture and packaging projects. The motivation for this work was to leverage the high degradability of this polymer in AM to create a new avenue for end-of-life disposal of the product. A composite feedstock was developed from PBS, polylactic acid (PLA), and wood flour (WF). While WF improves printability, reduces feedstock cost, and can improve stiffness, it also often results in a loss in tensile strength which can inhibit/limit applications of the material. To recover the loss in tensile strength suffered by including large quantities, up to 30 wt.% WF in the composite, a chain extender was used to compatibilize the different composite constituents with each other. The feedstock was used to 3D-print a demonstration mold for precast concrete catch basins and was tested by casting 750 lbs of concrete. The mold survived the demonstration. Aerobic degradation testing of the composite feedstock showed that it was highly degradable, 40% biodegradation after 90 days, compared to other incumbent biobased feedstocks.

100% biobased feedstock for large format additive

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium

Rapid Quantum Ground State Preparation via Dissipative Dynamics

Inspired by natural cooling processes, dissipation has become a promising approach for preparing low-energy states of quantum systems. However, the potential of dissipative protocols remains unclear beyond certain commuting Hamiltonians. This work provides significant analytical and numerical insights into the power of dissipation for preparing the ground state of noncommuting Hamiltonians. For quasi-free dissipative dynamics, including certain 1D spin systems with boundary dissipation, our results reveal a new connection between the mixing time in trace distance and the spectral properties of a non-Hermitian Hamiltonian, leading to an explicit and sharp bound on the mixing time that scales polynomially with system size. For more general spin systems, we develop a tensor network-based algorithm for constructing the Lindblad jump operator and for simulating the dynamics. Using this algorithm, we demonstrate numerically that dissipative ground state preparation protocols can achieve rapid mixing for certain 1D local Hamiltonians under bulk dissipation, with a mixing time that scales logarithmically with the system size. We then prove the rapid mixing result for certain weakly interacting spin and fermionic systems in arbitrary dimensions, extending recent results for high-temperature quantum Gibbs samplers to the zero-temperature regime. Together, these results show that dissipation can be a powerful tool for ground state preparation, with potential applications across condensed matter physics, quantum materials science, and beyond.

decoherence

Electron-Deficient Phenazines: Synthesis, Structures, Redox, and Optoelectronic Properties in the Gas Phase, in Solution, and in the Solid State

High-temperature gas-phase reactions of 1,4-C4F8I2 with phenazine (PHNZ) afforded new thermally- and air-stable derivatives with cyclo-C4F8 substituents, viz. 1,2-PHNZ(C4F8), 2,3-PHNZ(C4F8), 1,2:6,7-PHNZ(C4F8)2, 1,2:7,8-PHNZ(C4F8)2, and 1,2:8,9-PHNZ(C4F8)2. Reactions in the presence of Cu powder resulted in reductive-defluorination/aromatization of some of the cyclo-C4F8 substituents to produce 2,3-PHNZ(C4F4), 2,3-PHNZ(C4HF3), and 1,2:6,7-PHNZ(C4F8)(C4F4), which are formally tri- or tetrafluoro derivatives of benzo[a]phenazine. Single-crystal X-ray diffraction structures of these compounds demonstrate the planarity of 2,3-PHNZ(C4F4) and 2,3-PHNZ(C4HF3) and ordered p-stacking in 1,2-PHNZ(C4HF3), 1,2-PHNZ(C4F4), 1,2-PHNZ(C4F8), and 2,3-PHNZ(C4F8). Low-temperature gas-phase photoelectron spectra of the PHNZ(C4F8)1,2 compounds show that they are strong electron acceptors, with the three PHNZ(C4F8)2 isomers having electron affinities exceeding that of C60. Cyclic voltammetry in acetonitrile show that the five PHNZ(C4F8)1,2 compounds exhibit reversible one-electron reductions with large positive shifts in their reduction potentials relative to unsubstituted PHNZ. Spectroelectrochemical absorption and EPR spectra of PHNZ(C4F8)1,2- • anion radicals and extensive DFT calculations revealed the significant role that the different substitution patterns have on determining the optoelectronic properties of the PHNZ(C4F8)1,2 derivatives and PHNZ(C4F8)1,2- • anion radicals. Time-resolved microwave conductivity measurements and photoluminescence spectra of blends of the PHNZ(C4F8)1,2 derivatives with the donor poly(3-hexylthiophene) (P3HT) were recorded to assess the possible use of the new compounds in optoelectronic devices.

Balser, Sebastian

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID

Maintenance strategy, structural design, and site layout of the ST-E1 fusion power plant

An effective fusion reactor maintenance scheme enables safe operations and short downtimes. This in turn leads to high availability, which is critical to the commercial viability of a power-producing plant. In tokamak-based fusion power plants, the chosen maintenance approach has a significant impact on the spatial design of the tokamak, as well as the surrounding infrastructure, and therefore needs to be considered from the outset. Tokamak Energy has developed a pre-concept design of a fusion power plant, ST-E1. This work describes the major drivers and constraints that have been considered, presents the tokamak architecture and chosen maintenance regime, and discusses how this enables the plant’s two-phased approach to demonstrating commercial operations. It also shows the implications for the design of other systems areas, in particular the machine structural arrangement and bioshield and hot cell layout. The reactor core segmentation and removal scheme replaces entire toroidal segments radially through a large vacuum port, along a single axis only. The result is a change-tolerant machine and plant layout that can accommodate the evolving designs of the tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

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