Search NASA⌕ Search

SEARCH · Search NASA

Results for “Performance Fluid”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessment of Demographic, Genetic, and Imaging Variables Associated With Brain Resilience and Cognitive Resilience to Pathological Tau in Patients With Alzheimer Disease

Importance: Better understanding is needed of the degree to which individuals tolerate Alzheimer disease (AD)-like pathological tau with respect to brain structure (brain resilience) and cognition (cognitive resilience). Objective: To examine the demographic (age, sex, and educational level), genetic (APOE-ε4 status), and neuroimaging (white matter hyperintensities and cortical thickness) factors associated with interindividual differences in brain and cognitive resilience to tau positron emission tomography (PET) load and to changes in global cognition over time. Design, setting, an participants: In this cross-sectional, longitudinal study, tau PET was performed from June 1, 2014, to November 30, 2017, and global cognition monitored for a mean [SD] interval of 2.0 [1.8] years at 3 dementia centers in South Korea, Sweden, and the United States. The study included amyloid-β-positive participants with mild cognitive impairment or AD dementia. Data analysis was performed from October 26, 2018, to December 11, 2019. Exposures: Standard dementia screening, cognitive testing, brain magnetic resonance imaging, amyloid-β PET and cerebrospinal fluid analysis, and flortaucipir (tau) labeled with fluor-18 (18F) PET. Main outcomes and measures: Separate linear regression models were performed between whole cortex [18F]flortaucipir uptake and cortical thickness, and standardized residuals were used to obtain a measure of brain resilience. The same procedure was performed for whole cortex [18F]flortaucipir uptake vs Mini-Mental State Examination (MMSE) as a measure of cognitive resilience. Bivariate and multivariable linear regression models were conducted with age, sex, educational level, APOE-ε4 status, white matter hyperintensity volumes, and cortical thickness as independent variables and brain and cognitive resilience measures as dependent variables. Linear mixed models were performed to examine whether changes in MMSE scores over time differed as a function of a combined brain and cognitive resilience variable. Results: A total of 260 participants (145 [55.8%] female; mean [SD] age, 69.2 [9.5] years; mean [SD] MMSE score, 21.9 [5.5]) were included in the study. In multivariable models, women (standardized β = -0.15, P = .02) and young patients (standardized β = -0.20, P = .006) had greater brain resilience to pathological tau. Higher educational level (standardized β = 0.23, P < .001) and global cortical thickness (standardized β = 0.23, P < .001) were associated with greater cognitive resilience to pathological tau. Linear mixed models indicated a significant interaction of brain resilience × cognitive resilience × time on MMSE (β [SE] = -0.235 [0.111], P = .03), with steepest slopes for individuals with both low brain and cognitive resilience. Conclusions and relevance: Results of this study suggest that women and young patients with AD have relative preservation of brain structure when exposed to neocortical pathological tau. Interindividual differences in resilience to pathological tau may be important to disease progression because participants with both low brain and cognitive resilience had the most rapid cognitive decline over time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Electromagnetic Energy-Assisted Thermal Conversion of Fossil-Based Hydrocarbons to Low-Cost Hydrogen

