Search NASA⌕ Search

SEARCH · Search NASA

Results for “fluid flow systems”

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 55 records · Page 3

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↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗

Implications of rootless geothermal models: Missing processes, parameter compensation, and imposter convection

Numerical models of geothermal systems commonly capture only the top of a reservoir. Deeper areas of the reservoir are simplified to a boundary at the base of the model domain. Commonly, the basal boundary is given either a heat source or a source of mass and enthalpy. Here we developed and present simple numerical experiments which demonstrate that these approaches do not produce the correct model behavior in comparison to a model that captures the entire convecting domain with a heat flux only. We describe a variety of incorrect types of model behavior that arise directly from the choice of boundary condition, independent of the specific parameterization of the model. Heat sources are sensitive to the thickness of the domain and parameters take unphysical values to compensate for the reduced height. The combined mass/heat boundary can produce temperatures that show similarities to circulating geothermal systems, but with incorrect fluid flow and a strong boundary layer focused at the base of the simulated clay cap. These errors likely cause parameters to adjust their values to compensate for the incorrect physics. We highlight these issues and show an example from a developed reservoir model. Initial calibrations to natural state temperature were unsuitable for history matching. The 3D model required parameter adjustments to achieve a more realistic production model. Appropriate mitigation measures should be considered to reduce parameter compensation and improve decision-support models.

15 GEOTHERMAL ENERGY↗

Material coherence from trajectories via Burau eigenanalysis of braids

In this paper, we provide a numerical tool to study a material’s coherence from a set of 2D Lagrangian trajectories sampling a dynamical system, i.e., from the motion of passive tracers. We show that eigenvectors of the Burau representation of a topological braid derived from the trajectories have levelsets corresponding to components of the Nielsen–Thurston decomposition of the dynamical system. One can thus detect and identify clusters of space–time trajectories corresponding to coherent regions of the dynamical system by solving an eigenvalue problem. Unlike previous methods, the scalable computational complexity of our braid-based approach allows the analysis of large amounts of trajectories.

36 MATERIALS SCIENCE↗

Flow and pressure estimators in a waste heat recovery system

An apparatus includes a pump circuit structured to receive pump data indicative of an operating characteristic of a pump feeding a fluid to a waste heat recovery (WHR) system; a flow circuit structured to receive valve position data indicative of a position of a valve downstream of the pump, estimate a flow rate of the fluid exiting the pump, and estimate the flow rate of the fluid exiting the valve; and a pressure circuit structured to receive pressure data indicative of the pressure of the fluid exiting the valve, estimate a change in pressure of the fluid across the WHR system, and determine a pressure of the fluid in a hot section of the WHR system based on the pressure of the fluid exiting the valve and the change in the pressure of the fluid across the WHR system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Resilience and fault tolerance in high-performance computing for numerical weather and climate prediction

Progress in numerical weather and climate prediction accuracy greatly depends on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increased average frequency of hardware and software failures forces users to review their algorithms and systems in order to protect simulations from breakdown. This report surveys hardware, application-level and algorithm-level resilience approaches of particular relevance to time-critical numerical weather and climate prediction systems. A selection of applicable existing strategies is analysed, featuring interpolation-restart and compressed checkpointing for the numerical schemes, in-memory checkpointing, user-level failure mitigation and backup-based methods for the systems. Numerical examples showcase the performance of the techniques in addressing faults, with particular emphasis on iterative solvers for linear systems, a staple of atmospheric fluid flow solvers. The potential impact of these strategies is discussed in relation to current development of numerical weather prediction algorithms and systems towards the exascale. Trade-offs between performance, efficiency and effectiveness of resiliency strategies are analysed and some recommendations outlined for future developments.

54 ENVIRONMENTAL SCIENCES↗

Report on local data recovery approaches suitable for weather and climate prediction (Deliverable 1.3) (V.1.0)

