Atomistic computer simulations of dissolution of inorganic glasses.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Irradiated light water reactor (LWR) fuel contains material which is insoluble in nitric acid, and process solutions are unstable with respect to precipitation of compounds insoluble in nitric acid. Accordingly, solids will be encountered during fuel reprocessing, and operations must be conducted in a way to minimize formation of additional solids and to prevent their interference with reprocessing operations.
The dissolution rates of the aluminum alloys in the High Flux Isotope Reactor (HFIR) element carriers and the Material Test Reactor (MTR) L-bundles in the H-Canyon facility have been identified as the possible cause of extended dissolutions that result in significant time and financial expenditures. A study, carried out by Savannah River National Laboratory (SRNL) to determine relationships between the dissolution rates and the metallurgical properties of the aluminum alloy materials of construction of the HFIR carriers and the L-bundles, considered the dissolution rates of aluminum alloy (AA) series 1100, 6061, and 6063. The study determined that the aluminum alloy compositions played a principal role in the dissolution rate of the carrier/bundle components. Higher dissolution rates were correlated with lower concentrations of the minor element additions in the alloys and with specific element concentrations. Aluminum alloys 1100 and 6063 were found to have similar dissolution rates that were approximately two orders of magnitude (100X) greater than those of AA6061. Based on the results of the dissolution behavior study, a Technical Assistance Request (TAR) was first issued to determine if the replacement of AA6061-T6 with AA6063-T6 is feasible for the HFIR carrier lifting bails. A Technical Task Request was then issued to consider AA6063-T5 as well as other alloys to improve possible supply chain issues. The metallurgical properties of the L-bundle (specifically the end caps) were not evaluated in this report because L-Bundle drawings already allow for the use of AA6063-T6 in all structural components. The HFIR carriers are composed of thin-walled components with significant surface areas that allow for relatively quick overall dissolution times. Conversely, the carrier lifting bails and the supporting constituents are composed of solid bars and thick plate regions with relatively small surface areas that experience longer overall dissolution times. While the MTR L-bundle design includes allowances for the materials of construction to be either AA6061-T6 or AA6063-T6, the HFIR carriers are specified to be constructed fully with AA6061-T6 alloy. This report analyzes the recommendations of the dissolution behavior study to replace the materials of construction of the HFIR carrier lifting bails. The analysis considers the operational requirements of the lifting bail and its supporting structures. To decrease dissolution times, the analysis considers direct replacement of the material as well as reductions in the thicknesses of the components to decrease the mass of the elements. Material reductions are considered on options for using either AA6061 and/or AA6063. The calculations are based on specifications from the American Society of Mechanical Engineer (ASME) and The Aluminum Association, Inc. design codes. The analysis finds that direct replacement of the lifting bail material of construction with AA6063-T6, and AA6063-T5 as well as reductions in the dimensions of the lifting bail components are acceptable. Note that this study considers the structural suitability of the alloys. It does not consider their dissolution rates in the dissolvers.
The effect of immobile pore-water on gas transport in fractured rock has implications for numerical modeling of soil vapor extraction, methane leakage, gaseous CO 2 leakage from sequestration operations, radionuclide transport from underground nuclear explosions, and nuclear waste disposal. While the ability for immobile pore-water storage to effect gas transport has been recognized in the past, the details and specific scenarios leading to enhanced, retarded, or unaffected gas transport have not been explored. In this paper, we performed numerical investigations into the enhancement and retardation of gas transport due to immobile pore-water storage in order to identify implications for gas transport applications. To do this, we developed a numerical approach to model gas transport with a single-phase flow solution coupled to the advection–dispersion equation modified to account for immobile pore-water storage. Other than the immobility of pore water, the formulation contains all other physics included in two-phase formulations (advective and diffusive gas transport in fractures and rock matrix and dissolution in immobile pore water). The assumption of immobile pore water is valid here since for many applications involving transport of soluble gases in fractured rock, the rate of aqueous transport is insignificant compared to gas transport. We verify our modeling approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) barometric pumping of a fracture, and (4) gas transport with uniform fracture flow and diffusion into the matrix. We account for pore-water storage in our model by implementing a kinetic formulation of gas dissolution wherein the dissolved (aqueous) phase is considered an immobile constituent. Using this formulation, we model the effect of dissolution rate and saturation on the retardation of gas transport during pure diffusion and pure advection. We also demonstrate that although it is commonly believed that pore-water storage will always enhance gas transport in fractures during oscillatory flow (e.g., during reversing pressure gradients such as barometric pumping cycles), our simulations indicate that this may not always be the case. Our numerical investigations indicate that scenarios with lower effective diffusion coefficients ( ⪅<!-- <≈ --> 10 −<!-- - --> 5  m 2 / s ) and lower dissolution coefficients (i.e., the dissolution diffusion coefficient) ( ⪅<!-- <≈ --> 10 −<!-- - --> 11  m 2 / s ) will result in enhanced gas transport. Other combinations of gaseous diffusion and dissolution coefficients result in delayed gas transport or insignificant effect on gas transport. Based on these results, and given the fact that the free-air diffusion coefficients for many gases are slightly above 10 −<!-- - --> 5  m 2 / s (near the boundary of enhanced/non-enhanced gas transport), tortuosity would have to be significant for gas transport enhancement to occur. Similarly, if we consider the free-water diffusion coefficient of gases (around 10 −<!-- - --> 9  m 2 / s for many gases of interest) to be a maximum bound for the effective dissolution coefficients (i.e., the rate of diffusion into the water away from the air/water interface will limit the dissolution rate), the effective dissolution coefficient would have to be at least 2 orders of magnitude less for enhanced gas transport to occur. Otherwise, pore-water storage will delay or have negligible effect on the gas transport. The results also indicate that the assumption of instantaneous equilibrium, often invoked in numerical codes for dissolution processes (i.e., dissolution coefficient is effectively infinite), will fail to accurately capture important details of soluble gas transport. Finally, the implication of our results is that gas breakthrough times in fractured rock during oscillatory flow may be highly dependent on the effective gaseous diffusion and dissolution coefficients.
Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.
Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.
Epsilon metal (e-metal) is the metallic phase that forms as inclusions at the grain boundaries in the UO2 fuel during reactor operation. This metal is composed of Pd, Mo, Rh, Ru, and Tc. These metallic inclusions are insoluble in strong acid and remnants of these metallic inclusions have been found in the UO2 matrix that remains from the natural reactors in Gabon that were active 1.8 billion years ago, therefore e-metal should be an excellent waste form for the immobilization of the long-lived isotopes 107Pd (6.5 ? 106 a) and 99Tc (2.13 ? 105 a), with 99Tc being the isotope of interest for repository performance. Therefore, the chemical durability of this potential waste form is assessed in this study. Typically, corrosion rates for metallic materials are measured electrochemically because they are quick, inexpensive, and can reveal the mechanism by which a metal corrodes, at least initially. However, in a repository the waste form would be subjected to slowly flowing water without an applied electrical potential over long time periods. Therefore, the corrosion rates of e-metal specimens were measured with both electrochemical tests and the single-pass flow-through test (SPFT). Potentiodynamic and potentiostatic polarization results suggest that a thin passive film exists on the alloy surface, which seems to be responsible for its high corrosion resistance. Additionally, X-Ray photoelectron spectroscopic results suggest that Pd oxides are significantly enriched in the passive film Results from the SPFT show that the dissolution rates were weakly dependent on pH. Only Mo and Re were found in solution and were used for the calculation of the dissolution rates. In general, the electrochemically determined corrosion rates agree reasonably well with the initial dissolution rate measured with the SPFT test, but they are about one or two orders of magnitude higher than the steady state rates. The causes for this discrepancy are discussed.
A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.
The interfacial electron transfer (ET) between electron shuttling compounds and iron (Fe) oxyhydroxides plays a crucial role in the reductive dissolution of Fe minerals and the fate of surface-bound arsenic (As). However, the impact of exposed facets of highly crystalline hematite on reductive dissolution and As immobilization is poorly understood. In this study, we systematically investigated the interfacial processes of the electron shuttling compound cysteine (Cys) on various facets of hematite and the reallocations of surface-bound As(III) or As(V) on the respective surfaces. Our results demonstrate that the ET process between Cys and hematite generates Fe(II) and leads to reductive dissolution, with more Fe(II) generated on {001} facets of exposed hematite nanoplates (HNPs). Reductive dissolution of hematite leads to significantly enhanced As(V) reallocations on hematite. Nevertheless, upon the addition of Cys, a raipd release of As(III) can be halted by its prompt re-adsorption, leaving the extent of As(III) immobilization on hematite unchanged throughout the course of reductive dissolution. This is due to that Fe(II) can form new precipitates with As(V), a process that is facet-dependent and influenced by water chemistry. Electrochemical analysis reveals that HNPs exhibit higher conductivity and ET ability, which is beneficial for reductive dissolution and As reallocations on hematite. In conclusion, these findings highlight the facet-dependent reallocations of As(III) and As(V) facilitated by electron shuttling compounds and have implications for the biogeochemical processes of As in soil and subsurface environments.
In the event of a nuclear detonation, actinide and fission product-containing debris can be analyzed to provide information about the nature of the device. In-field sample dissolution prior to transport to a laboratory is critical for reducing the time between an event and delivery of information to decision makers. However, the best dissolution methods utilize liquid strong acid solutions, which are incompatible with field collection Standard Operating Procedure. What is required are solid materials that, when mixed with water in the field, can provide the strong digestion and dissolution properties of nitric and hydrofluoric acids. Fluorohydrogenate ionic liquids have shown a remarkable ability to dissolve a wide range of materials, from refractory actinide oxides to cementitious material such as calcium silicates. These strong dissolution properties are due to the neutral hydrofluoric acid that is bound to the anion component of the ionic liquid. Because these materials are liquids, they are also incompatible with Standard Operating Procedure. However, the physiochemical properties of ionic liquids, such as the melting point, can be readily tuned through minor modification of the molecular structure of the cations and anions. Thus, we hypothesize that by synthesizing modified cations, novel ionic compounds with melting points above ambient but with the same anion and bound hydrofluoric acid – and thus similar digestion and dissolution properties – as current fluorohydrogenate ionic liquids can be produced. These solid materials would be compatible with current Standard Operating Procedure, and would allow for rapid, in-field dissolution of samples, shortening the time to actionable information after an event.
Relative solvating power, that is, the ratio of the coordination ratios between a solvent and the reference solvent, was used to probe the quantitative structure–activity relationship of electrolyte solvents and the lithium polysulfide (LiPS) dissolution in lithium–sulfur batteries. Internally referenced diffusion-ordered nuclear magnetic resonance spectroscopy (IR-DOSY) was used to determine the diffusion coefficient and coordination ratio, from which the relative solvating power can be easily measured. The higher the relative solvating power of an ethereal solvent, the more severe will be the LiPS dissolution and the lower the coulombic efficiency of the lithium–sulfur battery. A linear relationship was established between the logarithm of relative solvating power of a solvent and the degree of LiPS dissolution, rendering relative solvating power an important parameter in choosing the electrolyte solvent for lithium–sulfur batteries.