SEARCH · Engineering Papers

Results for “constrained optimization”

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 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Skewing the CMB×LSS: a fast method for bispectrum analysis

Upcoming cosmic microwave background (CMB) lensing measurements and tomographic galaxy surveys are expected to provide us with high-precision data sets in the coming years, thus paving the way for fruitful cross-correlation analyses. Here, in this paper we study the information content of the weighted skew-spectrum, a nearly-optimal estimator of the angular bispectrum amplitude, as a means to extract non-Gaussian information on both bias and cosmological parameters from the bispectra of galaxies cross-correlated with CMB lensing, while gaining significantly on speed. Our results show that for the combination of the Planck satellite and the Dark Energy Spectroscopic Instrument (DESI), the difference in the constraints on bias and cosmological parameters from the skew-spectrum and the bispectrum is at most 17%. We further compare and find agreement between our theoretical skew-spectra and those estimated from N-body simulations, for which it is important to include gravitational non-linearities beyond perturbation theory and the post-Born effect for CMB lensing. We define an algorithm to apply the skew-spectrum estimator to the data and, as a preliminary step, we use the skew-spectra to constrain bias parameters and the amplitude of shot noise from the simulations through a Markov chain Monte Carlo likelihood analysis, finding that it may be possible to reach percent-level estimates for the linear bias parameter b 1 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗

DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗

DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

Views from EPOXI: Colors in Our Solar System as an Analog for Extrasolar Planets

The first visible-light studies of Earth-sized extrasolar planets will employ photometry or low-resolution spectroscopy. This work uses EPOCh medium-hand filter photometry between 150 and 950 nm obtained with the Deep Impact (DI) High Resolution Instrument (HRI) of Earth, the Moon, and Mars in addition to previous full-disk observations of the other six solar system planets and Titan to analyze the limitations of using photometric colors to characterize extrasolar planets. We determined that the HRI 350, 550, and 850 nm filters are optimal for distinguishing Earth from the other planets and separating planets to first order based on their atmospheric and surface properties. Detailed conclusions that can be drawn about exoplanet atmospheres simply from a color-color plot are limited due to potentially competing physical processes in the atmosphere. The presence of a Rayleigh scattering atmosphere can be detected by an increase in the 350-550 nm brightness ratio, but the absence of Rayleigh scattering cannot be confirmed due to the existence of atmospheric and surface absorbing species in the UV. Methane and ammonia are the only species responsible for strong absorption in the 850 nm filter in our solar system. The combination of physical processes present on extrasolar planets may differ from those we see locally. Nevertheless, a generation of telescopes capable of collecting such photometric observations can serve a critical role in first-order characterization and constraining the population of Earth-like extrasolar planets.

Crow, Carolyn A.↗

NASA NTRS2022

Solar Polar Flux Redistribution Based on Observed Coronal Holes

We explore the use of observed polar coronal holes (CHs) to constrain the flux distribution within the polar regions of global solar magnetic field maps in the absence of reliable quality polar field observations. Global magnetic maps, generated by the Air Force Data Assimilative Photospheric flux Transport (ADAPT) model, are modified to enforce field unipolarity thresholds both within and outside observed CH boundaries. The polar modified and unmodified maps are used to drive Wang–Sheeley–Arge (WSA) models of the corona and solar wind (SW). The WSA-predicted CHs are compared with the observations, and SW predictions at the WIND and Ulysses spacecraft are also used to provide context for the new polar modified maps. We find that modifications of the polar flux never worsen and typically improve both the CH and SW predictions. We also confirm the importance of the choice of the domain over which WSA generates the coronal magnetic field solution but find that solutions optimized for one location in the heliosphere can worsen predictions at other locations. Finally, we investigate the importance of low-latitude (i.e., active region) magnetic fields in setting the boundary of polar CHs, determining that they have at least as much impact as the polar fields themselves.

Solar coronal holes↗

DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmological constraints from galaxy clustering and galaxy-galaxy lensing using the MagLim lens sample

