Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring new frontiers in type 1 diabetes through advanced mass-spectrometry-based molecular measurements

Type 1 diabetes (T1D) is a devastating autoimmune disease for which advanced mass spectrometry (MS) methods are increasingly used to identify new biomarkers and better understand underlying mechanisms. For example, integration of MS analysis and machine learning has identified multimolecular biomarker panels. In mechanistic studies, MS has contributed to the discovery of neoepitopes, and pathways involved in disease development and identifying therapeutic targets. However, challenges remain in understanding the role of tissue microenvironments, spatial heterogeneity, and environmental factors in disease pathogenesis. Recent advancements in MS, such as ultra-fast ion-mobility separations, and single-cell and spatial omics, can play a central role in addressing these challenges. Here, in this work, we review recent advancements in MS-based molecular measurements and their role in understanding T1D.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Detailed space–time variations of the seismic response of the shallow crust to small earthquakes from analysis of dense array data

SUMMARY We compute high-resolution space–time variations of subsurface seismic properties from autocorrelation functions (ACF’s) of noise and local earthquakes, recorded by the Sage Brush Flat dense array deployed around the Clark branch of the San Jacinto fault. The resolved temporal changes are referred to as apparent velocity changes because they reflect both nonlinear response and variations of material properties such as cracking and damage. Apparent velocity changes are estimated at four frequency bands (10–15, 10–20, 15–30 and 20–40 Hz) for two local earthquake data sets. In one analysis, ACF’s from P- and S-wave windows of 31 small events with magnitudes below 3.1 are used to compute the apparent velocity variations with respect to the mean ACF of each phase, and we also use the mean ACF of noise data as reference to estimate the changes. In a further analysis, the temporal evolution of properties is computed using moving time windows in continuous waveform over one-hour long data with noise and earthquake signals. The apparent velocity changes and recovery times are frequency dependent and present a strong spatial variability across the array. The resolved changes are larger and recovery time shorter with data associated with higher frequencies. At frequencies larger than 15 Hz, and using the mean ACF of noise data as a reference, the apparent average velocity changes across the array during the passage of the P and S waves from the small local events are 2.5 per cent and 6 per cent, respectively. The apparent velocity changes decrease by one order of magnitude when the earthquake data are used as a reference. The relatively large changes in response to very low ground motion have important implications on nonlinear processes involving degradation and healing of the subsurface material during common earthquake shaking.

Geochemistry & Geophysics↗
DOE OSTI2024

The spatially resolved transcriptome signatures of glomeruli in chronic kidney disease

Here, we used digital spatial profiling (DSP) to describe the glomerular transcriptomic signatures that may characterize the complex molecular mechanisms underlying progressive kidney disease in Alport syndrome, focal segmental glomerulosclerosis, and membranous nephropathy. Our results revealed significant transcriptional heterogeneity among diseased glomeruli, and this analysis showed that histologically similar glomeruli manifested different transcriptional profiles. Using glomerular pathology scores to establish an axis of progression, we identified molecular pathways with progressively decreased expression in response to increasing pathology scores, including signal recognition particle–dependent cotranslational protein targeting to membrane and selenocysteine synthesis pathways. We also identified a distinct signature of upregulated and downregulated genes common to all the diseases investigated when compared with nondiseased tissue from nephrectomies. These analyses using DSP at the single-glomerulus level could help to increase insight into the pathophysiology of kidney disease and possibly the identification of biomarkers of disease progression in glomerulopathies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Feasibility of Adding Twitter Data to Aid Drought Depiction: Case Study in Colorado

