Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Materials simulation”

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 1,477 records · Page 82

Advances in the Application of Lasers in Materials Science, Chapter 5, Insights into Laser-Materials Interaction Through Modeling on Atomic and Macroscopic Scales

Computer simulations and theoretical analysis of laser-materials interactions are playing an increasingly important role in the advancement of modern laser technologies and broadening the range of laser applications. In this chapter, we first provide an overview of the current understanding of the laser coupling and transient variation of optical properties in metals, semiconductors and dielectrics, with the focus on the practical implications on the energy deposition and distribution in the irradiated targets. The continuum-level modeling of the dynamic evolution of laser-induced stresses, nonequilibrium phase transformations, and material redistribution within the laser spot are then discussed, and the need for the physical insights into the mechanisms and kinetics of highly nonequilibrium processes triggered by the laser excitation is highlighted. The physical insights can be provided by atomistic modeling, and several examples are discussed where large-scale molecular dynamics simulations are used for investigation of the mechanisms of the generation of crystal defects (vacancies, interstitials, dislocations, and twin boundaries) and the material redistribution responsible for the formation of laser-induced periodic surface structures in the single-pulse ablative regime. The need for the integrated computational approach fully accounting for the strong coupling between processes occurring at different time- and length-scales is highlighted.

Shugaev, Maxim V.↗

Additive Manufactured Materials Studies

The production of additively manufactured (AM) materials, is a growing technology that potentially offers more cost effective ways to produce metallic components for various uses. Mechanical properties of some AM metals have been studied, but more work is needed to better understand their dynamic properties as a function of orientation. Due to the melting and reheating of molten metal during fabrication, these materials tend to have internal porosity and slightly different microstructures along the weldments and normal to the weldments that may cause the materials to behave differently than their forged counterparts. The AM process builds up a part layer by layer, which provides opportunities to tailor geometry and optimize structural performance in addition to minimizing the amount of material waste for oddly shaped parts, thereby reducing cost. These advanced materials hold great promise for revolutionizing design, fabrication, and assembly of engineered components and could be used for many NNSA-specific applications including armoring complex elements of nuclear weapons, and customized diagnostic components in subcritical experiments. We have conducted materials studies in collaboration with the University of Nevada, Las Vegas (UNLV) as well as Los Alamos National Laboratory (LANL) under our Site-Directed Research and Development program and the Nevada Science Initiative. We are currently placing emphasis on quantifying how dynamic properties are affected by directionality and porosity of the material. A study of dynamic failure events in forged, layered, and additive manufactured (AM) titanium alloy subject to hypervelocity impact was conducted in collaboration with UNLV. Experiments were conducted using a two stage light-gas gun, 0.22-caliber Lexan projectiles, and different types of titanium target plates. A four-channel Photonic Doppler Velocimetry (PDV) system was used to study the free surface velocities to provide an understanding of the different failure mechanisms in these materials. The experimental measurements were used to validate computational simulations, which describe the behavior of these materials under shock. It was determined that AM, forged titanium, and multi-layered stacks produce similar velocity profiles during the early stage of impact, with the AM targets exhibiting spall at lower velocities and the multi-layered stacks exhibiting vibrations between plates. Simulations of solid and layered forged materials provide a good match to experimental data.

36 MATERIALS SCIENCE↗

Materials modeling efforts at LANL

This briefing describes some the materials theory, modeling and simulation capability at Los Alamos National Laboratory.

36 MATERIALS SCIENCE↗

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↗

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↗

Effect of electron correlation on the electronic structure and spin-lattice coupling of high- T c cuprates: Quantum Monte Carlo calculations

Electron correlation effects are particularly strong in high-temperature superconducting materials. Devising an accurate description of these materials has long been a challenge, with these strong correlation effects historically being considered impossible or impractical to simulate computationally. Using quantum Monte Carlo techniques, we have explicitly simulated electron correlations in several cuprate materials from first principles. These simulations accurately reproduce many important physical quantities of these materials, including the interaction-induced gap and superexchange coupling between copper spins, with no additional parameters beyond fundamental constants. We further investigate the dimensionless spin-lattice coupling parameter in the parent materials, showing that it varies dramatically, between 0.1 and 1.0, depending on the interlayer. Furthermore, this result indicates that the lattice and magnetic degrees of freedom are not independent in these systems, which may have ramifications for the origin of superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

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.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