The goal of this project was to develop and optimize catalysts for methane decomposition, particularly focusing on regeneration via an electromagnetic energy-assisted mechanism, to produce hydrogen more cost-effectively compared to electrolysis routes. To achieve this goal, the project pursued several key objectives. The project began with the preparation and testing of various catalysts. A nickel-silica based catalyst was identified as the most promising material for the pyrolysis of methane into carbon and hydrogen. Kinetic parameters for methane decomposition were determined, aiding in computational modeling efforts. Structured catalysts were investigated, highlighting the need for frequent cleaning or regeneration to maintain performance, with methane conversion rates exceeding 70% in tube furnace tests. Computational fluid dynamics modeling was employed to optimize reactor designs and electrode angles, leading the project team to propose a multi-compartment thermal conversion system for larger setups. This modeling work was important in understanding reaction characteristics, carbon deposition rates, and temperature profiles under various conditions. A bench-scale reactor system was assembled to evaluate catalyst regeneration using electromagnetic energy-assisted mechanisms. Experiments demonstrated the potential for carbon removal, though further optimization is needed. The carbon produced from the methane conversion process was evaluated for potential use in lithium-ion battery electrodes. The carbon exhibited properties similar to commercially available high-purity multi-walled carbon nanotubes and nanofibers, with a carbon content greater than 95%. Coin cell batteries assembled with this carbon showed that lower replacement levels (10% to 33%) outperformed the control group, improving specific capacity density and stability. However, higher replacement levels (100%) demonstrated poorer performance, suggesting that excessive carbon substitution negatively impacts battery performance. These findings indicate the potential marketability of the produced carbon as a component in lithium-ion batteries, though further testing is necessary to confirm long-term advantages and disadvantages associated with the use of the carbon product. These results however justified further technoeconomic assessments to determine if the process can provide low-cost hydrogen. The economic feasibility and technical performance of methane decomposition for hydrogen production were assessed, focusing on three plant configurations: 100E (electrically heated), 100C (combustion heated using produced hydrogen), and PE-Hybrid (a combination of pyrolysis (indicating decomposition) and electrolysis). The Levelized Cost of Hydrogen (LCOH) for the pyrolysis configurations was found to be approximately 25% lower than that of electrolysis. The 100E configuration had the lowest LCOH at $\$$3.12/kg. Including carbon product sales significantly improved the economics, with the 100C configuration achieving a negative LCOH of -$\$$0.35/kg. The PE-Hybrid configuration was not economically advantageous compared to pure pyrolysis plants due to its complexity and additional equipment requirements. Ultimately, methane pyrolysis presents a viable method for near carbon dioxide-free hydrogen production, with significant economic advantages over electrolysis, especially when considering the sale of carbon byproducts. The 100E and 100C configurations showed the most promise, with the choice between them ultimately depending on the prices of power and natural gas. In conclusion, this technology has the potential to lower hydrogen production costs by leveraging the methane decomposition process with the sale of valuable carbon byproducts. By optimizing catalyst performance and integrating electromagnetic energy-assisted regeneration, the process can achieve higher efficiency and economic viability, making it a competitive alternative to traditional hydrogen production methods.

08 HYDROGEN↗
DOE OSTI2024

A fast, matrix-based method to perform omnidirectional pressure integration

Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.

42 ENGINEERING↗
DOE OSTI2022

PeleC: An adaptive mesh refinement solver for compressible reacting flows

Reacting flow simulations for combustion applications require extensive computing capabilities. Leveraging the AMReX library, the Pele suite of combustion simulation tools targets the largest supercomputers available and future exascale machines. We introduce PeleC, the compressible solver in the Pele suite, and detail its capabilities, including complex geometry representation, chemistry integration, and discretization. We present a comparison of development efforts using both OpenACC and AMReX’s C++ performance portability framework for execution on multiple GPU architectures. We discuss relevant details that have allowed PeleC to achieve high performance and scalability. PeleC’s performance characteristics are measured through relevant simulations on multiple supercomputers. The success of PeleC’s design for exascale is exhibited through demonstration of a 160 billion cell simulation and weak scaling onto 100% of Summit, an NVIDIA-based GPU supercomputer at Oak Ridge National Laboratory. Our results provide confidence that PeleC will enable future combustion science simulations with unprecedented fidelity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of a Hybrid Single/Two-Phase Capillary-Based Micro-Cooler Using Copper Inverse Opals Wick with Silicon 3D Manifold for High-Heat-Flux Cooling Application: Preprint

Previously, two-phase capillary-based cooling from narrow (200-1,000 micrometer) heater-bridge copper inverse opal (CIO) wicks with heat flux levels exceeding 1,400 watts per square cm with a low superheat of approximately 10 degrees C was demonstrated. Here, we demonstrate the area scaling of the proposed technology to large-area micro-cooler for the high-heat-flux cooling of microprocessors and power electronics. We developed a hybrid single/two-phase micro-cooler that relies on capillary wicking in 25-micrometer-thick CIOs with an open channel silicon 3D-manifold for liquid delivery and vapor extraction, to achieve a high heat flux of approximately 400 watts per square cm over a heated area of 1 square cm. For the range of inlet water (21 degrees C water temperature) flowrates from 5 to 60 mL per min, we achieved total thermal resistances and vapor qualities of 0.68-0.2 square cm Kelvin per watt and 0.55-0.12, respectively. The high heat flux levels are achieved with flowrates that are 10 times smaller than conventional single- or two-phase microchannel cooling technology. The corresponding two-phase thermal resistances are in the range of 0.05 to 0.02 square cm Kelvin per watt with temperature superheat of 8 to 6 degrees C. While the overall performance of the large-area (10 by 10 square millimeter) capillary-based micro-cooler degraded compared to previous demonstration of the technology for a heated area of 5 by 5 square millimeter, preliminary computational fluid dynamics (CFD) modeling indicates that an improved manifold design will be able to achieve comparable performance.

