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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stacking disorder in α – RuCl 3 investigated via x-ray three-dimensional difference pair distribution function analysis

The van der Waals layered magnet α-RuCl 3 offers tantalizing prospects for the realization of Majorana quasiparticles. Efforts to understand this are, however, hampered by inconsistent magnetic and thermal transport properties likely coming from the formation of structural disorder during crystal growth, postgrowth processing, or upon cooling through the first order structural transition. Here, we investigate structural disorder in α-RuCl 3 using x-ray diffuse scattering and three-dimensional difference pair distribution function (3D-ΔPDF) analysis. We develop a quantitative model that describes disorder in α-RuCl 3 in terms of rotational twinning and intermixing of the high and low-temperature structural layer stacking. As a result, this disorder may be important to consider when investigating the detailed magnetic and electronic properties of this widely studied material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Magnetic properties of metastable honeycomb KCoAsO 4

We present comprehensive neutron scattering data on polycrystalline samples of a new metastable honeycomb material KCoAsO 4 . Below 𝑇 𝑁 = 14 K, the system orders into a zigzag antiferromagnetic state, with spins ordered into alternating ferromagnetic chains similar to the isostructural sister compound KNiAsO 4 . In the case of KCoAsO 4 , we find the moments are robustly canted out of plane closer to the crystallographic 𝑐 axis. A combination of weak interlayer coupling, lattice strain, and inhomogeneities lead to the coexistence of two magnetic ordering wave vectors 𝑘 1 = (1.5 0 0) and 𝑘 2 = (0.5 0 0.5), where the two structures differ only in their layer stacking. Inelastic data show the presence of a spin orbital mode at 24 meV, supporting a pseudospin $\tilde{𝑆}$ = 1/2 Kramer's doublet ground state of the Co 2+ ions. We model the low energy excitations using both a conventional XXZ Hamiltonian and generalized Kitaev Heisenberg Hamiltonian within a linear spin wave limit. While either model can qualitatively reproduce the observed spectra, the lack of fine features and observation of disorder prevent a clear-cut determination of the low energy Hamiltonian. In the case of an XXZ-type model, a large easy-axis anisotropy is necessary to reproduce the gapped spectra and canting of the magnetic moments. For the generalized Kitaev model, despite the large canting of the moments away from the honeycomb layers, we find a noticeable if nondominant Kitaev term persists. In conclusion, the contrast of the magnetic properties of KCoAsO 4 to other cobalt honeycombs highlights the sensitivity of the low energy magnetic properties of Co 2+ to fine details of its crystalline environment.

Honeycomb lattice↗
DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Remarks on the thesis of Aitor Bengoechea Fernández: Generalization of Doppler effect processing using S(α,β) thermal diffusion laws

The thesis titled “Generalization of Doppler effect processing using S(α,β) thermal diffusion laws” addresses a critical gap in the modeling of neutron scattering phenomena in nuclear materials, particularly where crystal lattice effects intersect with low-energy resonances. Conducted within the Doctoral School ED 352 at Aix-Marseille Université and hosted at CEA Cadarache, the research presents a theoretical and numerical framework for extending the applicability of Thermal Scattering Laws (TSLs) using higher-order correlation functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Thermodynamic Constraints on Electromicrobial Protein Production

Global consumption of protein is projected to double by the middle of the 21st century. However, protein production is one of the most energy intensive and environmentally damaging parts of the food supply system today. Electromicrobial production technologies that combine renewable electricity and CO 2 -fixing microbial metabolism could dramatically increase the energy efficiency of commodity chemical production. Here we present a molecular-scale model that sets an upper limit on the performance of any organism performing electromicrobial protein production. We show that engineered microbes that fix CO 2 and N 2 using reducing equivalents produced by H 2 -oxidation or extracellular electron uptake could produce amino acids with energy inputs as low as 64 MJ kg −1 , approximately one order of magnitude higher than any previous estimate of the efficiency of electromicrobial protein production. This work provides a roadmap for development of engineered microbes that could significantly expand access to proteins produced with a low environmental footprint.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Solving Critical Problems of the Muon Collider Higgs Factory: Optics, Magnets and their Protection, Detector Backgrounds

A low-energy medium-luminosity Muon Collider (MC) is being studied as a possible Higgs Factory (HF). Electrons from muon decays will deposit more than 300 kW in superconducting magnets of the HF collider ring. This imposes significant challenges to superconducting (SC) magnets used in the MC storage ring (SR) and interaction regions (IR). Magnet designs are proposed which provide high operating gradient and magnetic field in a large aperture to accommodate the large size of muon beams (due to low beta*), as well as a cooling system to intercept the large heat deposition from the showers induced by decay electrons. The distribution of heat deposition in the MC SR lattice elements requires large-aperture magnets in order to accommodate thick high-Z absorbers to protect the SC coils. Based on the developed MARS15 model and intensive simulations, a sophisticated radiation protection system was designed for the collider SR and IR to bring the peak power density in the superconducting coils below the quench limit and reduce the dynamic heat deposition in the cold mass by a factor of 100. The system consists of tight tungsten masks in the magnet interconnect regions and elliptical tungsten liners in the magnet aperture optimized individually for each magnet. These also reduce the background particle fluxes in the collider detector.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Interface-Driven Chiral Magnetism in Ultrathin Metallic Ferromagnets: Towards Skyrmion Spintronics

