Search NASA⌕ Search

SEARCH · Search NASA

Results for “EXTRAPOLATION”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A multi-scale calibration approach for process-oriented aggregated building energy demand models

Long-term energy planning relies largely on projections of future energy demand and hourly load profiles. Aggregate building models are increasingly being utilized to characterize the sensitivity of current and future building stocks to changes in climate, population, and building technology on city-to-regional-to-national scales. Due to challenges in the availability of data, those analyses have been limited to projection of energy demand at a single scale, usually state or country, while long-term planning power system models and production cost models might operate at different spatial scales such as energy regions which are focused on ensuring adequate generation infrastructure. In this work, we propose and evaluate a novel method to calibrate an aggregate building energy demand model (PNNL’s BEND model) against the best available data at the spatial scale of balancing authorities. This approach extends previous work on aggregated building energy demand by facilitating analysis of building energy demand across scales, in particular policy and operational decision-making scales. We show that the bias-corrected model estimates building electric loads reasonably well compared with estimates from a statistical model, but has the additional feature of flexibility across spatial scales. While the calibration approach is presently U.S.-centric and associated with U.S. energy regions, it can be extrapolated to other worldwide regions with similar scale challenges between policy and operational implementation decision making. We discuss the significant challenges involved in formulating and calibrating a complex physical model based on simulations of roughly 100,000 individual buildings against available aggregate regional electric load data and highlight areas for potential future work and improved data collection.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Futures for electrochromic windows on high performance houses in arid, cold climates

This study investigates high performance electrochromic windows used on a passive house and residential dwelling to IECC 2021 (i.e., IECC dwelling). In the lab, the electrochromic film switches transmitted solar heat gain coefficient (SHGC) from 0.09 to 0.7 and visible transmittance from 0.15 to 0.82 with power consumption of 1.23 W/m 2 during switching times less than 3 minutes. We extrapolate these results to a window assembly. Building energy models of the houses were evaluated in Santa Fe, New Mexico. A Monte Carlo analysis for 2020, 2040, 2060, and 2080 was conducted for Shared Socioeconomic Pathways 2-4.5, 3-7.0, and 5-8.5. Cases with and without the electrochromic windows and with and without electricity were used to determine energy use intensity and hours beyond thermal safety thresholds. The passive house showed 1.3-3.1% mean energy savings and the IECC dwelling 4.4-5.1% with electrochromic efficiency benefits growing into the future for both cases. Even so, overall savings decrease into the future for the passive house, due to growth in cooling load being dominant, conversely overall energy savings increase into the future for the IECC dwelling due to heating loads being dominant. For thermal resilience, the passive house exhibited a mean percent decrease of 0.02-0.31% hours in the extreme caution (i.e., > 32.2 °C, ≤ 39.4 °C) range while the IECC dwelling exhibited 0.38-4.38%. The study therefore shows that electrochromic windows will have smaller benefits for the passive house in comparison to the IECC dwelling. The relationship between electrochromic windows is shown to have a complex relationship between house efficiency and climate change by these results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Physics-informed hybrid modeling methodology for building infiltration

Infiltration is responsible for one-third to one-half of the space conditioning load of a typical residential home, but the modeling of infiltration for building energy modeling is either represented by over-simplified equations or dependent on over-generalized rules of thumb. Here, this paper develops a physics-informed data-driven methodology for modeling infiltration using building-specific empirical measurements. The developed hybrid methodology combines machine-learning categorization and grey-box sub-modeling to improve the accuracy and generalization of commonly used grey-box infiltration models. The developed methodology excels at predicting infiltration by improving the ability to predict infiltration under unseen environmental conditions using machine learning algorithms with physical significance. In a case study conducted using the iUnit, a modular studio apartment experimental test facility located at the National Renewable Energy Laboratory, we use empirical airtightness measurements to fit an infiltration model using the developed methodology. We find that the developed methodology can improve the overall model accuracy by 43% and improve extrapolation by 38%, compared with the model based on the common grey-box infiltration equation. We also notice that the selected features can improve the performance of a pure machine-learning model, indicating that our methodology identifies the features with the most physical significance to infiltration modeling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Image phase shift invariance based cloud motion displacement vector calculation method for ultra-short-term solar PV power forecasting

