SEARCH · Engineering Papers

Results for “global shape modeling”

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 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Global Magnetohydrodynamic Modeling of the Solar Corona

The coronal magnetic field defines the structure of the solar corona, the position of the heliospheric current sheet, the regions of fast and slow solar wind, and the most likely sites of coronal mass ejections. There are few measurements of the magnetic fields in the corona, but the line-of-sight component of the global magnetic fields in the photosphere have been routinely measured for many years (for example, at Stanford's Wilcox Solar Observatory, and at the National Solar Observatory at Kitt Peak). The SOI/MDI instrument is now providing high-resolution full-disk magnetograms several times a day. Understanding the large-scale structure of the solar corona and inner heliosphere requires accurately mapping the measured photospheric magnetic field into the corona and outward. Ideally, a model should not only extrapolate the magnetic field, but should self-consistently reconstruct both the plasma and magnetic fields in the corona and solar wind. Support from our NASA SR&T contract has allowed us to develop three-dimensional magnetohydrodynamic (MHD) computations of the solar corona that incorporate observed photospheric magnetic fields into the boundary conditions. These calculations not only describe the magnetic field in the corona and interplanetary spice, but also predict the plasma properties as well. Our computations thus far have been successful in reproducing many aspects of both coronal and interplanetary data, including the structure of the streamer belt, the location of coronal hole boundaries, and the position and shape of the heliospheric current sheet. The most widely used technique for extrapolating the photospheric magnetic field into the corona and heliosphere are potential field models, such as the potential field source-surface model (PFSS),and the potential field current-sheet (PFCS) model

Linker, Jon A.↗

NASA NTRS2023

Towards Soil Moisture Profile Estimation in the Root Zone Using L- and P-Band Radiometer Observations: A Coherent Modelling Approach

Precision irrigation management and crop water stress assessment rely on accurate estimation of root zone soil moisture. However, only the top 5cm soil moisture can be estimated using the two current passive microwave satellite missions, Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP), which operate at L-band (wavelength of ~21cm). Since the contributing depth of the soil to brightness temperature increases with observation wavelength, it is expected that a P-band (wavelength of ~40cm) radiometer could potentially provide soil moisture information from deeper layers of the soil profile. Moreover, by combining both L- and P- bands, it is hypothesized that the soil moisture profile can be estimated even beyond their individual observation depths. The aim of this study was to demonstrate the potential of combined L-band and P-band radiometer observations to estimate the soil moisture profile under flat bare soil using a stratified coherent forward model. Brightness temperature observations at L-band and P-band from a tower based experimental site across a dry (April 2019) and a wet (March 2020) period, covering different soil moisture profile shapes, were used in this study. Results from an initial synthetic study showed that the performance of a combined L-band and P-band approach was better than the performance of using either band individually, with an average depth over which reliable soil moisture profile information could be estimated (i.e. with a target root mean square error (RMSE) of less than 0.04 m3/m3) being 20cm for linear and 15cm for second-order polynomial functions. Other functions were also tested but found to have a poorer performance. Applying the method to the tower-based brightness temperature achieved an average estimation depth of 28cm (20cm) and 5cm (5cm) during the dry and wet periods respectively when using a second-order polynomial (linear) function. These findings highlight the opportunity of a satellite mission with L-band and P-band observations to accurately estimate the soil moisture profile to as deep as 30cm globally.

Foad Brakhasi↗

DOE OSTI2024

Data-driven Quasi-static Surrogate Models for Nuclear-powered Integrated Energy Systems

The integration of nuclear power into energy systems presents a promising avenue to address the growing global energy demands while minimizing greenhouse gas emissions. In this paper, we introduce a data-driven quasi-static surrogate model for nuclear-powered Integrated Energy Systems (IES) that comprises various components, including a small modular reactor (SMR), steam manifold, balance of plant (BOP), high-temperature steam electrolysis (HTSE), and district heating (DH) system. Traditional physics-based models for these components often entail significant computational resource and time consumption, necessitating the development of efficient surrogate models. The development of a complete surrogate model for the IES involves the creation of individual surrogate models for each component, leveraging machine learning techniques and simulated data. These isolated surrogate models are subsequently integrated, enabling a holistic view of the IES and reducing the computational burden associated with detailed physics-based simulations. This paper outlines the development process, validation, and the performance evaluation of the surrogate models. The exceptional performance, with low root-mean-squared errors and R-squared scores of at least 99.8% across all individual surrogate models, underscores their accuracy and practical applicability. These results demonstrate the potential of these models to expedite the analysis of nuclear-powered IES, offering insights that can shape future research and development efforts.

