Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials 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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulants for Testing the Blend Can Loading System

Eleven simulants have been provided to the Savannah River National Laboratory (SRNL) Mechanical Systems & Custom Equipment Development (MSCED) organization to further test the capabilities of the Blend Can Loading System (BCLS). The basis for the materials used in these simulants are based on SRNL-STI-2020-00503, “Simulant Recommendation for the Blend Can Loading System” (Ref. i). The 11 different simulants provided to MSCED are summarized in Table ES-1. This table provides a basis of why they were selected and what properties were specifically targeted, given the composition of the procured materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2001

Low Velocity Earth-Penetration Test and Analysis

Modeling and simulation of structural impacts into soil continue to challenge analysts to develop accurate material models and detailed analytical simulations to predict the soil penetration event. This paper discusses finite element modeling of a series of penetrometer drop tests into soft clay. Parametric studies are performed with penetrometers of varying diameters, masses, and impact speeds to a maximum of 45 m/s. Parameters influencing the simulation such as the contact penalty factor and the material model representing the soil are also studied. An empirical relationship between key parameters is developed and is shown to correlate experimental and analytical results quite well. The results provide preliminary design guidelines for Earth impact that may be useful for future space exploration sample return missions.

Fasanella, Edwin L.↗
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 OSTI2021

Dynamics of laser-generated magnetic fields using long laser pulses

In this work, we report on the experimental investigation of magnetic field generation with a half-loop gold sheet coil driven by long-duration (10 ns) and high-power (0.5 TW) laser pulses. The amplitude of the magnetic field was characterized experimentally using proton deflectometry. The field rises rapidly in the first 1 ns of laser irradiation, and then increases slowly and continuously up to 10 ns during further laser irradiation. The transient dynamics of current shape were investigated with a two-dimensional (2D) numerical simulation that included Ohmic heating of the coil and the resultant change of electrical resistivity determined by the coil material temperature. The numerical simulations show rapid heating at the coil edges by current initially localized at the edges. This current density then diffuses to the central part of the sheet coil in a way that depends both on normal current diffusion as well as temporal changes of the coil resistance induced by the Ohmic heating. The measured temporal evolution of the magnetic field is compared with a model that determines a solution to the coil current and voltage that is consistent with a plasma diode model of the drive region and a 2D simulation of current diffusion and dynamic resistance due to Ohmic heating in the laser coil.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Numerical Simulations of A Conceptual MSR-EES Shoulder Recession

The present study demonstrates our in-house material response solver, Icarus's capability to simulate the ablative conditions, including pyrolysis effects due to the interaction between hypersonic boundary layers and the thermal protection system (TPS). A conceptual aeroshell shoulder design, which undergoes mission-relevant flow and material conditions, is selected for the demonstration purpose. LAURA, a structured flow solver, is used to solve flow around the shoulder at several trajectory points of a flight path. The aerothermal dataset obtained from LAURA is used to enforce boundary conditions on the aeroshell wall to simulate ablation processes. The two-layered material system is stacked with HEEET (TPS) and Aluminum (actual material). A set of base parameters, such as the angle of material orientation with respect to the flow direction, convective heat-transfer co-efficient, and aerothermal boundary conditions on the rear end of the HEEET material layer, is selected to conduct a base-case simulation. A small amount of recession was observed, indicating that the design should be fine to go through the selected flow and material conditions. Furthermore, the parameters mentioned above are individually altered and are observed to affect the shoulder design's recession compared to the base case. Verification of our Icarus setup was also carried out using a one-dimensional grid and FIAT, a one-dimensional material-response solver, to build credibility for our results. As the current work is an uncoupled fluid-material-response simulation, a local sharp mesh deformation in the recessed surface near the shoulder corner is observed.

Prakash Shrestha↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Radioactive Waste Sludge Washing and Demonstration of the Nitric-Glycolic Acid Flowsheet for Sludge Batch 10 Qualification

For each sludge batch that is processed in the Defense Waste Processing Facility (DWPF), the Savannah River National Laboratory (SRNL) performs qualification testing to demonstrate that the sludge batch (SB) is processible. During processing of SB9, DWPF will be transitioning from the Nitric-Formic Acid (NFA) flowsheet to the Nitric-Glycolic Acid (NGA) flowsheet. Thus, the qualification of SB10 was requested to only be performed using the NGA flowsheet. In order to qualify the batch for the NGA flowsheet, Sludge Receipt and Adjustment Tank (SRAT) and Slurry Mix Evaporator (SME) cycles, designated SC-19, were performed using SB10 Tank 51 sample material. SRNL received Tank 51 material in the midst of Tank Farm washing. SRNL continued the washing in the SRNL Shielded Cells. The SRNL process included the addition of Sodium Reactor Experiment (SRE) material from H Canyon, simulating the transfer of SRE from H Canyon to Tank 51 during Tank Farm washing. The washed SB10 Tank 51 material, with SRE, was characterized prior to flowsheet qualification testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1987