Irradiance received on the earth's surface is the main factor that affects the output power of solar PV plants, and is chiefly determined by the cloud distribution seen in a ground-based sky image at the corresponding moment in time. It is the foundation for those linear extrapolation-based ultra-short-term solar PV power forecasting approaches to obtain the cloud distribution in future sky images from the accurate calculation of cloud motion displacement vectors (CMDVs) by using historical sky images. Theoretically, the CMDV can be obtained from the coordinate of the peak pulse calculated from a Fourier phase correlation theory (FPCT) method through the frequency domain information of sky images. The peak pulse is significant and unique only when the cloud deformation between two consecutive sky images is slight enough, which is likely possible for a very short time interval (such as 1?min or shorter) with common changes in the speed of cloud. Sometimes, there will be more than one pulse with similar values when the deformation of the clouds between two consecutive sky images is comparatively obvious under fast changing cloud speeds. This would probably lead to significant errors if the CMDVs were still only obtained from the single coordinate of the peak value pulse. However, the deformation estimation of clouds between two images and its influence on FPCT-based CMDV calculations are terrifically complex and difficult because the motion of clouds is complicated to describe and model. Therefore, to improve the accuracy and reliability under these circumstances in a simple manner, an image-phase-shift-invariance (IPSI) based CMDV calculation method using FPCT is proposed for minute time scale solar power forecasting. First, multiple different CMDVs are calculated from the corresponding consecutive images pairs obtained through different synchronous rotation angles compared to the original images by using the FPCT method. Second, the final CMDV is generated from all of the calculated CMDVs through a centroid iteration strategy based on its density and distance distribution. Third, the influence of different rotation angle resolution on the final CMDV is analyzed as a means of parameter estimation. Simulations under various scenarios including both thick and thin clouds conditions indicated that the proposed IPSI-based CMDV calculation method using FPCT is more accurate and reliable than the original FPCT method, optimal flow (OF) method, and particle image velocimetry (PIV) method.

14 SOLAR ENERGY↗
DOE OSTI2018

Assessment of long-term spatio-temporal radiofrequency electromagnetic field exposure

Highlights: • A year's worth of RF-EMF data is collected in an urban sensor network is analysed. • A full spatio-temporal model of the RF-EMF exposure in the study area is built. • Short-term variations of the RF-EMF exposure are collected in temporal profiles. • Models using data collected at different instants are improved by adjusting for time. As both the environment and telecommunications networks are inherently dynamic, our exposure to environmental radiofrequency (RF) electromagnetic fields (EMF) at an arbitrary location is not at all constant in time. In this study, more than a year's worth of measurement data collected in a fixed low-cost exposimeter network distributed over an urban environment was analysed and used to build, for the first time, a full spatio-temporal surrogate model of outdoor exposure to downlink Global System for Mobile Communications (GSM) and Universal Mobile Telecommunications System (UMTS) signals. Though no global trend was discovered over the measuring period, the difference in measured exposure between two instances could reach up to 42 dB (a factor 12,000 in power density). Furthermore, it was found that, taking into account the hour and day of the measurement, the accuracy of the surrogate model in the area under study was improved by up to 50% compared to models that neglect the daily temporal variability of the RF signals. However, further study is required to assess the extent to which the results obtained in the considered environment can be extrapolated to other geographic locations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2018

Assessing the inner core nucleation paradox with atomic-scale simulations

We investigate the conditions required to freeze liquid iron and iron alloys near the centre of Earth's core. It is usually assumed that inner core growth begins once the ambient core temperature falls below the melting temperature of the iron alloy at Earth's centre; however, additional (under)cooling is required to overcome the energy barrier associated with creating a solid–liquid interface. Predictions based on Classical Nucleation Theory (CNT) have estimated a required undercooling of ~1000 K, which cannot be reconciled with predicted core cooling rates of ~100 K Gyr -1 . This apparent contradiction has been called the ‘inner core nucleation paradox’. Here we address three major uncertainties in the application of CNT to inner core nucleation using atomic-scale simulations. First, we simulate freezing in Fe and Fe-O liquids at core conditions to self-consistently constrain all parameters required by the CNT equations. Second, we test the basic validity of CNT by directly calculating the waiting time to observe freezing events in Fe and Fe-O liquids. Third, we investigate the influence of wave-like forcings applied to the atomic simulations, which have been suggested as a means to significantly reduce the energy barrier. Our results are consistent with CNT in the computationally accessible parameter regime, though error estimates on the waiting time can reach 50% of the measurement at the largest undercooling temperatures. Using CNT to extrapolate to inner core conditions yields estimated undercooling of 730 ± 20 K for the pure iron system and 675 ± 35 K for the Fe-O system. Forcings corresponding to large pressure variations of O(10) GPa reduce these values by ~100 K. While our undercooling estimates are significantly lower than previous estimates they are not low enough to resolve the inner core nucleation paradox.

