Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Post-growth modulation doping by ion implantation

Modulation doping is a commonly adopted technique to create two-dimensional (2D) electrons or holes in semiconductor heterostructures. One constraint, however, is that the intentional dopants required for modulation doping are controlled and incorporated during the growth of heterostructures. Here, using undoped strained germanium quantum wells as the model material system, we show that modulation doping can be achieved post-growth of heterostructures by ion implantation and dopant-activation anneals. The carrier density is controlled ex situ by varying the ion fluence and implant energy, and an empirical calibration curve is obtained. While the mobility of the resulting 2D holes is lower than that in undoped heterostructure field-effect transistors built using the same material, the achievable carrier density is significantly higher. Potential applications of this modulation-doping technique are discussed.

42 ENGINEERING↗
DOE OSTI2021

Determining the electrical conductivity of metals using the 2 MA Thor pulsed power driver

We present the development of a pulsed power experimental technique to infer the electrical conductivity of metals from ambient to high energy density conditions. The method is implemented on Thor, a moderate scale (1–2 MA) pulsed power driver. The electrical conductivity of copper at elevated temperature (>4000 K) and pressure (>10 GPa) is determined, and a new tabular material model is developed, guided by density functional theory, which preserves agreement with existing experimental data. Minor modifications (<10%) are found to be necessary to the previous Lee–More–Desjarlais model isotherms in the vicinity of the melt transition in order to account for observed discrepancies with the new experimental data. An analytical model for magnetic direct drive flyer acceleration and Joule heating induced vaporization based on the Tsiolkovsky “rocket equation” is presented to assess sensitivity of the method to minor changes in electrical conductivity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Three-dimensional full-field velocity measurements in shock compression experiments using stereo digital image correlation

Shock compression plate impact experiments conventionally rely on point-wise velocimetry measurements based on laser-based interferometric techniques. This study presents an experimental methodology to measure the free surface full-field particle velocity in shock compression experiments using high-speed imaging and three-dimensional (3D) digital image correlation (DIC). The experimental setup has a temporal resolution of 100 ns with a spatial resolution varying from 90 to 200 μm/pixel. Experiments were conducted under three different plate impact configurations to measure spatially resolved free surface velocity and validate the experimental technique. First, a normal impact experiment was conducted on polycarbonate to measure the macroscopic full-field normal free surface velocity. Second, an isentropic compression experiment on Y-cut quartz–tungsten carbide assembly is performed to measure the particle velocity for experiments involving ramp compression waves. To explore the capability of the technique in multiaxial loading conditions, a pressure shear plate impact experiment was conducted to measure both the normal and transverse free surface velocities under combined normal and shear loading. The velocities measured in the experiments using digital image correlation are validated against previous data obtained from laser interferometry. Numerical simulations were also performed using established material models to compare and validate the experimental velocity profiles for these different impact configurations. The novel ability of the employed experimental setup to measure full-field free surface velocities with high spatial resolutions in shock compression experiments is demonstrated for the first time in this work.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Role of dislocations on martensitic transformation temperatures and microstructure: A molecular dynamics study

Microstructure and defects strongly affect martensitic transformations in metallic alloys. Significant progress has been made in understanding the atomic-level processes that control the role of grain boundaries and precipitates in these solid-to-solid phase transformations. Yet, the role of dislocations and their structures on martensitic transformation temperature and the resulting microstructure remains unclear. Therefore, we used large-scale molecular dynamics simulations to study the forward and reverse transformation of a martensitic material modeled after Ni63Al37 under cyclic thermal loading. The simulations reveal that dislocations in the austenite phase act as one-dimensional seeds for the martensite phase, which is present at temperatures significantly above the martensite start value. We find a reduction in the dislocation density during cyclic thermal loading, which results in the increase in martensite and austenite transition temperatures, in agreement with experiments. Importantly, we extracted a critical martensitic nuclei size for developing stable domains and found that relatively low dislocation densities are needed to grow independent martensitic variants resulting in a multi-domain structure.

Physics↗
DOE OSTI2026

Magnetic field dependence of the spin fluctuations in CeCu 5.8 ⁢Ag 0.2

