Search NASASearch

SEARCH · Search NASA

Results for “high energy density physics”

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,972 records · Page 45

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.

Aerodynamics

A Parametric Battery Model for the Conceptual Design of Electric Aircraft

Electric aircraft conceptual designs usually have an assumed specific power and specific energy for the aircraft’s energy storage system that dictates the total capacity and the peak power available. The specific energy and power are inherently linked to the state-of-charge and discharge rate from which they were derived; therefore, they will not directly correspond to the diverse circumstances encountered in various aircraft missions and flight segments. These parameters turn the battery system into a black box, disregard potential electrical restrictions, and disallow the aircraft and battery to be optimized as a system. Peering into this box, this study highlights the importance of incorporating a parametric battery model into the conceptual design workflow by splitting high-level terms such as power into voltage and current and investigating their variability during the discharge process. Through the modeling of these more detailed parameters, this methodology shows the feasibility of using low states-of-charge for contingency operations, including the end of the reserve mission, expanding the amount of usable capacity for electric aircraft. Stark differences in aircraft capabilities can arise between varying fidelity battery models due to late-mission, high-power flight operations. This parametric battery model effectively captures these differences by evaluating the limitations that arise within the individual battery cells and the aircraft powertrain. This paper shows that the unusable charge of a battery is set by the balked landing power requirement and can realistically range from 5% to 44% based on assumptions. This sets the analog to unusable fuel capacity in aircraft with liquid fuel systems. Needing only aircraft- and mission-level inputs and only seconds of run time, this model is a prime fit for the fast, accurate exploration of the electric aircraft conceptual design space.

Battery

Development of in-situ polymerized intrinsically conductive resin and low-cost carbon pigments offering high conductivity for sensing, EMI shielding and lighting protection

Electrically conductive composites are emerging across diverse industries such as electronic, automotive, aerospace, advanced air mobility, biomedical, infrastructure, defense and security offering static charge dissipation, electromagnetic interference shielding, lighting protection, sensing, dicing, corrosion monitoring, etc. Conductivity enhanced composites provide several advantages compared to conventional metals including weight reduction, corrosion resistance, energy efficient processability, tunable properties and multifunctionality. Polymers are typically insulating in nature and require conducting filler for electron transport. However, dispersion and polymer-filler interphases are critical and often disrupt conducting pathways. Besides, conductive fillers such as graphene, carbon nanotube, MXene, silver nanowire, etc. are expensive, limiting their wide adoption in composite industry. On the other hand, a limited number of intrinsically conductive polymers are available among which polyaniline (PANI) has been widely studied due to its high conductivity, thermal and chemical stability. However, PANI is difficult to process and exhibits weak mechanical properties. In brief, there is a significant demand for electrically conductive polymer formulation with cost-effective conducting fillers that offer processability in scale to expand the market of conductivity enhanced materials.