08 HYDROGEN↗

DOE OSTI2020

High-throughput mapping of the phage resistance landscape in E. coli

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2020

Object-Based Comparison of Data-Driven and Physics-Driven Satellite Estimates of Extreme Rainfall

The Global Precipitation Measurement (GPM) constellation of spaceborne sensors provides a variety of direct and indirect measurements of precipitation processes. Such observations can be employed to derive spatially and temporally consistent gridded precipitation estimates either via data-driven retrieval algorithms or by assimilation into physically based numerical weather models. We compare the data-driven Integrated Multisatellite Retrievals for GPM (IMERG) and the assimilation-enabled NASA-Unified Weather Research and Forecasting (NU-WRF) model against Stage IV reference precipitation for four major extreme rainfall events in the southeastern United States using an object-based analysis framework that decomposes gridded precipitation fields into storm objects. As an alternative to conventional ‘‘grid-by-grid analysis,’’ the object-based approach provides a promising way to diagnose spatial properties of storms, trace them through space and time, and connect their accuracy to storm types and input data sources. The evolution of two tropical cyclones are generally captured by IMERG and NU-WRF, while the less organized spatial patterns of two mesoscale convective systems pose challenges for both. NU-WRF rain rates are generally more accurate, while IMERG better captures storm location and shape. Both show higher skill in detecting large, intense storms compared to smaller, weaker storms. IMERG’s accuracy depends on the input microwave and infrared data sources; NU-WRF does not appear to exhibit this dependence. Findings highlight that an object-oriented view can provide deeper insights into satellite precipitation performance and that the satellite precipitation community should further explore the potential for ‘‘hybrid’’ data-driven and physics-driven estimates in order to make optimal usage of satellite observations.

extreme events↗

NASA NTRS2023

A Hyperspectral Inversion Framework for Estimating Absorbing Inherent Optical Properties and Biogeochemical Parameters in Inland and Coastal Waters

The simultaneous remote estimation of biogeochemical parameters (BPs) and inherent optical properties (IOPs) from hyperspectral satellite imagery of globally distributed optically distinct inland and coastal waters is a complex, unsolved, non-unique inverse problem. To tackle this problem, we leverage a machine-learning model termed Mixture Density Networks (MDNs). MDNs outperform operational algorithms by calculating the covariance between the simultaneously estimated products. We train the MDNs on a large ( N = 8237) dataset of co-aligned, in situ measured, hyperspectral remote sensing reflectance (R rs ), BPs, and absorbing IOPs from globally representative optically distinct inland and coastal waters. The estimated IOPs include absorption due to phytoplankton (a ph ), chromophoric dissolved organic matter (a cdom ), and non-algal particles (a nap ). The estimated BPs include chlorophyll-a, total suspended solids, and phycocyanin (PC). MDNs dramatically reduce uncertainty in the retrievals, relative to operational algorithms, when using a 50/50 dataset split, where the MDNs are trained on a randomly selected half of the in situ dataset and validated on the other half. Our model is shown to have higher, or equivalent, generalization performance than the calculated operational algorithms available for all BPs and IOPs (except PC) via a leave-one-out cross-validation assessment. The MDNs are sensitive to uncertainties in the hyperspectral satellite R rs , resulting from instrument noise and atmospheric correction; there is a difference of ~37.4–62.8% (using median symmetric accuracy) between the MDNs’ estimates derived from co-located satellite-derived R rs and in situ R rs . Of the IOPs, a cdom and a nap are less sensitive to uncertainties in hyperspectral satellite imagery relative to a ph , with remote estimates of a ph exhibiting incorrect spectral shape and magnitude relative to in situ measured IOPs. Despite the uncertainties in satellite derived R rs , the spatial distributions of BPs and IOPs in MDN-derived product maps of Lake Erie and the Curonian Lagoon, based on imagery taken with the Hyperspectral Imager for the Coastal Ocean (HICO) and PRecursore Iper-Spettrale della Missione Applicativa (PRISMA), are confirmed via co-aligned in situ measurements and agree with the literature’s understanding of these well-studied regions. The consistency and accuracy of the model on HICO and PRISMA imagery, despite radiometric uncertainties, demonstrate its applicability to future hyperspectral missions, such as the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission, where the simultaneous estimation model will serve as a key part of phytoplankton community composition analysis.