Quantum phase transitions are among the most intriguing phenomena that can occur when the electronic ground state of correlated metals are tuned by external parameters such as pressure, magnetic field, or chemical substitution. Such transitions between distinct states of matter are driven by quantum fluctuations, and can give rise to macroscopically coherent phases that are at the forefront of condensed matter research. However, the nature of the critical fluctuations, and thus the fundamental physics controlling many quantum phase transitions, remain poorly understood in numerous strongly correlated metals. Here we study the model material CeCu 5.8⁢ Ag 0.2 to gain insight into the implications of critical fluctuations originating from different regions in reciprocal space. By employing an external magnetic field along the crystallographic 𝑎 and 𝑐 axis as auxiliary tuning parameter, we observe a pronounced anisotropy in the suppression of the quantum critical fluctuations, reflecting the spin anisotropy of the long-range ordered ground state at larger silver concentration. Coupled with the temperature dependence of the quantum fluctuations, these results suggest that the quantum phase transition in CeCu 5.8⁢ Ag 0.2 is driven by three-dimensional spin-density wave fluctuations.

Boraley, Xavier [Paul Scherrer Inst. (PSI), Villig↗
DOE OSTI2024

Spin-flop coupling at L⁢a 0.5 ⁢S⁢r 0.5 ⁢Fe⁢O 3 /L⁢a 0.7⁢ S⁢r 0.3 ⁢Mn⁢O 3 interfaces

Antiferromagnetic (AFM) spintronics offer several benefits compared to their ferromagnetic (FM) counterparts, such as high storage capacity and faster processing speed, however, difficulties in manipulating and detecting the AFM moments impede their implementation. Spin-flop coupling, the interfacial perpendicular coupling between FM and AFM moments, can be utilized to control the orientation of AFM moments with the application of moderate magnetic fields on the scale of tenths of a Tesla. Here, in this work, epitaxial bilayers of AFM La 0.5 Sr 0.5 FeO 3 (LSFO)/FM La 0.7 Sr 0.3 MnO 3 (LSMO) with fixed LSMO thickness (~85 u.c.) and LSFO thicknesses varying from 10 to 50 u.c. were investigated to determine the effect of Sr doping and La 1-x Sr x FeO 3 magnetocrystalline anisotropy on the strength of spin-flop coupling. X-ray magnetic linear dichroism demonstrated that the spin-flop coupling strength decreased with increasing LSFO layer thickness, persisting at a thickness of 50 u.c. (~20 nm). Furthermore, photoemission electron microscopy revealed a domain-by-domain correlation between the FM and AFM domains consistent with the perpendicular orientation dictated by spin-flop coupling. These results demonstrate that LSFO/LSMO bilayers have the potential to serve as a model materials system for AFM spin transport measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Floquet engineering of polaritonic amplification in dispersive photonic time crystals

In this study, we investigate the dynamics of dispersive photonic time crystals (PTCs) and their potential applications for controlling light-matter interaction. Using the Lorentz-Drude as dispersive material model, we analyze the effects of periodic modulation of dispersion parameters through both theoretical and numerical analysis, revealing the emergence of hybrid bandgaps from the interaction of polaritonic branches with distinctive characteristics. Our study demonstrates that dispersive PTCs offer novel excitation channels and amplification possibilities, that require lower modulation frequencies compared to nondispersive systems thus alleviating experimental challenges for the realization of PTCs in the optical regime. These findings pave the way for advancements in polaritonic lasing and resonant Raman scattering.

Classical optics↗
DOE OSTI2026

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
DOE OSTI2021

A Practical Analysis Framework for Assessment of Printed Circuit Heat Exchangers in High-Temperature Nuclear Service

Abstract Printed circuit heat exchangers (PCHEs) have high thermal efficiency because of the numerous minuscule channels. These minuscule channels result in a high thermal exchange area per unit volume, making PCHE a top contender for an intermediate heat exchanger (IHX) in high-temperature reactors. Thousands of minuscule channels make finite element analysis of the PCHE computationally infeasible. A two-dimensional analysis is usually performed for the PCHE core, which cannot simulate the local channel level responses reasonably because of the absence of global constraint influence. At present, there is no analysis technique available in the ASME Code or literature that is computationally efficient and suitable for engineers to estimate PCHE local responses. A novel but practical two-step analysis framework is proposed for performing PCHE analysis. In the first step, the channeled core is replaced by orthotropic solids with similar stiffness to simulate the global thermomechanical elastic responses of the PCHE. In the second step, local submodel analysis with detailed channel geometry and loading is performed using the elastic-perfectly plastic (EPP) material model. The proposed two-step analysis technique provides a unique capability to estimate the channel corner responses to be used for PCHE performance assessment. This study first developed a methodology for calculating the elastic orthotropic properties of the PCHE core. Next, the two-step analysis is performed for a realistic size PCHE core, and different issues observed in the results are scrutinized and resolved. Finally, a practical finite element analysis framework for PCHEs in high-temperature nuclear service is recommended.

Engineering↗
DOE OSTI2021

Nonlinear ultrasonic technique for the characterization of microstructure in additive materials

This study employs nonlinear ultrasonic techniques to track microstructural changes in additively manufactured metals. The second harmonic generation technique based on the transmission of Rayleigh surface waves is used to measure the acoustic nonlinearity parameter, β. Stainless steel specimens are made through three procedures: traditional wrought manufacturing, laser-powder bed fusion, and laser engineered net shaping. The β parameter is measured through successive steps of an annealing heat treatment intended to decrease dislocation density. Dislocation density is known to be sensitive to manufacturing variables. In agreement with fundamental material models for the dislocation-acoustic nonlinearity relationship in the second harmonic generation, β drops in each specimen throughout the heat treatment before recrystallization. Geometrically necessary dislocations (GNDs) are measured from electron back-scatter diffraction as a quantitative indicator of dislocations; average GND density and β are found to have a statistical correlation coefficient of 0.852 showing the sensitivity of β to dislocations in additively manufactured metals. Moreover, β shows an excellent correlation with hardness, which is a measure of the macroscopic effect of dislocations.

42 ENGINEERING↗
DOE OSTI2022

A cellulose-derived supramolecule for fast ion transport

Supramolecular frameworks have been widely synthesized for ion transport applications. However, conventional approaches of constructing ion transport pathways in supramolecular frameworks typically require complex processes and display poor scalability, high cost, and limited sustainability. Here, we report the scalable and cost-effective synthesis of an ion-conducting (e.g., Na + ) cellulose-derived supramolecule (Na-CS) that features a three-dimensional, hierarchical, and crystalline structure composed of massively aligned, one-dimensional, and ångström-scale open channels. Using wood-based Na-CS as a model material, we achieve high ionic conductivities (e.g., 0.23 S/cm in 20 wt% NaOH at 25 °C) even with a highly dense microstructure, in stark contrast to conventional membranes that typically rely on large pores (e.g., submicrometers to a few micrometers) to obtain comparable ionic conductivities. This synthesis approach can be universally applied to a variety of cellulose materials beyond wood, including cotton textiles, fibers, paper, and ink, which suggests excellent potential for a number of applications such as ion-conductive membranes, ionic cables, and ionotronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

GOOSE (GRAPHITE OBJECT-ORIENTED SIMULATION ENVIRONMENT)

This package is a set of extensions to MOOSE (https://moose.inl.gov) adding anisotropic graphite material models for irradiation creep and swelling. The goal of the package is to provide a framework to model graphite components for use in future high temperature nuclear reactors.

MESSNER, MARKCHRISTIAN↗
DOE OSTI2022

In-situ thermomechanical measurements during burst testing of Zr claddings using digital image correlation and infrared thermography techniques

Burst testing has long been a relied-on method for assessing fuel cladding materials in light-water reactoraccident scenarios. However, burst testing historically provided minimal in-situ observations beyond a few thermocouple readings, pressure from transducers, and post-test deformation measurements. Acquisition of full cladding thermomechanical data throughout a burst test would be invaluable to fuel system qualification, material model development and fuel performance code validation, such as the BISON fuel performance code. Recently, methodologies for simultaneously applying 2-dimensional digital image correlation and infrared thermography techniques to quantify the relationship between temperature, internal cladding pressure, and deformation during burst testing were developed. This work details results from the application of these techniques to burst testing of Zircaloy-4 cladding segments subjected to loss-of-coolant accident transients in an air environment. Two types of tests were performed with: (a) cladding segments connected to a fixed pressure reservoir as well as (b) closed systems where internal pressure increased during transients. In-situ thermomechanical measurements including strain and thermal gradient data are presented, demonstrating the effectiveness of applying combined surface measurement techniques to burst testing

Bell, Sam↗
DOE OSTI2020

Sierra/SolidMechanics 4.56 User's Guide: Addendum for Shock Capabilities

This is an addendum to the Sierra/SolidMechanics 4.56 User's Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State's International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 4.56 User's Guide should be referenced for most general descriptions of code capability and use.

97 MATHEMATICS AND COMPUTING↗