Engineering PapersSearch

SEARCH · Engineering Papers

Results for “liquid phase ratio”

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 19 records

Better Constraining Supercooled Clouds Could Reduce Projected Warming Spread

The increase of climate sensitivity to rising greenhouse gas concentrations in the coupled model intercomparison project phase 6 (CMIP6) Earth system models (ESMs) compared to CMIP5 ESMs is primarily attributed to a larger extratropical cloud response to climate change, referred to as cloud feedback. The ratio of supercooled liquid cloud water relative to all cloud water, termed liquid phase ratio (LPR), which has also notably increased in many recent ESMs, is thought to be a primary driver of the extratropical cloud feedback increase. Unlike the preponderance of previous studies that compare native model LPR directly with observations, here we evaluate LPR over three ESM generations against Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) using an instrument simulator approach, which mimics instrument limitations and uses consistent cloud definitions and resolutions. We find that current coupled model intercomparison project (CMIP) ESMs collectively simulate greater LPR than previous CMIP generations and overestimate LPR compared to observations, contrary to previous findings, likely driven by past inconsistent comparisons of ESM outputs with CALIPSO observations. We further show that greater LPR in ESM present-day climate is unexpectedly correlated with a smaller extratropical cloud feedback, attributable to a decrease of cloud amount exceeding the increase of cloud optical depth. Finally, our results suggest that improving constraints on model LPR using our evaluation framework would likely reduce the spread in CMIP6 climate sensitivities, owing to its effect on extratropical cloud feedback from supercooled clouds.

greenhouse gas concentrations

Observational Constraint on A Feedback From Supercooled Clouds Reduces Projected Warming Uncertainty

The increase of carbon-dioxide-doubling-induced warming (climate sensitivity) in the latest climate models is primarily attributed to a larger extratropical cloud feedback. This is thought to be partly driven by a greater ratio of supercooled liquid-phase clouds to all clouds, termed liquid phase ratio. We use an instrument simulator approach to show that this ratio has increased in the latest climate models and is overestimated rather than underestimated as previously thought. In our analysis of multiple models, a greater ratio corresponds to stronger negative cloud feedback, in contradiction with single-model-based studies. We trace this unexpected result to a cloud feedback involving a shift from supercooled to warm clouds as climate warms, which corresponds to greater cloud amount and optical depth and weakens the extratropical cloud feedback. Better constraining this ratio in climate models – and thus this supercooled cloud feedback – impacts their climate sensitivities by up to 1 ˚C and reduces inter-model spread.

Climate and Earth system modelling

Impact of thermal energy storage properties on solar dynamic space power conversion system mass