58 GEOSCIENCES↗
DOE OSTI2022

Melting and phase relations of Fe-Ni-Si determined by a multi-technique approach

Many studies have suggested silicon as a candidate light element for the cores of Earth and Mercury. However, the effect of silicon on the melting temperatures of core materials and thermal profiles of cores is poorly understood, due to disagreements among melt detection techniques, uncertainties in sample pressure evolution during heating, and sparsity of studies investigating the combined effects of nickel and silicon on the phase diagram of iron. Here, in this study we develop a multi-technique approach for measuring the high-pressure melting and solid phase relations of iron alloys and apply it to Fe 0.8 Ni 0.1 Si 0.1 (Fe-11wt%Ni-5.3wt%Si), a composition compatible with recent estimates for the cores of Earth and Mercury. This approach combines results (20-83 GPa) from two atomic-level techniques: synchrotron Mossbauer spectroscopy (SMS) and synchrotron x-ray diffraction (XRD). Melting is independently detected by the loss of the Mossbauer signal, produced exclusively by solid-bound iron nuclei, and the onset of a liquid diffuse x-ray scattering signal. The use of a burst heating and background updating method for quantifying changes in the reference background during heating facilitates the determination of liquid diffuse signal onsets and leads to strong reproducibility and excellent agreement in melting temperatures determined separately by the two techniques. XRD measurements additionally constrain the hcp-fcc phase boundary and in-situ pressure evolution of the samples during heating. We apply our updated thermal pressure model to published SMS melting data on fcc-Fe and fcc-Fe 0.9 Ni 0.1 to precisely evaluate the effect of silicon on melting temperatures. We find that the addition of 10 mol% Si to Fe 0.9 Ni 0.1 reduces melting temperatures by ~250 K at low pressures (<60 GPa) and flattens the hcp-fcc phase boundary. Extrapolating our results, we constrain the location of the hcp-fcc-liquid quasi-triple point at 147±14 GPa and 3140±90 K, which implies a melting temperature reduction of 500 K compared with Fe 0.9 Ni 0.1 . The results demonstrate the advantages of combining complementary experimental techniques in investigations of melting under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state↗
DOE OSTI2021

Signatures of rare earth element distributions in fly ash derived from the combustion of Central Appalachian, Illinois, and Powder River basin coals

We report the distribution of Rare Earth elements (REE) in coal-derived fly ashes can have distinctive patterns when fly ashes are produced from different coals within or between basins, such as the Pennsylvanian Class F fly ashes from the Illinois and Central Appalachian basins. Both the Fire Clay coal and a blend of a number of eastern Kentucky coals show strong Gd peaks and an H-type distribution in the Upper Continental Crust-corrected plots. The Fire Clay coal-derived ash has a higher heavy REE concentration than the blended coal-derived ash. The Illinois Basin-derived fly as has an overall lower REE concentration than the latter ashes. Class C fly ash derived from Powder River Basin coals has, with the exception of an Eu peak, a flatter distribution of REE and an overall L-type or indistinct H- versus L-type distribution. The signatures of the REE in fly ashes may be useful in predicting their behavior in the extraction of the REE; simple extrapolations from the basic concentrations and the predicted extraction percentages for ashes from different basins are not necessarily indicative of the actual distribution of the extracted REE.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Model-based compositional predictions for a differential scanning calorimetry/thermogravimetric analysis-mass spectrometry system used for heat of vaporization measurements

A Differential Scanning Calorimetry/Thermogravimetric Analysis (DSC/TGA) procedure for determination of Heat of Vaporization (HOV) of fuel and surrogate fuel mixtures has been previously utilized for evaluation of oxygenate impacts on spark ignition fuels. This analysis has been further leveraged with mass spectrometry (MS) to elucidate vapor phase compositional changes during evaporation. In this study we focus on a simple system of three compounds (ethanol, n-hexane, and n-octane) to provide a model-based understanding of the evaporation process under these experimental conditions. An Aspen Plus model was set up to provide liquid and vapor phase concentration predictions during mixture evaporation using the initial mixture composition, and measured temperature and mass loss profiles from DSC/TGA experiments. Agreements between the trends of the MS ion count profiles and model composition predictions provide a degree of validation of this modeling approach that can allow the interpretation and extrapolation of DSC/TGA measurements for more complex fuel mixtures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2013

Systematic variations of argon diffusion in feldspars and implications for thermochronometry

