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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Feasibility Study of SSTO Base Heating Simulation in Pulsed-Type Facilities

A laboratory simulation of the base heating environment of the proposed reusable Single-Stage-To-Orbit vehicle during its ascent flight was proposed. The rocket engine produces CO2 and H2, which are the main combustible components of the exhaust effluent. The burning of these species, known as afterburning, enhances the base region gas temperature as well as the base heating. To determine the heat flux on the SSTO vehicle, current simulation focuses on the thermochemistry of the afterburning, thermophysical properties of the base region gas, and ensuing radiation from the gas. By extrapolating from the Saturn flight data, the Damkohler number for the afterburning of SSTO vehicle is estimated to be of the order of 10. The limitations on the material strengths limit the laboratory simulation of the flight Damkohler number as well as other flow parameters. A plan is presented in impulse facilities using miniature rocket engines which generate the simulated rocket plume by electric ally-heating a H2/CO2 mixture.

Park, Chung Sik↗
NASA NTRS2006

Development of Design Analysis Methods for C/SiC Composite Structures

The stress-strain behavior at room temperature and at 1100 C (2000 F) was measured for two carbon-fiber-reinforced silicon carbide (C/SiC) composite materials: a two-dimensional plain-weave quasi-isotropic laminate and a three-dimensional angle-interlock woven composite. Micromechanics-based material models were developed for predicting the response properties of these two materials. The micromechanics based material models were calibrated by correlating the predicted material property values with the measured values. Four-point beam bending sub-element specimens were fabricated with these two fiber architectures and four-point bending tests were performed at room temperature and at 1100 C. Displacements and strains were measured at various locations along the beam and recorded as a function of load magnitude. The calibrated material models were used in concert with a nonlinear finite element solution to simulate the structural response of these two materials in the four-point beam bending tests. The structural response predicted by the nonlinear analysis method compares favorably with the measured response for both materials and for both test temperatures. Results show that the material models scale up fairly well from coupon to subcomponent level.

Sullivan, Roy M.↗
NASA NTRS2009

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗
DOE OSTI2022

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

A Functional Comparison of Lunar Regoliths and Their Simulants

Lunar regolith simulants are essential to the development of technology for human exploration of the Moon. Any equipment that will interact with the surface environment must be tested with simulant to mitigate risk. To reduce the greatest amount of risk, the simulant must replicate the lunar surface as well as possible. To quantify the similarities and differences between simulants, the Figures of Merit were developed. The Figures of Merit software compares the simulants and regolith by particle size, particle shape, density, and bulk chemistry and mineralogy; these four properties dictate the majority of the remaining characteristics of a geologic material. There are limitations to both the current Figures of Merit approach and simulants in general. The effect of particle textures is lacking in the Figures of Merit software, and research into this topic has only recently begun with applications to simulants. In addition, not all of the properties for lunar regolith are defined sufficiently for simulant reproduction or comparison; for example, the size distribution of particles greater than 1 centimeter and the makeup of particles less than 10 micrometers is not well known. For simulants, contamination by terrestrial weathering products or undesired trace phases in feedstock material is a major issue. Vapor deposited rims have not yet been created for simulants. Fortunately, previous limitations such as the lack of agglutinates in simulants have been addressed and commercial companies are now making agglutinate material for simulants. Despite some limitations, the Figures of Merit sufficiently quantify the comparison between simulants and regolith for useful application in lunar surface technology. Over time, the compilation and analysis of simulant user data will add an advantageous predictive capability to the Figures of Merit, accurately relating Figures of Merit characteristics to simulant user parameters.

Rickman, D.↗
DOE OSTI2023

X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission

Abstract In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor or other means of deflection proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's X-ray energy deposition into the asteroid and the resulting material ablation. The relevant physics in these simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed X-ray energy deposition model developed using the radiation-hydrodynamics code Kull that can be used to initiate a nuclear mitigation mission calculation. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials at a range of porosities, two different source spectra, and a broad range of radiation fluences, source durations, and angles of incidence. Using blowoff momentum as the primary metric, the model-initiated simulation results match the full radiation-hydrodynamics results to within 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2017

Analyzing the Chemical and Spectral Effects of Pulsed Laser Irradiation to Simulate Space Weathering of a Carbonaceous Chondrite

Space weathering processes alter the chemical composition, microstructure, and spectral characteristics of material on the surfaces of airless bodies. The mechanisms driving space weathering include solar wind irradiation and the melting, vaporization and recondensation effects associated with micrometeorite impacts e.g., [1]. While much work has been done to understand space weathering of lunar and ordinary chondritic materials, the effects of these processes on hydrated carbonaceous chondrites is poorly understood. Analysis of space weathering of carbonaceous materials will be critical for understanding the nature of samples returned by upcoming missions targeting primitive, organic-rich bodies (e.g., OSIRIS-REx and Hayabusa 2). Recent experiments have shown the spectral properties of carbonaceous materials and associated minerals are altered by simulated weathering events e.g., [2-5]. However, the resulting type of alteration i.e., reddening vs. bluing of the reflectance spectrum, is not consistent across all experiments [2-5]. In addition, the microstructural and crystal chemical effects of many of these experiments have not been well characterized, making it difficult to attribute spectral changes to specific mineralogical or chemical changes in the samples. Here we report results of a pulsed laser irradiation experiment on a chip of the Murchison CM2 carbonaceous chondrite to simulate micrometeorite impact processing.

Thompson, M. S.↗
DOE OSTI2022