Ryan E. O'Shea↗

NASA NTRS2008

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

NASA NTRS2012

Airborne Lidar Measurements of Atmospheric Pressure Made Using the Oxygen A-Band

Accurate measurements of greenhouse gas mixing ratios on a global scale are currently needed to gain a better understanding of climate change and its possible impact on our planet. In order to remotely measure greenhouse gas concentrations in the atmosphere with regard to dry air, the air number density in the atmosphere is also needed in deriving the greenhouse gas concentrations. Since oxygen is stable and uniformly mixed in the atmosphere at 20.95%, the measurement of an oxygen absorption in the atmosphere can be used to infer the dry air density and used to calculate the dry air mixing ratio of a greenhouse gas, such as carbon dioxide or methane. OUT technique of measuring Oxygen uses integrated path differential absorption (IPDA) with an Erbium Doped Fiber Amplifier (EDF A) laser system and single photon counting module (SPCM). It measures the absorbance of several on- and off-line wavelengths tuned to an O2 absorption line in the A-band at 764.7 nm. The choice of wavelengths allows us to maximize the pressure sensitivity using the trough between two absorptions in the Oxygen A-band. Our retrieval algorithm uses ancillary meteorological and aircraft altitude information to fit the experimentally obtained lidar O2 line shapes to a model atmosphere and derives the pressure from the profiles of the two lines. We have demonstrated O2 measurements from the ground and from an airborne platform. In this paper we will report on our airborne measurements during our 2011 campaign for the ASCENDS program.

Riris, Haris↗

NASA NTRS2019

Informing Algorithms from Ground Validation: The GPM Combined Algorithm

Ground validation (GV) for the Global Precipitation Measurement (GPM) mission encompasses in-situ (e.g., gauge, disdrometer) measurements, ground radar products, and comprehensive datasets from dedicated airborne field campaigns. These datasets are used for direct validation of the precipitation products from GPM as well as to inform assumptions used by the algorithms that produce these products. This presentation will focus on current and potential such uses of GV data in the GPM combined radar-radiometer algorithm.The GPM combined algorithm, by virtue of using data from the Dual-frequency Precipitation Radar (DPR) and GPM Microwave Imager (GMI), is the most well-constrained instantaneous precipitation product from GPM. It also plays an important in role the passive microwave algorithms as a basis for the construction of brightness temperature-precipitation profile databases and is a calibrator for the multi- satellite gridded product IMERG. Since both radar and radiometer data are used as observational inputs, and even with these data the retrievals are underconstrained, microphysical properties of the hydrometeor profile that are relevant over the range of GMI and DPR wavelengths and incidence angles need to be assumed by the algorithm forward models. Such properties include the normalized intercept parameter (N w ) and shape parameter () of the particle size distribution, ice size-density relationship, and particle size-aspect ratio relationships. The sub-beam variability must also be prescribed in order to accurately simulate the observed radar reflectivity profiles without introducing significant biases. The sensitivity of the algorithm to these parameters, along with their mean values and variability, will be discussed. Examples of future directions include refinement with new data (particularly the assumptions related to ice), and seeking relationships between assumptions and observable or environmental data.

Munchak, Joe.↗

DOE OSTI2020

Solar Neutrino Observations at the Sudbury Neutrino Observatory

