Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluid System”

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 91 records · Page 5

Power Electronics Thermal Management

The 2017 Electrical and Electronics Technical Team Roadmap [11] proposes aggressive research and development targets aimed at improving power electronics technology to enable the mass-market penetration of electric-drive vehicles. Achieving these aggressive targets will require a decrease in cost (year 2025 cost target: $2.70/kW) and an increase in power density (year 2025 power density target: 100 kW/L) as compared with current on-road technology. Replacing traditional silicon device-based components with more efficient and higher-temperature wide-bandgap (WBG) semiconductor device-based components will enable increased power density. However, meeting the power density target will also require innovative thermal management solutions to increase the heat fluxes dissipated and allow for compact electronics packaging. This project evaluates, designs, and develops thermal management strategies that use dielectric fluid (single-phase heat transfer) as coolants.

ADVANCED PROPULSION SYSTEMS↗

Hierarchical Network Partitioning for Solution of Potential-Driven, Steady-State Nonlinear Network Flow Equations

The solution of potential-driven steady-state flow in large networks is a task which manifests in various engineering applications, such as transport of natural gas or water through pipeline networks. The resultant system of nonlinear equations depends on the network topology, and in general, there is no numerical algorithm that offers guaranteed convergence to the solution (assuming a solution exists). Some methods offer guarantees in cases where the network topology satisfies certain assumptions, but these methods fail for larger networks. On the other hand, the Newton-Raphson algorithm offers a convergence guarantee if the starting point lies close to the (unknown) solution. It would be advantageous to compute the solution of the large nonlinear system through the solution of smaller nonlinear sub-systems wherein the solution algorithms (Newton-Raphson or otherwise) are more likely to succeed. Here, this letter proposes and describes such a procedure, a hierarchical network partitioning algorithm that enables the solution of large nonlinear systems corresponding to potential-driven steady-state network flow equations.

42 ENGINEERING↗

Shape-anisotropy inverses the behavior of emergent vortices in active chiral fluids

Active colloidal fluids exhibit spontaneous emergence of correlated states, characterized by complex collective dynamics and self-organization. In geometrically confined systems, activity modulations trigger robust polar state reversals of a macroscopic vortex formed by colloidal rollers. Here, we reveal that the shape anisotropy of dielectric rollers transforms the system into a chiral active fluid comprised of rollers of arbitrary handedness. The chiral rollers self-assemble into multiple freestanding vortices with a spontaneously selected sense of rotation. We demonstrate that upon reactivation of the system after a complete cessation of activity beyond all relevant timescales, the vortices simultaneously restore their previous chiral states in striking contrast to the chiral state reversals demonstrated by spherical rollers. The analysis reveals that shape-anisotropy modifies the collective state memory of the ensemble. The findings provide insights into the emergence of complex collective behavior in chiral colloidal fluids governed by an interplay between shape anisotropy, chiral motion, and activity modulations.

Colloids↗

Combined TDLAS and chemiluminescence imaging in a flat flame burner operated with NH3/H2 blends

Ammonia is a promising hydrogen carrier due to its favorable storage and transport characteristics. However, its direct use in combustion systems is limited by low flammability and potential for high nitrogen oxide emissions. To better understand ammonia combustion, researchers conducted experiments using a flat flame burner and measured species profiles using tunable-diode-laser-absorption-spectroscopy and chemiluminescence imaging. They tested three flame conditions with different ammonia-hydrogen blends and oxygen levels. The results were compared to simulations using various kinetic mechanisms, including one that accounts for excited species chemistry. The goal is to provide direct information about species profiles in a simple system, isolating chemical kinetics from fluid dynamic effects, which can inform the development of more efficient and low-emission combustion systems using ammonia.

ammonia combustion↗

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiphysics Co-Optimization Design and Analysis of Double-Side Cooled Silicon Carbide-Based Power Module: Preprint

With the rapid growth of Electric Vehicles (EVs) and Hybrid Electric Vehicles (HEVs), much more rigorous design targets have been set for automotive power electronics, including high power density, high reliability, and low cost. Novel power module and inverter technologies based on wide bandgap (WEG) semiconductors have been developed to meet these design targets, while providing optimal power semiconductor operating temperature and promising thermomechanical performance. Compared with conventional cooling techniques which are normally applied only on one side of power module, double-side cooling approach is now believed to be the solution to enable high power density and low thermal resistance of WEG semiconductor-based power electronics. In this work, we develop a three-phase power module that is double-sided cooled using dielectric fluid jet impingement. In each phase, four silicon carbide (SiC) power semiconductors are bonded to copper busbars without electrical insulation layers. A finite element analysis (FEA) model is created for thermal and thermomechanical analysis. Based on FEA modeling results, we select particular dimensions for a parametric study to optimize thermal and mechanical performance. Using a multi-objective genetic algorithm (MOGA)-based optimization method, we have minimized the maximum junction temperature and thermal stresses within the power module. The multiphysics co-optimization approach has enabled an efficient design process of power modules with greatly reduced computational cost, as compared to conventional processes that rely on exhaustive numerical simulations and iterations.