A 16 parameter solar concentrator/heat receiver mass model is used in conjunction with Stirling and Brayton Power Conversion System (PCS) performance and mass computer codes to determine the effect of thermal energy storage (TES) material property changes on overall PCS mass as a function of steady state electrical power output. Included in the PCS mass model are component masses as a function of thermal power for: concentrator, heat receiver, heat exchangers (source unless integral with heat receiver, heat sink, regenerator), heat engine units with optional parallel redundancy, power conditioning and control (PC and C), PC and C radiator, main radiator, and structure. Critical TES properties are: melting temperature, heat of fusion, density of the liquid phase, and the ratio of solid-to-liquid density. Preliminary results indicate that even though overalll system efficiency increases with TES melting temperature up to 1400 K for concentrator surface accuracies of 1 mrad or better, reductions in the overall system mass beyond that achievable with lithium fluoride (LiF) can be accomplished only if the heat of fusion is at least 800 kJ/kg and the liquid density is comparable to that of LiF (1880 kg/cu m.

Juhasz, Albert J.

Impact of thermal energy storage properties on solar dynamic space power conversion system mass

A 16 parameter solar concentrator/heat receiver mass model is used in conjunction with Stirling and Brayton Power Conversion System (PCS) performance and mass computer codes to determine the effect of thermal energy storage (TES) material property changes on overall PCS mass as a function of steady state electrical power output. Included in the PCS mass model are component masses as a function of thermal power for: concentrator, heat receiver, heat exchangers (source unless integral with heat receiver, heat sink, regenerator), heat engine units with optional parallel redundancy, power conditioning and control (PC and C), PC and C radiator, main radiator, and structure. Critical TES properties are: melting temperature, heat of fusion, density of the liquid phase, and the ratio of solid-to-liquid density. Preliminary results indicate that even though overall system efficiency increases with TES melting temperature up to 1400 K for concentrator surface accuracies of 1 mrad or better, reductions in the overall system mass beyond that achievable with lithium fluoride (LiF) can be accomplished only if the heat of fusion is at least 800 kJ/kg and the liquid density is comparable to that of LiF (1800 kg/cu m).

Juhasz, Albert J.

Theory of Influence of a Low-Volatility, Soluble Impurity on Spherically-Symmetric Combustion of Fuel Droplets

Analyses are given for the evolution of liquid-phase mass fraction profiles and temperature during combustion of fuel droplets composed of binary miscible mixtures of low-volatility and high-volatility constituents, with small initial mass fractions of the low-volatility material, or slow introduction of low-volatility material into the liquid phase by absorption from the gas phase. The gas phase is assumed to remain quasi-steady and the liquid temperature spatially uniform. The ratio of the liquid-phase diffusion coefficient to the initial burning-rate constant is treated as a small parameter. Asymptotic analyses in this parameter are developed, and the conservation equations are integrated numerically in obtaining descriptions of the combustion history. It is shown that at the surface of the droplet a boundary layer arises in which the mass fraction of the low-volatility component increases with time. The results are used to explain qualitatively some observed conditions of flame contraction and liquid disruption in droplet combustion.

Shaw, B. D.

Probe measures gas and liquid mass flux in high mass flow ratio two-phase flows

Deceleration probe constructed of two concentric tubes with separator inlet operates successfully in flow fields where ratio of droplet flow rate to gas flow rate ranges from 1.0 to 20, and eliminates problems of local flow field disturbances and flooding. Probe is effective tool for characterization of liquid droplet/gas spray fields.

Burick, R. J.

Mechanistic Studies of Combustion and Structure Formation During Synthesis of Advanced Materials

Combustion in a variety of heterogeneous systems, leading to the synthesis of advanced materials, is characterized by high temperatures (2000-3500 K) and heating rates (up to 10(exp 6) K/s) at and ahead of the reaction front. These high temperatures generate liquids and gases which are subject to gravity-driven flow. The removal of such gravitational effects is likely to provide increased control of the reaction front, with a consequent improvement in control of the microstructure of the synthesized products. Thus, microgravity (mu-g) experiments lead to major advances in the understanding of fundamental aspects of combustion and structure formation under the extreme conditions of the combustion synthesis (CS) wave. In addition, the specific features of microgravity environment allow one to produce unique materials, which cannot be obtained under terrestrial conditions. The current research is a logic continuation of our previous work on investigations of the fundamental phenomena of combustion and structure formation that occur at the high temperatures achieved in a CS wave. Our research is being conducted in three main directions: 1) Microstructural Transformations during Combustion Synthesis of Metal-Ceramic Composites. The studies are devoted to the investigation of particle growth during CS of intermetallic-ceramic composites, synthesized from nickel, aluminum, titanium, and boron metal reactants. To determine the mechanisms of particle growth, the investigation varies the relative amount of components in the initial mixture to yield combustion wave products with different ratios of solid and liquid phases, under 1g and mu-g conditions; 2) Mechanisms of Heat Transfer during Reactions in Heterogeneous Media. Specifically, new phenomena of gasless combustion wave propagation in heterogeneous media with porosity higher than that achievable in normal gravity conditions, are being studied. Two types of mixtures are investigated: clad powders, where contact between reactants occurs within each particle, and mixtures of elemental powders, where interparticle contacts are important for the reaction; and 3) Mechanistic Studies of Phase Separation in Combustion of Thermite Systems. Studies are devoted to experiments on thermite systems (metal oxide-reducing metal) where phase separation processes occur to produce alloys with tailored compositions and properties. The separation may be either gravity-driven or due to surface forces, and systematic studies to elucidate the true mechanism are being conducted. The knowledge obtained will be used to find the most promising ways of controlling the microstructure and properties of combustion-synthesized materials. Low-gravity experiments are essential to create idealized an environment for insights into the physics and chemistry of advanced material synthesis processes.

