Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snapshot”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗
DOE OSTI2025

Filling data analysis gaps in time-resolved crystallography by machine learning

There is a growing understanding of the structural dynamics of biological molecules fueled by x-ray crystallography experiments. Time-resolved serial femtosecond crystallography (TR-SFX) with x-ray Free Electron Lasers allows the measurement of ultrafast structural changes in proteins. Nevertheless, this technique comes with some limitations. One major challenge is the quality of data from TR-SFX measurements, which often faces issues like data sparsity, partial recording of Bragg reflections, timing errors, and pixel noise. To overcome these difficulties, conventionally, large volumes of data are collected and grouped into a few temporal bins. The data in each bin are then averaged and paired with the mean of their corresponding jittered timestamps. This procedure provides one structure per bin, resulting in a limited number of averaged structures for the entire time interval spanned by the experiment. Therefore, the information on ultrafast structural dynamics at high temporal resolution is lost. This has initiated research for advanced methods of analyzing experimental TR-SFX data beyond the standard binning and averaging method. To address this problem, we use a machine learning algorithm called Nonlinear Laplacian Spectral Analysis (NLSA), which has emerged as a promising technique for studying the dynamics of complex systems. In this work, we demonstrate the power of this algorithm using synthetic x-ray diffraction snapshots from a protein with significant data incompleteness, timing uncertainties, and noise. Our study confirms that NLSA is a suitable approach that effectively mitigates the effects of these artifacts in TR-SFX data and recovers accurate structural dynamics information hidden in such data.

Trujillo, Justin (ORCID:0000000285505360)↗
DOE OSTI2020

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗
DOE OSTI2020

Reconstruction of the interatomic forces from dynamic scanning transmission electron microscopy data

We explore the possibility for reconstruction of the generative physical models describing interactions between atomic units in solids from observational electron microscopy data. In this work, scanning transmission electron microscopy (STEM) is used to observe the dynamic motion of Si atoms at the edge of monolayer graphene under continuous electron beam illumination. The resulting time-lapsed STEM images represent the snapshots of observed chemical states of the system. We use two approaches: potential of mean force calculation using a radial distribution function and a direct fitting of the graphene–Si interatomic pairwise potentials with force matching, to reconstruct the force fields in the materials. These studies lay the foundation for quantitative analysis of materials energetics from STEM data through the sampling of the metastable states in the chemical space of the system.

42 ENGINEERING↗
DOE OSTI2021

H- ion source research and development at the Oak Ridge National Laboratory

The U.S. Spallation Neutron Source (SNS) is a state-of-the-art neutron scattering facility delivering the world's most intense pulsed neutron beams to a wide array of instruments which are used to conduct investigations in many fields of science and engineering. Neutrons are produced by spallation of liquid Hg due to bombardment of short (~1µs), intense (~35 A) pulses of protons delivered at 60 Hz from a storage ring which is fed by a high-intensity, 1 GeV H- LINAC. This facility has operated almost continuously since 2006, with ion source performance increasing over those years, and currently providing 50-60 mA of H- ions with a duty-factor of 6% for maintenance-free runs of ~120 days with near 100% availability. Ion source research and development at ORNL has played a key role in enabling and supporting these achievements: this report provides some historical highlights of this effort, describes our current ion source testing facilities, discusses the goals of the program and finally provides a snapshot of some of our ongoing R&D activities. In particular, we have recently simplified and improved the reliability of the plasma ignition gun for the external antenna ion source which is discussed in detail.

Welton, Robert F.↗
DOE OSTI2022

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An LES study of secondary motion and wall shear stresses in a pipe bend

In this study, large-eddy simulations (LES) of a single-phase, turbulent flow in a 90 ° pipe bend are performed at three Reynolds numbers (5300, 27 000, and 45 000) to investigate the correlation between secondary flow motion and wall shear stresses, which is suspected to be a potential mechanism responsible for material erosion. The isothermal flows are validated against available experimental and numerical data first. The snapshot proper orthogonal decomposition (POD) is applied for the medium and high Reynolds number flows to identify the secondary flow motions and the oscillation of the Dean vortices that are found to cause swirl-switching. Distinguished frequencies of the POD time coefficients at Strouhal numbers of 0.25 and 0.28 are identified for Reynolds numbers at 27 000 and 45 000, respectively. Moreover, shear stress on the pipe wall and the associated power spectral density are obtained and shown to have the same oscillating frequency as the swirl-switching.