ADVANCED PROPULSION SYSTEMS↗

Historical and Future Windstorms in the Northeastern United States

Large-scale windstorms represent an important atmospheric hazard in the Northeastern US (NE) and are associated with substantial socioeconomic losses. Regional simulations performed with the Weather Research and Forecasting (WRF) model using lateral boundary conditions from three Earth System Models (ESMs: Geophysical Fluid Dynamics Laboratory (GFDL), Hadley Centre Global Environment Model (HadGEM) and Max Planck Institute (MPI)) are used to quantify possible future changes in windstorm characteristics and/or changes in the parent cyclone types responsible for windstorms. WRF nested within MPI ESM best represents important aspects of historical windstorms and the cyclone types responsible for generating windstorms compared with a reference simulation performed with the ERA-Interim reanalysis for the historical climate. The spatial scale and frequency of the largest windstorms in each simulation defined using the greatest extent of exceedance of local 99.9th percentile wind speeds (U > U999) plus 50-year return period wind speeds (U50,RP) do not exhibit secular trends. Projections of extreme wind speeds and windstorm intensity/frequency/geolocation and dominant parent cyclone type associated with windstorms vary markedly across the simulations. Only the MPI nested simulations indicate statistically significant differences in windstorm spatial scale, frequency and intensity over the NE in the future and historical periods. This model chain, which also exhibits the highest fidelity in the historical climate, yields evidence of future increases in 99.9th percentile 10 m height wind speeds, the frequency of simultaneous U > U999 over a substantial fraction (5–25%) of the NE and the frequency of maximum wind speeds above 22.5 ms−1. These geophysical changes, coupled with a projected doubling of population, leads to a projected tripling of a socioeconomic loss index, and hence risk to human systems, from future windstorms.

Pryor, Sara C. (ORCID:0000000348473440)↗

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗

Computational Thermal Hydraulics of a High-Performance Low-Enriched-Uranium Annular Target for HFIR Irradiation

Molybdenum-99 has historically been generated via isolation from fissioned highly enriched uranium (HEU) targets. Here, this isotope is in high demand due to its daily use across the world in radiopharmaceutical medical procedures. The primary objective of this work was to design and analyze an experimental target assembly containing one low-enriched-uranium (LEU) annular target for irradiation at the High Flux Isotope Reactor (HFIR). Efforts included incorporating spatially dependent energy sources from neutron and gamma interactions, quantifying thermal contact conductance at material interfaces, performing grid-independent studies, comparing turbulence models, and simulating various steady-state and transient scenarios relevant for irradiation qualification and eventual insertion. These models provide velocity, pressure, and temperature distributions in both space and time. Such results enable the selection of an appropriate irradiation location, fission rate density, and flow-limiting orifice size and demonstrate compliance with HFIR safety requirements such that insertion into the reactor can be approved. This analysis shows that across all scenarios, wetted surface temperatures remain below the coolant saturation temperature with no net vapor formation in the coolant. In every scenario, all components stay below 30% of the aluminum 6061 melting temperature. Computational fluid dynamics and system-level models predict peak target temperatures that agree within 4%, though the predicted axial location of the peak differs by about 10% of the heated length due to differences in flow development length. These results de-risk the irradiation of LEU (annular targets) and strengthen a domestic, HEU-independent 99 Mo supply by providing important fuel performance data to form the foundation for a robust licensing basis.

Molybdenum-99↗

