Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical equilibrium”

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

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Machine Learning Applications to Metal-Silicate Equilibria and their Insights into Core Formation

An extensive number of studies have experimentally investigated how elements distribute between metal and silicate phases, to better constrain core-mantle chemical equilibrium. Here, we present a new database compiling all (to our knowledge) experimental data on liquid metal-silicate partitioning from 118 peer-reviewed publications. We applied various machine learning techniques to gain further insights into these partitioning equilibria and their dependencies. We performed a network analysis to investigate the relationship between experiments and partition coefficients, which enables visualizing gaps in the experimental dataset and biases related to varying experimental conditions and analytical setup. In addition, semi-empirical thermodynamic models are commonly used to extrapolate these chemical reactions to the wide range of pressure, temperature and compositional conditions of planetary differentiation. These models are based on linear regressions that assume continuous relationship between partition coefficients and experimental variables. Here, we considered random forest regressions, which are algorithms based on ensembles of decision trees and does not consider continuous effects of each variable. The application of this regression significantly improves the prediction of metal-silicate partitioning for several elements including Ni, Si and Cr. We will show how this new approach improves our understanding of elemental exchange between metal and silicate and their implications for the Earth’s core formation.

siderophile element↗
NASA NTRS2016

Near-Infrared Spectroscopy of the Y0 WISEP J173835.52+273258.9 and the Y1 WISE J035000.32-565830.2: the Importance of Non-Equilibrium Chemistry

We present new near-infrared spectra, obtained at Gemini Observatory, for two Y dwarfs: WISE J035000.32−565830.2 (W0350) and WISEP J173835.52+273258.9 (W1738). A FLAMINGOS-2 R = 540 spectrum was obtained for W0350, covering 1.0 < λ micrometer < 1.7, and a cross-dispersed GNIRS R = 2800 spectrum was obtained for W1738, covering 0.993 - 1.087 micrometer, 1.191 - 1.305 micrometer, 1.589 - 1.631 micrometer, and 1.985 - 2.175 micrometer, in four orders. We also present revised Y JH photometry for W1738, using new NIRI Y and J imaging, and a re-analysis of the previously published NIRI H band images. We compare these data, together with previously published data for late-T and Y dwarfs, to cloud-free models of solar metallicity, calculated both in chemical equilibrium and with disequilibrium driven by vertical transport. We find that for the Y dwarfs the non-equilibrium models reproduce the near-infrared data better than the equilibrium models. The remaining discrepancies suggest that fine-tuning the CH4/CO and NH3/N2 balance is needed. Improved trigonometric parallaxes would improve the analysis. Despite the uncertainties and discrepancies, the models reproduce the observed near-infrared spectra well. We find that for the Y0, W1738, T(sub eff) = 425 +/- 25K and log g = 4.0 +/- 0.25, and for the Y1, W0350, T(sub eff) = 350 +/- 25K and log g = 4.0 +/- 0.25. W1738 may be metal-rich. Based on evolutionary models, these temperatures and gravities correspond to a mass range for both Y dwarfs of 3 - 9 Jupiter masses, with W0350 being a cooler, slightly older, version of W1738; the age of W0350 is 0.3 - 3 Gyr, and the age of W1738 is 0.15 - 1 Gyr.

two Y dwarfs↗
NASA NTRS1991

Oxygen production by pyrolysis of lunar regolith

Oxygen was identified as the most important product of initial lunar materials processing efforts. A source of oxygen on the Moon provides an alternative to the costly transport of propellant to the Moon or to low earth orbit. Pyrolysis, or vapor-phase reduction, involves heating a feedstock to temperatures sufficient to decompose the constituent metal oxides and release oxygen. The process relies on the vaporization of metal oxides in the form of reduced suboxides or atomic species. The reduced species must then be condensed without re-oxidizing, yielding oxygen in the gas phase. The feasibility of obtaining oxygen from common lunar minerals was demonstrated using solar furnace experiments. These results are discussed together with chemical equilibrium models which were extended to include the multicomponent oxides used in experiments. For the first time, both experiments and theoretical models dealt with the complex oxides that make up potential lunar feedstocks. Two major conclusions are drawn from this preliminary work. First, unbeneficiated regolith is a suitable feedstock for pyrolysis. Second, the process can operate at moderate temperatures, circa 2000 K, which could be supplied by direct solar or electrical energy. In addition to these advantages in choice of feedstock and energy source, the pyrolysis process requires no chemicals or reagents, making it an attractive process for lunar oxygen production.

