Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle simulations”

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.

At least 181 records · Page 10

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗

Revealing the Formation of the <20 MeV Inner Proton Radiation Belt at L ∼ 2 During the 10–11 May 2024 Superstorm

A new proton radiation belt was identified during the geomagnetic superstorm of 10–11 May 2024. To investigate its origin, we use an MHD-test particle simulation to model solar energetic proton (SEP) trapping and the evolution of the initial trapped proton population during the storm. The simulation weights the injected SEP population with interplanetary proton measurements and the initial trapped proton radial profile with energy-resolved Weather System Follow-on-Microwave (WSF-M) proton data, enabling quantitative comparison with observations. Results show that SEP contributions are minimal; the new ∼2–18 MeV belt at L ∼ 2 forms mainly through redistribution and energization of the initial trapped population within hours of the CME shock, driven by subsequent electric field impulses, rather than the initial shock alone. These results provide quantitative evidence based on realistic, plasma-dependent field dynamics and significantly advance understanding of inner radiation belt formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Generation of megatesla magnetic fields by intense-laser-driven microtube implosions

A microtube implosion driven by ultraintense laser pulses is used to produce ultrahigh magnetic fields. Due to the laser-produced hot electrons with energies of mega-electron volts, cold ions in the inner wall surface implode towards the central axis. By pre-seeding uniform magnetic fields on the kilotesla order, the Lorenz force induces the Larmor gyromotion of the imploding ions and electrons. Due to the resultant collective motion of relativistic charged particles around the central axis, strong spin current densities of ~ peta-ampere/cm 2 are produced with a few tens of nm size, generating megatesla-order magnetic fields. The underlying physics and important scaling are revealed by particle simulations and a simple analytical model. The concept holds promise to open new frontiers in many branches of fundamental physics and applications in terms of ultrahigh magnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental observation of a field-aligned ion beam produced by magnetic reconnection of two flux ropes

An ion beam field-aligned to the background guide field ($B$ 0 = $330$ G) was observed in a reconnection experiment on the Large Plasma Device and, to the authors' knowledge, is the first experimental observation of its kind. Two kink-unstable flux ropes ( L = 11 m, d = 7.6 cm) were made to collide, which allows magnetic reconnection to occur. Sub-Alfvénic ion beams with energies of up to 15 eV were then observed from measurements of the local ion energy distribution function. The beam ions do not appear to be heated. Further, they were correlated with the collision of the ropes and appear to be energized by magnetic reconnection. The results and interpretation of the measurements are supported by three-dimensional gyrokinetic particle simulations of the merging flux ropes and electric field measurements from previous experiments [W. Gekelman et al. , Astrophys. J. 853, 33 (2018)]. The mechanism behind the acceleration appears to be non-local.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding the electromagnetic topology during the ohmic breakdown in tokamaks considering self-generated electric fields

We report that the physical mechanisms of the ohmic breakdown in a tokamak have been understood based on the classical Townsend avalanche theory. However, a new systematic theory (Yoo et al 2018 Nat. Commun.9 3523) recently demonstrated that electron avalanches during the ohmic breakdown are completely different from the Townsend avalanche due to strong self-generated electric fields. In this study, we elucidate the multi-dimensional effects of the self-generated electric field on plasma dynamics during the ohmic breakdown. We also propose a novel electromagnetic topology analysis method that can easily predict the overall plasma behavior and where the main plasma is generated. The topology analysis method is validated by a state-of-art particle simulation for various magnetic configurations. New physical insights into the complex electromagnetic topology would facilitate designing more reliable and optimized ohmic breakdown scenarios in future tokamaks, such as ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗

Granular Flow

