Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials image simulations”

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 577 records · Page 32

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Image Distinguishability Analysis Testing Through Principal Components and Its Application to Hot Spot Scale Invariance

Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.

organic↗

Driving magnetic domains at the nanoscale by interfacial strain-induced proximity

Here, we investigate the local nanoscale changes of the magnetic anisotropy of a Ni film subject to an inverse magnetostrictive effect by proximity to a V 2 O 3 layer. Using temperature-dependent photoemission electron microscopy (PEEM) combined with X-ray magnetic circular dichroism (XMCD), direct images of the Ni spin alignment across the first-order structural phase transition (SPT) of V 2 O 3 were obtained. We find an abrupt temperature-driven reorientation of the Ni magnetic domains across the SPT, which is associated with a large increase of the coercive field. Moreover, angular dependent ferromagnetic resonance (FMR) shows a remarkable change in the magnetic anisotropy of the Ni film across the SPT of V 2 O 3 . Micromagnetic simulations based on these results are in quantitative agreement with the PEEM data. Direct measurements of the lateral correlation length of the Ni domains from XMCD images show an increase of almost one order of magnitude at the SPT compared to room temperature, as well as a broad spatial distribution of the local transition temperatures, thus corroborating the phase coexistence of Ni anisotropies caused by the V 2 O 3 SPT. We show that the rearrangement of the Ni domains is due to strain induced by the oxide layers’ structural domains across the SPT. Our results illustrate the use of alternative hybrid systems to manipulate magnetic domains at the nanoscale, which allows for engineering of coercive fields for novel data storage architectures.

36 MATERIALS SCIENCE↗

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗

Super-resolution Imaging of Plasmonic Near-Fields: Overcoming Emitter Mislocalizations

Plasmonic nano-objects have shown great potential in enhancing applications like biological/chemical sensing, light harvesting and energy transfer, and optical/quantum computing. Therefore, an extensive effort has been vested in optimizing plasmonic systems and exploiting their field enhancement properties. Super-resolution imaging with quantum dots (QDs) is a promising method to probe plasmonic near-fields but is hindered by the distortion of the QD radiation pattern. Here, we investigate the interaction between QDs and “L-shaped” gold nanoantennas and demonstrate both theoretically and experimentally that this strong interaction can induce polarization-dependent modifications to the apparent QD emission intensity, polarization, and localization. Based on FDTD simulations and polarization-modulated single-molecule microscopy, we show that the displacement of the emitter’s localization is due to the position-dependent interference between the emitter and the induced dipole, and can be up to 100 nm. Our results help pave a pathway for higher precision plasmonic near-field mapping and its underlying applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weak shock compaction on granular salt

This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-Speed Diagnostic and Simulation Capabilities for Reacting Hypersonic Reentry Flows (LDRD Final Report)

High-enthalpy hypersonic flight represents an application space of significant concern within the current national-security landscape. The hypersonic environment is characterized by high-speed compressible fluid mechanics and complex reacting flow physics, which may present both thermal and chemical nonequilibrium effects. We report on the results of a three-year LDRD effort, funded by the Engineering Sciences Research Foundation (ESRF) investment area, which has been focused on the development and deployment of new high-speed thermochemical diagnostics capabilities for measurements in the high-enthalpy hypersonic environment posed by Sandia's free-piston shock tunnel. The project has additionally sponsored model development efforts, which have added thermal nonequilibrium modeling capabilities to Sandia codes for subsequent design of many of our shock-tunnel experiments. We have cultivated high-speed, chemically specific, laser-diagnostic approaches that are uniquely co-located with Sandia's high-enthalpy hypersonic test facilities. These tools include picosecond and nanosecond coherent anti-Stokes Raman scattering at 100-kHz rates for time-resolved thermometry, including thermal nonequilibrium conditions, and 100-kHz planar laser-induced fluorescence of nitric oxide for chemically specific imaging and velocimetry. Key results from this LDRD project have been documented in a number of journal submissions and conference proceedings, which are cited here. The body of this report is, therefore, concise and summarizes the key results of the project. The reader is directed toward these reference materials and appendices for more detailed discussions of the project results and findings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Predictive models enhance feedstock quality of corn stover via air classification

