Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow regime”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

EXPERIMENTAL AND MODELLING ANALYSIS OF A LARGE-SCALE TWO-PHASE LOOP THERMOSYPHON

Liquid pumping requires large quantities of electrical energy, including about 7% of the energy of building heating, ventilation, and air conditioning (HVAC) systems. To replace pumped condenser-cooling tower water loops with a passive alternative system, we implemented a commercial-scale two-phase loop thermosyphon (TPLT). The unit consists of a 13 m riser integrated with a commercially available cooling tower and circulation water heater that simulates heat loads up to 25 kW. In addition to providing passive cooling capabilities, the cooling tower unit is also maintenance free, reliable, and can operate in both dry and wet modes. This study characterizes the performance (temperature difference between the evaporator and condenser) and the two-phase flow behavior of the loop under various refrigerant charges. Sight glasses installed throughout the loop are used to identify the operating flow regimes in the riser and downcomer. Over the range of operating conditions, we identified that there is an optimal refrigerant charge range for a specific heat load at which near-isothermal operation can be achieved. We further developed a model to predict the operating flow rate and gravitational height in the downcomer and compared it with the experimental data. The results show that the model agrees with the experimental data, in particular the threshold gravitational head height at which there will be subcooled liquid at the condenser exit, indicating that it can be used as a predictive tool for identifying the optimal loop charge for a given heat load.

Loop thermosyphon, passive two phase↗
DOE OSTI2023

Flow friction in a wide rectangular duct with a beveled inlet: Informal Notes of Fundamental Fluid Physics Studies

Fundamental streamwise pressure distributions were measured at Idaho State University to address a primary objective of obtaining insight into the apparent evolution of air flow regimes at moderate Reynolds numbers in thin, wide rectangular ducts with converging inlets. Ducts of 3.030 and 4.077 mm spacings s (called three and four mm), having non-dimensional lengths of (L/s) ≈ 500.8 and 372.2 and widths (W/s) ≈ 114.6 and 83.2, respectively, were compared. Flow rates ranged from Reynolds numbers based on hydraulic diameter (Re Dh ) of about 980 to 10,000 and corresponding non-dimensional downstream pressure gradient parameters (-K p ) covered the range 0.038 to 0.0061, respectively.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Overcoming persistent challenges in putting environmental flow policy into practice: a systematic review and bibliometric analysis

The implementation of environmental flows (e-flows) aims to reduce the negative impacts of hydrological alteration on freshwater ecosystems. Despite the growing attention to the importance of e-flows since the 1970s, actual implementation has lagged. Therefore, we explore the limitations in e-flows implementation, their systemic reasons, and solutions. We conducted a systematic review and a bibliometric analysis to identify peer-reviewed articles published on the topic of e-flows implementation research in the last two decades, resulting in 68 research and review papers. Co-occurrence of terms, and geographic and temporal trends were analyzed to identify the gaps in environmental water management and propose recommendations to address limitations on e-flows implementation. We identify the underlying causes and potential solutions to such challenges in environmental water management. The limitations to e-flow implementation identified were categorized into 21 classes. The most recognized limitation was the competing priorities of human uses of water (n = 29). Many secondary limitations, generally co-occurring in co-causation, were identified as limiting factors, especially for implementing more nuanced and sophisticated e-flows. The lack of adequate hydrological data (n = 24) and ecological data (n = 28) were among the most mentioned, and ultimately lead to difficulties in starting or continuing monitoring/adaptive management (n = 28) efforts. The lack of resource/capacity (n = 21), experimentation (n = 19), regulatory enforcement (n = 17), and differing authorities involved (n = 18) were also recurrent problems, driven by the deficiencies in the relative importance given to e-flows when facing other human priorities. In order to provide a clearer path for successful e-flow implementation, system mapping can be used as a starting point and general-purpose resource for understanding the sociohydrological problems, interactions, and inherited complexity of river systems. Secondly, we recommend a system analysis approach to address competing demands, especially with the use of coupled water-energy modeling tools to support decision-making when hydropower generation is involved. Such approaches can better assess the complex interactions among the hydrologic, ecological, socioeconomic, and engineering dimensions of water resource systems and their effective management. Lastly, given the complexities in environmental water allocation, implementation requires both scientific rigor and proven utility. Consequently, and where possible, we recommend a move from simplistic flow allocations to a more holistic approach informed by hydroecological principles. To ease conflicts between competing water demands, water managers can realize more 'pop per drop' by supporting key components of a flow regime that include functional attributes and processes that enhance biogeochemical cycling, structural habitat formation, and ecosystem maintenance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advanced Multi-Dimensional Capacitance Sensors Based Subsea Multiphase Mass Flow Meter to Measure and Monitor Offshore Enhanced Oil

