Search NASA⌕ Search

SEARCH · Search NASA

Results for “transportation 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 271 records · Page 15

Achievement of a high-density, high-confinement, and high-beta tokamak plasma regime in DIII-D, and implications for a lower-current path for ITER and FPP

Experiments on DIII-D have demonstrated a density-confinement synergy that enables sustainment of high performance in a previously unattained parameter regime of simultaneous very high energy confinement quality (H 98y2 ≥ 1.5), very high line-average density Greenwald fraction (ƒ Gr = πa 2 < n >/I P ≥ 1.4), and high toroidal beta (β T ≥ 3%). Tokamak operation in this regime is essential for a compact steady-state FPP, as well as for Q=10 with 500 MW of fusion power in ITER at I P << 15 MA. These experiments leveraged the knowledge that, in the high-poloidal-beta (β P ) regime, impurity and density gradients can enhance turbulence stabilization caused by high α MHD (α MHD ~(dβ P )⁄dr). This was described by theoretical predictions and gyrokinetic transport simulations [M.T. Kotschenreuther et al, 2024 Nucl. Fusion, 64 076033], and later confirmed by experiments on DIII-D [S. Ding et al, 2024 Nature 629 555]. To increase both β P and β T , the new experiments increased the ideal-wall stability β N -limit by using a smaller plasma-outer wall distance and higher triangularity in the plasma cross section (top/bottom average δ~0.9), enabled by the recent “shape & volume rise” (SVR) modification to the DIII-D divertor. The higher triangularity also contributed to achieving higher ƒ Gr by enabling higher pedestal density. At high density, the pedestal is ballooning limited and exhibits small and frequent ELMs, while the divertor is near detachment even without any impurity seeding. High plasma performance was attained and sustained reproducibly, with the eventual terminations brought about by an MHD mode destabilized as the current profile slowly continued to evolve. A path to stationary fully noninductive operation might include ECH injection to reduce both core impurity accumulation and the electron collisionality, thus increasing the bootstrap current. These experiments provide the first experimental demonstration of the ƒ Gr , H 98y2 , and β T values required simultaneously for ITER Q = 10 at I P < 10 MA, pointing to practical ways to improve the energy confinement in a fusion reactor.

