Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive 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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of Atmosphere and Ocean Horizontal Model Resolution on Tropical Cyclone and Upper-Ocean Response Forecasts in Four Major Hurricanes

A coupled atmosphere–ocean model is necessary for tropical cyclone (TC) prediction to accurately characterize ocean feedback on atmospheric processes within the TC environment. Here, the ECMWF coupled global model is run at horizontal resolutions from 9 to 1.4 km in the atmosphere, as well as 25 and 8 km in the ocean, to identify how resolution impacts forecast accuracy of four observed major TCs in the Atlantic: Irma, Florence, Teddy, and Ida. Most of the resolutions used here are unprecedented for global models. GOES-16 and synthetic aperture radar (SAR) satellite images and best track data are used for atmospheric validation. Salinity and temperature observations from Air-Launched Autonomous Micro-Observer (ALAMO) floats are used to validate modeled upper-ocean response, including mixed layer deepening, sea surface cooling, and near-inertial waves in the wakes of TCs. Increasing atmospheric resolution leads to more realistic TC structure and stronger winds, significantly improving TC intensity forecasts and modestly improving track errors. Ocean resolution impacts the upper-ocean response but does not influence atmospheric forecasts for the fast-moving TCs considered here. Stronger mixing, sea surface cooling, and near-inertial oscillations are found for both higher atmosphere and ocean resolutions, provided the initial upper-ocean state is the same for the two ocean resolutions. Whether this agrees better with the ALAMO observations also depends on the realism of the initial upper-ocean state in the model, emphasizing the importance of ocean initialization for the accurate upper-ocean response. Overall, the model at all resolutions correctly predicts stronger mixing, surface cooling, and near-inertial oscillation amplitudes to the right of a TC center, as observed by ALAMO floats.

Atmosphere-ocean interaction↗
DOE OSTI2024

Role of bath-induced many-body interactions in the dissipative phases of the Su-Schrieffer-Heeger model

The Su-Schrieffer-Heeger chain is a prototype example of a symmetry-protected topological insulator. Coupling it nonperturbatively to local thermal environments, either through the intercell or the intracell fermion tunneling elements, modifies the topological window. To understand this effect, we employ the recently developed reaction-coordinate polaron transform (RCPT) method, which allows treating system-bath interactions at arbitrary strengths. The effective system Hamiltonian, which is obtained via the RCPT, exposes the impact of the baths on the SSH chain through renormalization of tunneling elements and the generation of many-body interaction terms. By performing exact diagonalization and computing the ensemble geometric phase, a topological invariant, which is applicable even to systems at finite temperature, we distinguish the trivial band insulator (BI) from the topological insulator (TI) phases. Furthermore, through the RCPT mapping, we are able to pinpoint the main mechanism behind the extension of the parameter space for the TI or the BI phases (depending on the coupling scheme, intracell or intercell), which is the bath-induced, dimerized, many-body interaction. In conclusion, we also study the effect of on-site staggered potentials on the SSH phase diagram and discuss extensions of our method to higher dimensions.

1-dimensional systems↗
DOE OSTI2021

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

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

Transition to the Haldane phase driven by electron-electron correlations

One of the most famous quantum systems with topological properties, the spin S = 1 antiferromagnetic Heisenberg chain, is well-known to display exotic S = 1/2 edge states. However, this spin model has not been analyzed from the more general perspective of strongly correlated systems varying the electron electron interaction strength. Here, we report the investigation of the emergence of the Haldane edge in a system of interacting electrons – the two orbital Hubbard model—with increasing repulsion strength U and Hund interaction J H . We show that interactions not only form the magnetic moments but also form a topologically nontrivial fermionic many-body ground-state with zero-energy edge states. Specifically, upon increasing the strength of the Hubbard repulsion and Hund exchange, we identify a sharp transition point separating topologically trivial and nontrivial ground-states. Surprisingly, such a behaviour appears already at rather small values of the interaction, in a regime where the magnetic moments are barely developed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Investigation of the reaction mechanism of the hydrodeoxygenation of propionic acid over a Rh(1 1 1) surface: A first principles study

