Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy

Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.

58 GEOSCIENCES↗
DOE OSTI2023

Long-term sinking of nuclear waste canisters in salt formations by low-stress creep at high temperature

Abstract Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of $$-4.4\times 10^{-7}$$ - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus $$-2.1\times 10^{-2}$$ - 2.1 × 10 - 2 mm/year and $$-3.1\times 10^{-2}$$ - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

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

Self-organized criticality in multi-pulse gamma-ray bursts

The variability in multi-pulse gamma-ray bursts (GRBs) may help to reveal the mechanism of underlying processes from the central engine. To investigate whether the self-organized criticality (SOC) phenomena exist in the prompt phase of GRBs, we statistically study the properties of GRBs with more than 3 pulses in each burst by fitting the distributions of several observed physical variables with a Markov Chain Monte Carlo approach, including the isotropic energy E iso , the duration time T, and the peak count rate P of each pulse. Our sample consists of 454 pulses in 93 GRBs observed by the CGRO/BATSE satellite. The best-fitting values and uncertainties for these power-law indices of the differential frequency distributions are: α$^{d}_{E}$ = 1.54±0.09, α$^{d}_{T}$ = 1.82±$^{+0.14}_{0.15}$ and α$^{d}_{P}$ = 2.09$^{+0.18}_{–0.19}$, while the power-law indices in the cumulative frequency distributions are: α$^{c}_{E}$ = 1.44$^{+0.08}_{–0.10}$, α$^{c}_{T}$ = 1.75$^{+0.141}_{–0.13}$ and α$^{c}_{P}$ = 1.99$^{+0.16}_{–0.19}$. We find that these distributions are roughly consistent with the physical framework of a Fractal-Diffusive, Self-Organized Criticality (FD-SOC) system with the spatial dimension S = 3 and the classical diffusion β=1. The results of this study support that the jet responsible for the GRBs should be magnetically dominated and magnetic instabilities (e.g., kink model, or tearing-model instability) lead the GRB emission region into the SOC state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Chapter 2: Global Value Chain and Manufacturing Analysis on Geothermal Power Plant Turbines

The global geothermal power market has shown significant growth since the last decade and is expected to reach a total installed capacity of 18.4 gigawatts electric (GWe) by the end of 2021 (GEA, 2016). The global geothermal power plant turbine market is dominated by a small number of manufacturers. Between 2005 and 2015, 82% of the geothermal steam turbines were manufactured in Japan, and 74% of the geothermal binary cycle turboexpanders were manufactured in Israel. During this period, the United States played an important role in the global trade flow of fully assembled turbine units and turbine parts, with a high volume of imports and exports. Another significant growth area was in Italian turbine/turboexpander manufacturers, who have increased their market share in the last couple of years. One other important change in the manufacturing market was in Turkey, where the bonus on feed-in-tariff (FIT) for domestic hardware components boosted the national manufacturing sector between 2010 and 2020. When planning geothermal power projects, developers customize their power plant size to fit the available geothermal resource capacity. The turbine is designed and sized to optimize the efficiency and utilization of resource and revenue production. The rest of the power plant components such as heat exchangers (HX), water-cooled cooling towers (WCCT), or air-cooled condensers (ACC) are then chosen to complement the turbine size and design. These one-off manufacturing custom design turbines have relatively higher manufacturing set-up costs, longer lead times, and higher capital costs than the standard design turbines manufactured in larger volumes. However, turbines produced in standard increments and in larger manufacturing volumes could result in lower costs per turbine, but potentially lower efficiency. Based on pipeline projects and resource assessments, there is significant potential value in creating standard turbine sizes that could offer an economic advantage, as is done for modular microturbines.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Uranium integral fission product yields for a spectrally-shaped 14.1 MeV neutron source at the National Ignition Facility

This paper describes the experimental results for an energy tuning assembly created to modify the National Ignition Facility deuterium–tritium fusion neutron source into a notional thermonuclear and prompt fission neutron spectrum, which has applications in integral measurements, nuclear data benchmarks, and radiation effects on microelectronics. The Monte Carlo neutron transport utilized MCNP5 to estimate the ETA-modified fluence using the ENDF-B/VIII.0 and IRDFF-II continuous energy nuclear data libraries, and SCALE Sampler was used to estimate the systematic nuclear data covariance using ENDF-B/VII.1 and IRDFF-II in a 252-group structure. The experiment fielded eight activation foils and a highly enriched uranium sample. This provided fifteen reaction channels that were used in a forward-fit comparison to the modeled results and to unfold the neutron spectrum using STAYSL. Gamma-ray spectrometry was performed on the activation and highly enriched uranium foils, and the reduced Χ 2 between the modeled and experimental values was 1.21. The results from the STAYSL unfold, reduced Χ 2 = 1.62, indicated that the modeled neutron spectrum was achieved and the systematic nuclear data uncertainty associated with the neutron transport and activation product cross sections was representative of the experiment. Integral cumulative fission product yield data were collected for 37 mass chains with a combination of gamma-ray spectrometry and radiochemical analysis. Fission product analysis was generally in agreement with two models using a semi-empirical fit and the General Observables of Fission code, with the exception of mass chains 88, 109, 111, 112, 113, 129, 139, 142, 144, 151, and 156.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Cosmology Intertwined III: $f \sigma_8$ and $S_8$

