Search NASA⌕ Search

SEARCH · Search NASA

Results for “MULTIPLIER”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Risk assessment of wellbore leakage during underground hydrogen storage

The expansion of renewable energy sources would require large-scale energy storage options to overcome the intermittent nature of these sources. Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs offers a scalable and practical energy storage solution. These reservoirs are chosen for their availability and large capacity, but the unique properties of hydrogen raise concerns about potential leakage pathways, particularly through wellbores. In this study, we develop and apply, for the first time, reduced-order models (ROMs) specifically designed for efficient leakage risk prediction in UHS systems operating in depleted hydrocarbon reservoirs. Using 3,000 high-fidelity simulation scenarios, we examine the influence of 11 key parameters, including reservoir and aquifer depths, wellbore permeability and porosity, initial saturations of water, oil and gas fractions (hydrogen, light, intermediate, and heavy hydrocarbons), reservoir pressure multiplier, and the aquifer-to-reservoir volume ratio, to simulate leakage behavior over a 1,000-year timescale. We train ROMs using a two-step classification-regression approach, achieving R 2 values exceeding 99 % across all targets. These ROMs effectively capture the leakage evolution and identify critical controls of leakage, guiding the design of mitigation strategies. Results indicate that gas leakage occurs in about 27 % of scenarios as early as five years post-operation, reaching volumes of up to 106 ft3. Oil leakage is less frequent (~17 %) and typically begins decades later. Our findings also show that hydrogen often migrates first, owing to its smaller molecular size and higher buoyancy, followed by heavier hydrocarbons. Over time, these heavier components contribute significantly to the total leaked volume, reinforcing the need for targeted monitoring and remediation strategies. Our analysis highlights that deeper storage reservoirs, shallower aquifers, and low-permeability wellbores significantly reduce leakage risks. In conclusion, this work offers a robust framework for risk-informed UHS deployment, supporting energy security through reliable large-scale hydrogen storage while safeguarding environmental integrity.

08 HYDROGEN↗
DOE OSTI2026

Link statistics of dislocation network during strain hardening

Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.

Crystal plasticity↗
DOE OSTI2026

Distributed Tomographic Reconstruction with Quantization

Conventional tomographic reconstruction typically depends on centralized servers for both data storage and computation, leading to concerns about memory limitations and data privacy. Distributed reconstruction algorithms mitigate these issues by partitioning data across multiple nodes, reducing server load and enhancing privacy. However, these algorithms often encounter challenges related to memory constraints and communication overhead between nodes. In this paper, we introduce a decentralized Alternating Directions Method of Multipliers (ADMM) with configurable quantization. By distributing local objectives across nodes, our approach is highly scalable and can efficiently reconstruct images while adapting to available resources. To overcome communication bottlenecks, we propose two quantization techniques based on K-means clustering and JPEG compression. Numerical experiments with benchmark images illustrate the tradeoffs between communication efficiency, memory use, and reconstruction accuracy.

Miao, Runxuan↗
DOE OSTI2024

A safe reinforcement learning algorithm for supervisory control of power plants

Traditional control theory-based methods require tailored engineering for each system and constant fine-tuning. In power plant control, one often needs to obtain a precise representation of the system dynamics and carefully design the control scheme accordingly. Model-free Reinforcement learning (RL) has emerged as a promising solution for control tasks due to its ability to learn from trial-and-error interactions with the environment. It eliminates the need for explicitly modeling the environment’s dynamics, which is potentially inaccurate. However, the direct imposition of state constraints in power plant control raises challenges for standard RL methods. To address this, we propose a chance-constrained RL algorithm based on Proximal Policy Optimization for supervisory control. Our method employs Lagrangian relaxation to convert the constrained optimization problem into an unconstrained objective, where trainable Lagrange multipliers enforce the state constraints. In conclusion, our approach achieves the smallest distance of violation and violation rate in a load-follow maneuver for an advanced Nuclear Power Plant design.

constrained optimization↗
DOE OSTI2025

Accurate and rapid acoustic damage characterization in complex structures using sparse sensor networks and deep learning models