Senior, Constance L.↗
NASA NTRS1994

The energetics of Titan's ionosphere

We have developed a comprehensive model to study the dynamics and energetics of the ionosphere of Titan. We solved the one-dimensional, time-dependent, coupled continuity and momentum equations for several ion species, together with single ion and electron energy equations, in order to calculate density, velocity, and temperature profiles. Calculations were carried out for several cases corresponding to different local times and configurations of the Titan-Saturn system. In our model the effects of horizontal magnetic fields were assumed to be negligible, except for their effect on reducing the electron and ion thermal conductivities and inhibiting vertical transport in the subram region. The ionospheric density peak was found to be at an altitude of about 1100 km, in accordance with earlier model calculations. The ionosphere is chemically controlled below an altitude of about 1500 km. Above this level, ion densities differ significantly from their chemical equilibrium values due to strong upward ion velocities. Heat is deposited in a narrow region around the ionospheric peak, resulting in temperature profiles increasing sharply and reaching nearly constant values of 800-1000 deg K for electrons and 300 deg K for ions in the topside, assuming conditions appropriate for the wake region. In the subram region magnetic correction factors make the electron heat conductivities negligible, resulting in electron temperatures increasing strongly with altitude and reaching values in the order of 5000 deg K at our upper boundary located at 2200 km. Ion chemical heating is found to play an important role in shaping the ion energy balance in Titan's ionosphere.

Roboz, A.↗
NASA NTRS1972

The effect of phosphorus on the formation of the Widmanstaetten pattern in iron meteorites.

Use of a combination of a revised Fe-Ni-P phase diagram and laboratory cooling experiments on Fe-Ni-P alloys to determine the effect of P on the formation of the Widmanstaetten pattern. From the phase diagram results, two reaction paths were found for the formation of kamacite (1) gamma (taenite) yields alpha (kamacite) + gamma (taenite) and (2) gamma yields gamma + Ph /phosphide, (FeNi)3P/ yields alpha + gamma + Ph. The reaction path gamma yields alpha + gamma is preferred at low P contents, while at higher P contents and at Ni contents greater than 7.0 wt.%, the reaction path gamma yields gamma + Ph yields alpha + gamma + Ph controls the formation of kamacite. Above 7 wt.% Ni, the effect of P on the equilibrium nucleation temperature of kamacite is quite small, less than about plus or minus 30 C with respect to the Fe-Ni binary diagram. The addition of P (greater than 0.1 wt.%) to meteorites promotes nucleation of kamacite at higher temperatures and effectively lowers the amount of undercooling necessary to nucleate kamacite homogeneously. Ni has just the opposite effect, decreasing the temperature of nucleation and increasing the amount of undercooling. It is concluded that significant amounts of undercooling, 50 to 100 C, are necessary for the nucleation of the Widmanstaetten structure in meteorites, and that chemical equilibrium is maintained in the various phases of slowly cooled iron meteorites to 650 C and probably to 600 C.

Goldstein, J. I.↗
NASA NTRS2011

High Power Flex-Propellant Arcjet Performance

