Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear Models”

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

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Circularization of 23S rRNA but not 16S rRNA within archaeal ribosomes

Background Processing of archaeal 16S and 23S rRNAs is believed to involve excision of individual rRNAs from polycistronic precursors, circularization of excised rRNAs, and re-linearization before the incorporation into ribosomes. However, all the knowledge is derived from several isolated species, leaving open the possibility that different processes may occur in other archaeal groups. Results Here, we investigate rRNAs from diverse and mostly uncultivated archaea. Sequencing of total cellular RNA from eight phylum-level lineages indicates that archaeal circular 23S rRNA transcript abundances vastly exceed those of linear counterparts, and linear versions are often undetectable. As the majority of rRNAs derive from mature ribosomes, the data suggest that ribosomes contain circular 23S rRNAs. Thus, we directly sequence RNA extracted from isolated ribosomes of a model archaeon, Methanosarcina acetivorans, and confirm that the 23S rRNAs in the ribosomes are circular. Structural modeling places the 5′ and 3′ ends of the linear precursors of archaeal 23S rRNAs in close proximity to form a GNRA tetraloop (in which N is A, C, G, or U and R is A or G), consistent with their existence as circular molecules. We also confirm the existence of circular 16S rRNA intermediates in transcriptomes of most archaea, yet a circular form is not evident in some distinct archaeal groups, suggesting that certain archaea do not circularize 16S rRNA during processing. Conclusions Our findings uncover unexpected variations in the processing required to generate mature rRNAs and the conformation of functional molecules in archaeal ribosomes.

Archaea↗
DOE OSTI2025

Minimum GHG emissions and energy consumption of U.S. PET and polyolefin packaging supply chains in a circular economy

There is a wide agreement on the urgency of transforming linear management of plastics towards a circular economy model. However, no clear pathways exist as to required recycling technologies involved and system-wide environmental impacts. This study explores such pathways in the U.S. for the most commonly used packaging plastics through a combination of mechanical and emerging advanced recycling technologies. A system optimization model aimed at minimizing environmental impacts was developed to determine optimal end-of-life (EOL) management and locations of existing and emerging U.S. recycling infrastructures. Our study includes material flows from virgin resin production through semi-manufacturing processes to existing EOL disposal and recycling processes. An optimized circular plastics packaging system achieved greenhouse gas (GHG) emission savings of up to 28% and cumulative energy demand (CED) savings of up to 46%, compared to the linear economy. Moreover, these savings of GHG emissions and CED impacts represent a reduction of 0.16% and 0.49% compared to annual U.S. GHG emissions and energy consumption in 2022, respectively. The optimal recycling rates and systems-level circularity ranged from 78–99% and 57–75%, respectively. Increased energy savings led to increased GHG emissions showing a potential trade-off between GHG emissions and energy. Analysis of 40 scenarios showed the importance of material collection distances, blend limit of mechanically recycled resins, process yields, and mandated recycling rates for achieving a sustainable circular economy of plastics.

09 - BIOMASS FUELS↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2025

Effect of vector meson spin coherence on the observables for the chiral magnetic effect in heavy-ion collisions

The chiral magnetic effect (CME) in heavy-ion collisions reflects the local violation of P and CP symmetries in strong interactions and manifests as electric charge separation along the direction of the magnetic field created by the wounded nuclei. The experimental observables for the CME, such as the γ 112 correlator, the R Ψ$_2$ ⁡ (Δ⁢S) correlator, and the signed balance functions, however, are also subject to non-CME backgrounds, including those from resonance decays. A previous study showed that the CME observables are affected by the diagonal component of the spin density matrix, the ρ 00 for vector mesons. Here, in this work, we study the contributions from the other elements of the spin density matrix using a toy model and a multiphase transport model. We find that the real part of the ρ 1-1 component, Re ⁡ρ 1-1 , affects the CME observables in a manner opposite to that of the ρ 00 . All three aforementioned CME observables show a linear dependence on Re ⁡ρ 1-1 in the model calculations, supporting our analytical derivations. The rest elements of the spin density matrix do not contribute to the CME observables. The off-diagonal terms in the spin density matrix indicate spin coherence and may be nonzero in heavy-ion collisions due to local spin polarization or spin-spin correlations. Thus, Re ⁡ρ 1-1 , along with ρ 00 , could play a significant role in interpreting measurements in search of the CME.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗
DOE OSTI2024

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗
DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗
DOE OSTI2024

Kernel polynomial method for linear spin wave theory

