Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Structure”

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 1,171 records · Page 65

Compositionally-driven surface nanostructuring on refractory compositionally complex alloys under low energy helium bombardment

Additive manufacturing enables user-defined control of the compositional complexity, opening new design spaces for complex concentrated alloys (CCAs). Refractory CCAs may offer enhanced performance in the divertor region of a fusion reactor environment where plasma-facing materials will be subject to high temperatures, low energy He and D particles, and 14 MeV neutrons. In this work, specimens with the nominal composition of NbTaMoTi, NbTaMo, NbTaTi, and NbTa were fabricated via directed energy deposition (DED) then bombarded with 40 eV He ions to a fluence of 2x10 26 m −2 at ∼ 1000 K. Post irradiation, the surface morphology and composition were examined with electron microscopy and x-ray photoelectron spectroscopy to offer information on the spatial distribution of surface nano-structuring. Sub-surface He bubbles driving the nano-structuring were examined with electron microscopy techniques. Analysis indicates the local composition directly influences the He bubble and tendril size, while the region with the highest complexity showed the shortest nano-structuring. Molecular Dynamics simulations complement the experimental results, showing comparable helium bubble growth and migration as a function of compositional complexity. This work demonstrates the compositional dependence of surface nanostructure formation of compositionally complex alloys, important for future design of complex plasma facing materials in fusion reactors.

Complex concentrated alloys↗

Reactor Physics Simulations of the High Flux Isotope Reactor Permanent Beryllium Reflector Number 5 Design

This study evaluates reactor physics aspects of the proposed High Flux Isotope Reactor (HFIR) permanent beryllium reflector number 5 design. HFIR is a pressurized, light water–cooled, and beryllium reflected research reactor that operates at 85 MW for cycle lengths of approximately 24 to 26 days. A new permanent reflector design is highly desirable for increased versatility in irradiation experiments, to arrange the vertical experiment facilities (VXFs) to minimize their impact on neutron scattering if loaded with neutronabsorbing experiments, to enhance reflector thermal-structural-hydraulic performance, and to simplify the complex fabrication process. The current number 4 and concept number 5 reflectors have 22 and 28 VXFs, respectively, and are modeled and analyzed in this study with the MCNP5 and SCALE 6.1.3 nuclear simulation codes.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effect of nitrogen on tensile properties and structures of T-111 (tantalum, 8 percent tungsten, 2 percent hafnium) tubing

The effect of controlled nitrogen additions was evaluated on the mechanical properties of T-111 (Ta-8W-2Hf) fuel pin cladding material proposed for use in a lithium-cooled nuclear reactor concept. Additions of 80 to 1125 ppm nitrogen resulted in increased strengthening of T-111 tubular section test specimens at temperatures of 25 to 1200 C. Homogeneous distributions of up to 500 ppm nitrogen did not seriously decrease tensile ductility. Both single and two-phase microstructures, with hafnium nitride as the second phase, were evaluated in this study.

Buzzard, R. J.↗

Effect of energetic ions on edge-localized modes in tokamak plasmas

The most efficient and promising operational regime for the International Thermonuclear Experimental Reactor tokamak is the high-confinement mode. In this regime, however, periodic relaxations of the plasma edge can occur. These edge-localized modes pose a threat to the integrity of the fusion device. Here we reveal the strong impact of energetic ions on the spatio-temporal structure of edge-localized modes in tokamaks using nonlinear hybrid kinetic–magnetohydrodynamic simulations. A resonant interaction between the fast ions at the plasma edge and the electromagnetic perturbations from the edge-localized mode leads to an energy and momentum exchange. Energetic ions modify, for example, the amplitude, frequency spectrum and crash timing of edge-localized modes. The simulations reproduce some observations that feature abrupt and large edge-localized mode crashes. The results indicate that, in the International Thermonuclear Experimental Reactor, a strong interaction between the fusion-born alpha particles and ions from neutral beam injection, a main heating and fast particle source, is expected with predicted edge-localized mode perturbations. This work advances the understanding of the physics underlying edge-localized mode crashes in the presence of energetic particles and highlights the importance of including energetic ion kinetic effects in the optimization of edge-localized mode control techniques and regimes that are free of such modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SEM/EBSD Characterization of Baffle-Former Bolts Designated for High Energy X-ray Microscopy Studies at the Advanced Photon Source

