Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aggregation and fractionation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Mercury Removal from Contaminated Water by Wood-Based Biochar Depends on Natural Organic Matter and Ionic Composition

Biochars can remove potentially toxic elements, such as inorganic mercury [Hg(II)] from contaminated waters. However, their performance in complex water matrices is rarely investigated, and the combined roles of natural organic matter (NOM) and ionic composition in the removal of Hg(II) by biochar remain unclear. Here, we investigate the influence of NOM and major ions such as chloride (Cl – ), nitrate (NO 3 – ), calcium (Ca 2+ ), and sodium (Na + ) on Hg(II) removal by a wood-based biochar (SWP700). Multiple sorption sites containing sulfur (S) were located within the porous SWP700. In the absence of NOM, Hg(II) removal was driven by these sites. Ca 2+ bridging was important in enhancing removal of negatively charged Hg(II)-chloro complexes. In the presence of NOM, formation of soluble Hg-NOM complexes (as seen from speciation calculations), which have limited access to biochar pores, suppressed Hg(II) removal, but Cl – and Ca 2+ could still facilitate it. The ability of Ca 2+ to aggregate NOM, including Hg-NOM complexes, promoted Hg(II) removal from the dissolved fraction (<0.45 μm). Hg(II) removal in the presence of Cl – followed a stepwise mechanism. Weakly bound oxygen functional groups in NOM were outcompeted by Cl – , forming smaller-sized Hg(II)-chloro complexes, which could access additional intraparticle sorption sites. Therein, Cl – was outcompeted by S, which finally immobilized Hg(II) in SWP700 as confirmed by extended X-ray absorption fine structure spectroscopy. We conclude that in NOM containing oxic waters, with relatively high molar ratios of Cl – : NOM and Ca 2+ : NOM, Hg(II) removal can still be effective with SWP700.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of a Process to Separate Ensiled Corn Grain and Stover for Improved Utilization as a Biomass Feedstock

Corn grain and chopped stover were stored anaerobically at 37% to 50% (w.b.) aggregate moisture content to create a unique biomass feedstock. However, after storage, the two fractions must be separated to accommodate different conversion pathways. In this work, a modified cleaning system from a grain combine harvester was evaluated to investigate its effectiveness in separating the grain from the stover. While this system has been purposefully built for this task, there is no previous literature on separating grain at high material other than grain (MOG) to grain ratios. Using this system, material was separated into four fractions: Grain, Tailings, Heavy MOG, and Light MOG. Subsamples were collected, oven dried, and then hydrodynamically separated to quantify the grain content in the four fractions. Several different configurations of air velocities, mass flow rate, material MOG: grain ratios, and sieve types and openings were investigated. Grain capture was defined as the fraction of the total grain dry mass collected in the grain and tailings fractions. In five experiments, the grain capture effectiveness varied from 89.2% to 99.4% on a dry basis. The dry basis foreign matter (FM) in the clean grain was 0.6% to 10.8%. As a result, the use of a modified combine cleaning system has the potential to be one step in a system to fractionate corn grain from stover (i.e., MOG) in a biorefinery that uses both starch and cellulose as biomass feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2022

An aggregation-based nonlinear multigrid solver for two-phase flow and transport in porous media

A nonlinear multigrid solver for two-phase flow and transport in a mixed fractional-flow velocity-pressure-saturation formulation is proposed. The solver, which is under the framework of the full approximation scheme (FAS), extends our previous work on nonlinear multigrid for heterogeneous diffusion problems. The coarse spaces in the multigrid hierarchy are constructed by first aggregating degrees of freedom, and then solving some local flow problems. The mixed formulation and the choice of coarse spaces allow us to assemble the coarse problems without visiting finer levels during the solving phase, which is crucial for the scalability of multigrid methods. Specifically, a natural generalization of the upwind flux can be evaluated directly on coarse levels using the precomputed coarse flux basis vectors. We find that the resulting solver is applicable to problems discretized on general unstructured grids. The performance of the proposed nonlinear multigrid solver in comparison with the standard single level Newton's method is demonstrated through challenging numerical examples. It is observed that the proposed solver is robust for highly nonlinear problems and clearly outperforms Newton's method in the case of high Courant-Friedrichs-Lewy (CFL) numbers.

