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

Pore-Space-Partition-Enabled Exceptional Ethane Uptake and Ethane-Selective Ethane–Ethylene Separation

A modell material for C 2 H 6 /C 2 H 4 separation would simultaneously have the highest C 2 H 6 uptake capacity and the highest C 2 H 6 /C 2 H 4 selectivity. But such material is elusive. A benchmark material for ethane-selective C 2 H 6 /C 2 H 4 separation is peroxo-functionalized MOF-74-Fe that exhibits the best known separation performance due to its high C 2 H 6 /C 2 H 4 selectivity (4.4), although its C 2 H 6 uptake capacity is moderate (74.3 cm 3 /g). Here, we report a family of pore-space-partitioned crystalline porous materials (CPMs) with exceptional C 2 H 6 uptake capacity and C 2 H 6 /C 2 H 4 separation potential (i.e., C 2 H 4 recovered from the mixture) despite their moderate C 2 H 6 /C 2 H 4 selectivity (up to 1.75). The ethane uptake capacity as high as 166.8 cm 3 /g at 1 atm and 298 K, more than twice that of peroxo-MOF-74-Fe, was previously achieved even though the isosteric heat of adsorption (21.9–30.4 kJ/mol) for these CPMs is as low as about one-third of that for peroxo-MOF-74-Fe (66.8 kJ/mol). While the overall C 2 H 6 /C 2 H 4 separation potentials have not yet surpassed peroxo-MOF-74-Fe, these robust CPMs exhibit outstanding properties including high thermal stability (up to 450 °C) and aqueous stability, low regeneration energy, and a high degree of chemical and geometrical tunability within the same isoreticular framework.

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

Thermodynamics-guided alloy and process design for additive manufacturing

In conventional processing, metals go through multiple manufacturing steps including casting, plastic deformation, and heat treatment to achieve the desired property. In additive manufacturing (AM) the same target must be reached in one fabrication process, involving solidification and cyclic remelting. The thermodynamic and kinetic differences between the solid and liquid phases lead to constitutional undercooling, local variations in the solidification interval, and unexpected precipitation of secondary phases. These features may cause many undesired defects, one of which is the so-called hot cracking. The response of the thermodynamic and kinetic nature of these phenomena to high cooling rates provides access to the knowledge-based and tailored design of alloys for AM. Here, we illustrate such an approach by solving the hot cracking problem, using the commercially important IN738LC superalloy as a model material. The same approach could also be applied to adapt other hot-cracking susceptible alloy systems for AM.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Active site design enables industrial scale H 2 O 2 electrosynthesis with metal-free catalysts

The electrosynthesis of hydrogen peroxide (H 2 O 2 ) via a two-electron oxygen reduction reaction enables decentralized H 2 O 2 production. While metal-free carbon catalysts are sustainable and low-cost, their performance is hindered by poorly defined active sites and uncontrolled defect states. Here, we resolve these challenges through active site design and catalyst screening using fluorine (F) and nitrogen (N) codoped carbons as model materials. Statistical analysis combined with density functional theoretical calculations reveals that F-induced structural modification and defect passivation optimize OOH* binding, with F-doping and adjacent F atoms predominantly lowering abs ΔG(OOH*). Experimental results confirm that semi-ionic C–F bonds passivate defects in nitrogen-doped carbon, enhancing catalytic activity and durability. The resulting (N, F)-codoped carbon achieves nearly 100% H 2 O 2 selectivity at 0.5–0.65 V versus the reversible hydrogen electrode and maintains > 95% across 0.01–0.65 V versus the reversible hydrogen electrode. In an electrolyzer, (N, F)-codoped carbon exhibits an H 2 O 2 yield rate of 74.35 mol g cat. −1 h -1 and sustains 300 mA cm -2 for 105 hours with ~95% faradaic efficiency. Coupling the two-electron oxygen reduction reaction with methanol oxidation further reduces cell voltage and enhances productivity. This work provides a means to design efficient catalysts for industrial H 2 O 2 electrosynthesis.

H2O2 electrosynthesis↗
DOE OSTI2020

Intrinsic chemical reactivity of solid-electrolyte interphase components in silicon–lithium alloy anode batteries probed by FTIR spectroscopy

Herein we report the solid reaction products from the chemical reaction of aprotic battery electrolyte and three purported components of the Si-based anode SEI : SiO 2 nanoparticles (NPs), lithium silicate (Li x SiO y ) powders, and Si NPs. We use FTIR and classical molecular dynamics/density functional perturbation theory to assess the solid products remaining with these model materials after exposure to electrolyte. The absence of electrochemical bias provides a view of the chemical speciation resulting from early-stage chemical reactivity during battery assembly as well as under open circuit storage conditions. We believe these species represent the initial stages of SEI growth and predict they likely drive subsequent chemical and electrochemical reactions by controlling molecular interactons at the Si active material interface. We find that nominally equivalent materials react differently even before any electrochemistry is performed (e.g., acidic SiO 2 dissolves whereas alkaline SiO 2 is relatively robust), and derive new understanding of the chemical species that could and could not form stable SEI components in Si-based anodes. These results can be used to inform how to passivate Si anode surfaces and potentially generate an artificially engineered SEI that would be stable and enable next-generation battery anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