Microkinetic models based on first principles calculations have been used to study the vapor and liquid phase hydrodeoxygenation of propionic acid on a Rh(111) surface. Calculations suggest that both decarboxylation and decarbonylation do not occur at an appreciable rate under all reaction environments. Propanol and propionaldehyde are the main products on this surface and they are produced at similar rates in both vapor and liquid phase environments. While a condensed phase can shift the reaction rate, the dominant pathways and selectivity towards the various products are hardly affected. At 473 K, the turnover frequency is increased by about a factor of 1.5 in liquid water relative to the gas phase. In liquid 1,4-dioxane, the turnover frequency is also slightly increased relative to the gas phase. Here, given the uncertainty in Rh cavity radius in the liquid phase calculations, computations with different cavity radii have been performed. With larger Rh cavity radius, the promotional effect of the solvents on the turnover frequency becomes more significant, while practically no changes are observed for a smaller cavity radius.

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

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Correlation-assisted quantized charge pumping

Here, we investigate charge pumping in the vicinity of order-obstructed topological phases, i.e., symmetry-protected topological phases masked by spontaneous symmetry breaking in the presence of strong correlations. To explore this, we study a prototypical Su-Schrieffer-Heeger model with finite-range interaction that gives rise to orbital charge density wave order and characterize the impact of this order on the model's topological properties. In the ordered phase, where the many-body topological invariant loses quantization, we find that not only is quantized charge pumping still possible, but it is even assisted by the collective nature of the orbital charge density wave order. Remarkably, we show that the Thouless pump scenario may be used to uncover the underlying topology of order-obstructed phases.

1-dimensional systems↗
DOE OSTI2025

Neutrino-portal dark matter detection prospects at a future muon collider

With no concrete evidence for nongravitational interactions of dark matter to date, it is natural to wonder whether dark matter couples predominantly to the Standard Model (SM)’s neutrinos. Neutrino interactions (and the possible existence of additional neutrinophilic mediators) are substantially less understood than those of other SM particles, yet this picture will change dramatically in the coming decades with new neutrino sources. One potential new source arises with the construction of a high-energy muon collider (MuCol); due to muons’ instability, a MuCol is a source of high-energy collimated neutrinos. Importantly, since the physics of muon decays (into neutrinos) is very well-understood, this leads to a neutrino flux with systematic uncertainties far smaller than fluxes from conventional high-energy (proton-sourced) neutrino beams. In this work, we study the capabilities of a potential neutrino detector, “ MuCol ν ”, placed ∼ 100 m downstream of the MuCol interaction point. The MuCol ν detector would be especially capable of searching for a neutrinophilic mediator ϕ through the mononeutrino scattering process ν μ N → μ + ϕ X , exceeding searches from other terrestrial approaches for m ϕ in the ∼ few MeV- ten GeV range. Even with a 10 kg-yr exposure, MuCol ν is capable of searching for well-motivated classes of thermal freeze-out and freeze-in neutrino-portal dark matter. Published by the American Physical Society 2025

Adhikary, Jyotismita (ORCID:0000000253894170)↗
DOE OSTI2022

Differential methods for assessing sensitivity in biological models

Differential sensitivity analysis is indispensable in fitting parameters, understanding uncertainty, and forecasting the results of both thought and lab experiments. Although there are many methods currently available for performing differential sensitivity analysis of biological models, it can be difficult to determine which method is best suited for a particular model. In this paper, we explain a variety of differential sensitivity methods and assess their value in some typical biological models. First, we explain the mathematical basis for three numerical methods: adjoint sensitivity analysis, complex perturbation sensitivity analysis, and forward mode sensitivity analysis. We then carry out four instructive case studies. (a) The CARRGO model for tumor-immune interaction highlights the additional information that differential sensitivity analysis provides beyond traditional naive sensitivity methods, (b) the deterministic SIR model demonstrates the value of using second-order sensitivity in refining model predictions, (c) the stochastic SIR model shows how differential sensitivity can be attacked in stochastic modeling, and (d) a discrete birth-death-migration model illustrates how the complex perturbation method of differential sensitivity can be generalized to a broader range of biological models. Finally, we compare the speed, accuracy, and ease of use of these methods. We find that forward mode automatic differentiation has the quickest computational time, while the complex perturbation method is the simplest to implement and the most generalizable.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Modeling Framework for Data Center

