Search NASA⌕ Search

SEARCH · Search NASA

Results for “decomposition”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Predicting river turbidity in Pine Island Bayou using machine learning techniques coupled with variational mode decomposition

Elevated turbidity levels pose significant public health risks by facilitating the transport of harmful pollutants, including metals, organic compounds, and pathogenic microorganisms into the surface water. These conditions create serious challenges for public recreational water use and drinking water treatment, leading to economic losses and health risks. This study utilizes water monitoring data in Pine Island Bayou, Texas, and develops a Sequence-to-Sequence (S2S) model to predict turbidity using Attention-based Gated Recurrent Units with Encoder-Decoder (AT-GRU-ED) and Long Short-Term Memory (LSTM), coupled with Variational Mode Decomposition (VMD). Compared to the model without VMD, the model demonstrates satisfactory 72-hour turbidity prediction performance, achieving MAEs of 2.60 and 3.29 NTU (reductions of 53% and 58%), RMSEs of 21.08 and 31.49 NTU (reductions of 82% and 80%), and R² values of 0.96 and 0.84 on the validation and test sets, respectively. Feature importance analysis reveals that water temperature is the dominant factor influencing seasonal turbidity patterns, while real-time hourly rainfall significantly contributes to short-term variability. Turbidity typically peaks within 48 hours after rainfall events due to lagged effects from surface runoff and upstream flow. Findings suggest suspending recreational water use and water supply pumping for three days after heavy rainfall can benefit public health and improve water treatment processes. Discharges above 100 m3/s are found to accelerate sediment dilution and transport, reducing turbidity levels more quickly after the peak. In conclusion, the proposed model demonstrates reliable 72-hour turbidity prediction, supporting decision-making for water treatment plant operations and providing early warning for public recreational water use.

Deep learning↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Data-Driven Analysis of Multipactor Dynamics via Dynamic Mode Decomposition

Multipactor effect is a performance-limiting kinetic plasma effect that can occur in high-power microwave and radio frequency (RF) devices. Multipactor effect is of special concern in vacuum or near-vacuum conditions such as those in particle accelerators and spaceborne devices. In this work, we present a data-driven reduced-order model (ROM) based on dynamic mode decomposition (DMD) for modeling of multipactor effects. We study multipactor effects and the resulting nonlinear harmonic generation by processing high-fidelity data generated from electromagnetic particle-in-cell (EMPIC) simulations using the DMD algorithm. We also investigate time-delay embedding extensions of DMD with improved generalizability and accuracy for modeling the electron plasma current density behavior. Here, the results show that DMD provides valuable insights into multipactor phenomena by extracting relevant modal spatiotemporal patterns and frequencies. In addition, DMD offers the potential to time extrapolate EMPIC simulations at a minimal cost, thereby reducing overall simulation time.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗
DOE OSTI2024

Computing Sparse Tensor Decompositions via Chapel and C++/MPI Interoperability without Intermediate I/O

We extend an existing approach for efficient use of shared mapped memory across Chapel and C++ for graph data stored as 1-D arrays to sparse tensor data stored using a combination of 2-D and 1-D arrays. We describe the specific extensions that provide use of shared mapped memory tensor data for a particular C++ tensor decomposition tool called GentenMPI. We then demonstrate our approach on several real-world datasets, providing timing results that illustrate minimal overhead incurred using this approach. Finally, we extend our work to improve memory usage and provide convenient random access to sparse shared mapped memory tensor elements in Chapel, while still being capable of leveraging high performance implementations of tensor algorithms in C++.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

Rapid, wide-scale deployment of green process systems, such as carbon capture or water desalination systems, is essential for combatting climate change. Methods relying on traditional design or modularity fail to capture the benefits of both economies of numbers and economies of scale. We have proposed process family design, which designs a family of processes simultaneously exploiting opportunities for common elements. In previous work, we explored different optimization formulations to solve this problem. In this work, we develop a decomposition approach to tackle larger problems efficiently. We solve a water desalination case study, which is too large to solve within a reasonable timeframe with the discretization formulation. We exploit the block angular structure of the discretization problem to decompose and solve using Progressive Hedging (PH). We use the open-source Python package mpi-sppy to execute PH which allows us to leverage parallelization and a HPC cluster to further improve solution time.

Stinchfield, Georgia↗
DOE OSTI2026

