SEARCH · Engineering Papers

Results for “Shape Matching”

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗

DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗

DOE OSTI2023

An experiment-informed discrete element modelling study of knife milling for flexural biomass feedstocks

A discrete element method (DEM) based approach is used to study the relationships between material attributes (MAs), processing parameters (PPs), and quality attributes (QAs) for the knife milling of maize stalks. An approximate DEM shape model was conceptualized based on real maize stalks and calibrated based on experimental bending test data for flexural properties (elastic bending stiffness, elastic bending angle limit, elastoplastic ratio, etc.). DEM simulations of maize stalk comminution in a Jordan Reduction Solutions (“JRS”) knife mill were performed to investigate the relationships between the MAs (maize stalk size and breakage stress limit), PPs (impeller rotational speed), and QAs (mass throughput and output particle size distribution (PSD)). The DEM results suggest that stalk length has little influence on mass throughput and PSD, whilst stalks with larger cross sections tend to generate larger sizes of milled particles given the same breakage stress limit. Both the DEM and experimental results show that faster impeller rotation (or higher power) does not necessarily generate higher throughput or smaller output PSD, especially for maize stalks of higher breakage stress limit. The correlations between these MAs, PPs and QAs are found highly stochastic, though breakage stress limit dictates mass throughput, regardless of stalk size. The DEM-predicted output particle size tended to match the experimental data with coarse PSDs based on sieve size but showed weakened fidelity with finer material, indicating the potential for further model improvement.

09 BIOMASS FUELS↗

DOE OSTI2024

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Here, in this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Phyllosphere exudates select for distinct microbiome members in sorghum epicuticular wax and aerial root mucilage

Phyllosphere exudates create specialized microhabitats that shape microbial community diversity. We explored the microbiome associated with two sorghum phyllosphere exudates, the epicuticular wax and aerial root mucilage. We assessed the microbiome associated with the wax from sorghum plants over two growth stages, and the root mucilage additionally from nitrogen-fertilized and non-fertilized plants. In parallel, we isolated and characterized hundreds of bacteria from wax and mucilage, and integrated data from cultivation-independent and cultivation-dependent approaches to gain insights into exudate diversity and bacterial phenotypes. We found that Sphingomonadaceae and Rhizobiaceae families were the major taxa in the wax regardless of water availability and plant developmental stage to plants. The cultivation-independent mucilage-associated bacterial microbiome contained Erwiniaceae, Flavobacteriaceae, Rhizobiaceae, Pseudomonadaceae, Sphingomonadaceae, and its structure was strongly influenced by sorghum development but only modestly influenced by fertilization. In contrast, the fungal community structure of mucilage was strongly affected by the year of sampling but not by fertilization or plant developmental stage, suggesting a decoupling of fungal-bacterial dynamics in the mucilage. Our bacterial isolate collection from wax and mucilage had several isolates that matched 100% to detected amplicon sequence variants, and were enriched on media that selected for phenotypes including phosphate solubilization, putative diazotrophy, resistance to desiccation, capability to grow on methanol as a carbon source, and ability to grow in the presence of linalool and β-caryophyllene (terpenes in sorghum wax). This work expands our understanding of the microbiome of phyllosphere exudates and supports our long-term goal to translate microbiome research to support sorghum cultivation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

An acoustic investigation of the near-surface turbulence on Mars

