SEARCH · Engineering Papers

Results for “fluid flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Enhancing Fluid Flow Pressure and Saturation Prediction Accuracy and Reducing Uncertainty with Committee Machine – Illinois Basin Decatur Project (IBDP) as a Case Study

This is the conference paper accompanying an oral presentation at the 17th International Conference on Greenhouse Gas Control Technologies GHGT-17 held in Calgary, Canada, October 20-24, 2024. Carbon capture and storage (CCS) is a way to play a critical role in the global transition to a low-emission economy. Current progress is hampered by a number of factors, among which the lack of risk-informed design tools and decision support frameworks is seen as a major roadblock. Significant interest exists in using artificial intelligence to accelerate CCS site feasibility studies, as well as to facilitate the permit application process. Existing works commonly train a single deep learning model. This work investigates the feasibility of using a conventional ensemble learning (committee machine) technique to further improve prediction accuracy. Ensemble-based algorithms generally improve over individual base learners in terms of robustness and accuracy. Deep ensembles, however, are time-consuming to create and train. A pragmatic question is whether small-sized ensembles may lead to prediction improvement. Here we evaluated the efficacy of an ensemble learning technique using the latent spectral model (LSM), an efficient deep neural operator algorithm, as base learners. Preliminary results, obtained using the Illinois Basin-Decatur Project (IBDP) carbon sequestration data/model, show that small-sized ensembles can improve prediction over the base learners, achieving prediction accuracy of ~1.6 psi root mean square error (RMSE) on pressure (relative the average reservoir pressure of 3150 psi), and less than 1.3% for saturation.

Sun, Alexander↗

DOE OSTI2024

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗

DOE OSTI2024

Vacuum insulated glazing, valve, flexible seal, and method for making vacuum insulated glazing

A vacuum insulated glazing includes first and second spaced apart glass panes having perimeter portions defining a vacuum space between the glass panes. A vacuum valve is positioned at the perimeter portions of the glass panes and includes a valve body extending between the perimeter portions of the glass panes, and having first and second ends and a fluid conduit extending from the first end to the second end of the valve body and in fluid communication with the vacuum space. A one-way valve in the fluid conduit permits fluid flow from the vacuum space, and prevents fluid flow into the vacuum space. A suction fitting is provided for connecting the valve body to a suction device. An end seal between the perimeter portions of the glass panes is provided for hermetically sealing the vacuum space. A valve for a VIG and a method of making a VIG are also disclosed.

Bhandari, Mahabir↗

DOE OSTI2024

Taylor-Couette Flow Modeling (Rev.1)

The Taylor-Couette Flow is a form of viscous flow in which fluid flows between two cylindrical walls with a prescribed rotation. A diagram is shown on the right (Figure 1) where Ω is the rotational velocity, H is the distance between the two cylinders, R is the radius the outer cylinder, r is the radius of the inner cylinder and L is the length of the cylinders. When the rotating inner cylinder reaches a critical velocity, the flow of the liquid becomes unstable and creates vortices in between the layers. In this experiment, the Taylor-Couette cylinders were modeled so that the vortices would be created after reaching a critical velocity. After the initial experiment, the velocity of the inner cylinder was decreased and increased to look at how the results changed.

42 ENGINEERING↗

DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗

DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

DOE OSTI2024

Micro-architected flow through electrodes for energy storage

Disclosed are electrochemical reactors with electrodes that have variable porosity across the electrode. The electrodes are designed and micro-architected to have variable porosity and 3D flow. In one aspect, an electrochemical cell apparatus is disclosed. The apparatus includes an electrochemical vessel and an electrochemical fluid contained in the electrochemical vessel. The apparatus further includes a porous electrode submerged in the electrochemical fluid in the electrochemical vessel, the porous electrode having different porosities in different areas of the porous electrode. The different porosities inhibit electrochemical fluid flow and increase electrical conductivity in first areas of the porous electrode with decreased porosity compared to second areas, and enable increased electrochemical fluid flow and decrease electrical conductivity in the second areas of the porous electrode with increased porosity compared to the first areas.

Beck, Victor Alfred↗

DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2025

Undulating transducer for energy harnessing and propulsion apparatuses, methods and systems

