Search NASA⌕ Search

SEARCH · Search NASA

Results for “low order model”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulations of Quantum Circuits with Approximate Noise using qsim and Cirq

We introduce multinode quantum trajectory simulations with qsim, an open source high performance simulator of quantum circuits. qsim can be used as a backend of Cirq, a Python software library for writing quantum circuits. We present a novel delayed inner product algorithm for quantum trajectories which can result in an order of magnitude speedup for low noise simulation. We also provide tools to use this framework in Google Cloud Platform, with high performance virtual machines in a single mode or multinode setting. Multinode configurations are well suited to simulate noisy quantum circuits with quantum trajectories. Finally, we introduce an approximate noise model for Google's experimental quantum computing platform and compare the results of noisy simulations with experiments for several quantum algorithms on Google's Quantum Computing Service.

Isakov, Sergei V.↗
DOE OSTI2023

BEYONDPLANCK XIII. Intensity foreground sampling, degeneracies, and priors

We present the intensity foreground algorithms and model employed within the BEYONDPLANCK analysis framework. The BEYONDPLANCK analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the Planck Low Frequency Instrument (LFI) data analysis. Given the scope of the BEYONDPLANCK analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous Commander component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: (1) joint foreground spectral parameter and amplitude sampling, building on ideas from MIRAMARE; (2) component-based monopole determination; (3) joint spectral parameter and monopole sampling; and (4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current BEYONDPLANCK data set is the peak frequency of the anomalous microwave emission component, for which we find $ν$ p = 25.3 ± 0.5 GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects. This work will be organized within the Open Science-based COSMOGLOBE community effort.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Block–spiral magnetism: An exotic type of frustrated order

Significance Magnetic frustration in spin or electronic models emerges from the failure of the system to fulfill simultaneously conflicting local requirements. The latter typically arise from the lattice geometry or are induced by special spin–spin interactions in the system at various distances. Here we show that the competing energy scales of the seemingly nonfrustrated orbital-selective Mott phase of the low-dimensional multiorbital Hubbard model can originate a “block–chiral magnetism,” i.e., a state with rigidly rotating spin–magnetic islands. Furthermore, we show how such an exotic spin state influences the electronic properties of the system, revealing parity-breaking quasiparticles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2021

MFNets: data efficient all-at-once learning of multifidelity surrogates as directed networks of information sources

We present an approach for constructing a surrogate from ensembles of information sources of varying cost and accuracy. The multifidelity surrogate encodes connections between information sources as a directed acyclic graph, and is trained via gradient-based minimization of a nonlinear least squares objective. While the vast majority of state-of-the-art assumes hierarchical connections between information sources, our approach works with flexibly structured information sources that may not admit a strict hierarchy. The formulation has two advantages: (1) increased data efficiency due to parsimonious multifidelity networks that can be tailored to the application; and (2) no constraints on the training data—we can combine noisy, non-nested evaluations of the information sources. Finally, numerical examples ranging from synthetic to physics-based computational mechanics simulations indicate the error in our approach can be orders-of-magnitude smaller, particularly in the low-data regime, than single-fidelity and hierarchical multifidelity approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enumeration as a Tool for Structure Solution: A Materials Genomic Approach to Solving the Cation-Ordered Structure of Na 3 V 2 (PO 4 ) 2 F 3

While powder diffraction methods are routinely utilized to optimize structural models for compounds whose crystal structures are known, the determination of unknown structures is far more challenging. When the unknown structure is large, structure solution can become a virtually intractable problem using standard structure solution methodologies, especially when the space group cannot be unambiguously resolved. One such system is the promising Na-ion battery cathode material Na 3 V 2 (PO 4 ) 2 F 3 whose high temperature and room temperature structures were previously solved, but whose more complex low-temperature structure could not be determined. Here, a novel materials genomic approach is demonstrated for the solution of the unknown 100 K structure of Na 3 V 2 (PO 4 ) 2 F 3 in which enumeration methods are first used to generate a large number (~3,000) of trial structures based on plausible orderings of Na ions and then automated Rietveld refinements are carried out to optimize each of these trial structures. Based on both the analysis of the ensemble of optimized trial structures and the density functional theory energy minimization of selected trial structures, the 100 K structure of Na 3 V 2 (PO 4 ) 2 F 3 is best described as belonging to the space group A2 1 am with unit cell dimensions of a = 9.01928(4), b = 27.1379(1), c = 10.73307(5). The 100 K unit cell has a large volume of 2627.07(2) Å 3 with Z = 12 and 33 independent crystallographic sites (9 Na, 3 V, 3 P, 12 O, and 6 F) that is 3x and 6x larger than the room- and high-temperature polymorphs of this phase, respectively. Finally, the novel methods described here will be generally applicable for the solution of the complex cation-ordered structures that commonly occur for battery materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Low ozone dry deposition rates to sea ice during the MOSAiC field campaign: Implications for the Arctic boundary layer ozone budget

Dry deposition to the surface is one of the main removal pathways of tropospheric ozone (O 3 ). We quantified for the first time the impact of O 3 deposition to the Arctic sea ice on the planetary boundary layer (PBL) O 3 concentration and budget using year-round flux and concentration observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign and simulations with a single-column atmospheric chemistry and meteorological model (SCM). Based on eddy-covariance O 3 surface flux observations, we find a median surface resistance on the order of 20,000 s m –1 , resulting in a dry deposition velocity of approximately 0.005 cm s –1 . This surface resistance is up to an order of magnitude larger than traditionally used values in many atmospheric chemistry and transport models. The SCM is able to accurately represent the yearly cycle, with maxima above 40 ppb in the winter and minima around 15 ppb at the end of summer. However, the observed springtime ozone depletion events are not captured by the SCM. In winter, the modelled PBL O 3 budget is governed by dry deposition at the surface mostly compensated by downward turbulent transport of O 3 towards the surface. Advection, which is accounted for implicitly by nudging to reanalysis data, poses a substantial, mostly negative, contribution to the simulated PBL O 3 budget in summer. During episodes with low wind speed (<5 m s –1 ) and shallow PBL (<50 m), the 7-day mean dry deposition removal rate can reach up to 1.0 ppb h –1 . Our study highlights the importance of an accurate description of dry deposition to Arctic sea ice in models to quantify the current and future O 3 sink in the Arctic, impacting the tropospheric O 3 budget, which has been modified in the last century largely due to anthropogenic activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

System-size dependence of charged-particle suppression in ultrarelativistic nucleus-nucleus collisions

High-energy partons lose energy while propagating through the hot, strongly interacting medium produced in ultrarelativistic nucleus-nucleus collisions, leading to a suppression of particle production at high transverse momentum (p T ). The dependence of this energy loss on the size of the colliding nuclear system has yet to be firmly established experimentally. This Letter presents a systematic study of charged-particle suppression across four different nucleus-nucleus collision systems using nuclear modification factors (R AA ) measured by the CMS Collaboration at the CERN LHC. Previous CMS measurements of R AA in oxygen-oxygen, xenon-xenon, and lead-lead collisions are recast with identical p T intervals and are complemented by the first measurement of the charged-particle R AA in neon-neon collisions at $\sqrt{S_{NN}}$ = 5.36 Te V . The neon-neon data correspond to an integrated luminosity of 0.76 nb −1 . The R AA in all collision systems examined show similar qualitative trends as a function of p T , but have a magnitude which is ordered with the nucleon number A. The R AA feature a downward slope at low p T , a local minimum at around 5–7 Ge V , and an upward slope with increasing p T . The R AA are also compared in terms of A 1/3 , which is proportional to the nuclear radius. Models including only initial-state nuclear effects fail to reproduce the observed trends, whereas energy loss models reproduce the trends in the region p T > 9.6 Ge V .

CMS↗
DOE OSTI2023

The effect of pellet injection on turbulent transport in TJ-II

An improvement in confinement associated with the injection of pellets has been observed in TJ-II and other devices. By means of a simple model, we study here the modification of turbulent transport by the injection of pellets in TJ-II and how this modification affects the particle confinement time. The analysis of the results shows the relationship of the confinement improvement with the evolution of the shear flow due to turbulence, especially near low order rational surfaces. This work is based on pellet injections into the neutral beam injection heated phase of TJ-II plasmas. For this, single, double, or triple hydrogen pellet injection are considered and the post-injection evolution of selected plasma parameters are analyzed and compared.

Physics↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Cold neutron-deuteron capture and Wigner-SU(4) symmetry

We calculate the cold neutron-deuteron (nd) capture cross section,σ nd to next-to-next-to leading order (NNLO) using the model-independent approach of pionless effective-field theory [EFT(π)]. At leading order we find σ nd = 0.314 ± 0.217 mb, while the experimental result is 0.508(15) mb for a laboratory neutron velocity of 2200 m/s. At next-to-leading-order (NLO), we show that σnd is sensitive to the low-energy constant (LEC) $L$$^{(0)}_{1}$ of the two-nucleon isovector current appearing at NLO. A fit of $L$$^{(0)}_{1}$ at NLO to the triton magnetic moment yields a NLO prediction of σ nd = 0.393 ± 0.164 mb, where the error comes from propagating the error from the $L$$^{(0)}_{1}$ fit. At NNLO, we find that a new three-nucleon magnetic moment counterterm is required for renormalization-group invariance of both σnd and the triton magnetic moment. Fitting the NNLO correction to $L$$^{(0)}_{1}$ (denoted $L$$^{(1)}_{1}$) to cold neutron-proton capture (σnp) yields a NNLO prediction of σ nd = 0.447 ± 0.130 mb, where the error comes from propagating the error from the $L$$^{(1)}_{1}$ fit. We also study different fittings of $L$$^{(0)}_{1}$ and $L$$^{(1)}_{1}$ to σ np , σ nd , and/or the triton magnetic moment. For example, fitting $L$$^{(0)}_{1}$ simultaneously to σ np , σ nd , and the triton magnetic moment at NLO, and fitting $L$$^{(1)}_{1}$ simultaneously to σ np and σnd at NNLO, yields σ nd = 0.480 ± 0.114 mb and 0.511 ± 0.042 mb, respectively, where errors are naively estimated from EFT(π) power counting. Additionally, we discuss how Wigner SU(4) symmetry may alter the naive EFT(π) expansion of σ nd .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Magnetic dilution in the triangular lattice antiferromagnet NaYb 1−𝑥 ⁢Lu 𝑥 ⁢O 2

The delafossitelike compound NaYbO 2 hosts a triangular lattice of Yb 3+ moments and is a promising candidate for the realization of a quantum spin liquid ground state—an exotic, quantum-disordered magnetic phase featuring long-range entanglement of spins. Tuning this system away from this quantum-disordered regime toward classical order or spin freezing is a powerful approach to shed light on the nature of the parent ground state. Here we leverage the substitution of nonmagnetic Lu 3+ onto the Yb 3+ sites to study the effects of magnetic disorder in NaYbO 2 using low-temperature ac susceptibility, heat capacity, and muon spin relaxation (𝜇⁢SR) measurements. Here, our 𝜇⁢SR measurements reveal resilient, correlated magnetic fluctuations that persist to at least 15% dilution, precluding conventional spin freezing and magnetic inhomogeneity. Heat capacity and magnetic susceptibility resolve a rapid suppression of the field-induced “up-up-down” magnetic order upon dilution and a crossover in the power-law behavior of the low-temperature magnetic excitations associated with the zero-field quantum disordered ground state. Taken together, these results support the notion of a robust network of entangled moments in NaYbO 2 and provide experimental validation of several models of a Heisenberg triangular lattice antiferromagnet in the presence of disorder.

Muon spin relaxation & rotation↗
DOE OSTI2022

SBN-BD: O(10 GeV) Proton Beam Dump at Fermilab’s PIP-II Linac

Proton beam dumps are prolific sources of mesons enabling a powerful technique to search for vector mediator coupling of dark matter to neutral pion and higher mass meson decays. By the end of the decade the PIP-II linac will be delivering up to 1 MW of proton power to the FNAL campus. This includes a significant increase of power to the Booster Neutrino Beamline (BNB) which delivers 8 GeV protons to the Short Baseline Neutrino (SBN) detectors. By building a new dedicated beam dump target station, and using the SBN detectors, a greater than an order of magnitude increase in search sensitivity for dark matter relative to the recent MiniBooNE beam dump search can be achieved. This modest cost upgrade to the BNB would begin testing models of the highly motivated relic density limit predictions and provide novel ways to test explanations of the anomalous excess of low energy events seen by MiniBooNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

SBN-BD: $\mathcal{O}$(10 GeV) Proton Beam Dump at Fermilab&#x27;s PIP-II Linac

Proton beam dumps are prolific sources of mesons enabling a powerful technique to search for vector mediator coupling of dark matter to neutral pion and higher mass meson decays. By the end of the decade the PIP-II linac will be delivering up to 1 MW of proton power to the FNAL campus. This includes a significant increase of power to the Booster Neutrino Beamline (BNB) which delivers 8 GeV protons to the Short Baseline Neutrino (SBN) detectors. By building a new dedicated beam dump target station, and using the SBN detectors, a greater than an order of magnitude increase in search sensitivity for dark matter relative to the recent MiniBooNE beam dump search can be achieved. This modest cost upgrade to the BNB would begin testing models of the highly motivated relic density limit predictions and provide novel ways to test explanations of the anomalous excess of low energy events seen by MiniBooNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Learning Optimal Aerodynamic Designs

This project created a framework for efficient, accurate, and scalable deep neural network representations of design optimization problem solutions. The inputs to these DNN representations are the vector of design requirement parameters, the outputs are the optimal design variables, and the goal is to learn the map from inputs to outputs (i.e., inverse design). The team addressed the problem of the optimal shape design of aerodynamic lifting surfaces—in particular aircraft wings—using a Reynolds-Average Navier Stokes model to govern the CFD-based aerodynamic shape optimization. The inverse design map for such problems is very complex and high-dimensional, involving inputs and outputs on the order of 1000s. To approximate this inverse design map, the team developed algorithms to construct parsimonious DNN architectures, which automatically identify low-dimensional manifolds in which design requirements affect optimal shape parameters, and trained these architectures with multifidelity optimization methods. The resulting methodology accurately and automatically designs optimal aerodynamic lifting surfaces with very high accuracy (99%) at interactive speeds, of the order of milliseconds, resulting in factors of one million or more speedup relative to CFD-based design optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Towards Online Machine Learning in DUNE Data Acquisition

Processing the large volumes of data produced by liquid argon time projection chamber (LArTPC) experiments presents a significant challenge, especially those at the scale of DUNE. This is a particular challenge when aiming to trigger on low-energy neutrinos from core-collapse supernovae, which are typically buried in a high-rate radiological background. To enable real-time event selection suitable for such rare signals, we are developing machine learning based data filtering methods. In order to demonstrate the feasibility of this approach, we implemented such pipeline using the ICEBERG detector at Fermilab as a small-scale LArTPC, with a focus on identifying Michel electrons as a proxy for low-energy neutrino interactions. This poster will present the current status of integrating these machine learning models into the data acquisition (DAQ) system of this detector.

Dalager, Olivia [Fermilab]↗