Search NASA⌕ Search

SEARCH · Search NASA

Results for “structural dynamics testing”

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

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

At least 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Single particle impact explorations on the role of powder heat treatments in cold spray: Effects of strength and oxide structure

In cold spray, the optimization of the powder feedstock for improving the final properties of the coating and for improving the process itself are two different objectives. This paper explores the latter, studying the relationship between powder pre-heat treatment and the kinetic bonding between the individual particles and the substrate. Using a laser-induced particle impact test (LIPIT), we investigate the effect of different heat treatments on the critical adhesion velocity of Al 6061. Although heat treatments change the strength of the powders as well as the thickness of the native oxide layer each by as much as ∼50 %, the critical adhesion velocity remains constant to within only ∼10 %. A deeper analysis reveals that these two factors change in somewhat complementary ways for the heat treatments explored here, leading to counteracting effects. By quantitatively evaluating these two microstructural factors in tandem, we find that they are roughly equally important in determining the critical adhesion velocity. These observations thus provide insight for the development of new heat treatment regimens that might lower the adhesion velocity.

Al 6061↗
DOE OSTI2015

Hydrogen-bond Specific Materials Modification in Group IV Semiconductors

Semiconductor dielectric crystals consist of two fundamental components: lattice atoms and electrons. The former component provides a crystalline structure that can be disrupted by various defects or the presence of an interface, or by transient oscillations known as phonons. The latter component produces an energetic structure that is responsible for the optical and electronic properties of the material, and can be perturbed by lattice defects or by photo-excitation. Over the period of this project, August 15, 1999 to March 31, 2015, a persistent theme has been the elucidation of the fundamental role of defects arising from the presence of radiation damage, impurities (in particular, hydrogen), localized strain or some combination of all three. As our research effort developed and evolved, we have experienced a few title changes, which reflected this evolution. Throughout the project, ultrafast lasers usually in a pump-probe configuration provided the ideal means to perturb and study semiconductor crystals by both forms of excitation, vibrational (phonon) and electronic (photon). Moreover, we have found in the course of this research that there are many interesting and relevant scientific questions that may be explored when phonon and photon excitations are controlled separately. Our early goals were to explore the dynamics of bond-selective vibrational excitation of hydrogen from point defects and impurities in crystalline and amorphous solids, initiating an investigation into the behavior of hydrogen isotopes utilizing a variety of ultrafast characterization techniques, principally transient bleaching spectroscopy to experimentally obtain vibrational lifetimes. The initiative could be divided into three related areas: (a) investigation of the change in electronic structure of solids due to the presence of hydrogen defect centers, (b) dynamical studies of hydrogen in materials and (c) characterization and stability of metastable hydrogen impurity states under transient compression. This research focused on the characterization of photon and ion stimulated hydrogen related defect and impurity reactions and migration in solid state matter, which requires a detailed understanding of the rates and pathways of vibrational energy flow, of the transfer channels and of the coupling mechanisms between local vibrational modes (LVMs) and phonon bath as well as the electronic system of the host material. It should be stressed that researchers at Vanderbilt and William and Mary represented a unique group with a research focus and capabilities for low temperature creation and investigation of such material systems. Later in the program, we carried out a vigorous research effort addressing the roles of defects, interfaces, and dopants on the optical and electronic characteristics of semiconductor crystals, using phonon generation by means of ultrafast coherent acoustic phonon (CAP) spectroscopy, nonlinear characterization using second harmonic generation (SHG), and ultrafast pump-and-probe reflectivity and absorption measurements. This program featured research efforts from hydrogen defects in silicon alone to other forms of defects such as interfaces and dopant layers, as well as other important semiconducting systems. Even so, the emphasis remains on phenomena and processes far from equilibrium, such as hot electron effects and travelling localized phonon waves. This program relates directly to the mission of the Department of Energy. Knowledge of the rates and pathways of vibrational energy flow in condensed matter is critical for understanding dynamical processes in solids including electronically, optically and thermally stimulated defect and impurity reactions and migration. The ability to directly probe these pathways and rates allows tests of theory and scaling laws at new levels of precision. Hydrogen embedded in model crystalline semiconductors and metal oxides is of particular interest, since the associated local mode can be excited cleanly, and is usually well-separated in energy from the phonon bath. These basic dynamical studies have provided new insights for example into the fundamental mechanisms that control proton diffusion in these oxides. This area of materials science has largely fulfilled its promise to identify degradation mechanisms in electronic and optoelectronic devices, and to advance solid oxide proton conductors for fuel cells, gas sensors and proton-exchange membrane applications. It also provides the basis for innovations in materials synthesis involving atomic-selective diffusion and desorption.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting 3D Motions from Single-Camera Optical Test Data

