Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-phase 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 109 records · Page 6

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗

Theoretical analysis of a refrigerant injection-enhanced air source heat pump with thermal energy storage for extended cold-climate and defrosting operations

This study investigates a conceptual refrigerant injection-enhanced air source heat pump with integrated thermal energy storage (ITES-ASHP) that combines two-phase flash tank injection with phase change material (PCM) thermal energy storage (TES) to extend cold-climate operation and maintain indoor heating during defrost. A steady-state, mode-based thermodynamic cycle model is developed to evaluate heating and defrosting modes with and without refrigerant injection and TES integration over an outdoor temperature range of 8 to −25 °C. The framework is applied to refrigerants R-290, R-454B, R-410 A, and R-32. A parametric study is conducted for an 18-kW system rated for an outdoor temperature of 8 °C, with the compressor frequency scheduled to 150 Hz. The optimal PCM phase-change temperature is identified as 20 °C for R-290, R-454B, and R-410 A, and 15 °C for R-32. At the design condition of −25 °C outdoor temperature, the proposed concept achieves a combined COP of 2.22 for R-290, corresponding to a 60% increase relative to the theoretical single-stage baseline for combined heating and defrosting operation. Safety screening based on compressor discharge temperature (≤ 120 °C) and injection quality (≥ 0.65), together with performance evaluation, indicates that R-290 is the most suitable candidate, followed by R-410 A and R-454B. Under the assumptions and constraints used here, R-32 is found to be less suitable due to excessive discharge temperatures and a higher risk of wet compression.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Re-investigation of high-temperature phase equilibria in Fe-rich Sm–Fe–Ti alloys

Recently renewed attempts to develop high-performance rare-earth-lean permanent magnets based on the Sm(Fe,Ti) 12 compound have drawn attention to the limited knowledge about the high-temperature phase equilibria in the Sm–Fe–Ti system. Experimental investigation of equilibrated alloys with electron probe microanalysis, X-ray diffraction and thermomagnetic analysis revealed several inaccuracies in the currently accepted phase relations at 1000°C and allowed for a revision of the Fe-rich corner of the Sm–Fe–Ti phase diagram. The Sm(Fe,Ti) 12 and Sm 3 (Fe,Ti) 29 phases were found to have more extended Ti ranges of 5.1–9.7 at% and 2.8–6.9 at%, respectively. With increasing of the Ti content, the Curie temperature of the Sm(Fe,Ti) 12 remains nearly constant at 306–312°C, whereas that of the Sm 3 (Fe,Ti) 29 increases from 188°C to 207°C. The low-titanium Sm 3 (Fe,Ti) 29 phase equilibrates not only with the Sm(Fe,Ti) 12 and Sm2(Fe,Ti) 17 phases, but also with (α-Fe) solid solution. Newly demonstrated equilibrium between Sm 2 (Fe,Ti) 17 and TiFe 2 phases makes impossible the earlier reported equilibrium between the Sm 3 (Fe,Ti) 29 and Sm(Fe,Ti) 11 phases. Because of an invariant reaction at 1000 °C, the revised phase diagram also features a class II four-phase equilibrium Sm 3 (Fe,Ti) 29 + TiFe 2 + Sm(Fe,Ti) 12 + Sm 2 (Fe,Ti) 17 . Peritectic decomposition of the Sm(Fe,Ti) 11 phase, which occurs either at 1075°C or at 1087°C, was found to have among its products the Sm(Fe,Ti) 12 phase. Although no such equilibration was attempted, it must be possible to obtain above 1087°C a two-phase state composed of the Sm(Fe,Ti) 12 phase and a liquid – which is important for manufacturing of the Sm(Fe,Ti) 12 - based permanent magnets via the liquid-phase sintering.

36 MATERIALS SCIENCE↗

Assessment of binary eutectic Ni alloys for high-temperature applications via laser remelting

