Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Batch”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A direct transcription-based multiple shooting formulation for dynamic optimization

The growing need for fast and efficient solution techniques for solving dynamic optimization problems is driven by a broad spectrum of applications in scheduling and control. We suggest a novel framework for dynamic optimization that utilizes a multiple shooting “backbone” with discrete rather than continuous subproblems, thereby eliminating need for repeated time-integration. A Lagrangian relaxation (LR)-based decomposition scheme is proposed, which dualizes the state continuity requirements between subproblems and enables parallel solution of the problem. We demonstrate the applicability of our method on two case studies: the Van der Pol oscillator and a batch reactor.

42 ENGINEERING↗
DOE OSTI2024

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Universal method for the optimization of HDC coating uniformity on non-planar, non-stationary substrates for inertial confinement fusion targets

The thickness uniformity of chemical vapor deposited (CVD) diamond coatings on non-planar, non-stationary substrates depends on both the intrinsic instantaneous coating thickness distribution (ICTD) of the coating conditions used and, if applicable, on the frequency of substrate reorientation. While important for many CVD diamond applications, the relative impact of the ICTD and substrate reorientation on the coating thickness uniformity has not been studied. In this work, we systematically investigate the effect of these factors for microwave-plasma chemical vapor deposition (MPCVD) of diamond (referred to as high density carbon (HDC) in the inertial confinement fusion (ICF) community) coatings on spherical, rolling substrates. This coating technique is used to fabricate capsules for ICF experiments, which require extreme coating uniformity with <0.3 % thickness variation (so-called Mode 1 or M1) to ensure symmetric compression of imploding targets. To extract the otherwise unobservable reorientation timescale (Δt), Monte Carlo simulations were performed using experimental ICTD data as input. This combined approach confirms scaling relationships between the substrate reorientation timescale as well as coating thickness and coating uniformity, as expected from a 3D random walk. Simulations confirm that M1 is Rayleigh-distributed and scales as (Δt) 1/2 , consistent with the randomization of two angles that determine orientation of a sphere. We also demonstrate that, under the conditions studied, Δt is the dominant factor in determining thickness uniformity while the intrinsic ICTD has minimal impact. Finally, experiments show that Δt can be affected by total batch size under constant agitation conditions due to space constraints that limit the capsule reorientation kinetics. In conclusion, this study highlights the utility of a combined experiment-simulation approach as a general methodology for understanding and improving coating uniformity on non-planar, non-stationary substrates.

Capsule↗
DOE OSTI2025

Artificial intelligence driven laser parameter search: Inverse design of photonic surfaces using greedy surrogate-based optimization

Photonic surfaces designed with specific optical characteristics are becoming increasingly crucial for novel energy harvesting and storage systems. The design of these surfaces can be achieved by texturing materials using lasers. The optimal adjustment of laser fabrication parameters to achieve target surface optical properties is an open challenge. Thus, we develop a surrogate-based optimization approach. Our framework employs the Random Forest algorithm to model the forward relationship between the laser fabrication parameters and the resulting optical characteristics. During the optimization process, we use a greedy, prediction-based exploration strategy that iteratively selects batches of laser parameters to be used in experimentation by minimizing the predicted discrepancy between the surrogate model’s outputs and the user-defined target optical characteristics. This strategy allows for efficient identification of optimal fabrication parameters without the need to model the error landscape directly. We demonstrate the efficiency and effectiveness of our approach on two synthetic benchmarks and two specific experimental applications of photonic surface inverse design targets. By calculating the average performance of our algorithm compared to other state of the art optimization methods, we show that our algorithm performs, on average, twice as well across all benchmarks. Additionally, a warm starting inverse design technique for changed target optical characteristics enhances the performance of the introduced approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A directional electrode separator improves anodic biofilm current density in a well-mixed single-chamber bioelectrochemical system

