SEARCH · Search NASA

Results for “multidimensional”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗

DOE OSTI2025

Study of the muon component in the core-corona model using CONEX 3D

The discrepancy between models and data regarding the muon content in air showers generated by ultra-high energy cosmic rays still needs to be solved. The CONEX simulation framework provides a flexible tool to assess the impact of different interaction properties and thus address the muon puzzle. In this work, we present the multidimensional extension of CONEX and show its performance compared to CORSIKA by discussing muon-related air-shower features for three experiments: KASCADE, IceTop, and the Pierre Auger Observatory. We also implement an effective version of the core-corona model to demonstrate the impact of the core effect, as observed at the LHC, on the muon content in air showers produced by ultra-high energy cosmic rays. At a primary energy of E$_{0}$ = 10$^{19}$ eV, we obtain an increase of 15% to 20% in the muon content.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

AI-assisted detector design for the EIC (AID(2)E)

Artificial Intelligence is poised to transform the design of complex, large-scale detectors like ePIC at the future Electron Ion Collider. Featuring a central detector with additional detecting systems in the far forward and far backward regions, the ePIC experiment incorporates numerous design parameters and objectives, including performance, physics reach, and cost, constrained by mechanical and geometric limits. This project aims to develop a scalable, distributed AI-assisted detector design for the EIC (AID(2)E), employing state-of-the-art multiobjective optimization to tackle complex designs. Supported by the ePIC software stack and using G EANT 4 simulations, our approach benefits from transparent parameterization and advanced AI features. The workflow leverages the PanDA and iDDS systems, used in major experiments such as ATLAS at CERN LHC, the Rubin Observatory, and sPHENIX at RHIC, to manage the compute intensive demands of ePIC detector simulations. Tailored enhancements to the PanDA system focus on usability, scalability, automation, and monitoring. Ultimately, this project aims to establish a robust design capability, apply a distributed AI-assisted workflow to the ePIC detector, and extend its applications to the design of the second detector (Detector-2) in the EIC, as well as to calibration and alignment tasks. Additionally, we are developing advanced data science tools to efficiently navigate the complex, multidimensional trade-offs identified through this optimization process.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Trade can buffer climate-induced risks and volatilities in crop supply

Climate change is intensifying the frequency and severity of extreme events, posing challenges to food security. Corn, a staple crop for billions, is particularly vulnerable to heat stress, a primary driver of yield variability. While many studies have examined climate impact on average corn yields, little attention has been given to the climate impact on production volatility. This study investigates the future volatility and risks associated with global corn supply under climate change, evaluating the potential benefits of two key adaptation strategies: irrigation and market integration. A statistical model is employed to estimate corn yield response to heat stress and utilize NEX-GDDP-CMIP6 climate data to project future production volatility and risks of substantial yield losses. Three metrics are introduced to quantify these risks: Sigma (σ), the standard deviation of year-on-year yield change, which reflects overall yield volatility; Rho (ρ), the risk of substantial loss, defined as the probability of yield falling below a critical threshold; and Beta (β), a relative risk coefficient that captures the volatility of a region's corn production compared to the globally integrated market. The analysis reveals a concerning trend of increasing year-on-year yield volatility (σ) across most regions and climate models. This volatility increase is significant for key corn-producing regions like Brazil and the United States. While irrigated corn production exhibits a smaller rise in volatility, suggesting irrigation as a potential buffer against climate change impacts, it is not a sustainable option as it can cause groundwater depletion. On the other hand, global market integration reduces overall volatility and market risks significantly with less sustainability concerns. Furthermore, these findings highlight the importance of a multidimensional approach to adaptation in the food sector. While irrigation can benefit individual farmers, promoting global market integration offers a broader solution for fostering resilience and sustainability across the entire food system.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2024

SDSS-IV MaNGA: stellar rotational support in disc galaxies versus central surface density and stellar population age