Coupled information about the time-dependent production and temperature-dependent diffusion of radiogenic argon in feldspars can be used to constrain the thermal evolution attending a host of Earth and planetary processes. To better assess the accuracy of thermal models, an understanding of the mechanisms and pathways by which argon diffuses in feldspars is desirable. Here we present step-heating Ar diffusion experiments conducted on feldspars with diverse compositions, structural states, and microstructural characteristics. The experiments reveal systematic variations in diffusive behavior that appear closely related to these variables, with apparent closure temperatures for 0.1–1 mm grains of ~200–400 °C (assuming a 10 °C/Ma cooling rate). Given such variability, there is no broadly applicable set of diffusion parameters that can be utilized in feldspar thermal modeling; sample-specific data are required. Diffusion experiments conducted on oriented cleavage flakes do not reveal directionally-dependent diffusive anisotropy to within the resolution limits of our approach (approximately a factor of 2). Additional experiments aimed at constraining the physical significance of the diffusion domain are presented and indicate that unaltered feldspar crystals with or without coherent exsolution lamellae diffuse at the grain-scale, whereas feldspars containing hydrothermal alteration and/or incoherent sub-grain intergrowths do not. Arrhenius plots for argon diffusion in plagioclase and alkali feldspars appear to reflect a confluence of intrinsic diffusion kinetics and structural transitions that occur during incremental heating experiments. These structural transitions, along with sub-grain domain size variations, cause deviations from linearity (i.e., upward and downward curvature) on Arrhenius plots. An atomistic model for Arrhenius behavior is proposed that incorporates the variable lattice deformations of different feldspars in response to heating and compression. Furthermore, the resulting implications for accurately extrapolating laboratory-derived diffusion parameters to natural settings and over geologic time are discussed. We find that considerable inaccuracies may exist in published thermal histories obtained using multiple diffusion domain (MDD) models fit to Arrhenius plots for exsolved alkali feldspar, where the inferred Ar partial retention zones may be spuriously hot.

58 GEOSCIENCES↗
DOE OSTI2018

Catalytic generation of methane at 60–100 °C and 0.1–300 MPa from source rocks containing kerogen Types I, II, and III

Low temperature (60 and 100 °C) and long-term (6 months to 5 years) heating of pre-evacuated and sterilized shales and coals containing kerogen Types I (Mahogany Shale), II (Mowry Shale and New Albany Shale), and III (Springfield Coal and Wilcox Lignite) with low initial maturities (vitrinite reflectance Ro 0.39–0.62%) demonstrates that catalytically generated hydrocarbons may explain the occurrence of some non-biogenic natural gas accumulations where insufficient thermal maturity contradicts the conventional thermal cracking paradigm. Extrapolation of the observed rate of catalytic methanogenesis in the laboratory suggests that significant amounts of sedimentary organic carbon can be converted to relatively dry natural gas over tens of thousands of years in sedimentary basins at temperatures as low as 60 °C. Here, our laboratory experiments utilized source rock (shale and coal) chips sealed in gold and Pyrex® glass tubes in the presence of hydrogen-isotopically contrasting waters. Parallel heating experiments applied hydrostatic pressures from 0.1 to 300 MPa. Control experiments constrained the influence of pre-existing and residual methane in closed pores of rock chips that was unrelated to newly generated methane. This study’s experimental methane yields at 60 and 100 °C are 5–11 orders of magnitude higher than the theoretically predicted yields from kinetic models of thermogenic methane generation, which strongly suggests a contribution of catalytic methanogenesis. Higher temperature, longer heating time, and lower hydrostatic pressure enhanced catalytic methanogenesis. No clear relationships were observed between kerogen type or total organic carbon content and methane yields via catalysis. Catalytic methanogenesis was strongest in Mowry Shale where methane yields at 60 °C amounted to ~2.5 μmol per gram of organic carbon after one year of hydrous heating at ambient pressure. In stark contrast to the earlier findings of hydrogen isotopic exchange between water and thermogenic methane in hydrous pyrolysis experiments above 300 °C, the hydrogen isotopic composition of added water exerted limited influence on the δ2H value of methane generated catalytically at low temperatures. We hypothesize that the catalytic sites responsible for methanogenesis are located in hydrophobic microenvironments with limited access to water. The δ13CCH4 values of methane generated catalytically at 60–100 °C range from ~–57.6 to –41.4‰ and are thus similar to typical thermogenic methane (δ13CCH4 >–50‰) and microbially generated methane (<–55‰). Future studies need to evaluate the possibility that clumped isotope characteristics of catalytically generated methane can diagnose the low-temperature regime of catalytic methanogenesis. Furthermore, testing of freshly cored anoxic rocks is needed to determine whether the use of archived, oxygen-exposed rocks in geochemical maturation/catalysis studies introduces artifacts in experimental hydrocarbon yields.