Damage diagnosis in critical components is essential for ensuring the safety and reliability of operations across industries, spanning manufacturing, aerospace, and energy. Traditional acoustic nondestructive testing methods primarily focus on detecting defects through the direct scattering of single-mode incident waves from the damage, which limit their applicability to simple structures and small inspection areas. Our earlier research demonstrated that machine learning algorithms combined with sparse sensor networks can identify critical defect signatures even from multiply scattered, multi-mode acoustic signals, indicating the potential for improved defect inspection in complex, real-world structures. In this work, we demonstrate the successful implementation of this approach in a fixed sensor configuration to rapidly and accurately detect simulated defects in a geometrically complex, real-world structure, a brake rotor hub. Three different types of defects were physically simulated on the surface of the hub, and the collected data were used to train an autoencoder-based deep learning model. Two models were tested, one using single measurements and the other using multiple measurements taking advantage of the spatial distribution of the sensor network. After training, the multi-measurement model achieved 100 % accuracy in identifying, classifying, and locating unseen, unique damages. This work illustrates the potential of the proposed method for a wide range of industrial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Mu2e Digitizer ReAdout Controller (DiRAC): characterization and radiation hardness

The Mu2e experiment at Fermilab will search for the neutrino-less coherent conversion of a muon into an electron in the field of a nucleus. Mu2e detectors comprise a straw tracker, an electromagnetic calorimeter and a veto for cosmic rays. The calorimeter employs 1348 Cesium Iodide crystals readout by silicon photo-multipliers and fast front-end, and digitization electronics. The digitization board is named DiRAC (Digitizer ReAdout Controller) and 140 cards are needed for the readout of the full calorimeter. The DiRACs are hosted in crates located on the external surface of calorimeter disks, inside the detector solenoid cryostat and must sustain very high radiation and magnetic field so it was necessary to fully qualify it. Several version of prototypes were validated for operation in a high-vacuum (10−4 Torr) and under a 1T magnetic field. An extensive radiation hardness qualification campaign, carried out with photons, 14 MeV neutron beams, and 200 MeV protons, certified the DiRAC design to sustain doses up to 12 Krad, neutron fluences up to ∼ 1011 1 MeV neq/cm2, and very low occurrences of single-event effects. The qualification campaigns and quality assurance procedures will be reviewed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Fast spark-detection system for GEM detectors

The sPHENIX experiment is currently under commissioning at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL). The Time Projection Chamber (TPC) serves as a tracking detector for the experiment. The sPHENIX TPC uses a stack of four Gas Electron Multipliers (GEMs) as a gain stage in a reduced ion back-flow configuration. To mitigate the damaging effects of sparks in the GEMs, an online spark monitoring system was created. Once the system detects a spark in a GEM stack, the voltages across the GEMs in that stack can be lowered to prevent further sparking without affecting the gain and efficiency of the other modules. Spark signals are coupled out of the GEM stack by a pick-off capacitor attached to the bottom of the bottom GEM. Custom PCBs convert the oscillatory spark signal into a mono-polar pulse that is then digitized. The software then saves the waveform in a server and uses experimentally derived thresholds to determine how to react. As a result, the system has so far proven to be effective at improving the stability of the TPC and preventing damaging events while collecting cosmic ray data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

YSO implantation detector for beta-delayed neutron spectroscopy

A segmented-scintillator-based implantation detector was developed to study the energy distribution of β-delayed neutrons emitted from exotic isotopes. The detector comprises a 34 × 34 YSO scintillator coupled to an 8 × 8 anode Position-Sensitive Photo-Multiplier Tube (PSPMT) via a tapered light guide. The detector was used at RIBF, RIKEN for time-of-flight-based neutron spectroscopy measurement in the 78Ni region. The detector provides the position and timing resolution necessary for ion-beta correlations and time of flight measurements. The detector provides a high ~80% beta-detection efficiency and a sub-nanosecond timing resolution. This contribution discusses the details of the design, operation, implementation, and analysis developed to obtain neutron time-of-flight spectrum and the analysis methods in the context of neutron-rich nuclei in the 78Ni region.

Implantation detector↗
DOE OSTI2025