42 ENGINEERING↗
DOE OSTI2022

Studying the Richtmyer–Meshkov instability in convergent geometry under high energy density conditions using the Decel platform

The “Decel” platform at Sandia National Laboratories investigates the Richtmyer–Meshkov instability (RMI) in converging geometry under high energy density conditions [Knapp et al., Phys. Plasmas 27, 092707 (2020)]. In Decel, the Z machine magnetically implodes a cylindrical beryllium liner filled with liquid deuterium, launching a converging shock toward an on-axis beryllium rod machined with sinusoidal perturbations. The passage of the shock deposits vorticity along the Be/D2 interface, causing the perturbations to grow. Here, we present platform improvements along with recent experimental results. To improve the stability of the imploding liner to the magneto Rayleigh–Taylor instability, we modified its acceleration history by shortening the Z electrical current pulse. Next, we introduce a “split rod” configuration that allows two axial modes to be fielded simultaneously in different axial locations along the rod, doubling our data per experiment. We then demonstrate that asymmetric slots in the return current structure modify the magnetic drive pressure on the surface of the liner, advancing the evolution on one side of the rod by multiple ns compared to its 180° counterpart. This effectively enables two snapshots of the instability at different stages of evolution per radiograph with small deviations of the cross-sectional profile of the rod from the circular. Using this platform, we acquired RMI data at 272 and 157 μm wavelengths during the single shock stage. Finally, we demonstrate the utility of these data for benchmarking simulations by comparing calculations using ALEGRA MHD and RageRunner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of multi neutron-to-gamma converter array for measuring time resolved ion temperature of inertial confinement fusion implosions

The ion temperature varying during inertial confinement fusion implosions changes the amount of Doppler broadening of the fusion products, creating subtle changes in the fusion neutron pulse as it moves away from the implosion. A diagnostic design to try to measure these subtle effects is introduced—leveraging the fast time resolution of gas Cherenkov detectors along with a multi-puck array that converts a small amount of the neutron pulse into gamma-rays, one can measure multiple snapshots of the neutron pulse at intermediate distances. Further, precise measurements of the propagating neutron pulse, specifically the variation in the peak location and the skew, could be used to infer time-evolved ion temperature evolved during peak compression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Data Science Shows that Entropy Correlates with Accelerated Zeolite Crystallization in Monte Carlo Simulations

We have performed a data science study of Monte Carlo simulation trajectories to understand factors that can accelerate formation of zeolite nanoporous crystals, a process that can take days or even weeks. In previous work, Monte Carlo simulations predicted and experiments confirmed that using a secondary organic structure-directing agent (OSDA) accelerates crystallization of all-silica LTA zeolite, with experiments finding a three-fold speedup [PCCP 24, 142-148 (2022)]. However, it remains unclear what physical factors cause the speed-up. Here, we apply data science to analyze the simulation trajectories to discover what drives accelerated zeolite crystallization in Monte Carlo going from a one-OSDA synthesis (1OSDA) to a two-OSDA version (2OSDA). We encoded simulation snapshots using the Smooth Overlap of Atomic Positions approach, which represents all 2- and 3-body correlations within a given cutoff distance. Principal component analyses failed to discriminate datasets of structures from 1OSDA and 2OSDA simulations, while the Support Vector Machine (SVM) approach succeeded at classifying such structures with an area-under-curve (AUC) score of 0.99 (where AUC = 1 is a perfect classification) with all 3-body correlations, and as high as 0.94 with only 2-body correlations. SVM decision functions reveal relatively broad / narrow histograms for 1OSDA / 2OSDA datasets, suggesting that the two simulations differ strongly in information heterogeneity. Informed by these results, we performed pair (2-body) entropy calculations during crystallization, resulting in entropy differences that semi-quantitatively account for the speedup observed in the previous Monte Carlo simulations. We conclude that altering synthesis conditions in ways that substantially changes the entropy of labile silica networks may accelerate zeolite crystallization, and we discuss possible approaches for achieving such acceleration.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Observing 3-hydroxyanthranilate-3,4-dioxygenase in action through a crystalline lens