In a typical optical test, a stereo camera pair is required to measure the three-dimensional motion of a test article; one camera typically only measures motions in the image plane of the camera, and measurements in the out-of-plane direction are missing. Finite element expansion techniques provide a path to estimate responses from a test at unmeasured degrees of freedom. Treating the case of a single camera as a measurement with unmeasured degrees of freedom, a finite element model is used to expand to the missing third dimension of the image data, allowing a full-field, three-dimensional measurement to be obtained from a set of images from a single camera. The key to this technique relies on the mapping of finite element deformations to image deformations, creating a set of mode shape images that are used to filter the response in the image into modal responses. These modal responses are then applied to the finite element model to estimate physical responses at all finite element model degrees of freedom. The mapping from finite element model to image is achieved using synthetic images produced by a rendering software. Finally, the technique is applied first to a synthetic deformation image, and then is validated using an experimental set of images.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Sepsis reconsidered: Identifying novel metrics for behavioral landscape characterization with a high-performance computing implementation of an agent-based model

Objectives: Sepsis affects nearly 1 million people in the United States per year, has a mortality rate of 28–50% and requires more than $20 billion a year in hospital costs. Over a quarter century of research has not yielded a single reliable diagnostic test or a directed therapeutic agent for sepsis. Central to this insufficiency is the fact that sepsis remains a clinical/physiological diagnosis representing a multitude of molecularly heterogeneous pathological trajectories. Advances in computational capabilities offered by High Peformance Computing (HPC) platforms call for an evolution in the investigation of sepsis to attempt to define the boundaries of traditional research (bench, clinical and computational) through the use of computational proxy models. We present a novel investigatory and analytical approach, derived from how HPC resources and simulation are used in the physical sciences, to identify the epistemic boundary con- ditions of the study of clinical sepsis via the use of a proxy agent-based model of systemic inflammation. Design: Current predictive models for sepsis use correlative methods that are limited by patient heterogeneity and data sparseness. We address this issue by using an HPC version of a system-level validated agent-based model of sepsis, the Innate Immune Response ABM (IIRBM), as a proxy system in order to identify boundary conditions for the possible behavioral space for sepsis. We then apply advanced analysis derived from the study of Random Dynamical Systems (RDS) to identify novel means for characterizing system behavior and providing insight into the tractability of traditional investigatory methods. Results: The behavior space of the IIRABM was examined by simulating over 70 million sepsis patients for up to 90 days in a sweep across the following parameters: cardio-respiratory-metabolic resilience; microbial invasiveness; microbial toxigenesis; and degree of nosocomial exposure. In addition to using established methods for describing parameter space, we developed two novel methods for characterizing the behavior of a RDS: Probabilistic Basins of Attraction (PBoA) and Stochastic Trajectory Analysis (STA). Computationally generated behavioral landscapes demonstrated attractor structures around stochastic re- gions of behavior that could be described in a complementary fashion through use of PBoA and STA. The stochasticity of the boundaries of the attractors highlights the challenge for correlative attempts to characterize and classify clinical sepsis. Conclusions: HPC simulations of models like the IIRABM can be used to generate approximations of the behavior space of sepsis to both establish “boundaries of futility” with respect to existing investi- gatory approaches and apply system engineering principles to investigate the general dynamic properties of sepsis to provide a pathway for developing control strategies. Here, the issues that bedevil the study and treatment of sepsis, namely clinical data sparseness and inadequate experimental sampling of system behavior space, are fundamental to nearly all biomedical research, manifesting in the “Crisis of Reproducibility” at all levels. HPC-augmented simulation-based research offers an investigatory strategy more consistent with that seen in the physical sciences (which combine experiment, theory and simulation), and an opportunity to utilize the leading advances in HPC, namely deep machine learning and evolutionary computing, to form the basis of an iterative scientific process to meet the full promise of Precision Medicine (right drug, right patient, right time).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Wettability Alteration with Weak Acid-Assisted Surfactant Flood