Calculating dynamical spin correlations is essential for matching model magnetic exchange Hamiltonians to momentum-resolved spectroscopic measurements. A major numerical bottleneck is the diagonalization of the dynamical matrix, especially in systems with large magnetic unit cells, such as those with incommensurate magnetic structures or quenched disorder. In this paper, we demonstrate an efficient scheme based on the kernel polynomial method for calculating dynamical correlations of relevance to inelastic neutron scattering experiments. This method reduces the scaling of numerical cost from cubic to linear in the magnetic unit cell size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2025

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Boosting Barlow Twins Reduced Order Modeling for Machine Learning‐Based Surrogate Models in Multiphase Flow Problems

Abstract We present an innovative approach called boosting Barlow Twins reduced order modeling (BBT‐ROM) to enhance the reliability of machine learning surrogate models for multiphase flow problems. BBT‐ROM builds upon Barlow Twins reduced order modeling that leverages self‐supervised learning to effectively handle linear and nonlinear manifolds by constructing well‐structured latent spaces of input parameters and output quantities. To address the challenge of high contrast data in multiphase flow problems due to injection wells and faults, we employ a boosting algorithm within BBT‐ROM. This algorithm sequentially trains a set of weak models (i.e., inaccurate models), improving prediction accuracy through ensemble learning. To evaluate the performance of BBT‐ROM, we conduct three three‐dimensional multiphase flow problems, including waterflooding and geologic carbon storage (GCS), with varying numbers of input parameter cases and model domain features. The results demonstrate that BBT‐ROM excels at predicting non‐wetting phase saturation (e.g., oil or saturation) and fluid pressure, with average relative errors ranging from 0.5% to 3%. Importantly, BBT‐ROM showcases robustness when faced with limited input parameter space during GCS testing.

58 GEOSCIENCES↗
DOE OSTI2025

Digital Real-Time Simulation and Power Quality Analysis of a Hydrogen-Generating Nuclear-Renewable Integrated Energy System

This paper investigates the challenges and solutions associated with integrating a hydrogen-generating nuclear-renewable integrated energy system (NR-IES) under a transactive energy framework. The proposed system directs excess nuclear power to hydrogen production during periods of low grid demand while utilizing renewables to maintain grid stability. Using digital real-time simulation (DRTS) in the Typhoon HIL 404 model, the dynamic interactions between nuclear power plants, electrolyzers, and power grids are analyzed to mitigate issues such as harmonic distortion, power quality degradation, and low power factor caused by large non-linear loads. A three-phase power conversion system is modeled using the Typhoon HIL 404 model and includes a generator, a variable load, an electrolyzer, and power filters. Active harmonic filters (AHFs) and hybrid active power filters (HAPFs) are implemented to address harmonic mitigation and reactive power compensation. The results reveal that the HAPF topology effectively balances cost efficiency and performance and significantly reduces active filter current requirements compared to AHF-only systems. During maximum electrolyzer operation at 4 MW, the grid frequency dropped below 59.3 Hz without filtering; however, the implementation of power filters successfully restored the frequency to 59.9 Hz, demonstrating its effectiveness in maintaining grid stability. Future work will focus on integrating a deep reinforcement learning (DRL) framework with real-time simulation and optimizing real-time power dispatch, thus enabling a scalable, efficient NR-IES for sustainable energy markets.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗
DOE OSTI2025

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Velocity Dependence of Rate‐And‐State Friction in Granular Fault Gouge and Implications for Slow‐Slip Events

Abstract The standard rate‐and‐state friction (RSF) has extensively captured frictional behaviors, but it fails to explain the velocity dependence of frictional stability transition and widespread slow‐slip events (SSEs) in experiments and nature adequately. An alternative microphysical Chen‐Niemeijer‐Spiers (CNS) model can well describe the velocity dependence of frictional behaviors of granular gouges. Using the original CNS model, standard RSF parameters can be quantified microphysically. However, some micro‐parameters are not easy to estimate quantitatively, making it difficult to extrapolate to natural and experimental conditions. Here, we simplify the microphysically‐derived RSF parameters including direct effect a , evolution effect b , and critical slip distance D c , as well as equivalent values ( a eq , b eq , and D eq ). The simplified friction parameters directly illustrate their velocity dependence, namely the essentially constant a , a eq , and D c , negatively velocity‐dependent b and b eq , as well as varying D eq for different laws. They are roughly consistent with experimental results in various fault gouges. A modified CNS model is further derived from the original CNS model, establishing a direct link between the standard RSF and CNS models. The modified CNS model exhibits virtually identical frictional behaviors to the original CNS, but differs from the standard RSF at large velocity perturbations. Moreover, the linearized stability analysis indicates that the critical stiffness for the modified CNS model is velocity‐dependent. Compared with the standard RSF, the modified CNS model not only explains the velocity dependence of frictional stability transition, but also exhibits a more gradual transition for SSEs with a broader range of stiffness ratios.

Geochemistry & Geophysics↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