The Perseverance rover is carrying out an original acoustic experiment on Mars: the SuperCam microphone records the spherical acoustic waves generated by laser sparks at distances from 2 m to more than 8 m. These N-shaped acoustic waves scatter from the multiple local heterogeneities of the turbulent atmosphere. Therefore, large and random fluctuations of sound travel time and intensity develop as the waves cross the medium. The variances of the travel times and the scintillation index (normalized variance of the sound intensity) are studied within the mathematical formalism of the propagation of spherical acoustic waves through thermal turbulence to infer statistical properties of the Mars atmospheric temperature fluctuation field. The comparison with the theory is made by simplifying assumptions that do not include wind fluctuations and diffraction effects. Two Earth years (about one Martian year) of observations acquired during the maximum convective period (10:00–14:00 Mars local time) show a good agreement between the dataset and the formalism: the travel time variance diverges from the linear Chernov solution exactly where the density of occurrence of the first caustic reaches its maximum. Moreover, on average, waves travel faster than the mean speed of sound due to a fast path effect, which is also observed on Earth. To account for the distribution of turbulent eddies, several power spectra are tested and the best match to observation is obtained with a generalized von Karman spectrum with a shallower slope than the Kolmogorov cascade, φ(k)∝(1 + k 2 L 2 ) -4/3 . It is associated with an outer scale of turbulence, L, of 11 cm at 2 m above the surface and a standard deviation of 6 K over 9 s for the temperature. These near-surface atmospheric properties are consistent with a weak to moderate wave scattering regime around noon with little saturation. Overall, this study presents an innovative and promising methodology to probe the near-surface atmospheric turbulence on Mars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Molecular and structural characterization of a Bacillus cereus strain producing an anthrax-like capsule

Bacillus cereus is a ubiquitous Gram-positive, spore-forming, rod-shaped saprophytic bacterium, occasionally reported to cause food-borne illnesses. However, instances of B. cereus strains harboring anthrax toxin and capsule genes have elevated certain strains as formidable pathogens and biothreats. This study focuses on the genomic analysis and the structural characterization of capsular material produced by the virulent B. cereus PATH2418 strain, isolated from the wound of a traumatic open fracture patient. The genome was sequenced using Nanopore MinION sequencing, revealing a chromosome of 5,270,283 bp and three plasmids. One plasmid, pATH1, was found to encode an operon for the biosynthesis of a bacterial capsule. This operon had sequence homology to the Bacillus anthracis capBCADE operon, which encodes the poly-γ-D-glutamate (PDGA) capsule. The capsule production in B. cereus PATH2418 was influenced by temperature and CO 2 levels. Structural analysis of the capsular material using a combined approach of nuclear magnetic resonance (NMR) and high-performance liquid chromatography (HPLC) techniques confirmed the presence of a high-molecular-weight poly-γ-glutamate capsule, with an enantiomeric composition of approximately 67% D-glutamic acid and 33% L-glutamic acid, matching that of B. anthracis.

Bacillus cereus↗

DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Measurement of the 1-jettiness event shape observable in deep-inelastic electron-proton scattering at HERA

The H1 Collaboration reports the first measurement of the 1-jettiness event shape observable $τ^b_1$ in neutral-current deep-inelastic electron-proton scattering (DIS). The observable $τ^b_1$ is equivalent to a thrust observable defined in the Breit frame. The data sample was collected at the HERA ep collider in the years 2003–2007 with center-of-mass energy of $\sqrt{s}$ = 319 GeV, corresponding to an integrated luminosity of 351.1 pb -1 . Triple differential cross sections are provided as a function of $τ^b_1$, event virtuality $Q^2$, and inelasticity y, in the kinematic region $Q^2$ > 150 GeV 2 . Single differential cross section are provided as a function of $τ^b_1$ in a limited kinematic range. Double differential cross sections are measured, in contrast, integrated over $τ^b_1$ and represent the inclusive neutral-current DIS cross section measured as a function of $Q^2$ and y. The data are compared to a variety of predictions and include classical and modern Monte Carlo event generators, predictions in fixed-order perturbative QCD where calculations up to $\mathcal{O}$($α^3_s$) are available for $τ^b_1$ or inclusive DIS, and resummed predictions at next-to-leading logarithmic accuracy matched to fixed order predictions at $\mathcal{O}$($α^2_s$). These comparisons reveal sensitivity of the 1-jettiness observable to QCD parton shower and resummation effects, as well as the modeling of hadronization and fragmentation. Within their range of validity, the fixed-order predictions provide a good description of the data. Monte Carlo event generators are predictive over the full measured range and hence their underlying models and parameters can be constrained by comparing to the presented data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Advanced sorbents for modular oxygen production for REMS gasifiers