The use of social media, such as Twitter, has changed the information landscape for citizens’ participation in crisis response and recovery activities. Given that drought progression is slow and also spatially extensive, an interesting set of questions arise, such as how the usage of Twitter by a large population may change during the development of a major drought alongside how the changing usage facilitates drought detection. For this reason, contemporary analysis of how social media data, in conjunction with meteorological records, was conducted towards improvement in the detection of drought and its progression. The research utilized machine learning techniques applied over satellite-derived drought conditions in Colorado. Three different machine learning techniques were examined: the generalized linear model, support vector machines and deep learning, each applied to test the integration of Twitter data with meteorological records as a predictor of drought development. It is found that the integration of data resources is viable given that the Twitter-based model outperformed the control run which did not include social media input. Eight of the ten models tested showed quantifiable improvements in the performance over the control run model, suggesting that the Twitter-based model was superior in predicting drought severity. Future work lies in expanding this method to depict drought in the western U.S.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

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

Spatial Distribution of Ultraviolet Emission from Cometary Activity at 67P/Churyumov-Gerasimenko

The Alice ultraviolet spectrograph on board the Rosetta orbiter provided the first near-nucleus ultraviolet observations of a cometary coma from arrival at comet 67P/Churyumov-Gerasimenko in 2014 August through 2016 September. The characterization of atomic and molecular emissions in the coma revealed the unexpected contribution of dissociative electron impact emission at large heliocentric distances and during some outbursts. This mechanism also proved useful for compositional analysis, and Alice observed many cases that suggested elevated levels of the supervolatile O{sub 2}, identifiable in part to their emissions resulting from dissociative electron impact. In this paper, we present the first two-dimensional UV maps constructed from Alice observations of atomic emission from 67P during an increase in cometary activity on 2015 November 7–8. Comparisons to observations of the background coma and an earlier collimated jet are used to describe possible changes to the near-nucleus coma and plasma. To verify the mapping method and place the Alice observations in context, comparisons to images derived from the MIRO and VIRTIS-H instruments are made. The spectra and maps we present show an increase in dissociative electron impact emission and an O{sub 2}/H{sub 2}O ratio of ~0.3 for the activity; these characteristics have been previously identified with cometary outbursts seen in Alice data. Further, UV maps following the increases in activity show the spatial extent and emission variation experienced by the near-nucleus coma, informing future UV observations of comets that lack the same spatial resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Propagating neutronic uncertainties for FFTF LOFWOS Test #13

The safety evaluation conducted for licensing a Sodium-cooled Fast Reactor (SFR) may require modeling transients with best-estimate calculation tools that must first be validated against real-world measurements. To provide the community with a valuable benchmarking opportunity for validating SFR analysis tools and methods, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13. The impact of nuclear data uncertainties on neutronics parameters was previously investigated based on the COMMARA-2.0 covariance matrix. Since the transient simulation results are very sensitive to certain reactivity coefficients, it was decided to employ rigorous uncertainty propagation methods to quantify the impact of nuclear data uncertainties on the best-estimate predication of FFTF LOFWOS Test #13. The DAKOTA code is used to propagate neutronic uncertainties through SAS4A/SASSYS-1 transient simulations, while taking into account spatial and reaction-wise correlations within these uncertainties. This study shows that the remaining discrepancies observed between the Argonne National Laboratory (ANL) best-estimate results and the experimental measurements can be partly explained by the uncertainty associated with Gas Expansion Module (GEM) worth, which contributes the majority of the overall nuclear data uncertainty on the output from the FFTF LOFWOS Test #13 transient simulation. This study also confirmed the importance of including spatial and reaction-wise correlations of nuclear data uncertainties on feedback coefficients in the uncertainty propagation to avoid under-estimating their impact during the transient simulations.

42 ENGINEERING↗
DOE OSTI2022

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multiresolution convolutional autoencoders