A MW-class electrothermal arcjet based on a water-cooled, wall-stabilized, constricted arc discharge configuration was subjected to extensive performance testing using hydrogen and simulated ammonia propellants with the deliberate aim of advancing technology readiness level for potential space propulsion applications. The breadboard design incorporates alternating conductor/insulator wafers to form a discharge barrel enclosure with a 2.5-cm internal bore diameter and an overall length of approximately 1 meter. Swirling propellant flow is introduced into the barrel, and a DC arc discharge mode is established between a backplate tungsten cathode button and a downstream ringanode/ spin-coil assembly. The arc-heated propellant then enters a short mixing plenum and is accelerated through a converging-diverging graphite nozzle. This innovative design configuration differs substantially from conventional arcjet thrusters, in which the throat functions as constrictor and the expansion nozzle serves as the anode, and permits the attainment of an equilibrium sonic throat (EST) condition. During the test program, applied electrical input power was varied between 0.5-1 MW with hydrogen and simulated ammonia flow rates in the range of 4-12 g/s and 15-35 g/s, respectively. The ranges of investigated specific input energy therefore fell between 50-250 MJ/kg for hydrogen and 10-60 MJ/kg for ammonia. In both cases, observed arc efficiencies were between 40-60 percent as determined via a simple heat balance method based on electrical input power and coolant water calorimeter measurements. These experimental results were found to be in excellent agreement with theoretical chemical equilibrium predictions, thereby validating the EST assumption and enabling the utilization of standard TDK nozzle expansion analyses to reliably infer baseline thruster performance characteristics. Inferred specific impulse performance accounting for recombination kinetics during the expansion process implied nearly frozen flow in the nozzle and yielded performance ranges of 800-1100 sec for hydrogen and 400-600 sec for ammonia. Inferred thrust-to-power ratios were in the range of 30-10 lbf/MWe for hydrogen and 60-20 lbf/MWe for ammonia. Successful completion of this test series represents a fundamental milestone in the progression of high power arcjet technology, and it is hoped that the results may serve as a reliable touchstone for the future development of MW-class regeneratively-cooled flex-propellant plasma rockets.

Litchford, Ron J.↗
DOE OSTI2020

Calculating n-point charge correlations in evolving systems

Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Kinetic and equilibrium reactions on natural and laboratory generation of thermogenic gases from Type II marine shale

The phenomenon that laboratory pyrolysis experiments produce much wetter gases than those in natural reservoirs is a long-recognized and debated problem in the investigation of natural gases in sedimentary basins. In this study, we explore the discrepancy by pyrolyzing a type II kerogen from the Woodford Shale in Oklahoma, compared with the previous results on the produced natural gases from the Arkoma Basin generated from the same source rock (Liu et al., 2019) with the discussion of gas and isotopic compositions at bulk and position-specific (PS) levels. An improved GC-pyrolysis-GC IRMS method is applied for the determination of PS δ 13 C of propane produced in the pyrolysis of the Woodford Shale at Easy %R o from 0.76 to 3.27. Kinetic and thermodynamic considerations of the chemical and isotopic compositions of the natural and laboratory pyrolysis gases suggest that the generation of light hydrocarbons involves uni-directional cracking reactions, exchange reactions with water, and likely reversible reactions among light hydrocarbons and other H-containing volatiles. After the gas generation in the unconventional Woodford Shale reservoirs, the C 1 -C 4 gases might have approached close to chemical equilibrium of C 1 -C 3 and isotope equilibrium of C 2 -C 1 and C 3 -C 1 pairs at their peak temperatures. The capping H for the generation of C 1 -C 4 in the Woodford Shale gases appears to have experienced at least partial exchange with the water, while that in the pyrolysis gases is only originated from organic-bound compounds with large kinetic isotope effects (KIE). Our findings indicate that elevated compound-specific and PS δ 13 C values of propane in the wet-gas cracking stage are significantly influenced by the breakdown of the thermally stable compounds (e.g., remaining kerogen, residues). A first synthesis of PS δ 13 C and δ 2 H isotopic compositions of propane from this study and the literature data suggests relatively similar isotopic structures of propane precursors in kerogens. Finally, this study demonstrates that PS isotope analysis of propane can contribute to identifying various geological (e.g., maturation, wet-gas cracking, H exchange, diffusion) and biodegradation processes.

58 GEOSCIENCES↗
NASA NTRS1993

Carbonate compositions in CM and CI chondrites, and implications for aqueous alteration