ABSTRACT We investigate how the stellar rotational support changes as a function of spatially resolved stellar population age ($D_n4000$) and relative central stellar surface density ($\Delta \Sigma _1$) for MaNGA isolated/central disc galaxies. We find that the galaxy rotational support indicator $\lambda _{R_\mathrm{e}}$ varies smoothly as a function of $\Delta \Sigma _1$ and $D_n4000$. $D_n4000$ versus $\Delta \Sigma _1$ follows a ‘J-shape’, with $\lambda _{R_\mathrm{e}}$ contributing to the scatters. In this ‘J-shaped’ pattern rotational support increases with central $D_n4000$ when $\Delta \Sigma _1$ is low but decreases with $\Delta \Sigma _1$ when $\Delta \Sigma _1$ is high. Restricting attention to low-$\Delta \Sigma _1$ (i.e. large-radius) galaxies, we suggest that the trend of increasing rotational support with $D_n4000$ for these objects is produced by a mix of two different processes, a primary trend characterized by growth in $\lambda _{R_\mathrm{e}}$ along with mass through gas accretion, on top of which disturbance episodes are overlaid, which reduce rotational support and trigger increased star formation. An additional finding is that star-forming galaxies with low $\Delta \Sigma _1$ have relatively larger radii than galaxies with higher $\Delta \Sigma _1$ at fixed stellar mass. Assuming that these relative radii rankings are preserved while galaxies are star forming then implies clear evolutionary paths in central $D_n4000$ versus $\Delta \Sigma _1$. The paper closes with comments on the implications that these paths have for the evolution of pseudo-bulges versus classical bulges. The utility of using $\rm D_n4000$–$\Delta \Sigma _1$ to study $\lambda _{R_\mathrm{e}}$ reinforces the notion that galaxy kinematics correlate both with structure and with stellar-population state, and indicates the importance of a multidimensional description for understanding bulge and galaxy evolution.

Wang, Xiaohan (ORCID:0009000283121002)↗

DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗

DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗

DOE OSTI2026

Boosted decision tree reweighting of simulated neutrino interactions for O ( 1 ) GeV neutrino cross-section measurements

This paper illustrates a generic method for multidimensional reweighting of O ( 1 ) GeV neutrino interaction Monte Carlo samples. The reweighting is based on a boosted decision tree algorithm trained on high-dimensional space in detector final-state observables. This enables one generator’s events to be reweighted so that its reconstructed particle content and kinematics distributions, as well as detector efficiency, match those of a target model. The approach establishes an efficient way to reuse legacy Monte Carlo data, avoiding regeneration. As an example, we test its use in a measurement of transverse kinematic imbalance of the μ - and proton in charged-current quasielastic like ν μ events from the MINERvA experiment.

Lin, Z. [Rochester U.] (ORCID:0009000188903698)↗

DOE OSTI2026

Enhancements in Laser-Direct-Drive Nuclear Performance with Target Radius

Inertial confinement fusion is likely to require significant improvements in technology and design for large targets to implode at high driver energies. Here, to assess the potential for benefits, we report on nuclear performance as a function of target radius 𝑅 t in direct-drive cryogenic implosions as performed on the OMEGA laser. The neutron yield and areal density are found to increase as 𝑅 5.0±0.2 t and 𝑅 1.8±0.2 t , respectively, and exhibit a much stronger dependence on target radius than previous studies have assumed. As this paper demonstrates, these types of sensitivities should be expected when implosions are unstable, and have been degraded by a range of multidimensional effects in 3D. If the imperfections in a target and laser system are fixed in magnitude, it follows that larger implosions have higher relative quality, and approach criteria to ignite much more quickly than commonly appreciated.

inertial confinement fusion↗

DOE OSTI2025

Neural-network extraction of unpolarised transverse-momentum-dependent distributions

We present the first extraction of transverse-momentum-dependent distributions of unpolarized quarks from experimental Drell-Yan data using neural networks to parametrize their nonperturbative part. We show that neural networks outperform traditional parametrizations providing a more accurate description of data. This Letter establishes the feasibility of using neural networks to explore the multidimensional partonic structure of hadrons and paves the way for more accurate determinations based on machine-learning techniques.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Nonlocal Coherent Optical Nonlinearities of a Macroscopic Quantum System

The optical responses of solids are typically understood to be local in space. Whether locality holds for the optical response of a macroscopic quantum system has, however, been largely unexplored. Here, we use multidimensional coherent spectroscopy at the optical diffraction limit to demonstrate nonlocal optical nonlinearities in a semiconductor microcavity. In conclusion, these nonlocal optical responses are both coherent and quantum in nature, deriving from the macroscopic length scale of confined exciton-polariton wave functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

