SEARCH · Engineering Papers

Results for “response surface methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Impact of high-dose gamma-ray irradiation on electrical characteristics of N-polar and Ga-polar GaN p–n diodes

We investigate the impact of high-dose gamma-ray irradiation on the electrical performance of Ga-polar and N-polar GaN-based p-n diodes grown by metalorganic chemical vapor deposition. We compare the current density-voltage (J-V), capacitance-voltage (C-V), and circular transfer length method (CTLM) characteristics of the p-n diodes fabricated on Ga-polar and N-polar orientations before and after irradiation. The relative turn-on voltage increases for the Ga-polar diodes with increasing irradiation dose, while it increases initially and then starts to decrease for the N-polar diodes. The p-contact total resistance increases for Ga-polar and decreases for N-polar samples, which we attribute to the formation of point defects and additional Mg activation after irradiation. The J-V characteristics of most of the tested diodes recovered over time, suggesting the changes in the J-V characteristics are temporary and potentially due to metastable occupancy of traps after irradiation. X-ray photoelectron spectroscopy (XPS) and photoluminescence (PL) measurements reveal the existence of different types of initial defects and surface electronic states on Ga-polar and N-polar samples. Gallium vacancies (V Ga ) are dominant defects in Ga-polar samples, while nitrogen vacancies (V N ) are dominant in N-polar samples. The presence of a higher concentration of surface states on Ga-polar surfaces compared to N-polar was confirmed by calculating the band bending and the corresponding screening effect due to opposite polarization bound charge and ionized acceptors at the surface. The difference in surface stoichiometry in these two orientations is responsible for the different behavior in electrical characteristics after gamma-ray interactions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Anode Upcycling via Tailored Solvent Treatment

To achieve a truly closed-loop direct recycling process for lithium-ion batteries, all component materials must be recovered. To date, direct recycling method development has primarily focused on the high-value transition-metal cathode materials, while the inherently lower-value graphite has been challenging to recover in a cost-effective manner. However, end-of-life graphite contains a unique engineered value due to the presence of the solid electrolyte interphase (SEI). Growth of the SEI during the cell's active lifetime stabilizes the electronically reactive graphite surface through an irreversible consumption of Li, and thus necessitates both excess lithiation of the cathode and a costly and time-intensive formation procedure during manufacturing. An optimized pre-formed SEI that capitalizes on existing SEI components from end-of-life batteries has the potential to significantly reduce cathode lithiation requirements and eliminate the critical bottleneck of formation cycling during cell remanufacturing. Further, retaining Li at the anode obviates the need for a separate Li leaching and recovery step, improving the overall efficiency of the direct recycling line. In this work, we present a novel approach to "upcycling" spent graphite through use of tailored chemical treatment to remove adverse (i.e., highly resistive and/or poorly passivating) SEI species while retaining beneficially passivating components. We have explored a rational set of solvents spanning a range of polarity, proticity, and molecular size to evaluate structure-property-performance relationships between applied solvent(s), removed and remaining SEI species, and electrochemical response of the resulting graphite product. Further, we have developed and optimized a robust and holistic analysis procedure that couples symmetric-cell electrochemical testing, multi-modal materials characterization, and advanced electrochemical modeling. These analysis results inform a set of correlative metrics for graphite performance relative to both solvent properties and upcycled SEI composition. We demonstrate effective tunability in the residual SEI composition by varying solvent identity and concentration, and report on several promising solvent systems that achieve comparable or performance to pristine graphite.

anode recycling↗

DOE OSTI2026

Quantitative Infrared-to-Terahertz Nanospectroscopy of Semiconductors