Energy resolution and gain measurements in Argon-based gas mixtures: Exploring Ar:CF 4 for low energy measurements with TPCs

Time Projection Chambers (TPCs) are among the most advanced charged-particle detectors. Gas-filled TPCs have tracking capabilities that provide 3D-imaging of charged particles with a good energy resolution for spectroscopy. Different gas mixtures have different properties that determine the energy resolution as well as the spatial resolution. Therefore, optimization of operating conditions is required to simultaneously obtain adequate gain, energy resolution, spatial/track resolution, as well as higher drift velocities for high counting rates applications. Ar:CF 4 gas mixture has higher electron drift velocities and lower electron diffusion, which makes it an attractive candidate for TPC filling gas for low energy nuclear physics applications as compared to commonly used Ar:CH 4 and Ar:CO 2 gas mixtures, namely when tracking information is needed. However, other properties, including energy resolution and gain, remain largely unexplored in Ar:CF 4 especially at pressures and other operating conditions relevant for low-energy nuclear physics applications. Here, in this paper we report on gain and energy resolution measurements, using Gas Electron Multipliers (GEMs), in the less explored Ar:CF 4 mixture (Alfonsi et al., 2006), as well as in the more commonly used gas mixtures Ar:CH 4 and Ar:CO 2 . In addition to obtaining energy resolution and gain, we provide results from Garfield++ simulations for gain fluctuations, and their impact on energy resolution is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Investigating the effect of B T direction on W source-to-core pathways during the SAS-VW campaign on DIII-D

Experiments using the V-shaped closed slot tungsten (W) coated SAS-VW divertor in DIII-D studied the effects of the BT direction on core contamination of eroded tungsten from a closed slot divertor configuration. Core W content is inferred using soft-X ray tomography (SXR) and vacuum ultraviolet spectroscopy (SPRED), while W divertor erosion is inferred from visible spectroscopy of W emission (400.9 nm) measured by in-slot filterscopes (filtered photo-multipliers). Post-mortem analysis from the campaign discovered tile misalignment leading to suspected pronounced leading-edge erosion in the unfavorable BT direction (ion B ⃗×∇B ⃗ drift away from divertor) likely not captured by diagnostics. However, empirical findings show up to ~2-3x larger core contamination in the favorable BT direction even considering no additional W erosion from leading edges. A “source-to-core efficiency factor” is derived to estimate the effects of leading-edge erosion and compare W contamination for two pairs of H-mode discharges in opposite BT directions. While having differing absolute parameters, similar core impurity density gradients suggest comparable core impurity transport. These results show that favorable BT may have stronger source-to-core pathways for W impurities sourced from the outer divertor region. Possible explanations could include the effects of E ⃗×B ⃗ drifts on W transport in the scrape-off-layer (SOL) as well as previously determined fast SOL inner target directed flows in favorable BT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impact of nitrogen molecular breakup on divertor conditions in JET L-mode plasmas using SOLPS-ITER

SOLPS-ITER simulations of nitrogen-seeded, low-confinement mode plasmas in the Joint European Torus (JET) predict that the electron temperature in the low-field side (LFS) divertor leg is reduced locally by up to an order of magnitude when nitrogen is assumed to recycle as molecules (N 2 ) instead of atoms using a fixed nitrogen injection rate. The LFS divertor temperature reduction under the assumption of molecular recycling occurs due to a three-step mechanism: (1) the plasma penetration of nitrogen atoms is increased due to the strong triple bond of the N 2 molecule and the kinetic energy release in the dissociation event, both mechanisms contributing equally, (2) the abundance of (particularly multiply-charged) nitrogen ions in the divertor is increased and (3) the electron temperature is reduced due to the increase in radiation (by up to a factor of 4) from nitrogen ions. Setting the volume-integrated nitrogen radiated power to a constant value (0.6 MW) instead of the nitrogen injection rate, SOLPS-ITER predicts under the molecular nitrogen recycling assumption that the peak line-integrated N II, N III and N IV intensities in the LFS divertor are approximately within 15%, 35% and 5%, respectively, of the reference atomic nitrogen recycling case. The predicted peak N II, N III and N IV intensities under either assumption are within 30%, 65% and 5%, respectively, of measurements using the vertically viewing mirror-link divertor spectrometer (Meigs et al., 2010) in nitrogen-seeded JET L-mode plasmas (Lomanowski et al., 2019). ERO2.0 simulations using a constant nitrogen seeding rate on static background plasma solutions from EDGE2D-EIRENE (previously presented in Mäenpää et al., (2022), revised here to include fast reflections) predict that N II to N IV line emission is increased by 20% to 30% when nitrogen is assumed to recycle as molecules, demonstrating the importance of considering the effect of molecular dissociation reactions on the divertor plasma in a self-consistent manner.