This innovation is based on advanced multi-dimensional extensions of ECVT sensors that involve ECVT, Displacement Current Phase Tomography (DCPT), Maxwell-Wagner-Sillars polarization effect (MWS), and velocimetry which exploit the variation of electric properties between the oil, water, and gas coming out of the well. Capacitance sensors are embedded on the inside of a pipe spool and placed in-line with the subsea oil line piping on the extraction end of the well. The difference in dielectric, dielectric loss, surface polarization, and velocity of each phase are used in the multi-dimensional algorithm to measure the volume fraction, distribution, velocity, mass flow rate, and flow regime of the mixture moving through the pipe.

02 PETROLEUM↗
DOE OSTI2021

Observations of Stably Stratified Flow through a Microscale Gap

This paper reports the findings of a comprehensive field investigation on flow through a mountain gap subject to a range of stably stratified environmental conditions. This study was embedded within the Perdigão field campaign, which was conducted in a region of parallel double-ridge topography with ridge-normal wind climatology. One of the ridges has a well-defined gap (col) at the top, and an array of in situ and remote sensors, including a novel triple Doppler lidar system, was deployed around it. The experimental design was mostly guided by previous numerical and theoretical studies conducted with an idealized configuration where a flow (with characteristic velocity U 0 and buoyancy frequency N ) approaches normal to a mountain of height h with a gap at its crest, for which the governing parameters are the dimensionless mountain height G= (Nh)/U 0 and various gap aspect ratios. Modified forms of G were proposed to account for real-world atmospheric variability, and the results are presented in terms of a gap-averaged value G c . The nature of gap flow was highly dependent on G c , wherein a nearly neutral flow regime (G c <1), a transitional mountain wave regime (G c ~O(1)), and a gap-jetting regime ( G c >O(1)) were identified. The measurements were in broad agreement with previous numerical and theoretical studies on a single ridge with a gap or double ridge topography, although details vary. This is the first ever detailed field study reported on microscale (O(100 m)) gap flows, and provides useful data and insights for future theoretical and numerical studies.

Microscale, Gap, environmental fluid dynamics, win↗
DOE OSTI2021

Liquid–Liquid Microfluidic Flows for Ultrafast 5-Hydroxymethyl Furfural Extraction

We study the effects of liquid–liquid flow patterns on the extraction of 5-hydroxymethyl furfural from the aqueous phase into an organic phase (methyl isobutyl ketone, ethyl acetate, and 2-pentanol) at ultrashort residence times of 0.2–20 s in a coiled capillary microchannel. Slug, droplet, slug-droplet, dispersed, parallel, annular, and irregular two-phase flow patterns are generated. We evaluate the extraction efficiency and the volumetric mass transfer coefficient via experiments and computational fluid dynamics (CFD). The CFD predictions are in quantitative agreement with the experimental data. The convective and diffusive contributions to the mass transfer of different flow patterns are analyzed. For all solvent pairs, the extraction efficiency varies irregularly with increasing flow rate due to changes in the mass transfer mechanism caused by a transition between flow regimes. Overall, we demonstrate >90% extraction efficiencies with high volumetric mass transfer coefficients ranging from 0.006 to 2.17 s–1. Furthermore, a merit index is defined to assess the trade-off between the fast mass transfer and the pressure drop in biphasic microchannels. In general, the mass transfer rate increases significantly at high flow rates at the expense of a high-pressure drop. Furthermore, the ethyl acetate/water irregular flow pattern is optimal due to a large throughput, extraction efficiency, and mass transfer coefficient, with a minimal pressure drop.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2022

Scaling CFB risers: Maintaining microstructure dynamics

Over the years there have been numerous approaches to develop a scaleup methodology for circulating fluidized beds. Many of these past approaches have relied on macroscopic definitions of similitude to define dimensionless groups such as Reynolds number, Froude number, density ratios, and dimensional ratios involving different length scales. These approaches have demonstrated various levels of success, but often fail when used outside the bounds of experimental datasets. The issue with these methods is that they rely on macroscopic similitude with the hopes that those conditions would provide microscopic similitude as it does in single phase systems where the Buckingham Pi techniques have given reasonable success. Two of the main reasons for failure are 1) changing fluidization regimes and 2) changing material properties that result in the material falling into a different Geldart group. Both are certain to change the microstructure (local instantaneous volume solids fraction) and changes to the microstructure affect the interphase transport properties of the reactor system, namely the mass transfer and heat transfer. Here, this work presents a new approach towards scaling based upon maintaining similitude at the microscale. It reviews a novel dimensionless flow regime map developed from analysis of the deterministic chaos parameters and higher order moments to ensure microstructure similitude. The work then examines four cases from the literature to assess the scaling approach used: (1) same particle-different scale risers for Geldart Group B particles, (2) Glicksman scaling of atmospheric circulating fluidized bed combustor, (3) same particle-different scale risers for Geldart Group A particles and (4) Glicksman scaling for pressurized gasification. General findings are that none of the scaling approaches produced directly applicable microstructure similar behavior based upon the chaotic and higher moment analysis. The positive take away is that the approach can define scale conditions that will provide similar microstructure that can be us to predict transport behavior.