Semiconductor technology now employs few-nanometer features, necessitating tools probing electronic properties on the same length scale. While the concentration of free charge carriers is routinely measured, the scattering rate remains challenging to access at the nanoscale. Here, we present ultrabroadband (5–50 THz) synchrotron infrared nanospectroscopy as a quantitative metrology tool for semiconductors. This technique can determine both the charge carrier concentration and scattering rate with percent-level accuracy, and it is inherently capable of ∼10 nm spatial resolution. We study silicon with different doping levels and confirm the method’s accuracy by statistical analysis and comparison with established far-field infrared spectroscopy. Near-field measurements systematically reveal charge-carrier concentrations ∼30% lower than far-field values, consistent with increased surface sensitivity and surface depletion. Our work establishes synchrotron infrared nanospectroscopy as a precise tool for quantitative nanoscale semiconductor characterization and paves the way toward all-optical characterization of surface depletion effects.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Surface water nitrogen and sediment potential nitrate reduction rates, nutrient stocks, and stable isotopes from nine wetlands at the Tanglewood Biological Station, Alabama

This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation disturbances. Bimonthly surface water and sediment sampling events were conducted at nine wetland sites situated within the Tanglewood Biological Station in Alabama from April 2023 to February 2024. The contents included in the data package include surface water nitrogen (nitrogen oxides and ammonium) and sediment potential nitrate reduction rates (measured as potential denitrification and dissimilatory nitrate reduction to ammonium processing), nutrient stocks (total carbon, total nitrogen, and organic matter), and stable isotopes (carbon and nitrogen). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024) and related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025). In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) sediment potential denitrification rate, (2) sediment potential dissimilatory nitrate reduction to ammonium (DNRA) rate, (3) sediment total carbon and nitrogen content, (4) sediment stable isotopes (delta nitrogen-15 and delta carbon-13), (5) sediment percent organic matter, (6) surface water nitrous oxides, (7) surface water ammonium, and (8) methods codes. All files are .csv or .pdf.

Ammonium↗

DOE OSTI2023

Performance assessment for climate intervention (PACI): preliminary application to a stratospheric aerosol injection scenario

As the prospect of exceeding global temperature targets set forth in the Paris Agreement becomes more likely, methods of climate intervention are increasingly being explored. With this increased interest there is a need for an assessment process to understand the range of impacts across different scenarios against a set of performance goals in order to support policy decisions. The methodology and tools developed for Performance Assessment (PA) for nuclear waste repositories shares many similarities with the needs and requirements for a framework for climate intervention. Using PA, we outline and test an evaluation framework for climate intervention, called Performance Assessment for Climate Intervention (PACI) with a focus on Stratospheric Aerosol Injection (SAI). We define a set of key technical components for the example PACI framework which include identifying performance goals, the extent of the system, and identifying which features, events, and processes are relevant and impactful to calculating model output for the system given the performance goals. Having identified a set of performance goals, the performance of the system, including uncertainty, can then be evaluated against these goals. Using the Geoengineering Large Ensemble (GLENS) scenario, we develop a set of performance goals for monthly temperature, precipitation, drought index, soil water, solar flux, and surface runoff. The assessment assumes that targets may be framed in the context of risk-risk via a risk ratio, or the ratio of the risk of exceeding the performance goal for the SAI scenario against the risk of exceeding the performance goal for the emissions scenario. From regional responses, across multiple climate variables, it is then possible to assess which pathway carries lower risk relative to the goals. The assessment is not comprehensive but rather a demonstration of the evaluation of an SAI scenario. Future work is needed to develop a more complete assessment that would provide additional simulations to cover parametric and aleatory uncertainty and enable a deeper understanding of impacts, informed scenario selection, and allow further refinements to the approach.

climate intervention↗

DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

DOE OSTI2022

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves

Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Soil nitrogen mineralization rates, nutrient stocks, stable isotopes, and water volumetric measurements across terrestrial-aquatic interfaces from three wetlands at the Tanglewood Biological Station, Alabama

This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation events. Soil samples were collected across four sampling events along terrestrial-aquatic gradients at three wetland sites located within the Tanglewood Biological Station in Alabama from April 2024 to June 2025. The contents in this data package include soil in-situ nitrogen mineralization rates (measurements of net nitrification, net ammonification, and net mineralization), nutrient stocks (total carbon, total nitrogen, and organic matter), stable isotopes (carbon and nitrogen), and water volumetric measurements (water-filled pore space). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024), related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025), and related surface water sediment chemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3377325 (Molina Serpas et al., 2026). In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) net nitrification rate, (2) net ammonification rate, (3) areal net mineralization rate, (4) percent organic matter, (5) water-filled pore space, (6) total carbon content, (7) total nitrogen content, (8) stable carbon isotope (delta carbon-13), and (9) stable nitrogen isotope (delta nitrogen-15), and (10) methods codes. All files are .csv or .pdf.