Crystal plasticity including a phase-field deformation twinning model for the high-rate deformation of cyclotetramethylene tetranitramine

A new finite strain thermomechanical model for the high-rate deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) has been developed and applied to simulations of plate impact experiments. The crystal plasticity model is based on a model developed previously for RDX (Luscher et al. 2017), which is extended to incorporate deformation twinning. Twinning during normal plate impacts is simulated with a phase-field twin model. First, material parameters governing the kinetics of dislocation slip are calibrated on the subset of simulations which had negative Schmid factors for the twin system. Second, a parametric study of the twin material parameters was performed to find suitable values. Overall, the results of the simulations with the phase-field twinning model are reported for impacts on several crystal orientations. We find that the twin growth decreases with increasing distance from the impact surface because of dissipation of the shock front via dislocation-mediated plasticity, and that the simulated interface velocity with and without phase-field twinning do not show appreciable differences. These modeling results suggest that the significance of twinning in β-HMX cannot be determined with traditional loading configurations and diagnostics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal–Organic Framework

In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.

adsorption capacities↗
DOE OSTI2023

Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal

Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Etch Profile Simulation Using Level Set Methods

Etching and deposition of materials are critical steps in semiconductor processing for device manufacturing. Both etching and deposition may have isotropic and anisotropic components, due to directional sputtering and redeposition of materials, for example. Previous attempts at modeling profile evolution have used so-called "string theory" to simulate the moving solid-gas interface between the semiconductor and the plasma. One complication of this method is that extensive de-looping schemes are required at the profile corners. We will present a 2D profile evolution simulation using level set theory to model the surface. (1) By embedding the location of the interface in a field variable, the need for de-looping schemes is eliminated and profile corners are more accurately modeled. This level set profile evolution model will calculate both isotropic and anisotropic etch and deposition rates of a substrate in low pressure (10s mTorr) plasmas, considering the incident ion energy angular distribution functions and neutral fluxes. We will present etching profiles of Si substrates in Ar/Cl2 discharges for various incident ion energies and trench geometries.

Hwang, Helen H.↗
NASA NTRS2016

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗
DOE OSTI2023

Multiscale Development and Validation of the Stainless Steel Alloy Corrosion (SStAC) Tool for High Temperature Engine Materials

The goal of this project was to create a simulation tool to assist in part design, reducing costs by eliminating the need for conservative material selection, and enabling alloy optimization to improve corrosion resistance. Thus, it helps decrease costs for existing engines and assists in the material and engine design for the engines of the future. To meet our stated goal, we developed the Stainless-Steel Alloy Corrosion (SStAC) tool, which: (1) Models corrosion of engine valves fabricated from 21-2N valve steel in an engine environment at temperatures up to 800 °C. (2) Runs in 1D for a fast estimate of the corrosion rate, and in 2D or 3D for more detailed simulations that represent any part geometry and predict the precise location and rate of corrosion and its impact on the mechanical and thermal behavior. (3) Was implemented using the open-source Multiphysics Object Oriented Simulation Environment (MOOSE). (4) Couples the corrosion model with mechanics and thermal transport models. The SStAC tool can predict the corrosion rate of 21-2N engine valves with less than 10% error.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiative Heat Transfer and 2D Transition Metal Dichalcogenide Materials

Here we study the radiative heat transfer power in the family of transition metal dichalcogenide monolayers in their H- and T-symmetries. For this purpose, the electronic and optical properties computed from first-principles are used in effective models to understand the emerging scaling laws for metals and semiconductors as well as specific material signatures as control knobs for radiative heat transfer. Our combined approach of analytical modeling with properties from ab initio simulations can be used for other materials families to build a materials database for radiative heat transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles simulation of light-ion microscopy of graphene

The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10–200 keV light ions in free-standing graphene. Here we predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to three times higher signal and up to five times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. This predicted superiority of exit-side emission likely arises from anisotropic kinetic emission. The charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Organics Synthesized Using Iron-Grain Silicates

We use Fischer-Tropsch type (FTT) synthesis to produce hydrocarbons by hydrogenating carbon monoxide via catalytic reactions. The products of these reactions have been studied using 'natural' catalysts and calculations of the efficiency of FTT synthesis in the Solar Nebula suggest that these types of reactions could make significant contributions to the composition of material near three AU. We coat Fe-silicate grains with organic material using FTT synthesis to simulate the chemistry in the early Solar Nebula. In our experimental setup, we roughly model a nebular environment where grains are successively transported from hot to cold regions of the nebula. In other words, the starting gases and FTT products are continuously circulated through the grains at high temperature with intervals of cooling. Organics generated in this manner could represent the carbonaceous material incorporated in comets and meteorites. We analyze the resulting organics and present the results.

Johnson, N. M.↗
NASA NTRSDate not supplied

Simulation of Dust-Laden Flow in the DLR L2K Facility

An integrated CFD-particle trajectory code called US3D-DUST is used to simulate the dusty flow environments performed in the German Aerospace Center (DLR) L2K arcjet facility. The particle trajectory code uses a Lagrangian approach, where the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The L2K experiments measured mass loss to flat-face shaped samples of the Norcoat-Liège thermal protection system material. The US3D-DUST simulations were able to reproduce particle velocity and stagnation point heating rate data. When applied to the L2K dusty flow conditions, the US3D-DUST predictions of mass loss due to particle impacts and surface recession rates over-predicted the experimental values. It is recommended to perform additional dusty flow tests in the L2Kto provide and improve the surface damage model for Norcoat-Liège.

Thermal protection systems↗