LIGGGHTS-INL is a capability-extended version of LIGGGHTS. LIGGGHTS is an Open Source Discrete Element Method Particle Simulation Software (https://www.cfdem.com/media/DEM/docu/Manual.html). LIGGGHTS stands for LAMMPS improved for general granular and granular heat transfer simulations. LAMMPS is a classical molecular dynamics simulator. It is widely used in the field of Molecular Dynamics (https://lammps.sandia.gov/doc/Manual.html). Thanks to physical and algorithmic analogies, LAMMPS is a very good platform for DEM simulations. LAMMPS offers a GRANULAR package to perform these kinds of simulations. LIGGGHTS aims to improve those capability with the goal to apply it to industrial applications.

Xia, Yidong↗

High Intensity Gamma Ray Source Scintillation Attenuation Spectrometer and Filter Stack Monoenergetic Calibration

The Laser-Based X-ray Radiographic Imaging team at Los Alamos National Lab is looking for a monoenergetic MeV X-ray source to both verify our Monte Carlo N-particle simulations of detector performance and calibrate the instruments for future measurements. Funded by the Laboratory Directed Research and Development (LDRD) program, our overall goal is to improve the radiographic quality and reliability of laser-based X-ray sources for deployment at both dynamic and static radiography facilities. Laser-based X-ray sources have demonstrated smaller spot sizes to current electron accelerator based sources. The smaller laser spot size leads to significantly improved resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Utilizing HYSPLIT for Emergency Response Modeling at SRS

The threat of a disaster happening at the Savannah River Site (SRS) is low, but not zero. When an emergency arises, emergency responders utilize a variety of tools to detect, track, and mitigate the disaster. HYSPLIT is a community atmospheric dispersion model that is highly configurable and has been used operationally to simulate particle, chemical, and radionuclide releases. We have developed a preliminary user interface (UI) that configures HYSPLIT to run for a variety of Emergency Action Level (EAL) scenarios at SRS. Using a hypothetical EAL scenario of an accidental release of Pu-238 in H-area and K-area onsite, an analysis was done to show how a preliminary UI framework allows HYSPLIT to be run on a terminal without the need form an active network connection. The study was successful in multiple aspects. For the preliminary UI, users were able to input certain characteristics of the hypothetical EAL such as hours of emission (how long was Pu-238 being released), curies of release, release start time (when does the model run start), and type of radionuclide directly into the terminal. For HYSPLIT, the model was successful in being able to take the UI input, calculate the results from the inputs and output multiple figures of both total volumetric airborne concentrations and surface deposition. We used the Pu-238 Derived Intervention Level (DIL) deposition thresholds for produce, dairy, and beef for visualization of the surface deposition, showing how a release would impact ingestion pathways.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Central Density and Low-Mode Perturbation Control of Inertial Confinement Fusion Dynamic-Shell Targets

The dynamic-shell target is a new class of design for inertial confinement fusion (ICF). These targets address some of the target fabrication challenges prevalent in current ICF targets and take advantage of advances in manufacturing technologies. This study first examines how the dynamic-shell design can be used to control the density of the central region and therefore convergence ratio, thus expanding the design space for ICF. Additionally, the concern of low-mode perturbation growth is considered. A new class of high-performing beam configurations, based on icosahedral polyhedra and charged-particle simulations is proposed. These configurations achieve low levels of irradiation nonuniformity through selection of beam shapes that suppress the dominant symmetrical mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probability of Forming Gaps in the GD-1 Stream by Close Encounters of Globular Clusters

One of the most intriguing properties of the GD-1 stellar stream is the existence of three gaps. If these gaps were formed by close encounters with dark matter subhalos, the GD-1 stream opens an exciting window through which we can see the size, mass, and velocity distributions of the dark matter subhalos in the Milky Way. However, in order to use the GD-1 stream as a probe of the dark matter substructure, we need to disprove that these gaps are not due to the perturbations from baryonic components of the Milky Way. Here we ran a large number of test-particle simulations to investigate the probability that each of the known globular clusters (GCs) can form a GD-1-like gap, by using the kinematical data of the GD-1 stream and GCs from Gaia early data release 3 and by fully taking the observational uncertainty into account. We found that the probability that all of the three gaps were formed by GCs is as low as 1.7 × 10 -5 , and the expected number of gaps formed by GCs is only 0.057 in our fiducial model. Our result highly disfavors a scenario in which GCs form the gaps. Given that other baryonic perturbers (e.g., giant molecular clouds) are even less likely to form a gap in the retrograde-moving GD-1 stream, we conclude that at least one of the gaps in the GD-1 stream was formed by dark matter subhalos if the gaps were formed by flyby perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase-space Properties and Chemistry of the Sagittarius Stellar Stream Down to the Extremely Metal-poor ([Fe/H] ≲ -3) Regime

In this work, we study the phase-space and chemical properties of the Sagittarius (Sgr) stream, the tidal tails produced by the ongoing destruction of the Sgr dwarf spheroidal (dSph) galaxy, focusing on its very metal-poor (VMP; [Fe/H] < -2) content. We combine spectroscopic and astrometric information from SEGUE and Gaia EDR3, respectively, with data products from a new large-scale run of the StarHorse spectrophotometric code. Our selection criteria yield ~1600 stream members, including >200 VMP stars. We find the leading arm (b > 0°) of the Sgr stream to be more metal-poor, by ~0.2 dex, than the trailing one (b < 0°). With a subsample of turnoff and subgiant stars, we estimate this substructure's stellar population to be ~1 Gyr older than the thick disk's. With the aid of an N-body model of the Sgr system, we verify that simulated particles stripped earlier (>2 Gyr ago) have present-day phase-space properties similar to lower metallicity stream stars. Conversely, those stripped more recently (<2 Gyr) are preferentially akin to metal-rich ([Fe/H] > -1) members of the stream. Such correlation between kinematics and chemistry can be explained by the existence of a dynamically hotter, less centrally concentrated, and more metal-poor population in Sgr dSph prior to its disruption, implying that this galaxy was able to develop a metallicity gradient before its accretion. Finally, we identified several carbon-enhanced metal-poor ([C/Fe] > +0.7 and [Fe/H] ≤ -1.5) stars in the Sgr stream, which might be in tension with current observations of its remaining core where such objects are not found.

79 ASTRONOMY AND ASTROPHYSICS↗

Myths of German Graphite in World War II, with Original Translations

We re-examine the common narrative that a 1941 experimental error by physicists Walther Bothe and Peter Jensen led Germany to abandon graphite as a reactor moderator during World War II. We first detail the history of both German and American graphite experiments, noting that the Americans faced similar setbacks but succeeded only at the end of a costly 18-month graphite purification program. We then use Monte Carlo N-Particle simulations to reconstruct Bothe's 1941 experiment. We find the thermal absorption cross section of Bothe's Siemens electrographite to be 12.2 mb, in contrast to his reported 7.9 mb. This discrepancy arises because the neutrons in Bothe's experiment did not reach thermal equilibrium, leading to an underestimation of neutron absorption. Additionally, despite misconceptions that the Germans were unaware of boron impurities, we share evidence that Wilhelm Hanle accurately measured boron and cadmium impurities in the electrographite. To support our findings, we provide 9 excerpted and 3 complete English translations of classified wartime reports by Heisenberg, Joos, Bothe, Jensen, Höcker, Hanle, and Kremer. Our work aims to illuminate the rational constraints behind Germany's decision to forgo graphite moderation.

Park, Patrick J.↗

High Power Superconducting Electron Linear Accelerators For Industrial Application

Fermilab is executing a technology development program to develop a compact yet powerful electron accelerator. We are leveraging R&D breakthroughs in SRF cavities, cost-effective radio-frequency sources, modern cryo-coolers, and high average current electron guns. We show that a single accelerator module can deliver average beam power as high as 1 MW through detailed thermal, RF, and particle simulations.

Thangaraj, Jayakar Tobin↗

Measurement of the muon neutrino charged-current mesonless cross section in the NOvA near detector

NOvA is a long-baseline accelerator neutrino experiment at Fermilab. Its physics goals include precision neutrino oscillation measurements, neutrino interaction cross-section measurements and beyond Standard Model explorations. We present a measurement of muon neutrino charged-current cross section with zero mesons in the final state at the NOvA near detector. This measurement is performed as a function of the kinematics of the final state muon. Our chosen interaction channel is especially sensitive to quasielastic and meson exchange current interactions and it provides handles for constraining the cross section systematic uncertainties in oscillation analyses in present and future experiments. For particle identification, we use a convolutional neural network (CNN) trained on individual particles simulated in the NOvA Near detector. This allows us to select the desired signal while reducing the potential bias from neutrino interaction modeling. We study strategies for constraining the remaining charged-pion background via Michel electron information in a template fitting approach. The main experimental result is a two-dimensional differential cross section as a function of final-state muon kinetic energy and polar angle. The parameters of this measurement, including binning and unfolding, were optimized to reduce the expected systematic uncertainty in the total cross section. The final result shows good agreement with the main GENIE-based simulation framework that was independently fine-tuned in NOvA. We finally propose improvements and subsequent steps that build on this analysis and further dissect the final states of neutrino interactions. This work has been supported by US DOE grant DE-SC0015684.

Sánchez Falero, Sebastián Jesús [Iowa State U.]↗

The Ice Particle and Aggregate Simulator (IPAS). Part III: Verification and Analysis of Ice–Aggregate and Aggregate–Aggregate Collection for Microphysical Parameterization

Abstract The Ice Particle and Aggregate Simulator (IPAS) is used to theoretically represent the aggregation process of ice crystals. Aggregates have a variety of formations based on initial ice particle size, shape, and falling orientation, all of which influence water phase partitioning. Aggregate dimensional properties and density changes are calculated for monomer–monomer (MON–MON), monomer–aggregate (MON–AGG), and aggregate–aggregate (AGG–AGG) collection to be used by ice-microphysical models for improvement in aggregation parameterizations. Aggregates are chosen from a database of 9 744 000 preformed combinations to be further collected (see Part II). AGG–AGG collection results in more extreme and a smaller range of aggregate aspect ratios than MON–AGG collection. A majority of aggregates are closer to prolate than oblate spheroids, regardless of collection type, except for quasi-horizontally oriented particles that have extreme aspect ratios to begin with. MON–AGG collection frequently results in an increase in density upon collection, whereas MON–MON and AGG–AGG collection almost always result in particle density decreases, with extreme reductions near 99% for MON–MON collection. MON–MON collection results in the greatest decreases in density but then quickly becomes unaffected by the addition of more monomers due to inherent size differences between monomers and aggregates. Finally, a holistic analysis to in situ observations of cloud particle images is presented. IPAS 2D aspect ratios surround a median value of 0.6 and closely follow that of previous studies while varying by no more than ≈12% on average from observed aggregates.

54 ENVIRONMENTAL SCIENCES↗

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

We present two-dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein–Greene–Kruskal modes, which are exact nonlinear steady-state solutions of the Vlasov–Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 20482 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Physics↗