It has been recently demonstrated that eutectic alloys processed by additive manufacturing have excellent high-temperature mechanical properties. We suggest that nickel-base eutectic alloys may enable new combinations of structural and functional properties. To this end, we investigate the processability, microstructure, and thermal stability of five, binary near-eutectic Ni-X (X = B, Ce, La, Y, and Zr) alloys processed via surface laser-remelting. The microstructure of all alloys contain a two-phase lamellar eutectic microstructure consisting of γ-Ni and intermetallic phases; this microstructure is significantly finer (100–200 nm lamellar spacing) in the laser-remelted alloys than in the cast substrate (0.5–1.0 µm lamellar spacing). The microhardness of the laser-remelted alloys (550–770 HV) is 35–50% higher than that of the cast alloys (370–570 HV) due to this finer eutectic spacing. An anomalous eutectic microstructure appears at the meltpool boundaries, containing globular and lamellar γ-Ni phases. The alloys contain a high volume fraction (>50 vol%) of intermetallic phase which forms a continuous network, causing brittleness. Following laser-remelting trials, the alloys showed a high density of solid-state cracks, except for the Ni-Zr alloy which processed well. During thermal exposure at 700 and 900°C for up to 500 h, the eutectic microstructure coarsens. Coarsening occurs heterogeneously and initiates at the meltpool boundaries. This process occurs more slowly in the Ni-Zr and Ni-Y alloys, and more rapidly in the remaining alloys, resulting in greater microhardness retention in the Ni-Zr and Ni-Y alloys following thermal exposure at 700°C. Thus, among the five alloys, the Ni-Zr system exhibits a good combination of high-temperature mechanical properties and processability. Here, we conclude with recommendations for future work on designing additively manufactured alloys based on these eutectic Ni systems.

Additive manufacturing↗

Flow maldistribution in plate heat exchangers – Impact, analysis, and solutions

Plate heat exchangers are commonly used in various industrial applications, such as refrigeration, air conditioning, heat pumps, powerplants, and chemical industries. Plate heat exchangers are well known for their superior heat transfer performance, compactness, and low refrigerant charge. Despite offering several advantages, they suffer from flow maldistribution issues. The flow maldistribution can deteriorate both the heat transfer and pressure drop performance, ultimately resulting in a lower system efficiency where plate heat exchangers are deployed. The flow maldistribution issues become more pronounced when the heat exchanger size is relatively large and the number of plates is higher, limiting the deployment of plate heat exchangers in larger industrial systems. Consequently, analyzing and understanding the flow maldistribution behavior in plate heat exchangers and finding ways to mitigate flow maldistribution related issues become essential topics of interest. Further, this review aims to address the effect of flow maldistribution on plate heat exchanger characteristics. First, the experimental and numerical works on flow maldistribution under single-phase and two-phase conditions are detailed. Subsequently, the end-channel and end-plate effects are discussed. Then, the methods to mitigate flow maldistribution in plate heat exchangers are outlined. Finally, based on a thorough literature survey and industrial requirements, future research directions are recommended.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Anti-Pdc1p Nanobody as a Genetically Encoded Inhibitor of Ethanol Production Enables Dual Transcriptional and Post-translational Controls of Yeast Fermentations

Microbial fermentation provides a sustainable method of producing valuable chemicals. Adding dynamic control to fermentations can significantly improve titers, but most systems rely on transcriptional controls of metabolic enzymes, leaving existing intracellular enzymes unregulated. This limits the ability of transcriptional controls to switch off metabolic pathways, especially when metabolic enzymes have long half-lives. We developed a two-layer transcriptional/post-translational control system for yeast fermentations. Specifically, the system uses blue light to transcriptionally activate the major pyruvate decarboxylase PDC1, required for cell growth and concomitant ethanol production. Switching to darkness transcriptionally inactivates PDC1 and instead activates the anti-Pdc1p nanobody, NbJRI, to act as a genetically encoded inhibitor of Pdc1p accumulated during the growth phase. This dual transcriptional/post-translational control improves the production of 2,3-BDO and citramalate by up to 100 and 92% compared to using transcriptional controls alone in dynamic two-phase fermentations. Furthermore, this study establishes the NbJRI nanobody as an effective genetically encoded inhibitor of Pdc1p that can enhance the production of pyruvate-derived chemicals.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanical theory of nonequilibrium coexistence and motility-induced phase separation