Oil-wetness and heterogeneity are two key reasons for low oil recovery by waterflooding in carbonate reservoirs. Surfactants have been effective in altering the oil-wet matrix to a more water-wet condition and initiating spontaneous imbibition. Because it takes time for the surfactant to alter wettability, oil recovery from the tight matrix is slow and sometimes not economically feasible. Acids have the potential of dissolving minerals, which may alter wettability. In this study, the enhanced oil recovery (EOR) performance of an acid-assisted surfactant solution, a novel technique, was evaluated for low-temperature applications. A set of acids and their acetates were tested. Bulk rock-acid reaction, wettability alteration (WA) tests, and spontaneous imbibition measurements were conducted at reservoir temperature (35°C) to identify effective candidates. Coreflood tests were then performed to evaluate the selected acid-surfactant formulations. Before and after a coreflood test, the core was scanned using micro-computed tomography (CT) to investigate pore structure alteration. Bulk reaction measurements showed that weak acids, especially acetic acid (AA), have the desired low reaction rates at 35°C. WA tests showed that AA can remove the crude oil off the rock surface and alter wettability through mineral dissolution. The surfactant can reduce contact angles from 160° to 58°; adding acid into the surfactant can further reduce it to 52°. Spontaneous imbibition experiments showed the synergy between the acid and the surfactant; the AA-surfactant solution had the highest oil recovery (62.6%) among acid-surfactant formulations. The acid improves the WA efficiency by the surfactant through surface mineral dissolution and lower ζ-potential. The imbibition transports the acid-surfactant solutions into the matrix, which minimizes face dissolution. Coreflood tests show that the AA-surfactant flood can increase the oil recovery rate and recover about 8% more oil compared to the surfactant flood. Micro-CT showed that a few mineral particles were transported along the core and partially plugged pore throats, which reduced permeability and diverted flow leading to improved oil recovery. Here, the transport of the acid in reservoirs and the potential plugging issues have to be carefully evaluated in future studies.

58 GEOSCIENCES↗
DOE OSTI2022

Title Coherent X-ray Studies of Surface Growth and Patterning Processes

X-ray Photon Correlation Spectroscopy (XPCS) is being developed as a tool to study nanoscale dynamics of fluctuations during thin film growth and surface patterning. XPCS examines the evolution of the X-ray scattering speckle pattern in reciprocal space to reveal dynamics information not accessible through any other means. Unlike low-coherence conventional X-ray scattering, which incoherently averages over different regions of a sample, coherent X-ray scattering is sensitive to the detailed structure of a given sample at that moment in time. In thin film growth processes, heterodyning, which occurs due to coherent mixing of two scattered signals, can be used to investigate the relationship between surface growth velocity and defect propagation. And, for polycrystalline thin film growth, the spatial coherence of the X-ray beam can substitute for the missing spatial coherence of the growth process to track layered growth in detail, even in a growth regime in which there are no conventional growth oscillations of the X-ray intensity. This has opened the door for detailed dynamics studies of individual atomic layers during real-world growth and patterning, not just in perfect single-crystal growth cases, which greatly expands the applicability of in-situ X-ray scattering methods. Step-flow dynamics in mounded polycrystalline growth has been studied separately for two thin film organic semiconductors deposited by thermal deposition in a vacuum environment, C60 and diindenoperylene (DIP). Highly oriented polycrystalline thin films are readily obtained in both systems, where mounds are composed of crystalline monolayer-height steps and terraces in a so-called wedding cake morphology. The formation of mounds is understood to be due to significant Ehrlich-Schwoebel step edge barriers that inhibit molecules from hopping down from one layer to the one below. The mounds exhibit local step flow, which can be monitored using coherent X-ray scattering. This is made possible due to heterodyning between scattering from the average mounds and the moving steps, which becomes visible in XPCS analysis. The effect of desorption, i.e. re-evaporation of deposited molecules is found to be important for understanding these processes. The impact of this work is to enable testing of models of step dynamics, the shape of mounds, and merging of mounds to form continuous thin films. Highly ordered polycrystalline thin films deposited on inexpensive substrates have applications in thin film solar cells and other organic electronic devices. In a separate set of experiments, speckle analysis during self-organized ion beam nanopatterning reveals memory stretching back to the beginning of patterning in the early stages and enables measurement of the velocity of self-organized patterns across surfaces providing the possibility of stringent new tests of the theory of pattern formation and motion.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Large Eddy Simulation of turbulent flow in rod bundle with spacer grids