Integrated analysis and applications

A select overview is provided of ongoing research focusing on the development and verification of integrated structural analysis and optimal design capabilities for advanced aerospace propulsion and power systems. Subjects discussed include the following: (1) Composites - analytical models (composite mechanics), integrated computational methods, and characterization of composite structural response and durability for resin-, metal-, and ceramic-matrix systems; (2) Advanced inelastic analysis - algorithm/numerical methods for more accurate and efficient analysis; (3) Constitutive modeling - theoretical formulation and characterization of thermoviscoplastic material behavior; (4) Computational simulation - engine structures from components to assembly, and up to an entire engine system subjected to simulated test-stand and mission load histories; (5) Probabilistic structural analysis - quantification of the effects of uncertainty in geometry, material, loads, and boundary conditions on structural response for true reliability assessment; and (6) Interdisciplinary optimization - incorporation of mathematical optimization and multidisciplinary analyses to provide streamlined, autonomous optimal design systems.

Hopkins, Dale A.↗
NASA NTRSDate not supplied

PPW 2022 Short Course Introduction to Aerocapture, Entry, Descent, and Landing: Computational Materials

Computational materials science encompasses modeling and simulation techniques that are used to improve our understanding of materials from the atomic to macroscopic scales. This talk will focus on recent applications of computational materials science to ablative thermal protection materials, including heritage PICA and state-of-the-art woven materials like HEEET. Development of advanced molecular, microscopic, and mesoscopic models will be discussed.

phenolic↗
DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 MARK-I Station Blackout: Part I – Single Case Analysis

Here in this analysis, the two material interaction models available in the MELCOR code are benchmarked for a severe accident at a BWR under representative Fukushima Daiichi boundary conditions. This part of the benchmark investigates the impact of each material interaction model on accident progression through a detailed single case analysis. It is found that the eutectics model simulation exhibits more rapid accident progression for the duration of the accident. The slower accident progression exhibited by the interactive materials model simulation, however, allows for a greater degree of core material oxidation and hydrogen generation to occur, as well as elevated core temperatures during the ex-vessel accident phase. The eutectics model simulation exhibits more significant degradation of core components during the late in-vessel accident phase – more debris forms and relocates to the lower plenum before lower head failure. The larger debris bed observed in the eutectics model simulation also reaches higher temperatures, presenting a more significant thermal challenge to the lower head until its failure. At the end of the simulated accident scenario, however, core damage is comparable between both simulations due to significant core degradation that occurs during the ex-vessel phase in the interactive materials model simulation. A key difference between the two models’ performance is the maximum temperatures that can be reached in the core and therefore the maximum ΔT between any two components. When implementing the interactive materials model, users have the option to modify the liquefaction temperature of the ZrO 2 -interactive and UO 2 -interactive materials as a way to mimic early fuel rod failure due to material interactions. Through modification of the liquefaction of high melting point materials with significant mass, users may inadvertently limit maximum core temperatures for fuel, cladding, and debris components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Equivalent Fundamental-Mode Source Simulations for Spherical Uranium and Plutonium Systems [Abstract]

The equivalent fundamental-mode source refers to a “source” which is identically distributed in space, energy, and angle as that of a fundamental-mode fission source distribution. It has previously been shown that the system fixed source multiplication ($M_{fs}$) relates to the effective multiplication factor ($\kappa_{eff}$) as $M_{fs} = \frac{g^*}{1-k_{eff}}$ The term $\mathcal{g}$* therefore relates the system source distribution to the system equivalent fundamental-mode source. As previously shown, the $\mathcal{g}$* term is approximately 1 for deeply subcritical systems and then diverges as the system multiplication increases. For a system with a point source located in the center of an assembly (such a sphere), $\mathcal{g}$* will be greater than 1; this happens because there is a higher probability of induced fission when all of the starter neutrons originate from the center of the sphere. For a uniform source, $\mathcal{g}$* is less than 1, but does not diverge from 1 as dramatically as seen by a point source. One reason that $\mathcal{g}$* is an important parameter is that it is used in measurement methods to infer reactor kinetics parameters ($β_{eff}$ being one example). This work will discuss multiple ways to simulate $\mathcal{g}$* using MCNP®6.2, which include use of traditional methods versus a single input file method that was recently published. A discussion related to methods to simulate $M_{fs}$ will be presented. Spherical systems of Highly Enriched Uranium (HEU) will be compared with previous works. In addition, results of plutonium including multiple reflector materials will be investigated. Simulated systems will include both parametric studies as well as configurations used in recent experiments performed at the National Criticality Experiments Research Center (NCERC). These studies will include subtleties associated with system geometry that have not been previously shown. Last, this work will investigate the relationship between $\mathcal{g}$* and reactor kinetics parameters (such as $β_{eff}$).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Disorder effects in nitride semiconductors: impact on fundamental and device properties

