Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic analyses”

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

Agreement of Analytical and Simulation-Based Estimates of the Required Land Depth in Climate Models

Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ignition delay time and speciation of dibutyl ether at high pressures

Dibutyl ether (DBE) is a promising biofuel due to its high cetane number (~ 100) and high volumetric energy density (31.6 MJ/L). It could either be used directly in compression ignition engines or blended with other conventional or renewable fuels. Oxidation and pyrolysis kinetics of DBE are not well known, particularly at high pressures. In this work, we have experimentally investigated the chemical kinetics of DBE in three domains: (a) ignition delay time measurements in a rapid compression machine over T = 550–650 K, P = 10, 20, 40 bar, φ = 0.5, 1; (b) ignition delay time measurements in a shock tube over T = 900–1300 K, P = 20, 40 bar, φ = 0.5, 1; (c) laser-based carbon monoxide speciation measurements in a shock tube during DBE pyrolysis and oxidation over T = 1100–1400 K, P = 20 bar. Pressure time-histories measured in RCM experiments exhibited unique 3-stage and 4-stage ignition behavior predominantly at fuel-lean conditions. Experimental data were compared with the predictions of two recent chemical kinetic models of DBE. Sensitivity analyses were carried out to identify key reactions which may have caused the discrepancy between experiments and simulations. We found that the rate of decomposition of DBE may need to be revisited to improve the oxidative and pyrolytic predictions of DBE kinetic model.

09 BIOMASS FUELS↗
DOE OSTI2021

Ternary SiGeSn alloy nanocrystals via nonthermal plasma synthesis

We report on the synthesis of ternary SiGeSn nanocrystals (NCs) produced via nonthermal plasma synthesis from silane (SiH 4 ), germane (GeH 4 ), and tetramethylstannane (Sn(CH 3 ) 4 ) precursor sources. Detailed structural, chemical, and vibrational analyses show that all three elements are incorporated both on the NC surface and within the NC core. Incorporation of Sn into the NC core is realized using a secondary injection of SiH 4 and GeH 4 precursor gases in the after-glow region of the plasma, which kinetically traps Sn in the core. We demonstrate compositional tunability of the SiGeSn NCs in which the Si and Ge ratios can be varied broadly at low Sn concentrations. We also show tunability of the Sn content up to ~2 atomic percent as revealed by ICP-MS analysis. More generally, this report demonstrates how nonthermal plasma synthesis can be used to produce metastable ternary nanostructured alloys involving thermodynamically insoluble constituents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The temperatures of ejecta transporting in vacuum and gases

In this work, we measure continuous thermal radiance from evolving clouds of liquid metal fragments ejected into vacuum, nonreactive, and reactive gas. We implement a model for the thermalization of the ejecta and gas and use this to constrain the absolute temperature of the ejecta cloud. This model enables further analyses of ejecta thermal behavior under a variety of conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The FastEddy® Resident-GPU Accelerated Large-Eddy Simulation Framework: Moist Dynamics Extension, Validation and Sensitivities of Modeling Non-Precipitating Shallow Cumulus Clouds

Herein we describe the moist dynamics formulation implemented within the graphics processing unit-resident large-eddy simulation FastEddy® model, which includes a simple saturation adjustment scheme for condensation and evaporation processes. Two LES model intercomparison exercises for non-precipitating shallow cumulus clouds are simulated in order to validate this model extension, including a static forcing and a time-dependent forcing case. Overall, we find our dynamical, thermodynamical and microphysical quantities, along with turbulence variability and fluxes, to be commensurate with the corresponding model intercomparison results. In addition, sensitivities to specific model settings are investigated. Among these settings, it is shown that boundary layer and cloud layer structure and characteristics are sensitive to use of higher-order advection schemes impacting the vertical distribution of cloud content and associated turbulence statistics. Increasing the timescale of the saturation scheme leads to enhanced liquid water presence and decreases vertical velocity variance within the cloud deck. In some cases, these sensitivities agree with the model-to-model variability reported in the intercomparison exercises, highlighting the important role of specific model implementation choices in the context of shallow cumulus convection simulations. These analyses and findings also provide the basis for future extensions and applications of FastEddy® for modeling moist convection and precipitation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Revealing the role of redox reaction selectivity and mass transfer in current–voltage predictions for ensembles of photocatalysts