58 GEOSCIENCES↗
DOE OSTI2021

Effect of stereochemistry on nanoscale assembly of ABA triblock copolymers with crystallizable blocks

We report on the nanoscale assembly of poly(lactic acid)-b-poly(ethylene oxide)-b-poly(lactic acid) (PLA-PEO-PLA) triblock copolymers in water, focusing on the effect of stereochemistry, where the PLA blocks are statistical copolymers of l-lactide and d-lactide with l/d ratios of 100/0, 95/5, 90/10, 85/15, 75/25, and 50/50. Small-angle neutron scattering (SANS) shows a nearly constant d-spacing as concentration varies in triblock systems with l/d ratios of 90/10 and 95/5, which we attribute to inhomogeneity in the structure of these gels, supported by previous USANS and confocal microscopy studies. The SANS data fit well to a core-shell ellipsoid form factor model with a hard-sphere structure factor. Polymeric micelles with l/d ratios from 75/25 to 85/15 displayed very high aggregation numbers, consistent with a strong interaction between PLA chains and the enhanced storage modulus observed in rheological studies of these systems. While the 90/10 and 95/5 samples showed lower aggregation numbers, their SANS profiles shows close spacing between micelles, which may promote a high fraction of intermicellar bridging chains, also consistent with a higher storage modulus. Overall, these results provide insight into the micellar assembly behavior of block copolymers with a crystallizable block, and indicate that tuning stereochemistry of PLA-based block copolymers is an effective means of modifying micellar properties for specific applications.

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

Evidence linking calcium to increased organo-mineral association in soils

Geochemical indicators are emerging as important predictors of soil organic carbon (SOC) dynamics, but evidence concerning the role of calcium (Ca) is scarce. This study investigates the role of Ca prevalence in SOC accumulation by comparing otherwise similar sites with (CaCO 3 -bearing) or without carbonates (CaCO 3 -free). We measured the SOC content and indicators of organic matter quality (C stable isotope composition, expressed as δ 13 C values, and thermal stability) in bulk soil samples. We then used sequential sonication and density fractionation (DF) to separate two occluded pools from free and mineral-associated SOC. The SOC content, mass, and δ 13 C values were determined in all the fractions. X-ray photoelectron spectroscopy was used to investigate the surface chemistry of selected fractions. Our hypothesis was that occlusion would be more prevalent at the CaCO 3 -bearing site due to the influence of Ca on aggregation, inhibiting oxidative transformation, and preserving lower δ 13 C values. Bulk SOC content was twice as high in the CaCO 3 -bearing profiles, which also had lower bulk δ 13 C values, and more occluded SOC. Yet, contrary to our hypothesis, occlusion only accounted for a small proportion of total SOC (< 10%). Instead, it was the heavy fraction (HF), containing mineral-associated organic C, which accounted for the majority of total SOC and for the lower bulk δ 13 C values. Overall, an increased Ca prevalence was associated with a near-doubling of mineral-associated SOC content. Future investigations should now aim to isolate Ca-mediated complexation processes that increase organo-mineral association and preserve organic matter with lower δ 13 C values.

58 GEOSCIENCES↗
DOE OSTI2023

Vertical Resolution Impacts Explicit Simulation of Deep Convection

Abstract The aggregation of tropical convection greatly influences the mean‐state of the atmosphere, altering humidity distributions, total atmospheric radiative cooling, and cloud amounts. Although studies have demonstrated the sensitivity of convective aggregation to horizontal resolution and domain size, few studies have explored the impact of vertical resolution on convective aggregation. Here, we investigate the impact of vertical resolution on simulations of deep convection and convective aggregation using the System for Atmospheric Modeling convection resolving model. We analyze simulations of tropical radiative‐convective equilibrium with varying vertical levels (32, 64, 128, and 256) across small (100 km), medium (700 km) and large (1,500 km) domains. We demonstrate that relative humidity and cloud fraction decrease with increasing vertical resolution as a result of reduced turbulent mixing. Vertical resolution also influences the occurrence of, onset time, and equilibrium intensity of aggregated convection, and also appears to affect the sensitivity of convective aggregation to domain size. Understanding how simulated convection aggregates, as well as its simulated sensitivity to model formulation, is critical for making and interpreting future predictions of global climate change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Soil pore network response to freeze-thaw cycles in permafrost aggregates