The Sudbury Neutrino Observatory (SNO) is a 1000-tonne heavy water Cherenkov detector. Its usage of \dto as target allows the simultaneous measurements of the $\nu_e$ flux from $^8$B decay in the Sun and the total flux of all active neutrino species through the charged-current and the neutral-current interactions on the deuterons. Assuming the standard $^{8}$B shape, the $\nu_e$ component of the $^{8}$B solar neutrino flux is measured to be $\phi_e = 1.76^{+0.05}_{-0.05}{(stat.)}^{+0.09}_{-0.09} {(syst.)} x 10^6 {\rm cm}^{-2} {\rm s}^{-1}$ for a kinetic energy threshold of 5 MeV. The non-\nue component is found to be Φ μτ $\= 3.41^{+0.45}_{-0.45}{(stat.)}^{+0.48}_{-0.45} {(syst.)} x 10^6 {\rm cm}^{-2} {\rm s}^{-1}$. This $5.3\sigma$ difference provides strong evidence for $\nu_{e}$ flavor transformation in the solar neutrino sector. The total active neutrino flux is measured with the neutral-current reaction at a neutrino energy threshold of 2.2 MeV. This flux is determined to be $\phi_{NC} = 5.09^{+0.44}_{-0.43}{(stat.)}^{+0.46}_{-0.43} {(syst.)} x 10^6 {\rm cm}^{-2} {\rm s}^{-1}$, and is consistent with solar model predictions. Assuming an undistorted $^8$B spectrum, the night minus day rate is 14.0$\pm$6.3(stat.)$^{+1.5}_{-1.4}$(sys.)\% of the average rate in the charged-current channel. If the total active neutrino flux is constrained to have no asymmetry, the night-day asymmetry in the $\nu_e$ flux is found to be 7.0$\pm$4.9(stat.)$^{+1.3}_{1.2}$(sys.)\%. A global analysis of all the available solar neutrino data in terms of matter-enhanced oscillations of two active flavors strongly favors the Large Mixing Angle (LMA) solution.

Poon, A. W. P.↗

DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗

NASA NTRS1996

Transient Creep of a Composite Lower Crust: A Polymineralic Basis for Rapidly Evolving Postseismic Deformation Modes - 2

Postseismic horizontal strain and displacement following the June 28, 1992, Landers, California, earthquake (M(sub W) 7.3) is broad scale and cannot be explained solely by delayed afterslip located at the rupturing fault trace. Both the observed strain at Pifion Flat Observatory (PFO) and observed Global Positioning System receiver velocities evolve rapidly after the Landers-Big Bear earthquake sequence. The observed exponential decay of these motions, with timescales of 4-34 days, may reflect a soft creep rheology in the lower crust and brittle-ductile transition zone or even within the seismogenic crust itself. Here a simple model of a two-dimensional screw dislocation in a layered Maxwell viscoelastic Earth is employed in conjunction with a composite rheology to demonstrate that the short timescale transient response modes (approx. = 4-34 days) are consistent with the behavior of a biviscous lower crust. The lowest viscosity of this system is derivable from laboratory experimental data on the long-term creep of natural quarztites, and the highest viscosity is consistent with isostasy-related lower crustal flow in a continental extensional tectonic environment. The model predicts significant stress relaxation at the base of the seismogenic crust. Near the base of the seismogenic zone, and about 4 km away from the mainshock, the rate of predicted relaxation is of the order of 0.01 MPa/ d during the first 20 days of postseismic flow. Oblate spheroidal inclusions at 5% concentration levels that are both aligned and fairly flat in shape and that have a viscosity of 3-4 x 10(exp 15) Pa s are consistent with both the amplitude and decay time of horizontal crustal strain observed at PFO after the Landers mainshock. It is speculated that the structures exposed in cross sections and in seismic reflection profiles of the lower crust that have mylonitic associations are, in part, the cause of such rapid postseismic evolution in southeastern California. Unmylonitized quartz-rich rock at sufficiently elevated temperatures could also contribute to the rapid decay modes.

Ivins, Erik R.↗

NASA NTRS2011

Geochemical Mapping of 4 Vesta Begins