The synthesis of quinolinic acid from tryptophan is a critical step in the de novo biosynthesis of nicotinamide adenine dinucleotide (NAD + ) in mammals. Herein, the nonheme iron-based 3-hydroxyanthranilate-3,4-dioxygenase responsible for quinolinic acid production was studied by performing time-resolved in crystallo reactions monitored by UV-vis microspectroscopy, electron paramagnetic resonance (EPR) spectroscopy, and X-ray crystallography. Seven catalytic intermediates were kinetically and structurally resolved in the crystalline state, and each accompanies protein conformational changes at the active site. Among them, a monooxygenated, seven-membered lactone intermediate as a monodentate ligand of the iron center at 1.59-Å resolution was captured, which presumably corresponds to a substrate-based radical species observed by EPR using a slurry of small-sized single crystals. Other structural snapshots determined at around 2.0-Å resolution include monodentate and subsequently bidentate coordinated substrate, superoxo, alkylperoxo, and two metal-bound enol tautomers of the unstable dioxygenase product. These results reveal a detailed stepwise O-atom transfer dioxygenase mechanism along with potential isomerization activity that fine-tunes product profiling and affects the production of quinolinic acid at a junction of the metabolic pathway.

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

Mobile loop dynamics in adenosyltransferase control binding and reactivity of coenzyme B 12