The standard $\Lambda$ Cold Dark Matter cosmological model provides a wonderful fit to current cosmological data, but a few tensions and anomalies became statistically significant with the latest data analyses. While these anomalies could be due to the presence of systematic errors in the experiments, they could also indicate the need for new physics beyond the standard model. In this Letter of Interest we focus on the tension of the Planck data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models for solving this tension, and we discuss the importance of trying to fit with a single model a full array of data and not just one parameter at a time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Origin of the background absorption in carbon nanotubes: Phonon-assisted excitonic continuum

Excitonic effects in 1D semiconductors can be qualitatively different from those in higher dimensions. In particular, the Sommerfeld factor, the ratio of the above-band-edge excitonic continuum absorption to free electron-hole pair generation, has been shown to be less than 1 (suppressed) in 1D systems while it is larger than 1 (enhanced) in 2D and 3D systems. Strong continuum suppression indeed exists in semiconducting single-wall carbon nanotubes, a prototypical 1D semiconductor. However, absorption spectra for carbon nanotubes are typically fit with a combination of Lorentzians and a polynomial background baseline with little physical meaning. We performed absorption measurements in aligned single-chirality (6,5) carbon nanotube films. The obtained spectra were fit with our theoretical model obtained by solving the Boltzmann scattering equation involving fifty-nine different types of transitions among three different types of quasiparticles. Specifically, we took into account microscopic interactions between photons, phonons, and excitons, including their dispersions, which unambiguously demonstrated that the background absorption is due to phonon-assisted transitions from the semiconductor vacuum to finite-momentum continuum states of excitons. The excellent agreement we obtained between experiment and theory suggests that our numerical technique can be seamlessly extended to compute strongly out-of-equilibrium many-body dynamics and time-resolved spectra in low-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Soot-particle core-shell and fractal structures from small-angle X-ray scattering measurements in a flame

We have characterized soot particles measured in situ in a laminar co-flow ethylene-air diffusion flame using small-angle X-ray scattering (SAXS). The analysis includes temperature measurements made with coherent anti-Stokes Raman spectroscopy (CARS) and complements soot volume-fraction and maturity measurements made with laser-induced incandescence (LII). We compared the results of fits to the SAXS measurements using a unified model and a fractal core-shell model. Power-law parameters yielded by the unified model indicate that aggregates of primary particles are in the mass-fractal regime, whereas the primary particles are in the surface-fractal regime in the middle of the flame. Higher and lower in the flame, the primary-particle power-law parameter approaches 4, suggesting smooth primary particles. Furthermore, these trends are consistent with fits using the fractal core-shell model, which indicate that particles have an established core-shell structure in the middle of the flame and are internally homogeneous at higher and lower heights in the flame. Primary-particle size distributions derived using the fractal core-shell model demonstrate excellent agreement with distributions inferred from transmission electron microscopy (TEM) images in the middle of the flame. Higher in the flame, a second small mode appears in the size distributions, suggesting particle fragmentation during oxidation. Surface oxidation would explain (1) aggregate fragmentation and (2) loss of core-shell structure leading to smoother primary-particle surfaces by removal of carbon overlayers. SAXS measurements are much more sensitive to incipient and young soot particles than LII and demonstrate significant volume fraction from particles low in the flame where the LII signal is negligible.

42 ENGINEERING↗
DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

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

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2025

Spectrophotometric determination of the stability of La hydroxyl complexes at near neutral to alkaline pH from 25 to 75 °C

The hydrolysis of rare earth elements (REE) potentially controls their mobility during fluid-rock interaction in a broad range of pH and temperature conditions. However, there is still a lack of thermodynamic data for modeling accurately the stability of REE hydroxyl complexes in hydrothermal aqueous fluids. Here, in this study, UV–Vis spectrophotometric experiments were conducted from 25 to 75 °C in near-neutral to alkaline NaOH-bearing aqueous solutions with varying lanthanum (La) concentrations (0 to ∼0.23 mmol/kg). The color indicator m-cresol purple was used to determine in situ pH and derive the average OH− ligand number ($\overrightarrow{n}$) and formation constants for the La hydroxyl complexes (LaOH 2+ , La(OH) 2 + , and La(OH) 3 0 ). From 25 to 50 °C, $\overrightarrow{n}$ ranges between ∼1 and 2 at pH from 7.0 to 9.3. At 75 °C, $\overrightarrow{n}$ ranges between ∼1.5 and 3 at pH from 6.3 to 8.8. These results suggest the predominance of LaOH 2+ and La(OH) 2 + complexes from 25 to 50 °C, and an increased predominance of La(OH) 3 0 at 75 °C. The cumulative formation constants (β n °, n = 1 to 3) are derived for the reaction La 3+ + nOH − = La(OH) n 3-n , and fitted between 25 and 250 °C by combining the UV–Vis and literature solubility data. The resulting logβ n ° are expressed as function of temperature (T in Kelvin): logβ 1 ° = −1.786 + 0.0133 T + 1.049·10 3 /T; logβ 2 ° = −5.797 + 0.0267 T + 2.713·10 3 /T; logβ 3 ° = 6.435 + 0.0223 T + 512.7/T. A comparison between these new fits and existing extrapolations using the Helgeson-Kirkham-Flowers equation of state indicates significant differences in the predicted hydrolysis of La. The latter extrapolations should therefore be updated for the hydrolysis of REE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