Varma, A.

Reforming and decomposition of glucose in an aqueous phase

Exploratory experiments have been carried out to study the decomposition of glucose, a typical carbohydrate, in a high temperature-high pressure water reactor. The objective of the study was to examine the feasibility of such a process to decompose cellulosic waste materials in long-term space missions. At temperatures below the critical point of water, glucose decomposed to form liquid products and char. Little gas was noted with or without reforming catalysts present. The rate of the primary glucose reaction increased significantly with temperature. Partial identification of the liquid phase was made and the C:H:O ratios determined for both the liquid and solid products. One of the more interesting results from this study was the finding that when glucose was injected into a reactor held at the critical temperature (and pressure) of water, no solid products formed. Gas production increased, but the majority of the carbon was found in soluble furans (and furan derivatives). This significant result is now being investigated further.

Amin, S.

Sensitivity of Precipitation Processes to Microphysics and Resolution in a Cloud-Resolving Model

The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central U.S. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined. This will be mainly achieved by performing sensitivity tests using various horizontal (from 1- to 5-kilometers) and vertical resolutions (from 20- to 200-meters in the lower troposphere to 200- to 500-m in the middle and upper troposphere). Sensitivity tests using various microphysical schemes (warm rain only, and three ice with either graupel or hail) will also be performed. The PBL and diurnal variation of precipitation processes will also be evaluated. The budgets will be calculated for different regions (i.e., convective and stratiform regions).

Tao, Wei-Kao

Sensitivity of Precipitation Processes to Microphysics and Resolution in a Cloud-resolving Model

The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, South China Sea, eastern Atlantic, South America and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have a good sampling of cloud statistics.

Tao, Wei-Kuo

Sensitivity of Precipitation Processes to Microphysics and Resolution in a Cloud-Resolving Model

The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution on the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble of clouds. The model variables include horizontal and vertical velocities, potential temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulate squall lines that developed in the western Pacific, eastern Atlantic and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined.

Tao, Wei-Kuo

Sensitivity of Precipitation Processes of Microphysics and Resolution in a Cloud-Resolving Model

The Goddard Cumulus Ensemble (GCE) model is used to examine the impact of various microphysical schemes, and vertical and horizontal resolution ion the development, intensity and rainfall associated with mesoscale convective systems, idealized hurricanes and an ensemble f clouds. The model variables include horizontal and vertical velocities, potential temperatures, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water phases (vapor, liquid, and ice). The major characteristics of the GCE model are the explicit representation of warm rain and ice microphysical processes, and their complex interactions with solar and infrared radiative transfer processes, and with surface processes. For idealized hurricane, an axisymmetric version of the GCE model was developed and used successfully to simulate the tropical cyclogenesis process using both a Rankin vortex and saturated air within a specified radius as initial conditions. For mesoscale convective systems, the 3-D version of the GCE model was used to simulated squall lines that developed in the western Pacific, eastern Atlantic and central US. For the cloud ensemble, the GCE model was integrated for several days in order to have good sampling of cloud statistics. In this paper, the sensitivities of hurricane intensity to various microphysical processes and model grid resolution will be examined. This will be mainly achieved by performing sensitivity tests using various horizontal (from 1- to 5-km) and vertical resolutions (from 20- to 200-m in the lower troposphere to 200- to 500-m in the middle and upper troposphere). Sensitivity tests using various microphysical schemes (warm rain only, and three ice with either graupel or hail) will also be performed. The thermodynamic and water budget associated with various types of precipitation systems will also be evaluated. The budgets will be calculated for different regions (i.e., convective and stratiform regions).

Tao, Wei-Kuo

Modeling the Effects of Liquid-to-Gas Density Ratio on Slosh Dynamics