Pressurized water nuclear reactor fuel assemblies are made of fuel rods held by spacer grids, which are the core components of the nuclear power plants. Hot fuel rods are cooled by the coolant water flow, and the spacer grids with mixing vanes distort the mainstream, producing the swirling flow which can enhance the heat transfer to improve CHF (Critical Heat Flux). Heat transfer between fuel rods and coolant is a crucial part of thermal-hydraulic analysis. For the purpose to improve the thermal hydraulic performance of the fuel assemblies, many experiments have been conducted to understand the flow structures and three-dimensional behaviors of the turbulence in the rod bundle sub-channels and then to optimize the design of the spacer grid. In the experimental work, which requires lots of efforts to get reliable data, LDA (Laser Doppler Anemometry) and PIV (Particle Image Velocimetry) are generally used to measure the velocity profile and turbulent intensity in the sub-channels. Computational Fluid Dynamics (CFD) is another useful way to study the characteristics of the detailed flow condition and heat transfer phenomena in the subchannels. The main issue of the single-phase CFD is to accurately predict the flow and turbulence characteristics. RANS method is generally used to simulate the induced lateral velocity and the heat transfer process. With the development of computer technology, some more precise but previously time-consuming turbulent models, such as LES (Large Eddy Simulation), are gradually used for the simulation of rod bundle flow with current advanced CPU resources. Ikeno and Kajishima explained the decay of swirling turbulent flow in the rod bundle using LES method. Benhamadouche acquired the turbulence induced forces on the central rod of the rod bundle through performing the LES of the flow downstream the spacer grid. Chang and Tavoularis conducted hybrid simulations with SAS SST and LES of a 5 x 5 rod bundle, and several vortices were identified in the vicinity of the vanes and the generation mechanisms, evolution progress and influence on the velocities of these vortices were clarified. In the paper we present CFD simulations of turbulent flow in rod bundle sub-channels with spacer grids of two different experiments. One is MATiS-H benchmark experiment with scaled spacer grids, which is conducted by KAERI, and is used as a CFD benchmark test, and the other is a 5 x 5 rod bundle with typical spacer grids using PIV technology, with the cooperation of Harbin Engineering University. The main objectives of the calculations are to validate LES methodology by the benchmark experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

The Chemical Abundance Structure of the Inner Milky Way: A Signature of “Upside-down” Disk Formation