13-C↗

DOE OSTI2022

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2022

An approach to sulfate geoengineering with surface emissions of carbonyl sulfide

Sulfate geoengineering (SG) methods based on lower stratospheric tropical injection of sulfur dioxide (SO 2 ) have been widely discussed in recent years, focusing on the direct and indirect effects they would have on the climate system. Here a potential alternative method is discussed, where sulfur emissions are located at the surface or in the troposphere in the form of carbonyl sulfide (COS) gas. There are two time-dependent chemistry–climate model experiments designed from the years 2021 to 2055, assuming a 40 Tg-Syr -1 artificial global flux of COS, which is geographically distributed following the present-day anthropogenic COS surface emissions (SG-COS-SRF) or a 6 Tg-Syr -1 injection of COS in the tropical upper troposphere (SG-COS-TTL). The budget of COS and sulfur species is discussed, as are the effects of both SG-COS strategies on the stratospheric sulfate aerosol optical depth (~Δτ=0.080 in the years 2046–2055), aerosol effective radius (0.46 µm), surface SO x deposition (+8.9 % for SG-COS-SRF; +3.3 % for SG-COS-TTL), and tropopause radiative forcing (RF; ~-1.5 W m –2 in all-sky conditions in both SG-COS experiments). Indirect effects on ozone, methane and stratospheric water vapour are also considered, along with the COS direct contribution. According to our model results, the resulting net RF is –1.3 W m –2 , for SG-COS-SRF, and –1.5 W m –2 , for SG-COS-TTL, and it is comparable to the corresponding RF of –1.7 W m –2 obtained with a sustained injection of 4 Tg-Syr -1 in the tropical lower stratosphere in the form of SO 2 (SG-SO2, which is able to produce a comparable increase of the sulfate aerosol optical depth). Significant changes in the stratospheric ozone response are found in both SG-COS experiments with respect to SG-SO2 (~5 DU versus +1.4 DU globally). According to the model results, the resulting ultraviolet B (UVB) perturbation at the surface accounts for –4.3 % as a global and annual average (versus –2.4 % in the SG-SO2 case), with a springtime Antarctic decrease of –2.7 % (versus a +5.8 % increase in the SG-SO2 experiment). Overall, we find that an increase in COS emissions may be feasible and produce a more latitudinally uniform forcing without the need for the deployment of stratospheric aircraft. However, our assumption that the rate of COS uptake by soils and plants does not vary with increasing COS concentrations will need to be investigated in future work, and more studies are needed on the prolonged exposure effects to higher COS values in humans and ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

A New Method for High Resolution Surface Change Detection: Data Collection and Validation of Measurements from UAS at the Nevada National Security Site, Nevada, USA

The use of uncrewed aerial systems (UAS) increases the opportunities for detecting surface changes in remote areas and in challenging terrain. Detecting surface topographic changes offers an important constraint for understanding earthquake damage, groundwater depletion, effects of mining, and other events. For these purposes, changes on the order of 5–10 cm are readily detected, but sometimes it is necessary to detect smaller changes. An example is the surface changes that result from underground explosions, which can be as small as 3 cm. Previous studies that described change detection methodologies were generally not aimed at detecting sub-5-cm changes. Additionally, studies focused on high-fidelity accuracy were either computationally modeled or did not fully provide the necessary examples to highlight the usability of these workflows. Detecting changes at this threshold may be critical in certain applications, such as global security research and monitoring for high-consequence natural hazards, including landslides. Here we provide a detailed description of the methodology we used to detect 2–3 cm changes in an important applied research setting—surface changes related to underground explosions. This methodology improves the accuracy of change detection data collection and analysis through the optimization of pre-field planning, surveying, flight operations, and post-processing the collected data, all of which are critical to obtaining the highest output data resolution possible. We applied this methodology to a field study location, collecting 1.4 Tb of images over the course of 30 flights, and location data for 239 ground control points (GCPs). We independently verified changes with orthoimagery, and found that structure-from-motion, software-reported root mean square errors (RMSEs) for both control and check points underestimated the actual error. We found that 3 cm changes are detectable with this methodology, thereby improving our knowledge of a rock’s response to underground explosions.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Scaling carbon fluxes from eddy covariance sites to globe: synthesis and evaluation of the FLUXCOM approach