Magnetic skyrmions are localized chiral magnetic textures in the form of nanoscale vortices or bubbles that are topologically protected from being ‘unwound’. Their topological nature gives rise to rich behaviors including ordered lattice formation, emergent electrodynamics and robust current-driven displacement at remarkably low current densities. However, magnetic skyrmions have so far been restricted to just a few materials and observed only at low temperatures, limiting the experimental accessibility and technological application of these unique topological objects. This project aimed to realize magnetic skyrmions at room temperature in a new class of engineered chiral ferromagnets that exploit broken mirror symmetry at interfaces to generate helical magnetic order. Ultrathin magnetic heterostructures were engineered and patterned into laterally-confined nanostructures in which skyrmions were stabilized, manipulated, and detected for the first time under ambient conditions. Static and dynamic properties were examined using advanced x-ray imaging to provide a fundamental understanding of topological magnetic phases. The interactions between magnetic skyrmions and electron charge and spin currents were studied to identify fundamental new physics and enable device applications. These experimental studies, supported by detailed modeling, will help launch a new subfield of skyrmion-based spintronics in which individual skyrmions can be used to encode, store, and transport information in high-performance, low power memory and logic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Production of pions, kaons, and protons as a function of the relative transverse activity classifier in pp collisions at $ \sqrt{s} $ = 13 TeV

The production of π ± , K ± , and (p¯)p is measured in pp collisions at √s = 13 TeV in different topological regions of the events. Particle transverse momentum (p T ) spectra are measured in the “toward”, “transverse”, and “away” angular regions defined with respect to the direction of the leading particle in the event. While the toward and away regions contain the fragmentation products of the near-side and away-side jets, respectively, the transverse region is dominated by particles from the Underlying Event (UE). The relative transverse activity classifier, R T = N T / < N T >, is used to group events according to their UE activity, where N T is the measured charged-particle multiplicity per event in the transverse region and < N T > is the mean value over all the analysed events. The first measurements of identified particle p T spectra as a function of R T in the three topological regions are reported. It is found that the yield of high transverse momentum particles relative to the R T -integrated measurement decreases with increasing R T in both the toward and the away regions, indicating that the softer UE dominates particle production as R T increases and validating that R T can be used to control the magnitude of the UE. Conversely, the spectral shapes in the transverse region harden significantly with increasing R T . This hardening follows a mass ordering, being more significant for heavier particles. Finally, it is observed that the pT-differential particle ratios (p+p¯)/(π + +π – ) and (K + + K – )/(π + + π – ) in the low UE limit (R T → 0) approach expectations from Monte Carlo generators such as PYTHIA 8 with Monash 2013 tune and EPOS LHC, where the jet-fragmentation models have been tuned to reproduce e + e – results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Forward neutrinos from charm at the Large Hadron Collider

The currently operating FASER experiment and the planned Forward Physics Facility (FPF) will detect a large number of neutrinos produced in proton-proton collisions at the LHC. In this work, we estimate neutrino fluxes at these detectors from charm meson decays, which will be particularly important for the ν e and ν τ channels. We make prediction using both the next-to-leading order collinear factorization and the k T -factorization approaches to model the production of charm quarks as well as different schemes to model their hadronization into charm hadrons. In particular, we emphasize that a sophisticated modeling of hadronization involving beam remnants is needed for predictions at FASER and FPF due to the sensitivity to the charm hadron production at low transverse momenta and very forward rapidity. As example, we use the string fragmentation approach implemented in 8. While both standard fragmentation functions and 8 are able to describe LHCb data, we find that 8 predicts significantly higher rate of high energy neutrinos, highlighting the importance of using the correct hadronization model when making predictions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Spin-strain interactions under hydrostatic pressure in α-RuCl 3

We investigate the effects of hydrostatic pressure on 𝛼−RuCl 3 , a prototypical material for the Kitaev spin model on a honeycomb lattice with a possible spin-liquid ground state. Using ultrasound measurements at pressures up to 1.16 GPa, we reveal significant modifications of the acoustic properties and the 𝐻−𝑇 phase diagram of this material. Hydrostatic pressure suppresses the three-dimensional magnetic order and induces a dimerization transition at higher pressures. At low pressures, the sound attenuation exhibits a linear temperature dependence, while above 0.28 GPa, it becomes nearly temperature independent, suggesting a shift in the phonon scattering regime dominated by Majorana fermions. These findings provide new insights into spin-strain interactions in Kitaev magnets and deliver a detailed characterization of the 𝐻−𝑇 phase diagram of 𝛼−RuCl 3 under hydrostatic pressure.

Focused ion beam↗
DOE OSTI2023

Modified Virtual Oscillator Based Operation of Grid Forming Converters with Single Voltage Sensor

