Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

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

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗
DOE OSTI2020

Noncontact Flow Rate Using Laser Ultrasonics

Several types of advanced nuclear reactors are cooled with high-temperature liquid metal or molten salt flows. There is a critical need to measure flow velocity in flow channels for test purposes, and eventually in operational reactors. In conventional ultrasonic flow sensors, ultrasonic waves traveling in both the upstream and downstream directions are generated and detected by transducers that must contact the flow channels. A shift in the frequency or transit time between the two ultrasonic waves is measured to determine flow velocity. We describe here an initial effort to apply that sensing concept when the contact transducers are replaced by laser-based generation and detection instrumentation. This noncontact sensing avoids many practical problems associated with contact transducers when implemented on flow channels at high temperature. Laser-based flow monitoring can also be applied to hot-process piping in the geothermal energy, chemical-processing and petroleum-refining industries. Our effort has included theoretical simulation of noncontact laser-based flow monitoring, indicating capability of measuring flow velocities relevant to reactor cooling. It also included a room-temperature experimental demonstration using water as the flow liquid, and indicating capability of measuring flow velocity at a responsivity roughly consistent with simulation predictions. Plans were made for an experimental demonstration at Oregon State University using liquid metal at 110 degrees C.

02 PETROLEUM↗
DOE OSTI2016

Experimental Study for Two Phase Flow Instability in the PCHE Type Steam Generator in SMR

Conventional U-tube type heat exchangers have some limitations when applied to Small Modular Reactor (SMRs). In particular, the size of steam generators is important in the SMR, as it greatly affects the size, and thus the cost, of the reactor vessel.. Therefore, the circuit heat exchanger (PCHE) is suggested as an alternative. PCHE, a compact type heat exchanger manufactured by HEATRIC, has several advantages such as the reduced size, wide operation range until 900? and 900 bar and high power density due to the large heat transfer area over volume ratio. Recently, it is newly suggested as an application of a steam generator. MIT and KAIST have already conceptually suggested PCHE as a steam generator. Under boiling condition, flow instability is one of the critical issues-resulting in unstable mass flow rate, sudden temperature change and system control failure. According to the flow mechanism, there are various types of instabilities in the heat exchanger and each flow tendency is different. Typical examples are Ledinegg instability of the static instability and pressure-drop oscillation and density-wave oscillation of the dynamic instability. The characteristic of Ledinegg instability is that the flow under unstable condition changes into the stable flow. The characteristic of dynamic instability is that the flow fluctuates continuously. Pressure-drop oscillation has the period of tens of seconds and the system pressure also highly oscillates with the period of flow instability. In comparison, the density-wave oscillation has the period of few seconds and the system pressure is relatively stable. Previous studies on flow instability are focused on the single channel flow and helical tube type heat exchanger. There has been very limited experimental research on flow instability using PCHE. In KAIST, the test facility was constructed to study the thermal hydraulics in PCHE. The objectives of our test facility are: (1) to identify the occurrence of instability under two phase and super-heated steam flow in PCHE, and (2) to examine whether it can operate normally using throttling effect. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1945

Graphical and analytical methods for the determination of a flow of a compressible fluid around an obstacle

Chaplygin introduced the hodograph method in the theory of compressible fluid flows and developed a method for constructing stream functions of such flows. This method, which has been extensively used in investigation of compressible fluid flows, is limited in certain respects. The expression for the stream function obtained in this manner can represent only certain types of flow patterns. In general, flow patterns obtained in this way cannot represent the whole flow around an obstacle, but only a part of such a flow, and therefore several expressions are needed in order to obtain the whole flow. On the other hand, in many instances it is important to have a single expression representing the whole flow. Recently Von Karman and Tsien constructed more general types of stream functions, but only by replacing the true pressure density relation by the linear pressure-specific volume relation so that their method is essential limited to flows the maximum Mach number of which is not too large. In a companion report the author derived a new formula for stream functions based on the true pressure density relation. It is not subject to the limitations of the Chaplygin method. In the present report this formula is employed to construct two-dimensional subsonic compressible fluid flows around a body similar in shape to a given symmetric obstacle. The methods described in the report are illustrated by numerical examples.

FLOW, COMPRESSIBLE-MEASUREMENT↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

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

Pressure Drop-Flow Rate Curves for Single Phase Steam in Steam Generator U-tubes During Severe Accidents