Nonequilibrium phase transitions are routinely observed in both natural and synthetic systems. The ubiquity of these transitions highlights the conspicuous absence of a general theory of phase coexistence that is broadly applicable to both nonequilibrium and equilibrium systems. Here, in this study, we present a general mechanical theory for phase separation rooted in ideas explored nearly a half-century ago in the study of inhomogeneous fluids. The core idea is that the mechanical forces within the interface separating two coexisting phases uniquely determine coexistence criteria, regardless of whether a system is in equilibrium or not. We demonstrate the power and utility of this theory by applying it to active Brownian particles, predicting a quantitative phase diagram for motility-induced phase separation in both two and three dimensions. This formulation additionally allows for the prediction of novel interfacial phenomena, such as an increasing interface width while moving deeper into the two-phase region, a uniquely nonequilibrium effect confirmed by computer simulations. The self-consistent determination of bulk phase behavior and interfacial phenomena offered by this mechanical perspective provide a concrete path forward toward a general theory for nonequilibrium phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implementation of a Drift Flux Model into SAM with Development of a Verification and Validation Test Suite for Modeling of Noncondensable Gas Mixtures

The advanced thermal-hydraulic system code, System Analysis Module (SAM), was originally developed for the modeling of single-phase flow in advanced reactors. It has since been expanded to include a four-equation drift flux model for the modeling of two-phase flows containing a noncondensable gas. The model was expanded to support the modeling of molten salt reactor (MSR) designs in which the fuel is directly dissolved in the circulating coolant. These designs have shown that circulating gas bubbles can play an important role in the management of fission products and the operational behavior of the reactor. A drift flux model was implemented to more accurately capture the localized behavior of the void in the core and its impact on the mass transfer of fission products. A thorough assessment of the new model was performed by developing a verification and validation test suite. Verification problems were designed to test all major terms in the new governing equations. The new model converged to the correct solution at the expected order of accuracy for all verification cases. The validation cases included a wide range of flow and void conditions in different pipe geometries. Although higher void experiments show a slight underprediction of void by the drift flux model, experiments that aim to reproduce Molten Salt Reactor Experiment (MSRE) experimental conditions show good agreement with the model. The gas transport model was activated for a SAM model of the MSRE to demonstrate that it can be used in a more complex model. Finally, this gas transport model will be used along with an interfacial area transport equation being implemented in SAM for the prediction of mass transport behavior in MSR conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗

Three and Two Phase Rotating Field Inductive Couplers for Wireless Power Transfer with One Phase per Layer Windings

Multiphase inductive wireless charging coils have been proposed recently to improve coupler surface power density, reduce component stress and size, and provide near-constant power delivery to charge mobile electric systems. Several aspects for the fundamental characterization of multiphase coils are explored up to six phases including approximate mutual inductance with size and turn variation, induced voltage, and output power estimation. The relative component stress and size of passive components for resonant operation are compared between the multiphase variants. A combination of an experimentally validated 3D electromagnetic finite element analysis (FEA) and power electronic co-simulations are used to validate the estimated quantities approximated with a mixture of analytical equations and an artificial neural network model for mutual inductance. A novel three-phase transmitter, two-phase receiver coil pair is also proposed for electric vehicle charging to reduce the number of connections and compensation complexity on the vehicle-side with improved power output compared to a two-phase configuration.

Lewis, Donovin D. [University of Kentucky]↗

Understanding the Effect of Electrochemical Properties and Microstructure on the Microgalvanic Corrosion of Mg Alloys via Phase-Field Simulations

In this work, we apply a phase-field model to gain insights into the corrosion behavior of Mg alloys. Specifically, we study the effect of electrochemical properties of the materials and their environment and the spatial distribution of second phases in a magnesium alloy. To this end, we perform sensitivity analyses in which we separately vary the exchange current density and corrosion potential of the second phase, as well as the electrolyte conductivity. We describe the effects of these parameters on the corrosion current and the resulting corrosion-front morphology. We also confirm that the ratio between the Wagner length, calculated for the anodic phase, and the initial width of the anodic domain is a good predictor of the morphological characteristic of the corroding surface. In addition, we examine the effect of the microstructure on the corrosion behavior in both 2D and 3D systems by varying the shape, size, and distribution of the second-phase regions. Our findings provide insights into the influence of the two-phase microstructure on the corrosion behavior in a magnesium alloy.

Electrochemistry↗

