Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLE MULTIPLE METHOD”

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

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
DOE OSTI2025

Collective Modes in Multilayer Graphene/𝛼−RuCl 3 Heterostructures

Collective modes in multilayer graphene, such as plasmons and phonons, exhibit sensitivity to displacement fields and interlayer coupling, distinguishing them from their counterparts in single-layer graphene. Here, we engineer collective modes in charge-transfer heterostructures composed of multilayer graphene and 𝛼−RuCl 3 . In heterostructures with a single 𝛼−RuCl 3 interface, the charge transfer generates displacement fields up to 7 V/nm at the interface between 𝛼−RuCl 3 and the adjacent graphene layer—the highest value achieved through charge-transfer methods. As a result of the broken inversion symmetry, we discover enhanced nonlinear optical response and modified phonon selection rules. Conversely, we find that multilayer graphene sandwiched between two 𝛼−RuCl 3 flakes causes displacement fields to cancel. There, we achieve carrier densities as high as 8 ×10 13 cm −2 in multilayer graphene and restore the phonon selection rules to their unperturbed state. Meanwhile, we demonstrate that plasmonic properties derive from the depletion of multiple valence bands. As a result of the quasilinear band dispersion, these “Dirac multiband plasmons” are relatively unaffected by displacement fields. On the other hand, the inverted heterostructure sequence—two multilayer graphene sheets encapsulating 𝛼−RuCl 3 —activates significant alteration of the plasmons via interlayer plasmon-plasmon coupling. Hence, multilayer graphene and 𝛼−RuCl 3 heterostructures offer a gate-free platform for engineering collective modes derived from inversion symmetry and interlayer coupling.

