Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation and experiment”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗
DOE OSTI2022

Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock

Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.

58 GEOSCIENCES↗
DOE OSTI2020

Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures

Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ion transport in composites of binary electrolyte and single ion conductor—A chronoamperometry study

Composite electrolytes for lithium batteries typically combine materials with very different mechanical properties and ionic transport mechanisms and the degree to which these two phases affect each other is not well understood. In this work we used numerical simulations and experiments to investigate the transport in composite electrolytes consisting of polyethylene oxide (PEO) with Lithium bis-triuoromethanesulfonimide (LiTFSI) and Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 (LATP) lithium ion conducting glass-ceramic particles. In particular we are interested in how the introduction of a single ion conductor (SIC) changes the salt concentration gradients in the polymer electrolyte (PE) under applied potential. To study this, we performed numerical simulations and chronoamperometry experiments in electrolytes with different arrangements of the SIC and PE phases, i.e. layers and particulate composites. The results show that the particulate composites have the highest concentration gradients and take the longest time to reach steady state current. Furthermore, the high concentration gradient can be exacerbated by a high SIC/PE interfacial resistance. Finally, the best arrangement appears to have a layer of SIC impenetrable to anions in the polymer phase within the electrolyte membrane.

25 ENERGY STORAGE↗
DOE OSTI2024

Amantadine interactions with phase separated lipid membranes

Amantadine, a small amphilphic organic compound that consists of an adamantane backbone and an amino group, was first recognized as an antiviral in 1963 and received approval for prophylaxis against the type A influenza virus in 1976. Since then, it has also been used to treat Parkinson’s disease-related dyskinesia and is being considered as a treatment for corona viruses. Since amantadine usually targets membrane-bound proteins, its interactions with the membrane are also thought to be important. Biological membranes are now widely understood to be laterally heterogeneous and certain proteins are known to preferentially co-localize within specific lipid domains. Does amantadine, therefore, preferentially localize in certain lipid composition domains? To address this question, here, we studied amantadine’s interactions with phase separating membranes composed of cholesterol, DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine), POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine), and DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine), as well as single-phase DPhPC (1,2-diphytanoyl-sn-glycero-3-phos-phocholine) membranes. From Langmuir trough and differential scanning calorimetry (DSC) measurements, we determined, respectively, that amantadine preferentially binds to disordered lipids, such as POPC, and lowers the phase transition temperature of POPC/DSPC/cholesterol mixtures, implying that amantadine increases membrane disorder. Further, using droplet interface bilayers (DIBs), we observed that amantadine disrupts DPhPC membranes, consistent with its disordering properties. Finally, we carried out molecular dynamics (MD) simulations on POPC/DSPC/cholesterol membranes with varying amounts of amantadine. Consistent with experiment, MD simulations showed that amantadine prefers to associate with disordered POPC-rich domains, domain boundaries, and lipid glycerol backbones. Since different proteins co-localize with different lipid domains, our results have possible implications as to which classes of proteins may be better targets for amantadine.

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

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kick Signatures Through Multiphase Data

Oral presentation on NETL research delivered at the 2022 Resource Sustainability Meeting. This effort concerns advancing a concept for early kick detection through the development of synthetic (simulation and experiment) multiphase acoustic data corresponding to well characterized kick events. The simulation activities are the focus of the talk.

Carney, Janine↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Negative Muon Count for the FNAL Muon g-2 Experiment

The primary goal of the FNAL Muon g-2 Experiment, located in Batavia, IL, USA, is to measure the anomalous magnetic dipole moment, a μ = (g-2)/2, of a muon to a precision of 0.14 ppm. Initial results from the experiment confirm a discrepancy between the experimental value of the muon anomaly compared with the theoretical value first discovered by the Brookhaven Muon g-2 Experiment in 2001. This discrepancy hints at new physics. However, running the experiment with negative muons is essential in order to verify the experimental value of the muon anomaly. As neither experiments nor simulation studies pertaining to the muon anomaly have been performed with negative muons, the goal of this work is to estimate the rates of negative muons arriving at the FNAL Muon g-2 experiment so that scientists can estimate the time required to run the experiment with negative muons. In order to do this, we performed several simulations of the Muon Campus, which delivers muons to the FNAL Muon g-2 Experiment, using the G4Beamline v. 3.06 at the National Energy Research Computing Center. We performed separate simulations for different sections of the muon campus, and simulations were repeated twice, once for negative muons and once for positive muons. An initial 10 9 protons were shot on target for both the negative and positive muon simulations. We found that the rate of negative muons arriving at the experiment was 0.549 times the rate of positive muons arriving at the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