This chapter highlights the critical need for advanced modeling of data centers due to their rapidly increasing energy consumption and impact on grid reliability. Driven by the demand for AI applications, data centers are projected to consume a significant portion of US energy by 2028, putting stress on an already challenged power grid. The chapter emphasizes the importance of "fast" time-scale models to understand the dynamic interactions between data centers and the grid, especially given the rapid power fluctuations of AI workloads. It outlines a modeling framework that includes both offline and real-time EMT domain simulations, detailing the necessary representations for various components like utility interfaces, transformers, IT loads, UPS, cooling loads, Battery Energy Storage Systems (BESS), generators, protection systems, and higher-level control systems. While standard simulation tools like PSCAD offer basic models, custom development is often required to accurately capture the unique and fast-changing behaviors of modern data centers. The chapter also discusses key metrics and test cases for validating these models, focusing on transient load responses, protection relay coordination, and demand flexibility. Finally, it addresses the challenges of modeling large-scale data centers, such as computational complexity and the trade-off between model fidelity and practicality, suggesting hybrid modeling approaches as a solution. The overarching goal is to create a robust framework that helps assess data center impacts on grid stability, identify vulnerabilities, and inform the development of standards for reliable integration of these large loads into the bulk power system.

25 ENERGY STORAGE↗
DOE OSTI2025

The Atlantic Meridional Overturning Circulation’s Response to CO2 Increase: Assessing the Roles of Surface Flux and Oceanic Advection Feedbacks

Abstract The Atlantic meridional overturning circulation (AMOC) is projected to slow down in climate models due to greenhouse gas emissions. However, the physical mechanisms determining the rate of the projected AMOC slowdown remain unclear. Accordingly, this study isolates the roles of oceanic advection and surface flux feedbacks that might accelerate or decelerate the AMOC’s weakening using carbon dioxide (CO 2 ) quadrupling simulations in the CESM1.2 model. Surface flux feedbacks are isolated in partially coupled experiments in which either all surface flux components or the momentum flux responses to AMOC’s weakening that might provide feedback are suppressed, while a tracer decomposition of ocean density anomalies isolates the advection feedbacks. Comparing the ocean density components in the experiments shows that the AMOC’s response is initially determined by CO 2 -induced anomalous surface heat fluxes, and afterward, feedbacks determine its response. In the fully coupled case, surface heat flux feedback strongly promotes AMOC slowdown and causes its near shutdown, while a weaker but active AMOC is maintained when the surface flux feedback is inhibited in the partially coupled case. The positive surface heat flux feedback works by canceling out the negative oceanic heat advection feedback on deep-water formation in the subpolar North Atlantic (SPNA). With the heat advection feedback thus reduced, the positive salinity advection feedback becomes the dominant contributor to SPNA density changes and deep-water formation. In the partially coupled case, negative ocean heat advection feedback and the CO 2 -induced subtropical Atlantic saline anomalies imported into the SPNA play a stabilizing role. The results highlight the importance of SPNA salinity gradients and gyre circulation strength in determining the AMOC’s response rate or recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Implementation of INCL Nuclear Model in GENIE Generator