Herein we propose a multi-resolution convolutional autoencoder (MrCAE) architecture that integrates and leverages three highly successful mathematical architectures: (i) multigrid methods, (ii) convolutional autoencoders and (iii) transfer learning. The method provides an adaptive, hierarchical architecture that capitalizes on a progressive training approach for multiscale spatio-temporal data. This framework allows for inputs across multiple scales: starting from a compact (small number of weights) network architecture and low-resolution data, our network progressively deepens and widens itself in a principled manner to encode new information in the higher resolution data based on its current performance of reconstruction. Basic transfer learning techniques are applied to ensure information learned from previous training steps can be rapidly transferred to the larger network. As a result, the network can dynamically capture different scaled features at different depths of the network. The performance gains of this adaptive multiscale architecture are illustrated through a sequence of numerical experiments on synthetic examples and real-world spatial-temporal data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An adaptive synchronous extraction (ASE) method for estimating intensity and footprint of surface urban heat islands: A case study of 254 North American cities

The urban heat island (UHI) effect has attracted great attention due to its potential impacts on rapidly growing urban areas. Using remotely sensed estimates of land surface temperature (LST), a large number of studies have focused on the surface UHI (SUHI) effect, which can be characterized by its two fundamental properties: intensity and footprint. The SUHI intensity reflects the LST difference between the urban area and the background reference area (BRA), and the SUHI footprint indicates the spatial extent influenced by the heat island. Currently, numerous methods have been developed to estimate the SUHI intensity and footprint, but are still greatly challenged by three main issues. Namely, the discrepancy in BRA selection criterion brings great uncertainty to the estimated SUHI intensity, the estimation of SUHI footprint is largely constrained by the predefined models, and the quantification of SUHI effect is potentially influenced by several confounding factors. Here, we proposed an adaptive synchronous extraction (ASE) method, which is capable of adaptively selecting the most optimal BRA while removing the influence of confounding factors, and achieving synchronous estimation of SUHI intensity and footprint. We applied the ASE method to 254 North American cities and conducted an in-depth comparative analysis to discuss its applicability and benefits. The main results include: (1) The ASE method avoids the limitations of existing methods in BRA selection and model presetting, and shows resilience to parameter variations. This makes the ASE method highly applicable to quantify the SUHI intensity and footprint in cities with various thermal characteristics. (2) The ASE method can better highlight the spatial, seasonal and day-night contrasts in the estimated SUHI intensity. This superiority is particularly evident when comparing it to methods based on the equal-area buffer or the simplified urban-extent algorithm. (3) Confounding factors pose non-negligible impacts on the quantification of the SUHI effect. Typically, ignoring the influence of topographic relief or missing LST data can lead to an overall overestimation of the SUHI intensity, while not removing surrounding urban areas will cause some underestimation of the SUHI intensity. In conclusion, overall, the proposed ASE method provides a new generalizable tool for quantifying the SUHI effect, which has great potentials for future studies and urban climate assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Asymptotic hydrogen redistribution analysis in yttrium-hydride-moderated heat-pipe-cooled microreactors using DireWolf

Yttrium hydride (YH x ) is one of the materials being considered for moderating thermal and epithermal nuclear microreactors. One potential issue with YH x use is that the hydrogen redistributes in the hydride when thermal and concentration gradients are present. This hydrogen redistribution leads to spatial gradients in the hydrogen concentration, thus affecting neutron transport in the reactor. Here, by building upon observations in prior works, this paper aims to gain a better understanding of the reactivity feedback associated with such hydrogen redistributions. In particular, we wish to understand the sign (+/–) of the hydrogen redistribution neutronic feedback, its order of magnitude, and its underlying physical causes. To achieve this goal, the DireWolf multiphysics software driver was used to solve the coupled radiation transport, heat transfer, heat pipe two-phase flow, and hydrogen redistribution equations for the Simplified Microreactor Benchmark Assessment (SiMBA) problem, a full-core microreactor numerical benchmark developed at Idaho National Laboratory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Demographic analysis of invasible habitat fraction identifies context‐dependent roles of resource availability and biotic resistance in determining invasion success