A comparative study on cubic and tetragonal Ce-ZrO 2 supported Rh catalysts for N 2 O decomposition

Zirconium oxide (ZrO 2 ) exhibits strong synergy with cerium oxide (CeO 2 ), acting as a structural and electronic promoter during catalytic redox reactions. As a result, Ce-ZrO 2 composite oxides are widely used as supports in various catalytic systems. In our previous work, we demonstrated that the incorporation of Zr 4+ into the CeO 2 lattice significantly enhanced Rh dispersion, improved redox ability, and stabilized surface Rh species, which collectively boosted the de-N 2 O activity of Rh/Ce-ZrO 2 catalysts. Building on these findings, the present study emphasizes that the crystallographic phase of Ce-ZrO 2 , governed by the Ce/Zr ratio, plays a decisive role in tuning the physicochemical environment of Rh active sites and thereby optimizing catalytic performance. In conclusion, tailoring the Ce/Zr ratio to favor the cubic fluorite structure emerges as a promising strategy for the rational design of highly active and stable catalysts for N 2 O decomposition and potentially other redox-sensitive environmental applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stress intensity factor models using mechanics-guided decomposition and symbolic regression

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. Here, in this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

42 ENGINEERING↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
DOE OSTI2024

ZeroCAL: Eliminating Carbon Dioxide Emissions from Limestone’s Decomposition to Decarbonize Cement Production

Limestone (calcite, CaCO 3 ) is an abundant and cost-effective source of calcium oxide (CaO) for cement and lime production. However, the thermochemical decomposition of limestone (~800 °C, 1 bar) to produce lime (CaO) results in substantial carbon dioxide (CO 2(g) ) emissions and energy use, i.e., ~1 tonne [t] of CO 2 and ~1.4 MWh per t of CaO produced. Here, we describe a new pathway to use CaCO 3 as a Ca source to make hydrated lime (portlandite, Ca(OH) 2 ) at ambient conditions (p, T) while nearly eliminating process CO 2(g) emissions (as low as 1.5 mol. % of the CO 2 in the precursor CaCO 3 , equivalent to 9 kg of CO 2(g) per t of Ca(OH) 2 ) within an aqueous flowelectrolysis/ pH-swing process that coproduces hydrogen (H 2(g) ) and oxygen (O 2(g) ). Because Ca(OH) 2 is a zero-carbon precursor for cement and lime production, this approach represents a significant advancement in the production of zero-carbon cement. The Zero CArbon Lime (ZeroCAL) process includes dissolution, separation/recovery, and electrolysis stages according to the following steps: (Step 1) chelator (e.g., ethylenediaminetetraacetic acid, EDTA)-promoted dissolution of CaCO 3 and complexation of Ca 2+ under basic (>pH 9) conditions, (Step 2a) Ca enrichment and separation using nanofiltration (NF), which allows separation of the Ca-EDTA complex from the accompanying bicarbonate (HCO 3 – ) species, (Step 2b) acidity-promoted decomplexation of Ca from EDTA, which allows near-complete chelator recovery and the formation of a Ca-enriched stream, and (Step 3) rapid precipitation of Ca(OH) 2 from the Ca-enriched stream using electrolytically produced alkalinity. These reactions can be conducted in a seawater matrix yielding coproducts including hydrochloric acid (HCl) and sodium bicarbonate (NaHCO 3 ), resulting from electrolysis and limestone dissolution, respectively. Careful analysis of the reaction stoichiometries and energy balances indicates that approximately 1.35 t of CaCO 3 , 1.09 t of water, 0.79 t of sodium chloride (NaCl), and ~2 MWh of electrical energy are required to produce 1 t of Ca(OH) 2 , with significant opportunity for process intensification. This approach has major implications for decarbonizing cement production within a paradigm that emphasizes the use of existing cement plants and electrification of industrial operations, while also creating approaches for alkalinity production that enable cost-effective and scalable CO 2 mineralization via Ca(OH) 2 carbonation.

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

Diffusion Model-Guided Inverse Design of Bimetallic Catalysts for Ammonia Decomposition