Photocatalysts are conceptually simple reaction units where nanoscale semiconductors integrated with catalysts drive a pair of redox reactions on illumination. However, the proximity of reaction sites performing cathodic and anodic reactions poses dire challenges to realize large light-to-fuel conversion efficiencies. In this study, a powerful, yet straightforward, equivalent-circuit detail-balance modeling framework is developed and applied to evaluate the performance of photocatalytic systems featuring multiple light absorbers. Specifically, low bandgap iridium-doped strontium titanate is modeled as a Z-scheme photocatalyst to achieve desirable hydrogen evolution and iron-based redox shuttle oxidation reactions. Our model has unique capabilities to simulate competing redox reactions and address mass-transfer limitations. In a significant departure from state-of-the-art circuit models, our study develops tools to perform load-line analyses by incorporating a net electrochemical load curve that includes both desired and competing redox reactions. Consequently, reaction selectivity is predicted from equivalent circuit models for photocatalytic and photoelectrochemical systems. Our investigation into ensembles comprised of multiple, semi-transparent light absorbers reveals their potential to outperform a single, optically thick light absorber, particularly when operated under mass-transfer-limited conditions. However, this outcome hinges on minimizing mass-transfer rates of select redox species to prevent undesired reactions of hydrogen oxidation and/or redox shuttle reduction. Our findings demonstrate that reaction selectivity can be achieved by tuning asymmetry in redox species mass-transfer even with perfectly symmetric electrocatalytic charge-transfer coefficients. The influences of various kinetic, mass-transfer, and thermodynamic parameters are explored to offer crucial insights for synthesis of the next-generation of photocatalysts and selective coatings, and reactor designs.

25 ENERGY STORAGE↗
DOE OSTI2024

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

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

Harnessing dimethyl ether with ultra-low-grade heat for scaling-resistant brine concentration and fractional crystallization

Solvent-driven separations may enable scalable concentration of hypersaline brines, supporting a circular resource economy from the extraction of lithium and rare earth elements from spent battery and magnet leachates. This work analyses a novel solvent-driven water extraction (SDWE) system employing dimethyl ether (DME) and ultra-low-grade heat for brine concentration and fractional crystallization. SDWE exploits DME’s unique properties: (1) a low dielectric constant that promotes water solubility over charged solutes by a factor of 10 3 , and (2) a high volatility that facilitate efficient DME reconcentration with ultra-low-grade heat. The techno-economic viability of SDWE is assessed with a computational framework that encompasses a liquid–liquid separator and a solvent concentrator. We integrate the extended universal quasichemical model with the virial equation of state to predict the compositions of the complex three-phase DME-water mixture at vapor–liquid and liquid–liquid equilibrium. Subsequently, we optimize the thermodynamic and economic performance of SDWE, by controlling the interstage flash pressure, heat source temperature, and the number of concentrating stages. DME-based SDWE concentrates an input saline feed to 5.5 M and regenerates over 99 % of the DME using ultra-low-grade heat below 50 °C, with a DME/water selectivity ratio of 125. Here our calculations reveal that optimal performance is achieved at interstage flash pressures of 0.4 – 0.5 bar for heat source temperatures between 323–373 K, with improved exergetic efficiencies at lower temperatures. At a heat source temperature of 323 K and an interstage pressure of 0.489 bar, DME-driven SDWE achieves an optimal thermodynamic efficiency of 20.5 % and a projected specific cost of US$ 1.93 m -3 . These specific costs suggest that SDWE is competitive with commercialized thermal distillation technologies, while mitigating the traditional risks associated with scaling in heat and mass exchangers with hypersaline brines.

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

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗
DOE OSTI2025