Thermally induced steam generator tube rupture (TI-SGTR) can occur in pressurized water reactors (PWRs) with U-tube steam generators (SGs) during high pressure severe accident sequences resulting in containment bypass. Robust natural circulation patterns including in-vessel natural circulation and countercurrent flow to and from the SGs of superheated steam and other hot gases heats cooler structures in the reactor coolant system (RCS). The buoyancy driven flows return cooler and denser gases back to the core slowing the progression of some core degradation phenomena and failure of the lower head. As much as 50 percent of the decay energy and energy produced from core degradation oxidation reactions can be transported from the core to SGs resulting in the heat up and thermal creep of RCS pressure boundary components including hot leg nozzles, pressurizer surge line, and the SG tubes, many of which contain flaws. The creep failure of an individual tube is dependent on the transitory heating of the tube, a function of the mass flow rate through the tube and local temperature history of steam entering the tube. Buoyancy driven natural circulation flow rates are determined by changes in gravitational head, a combination of changes in fluid density due to heat transfer and change in elevation, and irreversible pressure drops resulting primarily from friction loss in the flow paths and other form losses. The nonuniform flow distribution in U-tubes of varying lengths and heights during single-phase liquid and two-phase natural circulation has been extensively studied in scaled test facilities and verified numerically in previous works. We hypothesize the countercurrent natural circulation of single-phase superheated steam in U-tube SGs is governed by the same characteristic pressure drop-flow rates curves of single-phase liquid and two-phase natural circulations in U-tubes connected to common plena derived in. Secondly, we propose flow distribution is a function of U-tube geometry - height, length, hydraulic diameter, etc., and the larger SGs of Combustion Engineering (CE) style 2 x 4 PWRs, which are the focus of this study, will show greater nonuniformity in tube by tube flow rates and ratio of forward flowing tubes to reverse flowing tubes. This work derives characteristic pressure drop-flow rate curves for single-phase superheated steam in CE-type SGs during a high pressure station blackout (SBO) sequence with failed auxiliary feedwater (AFW) and dry secondary side which is closely related to the TI-SGTR event. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Shear rate dependency on flowing granular biomass material

The commercialization of bioenergy has been significantly limited by various material handling issues due to the poor flowability of granular biomass materials. A good understanding of flow physics and robust constitutive models to predict flow behavior across multiple regimes are essential to address these issues. In this study, we investigated the multi-regime flow behavior of loblolly pine chips, a widely used bioenergy feedstock, through comprehensive inclined plane flow experiments and simulations. A quasi-static hypoplastic model and a cross-regime Drucker-Prager-µ(I) model were calibrated and validated against the physical experiments to investigate the quasi-static and dense flow behavior. The results show that for granular biomass, 1) plane flow (iso-thickness along the plane) exists within a smaller range of inclination while heap flow (varying thickness along the plane) exists within a broader range of inclination, as compared with conventional granular materials (e.g., glass beads); 2) the scaling law of granular biomass flowing on an inclined plane (Froude Number versus dimensionless thickness) forms a bi-linear trend with the turning point governed by the quasi-static and dense flow regimes; 3) the multi-regime DP-µ(I) model can capture the flow behavior in both regimes well at the cost of extra calibration. In conclusion, these findings advance the scientific understanding of the multi-regime flow behavior of granular biomass materials and shed light on formulating novel constitutive models to assist granular biomass handling in the bioenergy industry.

09 BIOMASS FUELS↗
DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2022

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗
DOE OSTI2018

Critical Review of Fluid Flow Physics at Micro‐ to Nano‐scale Porous Media Applications in the Energy Sector

While there is a consensus in the literature that embracing nanodevices and nanomaterials helps in improving the efficiency and performance, the reason for the better performance is mostly subscribed to the nanosized material/structure of the system without sufficiently acknowledging the role of fluid flow mechanisms in these systems. This is evident from the literature review of fluid flow modeling in various energy‐related applications, which reveals that the fundamental understanding of fluid transport at micro‐ and nanoscale is not adequately adapted in models. Incomplete or insufficient physics for the fluid flow can lead to untapped potential of these applications that can be used to increase their performance. This paper reviews the current state of research for the physics of gas and liquid flow at micro‐ and nanoscale and identified critical gaps to improve fluid flow modeling in four different applications related to the energy sector. The review for gas flow focuses on fundamentals of gas flow at rarefied conditions, the velocity slip, and temperature jump conditions. The review for liquid flow provides fundamental flow regimes of liquid flow, and liquid slip models as a function of key modeling parameters. The four porous media applications from energy sector considered in this review are (i) electrokinetic energy conversion devices, (ii) membrane‐based water desalination through reverse osmosis, (iii) shale reservoirs, and (iv) hydrogen storage, respectively. Review of fluid flow modeling literature from these applications reveals that further improvements can be made by (i) modeling slip length as a function of key parameters, (ii) coupling the dependency of wettability and slip, (iii) using a reservoir‐on‐chip approach that can enable capturing the subcontinuum effects contributing to fluid flow in shale reservoirs, and (iv) including Knudsen diffusion and slip in the governing equations of hydrogen gas storage.