This report documents the microstructural and surface characterization of two high-fluence baffle-former bolts (BFBs) harvested from a commercial pressurized water reactor (PWR) after nearly four decades of service. These bolts experienced extreme conditions including high neutron fluence, elevated temperatures, corrosive primary water, and sustained mechanical loads, all of which contribute to degradation mechanisms such as irradiation-assisted stress corrosion cracking (IASCC). Using scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), the study assessed bulk and near-surface microstructures, grain boundary characteristics, and in-service oxide layers. The results revealed significant differences between bolts exposed to 41 dpa and 29 dpa, including differences in recrystallization extent, deformation structures, and oxidation morphology. These findings establish critical baseline data for integration with advanced synchrotron-based high-energy X-ray microscopy (HEXM) and computational crystal plasticity models. Ultimately, the work supports improved mechanistic understanding and predictive capability regarding the long-term performance of reactor core internal components.

Lach, Tim [Oak Ridge National Laboratory (ORNL), O↗

Quantification of silicon carbide grain structure in TRISO fuel by BSE image analysis

The grain structure of the silicon carbide (SiC) layer in tristructural-isotropic (TRISO) fuels has commonly been qualified in fuel specifications based on a qualitative comparison to visual standards to ensure undesirable microstructure are not produced. This approach is inherently dependent on subjective judgement and provides limited data for analysis or comparison. An alternative method to rapidly quantify (i.e., provide a numeric measure of) the grain structure of the SiC layer in TRISO fuels has been developed. This method relies on the identification of silicon carbide grain boundaries based on brightness variation in backscattered electron (BSE) images through the application of automated image processing algorithms. The speed with which the requisite BSE imaging and analysis may be completed enables the practical use of this method on an industrial scale as a part of qualifying TRISO fuel to meet a given specification on silicon carbide grain structure, which is an important parameter for TRISO fuel performance. The method is benchmarked by comparison to prior Electron-Backscattered Diffraction (EBSD) results for TRISO particles from the Advanced Gas Reactor (AGR) program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Measuring thermal diffusivity and gap conductance in uranium nitride and Zircaloy relevant for microreactor applications

Heat transfer across nuclear fuels and structural interfaces is an important factor for evaluating the performance of nuclear power systems. Specifically, heat generated as nuclear fuel fissions must be transported through the cladding material and through the reactor to reach the steam turbine for power generation. As new microreactor designs emerge, maximizing the efficiency of this heat transfer process becomes crucial to make them commercially viable. This article examines thermal diffusivity and gap conductance in uranium nitride (UN) fuel and Zircaloy-4 (Zry4) cladding using light flash analysis (LFA). Thermal diffusivity measurements were made on monolithic UN pellets and Zry4 exposed to carbon at peak operating temperatures of microreactors and show that carbon ingress has a minimal effect on thermal diffusivity when compared with identical materials not exposed to carbon. Evaluation of gap conductance at the UN-Zry4 interface was done using one-dimensional two-layer thermal transport models as a function of applied pressure. Here the results show that increasing pressure on the UN-Zry4 interface leads to gains in gap conductance per unit area in fuel-cladding assemblies at microreactor operating temperatures. While many other variables are expected to influence UN-Zry4 interfacial gap conductance (e.g. contact surface roughness, porosity, localized heating, environmental gas pressure), the work offers a demonstration of using a conventional LFA apparatus to determine this parameter at elevated temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Materials Assessment for FLiBe Fusion Blankets

This evaluation was prepared at the request of DOE-FES and expands on an earlier evaluation prepared by Wilson, Wiffen and Keiser in 2015 (attached as Appendices A and B). Additional useful information for materials in fission molten salt reactor systems that use FLiBe can be found in a prior NRC report [Busby 2019]. The purpose of this report is to review the considerations that will determine the fusion blanket structural materials that will be suitable for use in a FLiBe breeder-coolant fusion power plant (FPP). Identifying a fusion-relevant structural material that is sufficiently compatible with FLiBe is a key issue to be addressed. The constraints and trade-offs made in structural materials selection include the usual power generation requirements such as mechanical properties, fabricability and durability, which assumes reasonable compatibility with the working fluid. For an FPP, additional constraints include neutron economy (breeding ratio), resistance to irradiation damage and low activation (i.e. waste disposal) characteristics.

36 MATERIALS SCIENCE↗

Silicon Halide-alkali Metal Flames as a Source of Solar Grade Silicon

A program is presented which was aimed at determining the feasibility of using high temperature reactions of alkali metals and silicon halides to produce low cost solar-grade silicon. Experiments are being conducted to evaluate product separation and collection processes, measure heat release parameters for scaling purposes, and determine the effects of the reactants and/or products on materials of reactor construction. During the current reporting period, the results of heat release experiments were used to design and construct a new type of thick-wall graphite reactor to produce larger quantities of silicon. A reactor test facility was constructed. Material compatibility tests were performed for Na in contact with graphite and several coated graphites. All samples were rapidly degraded at T = 1200K, while samples retained structural strength at 1700K. Pyrolytic graphite coatings cracked and separated from substances in all cases.

Olson, D. B.↗

Effects of irradiation damage on the elastic properties and hardness of complex transition metal diborides

Multi-principal component transition metal (TM) diborides represent a class of high-entropy ceramics (HECs) that have received considerable interest in recent years owing to their promising properties for extreme environment applications that include thermal/ environmental barriers, hypersonic vehicles, turbine engines, and next-generation nuclear reactors. While the addition of chemical disorder through the random distribution of TM elements on the cation sublattice has offered opportunities to tailor elastic stiffness and hardness, the effects of irradiationinduced structural damage on the physical properties of these complex materials have remained largely unexplored. To this end, changes in the hardness and elastic moduli of a high-entropy TM diboride (Hf0.2Nb0.2Ta0.2Ti0.2Zr0.2)B2 and three of its quaternary subsets following irradiation with 10 MeV gold (Au) ions to fluences of up to 6 × 10^15 Au cm-2 are investigated at the micrometer and sub-micrometer length-scales via the dispersion of laser-generated surface acoustic waves (SAW) and nanoindentation, respectively. The nanoindentation measurements show that the TM diborides exhibit an initial increase in hardness following irradiation with energetic Au ions, with a subsequent decrease in hardness following further irradiation. One quaternary composition, (Hf1/3Ta1/3Ti1/3)B2, exhibits a notable exception to the trend and continues to exhibit an increase in hardness with ion irradiation fluence. Although differences in the absolute values of the effective elastic moduli obtained from the measured SAW dispersion and nanoindentation are observed (and attributed to microstructural variations at the measurement length-scale), both techniques yield similar trends in the form of an initial reduction and subsequent saturation in the elastic modulus with increasing ion irradiation fluence. The quaternary TM diboride (Hf1/3Ta1/3Ti1/3)B2 again exhibits a departure from this trend. The high-entropy TM diboride (Hf0.2Nb0.2Ta0.2Ti0.2Zr0.2)B2 exhibits the greatest recovery in hardness and modulus when irradiated to high ion fluences following initial changes at low fluence, indicating superior resistance to radiation-induced damage over its quaternary counterparts. Opportunities for designing HECs with superior hardness and modulus for enhanced radiation resistance (compared to their single constituent counterparts) by tailoring chemical disorder and bond character in the lattice are discussed.

36 - MATERIALS SCIENCE↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Strength stability at high temperatures for additively manufactured alumina forming austenitic alloy

Several fast-spectrum nuclear reactors designed to generate high power (~450 MWe) rely on forced convection of media such as supercritical CO 2 , sodium, or liquid lead to cool the nuclear core, operating at temperatures up to 600 °C. Cost-effective, high-strength Fe-based alumina forming austenitic (AFA) alloys are a promising candidate for the fabrication of critical nuclear components. This study investigated laser powder bed fusion (LPBF) processing of an AFA alloy composition optimized for improved creep resistance. Electron microscopy revealed an elongated grain structure along the build direction with a fine sub-grain cellular structure decorated with (Cr,Fe,Nb) 23 C 6 carbide precipitates at the intercellular boundaries. Finally, at temperatures of 20–900 °C, the LPBF alloy's superior tensile properties compared to its arc-melted counterpart and other advanced steels (e.g., SS316) were attributed to the distribution of nano-sized carbide precipitates, whereas the high ductility was attributed to the LPBF alloy's elongated grain structure.

36 MATERIALS SCIENCE↗

Methods development towards automated, physics-informed, quantitative quality control of TRISO-SiC

Tristructural-isotropic (TRISO) fuel particles have been developed as a high-performance fuel for use in high-temperature gas-cooled reactor (HTGR) systems due to their high efficiency and stability under both normal and off-normal conditions. Broader deployment of this technology in advanced nuclear applications may benefit from quantitative quality assurance and quality control (QA/QC) methods that directly link TRISO properties to downstream performance. A key TRISO property is the SiC layer microstructure, which influences fission product retention during irradiation. However, existing QA/QC for the TRISO-SiC microstructure comprises only a qualitative visual inspection; therefore, there is a clear opportunity for the development of quantitative methods for TRISO QA/QC. Here, to this end, previous work has demonstrated an image processing approach to grain boundary (GB) identification and subsequent extraction of microstructural metrics; however, extensive twinning within the SiC layer complicates such analyses because twin GBs significantly influence microstructural metrics but are not expected to contribute to fission product transport. This study presents the initial development, training, and testing of an ML-based image segmentation algorithm designed to identify and remove twin GBs from standard backscattered electron micrographs, providing an industrially applicable, quantitative, and physically meaningful QA/QC approach for the TRISO-SiC microstructure. Although pixel-wise performance metrics for the twin predictions are low, the change in grain area and the number of GB pixels after twin removal predicted by the ML workflow are within 1% of the true values calculated using crystallographic data. This suggests that the model is well capable of predicting overall twin boundary structures and grain morphology, and continued advancement of this approach could enable automated, scalable, and physics-informed QA/QC for TRISO-SiC microstructures, supporting the reliable qualification of coated particle fuels for next-generation reactor systems.

Computer vision↗

Scattering of Radio Frequency Waves by Density Filaments

The edge region and the scrape-off layer of magnetically confined fusion devices, like tokamaks and stellarators, are replete with turbulent plasma that is a mixture of coherent, blob or filament like, structures [1] and incoherent fluctuations [2]. The variation in the density due to turbulence can be comparable to or greater than the ambient density [2]. As part of an overall effort to optimize the efficiency of operation, radio frequency (RF) waves are commonly used for heating fusion plasmas and, in tokamaks, for generating plasma current needed for confinement and controlling instabilities. The RF waves are excited by antenna structures that are placed near the wall of a fusion device. In order to deliver energy and momentum to charged particles in the core of fusion plasmas, RF waves have to propagate through the turbulent plasma. In present fusion devices, the scrap-off layer and the edge plasma region is of the order of a few centimeters. In reactor type devices, like ITER, this region is expected to be of the order of tens of centimeters. Since the efficiency of operation of a fusion reactor is of prime importance, it is imperative that we understand the effect of turbulence on RF waves. The fluctuations in density lead to changes in the plasma permittivity. As in conventional electrodynamics, the propagation of RF waves through different dielectric media is subject to reflection, refraction, and diffraction. In this paper, we summarize our theoretical and computational studies on the propagation of RF waves through filamentary structures present in the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of alkaline carbonate cooling to mitigate Ex-Vessel molten corium accidents

To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. Here the novel SM proposed is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO 2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of superlattice device structures

This report describes the investigation of growth properties, and the structure of epitaxial multilayer Si(Si(1x)Ge(x)) films grown on bulk Silicon Substrates. It also describes the fabrication and characterization of MOSFET and MESFET devices made on these epitaxial films. Films were grown in a CVD reactor using hydrides of Si and Ge with H2 and He as carrier gases. Growth temperatures were between 900 C and 1050 C with most films grown at 1000 C. Layer thickness was between 300A and 2000A and total film thickness was between 0.25 micro m and 7 micro m. The Ge content (X) in the alloy layers was between .05 and 0.2. N-type multilayer films grown on (100) p-type Si showed Hall mobility in the range 1000 to 1500 sq cm/v for an average carrier concentration of approx. 10 to the 16th power/cu cm. This is up to 50% higher than the Hall mobility observed in epitaxial Si films grown under the same conditions and with the same average carrier concentration. The mobility enhancement occurred in films with average carrier concentration (n) from 0.7 x 10 to the 16th power to 2 x 10 to the 17th power/cu cm, and total film thickness greater than 1.0 micro m. No mobility enhancement was seen in n-type multilayer films grown on (111) Si or in p-type multilayer films. The structure of the films was investigated was using SEM, TEM, AES, SIMS, and X-ray double crystal diffraction techniques. The film composition profile (AES, SIMS) showed that the transition region between layers is of the order of about 100A. The TEM examination revealed a well defined layered structure with fairly sharp interfaces and good crystalline quality. It also showed that the first few layers of the film (closest to the substrate) are uneven, most probably due to the initial growth pattern of the epitaxial film where growth occurs first in isolated islands that eventually growth and coalesce. The X-ray diffraction measurement determined the elastic strain and strain relief in the alloy layers of the film and the elastic strain in the intervening Si layers.

Gergis, I. S.↗

CaMn 0.9 Ti 0.1 O 3 based redox catalysts for chemical looping – Oxidative dehydrogenation of ethane: Effects of Na 2 MoO 4 promoter and degree of reduction on the reaction kinetics

Reduction kinetics and stability of 20 wt% Na 2 MoO 4 -promoted CaMn 0.9 Ti 0.1 O 3 were investigated for its applications in Chemical Looping – Oxidative Dehydrogenation (CL-ODH) of ethane, a potential alternative for ethylene production with higher efficiency and lower emissions. Here, the present work reports a kinetics model and parameters for a Na 2 MoO 4 -promoted, Ti-doped CaMnO 3 (CaMn 0.9 Ti 0.1 O 3 ) redox catalyst under H 2 and C 2 H 4 . A first-order reaction model provides the best fit for the reduction of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 , while the C 2 H 4 reduction is well described by an Avrami–Erofe’ev model. The activation energy for C 2 H 4 oxidation is approximately three times higher than that for H 2 conversion, showing that the activation of C 2 H 4 is significantly more difficult on the surface of the redox catalyst. The reduction rate of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 at 750 °C is more than two orders of magnitude greater than that under C 2 H 4 , while the reduction rate of unpromoted CaMn 0.9 Ti 0.1 O 3 is comparable under H 2 and C 2 H 4 , showing that the addition of Na 2 MoO 4 effectively suppresses C 2 H 4 combustion relative to H 2 oxidation. The kinetics results for Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 confirm its excellent selectivity towards hydrogen combustion, making it a promising candidate under CL-ODH. Additionally, the stability of the CaMn 0.9 Ti 0.1 O 3 @ Na 2 MoO 4 core-shell structure, which was the underlying reason for the excellent selectivity, was examined under both shallow and deep reductions. It was determined that deep reduction of the redox catalyst, e.g. higher than 80% solid conversion, would lead to loss of sodium and hence to decreased selectivity for hydrogen combustion. In contrast, the core-shell structure was well-maintained, exhibiting excellent performance after 50 redox cycles when deep reduction of the redox catalyst was avoided. This study offers a basis for both the CL-ODH reactor design and redox catalyst optimizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