High-precision phase control of an optical lattice with up to 50 dB noise suppression

An optical lattice is a periodic light crystal constructed from the standing-wave interference patterns of laser beams. It can be used to store and manipulate quantum degenerate atoms and is an ideal platform for the quantum simulation of many-body physics. A principal feature is that optical lattices are flexible and possess a variety of multidimensional geometries with modifiable band structure. An even richer landscape emerges when control functions can be applied to the lattice by modulation of the position or amplitude with Floquet driving. However, the desire of realizing high-modulation bandwidths while preserving extreme lattice stability has been difficult to achieve. In this paper, we demonstrate an effective solution that consists of overlapping two counterpropagating lattice beams and controlling the phase and intensity of each with independent acousto-optic modulators. Our phase controller mixes sampled light from both lattice beams with a common optical reference. This dual heterodyne locking method allows exquisite determination of the lattice position, while also removing parasitic phase noise accrued as the beams travel along separate paths. Here, we report up to 50 dB suppression in lattice phase noise in the 0.1–1 Hz band, along with significant suppression spanning more than 4 decades of frequency. The absolute phase diffusion of the lattice position is only 10 Å when integrated over 10 s. This method permits precise, high-bandwidth modulation (greater than 50 kHz) of the optical lattice intensity and phase. We demonstrate the efficacy of this approach by executing intricate time-varying phase profiles for atom interferometry.

Atom interferometry↗

DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Asymptotic-state prediction for fast flavor transformation in neutron star mergers

Neutrino flavor instabilities appear to be omnipresent in dense astrophysical environments, thus presenting a challenge to large-scale simulations of core-collapse supernovae and neutron star mergers (NSMs). Subgrid models offer a path forward, but require an accurate determination of the local outcome of such conversion phenomena. Focusing on “fast” instabilities, related to the existence of a crossing between neutrino and antineutrino angular distributions, we consider a range of analytical mixing schemes, including a new, fully three-dimensional one, and also introduce a new machine learning (ML) model. We compare the accuracy of these models with the results of several thousands of local dynamical calculations of neutrino evolution from the conditions extracted from classical NSM simulations. Our ML model shows good overall performance, but struggles to generalize to conditions from a NSM simulation not used for training. The multidimensional analytic model performs and generalizes even better, while other analytic models (which assume axisymmetric neutrino distributions) do not have reliably high performances, as they notably fail as expected to account for effects resulting from strong anisotropies. As a result, the ML and analytic subgrid models extensively tested here are both promising, with different computational requirements and sources of systematic errors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Loop-string-hadron approach to SU(3) lattice Yang-Mills theory: Hilbert space of a trivalent vertex

The construction of gauge-invariant states of SU(3) lattice gauge theories has garnered new interest in recent years, but implementing them is complicated by the need for SU(3) Clebsch-Gordon coefficients. In the loop-string-hadron (LSH) approach to lattice gauge theories, the elementary excitations are strictly gauge invariant, and constructing the basis requires no knowledge of Clebsch-Gordon coefficients. Originally developed for SU(2), the LSH formulation was recently generalized to SU(3), but limited to one spatial dimension. In this work, we generalize the LSH approach to constructing the basis of SU(3) gauge-invariant states at a trivalent vertex—the essential building block to multidimensional space. A direct generalization from the SU(2) vertex yields a legitimate basis; however, in certain sectors of the Hilbert space, the naive LSH basis vectors so defined suffer from being nonorthogonal. The issues with orthogonality are directly related to the “missing label” or “outer multiplicity” problem associated with SU(3) tensor products and may also be phrased in terms of Littlewood-Richardson coefficients or the need for a “seventh Casimir” operator. The states that are unaffected by the problem are orthonormalized in closed form. For the sectors that are afflicted, we discuss the nonorthogonal bases and their orthogonalization. A few candidates for seventh Casimir operators are readily constructed from the suite of LSH gauge-singlet operators. The diagonalization of a seventh Casimir represents one prescriptive solution toward obtaining a complete orthonormal basis, but a closed-form general solution remains to be found. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