58 GEOSCIENCES↗
DOE OSTI2018

An experimental study of the solubility and speciation of thorium in chloride-bearing aqueous solutions at temperatures up to 250 °C

In this study, the solubility and speciation of Th in Cl-bearing solutions saturated in ThO2 was investigated experimentally at 150, 175, 200, 225 and 250°C, and saturated water vapor pressure. At these experimental conditions (pH 25°C < 3.8, mNaCl < 2.5), Th-chloride complexation is insignificant, and the solubility of Th is controlled by hydrolysis. The two main species in solution were identified as Th(OH) 2 2+ and Th(OH) 4 ° and formed via the reactions Th 4+ + 2OH - ↔ Th(OH)$2+\atop{2}$ and Th 4+ + 4OH - ↔ Th(OH)$0\atop{4}$. The logarithms of the formation constants determined for the first reaction are 23.91 ± 0.05, 24.79 ± 0.07, 25.78 ± 0.05, and 25.85 ± 0.05 for 175, 200, 225 and 250°C, respectively, and 42.18 ± 0.06, 43.50 ± 0.07, 43.17 ± 0.08, and 43.74 ± 0.11 for the second reaction at the same temperatures. Formation constants could not be determined reliably for 150 °C. Extrapolations of our formation constants to 25°C are in excellent agreement with the values determined at 25 °C by Neck and Kim (2001) and Ekberg et al. (2000). Finally, for natural systems buffered by thorianite (ThO 2 ), the solubility determined in our study likely reflects a minimum, as the presence of additional ligands would lead to higher solubility.

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

Uranyl speciation in sulfate-bearing hydrothermal solutions up to 250 °C

The speciation of uranyl in sulfate-bearing solutions at temperatures up to 250 ºC has been investigated using in situ UV-Visible spectrophotometry. Formation constants for the reactions UO$^{2+}_{2}$ + SO$^{2–}_{4}$ $\Longleftrightarrow$ UO 2 SO$^{0}_{4}$ (logβ 1 = 3.38, 4.40, 5.44, 6.33, 7.74) and UO$^{2+}_{2}$ + 2SO$^{2–}_{4}$ $\Longleftrightarrow$ UO 2 (SO 4 )$^{2–}_{2}$ (logβ 2 = 4.10, 5.26, 6.83, 8.00, 9.70) derived at 25, 100, 150, 200 and 250 ºC respectively. These formation constants were fitted to the modified Ryzhenko – Bryzgalin (MRB) model to derive pK(298), A(zz/a) and B(zz/a) values for both complexes: 3.262, 2.212, -197.96 (UO 2 SO$^{0}_{4}$) and 4.189, 3.21, -473.14 UO 2 (SO 4 )$^{2–}_{2}$ respectively. Compared with values extrapolated using previously available data for 25, 70 and 75 ºC, our new data suggest higher stability of UO 2 SO$^{0}_{4}$ at temperatures above 150 ºC and significantly lower stability of UO 2 (SO 4 )$^{2–}_{2}$ at all temperatures above 25Cº. Molar absorbances of both sulfate species were also derived. At 25 Cº we found our molar absorbance for UO 2 SO$^{0}_{4}$ agreed well with those reported previously in the literature, however we report lower peak amplitudes for UO 2 (SO 4 )$^{2–}_{2}$ . We also noted significant temperature-dependent red shifts in the molar absorbance of UO 2 SO$^{0}_{4}$. We imply that these shifts could be explained by changes in sulfate bonding behaviour around the uranyl O04 ion - namely shifts in the distribution of complexes with sulfate bound in either a monodentate or bidentate configuration. Simple models calculated with our new data suggest that sulfate complexes may readily predominate over chloride complexes even in solutions containing upwards of 20 wt% chloride and only 100 ppm sulfate. In addition, moderate concentrations of sulfate (~1 wt%) can increase the solubility of uranium by an order of magnitude. Therefore, we imply that sulfate could play a vital role in mobilising uranium in hydrothermal systems, and that the removal of sulfate via precipitation of sulfate minerals may act as a means of depositing uranium under oxidizing conditions. This in turn might explain the associations between sulfate minerals and uranium mineralisation occasionally seen in large uranium ore deposits such as Olympic Dam and MacArthur River.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