Numerical weather and climate prediction rates as one of the scientific applications whose accuracy improvements greatly depend on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increasing average frequency of hardware and software failures forces users to review their algorithms and systems in order to protect simulations from breakdown. This report surveys approaches for fault-tolerance in numerical algorithms and system resilience in parallel simulations from the perspective of numerical weather and climate prediction systems. A selection of existing strategies is analyzed, featuring interpolation-restart and compressed checkpointing for the numerics, in-memory checkpointing, user-level failure mitigation-based and backup-based methods for the systems. Numerical examples showcase the performance of the techniques in addressing faults, with particular emphasis on iterative solvers for linear systems, a staple of atmospheric fluid flow solvers. The potential impact of these strategies is discussed in relation to current development of numerical weather prediction algorithms and systems towards the exascale. Trade-offs between performance, efficiency and effectiveness of resiliency strategies are analyzed and some recommendations outlined for future developments.

97 MATHEMATICS AND COMPUTING↗

Flow cell systems, flow cell batteries, and hydrogen production processes

Flow cell systems are provided. Example flow cell systems can include an H + /H 2 half-cell and a counterpart Fe 3+ /Fe 2+ or V 5+ /V 4+ half-cell. Flow cell systems can also include a half-cell in fluid communication with an electrolyte regeneration chamber. Embodiments of these flow cells systems can be configured to produce hydrogen through electrolysis. Flow cell battery systems are also disclosed. Example flow cell battery systems can include an H + /H 2 analyte; and a counterpart Fe 3+ /Fe 2+ or V 5+ /V 4+ catholyte. Processes for generating hydrogen are also disclosed. Example processes can include generating protons from a Fe 3+ /Fe 2+ or V 5+ /V 4+ electrolyte solution; and reacting the protons with H 2 O to form H 2 .

Wang, Wei↗

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗

Opportunities and Challenges in Passive Thermal-Fluid and Energy Systems

Abstract This article focuses on passive systems that are used in energy and thermal-fluid applications. These passive systems do not have moving parts and are reliable and cost-effective. Fluid motion in these passive devices could be driven by capillary force, gravity, osmotic pressure, and/or concentration gradient. The fundamental mechanisms and limitations of transport phenomena for passive systems are highlighted, followed by their applications in heat pipes, fuel cells, thermal energy storage, and desalination systems. The capabilities of the passive systems are limited by the balance between the driving force and transport resistance. Based on the fundamental understanding of fluid flow and phase change in passive systems, this study proposes associated transport phenomena and quantitative criteria to determine the maximum heat transfer rate, the transport distance, and minimum pore size of wick structures (when relevant) in these passive devices. This article concludes with the discussion of challenges and future opportunities of passive systems.

Engineering↗

Efficient Thermal Energy Storage with Radial Flow in Packed Beds

Thermal energy storage (TES) is a way to store excess heat in order to generate power later. These types of systems are valuable in the solar industry, where power production can still occur when sunlight is not available. This not only increases system flexibility, it drives down the cost of electricity. One option to store thermal energy is with a packed bed where the storage media resides in a cylindrical container. Flow from one end of the cylinder to the other deposits or recovers heat (axial flow). While a promising technology in terms of energy storage, it exhibits a high pressure drop that lowers overall system efficiency. This project seeks to flow the heat transfer fluid through the storage media in the radial direction. This system offers the potential to retain reasonable thermal performance while substantially lowering the pressure drop. The overall goal was to assess the performance of radial flow experimentally and numerically. Three main designs were considered in this work. The first system utilized a central pipe and four receiving pipes near the wall. The second system utilized a central pipe and an annulus with holes near the wall. The last approach considered segmenting the bed so only select zones would receive flow in the radial direction. In a high-aspect ratio system, the hole patterns in the piping are important for ensuring even flow into the bed, and a major design effort was testing holes patterns to promote even flow. The four pipe system was shown to be not feasible, as preferential flow paths occur that do not allow even and full heating of the packed bed. The annular system can lead to even flow, but the results show this is better accomplished through larger piping instead of considering variations in the hole pattern along the axial length of the pipe. Segmenting the bed can lead to similar exergetic performance when compared to axial or radial flow. However, the pressure drop in the segmented design is likely too high for practical implementation. A baseline commercial design showed exergetic efficiencies for axial, radial, and radial segments of 81.6%, 82.8%, and 80.2%, respectively. Pressure drop for the axial and radial results were 2.36 psi and 2.45 psi, respectively, with segments being nearly an order of magnitude higher. At these large scales, the aspect ratio of the system is important. From a baseline of 0.64, an aspect ratio of 0.32 for radial flow showed an exergetic efficiency of 86% but a pressure drop of 5.39 psi. All of these results provide new insights into packed bed thermal energy storage with radial flow. Competing effects must be considered when designing a radial system, and results show a radial design can show strong thermal results at the expense of system efficiencies from pressure drop. The COMSOL models used in this analysis are available to the public and other researchers. Public benefits include a potential increase in the efficiency of packed bed thermal energy storage. Higher efficiency in storage promotes green energy technologies while reducing costs.