The cosmic muon and detector simulation framework of the extreme energy events (EEE) experiment

Abstract This paper describes the simulation framework of the extreme energy events (EEE) experiment. EEE is a network of cosmic muon trackers, each made of three multi-gap resistive plate chambers (MRPC), able to precisely measure the absolute muon crossing time and the muon integrated angular flux at the ground level. The response of a single MRPC and the combination of three chambers have been implemented in a GEANT4-based framework (GEMC) to study the telescope response. The detector geometry, as well as details about the surrounding materials and the location of the telescopes have been included in the simulations in order to realistically reproduce the experimental set-up of each telescope. A model based on the latest parametrization of the cosmic muon flux has been used to generate single muon events. After validating the framework by comparing simulations to selected EEE telescope data, it has been used to determine detector parameters not accessible by analysing experimental data only, such as detection efficiency, angular and spatial resolution.

Abbrescia, M.↗
DOE OSTI2025

Modal Field Reconstruction in Resonant Cavities in the Fundamental and Undermoded Frequency Regimes

Theory, simulations, and experiments are presented that demonstrate reconstruction of electromagnetic fields in a cavity from sparse probe measurements. Such techniques are often referred to as virtual sensing, allowing fields at unobserved locations to be predicted. These methods are appropriate for the fundamental and undermoded regimes, providing the ability to estimate fields (and shielding effectiveness) throughout an arbitrarily shaped cavity from a few judiciously spaced probes. A modal simulation method is implemented that allows the response of arbitrarily shaped cavities to be rapidly computed with respect to varying probe locations and slot parameters, enabling statistical analysis of probe placement on reconstruction performance. A cylindrical vessel with numerous probe holes is developed for experiments, referred to as Perforated Vessel 2 (PV2). Experiments are performed on the vessel with and without a steel box inside, where transmit power is delivered into the vessel either through probes (probe injection) or through slots using an external antenna (slot excitation). Simulations and experiments illustrate that when the number of probes is minimal (equal to the number of mode coefficients to be estimated at each frequency), probe placement is critical to avoid missed peaks and to have acceptable reconstruction error. Probe placement becomes less important as the number of probes is increased, but care is still required to avoid probe locations giving poor performance.

42 ENGINEERING↗
DOE OSTI2020

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fuel Performance Analysis of Chromium-Coated Cladding under Burst Conditions

To reduce the oxidation of zirconium-based alloy cladding at high temperatures, accident tolerant fuel systems have been proposed. Of the concepts identified, chromium-coated cladding has been shown to slow oxidation without greatly impacting the fuel system geometry or neutronic performance. To determine how coated-cladding tubes will perform under high-temperature accident conditions, pressurized-tube burst tests have been performed using the Severe Accident Test Station at Oak Ridge National Laboratory. To begin modeling these tubes to better understand how the coating will impact cladding behavior, these burst tests were simulated with the BISON fuel performance code. Cladding tube surface temperatures for the burst test were developed by fitting thermocouple data into axial and azimuthal profiles, while pressure data were compared until cladding failure. The temperatures at failure and the pressure evolution show relatively good agreement between the simulation and experiment results. This is the first step of a larger effort to simulate the cladding deformation process under high-temperature transient conditions and assess the cladding margin to failure more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Low Purge and Mercury Recovery Testing with Sludge Batch 10 Tank 40 Simulant