The Li ge Intranuclear Cascade (INCL) model is a nuclear-physics model that simulates hadron (baryon, anti-baryon and meson) reactions on nuclei, for incident energies ranging from a few tens of MeV to 10-20 GeV. The INCL model has been well validated by hadron scattering data. In my work, I implement an interface in GENIE to use the INCL nuclear model in the simulations of both the initial state of the target nucleus and the Final State Interaction in neutrino-nucleus interaction. It has a consistent treatment of nuclear models in both neutrino interaction and hadron rescattering. A full event record including neutrino vertex and each vertex of hadron rescattering has been accomplished. Several processes, e.g. cluster production, Delta transportation and de-excitation will also be included as benefits of the implementation of the INCL model in GENIE. I will show some initial simulation results showcasing the new GENIE features and discuss plans for making them available for use in experimental analyses.

Liu, Liang [Fermilab]↗
DOE OSTI2025

Implementation of INCL nuclear model in GENIE Generator

The Liège Intranuclear Cascade (INCL) model is a nuclear-physics model that simulates hadron (baryon, anti-baryon and meson) reactions on nuclei, for incident energies ranging from a few tens of MeV to 10-20 GeV. The INCL model has been well validated by hadron scattering data. In my work, I implement an interface in GENIE to use the INCL nuclear model in the simulations of both the initial state of the target nucleus and the Final State Interaction in neutrino-nucleus interaction. It has a consistent treatment of nuclear models in both neutrino interaction and hadron rescattering. A full event record including neutrino vertex and each vertex of hadron rescattering has been accomplished. Several processes, e.g. cluster production, Delta transportation and de-excitation will also be included as benefits of the implementation of the INCL model in GENIE. I will show some initial simulation results showcasing the new GENIE features and discuss plans for making them available for use in experimental analyses.

Liu, Liang [Fermilab] (ORCID:000000026753925X)↗
DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Active learning emulators for nuclear two-body scattering in momentum space

In this work we extend the active learning emulators for two-body scattering in coordinate space with error estimation, recently developed by Maldonado et al. [Phys. Rev. C 112, 024002], to coupled-channel scattering in momentum space. Our full-order model (FOM) solver is based on the Lippmann-Schwinger integral equation for the scattering t-matrix as opposed to the radial Schrödinger equation. We use (Petrov-)Galerkin projections and high-fidelity calculations at a few snapshots across the parameter space of the interaction to construct efficient reduced-order models (ROMs), trained by a greedy algorithm for locally optimal snapshot selection. Both the FOM solver and the corresponding ROMs are implemented efficiently in Python using Google's JAX library. We present results for emulating scattering phase shifts in coupled and uncoupled channels and cross sections, and assess the accuracy of the developed ROMs and their computational speedup factors. We also develop emulator error estimation for both the t-matrix and the total cross section. The software framework for reproducing and extending our results is publicly available. Together with our recent advances in developing active-learning emulators for three-body scattering, these emulator frameworks set the stage for full Bayesian calibrations of chiral nuclear interactions and optical models against scattering data with quantified emulator errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Superconductivity of incoherent electrons in the Yukawa Sachdev-Ye-Kitaev model

We study a model of N fermions in a quantum dot, coupled to M bosons by a disorder-induced complex Yukawa coupling [Yukawa Sachdev-Ye-Kitaev (SYK) model], in order to explore the interplay between non-Fermi liquid and superconductivity in a strongly coupled, (quantum-)critical environment. We analyze the phase diagram of the model for an arbitrary complex interaction and arbitrary ratio of N/M, with special focus on the two regimes of non-Fermi-liquid behavior: an SYK-like behavior with a power-law frequency dependence of the fermionic self-energy and an impuritylike behavior with frequency independent self-energy. We show that the crossover between the two can be reached by varying either the strength of the fermion-boson coupling or the ratio M/N. We next argue that in both regimes the system is unstable to superconductivity if the strength of time-reversal-symmetry-breaking disorder is below a certain threshold. We show how the corresponding onset temperatures vary between the two regimes. We argue that the superconducting state is highly unconventional with an infinite set of minima of the condensation energy at T = 0, corresponding to topologically different gap functions. We discuss in detail similarities and differences between this model and the model of dispersion-full fermions tuned to a metallic quantum-critical point, with an effective singular dynamical interaction V(Ω) ∝ 1/|Ω| γ (the γ model).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