By December, the NASA Dawn spacecraft will have descended to a low altitude mapping orbit (LAMO), where the Gamma Ray and Neutron Detector (GRaND) will acquire global mapping data for up to four months. Measurements by GRaND will help answer elusive questions about how Vesta differentiated and the nature of processes that shaped Vesta s surface. The data will be analyzed to determine the abundances of Mg, Si, Fe, K, Th, and H at a spatial resolution of roughly 300 km full-width-at-half-maximum from a 465 km radius orbit. Thermal and fast neutron counting data will be analyzed to determine the neutron macroscopic absorption cross section and average atomic mass, providing constraints on additional elements, such as Ca and Al. GRaND will quantify the elemental composition of coarse spatial units identified by Dawn s Framing Camera (FC) and the Visible & Infrared Spectrometer (VIR). In addition, GRaND will map the mixing ratio of compositional end members selected from the howardite, eucrite and diogenite (HED) meteorites, determine the relative proportions of plagioclase and mafic minerals, and search for compositions that are absent or under-represented in the meteorite collection. While it is generally thought that Vesta s crust on a regional scale should be well-represented by linear mixing of HED whole-rock compositions, there are hints that Vesta may be more diverse than implied by this model. For example, the discovery of K-rich impact glasses in howardites suggests that K-rich rocks may be present on a portion of Vesta s surface, and the analysis of diogenites indicates considerable variability in the magmatic processes that formed them. The chemical composition of materials within Vesta s south polar structure may provide further clues to how it formed. An impact might have exposed mantle and lower crustal materials, which should have a distinctive compositional signature. We present the analysis of data acquired by GRaND from cruise through the descent to LAMO, including GRaND s sensitivity to different elements and geochemical processes.

Prettyman, Thomas H.↗

DOE OSTI2025

Online and Offline Data Quality Monitoring for the Mu2e Calorimeter

This thesis presents the design, implementation, and validation of a calorimeter Data Quality Monitoring (DQM) toolchain for the Mu2e experiment at Fermilab. Mu2e searches for charged lepton flavor violation via coherent muon-to-electron conversion in the field of an aluminum nucleus, $\mu^- Al \rightarrow e^-Al$, a process whose observation would constitute clear evidence of physics beyond the Standard Model. Achieving target sensitivity requires stringent control of detector performance and data integrity during acquisition, as subtle issues in readout configuration, data formatting, or electronics behavior can compromise reconstruction and bias downstream analyzes. To address these challenges, this work develops a multi-layer DQM approach spanning both raw data validation and reconstructed digi-level diagnostics. At the low level, a fragment analysis component performs word- and bit-field decoding of calorimeter readout blocks, enabling sanity checks of the expected structure and producing detailed error and integrity statistics useful for commissioning and troubleshooting. At the digi level, the CaloDigiDQM analyzer is implemented within the art framework and transforms each CaloDigiCollection into a structured hierarchy of ROOT histograms designed for fast drill-down diagnostics. The module generates coherent monitoring views at global, disk, board, and channel granularity, including occupancy, waveform-derived features (baseline, RMS, peak amplitude and position), and left-right sensor consistency metrics. Detector-aware channel-to-electronics mapping is performed through the conditions system (CaloDAQMap), ensuring that diagnostics remain aligned with hardware identifiers used in operations. For end-to-end testing without reliance on live DAQ data, a synthetic CaloDigi producer is developed to generate realistic waveforms with controlled noise and pulse shapes. The resulting system supports both offline ROOT-file production and online operation, including optional histogram streaming through otsdaq via ots::HistoSender. This toolchain provides a practical and scalable foundation for calorimeter commissioning and stable data collection, enabling early detection of anomalies and reducing operational risk for Mu2e.

Vakulenko, Mark [Drew U.] (ORCID:0009000276197818)↗

DOE OSTI2021

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗

DOE OSTI2020

Efficient method for estimation of fission fragment yields of $\textit{r}$-process nuclei

Background: More than half of all the elements heavier than iron are made by the rapid neutron capture process (or $\textit{r}$ process). For very-neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the $\textit{r}$-process endpoint, and the fission-fragment yields shape the final abundances of 110 ≤ $\textit{A}$ ≤ 170 nuclei. The knowledge of fission-fragment yields of hundreds of nuclei inhabiting very-neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. Purpose: In this study, we propose a model for the fast calculation of fission-fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. Methods: To generate realistic potential-energy surfaces and nucleonic localizations, we apply Skyrme density-functional theory. In this work, the distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. Results: We benchmark the method by studying the fission yields of 178 Pt, 240 Pu, 254 Cf, and 254,256,258 Fm and show that it satisfactorily explains the experimental data. We then make predictions for 254 Pu and 290 Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the $\textit{r}$-process nucleosynthesis occurring in neutron-star mergers. Conclusions: The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission-fragment distributions across the $\textit{r}$-process region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