capillary flow↗
DOE OSTI2025

Numerical Investigation of a Small Low-Flow Marine Hydrokinetic (MHK) Turbine for Small Autonomous Unmanned Mobile Recharge Stations (CRADA Final Report)

This project will improve the hydrodynamic performance of the undershot waterwheel turbine with a flow concentrator designed by the Participant, using a computational fluid dynamics study. The influence of different result parameters including the number of blades, turbine profile, submerged depth, tip speed ratio, concentrator geometry and wave + current conditions will be assessed to maximize power performance and minimize structural loads.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Fluid temperature predictions of geothermal borefields using load estimations via state observers

Fluid temperature predictions of geothermal borefields usually involves temporal superposition of its characteristic g-function, using load aggregation schemes to reduce computational times. Assuming that the ground has linear properties, it can be modeled as a linear state space system where the states are the aggregated loads. However, the application and accuracy of these models is compromised when the borefield is already operating and its load history is not registered or there are gaps in the data. This paper assesses the performance of state observers to estimate the borefield load history and obtain accurate fluid predictions. Results show that both Time-Varying Kalman Filter (TVKF) and Moving Horizon Estimator (MHE) provide predictions with average and maximum errors below 0.1 $^\circ C$ and 1 $^\circ C$ respectively. MHE outperforms TVKF in terms of n-step ahead output predictions and load history profile at the expenses of about 5 times more computational time.

Figueroa, Iago C.↗
DOE OSTI2024

General Relativistic Stability and Gravitational Wave Content of Rotating Triaxial Neutron Stars

Triaxial neutron stars can be sources of continuous gravitational radiation detectable by ground-based interferometers. The amplitude of the emitted gravitational wave can be greatly affected by the state of the hydrodynamical fluid flow inside the neutron star. In this work, we examine the most triaxial models along two sequences of constant rest mass, confirming their dynamical stability. We also study the response of a triaxial figure of quasiequilibrium under a variety of perturbations that lead to different fluid flows. Starting from the general relativistic compressible analog of the Newtonian Jacobi ellipsoid, we perform simulations of Dedekind-type flows. We find that in some cases the triaxial neutron star resembles a Riemann-S-type ellipsoid with minor rotation and gravitational wave emission as it evolves towards axisymmetry. The present results highlight the importance of understanding the fluid flow in the interior of a neutron star in terms of its gravitational wave content.

Luo, Yufeng (ORCID:0000000246230683)↗
DOE OSTI2023

ExaWind: Then and Now

The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present the current status of the ExaWind software stack in the context of the modeling and simulation capabilities when the project started in 2016.

computational fluid dynamics↗
DOE OSTI2025

Thermal Hydraulic Experimental Test Article: Second Year of Testing with Secondary Sodium System (Fiscal Year 2025 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is currently installed in the Mechanisms Engineering Test Loop (METL) 28” test vessel #4. Both the primary and secondary sodium systems remain online to facilitate continued testing. This fiscal year, work was performed using a COMSOL Multiphysics magnetohydrodynamic model to characterize flow more accurately in the secondary electromagnetic flowmeters. Experimental campaigns were then performed to study the thermal hydraulic differences between sodium and water as a surrogate fluid in the THETA geometry as well as a study to better characterize and understand temperature oscillations that exist at the outlet of the core to the hot pool. A peer-reviewed article was published in the ASME Journal of Nuclear Engineering and Radiation Science detailing the THETA facility and providing an overview of a test that was performed with the primary and secondary system online [1]. Work continues to develop a database to house experimental THETA data to better facilitate collaboration with industry and laboratory partners for their use of the data for code benchmarking/validation. THETA remains fully operational and is positioned for continued testing in fiscal year 2026.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
DOE OSTI2020

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗
DOE OSTI2024

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↗
DOE OSTI2020

Regional Assessment for the CO 2 Storage Potential in the Northern Niagaran Pinnacle Reef Trend

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the regional analysis of the Northern Niagaran Pinnacle Reef Trend (NNPRT) in Northern Michigan as part of the MRCSP programs Regional Assessment goals. The goal of this task was to perform an initial assessment of the geologic storage capacity and injectivity of the Northern Fairway of Michigan’s Niagaran Reef Trend. This was achieved through three main objectives: 1) development of a reef atlas, 2) understanding the regional trends in geology, and 3) estimating CO 2 and CO 2 -EOR resources. Data collected across multiple subtasks were integrated into a mappable database that could be accessed from many of the leading software. The reef atlas was used to map the reefs and associated characteristics across the entire NNPRT. Several methodologies were explored to estimate the CO 2 and CO 2 -EOR resources across the trend which included three scenarios: 1) storage only using fluid substitution and volumetric estimations, 2) CO 2 -EOR which applied proximity analysis concepts and measured performance metrics from the Core Energy reefs to predict performance at all oil reefs, and 3) enhanced storage scenario which combined CO 2 -EOR with fluid substitution to represent maximized storage after completion of CO 2 -EOR.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Intermediate Temperature Proton‐Conducting Solid Oxide Electrolysis Cells with Improved Performance and Durability