We present a model for the [α/Fe]--[Fe/H] distribution of stars in the inner Galaxy, 3 kpc<R<5 kpc, measured as a function of vertical distance |z| from the midplane by Hayden et al. (H15). Motivated by an “upside-down” scenario for thick disk formation, in which the thickness of the star-forming gas layer contracts as the stellar mass of the disk grows, we combine one-zone chemical evolution with a simple prescription in which the scale-height of the stellar distribution drops linearly from z{sub h}=0.8 kpc to z{sub h}=0.2 kpc over a timescale t {sub c}, remaining constant thereafter. We assume a linear-exponential star formation history, M-dot {sub ∗}(t)∝te{sup −t/t{sub s}{sub f}}. With a star formation efficiency timescale τ{sub ∗}=M{sub g}(t)/ M-dot {sub ∗}(t)=2 Gyr, an outflow mass-loading factor η= M-dot {sub out}(t)/ M-dot {sub ∗}(t)=1.5, t{sub sf}=3 Gyr, and t{sub c}=2.5 Gyr, the model reproduces the observed locus of inner disk stars in [α/Fe]--[Fe/H] and the metallicity distribution functions (MDFs) measured by H15 at |z|=0--0.5 kpc, 0.5--1 kpc, and 1--2 kpc. Substantial changes to model parameters lead to disagreement with the H15 data; for example, models with t{sub c}=1 Gyr or t{sub sf}=1 Gyr fail to match the observed MDF at high-|z| and low-|z|, respectively. The inferred scale-height evolution, with z {sub h}(t) dropping on a timescale t{sub c}∼t{sub sf} at large lookback times, favors upside-down formation over dynamical heating of an initially thin stellar population as the primary mechanism regulating disk thickness. The failure of our short-t {sub c} models suggests that any model in which thick disk formation is a discrete event will not reproduce the continuous dependence of the MDF on |z| found by H15. Our scenario for the evolution of the inner disk can be tested by future measurements of the |z|-distribution and the age–metallicity distribution at R=3--5 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

IDAES-PSE 2.5.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.5.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New diagnostics check for near-parallel variables and constraints. New diagnostics tools for identifying causes of infeasibility in models. New example for creating a custom model of a liquid-liquid extractor unit operation. Bug Fixes Fixed bug in Gibbs reactor that caused it to appear to have additional spurious degrees of freedom. Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Testing and Robustness Deployed the IDAES Diagnostics Toolbox to confirm that there are no structural or numerical issues in the core model libraries. Additional robustness tests for core model, and some associated improvements in the converge tester class. Fixed a number of issues that were causing unexpected warnings to be emitted during testing. Deprecations and Removals Removed examples for RIPE tool which has not been supported for a number of releases.

AS↗
DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗
DOE OSTI2019

Dark Energy Survey year 1 results: the relationship between mass and light around cosmic voids

What are the mass and galaxy profiles of cosmic voids? In this paper, we use two methods to extract voids in the Dark Energy Survey (DES) Year 1 redMaGiC galaxy sample to address this question. We use either 2D slices in projection, or the 3D distribution of galaxies based on photometric redshifts to identify voids. For the mass profile, we measure the tangential shear profiles of background galaxies to infer the excess surface mass density. The signal-to-noise ratio for our lensing measurement ranges between 10.7 and 14.0 for the two void samples. We infer their 3D density profiles by fitting models based on N-body simulations and find good agreement for void radii in the range 15–85 Mpc. Comparison with their galaxy profiles then allows us to test the relation between mass and light at the 10 per cent level, the most stringent test to date. We find very similar shapes for the two profiles, consistent with a linear relationship between mass and light both within and outside the void radius. We validate our analysis with the help of simulated mock catalogues and estimate the impact of photometric redshift uncertainties on the measurement. Our methodology can be used for cosmological applications, including tests of gravity with voids. Finally, this is especially promising when the lensing profiles are combined with spectroscopic measurements of void dynamics via redshift-space distortions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantifying Turbine-Level Risk to Golden Eagles Using a High-Fidelity Updraft Model and a Stochastic Behavioral Model