Singh, Harpreet (ORCID:0000000276185764)↗
DOE OSTI2022

MIDAS-DC-AC Tool: Fully Automating the Acquisition of AC Power Flow Solution (MIDAS-DC2AC)

MIDAS DC2AC is an automated tool for achieving a converged AC power flow solution from any dispatch, e.g., determined using DC power flow model based optimal power flow. The entire process is free of human interference. It is usually encountered in practice that, even with a DC power flow solution, acquiring the solved AC power flow solution, if exists, sometimes could be a challenging task, especially during the planning stage. It is also difficult to distinguish the unsolvable cases from diverging iterations. Manual adjustments to approach the desired power flow condition has been largely relied on in the past using lots of engineering heuristics. This tool provides a systematic way to first achieve a solvable AC power flow case by modifying the power flow condition, and then try to track the AC power flow solution while gradually removing the adopted changes. If all adopted changes can be completely removed, then the original AC power flow solution is obtained. Otherwise, insights for actionable controls are derived to help operation and planning. Currently, this tool has been implemented in Python using SIEMENS PTI PSS/E as power flow solver, where only adjusting generator terminal voltage set point is considered as an available means to try to turn an unsolved power flow to a solved one. In future, more means should be considered, including the operation of tap-changing transformers, switched shunts and redispatch of active power.

Tan, Jin↗
DOE OSTI2019

High Energy Storage Capacity Low Cost Iron Flow Battery