Under a DOE funded effort, Thermosolv LLC has been developing a sorbent-based oxygen production technology for gasification and oxy-combustion applications. The sorbents utilize oxygen-storage properties of certain perovskites to (1) selectively adsorb oxygen at moderate temperature from compressed air and (2) release the adsorbed oxygen into a vacuum or a sweep gas such as CO2 and/or steam. Through cyclic operations of multiple sorbent beds such a process can be made continuous. Pressure drop across the sorbent bed and other similar design considerations dictate that the sorbent be in the form of a high surface-area-to-volume pellet and yet still possess sufficient crush strength and integrity. Process efficiency is determined by the sorption and desorption kinetics and sorbent capacity. Typical process involves cycle times in the order of 100 seconds, a time too short to fully utilize the full volume of the sorbent and the sorption capacity of the sorbent pellets. As a part of this project, Thermosolv LLC undertook the development of oxygen sorbents as a high-surface area supported sorbent on a light-weight inert support to (1) reduce the cost, (2) increase the productivity, and (3) reduce the overall weight of the reactor. Using a previously developed perovskite (LSCF 1991) with its well-characterized performance, composite sorbent pellets each consisting of an inert core coated by a thin layer of the functional perovskite material were produced. Several low-density, inexpensive inert solids supports in the 1/8”-1/4” size were selected to keep the pressure drop across the sorbent bed in the acceptable range. A number of commercial technologies including spray coating, dip coating, incipient impregnation of porous substrate followed by thermal annealing at various temperatures and duration were employed and tested to produce robust composite supported sorbent pellets. A sense of optimum coating thickness was developed by testing sorbent pellets of various increasing diameter pellets. LSCF-1991extrudates ranging in size from 1/32” to 3/16” were tested in a TGA and in a fixed-bed pressure swing test set-up for sorption/desorption cycles of interest. The data show that for the operational conditions of interest the coating thickness for a supported sorbent pellet needs to be approximately 1/32” (about 0.8 mm). Subsequent work thereby concentrated on developing composite pellets of various substrates, shapes and sizes with coating layer of about this size range. Candidate support materials were chosen based on cost, inertness, mechanical strength, thermal expansion and chemical stability with sorbent material, and in a size range to give an acceptable pressure drop in fixed-bad reactor configurations. Coating application methods used for application of the sorbent onto the support included precipitation, spray coating and dip coating from sorbent slurries and sol gels. For all coated supports where we could successfully apply a uniform coating of desired thickness with an acceptable handling performance in terms of exfoliation, TGA-based cyclic sorption/desorption testing was performed to determine cycling capacity of oxygen. Among nearly twenty different support materials tested, best adhesion performance was obtained from stainless saddle supports. In the bench-scale fixed-bad tests stainless steel saddle support composite pellets with a coating thickness in the 0.6 mm or so range, the composite sorbent pellet performance approached up to 95% of that of the 2 mm parent material pellets. In a parallel approach to reducing the cost of the sorbent, alternate sorbent formulations replacing/reducing the amount of cobalt in the LSCF family were also investigated. Successful formulations that could match the performance of LSCF 1991 were identified based on TGA cyclic tests as LSCF 1919 and LSF 1910. In the range of operational envelop of cycle times and other relevant process operating conditions, the reduced Co formulations showed comparable performance in the bench-scale fixed-bad cyclic operations. As a part of this project, we also attempted substitution of La and Sr with Ca and Ba, and Co with Mn, Ni and Cu with little success, but the overall project goal of reducing sorbent cost in terms of raw material and manufacturing expenses was successful.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Quantifying Uncertainty in Ice Particle Velocity–Dimension Relationships Using MC3E Observations