This dataset contains data used for the paper "Pore network response to freeze-thaw cycles in permafrost aggregates". The Related References field will be updated with a full citation when available.Climate change in Arctic landscapes may increase freeze-thaw frequency within the active layer as well as newly thawed permafrost. A highly disruptive process, freeze-thaw can deform soil pores and alter the architecture of the soil pore network with varied impacts to water transport and retention, redox conditions, and microbial activity. Our objective was to investigate how freeze-thaw cycles impacted the pore network of newly thawed permafrost aggregates to improve understanding of what type of transformations can be expected from warming Arctic landscapes. We measured the impact of freeze-thaw on pore morphology, pore throat diameter distribution, and pore connectivity with X-ray computed tomography (XCT) using six permafrost aggregates with sizes of 2.5 cm3 from a mineral soil horizon (Bw; 28-50 cm depths) in Toolik, Alaska. Freeze-thaw cycles were performed using a laboratory incubation consisting of five freeze-thaw cycles (-10˚C to 20˚C) over five weeks. Our findings indicated decreasing spatial connectivity of the pore network across all aggregates with higher frequencies of singly connected pores following freeze-thaw. Water-filled pores that were connected to the pore network decreased in volume while the overall connected pore volumetric fraction was not affected. Shifts in the pore throat diameter distribution were mostly observed in pore throats ranges of 100 microns or less with no corresponding changes to the pore shape factor of pore throats. Responses of the pore network to freeze-thaw varied with aggregate, suggesting that initial pore morphology may play a role in driving freeze-thaw response. Our research suggests that freeze-thaw alters the microenvironment of permafrost aggregates during the incipient stage of deformation following permafrost thaw, impacting soil properties and function in Arctic landscapes undergoing transition. This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/Erin-Rooney/XCT-freezethaw (DOI: 10.5281/zenodo.5816355).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2026

While legume-based rotation influences the chemical composition of mineral-associated organic matter, tillage has little effect on its persistence

Legume-based rotations are key to enhancing soil carbon (C) storage in stable forms such as mineral-associated organic matter (MAOM), improving soil health, and helping maintain long-term C stock. This study investigates how plant inputs and tillage practices influence the distribution, composition, and persistence of MAOM in wheat cropping systems when legumes are incorporated into the rotation. Soil samples were collected from Walla Walla silty loam of northeast Oregon, including conventional tillage (CT) and no-till (NT) treatments in a long-term dryland winter wheat (Triticum aestivum L.) and spring pea (Pisum sativum L.) rotation experiment, with an adjacent grass pasture (GP) plot serving as the baseline. Bulk soil, particulate organic matter (POM), and MAOM fractions were scanned using mid-infrared (Mid-IR) spectroscopy and analyzed for total C and nitrogen (N) as well as radiocarbon age, aggregate stability, and mineralizable C. In WP systems, NT promoted greater C accumulation and stabilization in MAOM than CT, reaching levels similar to those under native GP conditions. Furthermore, our findings indicate that decades of continuous legume-based rotation were insufficient to restore total bulk C content to levels observed in the undisturbed reference site (GP) in the 0–30 cm layer, with increases limited to the surface layer. The Mid-IR data provided insights into how cultivation may alter the chemical composition of MAOM. An increase in amide bands was observed in GP, while legume-based rotation systems favored an abundance of aliphatic (C-H) components in the MAOM fraction. In cultivated soils, C in MAOM exhibited longer residence times under arable conditions than in grassland, while POM remained an actively cycled pool across all treatments. Despite variations in the chemical composition of MAOM across treatments, mineralizable C emissions from MAOM fractions remained statistically unchanged during the 96-hour incubation period. Our findings suggest that soil physical protection, enhanced by reduced tillage, is the primary factor influencing MAOM cycling.

Carbon stabilization↗
DOE OSTI2021

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Relating Molecular Properties to the Persistence of Marine Dissolved Organic Matter with Liquid Chromatography–Ultrahigh-Resolution Mass Spectrometry