Abstract Semiconductor structures used for fundamental or device applications most often incorporate alloy materials. In “usual” or “common” III–V alloys, based on the InGaAsP or InGaAlAs material systems, the effects of compositional disorder on the electronic properties can be treated in a perturbative approach. This is not the case in the more recent nitride-based GaInAlN alloys, where the potential changes associated with the various atoms induce strong localization effects, which cannot be described perturbatively. Since the early studies of these materials and devices, disorder effects have indeed been identified to play a major role in their properties. Although many studies have been performed on the structural characterization of materials, on intrinsic electronic localization properties, and on the impact of disorder on device operation, there are still many open questions on all these topics. Taking disorder into account also leads to unmanageable problems in simulations. As a prerequisite to address material and device simulations, a critical examination of experiments must be considered to ensure that one measures intrinsic parameters as these materials are difficult to grow with low defect densities. A specific property of nitride semiconductors that can obscure intrinsic properties is the strong spontaneous and piezoelectric fields. We outline in this review the remaining challenges faced when attempting to fully describe nitride-based material systems, taking the examples of LEDs. The objectives of a better understanding of disorder phenomena are to explain the hidden phenomena often forcing one to use ad hoc parameters, or additional poorly defined concepts, to make simulations agree with experiments. Finally, we describe a novel simulation tool based on a mathematical breakthrough to solve the Schrödinger equation in disordered potentials that facilitates 3D simulations that include alloy disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2014

A Large Lunar Surface Testbed from Low Cost Material

For users needing to simulate the lunar surface, several distinct avenues have been used. Numerous volcanic areas, including Hawaii, have been used. While providing very large areas and scenic interest, field parties to such an area is expensive and limits testing time. An alternative is to build test facilities locally. This has been done many ways, contrast GRC-1, GSC-1, BP-1 and the KSC Morpheus facility [1-4]. GRC-1 is a mixture of sand and clay; GSC-1 and BP-1 are waste materials created in the process of crushing basaltic rock. The Morpheus field used salvaged concrete and crushed quartz rock [5]. Here I report about a 30 m X 30 m test area at MSFC which was both low cost and relatively high fidelity [6].

Rickman, Douglas↗
DOE OSTI2023

Frustration- and doping-induced magnetism in a Fermi–Hubbard simulator

Geometrical frustration in strongly correlated systems can give rise to a plethora of novel ordered states and intriguing magnetic phases, such as quantum spin liquids. Promising candidate materials for such phases can be described by the Hubbard model on an anisotropic triangular lattice, a paradigmatic model capturing the interplay between strong correlations and magnetic frustration. However, the fate of frustrated magnetism in the presence of itinerant dopants remains unclear, as well as its connection to the doped phases of the square Hubbard model. Here we investigate the local spin order of a Hubbard model with controllable frustration and doping, using ultracold fermions in anisotropic optical lattices continuously tunable from a square to a triangular geometry. At half-filling and strong interactions U/t ≈ 9, we observe at the single-site level how frustration reduces the range of magnetic correlations and drives a transition from a collinear Néel antiferromagnet to a short-range correlated 120° spiral phase. Away from half-filling, the triangular limit shows enhanced antiferromagnetic correlations on the hole-doped side and a reversal to ferromagnetic correlations at particle dopings above 20%, hinting at the role of kinetic magnetism in frustrated systems. Furthermore, this work paves the way towards exploring possible chiral ordered or superconducting phases in triangular lattices and realizing t–t' square lattice Hubbard models that may be essential to describe superconductivity in cuprate materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

CRADA NO. NFE-22-00001 Between CNS and Northrop Grumman Systems Corporation

Consolidated Nuclear Security, LLC (CNS) and Northrop Grumman Systems Corporation (NGC) have established a partnership through the subject Cooperative Research and Development Agreement (CRADA) in three areas of collaboration: (1) Additive Manufacturing (AM) of manufactured carbon-carbon composite materials, (2) modeling and simulation, and (3) specialized chip and circuit board manufacturing for supply chain assurance.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Simulating AOTV heating environment in an arc jet

The application of arc-jet facilities to the simulation of the combined radiative and convective heating environments calculated for AOTV return missions is presented. The AOTV mission requirement and design configuration determine the magnitude of the radiative and convective heating rates. This paper considers methods for obtaining radiation from the high temperature plasma in arc-jet flows. A computer code that predicts arc-heater performance was used to determine the effect of gas mixtures on arc-heater performance, with the intent of increasing the gas temperature and resultant radiative emission. Numerical radiation solutions were compared with direct measurements of the radiative intensity from the gas cap. Convective heating measurements were made as well. The arc-jet was operated at AOTV enthalpy and pressure levels for these measurements to assess the simulation capabilities for AOTV material development/testing. A parametric study is presented that defines the capability of the Interaction Heating Facility arc jet at Ames to simulate both the radiative and the convective heating for current AOTV configurations.

Balboni, J.↗