14 SOLAR ENERGY↗

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

Using Machine Learning Tools to Predict Compressor Stall

Clean energy has become an increasingly important consideration in today’s power systems. As the push for clean energy continues, many coal-fired power plants are being decommissioned in favor of renewable power sources such as wind and solar. However, the intermittent nature of renewables means that dynamic load following traditional power systems is crucial to grid stability. With high flexibility and fast response at a wide range of operating conditions, gas turbine systems are poised to become the main load following component in the power grid. Yet, rapid changes in load can lead to fluid flow instabilities in gas turbine power systems. These instabilities often lead to compressor surge and stall, which are some of the most critical problems facing the safe and efficient operation of compressors in turbomachinery today. Although the topic of compressor surge and stall has been extensively researched, no methods for early prediction have been proven effective. This study explores the utilization of machine learning tools to predict compressor stall. The long short-term memory (LSTM) model, a form of recurrent neural network (RNN), was trained using real compressor stall datasets from a 100 kW recuperated gas turbine power system designed for hybrid configuration. Two variations of the LSTM model, classification and regression, were tested to determine optimal performance. The regression scheme was determined to be the most accurate approach, and a tool for predicting compressor stall was developed using this configuration. Overall, results show that the tool is capable of predicting stalls 5–20 ms before they occur. With a high-speed controller capable of 5 ms time-steps, mitigating action could be taken to prevent compressor stall before it occurs.

42 ENGINEERING↗

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗

Effects of flow on uranium speciation in soils impacted by acidic waste fluids

Radioactive acidic liquid waste is a common byproduct of uranium (U) and plutonium (Pu) enrichment and recycling processes whose accidental and planned release has led to a significant input of U into soils and sediments across the world, including at the U.S. DOE's Hanford site (WA, USA). Because of the particularly hazardous nature of U, it is important to predict its speciation when introduced into soils and sediments by acidic waste fluids. Of fundamental importance are the coupled effects of acid-driven mineral transformation and reactive transport on U speciation. To evaluate the effect of waste-fluid residence time and co-associated dissolved phosphate concentrations on U speciation in impacted soils and sediments, uncontaminated surface materials (from the Hanford Site) were reacted with U-containing synthetic acidic waste fluids (pH 2) amended with dissolved phosphate concentrations in both batch (no flow) and flow-through column systems for 7–365 days. By comparing dissolved U behavior and solid phase speciation as a function of flow regimen, we found that the availability of proton-promoted dissolution products (such as Si) to sequester U into uranyl silicates was dependent on waste fluid-sediment contact time as uranyl silicates were not detected in short contact time flow-through systems but were detected in no-flow, long contact time, reactors. Moreover, the dominance of uranyl phosphate as neoprecipitate U scavenger (principally in the form of meta-ankoleite) in phosphate amended systems confirmed the importance of phosphate amendments for an efficient sequestration of U in the soils and sediments. Overall, our experiments suggest that the formation of uranyl silicates in soils impacted by acidic waste fluids is likely to be limited unless reaction products are allowed to accumulate in soil pores, highlighting the importance of investigating soil U speciation in flow-through, transport-driven systems as opposed to no-flow, batch systems. Finally, this study provides insights into uranium speciation and its potential changes under acidic conditions for better prediction of risks and subsequent development of efficient remediation strategies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