To minimize the effects of wind farms on Golden Eagle (Aquila chrysaetos) populations while enabling sustainable development of renewable energy resources, it is important to understand how eagles interact with atmospheric flows, terrain features, and anthropogenic structures. Models that predict migratory flight paths provide one tool that helps us grasp how the location of wind farms may influence interactions and impacts on migrating Golden Eagles. The current state-of-the-art in predicting migratory flight paths uses a deterministic fluid-flow analogy to predict eagle trajectory using only an orographic updraft potential computed from topographical features. This model does not take into account variables, such as thermal updrafts and time varying atmospheric conditions that are known to influence migratory behavior. In this work, we improve on the model with the objective of developing tools that advance our understanding of how atmospheric flows and terrain features affect migratory eagle behavior and their interactions with wind farms. Specifically, we 1) incorporate both orographic and thermal updraft information in simulating eagle flight paths; 2) incorporate stochasticity into eagle travel patterns to better capture the influence of exogenous factors on, and the inherent stochasticity of eagle behavior; 3) consider spatio-temporal atmospheric data at wind-farm-scale when computing updraft potential; and 4) account for how atmospheric conditions and the direction of migration change seasonally and how these changes affect eagle migratory flight behavior. We tested the model using a 50km by 50km region with 50 m resolution in the western United States. We simulated 900 independent, probabilistic eagle tracks during southerly and northerly migration, assuming eagles solely rely on orographic updrafts. The preliminary results indicate that the inclusion of finer resolution atmospheric data allows for the inclusion of realistic conditions that an eagle experiences. The stochasticity in eagle tracks provides a platform to include uncertainty in eagle decision making and help produce robust eagle presence maps. We will deploy updraft and downdraft velocities computed using a high-fidelity, wind farm scale, computational fluid dynamics solver under development at National Renewable Energy Laboratory. This work is a first step in the development of a predictive and generalizable eagle behavior model at the wind farm scale that does not rely on empirical data collection. Although the current model is intended for migratory eagles, we will extend and refine this model to inform the development of additional behavioral modes, including resident eagle behavior. This modeling approach improves our ability to understand eagle use of the landscape at a fine scale, and it is our hope that this work will ultimately help advance strategies that minimize the impact of wind development on Golden Eagle populations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2024

Pinus ponderosa physiology data from the 2023 girdling manipulation field experiment at Game Ridge, Missoula County, Montana, USA

This dataset contains physiology data on Pinus ponderosa saplings, including: leaf water potential (MPa, megapascals), leaf osmotic potential (MPa, megapascals), leaf pressure potential (MPa, megapascals), leaf relative water content (% saturated mass), as well as soluble sugars (% dry mass),starch (% dry mass), and total nonstructural carbohydrates (% dry mass) found in branch phloem and xylem. Data files and data dictionary(ies) are uploaded as .csv files and .xlsx files. The Users Guide is a .pdf file. Location data can be found in the Google Earth file GameRidge_SitePlotCoordinates.kmz.kml included here. These datasets were collected for Plant Carbohydrate Depletion, Mycorrhizal Networks, and Vulnerability to Drought: An Experimental Test in the Field. This experiment examined the interdependency between plant hydraulics and carbohydrate availability and sought to develop ways to incorporate interactions with below ground symbiotic organisms to better model and quantify forest response to drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Uncertainty Quantification in Vibro-Acoustics Diagnosis of Alkali-Silica Reaction Degradation in Medium-Sized Concrete Samples

The objectives of this ongoing research project focus on health monitoring and data analytics of concrete slabs containing reactive aggregates and thus subjected to degradation due to alkali-silica reaction (ASR). A controlled concrete slab with four pockets of reactive aggregates (pure silica, wells, placitas, and spratt) was cast at Vanderbilt University and cured in representative conditions to accelerate degradation due to ASR. A set of four concrete samples were also cast and cured at the University of Alabama for ASR testing. Of these four samples, two slabs contained reactive aggregates while the other two had the non-reactive aggregate counterparts mixed throughout the samples. Vibro-acoustic testing was used on these slabs to locate ASR damage within the reactive samples. Vibro-acoustic modulation (VAM) is a vibration-based NDE method that utilizes signatures of nonlinear dynamic interactions on contact surfaces of crack or delamination damage to detect and localize the damage. VAM analysis was conducted on both the Vanderbilt and Alabama samples using multiple variables for damage detection and localization. The results from the data analysis of the vibro-acoustic testing on concrete slabs cured at Vanderbilt University and the University of Alabama are discussed in detail in this report. Results for damage localization are dependent on multiple variables used in the vibro-acoustic modulation experiments. A major focus of this report is to quantify the uncertainty in the diagnosis due to multiple factors and uncertainty sources. The uncertainty quantification methodology presented in this report was applied to VAM-based diagnosis and prognosis. However, the methodology is general, and is capable of being applied to multiple technique that collect spatially distributed data. Future work needs to investigate the incorporation of uncertainty quantification in developing a robust Prognostics and Health Management framework. Digital image correlation is a three-dimensional, full-field, optical NDE technique to measure contour, deformation, vibration, and strain. The application of the digital image correlation technique to study ASR-related degradation on a large concrete specimen at the University of Tennessee is discussed and observations are presented in this report.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2018