42 ENGINEERING↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
DOE OSTI2020

Development and performance of a 2.9 Tesla dipole magnet using high-temperature superconducting CORC® wires

Although the high-temperature superconducting (HTS) REBa 2 Cu 3 O x (REBCO, RE-rare earth elements) material has a strong potential to enable dipole magnetic fields above 20 T in future circular particle colliders, the magnet and conductor technology needs to be developed. As part of an ongoing development to address this need, here we report on our CORC® canted cosθ magnet called C2 with a target dipole field of 3 T in a 65 mm aperture. The magnet was wound with 70 m of 3.8 mm diameter CORC® wire on machined metal mandrels. The wire had 30 commercial REBCO tapes from SuperPower Inc. each 2 mm wide with a 30 µm thick substrate. The magnet generated a peak dipole field of 2.91 T at 6.290 kA, 4.2 K. The magnet could be consistently driven into the flux-flow regime with reproducible voltage rise at an engineering current density between 400-550 A mm -2 , allowing reliable quench detection and magnet protection. The C2 magnet represents another successful step towards the development of high-field accelerator magnet and CORC® conductor technologies. The test results highlighted two development needs: continue improving the performance and flexibility of CORC® wires and develop the capability to identify locations of first onset of flux-flow voltage.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗
DOE OSTI2022

Flow Loop Study of ECT-Based Volume Fraction Monitoring in Oil–Water Two-Phase Flows

Measurement of phase volume fractions in water-containing multiphase flows is necessary for the optimization of a host of industrial flow processes. Many water-containing multiphase flows can be classified as either water-dispersed or water-continuous mixtures. A recently developed approach based on Hanai’s mixture formula and utilizing electrical capacitance tomography (ECT) sensors have shown good potential for obtaining water volume fraction estimates in two-phase water-containing flows with different water salinity levels. However, the proposed approach was investigated via controlled experiments restricted to static configurations while, in practice, multiphase flows can be dynamic and unpredictable. In this work, we perform a flow loop study of the proposed ECT-based method for volume fraction estimation in oil–water two-phase flows. Here, we evaluate the performance of the proposed method in both water-dispersed and water-continuous flow regimes by employing different types of capacitive sensors in cylindrical arrangements and parallel-plate rectangular arrangements.

02 PETROLEUM↗
DOE OSTI2023

Toward local core outlet temperature monitoring in gas-cooled nuclear reactors using distributed fiber-optic temperature sensors

Gas-cooled nuclear reactors operate at temperatures up to 950 °C with significant spatial variations in power generation and coolant flow through hundreds of parallel flow channels, resulting in complex mixing at the core outlet. The high temperatures and complex mixing can damage downstream components, challenge reactor calorimetry for power determination, and result in significant conservatism in calculated peak fuel temperatures, which ultimately limits the total power output. Directly measuring each gas stream individually is unrealistic using single point thermocouples, requiring more robust sensing techniques. Fiber-optic sensors are resilient to high temperatures (up to 1,000 °C) and radiation damage. Moreover, distributed measurements can be made along the length of one fiber, making them potential candidates for monitoring local core outlet temperatures to improve core calorimetry and identify hot channels. Here, to the authors’ knowledge, this manuscript is the first to report spatially distributed fiber optic temperature measurements to quantify mixing at the outlet of a relevant orifice plate under prototypic temperature and flow regimes to assess the feasibility of using fiber-optic sensors for distributed measurements of coolant temperature in gas-cooled reactors. A single optical fiber captured dynamic changes in local temperatures, whereas the mixed outlet temperature exhibited a muted response and could not identify which channels were responsible for the change in the mixed outlet temperature. The discussion focuses on potential challenges for deploying distributed optical fibers in gas-cooled reactors, including the effects of vibrations, radiation-induced signal attenuation and drift, and routing of the fiber while minimizing flow obstructions.

42 ENGINEERING↗
DOE OSTI2020

Projecting spatiotemporally explicit effects of climate change on stream temperature: A model comparison and implications for coldwater fishes