Cobalamin is a complex organometallic cofactor that is processed and targeted via a network of chaperones to its dependent enzymes. AdoCbl (5'-deoxyadenosylcobalamin) is synthesized from cob(II)alamin in a reductive adenosylation reaction catalyzed by adenosyltransferase (ATR), which also serves as an escort, delivering AdoCbl to methylmalonyl-CoA mutase (MCM). The mechanism by which ATR signals that its cofactor cargo is ready (AdoCbl) or not [cob(II)alamin] for transfer to MCM, is not known. In this study, we have obtained crystallographic snapshots that reveal ligand-induced ordering of the N terminus of Mycobacterium tuberculosis ATR, which organizes a dynamic cobalamin binding site and exerts exquisite control over coordination geometry, reactivity, and solvent accessibility. Cob(II)alamin binds with its dimethylbenzimidazole tail splayed into a side pocket and its corrin ring buried. The cosubstrate, ATP, enforces a four-coordinate cob(II)alamin geometry, facilitating the unfavorable reduction to cob(I)alamin. The binding mode for AdoCbl is notably different from that of cob(II)alamin, with the dimethylbenzimidazole tail tucked under the corrin ring, displacing the N terminus of ATR, which is disordered. In this solvent-exposed conformation, AdoCbl undergoes facile transfer to MCM. The importance of the tail in cofactor handover from ATR to MCM is revealed by the failure of 5'-deoxyadenosylcobinamide, lacking the tail, to transfer. In the absence of MCM, ATR induces a sacrificial cobalt–carbon bond homolysis reaction in an unusual reversal of the heterolytic chemistry that was deployed to make the same bond. Finally, the data support an important role for the dimethylbenzimidazole tail in moving the cobalamin cofactor between active sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Lithium superoxide encapsulated in a benzoquinone anion matrix

Lithium peroxide is the crucial storage material in lithium–air batteries. Understanding the redox properties of this salt is paramount toward improving the performance of this class of batteries. Lithium peroxide, upon exposure to p –benzoquinone (p–C 6 H 4 O 2 ) vapor, develops a deep blue color. This blue powder can be formally described as [Li 2 O 2 ] 0.3 • [LiO 2 ] 0.7 • {Li[p–C 6 H 4 O 2 ]} 0.7 , though spectroscopic characterization indicates a more nuanced structural speciation. Infrared, Raman, electron paramagnetic resonance, diffuse-reflectance ultraviolet-visible and X-ray absorption spectroscopy reveal that the lithium salt of the benzoquinone radical anion forms on the surface of the lithium peroxide, indicating the occurrence of electron and lithium ion transfer in the solid state. As a result, obligate lithium superoxide is formed and encapsulated in a shell of Li[p–C 6 H 4 O 2 ] with a core of Li 2 O 2 . Lithium superoxide has been proposed as a critical intermediate in the charge/discharge cycle of Li–air batteries, but has yet to be isolated, owing to instability. The results reported herein provide a snapshot of lithium peroxide/superoxide chemistry in the solid state with redox mediation.

25 ENERGY STORAGE↗
DOE OSTI2021

A conserved epitope III on hepatitis C virus E2 protein has alternate conformations facilitating cell binding or virus neutralization

Epitope III, a highly conserved amino acid motif of 524 APTYSW 529 on the hepatitis C virus (HCV) E2 glycoprotein, resides in the critical loop that binds to the host receptor CD81, thus making it one of the most important antibody targets for blocking HCV infections. Here, we have determined the X-ray crystal structure of epitope III at a 2.0-Å resolution when it was captured by a site-specific neutralizing antibody, monoclonal antibody 1H8 (mAb1H8). The snapshot of this complex revealed that epitope III has a relatively rigid structure when confined in the binding grooves of mAb1H8, which confers the residue specificity at both ends of the epitope. Such a high shape complementarity is reminiscent of the “lock and key” mode of action, which is reinforced by the incompatibility of an antibody binding with an epitope bearing specific mutations. By subtly positioning the side chains on the three residues of Tyr 527 , Ser 528 , and Trp 529 while preserving the spatial rigidity of the rest, epitope III in this cocrystal complex adopts a unique conformation that is different from previously described E2 structures. With further analyses of molecular docking and phage display–based peptide interactions, we recognized that it is the arrangements of two separate sets of residues within epitope III that create these discrete conformations for the epitope to interact selectively with either mAb1H8 or CD81. These observations thus raise the possibility that local epitope III conformational dynamics, in conjunction with sequence variations, may act as a regulatory mechanism to coordinate “mAb1H8-like” antibody-mediated immune defenses with CD81-initiated HCV infections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Rippled metamaterials with scale-dependent tailorable elasticity

Thermally induced ripples are intrinsic features of nanometer-thick films, atomically thin materials, and cell membranes, significantly affecting their elastic properties. Despite decades of theoretical studies on the mechanics of suspended thermalized sheets, controversy still exists over the impact of these ripples, with conflicting predictions about whether elasticity is scale-dependent or scale-independent. Experimental progress has been hindered so far by the inability to have a platform capable of fully isolating and characterizing the effects of ripples. This knowledge gap limits the fundamental understanding of thin materials and their practical applications. Here, we show that thermal-like static ripples shape thin films into a class of metamaterials with scale-dependent, customizable elasticity. Utilizing a scalable semiconductor manufacturing process, we engineered nanometer-thick films with precisely controlled frozen random ripples, resembling snapshots of thermally fluctuating membranes. Resonant frequency measurements of rippled cantilevers reveal that random ripples effectively renormalize and enhance the average bending rigidity and sample-to-sample variations in a scale-dependent manner, consistent with recent theoretical estimations. The predictive power of the theoretical model, combined with the scalability of the fabrication process, was further exploited to create kirigami architectures with tailored bending rigidity and mechanical metamaterials with delayed buckling instability.

Applied Physical Sciences↗
DOE OSTI2023

On the Language of Reliability: A System Engineer Perspective

In its classical definition, risk is defined by three elements: what can go wrong, what are its consequences, and how likely is it to occur? While this definition makes sense in a regulatory-based framework where for the current fleet of operating light water reactors (LWRs), the risks associated with nuclear power plants typically are characterized in terms of core damage and large early release frequency (LERF), this approach does not provide a useful snapshot of the health of the plant from a broader perspective. This is due to the very narrow context in which the term “risk” typically is defined as nuclear safety aspects that have the potential to impact public health. In this paper, we take the viewpoint of nuclear safety that is reflective of the current fleet of operating LWRs for which core damage frequency and LERF are appropriate metrics. For other advanced reactor designs, other more applicable technology neutral metrics of reactor safety metrics would be specified. A possible alternate path would start by redefining the word risk with a broader meaning that better reflects the needs of a system health and asset management decision-making process. Rather than asking how likely an event could occur (in probabilistic terms), we can ask how far this event is from occurring. Our approach starts by defining and quantifying component and system health in terms of a “distance” between its actual and limiting conditions, i.e., determination of the margin that exists between the current state/condition and the state where the component/system is no longer capable of achieving its intended function. A margin is a measure that is more reflective of the current state or performance of a component, and therefore more closely tied to decisions that are made on an ongoing basis. Further, we will show how, given the data available from plant equipment reliability and monitoring (e.g., pump vibration data) and prognostic (e.g., component remaining useful life estimation) data, a margin can be described and determined for all types of maintenance approaches (e.g., corrective or predictive maintenance). We show how classical reliability models (e.g., fault trees) can be used to quantify the system margin provided component margin values. In the approach described in this paper, the propagation of margin values through classical reliability models are not performed using classical probabilistic calculations applied to sets (as performed in a typical plant probabilistic risk assessment). Instead, we show how it is possible to propagate margin values through Boolean logic gates (i.e., AND and OR operators) through distance-based operations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Offline Maximizing Minimally Invasive Proper Orthogonal Decomposition for Reduced-Order Modeling of S n Radiation Transport

Deterministic solutions to the Sn radiation transport equation can be computationally expensive to calculate. Reduced-order modeling enables efficient approximation of the full-order model (FOM) solution. We propose a novel method for constructing reduced-order models (ROMs) of the S n radiation transport equation, offline maximizing minimally invasive (OMMI) proper orthogonal decomposition (POD). POD uses the method of snapshots to create a reduced-order basis for constructing an ROM. Minimally invasive POD leverages the sweep infrastructure existing in deterministic transport codes to create a POD-based ROM, even when infeasible by traditional methods. Offline maximizing minimally invasive proper orthogonal decomposition (OMMI-POD) extends minimally invasive POD by performing sweeps offline, therefore maximizing the potential speedup. OMMI-POD does so by creating a library of reduced systems from a training set. This library of reduced systems is then interpolated to provide a rapid approximate solution of the S n radiation transport equation. The model is evaluated on a set of test problems, achieving a low error with a 466 times speedup over the FOM. Also presented is a study of the effect of sampling method on the performance of OMMI-POD, specifically comparing naive uniform sampling to the more accurate and computationally expensive greedy sampling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Metalloproteins in an Era of Modern Crystallography and Why the Beamline Matters

The Structural Molecular Biology (SMB) macromolecular crystallography (MC) group at the Stanford Synchrotron Radiation Lightsource (SSRL) have developed state-of-the-art capabilities tailored for metalloenzyme structural analysis. Metalloproteins sit at the center of biology’s most audacious chemistry. From multi-electron redox catalysis to radical rearrangements and light-driven transformations, metal sites give proteins access to reaction landscapes that would otherwise be inaccessible under ambient conditions. Yet their study presents a fundamental paradox for MC studies: the very X-rays we use to reveal atomic structure can alter the electronic states we seek to understand. As the field moves beyond static snapshots toward mechanistic insight, success increasingly depends on our ability to maintain metal centers in their native state throughout the experiment. The SSRL SMB-MC beamlines integrate a suite of capabilities specifically designed to address these challenges. By combining in situ spectroscopic verification, intelligent dose management, controlled reaction initiation, optimized anomalous diffraction, and real-time crystallographic diffraction analysis, these tools enable researchers to interrogate metalloproteins with unprecedented rigor. This article explores how these complementary approaches are reshaping our ability to capture metalloprotein chemistry, and what this means for mechanistic studies at synchrotron beamlines.

Maggiolo, Ailiena O. [SLAC National Accelerator La↗