A new flow battery was proposed that utilizes low cost materials: iron as the only active element, cheap aqueous electrolytes, and inexpensive separators. During charging, ferrous iron (Fe 2+ ) is oxidized to ferric iron (Fe 3+ ) at the positive electrode while it (Fe 2+ ) is reduced to form iron metal (Fe 0 ) at the negative electrode. Because iron is plated at the negative electrode during charging, conventional electrode structures couple the energy storage capacity and the power rating of the battery. In order to decouple the energy and power ratings and regain the economic advantages of a flow battery, a slurry electrode design was proposed for use as the negative electrode. The slurry electrode is made by flowing electrically conductive particles in an electrolyte containing the dissolved iron species. On charging, iron is plated onto the particles. The particles can then carry the iron metal out of the cell to be stored in external reservoirs. This project looked to develop a slurry for the negative electrode of the all-iron flow battery that can operate at 200 mA/cm 2 with <100 mV overpotential, be pumped with parasitic energy costs < 10%, and promote plating on the particles. A scalable and stackable cell design of 1100 cm 2 was developed. The project was initially under contract in mid-January 2013. The project was converted from a one-year proof of concept seedling to an additional two-year program, then extended with an 18-month plus-up, and finally with an additional 6-month approved cost-overrun. During the seeding project, the concept of plating on carbon particles in a slurry was demonstrated. However, this was accomplished with high-cost multi-walled carbon nanotubes as the slurry material. The following two-year program discovered several low-cost carbon blacks that could also be used in this application. These findings were supported by an analysis that demonstrated the effectiveness of the carbon that takes into account not just the effective electronic conductivity of the slurry, but also the surface area per volume of the slurry. It was demonstrated that stable, flowable slurries could be created with these carbons, even when the carbons where fully loaded with iron during charging. During the two-year program, laboratory-scale 50 cm 2 single-cells demonstrated charge-discharge cycling at current densities of 150 mA/cm 2 and energy efficiencies of >50%. Meanwhile, cost modeling continued to show promise for this energy storage approach. Following the conversion period, the overall objective of the nominally 18-month plus up period was to close the technical gap of a 25X scale-up and reduce risk to enable commercial investment required to develop a slurry iron flow battery product. A specific objective was to demonstrate 150 cycles of performance exceeding 70% energy efficiency with a 10-cell stack of 1156 cm 2 cells operating at 100 mA/cm 2 . The 10-cell short stack was designed to have a power delivery rating of approximately 1 kW, which represents a size commensurate with the requirements of our commercialization partner at that time. A complete battery support system for the full-size cell was designed, fabricated, and tested. This includes pumps (multiple types were considered, and constant volume hose pumps were found to be acceptable for performance and cost), electrolyte rebalance reactor, and sensors and control systems. The support system proved to be well-behaved and worked well. A new low-cost membrane was developed for this system based on a microporous polymer support with a layer of hydrogel polymer. The membrane had low resistivity, and significantly reduced pressure driven fluid flow across the membrane (as compared to the un-coated microporous support). This membrane approach was scaled-up in conjunction with a commercial toll coating company, and a large amount of membrane was manufactured to support large cell and stack testing. A new concept for an electrolyte rebalance reactor to react hydrogen gas from the negative electrolyte (from plating inefficiency) with excess ferric ion in the positive electrolyte was developed. This design was very simple, low cost, and effective. The full scale 1156 cm 2 single-cell was designed, modeled, fabricated and tested. Modeling showed that uniform flow-distribution should be obtained with the designed flow distributor. The design took into consideration the need to stack the cells, and components were fabricated for a ten-cell stack. Testing of the cell for charge and discharge demonstrated current densities of 50 mA/cm 2 with efficiencies of 50%. However, optimization and long cycling could not be accomplished because the slurry electrode would plug within the cell and block the flow. This would shut down the cell after several charge-discharge cycles. Various approaches were pursued to address this problem with some progress. However, there was not adequate time to find solutions to the problem although several promising approaches were identified for future work. The final deliverable for this proposed 18-month effort was to be a 1 kW, 6 kWh slurry flow battery, with complete balance of plant, ready for testing at an ARPA-E ‘Charges’ site with an estimated capital cost less than $\$ $45/kWh for a six-hour system (excluding power electronics). The deliverable milestone was not achieved for reasons mentioned above related to the slurry plugging. At the completion of this project, it was determined that the cost targets are achievable. The cost estimate for a 1 MW system with a one-hour capacity was estimated at $\$ $140/kWh (including power electronics) and for a six-hour system was estimated at $\$ $40/kWh (including power electronics). These costs include the battery, the balance of plant, the power conditioning and a concrete pad. This program resulted in three patents including the slurry-based iron flow battery, the rebalance reactor, and the membrane separator. At the time of this report, slurry-based iron flow battery patents have been issued in the US, China, Japan, Korea, and Europe. The other patent applications were under various stages of examination. A licensee was identified and the company initiated a scale-up and development program to take this technology to the marketplace, and CWRU staff have transferred their know-how and knowledge base to the company. This technology will enhance the economic and energy security of the United States by enabling the use of intermittent renewable energy technology, such as wind and solar power, as dispatchable resources. This report summarizes the results of this research and development program.

25 ENERGY STORAGE↗
NASA NTRS1975

Experimental and Theoretical Studies of Pulsating Turbulent Flow

The objective of this investigation was to study the effects of small amplitude sinusoidal pulsations on fully developed turbulent flow in a tube from both experimental and theoretical viewpoints. Theoretical models for the macroscopic behavior of pulsating turbulent tube flow were developed for the two cases of very low and very high pulsation frequencies. The models are based on assumptions of quasi-steady and frozen eddy viscosity flow behavior, respectively. The models successfully predict unsteady velocity profiles, thereby supporting the currently proposed definitions of frequency regimes in pulsating turbulent flow. Experimental measurements were made of the time-dependent pressure drop and velocity profiles over the range of frequency-to-Reynolds number ratios from 0.0095 to 0.24. The two macroscopic models developed in this study predict unsteady velocity profiles which are in moderately good agreement with the experiments in their respective frequency regimes, and a previously developed quasi-steady model is found to predict experimental velocity profiles well in both the quasisteady and the frozen eddy viscosity frequency regimes. The effect of flow pulsations on the dissipation of turbulence energy in the vicinity of the wall was measured in the lower transition frequency regime. The long-time averaged dissipation was observed to be unchanged from the steady flow dissipation, within the accuracy of the experiment. A theoretical model of the periodic viscous sublayer was also developed and applied to pulsating flow in a tube, in order to investigate the effects of flow pulsations on the rate of production of turbulence in the region of the wall. The periodic viscous sublayer model predicts sublayer growth periods in steady flow which agree with the published experimental data. When the model is applied to pulsating flow, the response of the sublayer growth period falls into three frequency regimes, the parameters of which are in approximate agreement with the frequency regimes which are defined on the basis of macroscopic flow behavior. The sublayer renewal cycle exhibits quasi-steady flow behavior when the sublayer growth period is much less than the pulsation period, transition behavior when these two periods are approximately equal, and frozen eddy viscosity behavior when the sublayer period is much longer than the pulsation period. The effect of the sublayer growth and renewal cycle on the level of turbulence was investigated by two methods. The velocity fluctuations seen by a point velocity probe located close to the wall were predicted from the model in one method and the rate of turbulence production was estimated from the frequency of sublayer renewal events in the other.