In this study, a directional electrode separator (DES) was designed and incorporated into a single-chamber bioelectrochemical system (BES) to reduce migration and reoxidation of hydrogen. This issue arises when H 2 , generated at the cathode, travels to the anode where anodic biofilms use H 2 . To test the feasibility of our design, a 3D-printed BES reactor equipped with a DES was inoculated with anaerobic digestor granules and operated under fed-batch conditions using fermented corn stover effluent. The DES equipped reactor achieved significantly higher current densities (~53 A/m²) compared to a conventional single-chamber BES without a separator (~16 A/m²), showing a 3.3 times improvement. Further, control abiotic electrochemical experiments revealed that the DES exhibited significantly higher proton conductivity (456±127 µS/mm) compared to a proton exchange membrane (67±21 µS/mm) with a statistical significance of P=0.03. The DES also effectively reduced H 2 migration to the anode by 21-fold relative to the control. Overall, incorporating a DES in a single-chamber BES enhanced anodic current density by reducing H 2 migration to the anode.

3D printed BES↗
DOE OSTI2022

Investigation of two concrete thermal energy storage system configurations for continuous power production

Two Modelica concrete thermal energy storage (CTES) models are built to analyze potential CTES system designs. One design is the single pipe network design wherein a heat transfer fluid (HTF) flows in one direction during heat deposition and the opposite direction during heat removal. All pipes in the network are used for either deposition (charging) or removal (discharging). The second design is a dual pipe network design in which a network of HTF pipes carries charging fluid and a separate piping network carries discharging HTF. Operation of these two configurations to produce constant steam is evaluated. Here, results indicate that the designs are likely appropriate for different applications with the single pipe network appearing appropriate for batch energy applications and the dual pipe network for continuous energy applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Investigation of the structure, filler interaction and degradation of disulfide elastomers made by Reversible Radical Recombination Polymerization (R3P)

In this paper we report the synthesis and analysis of cyclic polydisulfides, their redox depolymerization as well as their interaction with carbon black. Cyclic poly(3,6-dioxa-1,8-octanedithiol)s (pDODT) were synthesized by Reversible Radical Recombination Polymerization (R3P). R3P is a scalable “green” polymerization process using triethylamine (TEA), H 2 O 2 and air for the polymerization of dithiol monomers. pDODTs were synthesized in 5 and 20 g batches. The effect of peroxide concentration (3, 5, 10, 20, 30%) was investigated. It was found that the molecular weights increased exponentially with increasing peroxide concentrations. At 30% peroxide concentration pDODTs with M n up to 400,000 g/mol and polydispersity of 2 were obtained. 700 and 800 MHz NMR were used for investigating the chemical structure of the polymers. The absence of thiol end group signals in polymers made with > 10wt % H 2 O 2 and M n < 100,000 g/mol, including a sample of a fractionated polymer with starting Mn = 600,000 g/mol, verified cyclic structures. Dithiothreitol reduced the polymers to monomer and a small percentage of oligomers. Carbon black was shown to increase the moduli of high molecular weight pDODTs. Swelling tests in chloroform showed that polymers made with 20 and 30w% peroxide swelled 14 and 10 times of their starting weight without dissolving, indicating chemical interaction between the polymers and carbon black. In conclusion, further analysis of these interesting elastomers, available at the 20 g scale, is in progress.

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

Catalytic monoterpenes conversion over Pd and zeolite to sustainable aviation fuel

Waste terpenes from wood drying are typically processed using a Regenerative Thermal Oxidizer (RTO) and converted to CO2 which is then released to the atmosphere. A novel alternative to this approach is to capture these terpenes using Fluidized Bed Concentrator (FBC) technology developed by CaptisAire which provides a useful waste carbon feedstock for subsequent upgrading. We have developed a catalytic approach where these waste terpenes (comprised of pinenes, limonene, and camphene), are converted to aromatic, and cycloalkane hydrocarbons suitable for use as a Sustainable Aviation Fuel (SAF) blendstock without the need for additional hydrogen. This provides a renewable source for aviation fuel aromatics and increases the potential blend content of renewably sourced alkanes from Hydroprocessed Esters and Fatty Acids (HEFA) and Alcohol to Jet (ATJ) processes. Terpenes were converted in a batch reactor at 260 °C and 120 psi N2 using a combination of Pd/C and H-BEA zeolite. Experimental results infer the role of the zeolite Brønsted acid to facilitate isomerization of pinenes and camphene to limonene, which then undergoes subsequent dehydro-aromatization to p-cymene using Pd/C. Hydrogen release from the dehydro-aromatization process, enables inert conversion of monoterpenes to cycloalkanes without additional hydrogen. Monoterpenes conversion to C10 aromatics (60%) and C10 cycloalkanes (40%) in an inert environment, provides a viable route for SAF blendstock sourced directly from captured waste terpenes. Transition to commercial off-the-shelf catalysts, facilitates scale-up of the terpene conversion process and improves commercial viability. We subsequently created hydrocarbon blends using HEFA or ATJ and converted terpenes, providing a 100% SAF blend with potential to replace traditional Jet-A.