Role of Histidine‐Containing Peptoids in Accelerating the Kinetics of Calcite Growth

Carbonate mineralization, the conversion of CO 2 into stable, thermodynamically favorable carbonate minerals, offers a promising strategy for permanent and environmentally friendly carbon storage, with minimal risk of long-term leakage and minimal monitoring requirements. Drawing inspiration from carbonic anhydrase (CA), a family of zinc-containing metalloenzymes that catalyze the hydration of CO 2 to bicarbonate and promote carbonate precipitation, a class of histidine-containing peptoids was designed that is capable of coordinating with Zn 2+ ions to act as CA mimetics for accelerating calcite step growth. In situ atomic force microscopy (AFM) measurements reveal that these peptoids significantly enhance step advancement, with a more pronounced effect observed when combined with Zn 2+ ions and under higher calcium-to-carbonate activity ratios, indicating that peptoids facilitate the incorporation of CO 3 2− ions at step edges. Solution NMR and 3D atomic force microscopy (3D AFM) analyses show that the coordination of peptoids with Zn 2+ promotes both the deprotonation of HCO 3 − to CO 3 2− and restructures the interfacial hydration layers of calcite, collectively lowering the activation barrier for step growth. These findings establish a design framework for sequence-defined polymers to regulate carbonate mineralization, offering promising applications in CO 2 capture and long-term storage.

CO2 mineralization↗
DOE OSTI2022

Size-Dependent Characteristics of Surface-Rooted Three-Dimensional Convective Objects in Continental Shallow Cumulus Simulations

A segmentation algorithm is applied to high resolution simulations of shallow continental convection to identify individual convective 3D objects within the convective boundary layer, in order to investigate which properties of the objects vary with the object width. The study analyses the geometry of the objects, along with their profiles of vertical velocity and total water, to assess various assumptions often used in spectral mass-flux convection schemes. The methodology of this paper is unique in that we use (a) a novel application of the watershed algorithm to detect individual objects in the constantly evolving continental boundary layer efficiently, and (b) an unprecedentedly large number of simulations being analyzed. In total, 26 days of LASSO simulations at the Atmospheric Radiation Measurement Southern Great Plains site are analyzed, yielding roughly one million objects. Plume-like surface-rooted objects are found to be omnipresent, the vertical extent of which is strongly dependent on the object width. The vertical velocity and moisture anomaly profiles of the widest objects are roughly consistent with the classic buoyancy-driven rising plume model. The kinematic and thermodynamic properties of the objects vary with object width. This width dependence is strongest above cloud base, but much weaker below. Finally the impact of neglecting the contribution of covariances to the vertical moisture flux, which is commonly used in mass-flux parameterizations, is investigated. The average effect of neglecting covariances increases linearly with object width, leading to a 20% flux underestimation for 2 km wide objects. Implications of the results for spectral convection scheme development are briefly discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Group additivity values for entropy and heat capacities of C 2 –C 8 alkanes, alkyl hydroperoxides, and their radicals

Group additivity values for the thermodynamic properties of oxygenated radicals are poorly determined due to the absence of high quality reference data. Here, a set of 58 group additive values (GAV) for the standard entropy and heat capacity of relevance to alkane oxidation is derived from fits to an extensive and accurate database of standard entropies (298.15 K) and heat capacities (300–3000 K) recently calculated with the “STAR-1D” formalism. The 192 species in this database represent the alkanes (RH), alkyl radicals ($\dot{R}$), alkyl hydroperoxides (RO 2 H), alkyl-peroxy (R$\dot{O}$ 2 ) and hydroperoxy-alkyl ($\dot{Q}$OOH) radicals for all of the isomers of C 2 –C 5 alkane fuels and a select number of isomers of C 6 –C 9 species. The STAR-1D thermochemical data for this set of species was previously obtained from a coupling of scaled B2PLYPD3/cc-pVTZ vibrational analyses with scaled ωB97X-D/cc-pVTZ one-dimensional hindered rotor corrections. The 2σ uncertainties in the GAV results relative to the STAR-1D data set are 2.4 cal K –1 mol –1 for the entropies and at most 2.0 cal K –1 mol –1 for the heat capacities in the temperature range 500–800 K. The 2σ fitting uncertainties in the heat capacity gradually reduce at higher temperatures reaching a value of 0.8 cal K –1 mol –1 at 2000 K decreasing to only 0.5 cal K –1 mol –1 at 3000 K. The high degree of accuracy for the GAV representations is obtained through the introduction of various new group terms, together with the re-optimization of existing group terms. Among the full set of 58 GAV terms, 25 include non-next-nearest neighbor interactions (NNI) and β-corrections. The updated GAVs can be applied in the prediction of entropies and heat capacities for a wide range of hydrocarbons and hydroperoxide species and their radicals, which is important to the accurate prediction of fuel reactivity at low-temperatures in the range 600–1000 K.

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