Report on Initial Sodium Testing on the Thermal Hydraulic Experimental Test Article (THETA) (Fiscal Year 2024 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that is used to develop sodium components and instrumentation as well as to acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility simulates nominal thermal hydraulic conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field may be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility was designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary coolant and secondary coolant system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. In fiscal year 2023, thermal stratification tests were completed with the primary system online, while the secondary system was being constructed [1]. These tests had shown that the core barrel and intermediate heat exchanger (IHX) outlet required increased thermal insulation. The THETA primary system was removed from METL, cleaned, thermal insulators installed, and then inserted into METL Test Vessel 4. At the time of this writing the THETA primary and secondary system are operational. During this fiscal year 100+ hours of testing was completed to characterize thermal hydraulic phenomena associated with steady state and transient conditions in a pool type liquid metal cooled reactor. A majority of the testing campaign was completed to satisfy the experimental data acquisition requirements for the GAIN Voucher with Oklo, CRADA 2021-21121. THETA is still operational at the time of this publication and future testing is planned for fiscal year 2025. Work is underway to publish existing and future data to an online database to facilitate collaboration with SFR engineers looking to validate their systems code or computational fluid dynamics models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Editorial: Subsurface microbiology within hydrocarbon resources or stored gases

A Research Topic on the microbiology of hydrocarbon and gas storage reservoirs has far reaching industrial applications. In recent decades, there has been a growing interest in understanding microbial communities in subsurface energy reservoirs, such as coal, oil, and shale beds. This area of research has broadened to include gas storage reservoirs for hydrogen and CO 2 . Scientists are beginning to unravel the unexpected impact microorganisms have on these systems, through changing the fluid geochemistry, the gas content, and even the permeability. By recognizing the influence of these tiny organisms on our engineered environments, we can develop better risk assessments, target mitigation strategies, expand energy production, and refine operational guidance, ultimately contributing to a more sustainable energy future.

59 BASIC BIOLOGICAL SCIENCES↗

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A buffering heat exchanger/thermal energy storage system for desalination applications

A buffering heat exchanger/thermal energy storage (BHXTES) system was designed using a fast-response phase-change material/graphite foam (PCM/GF) medium for desalination and potentially other industrial applications. Unlike a store-now-and-use-later TES system, this study focused on developing a dynamically balanced, continuously used, system that optimized three functions. A lab-scale prototype was designed and fabricated with a PureTemp 151/GF storage medium. The experimental data showed excellent thermal performance repeatability indicating minimal effects of PCM expansion and contraction on the GF as well as negligible PCM redistribution effects. The experimental data also served as a validation of the numerical model of the BHXTES system, after which optimization of the BHXTES system was conducted through numerical simulations for both PureTemp 151/GF and solar salt/GF as the storage medium. The results show that, with a Therminol 55 heat transfer fluid, continuous operation of the system can be realized with a solar energy source through cycles of 8-h changing for desalination and energy storage with the supplied heat followed by 16-h discharging for desalination with the stored energy. As a result, the working fluid outlet temperatures are generally in the ranges of approximately 128–160 °C for a PureTemp 151/GF medium and 198–220 °C for a solar salt/GF medium.

Buffering↗

Experimental Investigation of Steel-Borne Acoustic Pulses for Fault Pinpointing in Pipe-Type Cable Systems: A Scaled-Down Model Approach

Pipe-type cable systems, including high-pressure fluid-filled (HPFF) and high-pressure gas-filled cables, are widely used for underground high-voltage transmission. These systems consist of insulated conductor cables within steel pipes, filled with pressurized fluids or gases for insulation and cooling. Despite their reliability, faults can occur due to insulation degradation, thermal expansion, and environmental factors. As many circuits exceed their 40-year design life, efficient fault localization becomes crucial. Fault location involves prelocation and pinpointing. Therefore, a novel pinpointing approach for pipe-type cable systems is proposed, utilizing accelerometers mounted on a steel pipe to capture fault-induced acoustic signals and employing the time difference of arrival method to accurately pinpoint the location of the fault. The experimental investigations utilized a scaled-down HPFF pipe-type cable system setup, featuring a carbon steel pipe, high-frequency accelerometers, and both mechanical and capacitive discharge methods for generating acoustic pulses. The tests evaluated the propagation velocity, attenuation, and pinpointing accuracy with the pipe in various embedment conditions. The experimental results demonstrated accurate fault pinpointing in the centimeter range, even when the pipe was fully embedded, with the acoustic pulse velocities aligning closely with the theoretical values. These experimental investigation findings highlight the potential of this novel acoustic pinpointing technique to improve fault localization in underground systems, enhance grid reliability, and reduce outage duration. Further research is recommended to validate this approach in full-scale systems.

Moutassem, Zaki↗

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗

Recurrent convolutional neural networks for modeling nonadiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time evolution in fluid dynamics, weather predictions, and even chaotic systems due to their ability to capture temporal dependencies and sequential patterns in data. Here we present an RNN model based on convolutional neural networks for modeling the nonlinear nonadiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The Physics-Aware Recurrent Convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. Here, we apply our RNN approach to learn the space-time evolution of a one-dimensional semiclassical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

Holstein model↗