Mechanical models are commonly used in Guidance, Navigation and Controls (GN&C) system-level models to represent propellant slosh forces on launch vehicles and spacecraft. The slosh model parameters are typically predicted using semi-empirical analytical methods that are based on experimental data for high liquid-to-gas density ratios (such as water and air at standard sea level conditions). Model parameter calculations typically neglect any contribution to slosh dynamics from the gas phase. However, some cryogenic propellant systems operate in conditions where the density ratio can be orders of magnitude smaller. Analytical and computational modeling was used in this effort to investigate the effects of liquid-to-gas density ratio on slosh dynamics. No experimental data nor previous studies are available at this time that quantify the effect of liquid-to-gas density ratio on slosh dynamics. However, in a recent CFD study, results showed that density ratio can have a significant effect on slosh model parameters. In this study, a dual pendulum slosh model was created that distinctly represents both liquid and gas phase dynamics. Solution of the model, which invokes the Euler-Lagrange equations of motion, demonstrates that slosh frequency is a combination of liquid slosh frequency and gas phase slosh frequency. Additionally, it shows that slosh mass is reduced due to the opposing motion of the gas phase with respect to that of the liquid phase. The slosh model parameter trends were verified using Computational Fluid Dynamics (CFD) analysis results for various liquid-to-gas density ratios, tanks, and fill levels.

Christopher D Moore

Partition between supercooled liquid droplets and ice crystals in mixed-phase clouds based on airborne in situ observations

Abstract. The onset of ice nucleation in mixed-phase clouds determines the lifetime and microphysical properties of ice clouds. In this work, we develop a novel method that differentiates between various phases of mixed-phase clouds, such as clouds dominated by pure liquid or pure ice segments, compared with those having ice crystals surrounded by supercooled liquid water droplets or vice versa. Using this method, we examine the relationship between the macrophysical and microphysical properties of Southern Ocean mixed-phase clouds at −40 to 0 °C (e.g. stratiform and cumuliform clouds) based on the in situ aircraft-based observations during the US National Science Foundation Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) flight campaign. The results show that the exchange between supercooled liquid water and ice crystals from a macrophysical perspective, represented by the increasing spatial ratio of regions containing ice crystals relative to the total in-cloud region (defined as ice spatial ratio), is positively correlated with the phase exchange from a microphysical perspective, represented by the increasing ice water content (IWC), decreasing liquid water content (LWC), increasing ice mass fraction, and increasing ice particle number fraction (IPNF). The mass exchange between liquid and ice becomes more significant during phase 3 when pure ice cloud regions (ICRs) start to appear. Occurrence frequencies of cloud thermodynamic phases show a significant phase change from liquid to ice at a similar temperature (i.e. −17.5 °C) among three types of definitions of mixed-phase clouds based on ice spatial ratio, ice mass fraction, or IPNF. Aerosol indirect effects are quantified for different phases using number concentrations of aerosols greater than 100 or 500 nm (N>100 and N>500, respectively). N>500 shows stronger positive correlations with ice spatial ratios compared with N>100. This result indicates that larger aerosols potentially contain ice-nucleating particles (INPs), which facilitate the formation of ice crystals in mixed-phase clouds. The impact of N>500 is also more significant in phase 2 when ice crystals just start to appear in the mixed phase compared with phase 3 when pure ICRs have formed, possibly due to the competing aerosol indirect effects on primary and secondary ice production in phase 3. The thermodynamic and dynamic conditions are quantified for each phase. The results show stronger in-cloud turbulence and higher updraughts in phases 2 and 3 when liquid and ice coexist compared with pure liquid or ice (phases 1 and 4, respectively). The highest updraughts and turbulence are seen in phase 3 when supercooled liquid droplets are surrounded by ice crystals. These results indicate both updraughts and turbulence support the maintenance of supercooled liquid water amongst ice crystals. Overall, these results illustrate the varying effects of aerosols, thermodynamics, and dynamics through various stages of mixed-phase cloud evolution based on this new method that categorizes cloud phases.

54 ENVIRONMENTAL SCIENCES

Mixing of an Airblast-atomized Fuel Spray Injected into a Crossflow of Air