The cosmological information extracted from photometric surveys is most robust when multiple probes of the large scale structure of the Universe are used. Two of the most sensitive probes are the clustering of galaxies and the tangential shear of background galaxy shapes produced by those foreground galaxies, so-called galaxy-galaxy lensing. Combining the measurements of these two two-point functions leads to cosmological constraints that are independent of the way galaxies trace matter (the galaxy bias factor). The optimal choice of foreground, or lens, galaxies is governed by the joint, but conflicting requirements to obtain accurate redshift information and large statistics. We present cosmological results from the full 5000 deg 2 of the Dark Energy Survey’s first three years of observations (Y3) combining those two-point functions, using for the first time a magnitude-limited lens sample (MagLim) of 11 million galaxies, especially selected to optimize such combination, and 100 million background shapes. We consider two flat cosmological models, the Standard Model with dark energy and cold dark matter ( Λ CDM ) a variation with a free parameter for the dark energy equation of state ( w CDM ). Both models are marginalized over 25 astrophysical and systematic nuisance parameters. In Λ CDM we obtain for the matter density Ω m = 0.320 - 0.034 + 0.041 and for the clustering amplitude S 8 ≡ σ 8 ( Ω m / 0.3 ) 0.5 = 0.778 - 0.031 + 0.037 , at 68% C.L. The latter is only 1 σ smaller than the prediction in this model informed by measurements of the cosmic microwave background by the Planck satellite. In w CDM we find Ω m = 0.32 - 0.046 + 0.044 , S 8 = 0.777 - 0.051 + 0.049 and dark energy equation of state w = - 1.031 - 0.379 + 0.218 . We find that including smaller scales, while marginalizing over nonlinear galaxy bias, improves the constraining power in the Ω m - S 8 plane by 31% and in the Ω m - w plane by 41% while yielding consistent cosmological parameters from those in the linear bias case. These results are combined with those from cosmic shear in a companion paper to present full DES-Y3 constraints from the three two-point functions ( 3 × 2 pt ).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

A Unified Beam-Dynamics and Hardware Design Framework for Hybrid Nonlinear-Kicker Injection in NSLS-II

Nonlinear kickers (NLKs) enable off-axis injection in ultralow-emittance storage rings by providing a strong kick to the injected beam while remaining nearly transparent to the stored beam. In hybrid schemes, a conventional four-kicker bump defines the injected trajectory, and the NLK reduces the first-turn action under constrained beam offset and optics conditions. Effective operation additionally requires stable and reproducible first-turn injection trajectories. We develop a compact action–angle framework that expresses NLK dynamics in terms of Courant–Snyder invariants and yields an analytical bound on achievable action reduction. This formulation provides direct design rules for NLK placement, phase advance, injected-beam offset, and kicker field profile. Within this framework, we identify the 8-wire NLK as a practical baseline and extend its design by relaxing the square-geometry constraint, enabling inward shifting of the off-axis field peak while preserving on-axis field and gradient cancellation. Application to the NSLS-II lattice shows how aperture, pulsed-power, and mechanical constraints combine to determine a coupled design solution. Multi-turn tracking confirms that candidate NLK locations maintain sufficient stay-clear (aperture-clearance) margin, while the optimized wire geometry reduces the required current and Lorentz force load. The results establish a unified approach for NLK-assisted injection design and provide a practical pathway for upgrades in diffraction-limited storage rings.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Low‐Temperature and High‐Rate Rechargeable Aluminum Batteries Enabled by Ternary Eutectic Electrolytes

Abstract Rechargeable aluminum batteries (RABs) have garnered extensive scientific attention as a promising alternative chemistry due to the inherent advantages associated with aluminum (Al) metal anodes, including their high theoretical capacities, cost‐effectiveness, environmental friendliness, and inherent non‐flammable properties. Nonetheless, the practical energy density of RABs is constrained by the electrolytes that support lower operational voltage windows. Herein, we report a ternary eutectic electrolyte composed of 1‐ethyl‐3‐methylimidazolium chloride ([C 2 C 1 im]Cl):1‐butyl‐3‐methylimidazolium chloride ([C 4 C 1 im]Cl):aluminum chloride (AlCl 3 ) for the application of RABs. The electrolyte exhibits a high operational potential window (~3 V vs. Al/Al 3+ on SS 316) and high ionic conductivity (~8.3 mS cm −1 ) while exhibiting only a low temperature glass transition at −65 °C suitable for all‐climate conditions. Al||graphene nanoplatelets cell delivers a high capacity of ~117 mAh/g, and ~43 mAh/g at a very high current densities of 1 A/g and 5 A/g, respectively. The cells render a reversible capacity of 20 mAh/g at −20 °C and 17 mAh/g at −40 °C, indicating their suitability for operation under extreme environmental conditions. We comprehensively evaluated the design and optimization of carbon paper‐based battery systems. The ternary eutectic electrolyte demonstrates exceptional electrochemical performance, thus signifying its substantial potential for utilization in high‐performance energy storage systems in all climates.

Vadthya, Raju↗

DOE OSTI2026

The solubility of NdPO 4 and DyPO 4 and stability of Nd and Dy chloride and hydroxyl complexes as a function of pH and salinity to 450 °C