Feedstock heterogeneity is a fundamental obstacle to cost-competitive biobased products. Agricultural products like corn stover have anatomical components that vary in their chemical composition, mechanical properties, structure, and response to chemical and biological treatments. A technique that can enrich streams in select anatomical fractions would allow a tailored deconstruction approach to increase overall process efficiency. Air classification can be leveraged for such refining, however, fundamental characterization and understanding of the particle properties that underly the physics of air classification are only modestly documented. Here, we determine fundamental particle properties including mass-to-area ratio, drag coefficient, and partition velocity that describe how anatomical tissues of corn stover behave during air classification. In this work, mass-to-area ratios of anatomical tissues vary by nearly two orders of magnitude from 2.3 mg/mm 2 for cob to 0.04 mg/mm 2 for leaf. Drag coefficients of longer, fibrous materials (i.e., rind, husk, and sheath) are shown to correlate with particle area (p-value < 0.001) whereas granular tissues (i.e., cob, pith, and leaf) correlate better with mass-to-area ratio (p-values < 0.001). When compared to experimental observations, a simulated two-stage air classification and size reduction scenario predicts the overall partitioning of anatomical tissues within 15% for pith, husk, rind, and cob tissues. The model predicts an air-classified fraction preferentially enriched in cob (purity = 20%), rind (purity = 74%), and pith (purity = 4.5%) with a mass yield of 47%. Empirical relations for these properties can be used to predict the partitioning of corn stover during air classification based on anatomical type and size.

09 BIOMASS FUELS↗

Microwave Imaging of Edge Conductivity in Graphene at Charge Neutrality and Quantum Hall States

Here, we report local conductivity imaging of edge states in monolayer graphene by millikelvin microwave impedance microscopy. At the charge-neutrality point, as the magnetic field increases, the local conductivity at the edge drops to zero more slowly than in the bulk. This behavior is consistent with the calculated spatial profile of the charge gap in the canted antiferromagnetic phase. For comparison, we also perform microwave imaging of integer quantum Hall states away from neutrality, which host dissipationless chiral edge channels. The evolution of the edge signal as a function of the bulk gap is fundamentally different between the Landau level filling factor 𝜈 = 0 and |𝜈| ≥ 1 integer quantum Hall states, which can be qualitatively explained by numerical simulations and theoretical analysis. Our results provide a comprehensive microscopic picture of the edge and bulk states as the Fermi level moves across the unique Landau-level spectrum of graphene.

Hall effect↗

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Topology-driven compressive behavior of Inconel 718 lattice structures with Z-strut reinforcement fabricated by laser powder bed fusion

This study investigates the compressive deformation behavior and mechanical performance of Inconel 718 lattice structures fabricated by laser powder bed fusion (LPBF). Four unit-cell topologies—BCC, FCC, BCCZ, and FCCZ—were designed with a unit-cell size of 3 mm and fabricated under identical process conditions to isolate the effect of topology. Measured relative densities ranged from 14.65% to 17.72%. Compressive testing showed that Z-strut-reinforced topologies (BCCZ: 54.6 MPa; FCCZ: 80.2 MPa) exhibited higher strength than their unreinforced counterparts, which may be associated with mixed-mode deformation behavior enabled by the vertically aligned Z-struts. Finite element simulations and Digital Image Correlation (DIC) analysis support the observation of a transition from node-dominated deformation in BCC/FCC to mixed-mode deformation in BCCZ/FCCZ. These findings suggest that unit-cell topology is a key design variable for tailoring deformation mechanisms in LPBF lattice structures.

36 MATERIALS SCIENCE↗

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-Efficient Generation of Protein Conformational Ensembles with Backbone-to-Side-Chain Transformers

Excitement at the prospect of using data-driven generative models to sample configurational ensembles of biomolecular systems stems from the extraordinary success of these models on a diverse set of high-dimensional sampling tasks. Unlike image generation or even the closely related problem of protein structure prediction, there are currently no data sources with sufficient breadth to parametrize generative models for conformational ensembles. To enable discovery, a fundamentally different approach to building generative models is required: models should be able to propose rare, albeit physical, conformations that may not arise in even the largest data sets. Here, in this work, we introduce a modular strategy to generate conformations based on “backmapping” from a fixed protein backbone that (1) maintains conformational diversity of the side chains and (2) couples the side-chain fluctuations using global information about the protein conformation. Our model combines simple statistical models of side-chain conformations based on rotamer libraries with the now ubiquitous transformer architecture to sample with atomistic accuracy. Together, these ingredients provide a strategy for rapid data acquisition and hence a crucial ingredient for scalable physical simulation with generative neural networks.