Abstract Ice particle terminal fall velocity ( V t ) is fundamental for determining microphysical processes, yet remains extremely challenging to measure. Current theoretical best estimates of V t are functions of Reynolds number. The Reynolds number is related to the Best number, which is a function of ice particle mass, area ratio ( A r ), and maximum dimension ( D max ). These estimates are not conducive for use in most models since model parameterizations often take the form , where ( α , β ) depend on habit and D max . A previously developed framework is used to determine surfaces of equally plausible ( α , β ) coefficients whereby ice particle size/shape distributions are combined with V t best estimates to determine mass- ( V M ) or reflectivity-weighted ( V Z ) velocities that closely match parameterized V M ,SD or V Z ,SD calculated using the ( α , β ) coefficients using two approaches. The first uses surfaces of equally plausible ( a , b ) coefficients describing mass ( M )–dimension relationships (i.e., ) to calculate mass- or reflectivity-weighted velocity from size/shape distributions that are then used to determine ( α , β ) coefficients. The second investigates how uncertainties in A r , D max , and size distribution N ( D ) affect V M or V Z . For seven of nine flight legs flown on 20 and 23 May 2011 during the Mesoscale Continental Convective Clouds Experiment (MC3E), uncertainty from natural parameter variability—namely, the variability in ice particle parameters in similar meteorological conditions—exceeds uncertainties arising from different A r assumptions or D max estimates. The combined uncertainty between A r , D max , and N ( D ) produced smaller variability in ( α , β ) compared to varying M ( D ), demonstrating M ( D ) must be accurately quantified for model fall velocities. Primary sources of uncertainty vary considerably depending on environmental conditions. Significance Statement Ice particle fall velocity is fundamental for numerous processes within clouds, and hence is a critical property that must be accurately represented in weather and climate models. Using aircraft observations of ice particle shapes and sizes obtained in clouds behind midlatitude thunderstorms, this work develops a new framework for estimating ice particle fall velocities and their uncertainty, including quantifying the importance of different uncertainty sources from cloud microphysics measurements. Natural parameter variability contributes the most uncertainty in ice particle fall velocity estimates, although other sources can also be important contributors to uncertainty in certain conditions. Additional work examining ice particle data is needed to further understand how dependent uncertainty in certain ice particle properties are to local environmental conditions.

Meteorology & Atmospheric Sciences↗

DOE OSTI2026

Investigating the Effects of Bars on Star Formation and Nuclear Activity of Galaxies Using DESI Survey Data

We present a statistical analysis of the connections between galactic bars, star formation, and active galactic nucleus (AGN) activity using 33,201 disk galaxies (0.01 < z < 0.05) from Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) cross-matched with Galaxy Zoo DESI. Based on morphological classifications, we identify 3508 strongly barred and 8335 weakly barred systems. We find that barred galaxies exhibit a clear bimodal distribution in color–mass space: weak bars are preferentially found in bluer, lower-mass disks, whereas strong bars are more common in massive, redder systems. Strongly barred galaxies are on average more massive and metal rich than unbarred systems. In addition, strong bars enhance central star formation rates (SFRs) in low-mass galaxies but reduce specific SFRs in massive systems, reflecting a dual role where bars initially trigger central star formation but eventually promote quenching by accelerating gas consumption. In terms of nuclear activity, barred galaxies display a higher incidence of AGN activity. The presence of a bar is also associated with an increased fraction of powerful AGN, with the highest proportions found in strongly barred systems. However, the correlations between AGN activity and detailed bar structural parameters are weak, suggesting that the link between bars and nuclear activity is indirect and regulated by multiple factors. Overall, our results support a scenario in which bars facilitate angular momentum transport and gas inflow, thereby driving central star formation and fueling supermassive black hole accretion while operating alongside other processes that shape galaxy evolution.

Liu, Jianfei [Chinese Academy of Sciences (CAS), B↗

DOE OSTI2023

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel’dovich effect (kSZ). We precisely measure this effect (𝑆/𝑁 ≈ 10) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over ≥ 4,000 deg 2 . We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in 𝑔, 𝑟, and 𝑧 bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed (PTE)⁢ = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris 𝑧 = 0.5, PTE⁢ = 0.37) over weaker-feedback cases (IllustrisTNG 𝑧 = 0.8, PTE⁢ = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast 𝑆/𝑁 ∼ 50 for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. Furthermore, these measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Characterization and controllability of radiated power via extrinsic impurity seeding in strongly negative triangularity plasmas in DIII-D

Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