Lawal, Ajibola [ORNL] (ORCID:0000000187266570)↗
DOE OSTI2020

Production of levulinic acid and biocarbon electrode material from corn stover through an integrated biorefinery process

To overcome the inefficient biomass conversion, waste generation, and lack of co-production in biorefineries, an integrated process was proposed for the conversion of corn stover into levulinic acid and biocarbon electrode material. Corn stover was pretreated through hydrothermal process using 0.45 wt% K 2 CO 3 which removed 76 wt% lignin and 85 wt% xylan while preserving 83 wt% glucan. This was followed by acid hydrolysis to produce levulinic acid at varying H 2 SO 4 concentrations and reaction time in a batch reactor at 190 °C. At a reaction time of 5 min in 2 wt% H 2 SO 4 , 35.8 wt% and 30 wt% glucan in raw and pretreated corn stover was converted to levulinic acid, respectively. The residue from acid hydrolysis was converted into biocarbon for supercapacitor electrodes via a two-step thermal activation process which showed a specific capacitance of 120 F g -1 . The proposed integrated biorefinery concept provides multiple value-added products for a greater financial and environmental sustainability.

09 BIOMASS FUELS↗
DOE OSTI2021

Labile Fe(III) supersaturation controls nucleation and properties of product phases from Fe(II)-catalyzed ferrihydrite transformation

We report Fe(II)-catalyzed ferrihydrite (Fh) transformation to more crystalline iron (oxyhydr)oxide phases is a widely occurring geochemical process which has been extensively studied as a function of Fe(II)/Fh ratios at fixed Fh loadings. However, recent isolation of an intermediate Fe(III) species resulting from Fe(II)-Fh contact that facilitates transformation by dissolution/reprecipitation suggests that the kinetics and properties of product phases will instead depend mostly on its rate of accumulation to a critical concentration, consistent with principles in the classical nucleation theory (CNT). This suggests a dependence both on the loading of Fe(II) on the surface, which controls the rate of labile Fe(III) formation, as well as the available volume of solution, which also impacts how fast it can achieve its critical concentration to nucleate product phases. To specifically examine the latter effect, here we studied transformation of 15 mg Fh in 1 mM FeSO 4 solutions at pH 7.2 in batch suspensions of 30 mL, 150 mL, and 450 mL volumes. Time-dependent concentrations of aqueous Fe(II), surface-associated Fe(II), and resulting labile Fe(III) were monitored along with bulk solids characterization as a function of time. Transmission electron microscopy (TEM) was used to visualize the evolution of phases at identical locations on TEM grids. The collective results show that the rates of Fh loss and emergence of product lepidocrocite (Lp) and goethite (Gt) as well as their phase proportions, nucleation mode and morphological properties depend directly on the rate of accumulation of the labile Fe(III) precursor to its critical concentration, which in our experiments was controlled simply by varying the available volume of solution into which it enters. Statistical analyses of TEM image data suggest that while both heterogeneous and homogeneous nucleation occurred in all experiments, the former was increasingly favored at lower Fh/solution ratio due to its lower nucleation barrier being more favorable at attendant lower supersaturations of Fe(III).

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

Mechanisms and kinetics of citrate-promoted dissolution of a uranyl phosphate mineral