Conservation planners and resource managers seek information about how the availability and locations of cold-water habitats will change in the future and how these predictions vary among models. In this work, we used a physical process-based model to demonstrate the implications of climate change for streamflow and water temperature in two watersheds with distinctive flow regimes: the Snoqualmie watershed (WA) and Siletz watershed (OR), USA. Our model incorporated a downscaled ensemble of global climate model outputs and was calibrated with in situ and remotely sensed water temperatures. Furthermore, we compared predictions from our processed-based model to those from a publicly available and widely used statistical model. The process-based model projected greater changes in summer maximum water temperatures for the mixed-rain-snow Snoqualmie watershed than for the rain-dominated Siletz watershed as a result of the near-complete loss of winter snowpack and significant reduction in summer flow in the Snoqualmie watershed expected by the 2080s. Both models projected generally similar future spatial patterns of maximum water temperature in the two rivers, with cool reaches distributed farther upstream and fewer in number. However, the process-based model projected higher spatial heterogeneity in water temperature due to our spatially explicit simulation of streamflow and because we calibrated the model with spatially continuous remotely sensed water temperature data. We used stream temperature projections to assess the vulnerability of Pacific salmon and trout to changes in the spatial distribution of cold-water habitats during August by the 2080s. Results suggest that salmonids may have fewer summertime cold-water habitats in both watersheds. Projected stream warming may further limit particular species and life stages, especially in the Snoqualmie watershed. Our comparison of models highlights the importance of considering what might be gained by using a process-based model for evaluating and prioritizing management actions that mitigate climate impacts on cold-water habitats for stream fishes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Influence of Oxygen Flow and Stoichiometry on Optical Properties and Damage Resistance of Hafnium Oxide Thin Films

Hafnium oxide (HfO 2 ) is predominantly used as a high-index material in multi-layer dielectric coatings for high-peak- and high-average-power lasers, but laser damage often initiates within the HfO 2 layers despite their wide bandgap. Oxygen deficiency during deposition can introduce vacancy-related sub-bandgap states and absorptive defects, lowering damage resistance. This study investigates how oxygen flow during HfO 2 deposition with ion beam sputtering (IBS) affects its stoichiometry, defect formation, and nanosecond laser-induced damage threshold (LIDT) and whether single-layer trends predict multilayer performance. Single layers were deposited at varying oxygen flows, characterized for optical and structural properties, and tested for the LIDT at 1064 nm and 355 nm. Increasing oxygen flow drove the layer toward near-stoichiometric HfO 2 , reduced the refractive index, and altered the density of surface pinhole-like features. The single-layer LIDT at 355 nm increased with oxygen, whereas the 1064 nm LIDT was comparatively less sensitive to oxygen flow, consistent with the wavelength-dependent roles of absorptive precursors and microstructural defects. In contrast, a HfO 2 -based high-reflector (HR) showed a higher LIDT at lower oxygen flow, indicating that the family of damage precursors changes between single layers and multilayers; in stacks, structural properties such as stress, gas entrapment and thermal dissipation may outweigh the isolated absorptive defects found in single layers. These results demonstrate that the optimal oxygen flow condition depends on both LIDT wavelength and film architecture. We identified, for single layers, a 15–35 sccm window for maximizing the 1064 nm LIDT and a high-flow optimum (45 sccm) for the 355 nm LIDT and, for 355 nm HR stacks, a distinct lower-flow regime (~10 sccm).

Optics and optical instruments↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2021

Performance evaluation of a Terry GS-2 steam impulse turbine with air-water mixtures

Terry steam turbines are widely used in various industries because of their robust design. Within the nuclear power generation industry, they are used in the Reactor Core Isolation Cooling System to remove decay heat during reactor isolation events. During the Fukushima Daiichi nuclear power station disaster in Japan in 2011, the Reactor Core Isolation Cooling System and associated Terry turbine operated for over 70 hours in Unit 2; this runtime is well beyond the expected operating duration. Theories suggest the turbine was subjected to a two-phase inlet flow, which could degrade the turbine performance. In this work, an experimental test rig was constructed to test a full-scale Terry model GS-2 steam turbine under two-phase air/water flows. Steady-state efficiency and torque performance maps of the turbine were developed over a range of turbine inlet pressures (1.38–4.83 bar or 20–70 psia), air mass fractions (0.05–1.0) and rotational speeds up to 4000 RPM. Furthermore, turbine performance followed expected trends with torque varying linearly and efficiency varying quadratically with rotational speed. In addition, high-speed images of the two-phase flow entering the turbine were also analyzed to understand how changes in inlet pressure and air mass fraction affect the flow regime and homogenization. The present tests with air–water two-phase mixtures are an important step towards providing an understanding of the full-scale Terry turbine’s behavior and performance curves under two-phase conditions. The results of this work will be combined with air/water and steam/water data gathered using a small-scale Terry ZS-1 steam turbine in order to understand the scaling relationship between large and small size Terry turbines and fluid pairs. The combined data set will enable further development of analytical models over a wide range of conditions and may be used to provide technical justification for expanded use of the Terry turbines in nuclear power plant safety systems and other systems.

42 ENGINEERING↗