The injection of a spray of fuel droplets into a crossflow of air provides a means of rapidly mixing liquid fuel and air for combustion applications. Injecting the liquid as a spray reduces the mixing length needed to accommodate liquid breakup, while the transverse injection of the spray into the air stream takes advantage of the dynamic mixing induced by the jet-crossflow interaction. The structure of the spray, formed from a model plain-jet airblast atomizer, is investigated in order to determine and understand the factors leading to its dispersion. To attain this goal, the problem is divided into the following tasks which involve: (1) developing planar imaging techniques that visualize fuel and air distributions in the spray, (2) characterizing the airblast spray without a crossflow, and (3) characterizing the airblast spray upon injection into a crossflow. Geometric and operating conditions are varied in order to affect the atomization, penetration, and dispersion of the spray into the crossflow. The airblast spray is first characterized, using imaging techniques, as it issues into a quiescent environment. The spray breakup modes are classified in a liquid Reynolds number versus airblast Weber number regime chart. This work focuses on sprays formed by the "prompt" atomization mode, which induces a well-atomized and well-dispersed spray, and which also produces a two-lobed liquid distribution corresponding to the atomizing air passageways in the injector. The characterization of the spray jet injected into the crossflow reveals the different processes that control its dispersion. Correlations that describe the inner and outer boundaries of the spray jet are developed, using the definition of a two-phase momentum-flux ratio. Cross-sections of the liquid spray depict elliptically-shaped distributions, with the exception of the finely-atomized sprays which show kidney-shaped distributions reminiscent of those obtained in gaseous jet in crossflow systems. A droplet trajectory analysis overpredicts the liquid mass penetration, and indicates a need for a more rigorous model to account for the three-dimensional mixing field induced by the jet-crossflow interaction. Nonetheless, the general procedures and criteria that are outlined can be used to efficiently assess and compare the quality of sprays formed under different conditions.

Leong, May Y.

Scalable Enrichment of 48 Ca at the Solid/liquid Interface by Chemical and Electrochemical Methods

This award targets to develop methods to enrich 48 Ca, which is a critical isotope for synthesizing superheavy elements and testing the standard model through neutrinoless double beta decay. The team first tested chemical exchange-based separation between solids and liquids, which is based on the free energy change due to the different vibrational frequencies caused by Ca isotopes in a material. However, the separation factor (alpha), which is defined as the ratio of 40 Ca/ 48 Ca ratios in the two phases, only reach ~1.01. The team then developed liquid centrifugation-based isotope separation, where a Ca salt aqueous solution is centrifuged at a speed of ~60 kRPM, and 48 Ca is enriched at the bottom of a centrifuge tube due to its larger mass. A high α of ~1.2-1.4 is achieved for 40 Ca/ 48 Ca at 40 °C. This method is further approved to be generic for any isotope that can be dissolved in a liquid solution or form liquid chemicals near room temperature. The experimental results also align well with modeling prediction. The team further develop a model to evaluate isotope separation in countercurrent liquid centrifugation. The team found that the countercurrent configuration can also enhance isotope separation in liquids, similar with gas centrifugation, which boost separation for isotopes which are difficult to be gasified near room temperature.

07 ISOTOPE AND RADIATION SOURCES

Predicting Liquid–Liquid Phase Separation of Submicrometer Proxies for Atmospheric Secondary Aerosol

Liquid–liquid phase separation (LLPS) of atmospheric aerosols can significantly impact climate, air quality, and human health. However, their complex composition, small size, and history-dependent properties result in great uncertainty in the modeling of aerosol phase state and atmospheric processes. Herein, using cryogenic transmission electron microscopy (cryo-TEM), we examined model submicron aerosols composed of organic compounds and ammonium sulfate, and established a parameterization for the separation relative humidity (SRH) that accounts for chemical composition, particle size, and equilibration time. We evaluated different variables that describe chemical composition: O/C ratio, partition coefficient, solubility, molar mass, and polarizability. The O/C ratio fits the SRH of micrometer droplets best, and by using a scaling factor to translate the micrometer SRH parameterization to submicron aerosols, we incorporate the effects of size and equilibration time. The measured scaling factor for the submicron mean SRH (30nm – 1μm, 20 min equilibration times) is 0.80, the factor becomes 1 with equilibration time over 1 hour, and is equal to 0, meaning that SRH is absent, when the aerosol dry diameter is smaller than 30 nm. Furthermore, our parameterization will aid in universal SRH modeling, potentially leading to more accurate predictions of aerosol mass, optical properties, hygroscopicity, and heterogeneous chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W