Data for "Anti-Pdc1p Nanobody as a Genetically Encoded Inhibitor of Ethanol Production Enables Dual Transcriptional and Post-translational Controls of Yeast Fermentations"

Microbial fermentation provides a sustainable method of producing valuable chemicals. Adding dynamic control to fermentations can significantly improve titers, but most systems rely on transcriptional controls of metabolic enzymes, leaving existing intracellular enzymes unregulated. This limits the ability of transcriptional controls to switch off metabolic pathways, especially when metabolic enzymes have long half-lives. We developed a two-layer transcriptional/post-translational control system for yeast fermentations. Specifically, the system uses blue light to transcriptionally activate the major pyruvate decarboxylase PDC1 , required for cell growth and concomitant ethanol production. Switching to darkness transcriptionally inactivates PDC1 and instead activates the anti-Pdc1p nanobody, NbJRI, to act as a genetically encoded inhibitor of Pdc1p accumulated during the growth phase. This dual transcriptional/post-translational control improves the production of 2,3-BDO and citramalate by up to 100 and 92% compared to using transcriptional controls alone in dynamic two-phase fermentations. This study establishes the NbJRI nanobody as an effective genetically encoded inhibitor of Pdc1p that can enhance the production of pyruvate-derived chemicals.

metabolic engineering↗

Parametric and Sensitivity Analysis of a Steam Generator Model Using Python and Machine-Learning Tools

For this study, we used Python and machine-learning tools to perform a comprehensive parametric and sensitivity analysis on a steam generator (SG) model. (The Python model was based on a previously completed MATLAB framework for the Holtec SMR-160 SG.) We investigated the influence of various input parameters (e.g., heat transfer coefficient [HTC], Nusselt number, and heat exchanger effectiveness) on the system’s output. With machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN), which was developed at Idaho National Laboratory, we were then able to perform an automated analysis of the SG inputs’ effect on the HTC. The analysis results give valuable insights into the performance and optimization of SG systems. We found the inlet mass flow rate (MFR) to have the greatest impact on the HTC, followed closely by the inlet temperature, and then pressure. Shifting of the input parameters causes the location of the maximum HTC along the SG length to change incrementally. The cold leg (CL) MFR was also found to impact the HTC magnitude as well as the location of the maximum HTC. At between 0.4–0.9 of the total SG length, the input parameters experience maximum impact on the HTC, leading us to suggest that sensors be efficiently placed on the SG so as to closely and effectively monitor thermal-hydraulic properties during reactor operation. We also found that the sensitivity data calculated manually agrees with the RAVEN – based data, confirming the same range of maximum sensitivity. However, the RAVEN-based analysis showed that cold leg pressure and hot leg temperature have a greater impact on the heat transfer coefficient than the mass flow rate, implying that a manual sensitivity study taking only two samples is not accurate.

20 FOSSIL-FUELED POWER PLANTS↗

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗

SULI Report - Development of a Molten Salt Circulation Loop for in-situ Spectroscopy

This project supports the development of real time optical monitoring capabilities for molten salt reactor (MSR) environments by designing, testing, and refining a molten salt circulation loop suitable for combined laser induced breakdown spectroscopy (LIBS) and ultraviolet visible (UV Vis) absorption measurements. Online spectroscopic monitoring is increasingly important for nuclear safeguards, corrosion tracking, and material accountancy, yet MSR process fluids present substantial challenges due to their chemical complexity and hazards such as high temperatures and radiation. To address these needs, this work focuses on Phase I, the development of a room temperature aqueous circulation loop that serves as a surrogate platform for evaluating flow behavior, optical access, and component performance prior to high temperature salt operation. Initial testing identified several practical issues—including leaks, obstructions, and two-phase flow through the absorption cell—that were systematically resolved through hardware replacement, flow path redesign, and venturi pressure optimization. Relocating the flow cell upstream of the primary venturi enabled periods of stable single-phase flow, demonstrating the feasibility of integrating optical diagnostics into a circulation system. The results of Phase I provide essential design insight for Phase II, which will incorporate furnace compatible materials and LiCl KCl eutectic salt. Completion of the molten salt system will deliver a reusable testbed for evaluating multimodal spectroscopic techniques, advancing nondestructive, real time monitoring tools for future MSR and nuclear fuel cycle applications.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