FLUXNET comprises globally distributed eddy-covariance-based estimates of carbon fluxes between the biosphere and the atmosphere. Since eddy covariance flux towers have a relatively small footprint and are distributed unevenly across the world, upscaling the observations is necessary to obtain global-scale estimates of biosphere–atmosphere exchange. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVMs), here we provide a systematic assessment of the latest upscaling efforts for gross primary production (GPP) and net ecosystem exchange (NEE) of the FLUXCOM initiative, where different machine learning methods, forcing data sets and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent across FLUXCOM and DGVM ensembles (R 2 >0.94 at 1° spatial resolution) while the majority of DGVMs show, for 70 % of the land surface, values outside the FLUXCOM range. Global mean GPP magnitudes for 2008–2010 from FLUXCOM members vary within 106 and 130 PgC yr -1 with the largest uncertainty in the tropics. Seasonal variations in independent SIF estimates agree better with FLUXCOM GPP (mean global pixel-wise R 2 ~0.75) than with GPP from DGVMs (mean global pixel-wise R 2 ~0.6). Seasonal variations in FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R 2 >0.92). Interannual variability of global NEE in FLUXCOM is underestimated compared to inversions and DGVMs. The FLUXCOM version which also uses meteorological inputs shows a strong co-variation in interannual patterns with inversions (R 2 =0.87 for 2001–2010). Mean regional NEE from FLUXCOM shows larger uptake than inversion and DGVM-based estimates, particularly in the tropics with discrepancies of up to several hundred grammes of carbon per square metre per year. These discrepancies can only partly be reconciled by carbon loss pathways that are implicit in inversions but not captured by the flux tower measurements such as carbon emissions from fires and water bodies. We hypothesize that a combination of systematic biases in the underlying eddy covariance data, in particular in tall tropical forests, and a lack of site history effects on NEE in FLUXCOM are likely responsible for the too strong tropical carbon sink estimated by FLUXCOM. Furthermore, as FLUXCOM does not account for CO 2 fertilization effects, carbon flux trends are not realistic. Overall, current FLUXCOM estimates of mean annual and seasonal cycles of GPP as well as seasonal NEE variations provide useful constraints of global carbon cycling, while interannual variability patterns from FLUXCOM are valuable but require cautious interpretation. Exploring the diversity of Earth observation data and of machine learning concepts along with improved quality and quantity of flux tower measurements will facilitate further improvements of the FLUXCOM approach overall.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Real-time characterization of particles produced by laser ablation for analysis by inductively coupled plasma mass spectrometry

Not Available

Suski, Kaitlyn J.↗

DOE OSTI2025

Charge regulation effects on colloidal mixture nanoparticles

Changes in pH within a system containing dissociable sites affect the protonation and deprotonation of these groups, thereby influencing their physical properties. In response, the system modifies their surface charge, affecting electrostatic interactions, aggregation, stability, and structural behavior. Although the pH can be tuned in experiments, it is difficult to model this phenomenon using simulations or theoretical approaches. Here, we perform hybrid Monte Carlo-molecular dynamics simulations to model charge regulation effects in an equimolar colloidal charged system. We compare charge regulation effects with those of a system in which the charges of colloidal nanoparticles are not dissociable. The comparison between the two cases modifies the phase diagram, and it changes the volume fraction where a percolation network of nanoparticles is found. Charge regulation is found to destroy network formation, as the charge in the nanoparticles is modified because of the cooperativity dependency of the degree of charge dissociation sites among the nanoparticles favoring cluster formation. Furthermore, our work suggests that the ionic and/or electronic conductivity in functionalized nanoparticles can be modified by changing pH values. It also guides the experimental design of oppositely charged nanoparticles as inks for 3D printing processes.

Classical statistical mechanics↗

DOE OSTI2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport↗