The Role of Dimensionality in Understanding Granuloma Formation

Within the first 2–3 months of a Mycobacterium tuberculosis (Mtb) infection, 2–4 mm spherical structures called granulomas develop in the lungs of the infected hosts. These are the hallmark of tuberculosis (TB) infection in humans and non-human primates. A cascade of immunological events occurs in the first 3 months of granuloma formation that likely shapes the outcome of the infection. Understanding the main mechanisms driving granuloma development and function is key to generating treatments and vaccines. In vitro, in vivo, and in silico studies have been performed in the past decades to address the complexity of granuloma dynamics. This study builds on our previous 2D spatio-temporal hybrid computational model of granuloma formation in TB (GranSim) and presents for the first time a more realistic 3D implementation. We use uncertainty and sensitivity analysis techniques to calibrate the new 3D resolution to non-human primate (NHP) experimental data on bacterial levels per granuloma during the first 100 days post infection. Due to the large computational cost associated with running a 3D agent-based model, our major goal is to assess to what extent 2D and 3D simulations differ in predictions for TB granulomas and what can be learned in the context of 3D that is missed in 2D. Our findings suggest that in terms of major mechanisms driving bacterial burden, 2D and 3D models return very similar results. For example, Mtb growth rates and molecular regulation mechanisms are very important both in 2D and 3D, as are cellular movement and modulation of cell recruitment. The main difference we found was that the 3D model is less affected by crowding when cellular recruitment and movement of cells are increased. Overall, we conclude that the use of a 2D resolution in GranSim is warranted when large scale pilot runs are to be performed and if the goal is to determine major mechanisms driving infection outcome (e.g., bacterial load). To comprehensively compare the roles of model dimensionality, further tests and experimental data will be needed to expand our conclusions to molecular scale dynamics and multi-scale resolutions.

2D vs. 3D↗
DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

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

A systematic study of lanthanide titanates (A 2 Ti 2 O 7 ) chemical durability: corrosion mechanisms and control parameters

Dense lanthanide titanates ranging from La 2 Ti 2 O 7 to Yb 2 Ti 2 O 7 were synthesized by spark plasma sintering with similar micro-sized microstructure, and their chemical durability and degradation mechanisms were investigated by semi-dynamic leaching tests. Different corrosion mechanisms from incongruent to congruent dissolution were observed with preferential elemental releases of lanthanides in A 2 Ti 2 O 7 (A = La to Gd) and Ti in Er 2 Ti 2 O 7 and Yb 2 Ti 2 O 7 . Initial leaching was controlled by surface dissolution and long-term dissolution shifted to a diffusion-controlled mechanism for titanate pyrochlores from Sm 2 Ti 2 O 7 to Yb 2 Ti 2 O 7 . Both short-term and long-term leaching rates decreased from La 2 Ti 2 O 7 to Yb 2 Ti 2 O 7 , and materials became more corrosion resistant with decreasing ionic radius ratios of the rare-earth elements and titanium. Amorphous passivation films of several nm-thick enriched with Ti formed on the surfaces of the corroded materials, and the formation mechanisms were governed either by surface reorganization with preferential release of lanthanides or dissolution-precipitation from saturated bulk solutions with preferential release of titanium. Finally, these results demonstrate that the ionic size plays a dominant role in determining chemical durability of lanthanide titanates and a strong correlation is established among chemical composition/structural variation, passivation film formation and corrosion mechanisms, and materials chemical durability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