The mobility of rare earth elements (REE) in geological systems is often controlled by the stability of monazite and xenotime, which are important hosts for the light and heavy REE. These REE phosphates constitute important resources and provide information on ore formation conditions and timing due to their uses as geothermometers and geochronometers. While the thermodynamic properties for these minerals are well-established up to high temperature and pressure, the properties for the aqueous REE complexes are not well constrained with limited information up to 300 °C. In this study, the solubility of synthetic NdPO 4 and DyPO 4 endmembers were measured in hydrothermal sub- to supercritical NaCl-bearing aqueous solutions from 350 at saturated water vapor pressure to 450 °C at 700 bar. The speciation of Nd and Dy was determined in acidic to alkaline solutions (pH 25 °C values from 2 to 10), and indicates that the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 predominate at low salinity (0.01 mol/kg NaCl) with some contribution of REE chloride species REECl 2+ and REECl 2 + in acidic fluids. The REE phosphate solubility measured in these experiments is up to two orders of magnitude higher than predicted using existing thermodynamic properties from literature for aqueous species extrapolated from lower temperature data to supercritical conditions. To address this discrepancy, the thermodynamic properties of the REE aqueous species were optimized using GEMSFITS to derive the formation constants for the REE chloride (β Cl ) and hydroxyl (β OH ) complexes in the studied temperature range. This study highlights the importance of the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 over a wide range of pH and temperature. As a result, the revised thermodynamic properties for Dy and Nd chloride and hydroxyl species more accurately predict the solubility of NdPO 4 and DyPO 4 and provide insight into the speciation and mobility of the light and heavy REE in hydrothermal supercritical fluids in the crust.

58 GEOSCIENCES↗

DOE OSTI2024

Super Resolving Unrolled Neural Networks for Remote Sensing

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Opportunities for Treatment and Reuse of Agricultural Drainage in the United States

Irrigation accounts for 42% of the total freshwater withdrawals in the United States. Climate change, the pressure of a growing population, degrading water quality, and increased competition from other sectors could constrain continuous supply to meet future agricultural water demand. This study presents an evaluation framework to assess the potential reuse of agricultural drainage water for crop irrigation. Using a regional approach, we review the current state of agricultural drainage treatment and reuse and the institutional, economic, and other barriers that can influence the reuse decision. In the 31 eastern states, agricultural drainage contains valuable nutrients that can be reused for irrigation with minimal treatment, while the 17 western states struggle with large volumes of saline drainage that can contain constituents of concern (e.g., selenium), preventing reuse without treatment. Using a new decision-support tool called WaterTAP3, a potential treatment train for saline agricultural drainage was analyzed to identify treatment challenges, research needs, and the potential implementation at a larger scale. As demonstrated by our case study, desalination of agricultural drainage is costly and energy intensive and will require sizable investments to fully develop and optimize technologies as well as manage the generated waste and brine.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Gravitational waves and galaxies cross-correlations: a forecast on GW biases for future detectors

ABSTRACT Gravitational waves (GWs) have rapidly become important cosmological probes since their first detection in 2015. As the number of detected events continues to rise, upcoming instruments like Einstein Telescope (ET) and Cosmic Explorer (CE) will observe millions of compact binary mergers. These detections, coupled with galaxy surveys by instruments such as the Dark Spectroscopic Energy Instrument (DESI), Euclid, and the Vera Rubin Observatory, will provide unique information on the large-scale structure of the universe by cross-correlating GWs with the distribution of galaxies hosting them. In this paper, we focus on how cross-correlations constrain the clustering bias of GWs emitted by the coalescence of binary black holes (BBHs). This parameter links BBHs to the underlying dark matter distribution, hence informing us how they populate galaxies. Using a multitracer approach, we forecast the precision of these measurements under different survey combinations. Our results indicate that current GW detectors will have limited precision, with measurement errors as high as $\displaystyle \sim 50~{{\ \rm per\ cent}}$. However, third-generation detectors like ET, when cross-correlated with Legacy Survey of Space and Time (LSST) data, can improve clustering bias measurements to within 2.5 per cent. Furthermore, we demonstrate that these cross-correlations can enable a per cent-level measurement of the magnification lensing effect on GWs. Despite this, there is a degeneracy between magnification and evolution biases, which hinders the precision of both. This degeneracy is most effectively addressed by assuming knowledge of one bias or targeting an optimal redshift range of $\displaystyle 1 \lt z \lt 2.5$. Our analysis opens new avenues for studying the distribution of BBHs and testing the nature of gravity through large-scale structure.

Zazzera, Stefano (ORCID:0000000158979221)↗

DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗

DOE OSTI2021

Spatial biases of information influence global estimates of soil respiration: how can we improve global predictions?