Marine dissolved organic matter (DOM) contains a complex mixture of small molecules that eludes rapid biological degradation. Spatial and temporal variations in the abundance of DOM reflect the existence of fractions that are removed from the ocean over different time scales, ranging from seconds to millennia. However, it remains unknown whether the intrinsic chemical properties of these organic components relate to their persistence. Here, we elucidate and compare the molecular compositions of distinct DOM fractions with different lability along a water column in the North Atlantic Gyre. Our analysis utilized ultrahigh-resolution Fourier transform ion cyclotron resonance mass spectrometry at 21 T coupled to liquid chromatography and a novel data pipeline developed in CoreMS that generates molecular formula assignments and metrics of isomeric complexity. Clustering analysis binned 14 857 distinct molecular components into groups that correspond to the depth distribution of semilabile, semirefractory, and refractory fractions of DOM. The more labile fractions were concentrated near the ocean surface and contained more aliphatic, hydrophobic, and reduced molecules than the refractory fraction, which occurred uniformly throughout the water column. These findings suggest that processes that selectively remove hydrophobic compounds, such as aggregation and particle sorption, contribute to variable removal rates of marine DOM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

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

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

42 ENGINEERING↗
DOE OSTI2021

Changing chemistry of particulate manganese in the near- and far-field hydrothermal plumes from 15°S East Pacific Rise and its influence on metal scavenging

Dissolved Mn(II) in the hydrothermal plume is known to be microbially oxidized to form Mn(III/IV) oxides, and the Mn oxides scavenge other trace elements in seawater. In the GEOTRACES GP16 cruise, dissolved Mn (dMn) and particulate Mn (pMn) were found to be transported over 4000 km westwards from the Southern Eastern Pacific Rise. Previous studies in this plume showed different removal rates of dMn and pMn as well as pMn size distribution between the near-field (<80 km from the ridge axis) and far-field (>80 km) plumes. In order to understand Mn cycling in these plumes, spatial distribution, oxidation states, and mineral structures of Mn in small size fraction (SSF; 0.8–51 μm) and large size fraction (LSF; >51 μm) particles from the near-field and far-field plumes were examined using micro X-ray fluorescence spectrometry (μ-XRF), X-ray absorption near-edge structure spectroscopy (XANES), chemical species mapping, and extended X-ray absorption fine-structure spectroscopy (EXAFS). In the near-field plume, pMn in the SSF is dominated by oxidized Mn with Mn(III) fractions of ~30%. They are a mixture of δ-MnO 2 and triclinic birnessites that is known to be formed as a result of autocatalytic Mn(II) oxidation at the surface of freshly-formed δ-MnO2, suggesting that both microbial and autocatalytic Mn oxidation occur in the near-field plume. The LSF pMn in the near-field plume is also oxidized and often found in large aggregates several hundreds of μm in size. These aggregates settle out in the near-field and during transport, and are not found in the far-field plume. In the far-field plume where Mn oxides are not newly formed, pMn in the SSF is oxidized, but their Mn(III) fractions are smaller than in the near-field pMn. Unlike the SSF, the far-field plume LSF pMn is dominated by reduced Mn, implying very slow aggregation of pMn in the far-field plume. The different characteristics of pMn between the near-field and far-field plumes affect its scavenging of other trace elements. In the near-field plume, Co, Mo, 231 Pa are associated with pMn, but not in the far-field plume. 231 Pa is adsorbed to pFe rather than pMn in the far-field plume, and Pb is adsorbed to pFe in the entire plume. The result shows that freshly-formed Mn oxides in the near-field plume have higher scavenging capacity than the far-field plume pMn. Our findings suggest that the mineralogical age of Mn oxides may be an important parameter that controls the scavenging of many other trace elements and isotopes.

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

Quantum computing approach for building surface sunlit in urban-scale energy modeling