Abstract Theories of plant invasions predict that plant communities should be more easily invaded when resources increase and/or competition decreases. We tested this with an experimentally introduced plant population by manipulating precipitation and resident community biomass. We used a spatially explicit demographic approach to develop a new population‐level metric of invasibility that quantifies the invasible habitat fraction (IHF) across the landscape. The existing community was essentially uninvasible (median IHF ≈ 0%), but experimental manipulations greatly increased the range of outcomes, with maximum observed IHF values over 50%. However, changes in invasibility were often context‐dependent, resulting in some outcomes that aligned with existing theory, and others that were not readily predicted. Moreover, variation in invasibility was often driven by specific sets of invader demographic vital rates. Removing competitors revealed the capacity for strong biotic resistance, but this interacted with precipitation such that little biotic resistance was detected under drought conditions. Adding precipitation typically had little positive effect on invasibility, and moderate drought relief led to relatively high invasibility. However, the latter was driven to a large extent by interactions with mammal herbivory that otherwise inhibited invasion in one year. Synthesis . Our findings show that interactions between abiotic and biotic factors, as well as legacy effects, can strongly mediate invasibility. This study also highlights the importance of incorporating spatial heterogeneity into population‐level assessments of invasion, as initial population declines do not necessarily indicate resistance to invasion.

Griffith, Alden B.↗
DOE OSTI2022

GROWdb US River Systems - Samples