Carbonate minerals in fourteen CM chondrites and two CI chondrites have been analyzed by electron microprobe to provide a better understanding of the aqueous processes that affected carbonaceous chondrite parent bodies. Calcites in CM chondrites and dolomites and magnesites in CI chondrites display the compositions expected of stable phases formed at low temperatures. Dolomites in CM chondrites, identified here for the first time in five members of the group, have small amounts of excess Ca which may reflect metastable growth. The distribution of Fe between dolomite and coexisting serpentine differs in the two chondrite groups. If the distributions reflect an approach to chemical equilibrium, then the difference implies higher alteration temperatures for the CI group than the CM group in agreement with the results of previously published oxygen isotope thermometry and mineral solubility modeling of the alteration process. Dolomite Fe contents are relatively uniform in the two chondrite groups. Dolomite Mn contents, by contrast, vary widely. The variations may reflect transport-controlled coprecipitation of Mn resulting from a heterogeneous distribution of the element in the anhydrous precursor material. If this interpretation is correct, then the altering fluids were essentially immobile white hydration reactions proceeded on the meteorite parent bodies. The near closed-system character of the alteration process, long known from bulk chemical analyses of the meteorites, is a direct consequence of the limited mobility of dissolved species.

Johnson, Craig A.↗
NASA NTRS1990

Phase compatibilities of YBa2Cu3O(9-delta) type structure in quintenary systems Y-Ba-Cu-O-X (impurity)

Electrical transport properties of the oxidic high T(sub c) superconductors are significantly affected by the presence of minor amounts of various elements adventing as impurities, e.g., from the chemical environment during manufacturing. YBa2Cu3O(9-delta) is prone to an extinction of the superconductivity on (partial) substitution of all four elemental components. E.g., Pr (for Y), La (for Ba), Zn (for Cu) or peroxygroup (for O) substituents will alter some of the superconductivity preconditions, like mixed valence state in Cu3O7/O(9-delta) network or structural distortion of the network. Although various pseudoternary chemical equilibrium phase diagrams of the Y(O)-Ba(O)-Cu(O) system now are available, no consensus is generally shown, however, this is partly due to lack of compatible definitions of the equilibrium conditions. Less information is available about the phase compatibilities in the appropriate quaternary phase diagram (including oxygen) and virtually no information exists about any pentenary phase diagrams (including one impurity). Unfortunately, complexity of such systems, stemming both from number of quaternary or pentenary compounds and from visualizing the five-component phase system, limits this presentation to more or less close surroundings of the YBa2Cu3O(9-delta) type phase in appropriate pseudoquaternary or pseudopseudoternary diagrams, involving Y-Ba-Cu and O, O-CO2, alkaline metals, Mg and alkaline earths, and Sc and most of the 3-d and 4-f elements. The systems were investigated by means of x ray diffraction, neutron diffraction and chemical analytical methods on samples prepared by sol-gel technique from citrates. The superconductivity was characterized by measuring the diamagnetic susceptibility by SQUID.

Karen, P.↗
DOE OSTI2023

High-temperature chemical expansion of Pr0.1Ce0.9O2-δ thin films determined by Differential Laser Doppler Vibrometry

The chemical expansion of thin-film praseodymium-cerium oxide solid solutions (Pr 0.1 Ce 0.9 O 2-δ , PCO10) is investigated at temperatures ranging from 500 °C to 800 °C and oxygen activities down to about 10 -22 . Differential Laser Doppler Vibrometry (D-LDV) in combination with periodic electrochemical pumping of oxygen is used to detect very small displacements of the sample surface. With decreasing pumping frequency, the film approaches chemical equilibrium, which is nearly reached at 700 °C below 0.01 Hz. D-LDV enables direct measurement of film thickness changes as well as the bending of the underlying substrates due electrochemically induced lateral mechanical stress. The substrate bending is found to be more than an order of magnitude larger than the film thickness change, increasing the detection limit by the same factor. The thickness change of the thin film is calculated from the measured displacement of the sample surface and compared with directly measured values. Thin-film literature data available for low oxygen activities up to 650 °C agree with our measurements. In addition, the thickness change of the film is determined from the displacement at different sample positions, which validates previous measurements. Finally, thin film chemical expansion is found to be larger by a factor of about 2 than the corresponding bulk value in agreement with predictions from a mechanical model that the out-of-plane expansion of a constrained film is increased.

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

Nonequilibrium effects on stability of hybrid stars with first-order phase transitions

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