36 MATERIALS SCIENCE↗

Data-driven prediction of geometry- and toolpath sequence-dependent intra-layer process conditions variations in laser powder bed fusion

Geometrical features and toolpath sequence are two important factors that cause process condition variations, such as variations in the meltpool temperature or meltpool size, that might lead to undesired material properties in the laser powder bed fusion (LPBF) process. Due to the high dynamics and complex physics of the LPBF process, it is difficult to predict variations in process conditions with simulations alone. Advances in measurement technology and computational technologies open up new possibilities for smart manufacturing. In this paper, a data-driven method to predict intra-layer variations in the processing conditions that source from the toolpath sequence and part geometry is presented. The approach is demonstrated using two-color on-axis pyrometer measurements. Three demonstration cases are presented in which it is demonstrated (1) how the trained predictive model can be used as a filter to ease the interpretation of process variations and discover patterns related to toolpath and part geometry, and (2) how to generate predictions that can be used for feedforward control, i.e., for adjusting laser power or scanning speed along the toolpath using a meltpool temperature prediction model generated based on on-axis measurements. Results show that the developed prediction model is able to meaningfully predict process variations resulted from toolpath sequence and geometry. Predictions are aligned with the results from the related work of others and for the case of 180° laser path turnarounds in our high-speed X-ray imaging experiments. In conclusion, the potential issues related to the current maturity status of the process and measuring equipment that could in practice affect the performance of the proposed solutions are also discussed.

42 ENGINEERING↗

Measurement of mix at the fuel–ablator interface in indirectly driven capsule implosions on the National Ignition Facility

The interface between the capsule ablator and fuel ice layer is susceptible to hydrodynamic instabilities. The subsequent mixing of hot ablator material into the ice reduces fuel compression at stagnation and is a candidate for reduced capsule performance. The ability to diagnose ice–ablator mix is critical to understanding and improving stability at this interface. Combining the crystal backlighter imager with the single line of sight camera on the National Ignition Facility (NIF) allows direct measurement of ice–ablator mix by providing multiple quasi-monochromatic radiographs of layered capsule implosions per experiment with high spatial (∼12 μm) and temporal (∼35 ps) resolution. The narrow bandwidth of this diagnostic platform allows radiography of the inner edge of the capsule limb close to stagnation without capsule self-emission contaminating the data and removes opacity uncertainties typically associated with the spectral content of the radiograph. Analysis of radiographic data via a parameterized forward-fitting Abel inversion technique provides measurements of the distribution of mix mass inwards from the ice–ablator interface. The sensitivity of this mix measurement technique was demonstrated by applying it to layered experiments in which the stability of the ice–ablator interface was expected to vary significantly. Additional experiments suggest that high-density carbon capsules that employ a buried-layer dopant profile suffer from mixing at the innermost doped–undoped interface. Data from these experiments suggest that opacity models used in hydrodynamic simulations of NIF experiments can potentially over-predict the opacity of doped capsules. LLNL-JRNL-850535-DRAFT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

In‐Operando Spatiotemporal Imaging of Coupled Film‐Substrate Elastodynamics During an Insulator‐to‐Metal Transition

Abstract The drive toward non‐von Neumann device architectures has led to an intense focus on insulator‐to‐metal (IMT) and the converse metal‐to‐insulator (MIT) transitions. Studies of electric field‐driven IMT in the prototypical VO 2 thin‐film channel devices are largely focused on the electrical and elastic responses of the films, but the response of the corresponding TiO 2 substrate is often overlooked, since it is nominally expected to be electrically passive and elastically rigid. Here, in‐operando spatiotemporal imaging of the coupled elastodynamics using X‐ray diffraction microscopy of a VO 2 film channel device on TiO 2 substrate reveals two new surprises. First, the film channel bulges during the IMT, the opposite of the expected shrinking in the film undergoing IMT. Second, a microns thick proximal layer in the substrate also coherently bulges accompanying the IMT in the film, which is completely unexpected. Phase‐field simulations of coupled IMT, oxygen vacancy electronic dynamics, and electronic carrier diffusion incorporating thermal and strain effects suggest that the observed elastodynamics can be explained by the known naturally occurring oxygen vacancies that rapidly ionize (and deionize) in concert with the IMT (MIT). Fast electrical‐triggering of the IMT via ionizing defects and an active “IMT‐like” substrate layer are critical aspects to consider in device applications.

36 MATERIALS SCIENCE↗

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