Kumar, Vipin [Oak Ridge National Laboratory (ORNL)

Experimental and analytical ion thruster research

The results of further spectroscopic studies on the plume from a 3 cm ion source operated on an argon propellant is reported on. In particular, it is shown that it should be possible to use the spectroscopic technique to measure the plasma density of the ion plume close to the grids, where it is difficult to use electrical probe measurements. How the technique, along with electrical probe measurements in the far downstream region of the plume, can be used to characterize the operation of a three-grid, 15 cm diameter thruster from NASA JPL is outlined. Pumping speed measurements on the Vacuum Research Facility have shown that this facility should be adequate for testing the JPL thruster at pressures in the low 10(exp -5) Torr range. Finally, we describe a simple analytical model which can be used to calculate the grid impingement current which results from charge-exchange collisions in the ion plume.

Wilhelmus M. Ruyten

Lunar and Planetary Science XXXV: Icy Worlds: Moving and Grooving

Reports from the conference session entitled Icy Worlds: Moving and Grooving, include:Mass Anomalies on Ganymede; Europan Chaos and Lenticulae: A Synthesis of Size, Spacing, and Areal Density Analyses; Thermal and Topographic Tests of Europa Chaos Formation Models; Flexure of Europa s Lithosphere Due to Ridge-Loading; Ridges on Europa: Origin by Incremental Ice-Wedging ; Convergent Boundaries on Europa: a Numerical Approach to Euler Pole Analysis and Its' Implications for Plate Reconstruction; Numerical Simulations of Subsolidus Convection in the Ice Shell of Europa: Implications for the Thermal Evolution and Present State; Effects of Plasticity on Convection in an Ice Shell: Implications for Europa; Non-Newtonian Convection and Compositional Buoyancy: Advances in Modeling Convection and Dome Formation on Europa; Convective Instability in Ice I: Application to Callisto and Ganymede; Crater Size Distributions on Callisto: A Galileo SSI Summary; Neutron Diffraction Studies of Planetary Ices; and H2O2 Synthesis Induced by Irradiation of H2O with Energetic H+ and Ar+ Ions at Various Temperatures.

Source record

Analytical modeling of battery cycle life

An analytical model related to the physical and chemical processes involved in battery wear and failure is discussed. The model is described using the data from the Crane tests.

Gross, S.

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION

Monitoring Depolymerization in Mesopores Using Dynamic Properties of Polymeric Melt Accessed via Dielectric Spectroscopy

Traditional design principles for heterogeneous catalysis guide the use of catalytic particles with mesosized (∼2–50 nm) pores to increase the number of surface-active sites by way of an increased surface area. However, the entry of long-chain polymers into such pores may be significantly limited by the size and entanglement of polymers in the melt state, thereby decreasing the number of accessible sites. Assessment of catalyst performance from traditional reactor-based studies averages over intrapore reaction events as well as reactions on the surface of a particle, resulting in an inability to distinguish between differences in site accessibility and activity. Techniques that assess the intrapore performance can inform the design of future heterogeneous catalysts for polymer upcycling. In this work, we demonstrate the use of broadband dielectric spectroscopy to monitor depolymerization of a polymer melt within mesopores via changes in the segmental relaxation time scale of amorphous polymer chains. In particular, we highlight the use of an anodic aluminum oxide (AAO) membrane as a readily available model for catalyst pores with a well-characterized pore morphology. The decrease in the segmental relaxation (α-relaxation) time of the melt with increasing chain scission emerges as a measure of the extent of polymer deconstruction inside mesopores. To demonstrate the utility of this technique, we demonstrate the decomposition of two commercial poly(propylene carbonate) polymers with different decomposition rates within mesopores. As the polymers depolymerize, their segmental relaxation time decreases as the molecular weight decreases (as predicted by the Fox–Flory equation). The BDS-measured change in segmental relaxation time mirrors the expected trend based on change in molecular weight measured by size exclusion chromatography.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Nuclear Safety [Vol. 33, No. 4, October-December 1992]

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: 499 Technical Note: Nuclear Regulatory Commission’s Examination Question Bank System, C. L Jensen and N. L. Skinner; 506 Report on the Specialists’ Meeting on Passive and Active Safety Features of Liquid-Metal Fast Breeder Reactors Organized by the International Atomic Energy Agency at Oarai Engineering Centre of Power Reactor and Nuclear Development Corporation, Japan, November 5-7,1991, S. R. Paranjpe; ACCIDENT ANALYSIS: 514 A Review of Hydrogen Production During Melt/Water Interaction in LWRs, D. F. Fletcher, B. D. Turland, and S. P. A. Lawrence; 535 A Review of the U 3 0 8 -Al Reaction as a Potential Heat Source in Research and Test Reactor Accidents, J. L. Snelgrove; CONTROL AND INSTRUMENTATION: 543 Multilevel Interfaces for Power Plant Control Rooms II: A Preliminary Design Space, K. J. Vicente; DESIGN FEATURES: 549 Twenty-Second DOE/NRC Nuclear Air Cleaning and Treatment Conference, R. R. Bellamy, D. W. Moeller, and M. W. First; 562 Book Review: Safety Assessment of Unit 5 (WWER-440/W-213) of the Greifswald Nuclear Power Station, M. H. Fontana; WASTE AND SPENT FUEL MANAGEMENT: 576 Activities Related to Waste and Spent Fuel Management, M. D. Muhlheim and E. G. Silver; OPERATING EXPERIENCES: 586 Analysis of Corrective Actions Applied to Nuclear Power Plant Operations, A. E. Cross; 593 Reactor Shutdown Experience, Compiled by J. W. Cletcher; 596 Operating U.S. Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver; 609 Non-Commercial Power Operating Experience; 614 World Nuclear Performance, Compiled by M. D. Muhlheim; RECENT DEVELOPMENTS: 623 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 635 Reports, Standards, and Safety Guides, D. S. Queener; 639 Proposed Rule Changes as of June 30, 1992; ANNOUNCEMENTS: 513 Fourth International Conference on Nuclear Fuel Reprocessing and Waste Management RECOD '94; 585 International Topical Meeting on Advanced Reactors Safety; 622 International Workshop on Implementation of ALARA at Nuclear Power Plants; 648 Short Course on Nuclear Criticality Safety; 648 8th Topical Meeting on Problems in Reactor Physics MEPhi-VOLGA-93; 643 The Authors; 646 Reviewers of Nuclear Safety, Vol. 33

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Pinned in Microgravity: How Colloidal Suspensions Inhibit Droplet Motion

Understanding the fundamental behavior of colloidal fluids in microgravity is critical for advancing material science, manufacturing, and biological processes. In this study we investigate the dynamics of a droplet containing colloidal suspension when ejected in short term microgravity. These conditions are generated using a 2.2 second drop tower. A custom droplet generator was designed to produce consistent size-controlled droplets containing passive colloidal particles. The particles are suspended in water and subsequently ejected out of a hydrophobic wedge under microgravity conditions, where the droplet ejection time can be predicted [Torres & Weislogel, 2021]. The ejection times will be compared with de-ionized water ejected from the same system to investigate the effects of change in density and surface tension caused by the colloids. Successful passive colloidal ejection will lay the foundation for active colloid research in the future.

colloids

A Battery Charger and State of Charge Indicator

A battery charger which has a full wave rectifier in series with a transformer isolated 20 kHz dc-dc converter with high frequency switches, which are programmed to actively shape the input dc line current to be a mirror image of the ac line voltage is discussed. The power circuit operates at 2 kW peak and 1 kW average power. The BC/SCI has two major subsystems: (1) the battery charger power electronics with its controls; and (2) a microcomputer subsystem which is used to acquire battery terminal data and exercise the state of charge software programs. The state of charge definition employed is the energy remaining in the battery when extracted at a 10 kW rate divided by the energy capacity of a fully charged new battery. The battery charger circuit is an isolated boost converter operating at an internal frequency of 20 kHz. The switches selected for the battery charger are the single most important item in determining its efficiency. The combination of voltage and current requirements dictate the use of high power NPN Darlington switching transistors. The power circuit topology is a three switch design which utilizes a power FET on the center tap of the isolation transformer and the power Darlingtons on each of the two ends. An analog control system is employed to accomplish active input current waveshaping as well as the necessary regulation.

Latos, T. S.

Getting Warmer: IceCube Nears Freeze Out

IceCube has recently detected a diffuse population of high-energy neutrinos arising from the Milky Way. We use this high-significance detection to place new limits on dark matter (DM) annihilation to neutrinos with two complementary approaches. The first method uses the background-subtracted Galactic longitude distribution of shower events to place a conservative bound on the DM annihilation cross section that does not rely on any assumed Galactic cosmic ray emission model; the resulting limits on the velocity-averaged annihilation cross section improve upon existing bounds by factors of a few. The second method uses the template-dependent neutrino energy spectra from the Inner Galaxy, inferred under different Galactic cosmic ray emission models. This complementary approach shows that the inferred Galactic neutrino intensities are already sensitive to DM contributions near the thermal-relic benchmark for a range of TeV-scale DM masses, though this comparison is more model-dependent. Our results demonstrate that measurements of diffuse Galactic neutrino emission can be used as a powerful probe of DM annihilation into neutrinos. Future observations with IceCube-Gen2 and KM3NeT will substantially extend this sensitivity, potentially allowing a decisive test of the thermal freeze-out mechanism with Galactic neutrino observations.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP;

Characterization of Ternary NiTiPt High-Temperature Shape Memory Alloys

Pt additions substituted for Ni in NiTi alloys are known to increase the transformation temperature of the alloy but only at fairly high Pt levels. However, until now only ternary compositions with a very specific stoichiometry, Ni 50-x Pt x Ti 50 , have been investigated and then only to very limited extent. In order to learn more about this potential high-temperature shape memory alloy system, a series of over twenty alloys along and on either side of a line of constant stoichiometry between NiTi and TiPt were arc melted, homogenized, and characterized in terms of their microstructure, transformation temperatures, and hardness. The resulting microstructures were examined by scanning electron microscopy and the phase compositions quantified by energy dispersive spectroscopy. "Stoichiometric" compositions along a line of constant stoichiometry between NiTi to TiPt were essentially single phase but any deviations from a stoichiometry of (Ni,Pt) 50 Ti 50 resulted in the presence of at least two different intermetallic phases, depending on the overall composition of the alloy. Essentially all alloys, whether single or two-phase, still under went a martensitic transformation. It was found that the transformation temperatures were depressed with initial Pt additions but at levels greater than 10 at.% the transformation temperature increased linearly with Pt content. Also, the transformation temperatures were relatively insensitive to alloy stoichiometry within the range of alloys examined. Finally, the dependence of hardness on Pt content for a series of Ni 50-x Pt x Ti 50 alloys showed solution softening at low Pt levels, while hardening was observed in ternary alloys containing more than about 10 at.% Pt. On either side of these "stoichiometric" compositions, hardness was also found to increase significantly.

Transformation Temperatures

Scattering theory of frequency-entangled biphoton states facilitated by cavity polaritons

The use of quantum light to probe exciton properties in semiconductor and molecular nanostructures typically occurs in the low-intensity regime. A substantial enhancement of exciton-photon coupling can be achieved with photonic cavities, where excitons hybridize with cavity modes to form polariton states. Here, to provide a theoretical framework for interpreting emerging experimental efforts in this direction, we develop a scattering theory describing the interaction of frequency-entangled photon pairs with cavity polariton and bipolariton states under various coupling regimes. Employing the Tavis-Cummings model in combination with our scattering approach, we present a quantitative analysis of how the interaction of the entangled photon pair with the polariton or bipolariton modifies its joint spectral amplitude (JSA). Specifically, we examine the effects of the cavity-mode steady-state population, exciton-cavity coupling strength, and different forms of the input photon JSA. Our results show that the entanglement entropy of the scattered photons is highly sensitive to the interplay between the input JSA and the spectral line shapes of the polariton resonances, emphasizing the cavity filtering effects. We suggest that biphoton-scattering quantum light spectroscopy best serves as a sensitive probe of polariton and bipolariton states in the photon-vacuum cavity state. Our approach is not only robust to various regimes of cavity-exciton coupling, but also amenable to extensions beyond the Tavis-Cummings model, enabling the representation of a broad class of molecular systems and solid state quantum materials.

36 MATERIALS SCIENCE

Toward an AI-Powered Software Pipeline for Real-Time Tracking and Analysis of Wildfire and Smoke

Real-time tracking of wildfires and smoke is crucial for effective response, minimizing damage, protecting lives, and efficiently managing resources during fire emergencies. We develop a web-based AI-powered pipeline that detects wildfires in aerial video and estimates deployment-relevant behavior metrics, including cumulative burned area, burned-area growth rate, fire spread direction, and smoke dispersion. The system combines a YOLO-based detector with YCbCr-based fire segmentation, HSV-based smoke segmentation, Farneback optical flow, and centroid-based spatiotemporal tracking. Using ground sampling distance (GSD), pixel-level fire masks are converted to physical burned-area measurements by correlating fire pixel counts with camera altitude and tilt angle. We benchmark YOLO variants and non-YOLO baselines (GoogLeNet, CNN, DBN, Autoencoder, U-Net, and AlexNet) on the IEEE FLAME dataset and a newly created aerial frame dataset, Wildfire-DB. Cross-dataset evaluation uses a strict threshold-transfer protocol: decision thresholds are selected on FLAME validation and transferred unchanged to Wildfire-DB to quantify generalization under domain shift. YOLOv6 achieves the strongest cross-dataset frame-level fire detection on Wildfire-DB (ROC-AUC 0.8200, PR-AUC 0.8044, and transferred-threshold F1 0.7596). For tracking-oriented deployment requiring oriented localization, YOLO11-OBB provides the most reliable cross-dataset behavior among OBB-capable models while remaining computationally feasible. To analyze the feasibility of UAV deployment, we further measure inference efficiency using synchronized GPU and CPU power logs on a fixed workload of 1569 frames. YOLO-family models process the video in 5.73–12.47 seconds with net energy of 1247.28–1775.39 J, substantially lower latency and energy than heavier classification and reconstruction baselines. Overall, model optimality depends on operational objectives: YOLOv6 is best for cross-dataset detection robustness, whereas YOL...

Color segmentation

Monolayer Films of Binary Metal Oxide Nanoparticles for Carbon Nanotube Synthesis

Binary metal oxide nanoparticles (biMO-NPs) combining transition and main-group metals with well-organized atomic architectures have unique potential as robust and economical heterogeneous catalysts. Herein, metal oxide nanoparticles (CoAlxOy) using cobalt and aluminum were synthesized, and their structure, morphology, and chemical composition were investigated using various characterization techniques. The biMO-NPs, which had a Gaussian-type size distribution (∼6 nm), were assembled on plain silicon substrates modified chemically with linker molecules bearing suitable end groups. Surface analysis revealed an ordered, uniform, close-packed arrangement of biMO-NPs forming a monolayer architecture with nanoscale thickness (∼12 nm), high surface uniformity, and negligible defects. The nanoparticle monolayer film (NP-MF) was employed as a catalyst to grow carbon nanotube (CNT) forests via the thermal chemical vapor deposition (CVD) method, similar to those grown from catalyst thin films deposited by physical vapor deposition. Structural characterization confirmed the growth of a dense, uniform, high-quality vertically aligned carbon nanotube (VA-CNT) forest with a length of ∼125 µm, which is comparable to the VA-CNTs grown from expensive thin films. Thus, CNT growth was successfully catalyzed by the biMO-NPs, highlighting a simple and inexpensive alternative for catalyst deposition. In this innovative wet-chemistry approach, the catalyst and catalyst support are combined within a single nanoparticle before NP-MF assembly. Interestingly, this approach provides a scalable and cost-effective pathway for CNT growth, eliminating the need for expensive thin deposition techniques.

Regmi, Bishow [University of Cincinnati]

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE

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