On the Composition of Young, Directly Imaged Giant Planets

The past decade has seen significant progress on the direct detection and characterization of young, self-luminous giant planets at wide orbital separations from their host stars. Some of these planets show evidence for disequilibrium processes like transport-induced quenching in their atmospheres; photochemistry may also be important, despite the typically large orbital distances. Disequilibrium chemical processes such as these can alter the expected composition, spectral behavior, thermal structure, and cooling history of the planets, and can potentially confuse determinations of bulk elemental ratios, which provide important insights into planet-formation mechanisms. Using a thermo/photochemical kinetics and transport model, we investigate the extent to which disequilibrium chemical processes affect the composition and spectra of directly imaged giant exoplanets. Results for specific "young Jupiters" such as HR 8799 b and c and 51 Eri b are presented, as are general trends as a function of planetary effective temperature, surface gravity, incident ultraviolet flux, and strength of deep atmospheric convection. We find that quenching is very important on young Jupiters, leading to CO/CH4 and N2/NH3 ratios much greater than; and H2O mixing ratios a factor of a few less than chemical equilibrium predictions. Photochemistry can also be important on such planets, with CO2 and HCN being key photochemical products. Carbon dioxide becomes a particularly major constituent when stratospheric temperatures are low and recycling of water following H2O photolysis becomes stifled. Young Jupiters with effective temperatures less than 700 degrees Kelvin are in a particularly interesting photochemical regime that differs from both transiting hot Jupiters and our own solar-system giant planets.

photochemistry↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

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

From Langmuir isotherm to Brunauer–Emmett–Teller isotherm

Abstract Developed for multilayer adsorption, the Brunauer–Emmett–Teller ( BET) isotherm considers the adsorption of the first layer as an equilibrium chemical reaction between adsorbate molecules and adsorption sites and the adsorption of the second and subsequent layers as a condensation–evaporation process. Following the recent development of an activity‐based formulation for the Langmuir isotherm for monolayer adsorption, we present an activity‐based formulation for the BET isotherm in which species concentrations are replaced with species activities. Capturing the adsorbent surface heterogeneity for the adsorption of the first layer, the resulting thermodynamic BET isotherm is shown to accurately represent pure component adsorption isotherms over the relative pressure range of zero to unity or prior to the onset of capillary condensation. The thermodynamic BET isotherm should facilitate accurate estimation of monolayer adsorption capacity and the corresponding adsorbent surface area.

Vyawahare, Pradeep↗
DOE OSTI2021

Repetitively pulsed nanosecond discharge plasma decay in propane–oxygen gas mixture in the presence of a heating electric field

Plasma decay in the afterglow of a repetitively pulsed nanosecond discharge in a stoichiometric propane–oxygen mixture was experimentally investigated when a weak heating DC electric field was applied and in its absence. The discharge was ignited at room gas temperature and a pressure of 1–2 Torr and was characterized by low specific energy inputs (<0.004 eV per molecule in one pulse). Using microwave interferometry, the temporal evolution of the electron density during plasma decay was studied, and the effective recombination coefficients were obtained from data processing. It was shown that the rate of plasma decay behaved in a non-monotonic manner with increasing degree of propane oxidation; at first the decay rate grew, then passed through a maximum, fell and saturated in the limit of a large (~2000) number of pulses. In this limit, the effect of the heating DC electric field on the plasma decay decreased with approaching chemical equilibrium. Numerical simulation of the observed effects was performed for low and high oxidation degrees of propane taking into account changes in the composition of positive ions in the plasma. Good agreement was obtained between measurements and calculations of the electron density during plasma decay in these cases. Here it was shown that the formation of cluster ions in the discharge afterglow plays a fundamental role. The plasma decay was controlled by electron recombination with hydrocarbon cluster ion at low oxidation degree of propane and with water cluster (hydrated) ions at high oxidation degree. A hypothesis was proposed to explain the observed nonmonotonic behavior of the plasma decay rate with an increase in the propane oxidation in the discharge, based on the formation of hydrated hydrocarbon ions C x H y + (H 2 O) k at moderate oxidation degrees.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