Computation of the Thermal Expansion Coefficient of Graphene with Gaussian Approximation Potentials

Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this work, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

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

Sulfur speciation in dacitic melts using X-ray absorption near-edge structure spectroscopy of the S K -edge (S-XANES): Consideration of radiation-induced changes and the implications for sulfur in natural arc systems

Abstract The synchrotron technique of micro X-ray absorption near-edge structure spectroscopy at the sulfur K-edge (S-XANES) provides a unique opportunity to measure the proportion of different oxidation states of sulfur (S) in silicate glasses. Although applied extensively in the analysis of basaltic silicate glasses, few S-XANES studies have investigated variations in S oxidation states with fO2 in felsic silicate glasses. In addition, no study has systematically compared the S-XANES results obtained from the same samples at different photon flux densities to quantify the relationship between exposure time and changes in S speciation in silicate glass, as has been done for Fe and V. This study evaluates observed differences in S speciation measured in experimentally produced H2O-saturated dacitic glasses over a range of reducing to oxidizing conditions (from log fO2 = ΔFMQ-0.7 to ΔFMQ+3.3; FMQ is the fayalite-magnetite-quartz mineral redox buffer) and equilibrated at 1000 °C and 300 MPa. S-XANES spectra were collected at three different photon flux densities using three microspectroscopy beamlines. As is observed in S-XANES analyses of basaltic silicate glasses, beam-induced changes to the S6+/ΣS are observed as a function of photon flux density and beam exposure time. Our results demonstrate that silicate glasses of dacitic composition undergo beam-induced photo-reduction in samples equilibrated at ΔFMQ > +1.75 and photo-oxidation if equilibrated at ΔFMQ < +1. The time required to observe beam-induced changes in the spectra varies as a function of flux density, and our study establishes an upper photon density limit at ~1.0 × 1012 photons/µm2. The S6+/ΣS calculated from spectra collected below this absorbed photon limit at intermediate flux densities (~1–4 × 109 photons/s per µm2) are affected by beam damage, as no conditions were found to be completely free of beam-induced changes. However, the S6+/ΣS ratios calculated below the limit at intermediate flux densities are consistent with thermodynamic constraints, demonstrating that S6+/ΣS ratios calculated from S-XANES spectra can be considered reliable for estimating the oxygen fugacity. Our results carry important implications for the S budget of felsic magmas and dissolution mechanisms in evolved melts. While our results from all three flux densities show the presence of S4+ dissolved in relatively oxidized (ΔFMQ > +1.75) dacitic glass, even in the spectra exposed to the lowest photon densities, we are unable to rule out the possibility that the S4+ signal is the result of instantaneous X-ray irradiation induced beam damage using S-XANES alone. When our spectra are compared to S-XANES spectra from basaltic silicate glasses, important differences exist in the solubility of S2– and S6+ between dacitic silicate melts, pointing to differences in solubility mechanisms as melt composition changes. This study highlights the need for further investigation into beam damage systematics, presence of S4+, and the solubility mechanisms of different oxidation states of S as silicate melt composition changes.

Geochemistry & Geophysics↗