Kingston, G. C.↗
NASA NTRS1987

Relative ages of lava flows at Alba Patera, Mars

Many large lava flows on the flanks of Alba Patera are astonishing in their volume and length. As a suite, these flows suggest tremendously voluminous and sustained eruptions, and provide dimensional boundary conditions typically a factor of 100 larger than terrestrial flows. One of the most striking features associated with Alba Patera is the large, radially oriented lava flows that exhibit a variety of flow morphologies. These include sheet flows, tube fed and tube channel flows, and undifferentiated flows. Three groups of flows were studied; flows on the northwest flank, southeast flank, and the intracaldera region. The lava flows discussed probably were erupted as a group during the same major volcanic episode as suggested by the data presented. Absolute ages are poorly constrained for both the individual flows and shield, due in part to disagreement as to which absolute age curve is representative for Mars. A relative age sequence is implied but lacks precision due to the closeness of the size frequency curves.

Schneeberger, Dale M.↗
NASA NTRS1986

Flow-separation patterns on symmetric forebodies

Flow-visualization studies of ogival, parabolic, and conical forebodies were made in a comprehensive investigation of the various types of flow patterns. Schlieren, vapor-screen, oil-flow, and sublimation flow-visualization tests were conducted over an angle-of-attack range from 0 deg. to 88 deg., over a Reynolds-number range from 0.3X10(6) to 2.0X10(6) (based on base diameter), and over a Mach number range from 0.1 to 2. The principal effects of angle of attack, Reynolds number, and Mach number on the occurrence of vortices, the position of vortex shedding, the principal surface-flow-separation patterns, the magnitude of surface-flow angles, and the extent of laminar and turbulent flow for symmetric, asymmetric, and wake-like flow-separation regimes are presented. It was found that the two-dimensional cylinder analogy was helpful in a qualitative sense in analyzing both the surface-flow patterns and the external flow field. The oil-flow studies showed three types of primary separation patterns at the higher Reynolds numbers owing to the influence of boundary-layer transition. The effect of angle of attack and Reynolds number is to change the axial location of the onset and extent of the primary transitional and turbulent separation regions. Crossflow inflectional-instability vortices were observed on the windward surface at angles of attack from 5 deg. to 55 deg. Their effect is to promote early transition. At low angles of attack, near 10 deg., an unexpected laminar-separation bubble occurs over the forward half of the forebody. At high angles of attack, at which vortex asymmetry occurs, the results support the proposition that the principal cause of vortex asymmetry is the hydrodynamic instability of the inviscid flow field. On the other hand, boundary-layer asymmetries also occur, especially at transitional Reynolds numbers. The position of asymmetric vortex shedding moves forward with increasing angle of attack and with increasing Reynolds number, and moves rearward with increasing Mach number.

Keener, Earl R.↗
NASA NTRS1987

Fluid particle motion and Lagrangian velocities for pulsatile flow through a femoral artery branch model

A flow visualization study using selective dye injection and frame by frame analysis of a movie provided qualitative and quantitative data on the motion of marked fluid particles in a 60 degree artery branch model for simulation of physiological femoral artery flow. Physical flow features observed included jetting of the branch flow into the main lumen during the brief reverse flow period, flow separation along the main lumen wall during the near zero flow phase of diastole when the core flow was in the downstream direction, and inference of flow separation conditions along the wall opposite the branch later in systole at higher branch flow ratios. There were many similarities between dye particle motions in pulsatile flow and the comparative steady flow observations.

Cho, Y. I.↗