In the past decade, artificial intelligence and deep learning have played increasingly prominent roles in materials design and discovery. Among these, generative AI models, known for their ability to create unique and complex structures, have emerged as state-of-the-art tools for materials screening due to their high efficiency and low computational cost. In catalysis, one of the major challenges is identifying promising material candidates within an immense chemical space. This challenge can be addressed using generative approaches, such as diffusion-based inverse design models. In this study, we present a machine learning-guided workflow that employed a diffusion model for the inverse design of bimetallic alloy catalysts for low-carbon ammonia decomposition, a key reaction for ammonia emission control and sustainable hydrogen production. Catalyst candidates were evaluated using nitrogen adsorption energy as the key descriptor, inspired by multiscale modeling. The proposed workflow identified low-cost, environmentally friendly catalysts with excellent catalytic performance, which have been validated theoretically and experimentally. Our framework decoupled the generative and property-prediction components, enhancing both flexibility and accuracy in the catalytic material design process.

Adsorption↗
DOE OSTI2024

Lagrangian Decomposition of the Atlantic Ocean Heat Transport at 26.5°N

Abstract The Atlantic Meridional Overturning Circulation (AMOC) plays a critical role in the global climate system through the redistribution of heat, freshwater and carbon. At 26.5°N, the meridional heat transport has traditionally been partitioned geometrically into vertical and horizontal circulation cells; however, attributing these components to the AMOC and Subtropical Gyre (STG) flow structures remains widely debated. Using water parcel trajectories evaluated within an eddy‐rich ocean hindcast, we present the first Lagrangian decomposition of the meridional heat transport at 26.5°N. We find that water parcels recirculating within the STG account for 37% (0.36 PW) of the total heat transport across 26.5°N, more than twice that of the classical horizontal gyre component (15%). Our findings indicate that STG heat transport cannot be meaningfully distinguished from that of the basin‐scale overturning since water parcels cooled within the gyre subsequently feed the northward, subsurface limb of the AMOC.

Tooth, Oliver J.↗
DOE OSTI2024

Soil Moisture‐Cloud‐Precipitation Feedback in the Lower Atmosphere From Functional Decomposition of Satellite Observations

Abstract The feedback of topsoil moisture (SM) content on convective clouds and precipitation is not well understood and represented in the current generation of weather and climate models. Here, we use functional decomposition of satellite‐derived SM and cloud vertical profiles (CVP) to quantify the relationship between SM and the vertical distribution of cloud water in the central US. High‐dimensional model representation is used to disentangle the contributions of SM and other land‐surface and atmospheric variables to the CVP. Results show that the sign and strength of the SM‐cloud‐precipitation feedback varies with cloud height and time lag and displays a large spatial variability. Positive anomalies in antecedent 7‐hr SM and land‐surface temperature enhance cloud reflectivity up to 4 dBZ in the lower atmosphere about 1–3 km above the surface. Our approach presents new insights into the SM‐cloud‐precipitation feedback and aids in the diagnosis of land‐atmosphere interactions simulated by weather and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2024

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal decomposition of neptunyl ammonium nitrate: mechanistic insights and structural characterization of the Np 2 O 5 intermediate phase

Neptunium (Np) possesses a rich and unique chemistry that often diverges from other actinide elements yet remains relatively underexplored compared with the other light actinides. A resurgence of interest in Np has been spurred by the application of 237 Np for plutonium-238 ( 238 Pu) production for use in radioisotope thermoelectric generators (RTGs), necessitating evaluation of Np chemical reactions and materials. The work presented here studied the thermal decomposition of neptunyl ammonium nitrate (NH 4 Np VI O 2 (NO 3 ) 3 ) for synthesis of neptunium dioxide (NpO 2 ), which is the target material used for production of 238 Pu. Additionally, structural characterization of the intermediate solid Np pentoxide (Np 2 O 5 ) was performed. Advanced solid-state characterization techniques, including simultaneous thermal analysis (STA), powder X-ray diffraction (pXRD), Raman spectroscopy, and density functional theory (DFT) modeling have been combined to study the reaction pathways. Analysis revealed that NH 4 Np VI O 2 (NO 3 ) 3 thermally decomposes to a proposed neptunyl nitrate intermediate, followed by Np 2 O 5 and finally NpO 2 , all within the temperature range of 150 °C–600 °C. Further characterization of the pentoxide intermediate provided the first Raman spectra of pure-phase Np 2 O 5 and associated DFT modeling confirmed Raman peak assignments for this phase. These findings provide mechanistic information to advance production of the critical radioisotope 238 Pu and advance the state of knowledge on Np materials chemistry using modern characterization techniques.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