The formation of uranyl phosphate precipitates is a remediation strategy to limit the mobility of uranium in contaminated soils. However, exposure to organic ligands, like the plant exudate citrate, can remobilize the uranium. The purpose of this study is to provide a more thorough comprehension of citrate-promoted dissolution of a uranyl phosphate mineral, chernikovite [(H 3 O)(UO 2 )(PO 4 )·3H 2 O], by determining the extent of uranium release from chernikovite at a wide range of citrate concentrations. Here we have quantified the kinetics of dissolution and proposed potential mechanisms of chernikovite dissolution to gain a better understanding of the fate of uranyl phosphate precipitates in the environment. Batch dissolution and continuously stirred tank reactor (CSTR) experiments indicate that increasing citrate concentrations from 0.1 mM to 50 mM increases the concentration of dissolved uranium. However, at citrate concentrations of 10 mM and greater, the effectiveness of the ligand to enhance the dissolution decreases. Density functional theory modeling and Raman spectroscopy indicate a strong interaction between citrate and the uranyl ions at the surface of the uranyl phosphate. The interaction between citrate and uranyl suggests the formation of an alteration layer of uranyl-citrate surface complexes, which could impede uranium release at citrate concentrations of 10 mM or higher. Together, these results indicate that the citrate-promoted dissolution of chernikovite is a fast reaction that is hindered by a combination of surface saturation and a secondary-phase precipitation reaction at higher concentrations of citrate and reaction times greater than or equal to 72 min based on the hydraulic residence time of CSTR experiments. These findings improve our understanding of the stability of uranyl phosphate for use in sequestration of uranium groundwater plumes at contaminated sites.

58 GEOSCIENCES↗
DOE OSTI2023

Influence of clay mineral weathering on green rust formation at iron-reducing conditions

Green rusts (GR) are important drivers for trace metal and nutrient cycling in suboxic environments. We investigated whether green rusts would incorporate aluminum (Al) or other elements from naturally-formed clay minerals containing easily-weatherable clay minerals (e.g. mica, interlayered clays). We isolated the clay minerals from a Matapeake silt loam soil by removal of silt and sand, organic matter, and reducible oxides to study mechanisms of interaction between Fe(II) and soil-sourced clay minerals. We conducted batch Fe(II) sorption experiments at multiple near-neutral pHs (6.5–7.5) and reaction times (2h–365 days). Mineral transformations were characterized by selective extractions, X-ray diffraction (XRD), and Fe X-ray absorption spectroscopy (XAS) analyzed by shell-fitting and linear combination fitting (LCF) with natural and synthetic standards. Clay mineral fraction contained a mixture of quartz, kaolinite, interlayered vermiculite, mica, and chlorite with significant structural Fe (2.6% wt). Uptake of Fe(II) increased with pH and kinetics were rapid until 5 days, followed by slow continuous Fe(II) uptake. Citrate-bicarbonate desorption kinetics from Fe(II) sorbed clay released more Al and silicon (Si) compared with unreacted soil clay fraction whereas magnesium (Mg) and potassium (K) were unaffected. Citrate-bicarbonate extracted Fe contained more Fe(II) than an ideal GR with an Fe(II)/Fe(III) molar ratio of 5.50. Analysis of the Fe EXAFS by both LCF and shell fitting was best modeled as a combination of Fe(III)-clay reduction to Fe(II) and precipitation of GR and Fe(II)-Al LDH. After 7 days of Fe(II) sorption, LCF identified 55.2% total Fe in clay, 33.4% GR(Cl) and 11.4% Fe(II)-Al LDH. Further, these results provide novel evidence of Fe(II)-Al LDHs precipitating on naturally-formed soil clay minerals as a minor phase to GR. The geochemical implications are that GRs formed in soils and sediments should be considered to have Al and Si as well as Mg substitutions affecting their structure and reactivity.

58 GEOSCIENCES↗
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2026

Comparative Dissolution of Iron-Bearing Minerals by Catecholate and Hydroxamate Siderophores under Oxic and Anoxic Conditions