Embodiments disclosed may serve the function of a pump, a generator, or a propulsion system. Two flexible fins are coupled via their longitudinal inner edges to fin coupling modules. An external fluid flow causes traveling fin undulations. The angle of the fins relative to each other changes as undulations travel along the fins which can be translated through the fin coupling modules into mechanical work, powered by the external fluid flow. The mechanical work of the fin coupling modules may be the cyclical extension and contraction of pistons or bellows incorporated into the fin coupling modules, which may draw ambient fluid during their extension cycles and compress fluid during their contraction cycles to create a pump. The pistons or bellows may also be sealed from ambient fluid but connected through conduits through which contained fluid is exchanged in a continuous circulation of fluid to power a pump or a turbine.

Filardo, Benjamin Pietro↗

DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗

DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

DOE OSTI2025

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Effect of CO 2 on the water slip flow at silica surfaces for nanometer slit pores of talc

In our previous study, Molecular Dynamics (MD) simulated water flow with no dissolved gas in nanometer slit pores of the hydrophobic talc surfaces revealed a slip length of about 0.5 nm, which is close to the size of “water exclusion zone” but much less than experimental slip length. Atomic Force Microscopy (AFM) imaging confirmed the presence of pancake shape nanobubbles at the hydrophobic talc (001) surface (size from tens to hundreds of nanometers). The effect of CO 2 on water slip flow in talc nanopores is reported to reveal the importance of dissolved gas for confined fluid flow. MD simulated 7 nm CO 2 nanobubbles were found to attach and spread at the talc (001) surface, consistent with the AFM observation. Simulated CO 2 saturated water flow in the slit pores of talc (001) surfaces predicted the generation of CO 2 nanobubbles at the end of the slit pore, due to the stabilization of a CO 2 film at the talc (001) surfaces and a water barrier at the edge surfaces of talc. A “critical thickness” of about 1.5 nm was found for CO 2 nanobubbles to be stabilized at the end of the talc (001) surface. With 1.5 nm CO 2 nanobubbles at the talc (001) surface of a 6 nm slit pore, a water slip length of 1.8 nm was determined. Furthermore, the nanometer slit pore simulations with dissolved CO 2 in water revealed an increased slip length. It is expected that an increased slip length value would be found with an increased CO 2 film thickness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗

DOE OSTI2024

Dataset for scientific paper "Simulated plant‑mediated oxygen input has strong impacts on fine‑scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands", a modeling study based on field observation at the tidal salt marshes of the Parker River Estuary, Massachusetts, United States

This dataset is the raw and processed data for the paper "Simulated plant ‑ mediated oxygen input has strong impacts on fine ‑ scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands". This study investigated how plant-mediated oxygen input affects subsurface biogeochemical reactions of organic carbon degradation and the resulting methane emissions of coastal wetlands by model simulation. We used the subsurface geochemical simulator PFLOTRAN for the modeling, which produced the simulated changes in porewater chemical substances and methane emissions over 10 days under different scenarios of plant-mediated oxygen input.Specifically, this dataset contains: 1) the input files for PFLOTRAN of all simulation runs conducted in this study. Those files are with an extension of ".in", containing information of the biogeochemical reaction network (stoichiometry, reaction rate, Monod constants, etc), fluid flow rate and oxygen concentration in the fluid which together simulated the plant-mediated oxygen input, the configuration of artificial reactions that simulated the methane fluxes, etc. The PFLOTRAN input files are text files, which can be opened by NotePad, but running these input files will require proper installation of PFLOTRAN (instruction: https://documentation.pflotran.org/user_guide/how_to/installation/installation.html). 2) the raw and processed model output from PFLOTRAN of all simulation runs, and 3) the python scripts used to process the raw model output, including random allocation of root cells, converting raw data into organized formats, calculating the methane fluxes based on the model output, data visualization, etc. The raw and processed model output from PFLOTRAN are in .spydata format, which can be viewed with Python. and 3) the python scripts for data processing and analysis are programming scripts, which can be opened with Python.This modeling work, in particular the model parameterization of root density and initial conditions of porewater concentrations of biogeochemical substances, was based on field measurements at the salt marsh of the Upper Parker River Estuary, Massachusetts, United States.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗

DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