Researchers at the Savannah River National Laboratory (SRNL) have completed testing as requested by Savannah River Mission Completion (SRMC) to perform experiments to determine the impact of using a lower or inert purge in Sludge Batch (SB) 10 processing under the Nitric-Glycolic Acid (NGA) flowsheet. A key objective of this testing was also to determine the mercury speciation and recovery during each experiment. The testing was performed as part of the SB10 Technical Task Request (TTR) and Task Technical and Quality Assurance Plan (TTQAP). Two sets of tests were performed, and a Run Plan was approved prior to each set of experiments to document the planned testing. Three initial experiments were completed to determine whether a low air purge would be beneficial to CPC processing at higher acid stoichiometry (110%) based on the Koopman minimum acid (KMA) equation (116% Hsu). One of the tests, an inert nitrogen purge experiment, was also completed to demonstrate that excluding oxygen did not introduce any new hazards. The experiments were designed to be identical except for the change in purge gas and purge flowrate from run to run. After reviewing the results from the initial three experiments, six additional tests were proposed by SRNL to support the lower purge study and to look for processing alternatives for improving mercury recovery. These additional six tests were all completed at a very low acid stoichiometry to mimic the pH experienced during processing in DWPF (~7). DWPF is processing SB10 sludge at an acid stoichiometry of 90% based on the Hsu equation. The additional SRNL experiments were performed at an acid stoichiometry of 62.5% KMA stoichiometry (66.3% Hsu) to produce a Sludge Receipt and Adjustment Tank (SRAT) product with a pH of about 7. All experiments used simulants of both SWPF streams, although no entrained solvent was added during any of the experiments

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Coupled physics simulation of fracture in nuclear fuel pellets induced by resistive heating

While it is well known that the UO 2 fuel used in light-water reactors cracks early in its life in the reactor, it is very difficult to collect in-reactor data regarding fuel cracking initiation and propagation to validate fracture models in fuel performance codes. Most fracture data comes from the post-irradiation examination of fuel that has been subjected to long and complex power histories. Replicating the thermal gradients experienced in this fuel to study crack initiation and propagation in a laboratory environment is challenging. This is primarily because the fuel is heated volumetrically by fission, which is difficult to replicate. One way to approximate this volumetric heating out of the reactor is to use electrical resistance (or Joule heating). An electrical conductivity model has been developed and coupled with the existing thermomechanical modeling capability in the BISON nuclear fuel performance code. This will permit simulation of experiments that use resistive heating for the purpose of validating fracture models in BISON. The modeling approach developed here is applied to simulate both historical and current resistive heating experiments on ceramic nuclear fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Testing to Evaluate Processes Expected to Occur during MSR Salt Spill Accidents

Obtaining a license for a new nuclear reactor requires the identification and assessment of the potential consequences of specified accident scenarios, which are achieved using accident progression models. Those models need to be parameterized and validated using experimental data, but existing experimental data addressing processes relevant to molten salt reactor accidents are sparse. Specifically, experimental data that quantify the sensitives of important processes to the initial conditions of the spill, the ambient environment, and the containment features are needed to parameterize individual process models. Integrated experiments that simulate accident scenarios are also needed to provide data for model validation, but these experiments will require the use of proven methods to quantify the processes under evaluation. The overarching objectives of this work are to develop the methods for simulating the targeted processes, to determine the effectiveness of the methods in producing the data required for model development, to generate data that can be used to parameterize individual process models, and to provide key insights into the behavior of spilled molten salt that should be considered in models. Experimental methods were designed to quantify aspects of individual processes expected to occur during or after a molten salt spill accident that will affect the fate of spilled molten salt and the radionuclides within. These processes include 1) molten salt spreading and heat transfer, 2) molten salt flowing and freezing in tubing, 3) stainless steel corrosion kinetics in molten salt, and 4) molten salt splashing and aerosol generation. The initial tests described in this report were conducted using eutectic FLiNaK to demonstrate the test methods, the data that are generated, and the analyses of the data to derive values needed for modeling. The primary variables that were tested include initial salt temperature and the presence of volatile surrogate fission products (e.g., cesium and iodine). The developed methods are shown to be effective in quantifying the desired processes and can be applied to study more complex salt compositions of interest to molten salt reactor developers, a wide range of environmental conditions of interest to modelers, and additional variables relevant to salt spill accidents. The developed methods and insights gained from laboratory tests can also be incorporated in future large-scale integrated tests used to simulate molten salt spill accidents.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

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]↗