Siderophores play a crucial role in biological iron (Fe) acquisition and mobilization by promoting Fe mineral and rock weathering. While the effects of the hydroxamate siderophore desferrioxamine B (DFOB) have been extensively studied, the role of catecholates in the dissolution of Fe mineral and rock under varying redox conditions remain limited. Moreover, despite Fe being one of the most redox-active metals, the siderophore-mediated redox transformations of solid-phase Fe and their impact on mineral dissolution are not well understood. Herein, this study systematically investigated dissolution behavior of both Fe(II) and Fe(III)-bearing minerals and rocks (Fe(II)- bearing basalt and olivine, and Fe(III)- bearing nontronite and goethite), by two structurally distinct biological siderophores, catecholate protochelin and hydroxamate DFOB, under both oxic and anoxic conditions. Batch experiments quantified Fe and transition metals (Mn, Ni, Zn, Mo) released from the tested minerals in the presence of siderophores. Throughout the dissolution experiments, siderophore concentrations and Fe-siderophore complexation in solutions were measured using UV-vis spectrometry. Fe redox changes under oxic and anoxic conditions by siderophores were determined, and mineral surface alterations following siderophore treatments were characterized through scanning electron microscopy (SEM) and time-of-flight secondary ion mass spectroscopy (ToF-SIMS). Results revealed distinct interplays of dissolution mechanisms, including ligand-exchange promoted dissolution and reductive dissolution, along with Fe(III) reduction and Fe(II) oxidation, contributing to the Fe(II)- and Fe(III)-rich mineral weathering under varying redox conditions. Under oxic conditions, both protochelin and DFOB enhanced Fe release from Fe(II)-rich basalt and olivine more effectively than from Fe(III)-rich minerals. DFOB induced greater Fe(II)-mineral dissolution than protochelin. This difference was attributed to a higher level of Fe(II) oxidation by DFOB in contrast to protochelin, as well as the protochelin oxidation with the loss of Fe binding groups (catechols). Under anoxic conditions, both siderophores significantly reduced Fe(III) in nontronite and goethite, with protochelin demonstrating substantially stronger Fe(III) reduction capacity than DFOB. However, Fe(III) reduction negatively impacted Fe release from Fe(III)-rich minerals. Instead, Fe release from Fe(II)-rich minerals was enhanced under anoxic conditions due to the absence of Fe(II) oxidation and labile Fe(II). Variations in siderophore adsorption also contributed to dissolution differences. Beyond iron, the release of transition metals from tested minerals was influenced by elemental contents, mineral compositions, and siderophore types, revealing distinct metal- and siderophore-dependent patterns. This systematic investigation highlights the roles of siderophores in Fe redox cycling and metal mobilization, enhancing our understanding of different siderophore behaviors in siderophore-mediated microbial metal acquisition within redox-dynamic environments, with implications for bioleaching industries and applications in agriculture and climate change mitigation.

Guo, Dongyi↗
DOE OSTI2026

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↗
DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2021

Temperature-responsive smart tracers for field-measurement of inter-well thermal evolution: Heterogeneous kinetics and field demonstration

Temperature-responsive smart tracers enable advanced warning of a “premature thermal breakthrough” and may therefore improve reservoir management and reduce financial risk. A successful calculation of inter-well fluid temperatures requires that temperature-dependent kinetics are known, which typically result from homogeneous batch reactor experiments. However, recent meso-scale field experiments at the Altona Field Laboratory involving fluid flow in a discretely-fractured reservoir suggest that silica-fluid interfaces may accelerate hydrolysis kinetics. Here, the Arrhenius parameters of phenyl acetate hydrolysis are investigated under heterogeneous reaction conditions in packed-bed column experiments. The breakthrough curve of the reaction product is compared to an inert reference tracer (carbon-cored nanoparticles). Temperatures experienced during field testing ranging from 10 to 40 °C were studied using a phosphate buffer solution to maintain a near-neutral pH, as found in the reservoir. The empirically-determined pre-exponential factor and activation energy values are subsequently used in a first-order kinetic model to improve calculations of effective reservoir temperatures for these meso-scale field tests. Here, the results suggest that hydrolysis reactions can directly monitor inter-well reservoir temperatures, but only if the influence of solid-fluid interactions are carefully accounted for.

15 GEOTHERMAL ENERGY↗