Soil respiration (Rs), the efflux of CO2 from soils to the atmosphere, is a major component of the terrestrial carbon cycle, but is poorly constrained from regional to global scales. The global soil respiration database (SRDB) is a compilation of in-situ Rs observations from around the globe that has been consistently updated with new measurements over the past decade. It is unclear whether the addition of data to new versions of the SRDB has produced better-constrained global Rs estimates. We compared two versions of the SRDB (v3.0 n=5173 and v5.0 n=10366) to determine how additional data influenced the global sum and uncertainty of Rs as predicted by a quantile regression forest model. A model parameterized using SRDBv3 yielded a global Rs sum of 88.6 Pg C yr-1, and associated uncertainty of 29.9 (mean absolute error) and 57.9 (standard deviation) Pg C yr-1, whereas parameterization using SRDBv5 yielded 96.5 Pg C yr-1 and associated uncertainty of 30.2 (mean average error) and 73.4 (standard deviation) Pg C yr-1. Empirically estimated global heterotrophic respiration (Rh) from v3 and v5 were 49.9-50.2 (mean 50.1) and 53.3-53.5 (mean 53.4) Pg C yr-1, respectively. Interestingly, SRDBv5’s inclusion of new data from underrepresented regions (e.g., Asia, Africa, South America) thus resulted in overall higher model uncertainty. The SRDBv5 is still biased towards northern latitudes and temperate zones, so we tested that an optimized global distribution of Rs measurements resulted in a global sum of 96.4 ± 21.4 Pg C yr-1 and lower model uncertainty. These results support current global estimates of Rs but highlight spatial inconsistencies that influence model parameterization and interpretation that should be addressed to improve regional- to global-scale predictions.

Stell, Emma↗

NASA NTRS2012

Assessment of Mars Atmospheric Temperature Retrievals from the Thermal Emission Spectrometer Radiances

Motivated by the needs of Mars data assimilation. particularly quantification of measurement errors and generation of averaging kernels. we have evaluated atmospheric temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) radiances. Multiple sets of retrievals have been considered in this study; (1) retrievals available from the Planetary Data System (PDS), (2) retrievals based on variants of the retrieval algorithm used to generate the PDS retrievals, and (3) retrievals produced using the Mars 1-Dimensional Retrieval (M1R) algorithm based on the Optimal Spectral Sampling (OSS ) forward model. The retrieved temperature profiles are compared to the MGS Radio Science (RS) temperature profiles. For the samples tested, the M1R temperature profiles can be made to agree within 2 K with the RS temperature profiles, but only after tuning the prior and error statistics. Use of a global prior that does not take into account the seasonal dependence leads errors of up 6 K. In polar samples. errors relative to the RS temperature profiles are even larger. In these samples, the PDS temperature profiles also exhibit a poor fit with RS temperatures. This fit is worse than reported in previous studies, indicating that the lack of fit is due to a bias correction to TES radiances implemented after 2004. To explain the differences between the PDS and Ml R temperatures, the algorithms are compared directly, with the OSS forward model inserted into the PDS algorithm. Factors such as the filtering parameter, the use of linear versus nonlinear constrained inversion, and the choice of the forward model, are found to contribute heavily to the differences in the temperature profiles retrieved in the polar regions, resulting in uncertainties of up to 6 K. Even outside the poles, changes in the a priori statistics result in different profile shapes which all fit the radiances within the specified error. The importance of the a priori statistics prevents reliable global retrievals based a single a priori and strongly implies that a robust science analysis must instead rely on retrievals employing localized a priori information, for example from an ensemble based data assimilation system such as the Local Ensemble Transform Kalman Filter (LETKF).

Hoffman, Matthew J.↗

DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nanotomography for Quantitative 3D Particle Reconstruction

Particulates are ubiquitous across fuel cycle operations and carry critical information about particle formation, processing, and potential proliferation-related activities. Traditional analytical techniques, including micro-Raman spectroscopy and standard electron microscopy, are often limited in spatial resolution or dimensionality, particularly when used to examine metallic or submicron-scale features. Understanding particle morphology, phase distribution, and internal porosity is essential for constraining formation conditions, thermodynamic environments, and material transport behavior. In this report, we demonstrate the application of plasma focused ion beam nanotomography to reconstruct micron-scale particulates at nanoscale resolution. Using high-resolution backscattered electron imaging and Avizo software, we obtained 3D reconstructions that enabled quantitative analysis of particle morphology, phase composition, and internal voids. Representative examples include a Ta particle with a large central void and a composite particle with embedded tetrahedral crystalline structures. These reconstructions reveal structural and compositional details that are inaccessible through conventional 2D imaging. The results demonstrate that nanotomography provides both qualitative and quantitative insights into particle formation and behavior. Using nanotomography, porosity and phase distributions can be quantified to inform models of particle density, transport, and solidification conditions. Beyond technical insights, the workflow developed here establishes a transferable capability for analyzing heterogeneous particles and has potential applications in bulk materials studies via x-ray computed tomography or other volumetric imaging modalities. Ongoing efforts are focused on optimizing the workflow to process multiple particles simultaneously, increasing throughput and statistical robustness. Overall, this work illustrates the power of nanotomography as a tool for connecting particulate morphology to formation mechanisms, composition, and transport, thereby strengthening analytical capabilities for nuclear forensics, fuel cycle analysis, and related scientific investigations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