Technical understanding: A high-performance and stable intermediate temperature proton-conducting solid oxide electrolysis cell (PCEC) will add significantly to the hydrogen roadmap of EERE. However, the low electrocatalytic activity of the air electrode and the instability of the electrolyte are technical barriers. In this research, we utilized computational fluid dynamic modeling and materials/microstructure engineering to investigate the complex reaction chain and to improve the performance and reliability of the protonic ceramic electrochemical cells. An active and stable Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ (BZCYYb) scaffold conformal coating electrode was rationally designed where the robust Pr 1.8 Ba 0.2 NiO 4+δ functions both as a catalytic active and a protective layer. The PCECs with porosity-optimized scaffold conformal coating electrodes demonstrate 1.46 Acm -2 under 1.3 V at 600°C, which is among the highest in recent years. The final optimization by computational fluid dynamic modeling and conformal coating technique was applied to evolve the PCEC air electrode from Pr 2 NiO 4+δ to Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ scaffold conformal coating. The electrolysis current was boosted from 0.445A/cm 2 to 1.46A/cm 2 at 600°C and 1.3V. A superior low degradation rate of 1mV/kh was achieved in a 5000-hour of operation. Technical effectiveness and economic feasibility: The effectiveness of this TCO conformal coating on electrolyte backbone as air electrode has been proven by long-term demonstration. Based on the highly conductive BZCYYb electrolyte and the TCO steam electrode PNO, this project rationally took advantage of the existent materials to structure a PCEC device with high TRL and low risk. The solution infiltration method used to structure the conformal coating is a standard method used by the solid oxide cell community. Very minimal adoption is needed in the current SOC fabrication process. It is a cost-effective, scalable method for manufacturing, and brings up remarkable impact on the current technical status of PCEC. Benefit to the public: The findings in the research showed the industry that the severe challenge in PCEC can be overcome in a practical way. A various energy-hydrogen ecosystem is not beyond reach. The success of this technology provides sustainable and affordable H 2 to society. It will stimulate investment from the governmental and private sectors in response to the device manufacture and deployment, and the related infrastructure construction. It will generate high-pay domestic jobs and environmental benefits to the general public in the United States.

08 HYDROGEN↗
DOE OSTI2024

Exploring the Feasibility of INCONEL® ALLOY 740H® for Power Plant Headers: Integrating Machine Learning with Computational Fluid Dynamics (CFD)

This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

A space-time parallel algorithm with adaptive mesh refinement for computational fluid dynamics

This work describes a space-time parallel algorithm with space-time adaptive mesh refinement (AMR). AMR with subcycling is added to multigrid reduction-in-time (MGRIT) in order to provide solution efficient adaptive grids with a reduction in work performed on coarser grids. This algorithm is achieved by integrating two software libraries: XBraid (Parallel time integration with multigrid. https://computation.llnl.gov/projects/parallel-timeintegration-multigrid) and Chombo (Chombo software package for AMR applications—design document, 2014). The former is a parallel time integration library using multigrid and the latter is a massively parallel structured AMR library. Employing this adaptive space-time parallel algorithm is Chord (Comput Fluids 123:202–217, 2015), a computational fluid dynamics (CFD) application code for solving compressible fluid dynamics problems. For the same solution accuracy, speedups are demonstrated from the use of space-time parallelization over the time-sequential integration on Couette flow and Stokes’ second problem. On a transient Couette flow case, at least a 1.5× speedup is achieved, and with a time periodic problem, a speedup of up to 13.7× over the time-sequential case is obtained. In both cases, the speedup is achieved by adding processors and exploring additional parallelization in time. The numerical experiments show the algorithm is promising for CFD applications that can take advantage of the time parallelism. Future work will focus on improving the parallel performance and providing more tests with complex fluid dynamics to demonstrate the full potential of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