GROW Overview We developed the Genome Resolved Open Watersheds database (GROWdb), which aims to increase genomic sampling and understanding of global river microbiomes. An emphasis of GROWdb is to create a publicly available and ever-expanding microbial genome database that is focused on rivers while being interoperable with databases from other ecosystems. GROWdb is based on a network-of-networks approach to move beyond a small collection of well-studied rivers, towards a spatially distributed, global network of systematic observations. GROWdb represents the first microbial, river-focused resource parsed at various scales from genes to MAGs to community level including expression and potential based measurements that will be of interest to microbiologists, ecologists, geochemists, hydrologists, and modelers. Dataset Acknowledgement GROWdb contains data from various research campaigns, please acknowledge the following data generators, as appropriate: WHONDRS derived genomes or samples - include this statement in your acknowledgements: “This study used data from the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) under the River Corridor Science Focus Area (SFA) at the Pacific Northwest National Laboratory (PNNL) that was generated at the U.S. Department of Energy (DOE) Joint Genome Institute User Facility. PNNL is operated by Battelle Memorial Institute for the U.S. DOE under Contract No. DE-AC05-76RL01830. The SFA is supported by the U.S. DOE, Office of Biological and Environmental Research (BER), Environmental System Science (ESS) Program.” Total Samples loaded onto this Narrative: 178 Note: Not all GROW samples may be loaded into KBase Data Availability The data underlying GROWdb are accessible across various platforms to ensure all levels of data structure are widely available. First, all reads and MAGs are publicly hosted on National Center for Biotechnology (NCBI) under Bioproject PRJNA946291. Second, all data related data presented here including MAG annotations, extended data tables, phylogenetic tree files, antibiotic resistance gene database files, and MAG abundance tables are available in Zenodo (link). Beyond the flat database files listed above, our aim for GROWdb was to maximize data use by making the data available in searchable and interactive platforms including the National Microbiome Data Collaborative (NMDC) data portal, the Department of Energy’s Systems Biology Knowledgebase (KBase), and a GROW specific user interface released here, GROWdb Explorer. Each platform provides different ways to interact with GROWdb: NMDC GROWdb formed a pilot project for the NMDC. Specifically, individual GROWdb datasets (metagenomes, metatranscriptomes, etc) are easily accessible and searchable through the NMDC data portal, where they are systematically connected to each other and to a rich suite of sample information and standard analysis results, following Findable, Accessible, Interoperable, and Reusable (FAIR) data practices. KBase GROWdb is publicly available within KBase, including samples (this Narrative), MAGs, and corresponding genome scale metabolic models. Access within KBase allows for immediate access and reuse of data, including comparison to private data using KBase’s 500+ analysis tools. Other linked narratives in KBase: GROW Metagenome Assembled Genomes (MAGs) GROW Metabolic Models GROWdb Explorer GROWdb data is also explorable through a graphical user interface built through the Colorado State University Geospatial Centroid (https://geocentroid.shinyapps.io/GROWdatabase/), allowing users to search and graph microbial and spatial data simultaneously. In summary, this microbial genome resource represents the first publicly available genome collection from rivers and offers data that can be leveraged across microbiome studies. GROWdb is an expanding repository to incorporate and unify global river multi-omic data for the future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
DOE OSTI2025

Benchmarking image processing techniques for porosity measurement in polymer additive manufacturing: Review and experimental analysis

An image processing workflow is proposed for porosity measurement in polymer additive manufacturing. Various techniques, including global and local thresholding, region growing, and K-means clustering, were applied to microscopic images of carbon fiber reinforced acrylonitrile butadiene styrene (CF-ABS) and benchmarked for their ability to accurately measure porosity. Global methods included Otsu, minimum error, iterative, and entropy-based thresholding, while local methods included Niblack, Bernsen, Sauvola, and Bradley-Roth algorithms. Artificial uneven illumination was introduced to test local adaptive thresholds. Results showed significant differences in porosity values across methods. Otsu, region growing, and K-means clustering excelled under uniform illumination, while Sauvola and Bradley-Roth performed better with uneven illumination. Comparison with X-ray computed tomography (XCT) revealed slightly lower porosity values (2.55 %) than optimized methods (2.73–2.79 %) due to XCT's lower resolution excluding smaller pores. While XCT offers finer pore detection, it limits sample volume and underestimates porosity due to spatial variation. Validation using artificial grayscale images with 5 % porosity confirmed that Otsu, Bradley-Roth, region growing, and Sauvola algorithms produced accurate results. Although tested on a single material system, these methods can be adapted to others with optimization. In conclusion, given XCT's high computational and time costs, this study highlights suitable image processing techniques as cost-effective alternatives for porosity analysis in polymer composites.

Additive manufacturing↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2023

An empirical calibration of the serpentine-water oxygen isotope fractionation at T = 25–100°C

The ability to precisely constrain the physical–chemical conditions of serpentinization, such as temperature and fluid source, is limited by the accuracy of calibrations for oxygen isotope fractionation between serpentine and water – i.e., 1000 lnα( serpentine-water ) – which disagree by up to 20‰ when extrapolated to T<200°C. In this study, we present a new empirical calibration of 1000 lnα( serpentine-water ) aiming to improve applications of oxygen isotope geochemistry to very low-T serpentinization (T<100°C). We used the high-spatial resolution capabilities of Secondary Ion Mass Spectrometry (SIMS) to analyze oxygen isotope ratios in mineral pairs of calcite + serpentine, quartz + serpentine and talc + serpentine co-crystallized at scales≤50μm in six serpentinite samples from the Samail ophiolite (Oman). SIMS analysis shows that the mineral pairs are relatively homogenous in oxygen isotope ratios with variability in δ 18 O values ≤2‰ (2s). Clumped isotope thermometry and petrological constraints indicate crystallization temperatures from 25 to 100°C for the investigated samples. These independent constraints on temperature allowed us to derive 1000 lnα( serpentine-water ) by combining Δ 18 O( quartz-serpentine ), Δ 18 O( talc-serpentine ) and Δ 18 O( calcite-serpentine ) measured by SIMS in the investigated samples with Δ 18 O( quartz-water ), Δ 18 O( talc-water ) and Δ 18 O( calcite-water ) calibrations available in the literature. Our empirical calibration of 1000 lnα( serpentine-water ) = 1.04±0.20 (2SE) × 10 6 /T 2 (T in K), from T = 25 to 100°C, is within uncertainty of former high-temperature empirical calibrations extrapolated to T<100°C and matches experimental calibrations when extrapolated to T=250°C. In conclusion, the new serpentine-water calibration enables more accurate reconstructions of fluid-rock interactions occurring during low-temperature serpentinization processes in various tectonic settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