Solar shadow calculations are needed in building energy modeling and performance simulation of PV systems installed on roofs or facades of buildings. We present a quantum computing approach for calculation of building surface sunlit fractions by recasting solar visibility as a binary optimization problem solved by quantum annealing. Each triangulated surface centroid is encoded as a binary qubit indicating sunlit or shaded status. Geometric visibility constraints are derived from the Möller-Trumbore intersection algorithm and converted into a constrained quadratic binary model compatible with contemporary quantum annealers. The coefficients were embedded to D-Wave quantum computer. To demonstrate feasibility, we conducted a case study in San Francisco for a target building with 52 triangles and roughly 2700 nearby triangles within 50 m evaluated at representative winter and summer solar positions. The results demonstrated that quantum annealing can reliably calculate and distinguish sunlit from shaded surfaces. Quantum samples achieved average accuracy exceeding 92.4 %, with the aggregate surface-level agreement approaching 99.9 %. The outputs of quantum computers agreed closely with classical algorithms, indicating practical feasibility and promising scalability. Finally, the hourly sunlit fractions of building surfaces can be obtained for urban energy modelling. This is the first study to apply quantum computing to the solar shadow and building surface sunlit calculation. It introduces a new paradigm that differs fundamentally from traditional approaches.

Deng, Zhipeng↗
DOE OSTI2022

On the provenance of GEMS, a quarter century post discovery

The provenance of GEMS (glass with embedded metal and sulfides) in cometary type interplanetary dust particles is investigated using analytical scanning transmission electron microscopy and secondary ion mass spectrometry. We review the current state of knowledge and closely examine the densities, elemental compositions and distributions, iron oxidation states, and isotopic compositions of a subset of GEMS in chondritic porous interplanetary dust. We find that GEMS are underdense with estimated densities that are 35–65% of compositionally equivalent crystalline aggregates. GEMS low densities result in a lower contribution to the bulk compositions of IDPs than has been assumed based on their volume fraction. We also find that element/Si ratios, assumed to be primary (indigenous), are instead perturbed by contamination and secondary alteration, including pulse heating during atmospheric entry. Fe in pyrrhotite inclusions was oxidized and Mg, S, Ca, and Fe were depleted relative to lithophile Al and Si, resulting in reduction in element/Si ratios. Because they trap outgassing elements, Fe-rich oxide rims that formed on the surface of GEMS are serendipitous “witness plates” to the changes in composition that accompany atmospheric entry. As a result of alteration, GEMS elemental compositions cannot reliably inform about their provenance. Except for highly anomalous oxygen isotope ratios measured in some large GEMS grains that indicate a contribution from circumstellar dust, oxygen isotope compositions are generally poor indicators of provenance. Prior work indicates that most GEMS fall close to the terrestrial oxygen isotope composition, which, however, does not exclude a presolar interstellar origin. Nitrogen isotopic compositions are more diagnostic. Elevated 15 N/ 14 N ratios indicate that GEMS accreted in conjunction with formation of organic matter by ion–molecule reactions in a cold (<50 K) presolar environment like the extreme outer nebula or interstellar medium. Considering all observations, we conclude that GEMS are most likely processed interstellar silicates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗
DOE OSTI2023

Linking water quality, fouling layer composition, and performance of reverse osmosis membranes

Fouling of polyamide membranes during reverse osmosis (RO) is a major challenge for adopting membrane technologies to treat highly contaminated waters, especially those containing organic foulants (e.g., natural organic matter (NOM), polysaccharides) and dominant cations (e.g., sodium, magnesium, calcium). Here, this work combines bench-scale membrane fouling experiments with detailed characterization of feedwater chemistry and fouling layer composition/morphology to reveal fundamental mechanisms of (in)organic fouling during RO. Divalent cations are shown to promote fouling by hydrophobic NOM containing aromatic and carboxyl groups, while NOM fouling in the presence of a monovalent cation, sodium, occurs by smaller fulvic acids containing larger fractions of carboxyl groups and other oxygen-rich moieties. Calcium-carboxyl bridging occurs in solution and near the membrane surface to induce NOM aggregation on nanometer length scales. In complex waters containing foulant mixtures, co-fouling by calcium-carboxyl bridging and CaCO 3 precipitation influence membrane performance at longer timeframes. However, the flux decline observed for the co-fouling mechanism was less significant than the sum of its parts, suggesting both synergistic and antagonistic fouling mechanisms should be considered in membrane design/operation. These results encourage the design of pretreatment processes to reduce concentrations of multivalent ions and hydrophobic NOM in RO feedwaters, and of membrane materials to limit attachment/deposition of aggregates to/on polyamide surfaces.

42 ENGINEERING↗