Moore, Samuel L. [Columbia University, New York, N↗
DOE OSTI2025

Visualizing Stereodynamics in Cold Collisions through Shape Resonance Wavefunctions

Shape resonances in cold collisions are often strongly affected by stereodynamics. These resonances can sometimes be assigned by single-channel quantization along the scattering coordinate. However, sensitive steric control of collision implies strong anisotropy of the interaction potential energy surface, which usually leads to coupling among multiple scattering channels in the strongly interacting region. Hence, the resonances might be the result of quantization in a multidimensional space. Using the cold rotationally inelastic collision between para-H 2 (v 1 = 1, j 1 = 2) and HF (v 2 = 0, j 2 = 0) as an example, we analyze four low-lying shape resonances via diagonalization of the full-dimensional Hamiltonian with a stabilization method. While some resonances can indeed be assigned with a single partial wave, others apparently involve more than one scattering channel. Furthermore, a new model based on these resonance wavefunctions is developed to better understand the stereodynamics through shape resonances in cold collisions.

Collisions↗
DOE OSTI2025

Data Assimilation for Robust UQ Within Agent-Based Simulation on HPC Systems

Agent-based simulation provides a powerful tool for in silico system modeling. However, these simulations do not provide built-in methods for uncertainty quantification (UQ). Within these types of models a typical approach to UQ is to run multiple realizations of the model then compute aggregate statistics. This approach is limited due to the compute time required for a solution. When faced with an emerging biothreat, public health decisions need to be made quickly and solutions for integrating near real-time data with analytic tools are needed. We propose an integrated Bayesian UQ framework for agent-based models based on sequential Monte Carlo sampling. Given streaming or static data about the evolution of an emerging pathogen this Bayesian framework provides a distribution over the parameters governing the spread of a disease through a population. These estimates of the spread of a disease may be provided to public health agencies seeking to abate the spread. By coupling agent-based simulations with Bayesian modeling in a data assimilation, our proposed framework provides a powerful tool for modeling dynamical systems in silico. We propose a method which reduces model error and provides a range of realistic possible outcomes. Moreover, our method addresses two primary limitations of ABMs: the lack of UQ and an inability to assimilate data. Our proposed framework combines the flexibility of an agent-based model with UQ provided by the Bayesian paradigm in a workflow which scales well to HPC systems. We provide algorithmic details and results on a simulated outbreak with both static and streaming data.

Spannaus, Adam [ORNL] (ORCID:0000000225213657)↗
DOE OSTI2025

High-throughput micro-scale bandgap mapping for perovskite-inspired materials with complex composition space

Abstract To realize the full promise of high-throughput experimental workflows, the rate of sample synthesis must be matched by that of characterization. Of growing interest are contactless optical techniques that can rapidly measure material homogeneity and properties. Here, we present a hyperspectral imaging method to measure local optical bandgap distributions within samples, utilizing spatially-resolved reflectance spectra coupled with automated data analysis. We collect approximately one million optical bandgap data across the compositional space of Cs 3 (Bi x Sb 1-x ) 2 (Br y I 1-y ) 9 perovskite-inspired materials. Our results show non-monotonic bandgap variations (i.e., bandgap bowing) along six composition gradient sequences, in addition to identifying samples with multiple bandgaps in statistics. High-throughput transient absorption spectroscopy reveals that within these compositions, the depletion of the ground state carriers to excited states occurred at discrete energy levels with independent carrier dynamics, consistent with the bandgap observation and indicative of phase separation. This work demonstrates the potential for rapid optical measurements to assess material quality and homogeneity in a high-throughput experimental setting, supporting screening and recipe optimization of optoelectronic material candidates with desired carrier dynamics and optical properties.

Science & Technology - Other Topics↗
DOE OSTI2025

Emergent Nonthermal Fluid from Jets in the Massive Schwinger Model Using Tensor Networks

We analyze the correlation between the energy, momentum, and spatial entanglement produced by two luminal jets in the massive Schwinger model. Using tensor network methods, we show that for m/g > 1/π, in the vicinity of the strong- to weak-coupling transition, a nearly perfect and chargeless effective fluid behavior appears around the midrapidity region with a universal energy-pressure relationship. The evolution of energy and pressure is strongly correlated with the rise of the spatial entanglement entropy, indicating a key role of quantum dynamics. Some of these observations may be used to analyze high multiplicity jet fragmentation events, energy-energy and energy-charge correlators at current collider energies.

Ab initio calculations↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Three-particle formalism for multiple channels: the ηππ + $ K\overline{K}\pi $ system in isosymmetric QCD

We generalize previous three-particle finite-volume formalisms to allow for multiple three-particle channels. For definiteness, we focus on the two-channel ηππ and $ K\overline{K}\pi $ system in isosymmetric QCD, considering the positive G parity sector of the latter channel, and neglecting the coupling to modes with four or more particles. The formalism we obtain is thus appropriate to study the b 1 (1235) and η(1295) resonances. The derivation is made in the generic relativistic field theory approach using the time-ordered perturbation theory method. We study how the resulting quantization condition reduces to that for a single three-particle channel when one drops below the upper ($ K\overline{K}\pi $) threshold. We also present parametrizations of the three-particle K matrices that enter into the formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Using the optimal combined index weight ratio to improve the probability of anomaly detection in big area additive manufacturing

Big Area Additive Manufacturing (BAAM) of composites requires significant time, energy, and material, so it is critical to reduce production inefficiencies to make functional parts without multiple iterations. Statistical process control coupled with Principal Component Analysis (PCA) is a powerful technique that provides a quick, computationally inexpensive, and intuitive way for operators to detect defects that form in a manufacturing process without massive datasets. Recently, a combined index that is a weighted sum of the Hotelling's T 2 and squared residual error statistics has been proposed that can be monitored in one chart, improving interpretation accuracy and simplicity. However, the literature does not offer a formal method to optimise the weights. Here, we introduce two new approaches to the traditional weight selection approach using simulated and BAAM image data. Approach 1 uses a theoretically motivated optimum inspired by probabilistic principal component analysis. Approach 2 systematically varies the ratio of the weights to find the optimum. We show that approach 1 delivers optimal anomaly detection performance in select cases while approach 2 fares better in practice. Surprisingly, we also show that choosing a more complex PCA model has a minimal negative impact on anomaly detection performance compared to a more simplistic model.

3-dimensional printing↗
DOE OSTI2026

Monte Carlo Event Generation with Continuous Normalizing Flows

We apply continuous normalizing flows trained with the flow matching method to the problem of phase-space sampling in Monte Carlo event generation for high-energy collider physics. Focusing on lepton-pair and top-quark pair production with multiple jets, the two computationally most expensive processes at the Large Hadron Collider, we train helicity-conditioned continuous normalizing flows to remap the random numbers used in matrix element evaluation. Compared to standard methods, we achieve unweighting efficiency improvements by factors of up to 184 and 25 for the two processes at their respective highest jet number, at the cost of an increased evaluation time. When combining the advantages of continuous normalizing flows with the fast evaluation times of coupling-layer-based flows, using the RegFlow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest jet numbers. These substantial gains highlight the promise of samplers based on machine learning for next-generation collider experiments.

Bothmann, Enrico [CERN; Gottingen U.] (ORCID:00000↗
DOE OSTI2026

Test of Harmonic Coil Accuracy Using a Multi Vertex Probe

Magnets are integral components of particle accelerators, with different multipoles used to perform various beam conditioning functions, like steering and focusing. Testing accelerator magnets to ensure compliance with specification is a crucial step of fabrication. One of the most popular techniques for verifying magnetic field strength and integrity is the harmonic coil probe. Harmonic coil probes are comprised of passive loops that are rotated around the aperture of a magnet and are used to measure magnetic flux. Oftentimes the measurements of the harmonic coil probe are affected by transverse vibrations along the axis of rotation and torsional vibrations of the position closure sensor, both causing spurious harmonics. To make matters challenging, it is difficult to distinguish between the real harmonics of the magnet and the spurious harmonics during analysis. This research explores the use of a novel configuration of harmonic coil probe, Multi-Vertex Probe (MVP), to address multiple factors that affect the accuracy and repeatability of harmonic coil probe measurement. Most of the results are based on simulation data, in which artificial torsional distortions are introduced and removed with the proposed processing method. In addition to implementing torsional vibration removal in simulation, laboratory analysis regarding the source of vibrations was performed using a pre-fabricated MVP. More specifically, various mechanical conditions were used to explore the standard deviation of main harmonic’s flux amplitude. The conditions included: coupling type (rigid or soft), rotation speed, probe starting position within magnetic field, and cabling location inside the magnetic field. Along with this, the MVP’s measured noise floor in and outside a magnetic field was quantified.

Lofquist, Claire [Northern Illinois U.]↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2025

Enhancing Security and Resiliency in Operational Technology Environments Through Network Slicing and Federated Learning

The growing convergence of Information Technology (IT) and Operational Technology (OT) within Industry 4.0 environments has introduced new demands on industrial network infrastructure. As cyber-physical systems become increasingly interconnected, ensuring the secure, timely, and efficient exchange of critical data is essential. This thesis explores how network slicing, a method of creating isolated virtual network segments, can be applied within OT environments to address challenges such as latency, security, and resource allocation. The first research question addressed in this thesis is: How can OT networks take advantage of NFV and SDN technology to become cyber resilient? This study examines the operational, security, and architectural implications of introducing network slicing into traditionally static OT infrastructures such as Industrial Control Systems (ICS) and SCADA. Through simulated deployments and case studies, the research demonstrates how slicing enables better isolation between critical and non-critical services, thereby improving response time, throughput, and security in sensitive environments. The second question considers: How to dynamically implement network slicing and take advantage of network resources towards integrating decentralized machine learning? In response, this thesis proposes a framework that combines Software-Defined Networking (SDN), Network Function Virtualization (NFV), and Federated Learning (FL) to enable real-time analytics while maintaining data locality. The proposed approach reduces the burden on centralized infrastructure and minimizes privacy risks by supporting on-site training of models across distributed OT nodes, coordinated through dynamically allocated network slices. The third focus explores: How slicing helps to increase the resiliency of OT networks through the orchestration of a dynamic DMZ? To answer this, the thesis presents a method for creating and managing Dynamic Demilitarized Zones (DMZs) using network slicing. This enables flexible and automated isolation of sensitive subsystems during threat scenarios or high-risk operations. Coupled with intelligent orchestration and containerized security services, the dynamic DMZ significantly enhances the system's ability to respond to cyber incidents without halting production. Ultimately, this thesis contributes a comprehensive architecture that blends network slicing with machine learning, secure segmentation, and automation, paving the way for resilient, adaptive, and intelligent OT environments. Performance evaluations across multiple scenarios show improvements in system reliability, threat response time, model accuracy, and resource utilization, providing a strong foundation for future industrial automation systems.

Rodiles Delgado, Brian G↗
DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗
DOE OSTI2025

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