Divertor↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2024

Accelerated oxidation of epoxy thermosets with increased O 2 pressure

Polymer oxidation is usually accelerated with temperature, which is therefore applied in nearly every experimental approach dealing with predictive materials aging. Because of this simple approach, we may tend to neglect that the effective concentration of oxygen also acts as a rate multiplier for oxidation. Increasing the oxygen partial pressure in an aging environment accelerates oxidation, and while there is often a near proportional increase initially, the effect of additional oxygen usually transitions to a saturation level at some elevated pressure. This has been theoretically described in the general autoxidation scheme and is well recognized. However, for many materials the exact rate behavior under moderately increased oxygen concentration remains to be established. We therefore review epoxy oxidation and offer a broader overview of its behavior under increased O 2 partial pressure. Experimental data are given for a few thermoset materials demonstrating their rate behavior under O 2 partial pressure up to 4 atm, meaning approximately 20 times more than under standard atmospheric conditions. Confirmative evidence suggests that epoxy materials will reach saturation oxidation rates only at significantly higher O 2 partial pressure. In such a high-pressure regime it is theoretically possible to not only accelerate oxidation, but to transition into a condition where O 2 diffusion can be increased without further accelerating the oxidation rate. Finally, this can reduce diffusion limited oxidation effects under specific accelerated aging conditions as a combination of temperature and O 2 partial pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Asymptotic coefficients of the attached-eddy model derived from an adiabatic atmosphere

The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘−1’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient (A 1 ) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than 4 × 10 7 , the inferred A 1 converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘−1’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to A 1 .

atmospheric flows↗
DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2026

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗
DOE OSTI2026

Cyclotron breaking: a mechanism for parallel ion cyclotron waves to heat the fast solar wind

The Parker Solar Probe mission has observed near-continuous power in parallel ion cyclotron waves (PICWs) in the young, fast solar wind. These waves are unlikely to be directly produced by the turbulent cascade and are likely born of a local instability; yet, they are observed to both cool – and heat – the plasma. We propose that these observations can be self-consistently explained as the natural consequence of PICWs propagating in the inhomogeneous solar wind after they have been driven unstable. In this work, we argue that strong proton heating by a turbulent cascade of oblique ICWs will result in PICWs being driven unstable in a process known as quasi-linear focusing. Because the power in the turbulent cascade is concentrated at scales above the turbulent transition region, PICWs will be driven unstable within a range of wavenumbers parallel to the background magnetic field, 𝑘 ∥ , that is bounded from above by 𝑘$^{∗}_{∥P}$, corresponding to the start of the transition region. As unstable PICWs propagate away from the Sun to regions of lower proton density, their 𝑘 ∥ , multiplied by the proton inertial length 𝑑 p , increases. Eventually, 𝑘$^{∗}_{∥P}$ of the PICWs becomes larger than 𝑘$^{∗}_{∥P}$⁢𝑑 p and the waves damp, heating the solar wind. We call this effect ‘cyclotron breaking’, in analogy with ocean waves breaking on the shore. We then discuss the testable predictions of the theory, including a distinct heating signature in which PICWs cool fast protons and heat slow protons at any given heliocentric distance 𝑟. Finally, we conjecture that cyclotron breaking can lead to net heating by PICWs if the power emitted as PICWs decreases sufficiently rapidly with 𝑟 that local emission of PICWs is overwhelmed by the local damping of PICWs generated closer to the Sun.

plasma heating↗