Garofalo, Andrea M. [General Atomics, San Diego, C↗

Superconducting proximity effect in a transparent van der Waals superconductor-metal junction

Here, we report on Andreev reflections at clean NbSe 2 -bilayer graphene junctions. The high transparency of the junction, which manifests as a large conductance enhancement of up to 1.8, enables us to see clear evidence of a proximity-induced superconducting gap in bilayer graphene and two Andreev reflections through a vertical NbSe 2 -graphene and a lateral graphene-graphene junction, respectively. Quantum transport simulations capture the complexity of the experimental data and illuminate the impact of various microscopic parameters on the transmission of the junction. Our work establishes the practice and understanding of an all-van-der-Waals, high-performance superconducting junction. The realization of a highly transparent proximized graphene-graphene junction opens up possibilities to engineer emergent quantum phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Giant spin transport anisotropy in magnetic topological insulators

Here, we report on exceptionally long spin transport and giant spin lifetime anisotropy in the gapped surface states of three-dimensional (3D) magnetic topological insulators (MTIs). We examine the properties of these states using the three-dimensional Fu-Kane-Mele Hamiltonian model in the presence of a magnetic exchange field. The corresponding spin textures of surface states, which are well reproduced by an effective two-band model, hint at a considerable enhancement of the lifetime of out-of-plane spins compared to in-plane spins. This is confirmed by large-scale spin transport simulations for 3D MTIs with disorder. The energy dependence of the spin lifetime anisotropy arises directly from the nontrivial spin texture of the surface states, and is correlated with the onset of the quantum anomalous Hall phase. Our findings suggest different spin filtering capabilities of the gapped surface MTI states, which could be explored by nonlocal spin valve measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spectroscopic Measurements and Models of Energy Deposition in the Substrate of Quantum Circuits by Natural Ionizing Radiation

Naturally occurring background radiation is a potential source of correlated decoherence events in superconducting qubits that will challenge error-correction schemes. In order to characterize the radiation environment in an unshielded laboratory representative of superconducting qubits’ environments, we performed broadband, spectroscopic measurements of background radiation events inside a millikelvin refrigerator. The spectrometer was designed to mimic the size and composition of a quantum circuit. Specifically, we measured the background radiation spectra in silicon substrates of two thicknesses, 500 and 1500 µm, and one area, 25 mm 2 . The observed spectra span energies from a few kilo-electron-volts up to nearly 10 MeV, are nearly featureless, and decrease in intensity by a factor of 40 000 between 100 keV and 3 MeV for the 500-µm substrate. We integrate the spectra to obtain the average event rates and deposited power levels. These quantities correspond to a rate of 0.023 events per second and a power of 4.9 keV s -1 , when counting events that deposit at least 40 keV for the 500-µm-thick substrate. We find that the cryogenic measurements are in good agreement with predictions based on simple measurements of the terrestrial gamma-ray flux outside the refrigerator, published models of cosmic-ray fluxes, a crude model of the cryostat, and radiation-transport simulations. This model requires no free parameters to predict the background radiation spectra in the silicon substrates. The agreement between measurements and predictions demonstrates that the model we present can be used to assess the relative contributions of terrestrial and cosmic-ray sources to background radiation interactions in silicon substrates of varying thickness. These spectroscopic measurements are performed with a novel combination of superconducting microresonators located on micromachined silicon islands that define the interaction volume with background radiation. The resonators transduce deposited energy to a readily detectable electrical signal. Microresonator readout closely resembles dispersive superconducting qubit readout, so similar devices—with or without micromachined islands—are suitable for integration with superconducting quantum circuits as detectors for background radiation events. For our specific laboratory conditions, we find that gamma-ray emissions from radioisotopes are responsible for the majority of events that deposit E < 1 ⁢Me⁢V. We present results demonstrating that the background radiation spectrum contains relevant contributions from cosmic-ray particles other than muons, particularly a tail of multi-mega-electron-volt events due to protons and neutrons. These observations suggest several paths to reducing the impact of background radiation on quantum circuits, supported by an empirically validated model for generating reliable predictions of radiation interactions with silicon substrates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pricing Strategy of Electric Vehicle Aggregators Based on Locational Marginal Price to Minimize Photovoltaic (PV) Curtailment

The global climate crisis demands urgent action to mitigate global warming. Using renewable energy sources, such as solar and wind power, for electricity generation is crucial. This shift from centralized to distributed power systems, however, brings challenges, including voltage fluctuations and renewable energy curtailment. The rapid growth of the electric vehicle (EV) industry adds complexity, increasing overall electricity demand and straining the power supply during peak charging times. This paper proposes a scheduling strategy for EV aggregators to reduce renewable energy curtailment and stabilize grid operation by strategically scheduling EV charging. Using Multi -Agent Transport Simulation (MATSim), a traffic simulation tool, EV driving data in Denver, Colorado, USA, were modeled. The EV aggregator adjusts charging fees based on locational marginal prices, encouraging EVs to charge at different stations according to pricing. Simulations on an IEEE 33-bus system with distributed energy resources and EV charging stations validate the proposed algorithm, demonstrating its effectiveness in reducing curtailment by 12.55% and stabilizing grid operation.

33 ADVANCED PROPULSION SYSTEMS↗

Testing the Activation Analysis for Fusion in OpenMC

OpenMC is a community-developed Monte Carlo neutron and photon transport simulation code. It can perform fission simulations such as fixed-source, k-eigenvalue, and subcritical multiplication calculations on models built using either a constructive solid geometry or CAD representation. To explore the use of OpenMC for fusion activation analysis, a detailed model of the Fusion Neutronics Science Facility (FNSF) was first developed for comparisons against an existing SERPENT model. A 90-degree model of FNSF in Standard-Triangle-Language (STL) CAD format was converted to Constructive Solid Geometry (CSG) using each code's built-in functions, and the geometries were validated by ensuring no cells overlapped and no particles were lost during simulations. The neutron fluxes were calculated and compared for multiple components close to the plasma. The results show differences mostly below 1% in fluxes and averaged 8% for activity and decay heat. Here, the work described in this study tests the CAD-based geometry using the DagMC toolkit in OpenMC and compares the activation analysis of OpenMC to SERPENT code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Aggregated EV Flexibility With TSO-DSO Coordination

Electric vehicles (EVs) are becoming a promising source of grid ancillary services due to the temporal and spatial charging flexibility, quick response and storage capability. Such advantages are increasing with government policy promotion and technology improvement. However, the exploration of EV flexibility requires the coordination of both transmission system operators (TSOs) and distribution system operators (DSOs), to ensure the safe and reliable operation of power network. In this paper, we propose a coordinated evaluation method that determines the optimal utilization of EV temporal flexibility without compromising EV owners’ usage. At the distribution level, DSOs first evaluate EV aggregators’ operational boundaries to exploit distribution level services. At the transmission level, TSOs then determine EV charging schedules and ancillary service capacity simultaneously, taking into account the requirement from DSOs. Here, we validate the model in a case study using the IEEE 123 node test feeder and EV charging sessions obtained from a transportation simulation tool that uses real-world data.

33 ADVANCED PROPULSION SYSTEMS↗

The high-density equation of state in heavy-ion collisions: constraints from proton flow

Abstract A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $$m_T$$ m T . We also show that an increase in $$v_2$$ v 2 is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $$n_{\text {CE}}<4 n_0$$ n CE < 4 n 0 , show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $$4n_0$$ 4 n 0 .

74 ATOMIC AND MOLECULAR PHYSICS↗

Adapting CLUTCH methodology to multigroup TSUNAMI-3D for eigenvalue sensitivity calculations

The sensitivity of the eigenvalue to uncertainties in nuclear data and its evaluation are important for nuclear criticality safety. TSUNAMI-3D sequences within the SCALE code system offer several options to the user community for calculating eigenvalue sensitivity coefficients with multigroup (MG) and continuous energy (CE) 3D transport capabilities. TSUNAMI-3D sequences implement the adjoint-based perturbation theory with MG KENO code, the Contributon Linked eigenvalue sensitivity/Uncertainty estimation via Track length importance CHaracterization (CLUTCH) method with CE KENO code, and the Iterated Fission Probability (IFP) method with CE KENO and Shift codes. Each method has benefits and limitations depending on the problem that is run. The work presented here aims to adapt the CLUTCH method, which enables the Contributon method's mesh-free, memory-efficient approach for calculating adjoint-weighted tallies for sensitivity calculations, to the MG TSUNAMI-3D sequence. This application would eliminate the explicit adjoint KENO calculation, as well as the memory-consuming mesh flux moment tallies required by the conventional MG TSUNAMI-3D. Smaller memory footprints in the CLUTCH methodology and relatively shorter runtimes in MG KENO transport can make MG TSUNAMI-3D a viable method for some complex problems. Moreover, this adaptation allows MG sensitivity calculations with Shift, ORNL's next-generation high-performance Monte Carlo transport code, which currently does not offer any sensitivity capabilities with MG particle transport simulations. Initial implementation of the new MG TSUNAMI-3D sequence and its preliminary results with a selected critical benchmark experiment in the Verified, Archived Library of Inputs and Data (VALID) are presented in this study.

KENO↗

Polaris Studio

SF-24-040 Polaris-studio is the Python entry point/front end for using the Polaris transportation simulator. The package contains the code used to create new models from open and user-provided data, an extensive range of convenience tools for data preparation and result analysis and the ability of designing, launching and tracking studies combining multiple individual Polaris simulations from a Python terminal. The Polaris Studio software contains several modules for accomplishing the above tasks, including the PolarisLib - the Polaris object model, Polaris Manager - code for setting up and running Polaris studies, QPOLARIS - an add-on for the QGis software to allow model files to be created and edited, PolarisVis - a juptyer-based postprocessing library, and Polaris Tools - a useful collection of supplemental tools for interacting with Polaris files.

Auld, Joshua↗

pymcnp v1

This is a software package that makes it easier to handle MCNP, a neutron transport simulation toolkit, simulations. It helps with reading and writing output from MCNP and modifying the input files programmatically. Furthermore, it also allows the visualization of input geometry and of the results. Pymcnp provides access to these features using the Python programming language.

Persaud, Arun↗

Local models for scatter estimation and descattering in polyenergetic X-ray tomography

We propose a new modeling approach for scatter estimation and descattering in polyenergetic X-ray computed tomography (CT) based on fitting models to local neighborhoods of a training set. X-ray CT is widely used in medical and industrial applications. X-ray scatter, if not accounted for during reconstruction, creates a loss of contrast in CT reconstructions and introduces severe artifacts including cupping, shading, and streaks. Even when these qualitative artifacts are not apparent, scatter can pose a major obstacle in obtaining quantitatively accurate reconstructions. Our approach to estimating scatter is, first, to generate a training set of 2D radiographs with and without scatter using particle transport simulation software. To estimate scatter for a new radiograph, we adaptively fit a scatter model to a small subset of the training data containing the radiographs most similar to it. We compared local and global (fit on full data sets) versions of several X-ray scatter models, including two from the recent literature, as well as a recent deep learning-based scatter model, in the context of descattering and quantitative density reconstruction of simulated, spherically symmetrical, single-material objects comprising shells of various densities. Our results show that, when applied locally, even simple models provide state-of-the-art descattering, reducing the error in density reconstruction due to scatter by more than half.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling the impact of riparian hollows on river corridor nitrogen exports, Frontiers in Water: Dataset

This package contains all reactive transport simulation outputs used in "Modeling the impact of riparian hollows on river corridor nitrogen exports" (Rogers, D.B., et al., 2021). The simulation outputs reflect hydrological and biogeochemical behavior of a single, generalized riparian hollow located in the intra-meander region of a sub-alpine floodplain under transient hydrological conditions. The context for, purpose of, and conclusions drawn from this model are detailed in Rogers, D.B., et al., 2021.The .zip file contains simulation outputs for each of the five time periods analyzed (see metadata for specific dates), with each folder representing one continuous simulation period. Within each folder are the raw simulation outputs, which are comprised of hundreds of files. Files ending in ".vtk" are best visualized and analyzed using ParaView (https://www.paraview.org). All other files are best viewed or analyzed using R. For a detailed description of the meaning of each file, refer to the ".fls" file located in each folder.

54 ENVIRONMENTAL SCIENCES↗

Design of a 10 MeV beamline at the Upgraded Injector Test Facility for e-beam irradiation

Electron beam irradiation near 10 MeV is suitable for wastewater treatment. The Upgraded Injector Test Facility (UITF) at Jefferson Lab is a CW superconducting linear accelerator capable of providing an electron beam of energy up to 10 MeV and up to 100 µA current. This contribution presents the beam transport simulations for a beamline to be used for the irradiation of wastewater samples at the UITF. The simulations were done using the code General Particle Tracer with the goal of obtaining an 8 MeV electron beam of radius (3-σ) of ~2.4 cm. The achieved energy spread is ~74.5 keV. The space charge effects were investigated when the bunch charge is varied to be up to 1000 times and the results showed that they do not affect the beam quality significantly.

Li, X.↗

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗

Effects of delta ray electrons on measurement uncertainties of harp system

A harp system, which is a multi-wire beam profile monitoring (MWPM) system, is planned upstream of the spallation target to make in situ calibration of beam current density configuration on the target along with beam imaging from luminescent coating on the beam entrance window at the Second Target Station (STS) of the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL). This beam interception-based beam diagnostics system on the target will be used to ensure that the maximum beam loads on the target are within the design range during neutron production. Current design of the harp consists of three layers of measurement wires each of which is sandwiched between voltage biasing wire planes. The signal obtained from each measurement wire layer is disturbed by secondary electrons (SE) and delta rays produced by beam-matter interactions in neighboring wires and ionization of residual gases in accelerator vacuum. While the backgrounds from SE can be suppressed by voltage biasing, the delta-ray electrons with kinetic energies above keV ranges overcome the electric potential bias. In this paper, we study the effects of delta-rays on the measurement uncertainties of MWPM using the particle transport simulation code FLUKA. Furthermore, the cases where the harp system is installed in the proximity of a large delta ray sources such as proton beam window or in the core vessel filled with sub-atmospheric gas have been studied.

Lee, Yong Joong↗

Physics Considerations for a Harp System Design at the Second Target Station of the Spallation Neutron Source

A harp system is being developed for monitoring proton beam profile direct upstream of the proton beam window at the Second Target Station of the Spallation Neutron Source, Oak Ridge National Laboratory. It consists of two sensor planes which have arrays of thin conducting wires aligned vertically and horizontally, respectively. It monitors beam profiles in two transverse directions to the beam axis by measuring the net-charge depositions in the sensor wires, which are caused by ejection of secondary electrons and delta rays driven by electromagnetic interactions with high energy protons. The net charge deposition in a sensing wire linearly correlates with the number of incident protons on it. This correlation is perturbed when the wire interacts with secondary electrons and delta rays originating from beam matter interactions in neighboring wires, PBW and residual gases. In this paper, we analyze the physical phenomena that affects the measurement uncertainties of the harp using particle transport simulations.

Lee, Yong Joong [ORNL] (ORCID:0000000298381723)↗

An Analytic Benchmark for the Solution to the Isotopic Fission Spectrum Mixture Problem

In neutron transport calculations, it is common to specify material compositions in terms of constituent isotopes. Material compositions may be described by isotope number densities and associated microscopic cross sections. For general reaction cross sections, the macroscopic cross sections of a composition are simply the summation of the sum of the products of isotope number densities and microscopic cross sections. A notable exception is the mixture of the neutron fission spectrum in fissile material. To demonstrate proper and improper mixture of the neutron fission spectrum, a zero-dimensional test problem is developed. Using the test problem, it is demonstrated that the improper mixing of the fission spectrum results in an eigenvalue error of approximately 65 pcm. Though this may seem small, any error in such a simple calculation is unacceptable. In similar infinite-homogeneous test problems, eigenvalue errors of more than 300 pcm have been observed. It is also shown that the error due to fission spectrum mixing can be exaggerated for computer program verification purposes. It is concluded that the proper mixture of the neutron fission spectrum is essential for accurate neutron transport simulations. The effect of improper neutron fission spectrum mixture is demonstrated and quantified and this test problem may be used for verification purposes in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