Virtual oscillator based control architecture for grid forming converter with inner current control loop for additional objective accomplishments is proposed in this paper. In general, virtual oscillator is based on Lienard type of oscillator equation with a cubic nonlinearity. Virtual oscillator alone cannot accommodate unbalanced or harmonic mitigation while operating converters in grid forming mode. Therefore, to accomplish better tracking performance, incorporating the provision for nonlinear and/or unbalanced loads, an inner current control loop is utilized which is based upon Lyapunov energy function type control architecture due to its advantages. To reduce the overall cost of implementation and ensure proper phase difference in the generated voltages, the implemented virtual oscillator operates with only one voltage feedback from the point of common coupling as presented in this paper. To verify the effectiveness of the proposed approach, the overall system has been modeled in MATTLAB/Simulink and PLECS domain and case studies showing the successful production of low harmonic load voltages under nonideal loading condition along with other important results are presented in this paper.

dual second order generalized integrator↗
DOE OSTI2022

Data used for Figure 3 of the Nature Reviews Earth and Environment (NREE) paper: "A low-to-no snow future and its impacts on water resources in the western United States

In order to synthesize western United States snowpack projections, this dataset contains the results of 18 peer-review journal articles over three periods of interest (2025-2049, 2050-2074, and 2075-2099) and over 4 mountain ranges (Cascades, Sierra Nevada, Rockies, Wasatch/Uinta) in addition to western-US wide projections. Distinction is made by model type (Earth System Models, bias-corrected statistically downscaled Earth System Models, and regional climate models). RCP4.5 and RCP8.5 emission scenarios are considered. Percent snow water equivalent (SWE) loss considers 1 April, peak SWE and/or seasonal SWE. Heterogeneity in projected snowpack changes exists across mountain ranges and for different modeling approaches, but generally indicate agreement in decreases by the end of the century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Space-Time Reduced-Order Modeling for Uncertainty Quantification

This work focuses on the space-time reduced-order modeling (ROM) method for solving large-scale uncertainty quantification (UQ) problems with multiple random coefficients. In contrast with the traditional space ROM approach, which performs dimension reduction in the spatial dimension, the space-time ROM approach performs dimension reduction on both the spatial and temporal domains, and thus enables accurate approximate solutions at a low cost. We incorporate the space-time ROM strategy with various classical stochastic UQ propagation methods such as stochastic Galerkin and Monte Carlo. Numerical results demonstrate that our methodology has significant computational advantages compared to state-of-the-art ROM approaches. By testing the approximation errors, we show that there is no obvious loss of simulation accuracy for space-time ROM given its high computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2020

A control-oriented combustion model framework for compression ignition engines operating on low-reactivity fuel

This work focuses on zero-dimensional modeling of the heat release rate in a compression ignition engine operating on gasoline-like fuels. Due to the properties of gasoline, such as high volatility and longer ignition delay than diesel, the injection strategies can vary significantly from the operation with conventional diesel fuel. Different injection strategies are commonly used to achieve varying degrees of in-cylinder stratification in order to shape the combustion event and maximize efficiency. The proposed zero-dimensional combustion model was developed to account for the different stages in combustion caused by the fuel stratification. As the ignition delay model is an integral part of the entire combustion process and significantly affects the prediction accuracy, special attention has been paid to local phenomena influencing ignition delay. A one-dimensional spray model by Musculus and Kattke was employed in conjunction with a Lagrangian tracking approach in order to estimate the local air–fuel ratio within the spray tip, as a proxy for reactivity. The local air–fuel ratio, in-cylinder temperature and pressure were used in an integral fashion to estimate the ignition delay. Heat release rates were modeled using first-order non-linear differential equations. The proposed combustion model was validated against experimental data of a heavy-duty compression ignition engine with up to three injection events at mostly 1038 r/min and 14 bar brake mean effective pressure. Further validation of the model was carried out at other engine loads and speeds. Model prediction errors in CA50 of less than 1 °CA across all conditions were found. Modeling results of other combustion metrics such as combustion duration and indicated mean effective pressure are also highly satisfactory. In addition, the model has been shown to be capable of estimating the ringing intensity for most conditions.

Pamminger, Michael↗
DOE OSTI2025

Plausible indication of gamma-ray absorption by dark matter in NGC 1068

NGC 1068 is the brightest extragalactic source in high-energy neutrinos as seen by IceCube, yet the accompanying gamma-ray flux is orders of magnitude weaker. It has been argued that this indicates that the bulk of neutrinos and gamma rays are emitted in the innermost vicinity of the central supermassive black hole, which is transparent to neutrinos, but opaque to gamma rays. Even in such extreme scenarios for the acceleration of cosmic rays, astrophysical models typically overestimate the low-energy gamma-ray flux and/or require some fine-tuning in the physical parameters. Here, in this work, we suggest instead that the dark matter surrounding the supermassive black hole may absorb the gamma rays, inducing the observed deficit. We show that for a dark matter–photon scattering cross section in the range σ DM−γ /m DM ≃ 10 −28 −10 −30 cm 2 /GeV, Fermi-LAT measurements can be well reconciled with IceCube data. We also present some simple particle physics examples that achieve the correct spectral energy dependence while respecting complementary constraints.

Active Galactic Nuclei↗