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

Progressive Damage Analyses of Skin/Stringer Debonding

The debonding of skin/stringer constructions is analyzed using a step-by-step simulation of material degradation based on strain softening decohesion elements and a ply degradation procedure. Decohesion elements with mixed-mode capability are placed at the interface between the skin and the flange to simulate the initiation and propagation of the delamination. In addition, the initiation and accumulation of fiber failure and matrix damage is modeled using Hashin-type failure criteria and their corresponding material degradation schedules. The debonding predictions using simplified three-dimensional models correlate well with test results.

Daville, Carlos G.↗
NASA NTRS2005

Present Status of Lunar Regolith Simulants, Workshop Overview and Objectives

As NASA embarks on a new phase of Space Exploration aiming at establishing a long- term human presence on the Moon and Mars, we face some immediate challenges in the development of tbe technologies that will take us there. Among these challenges, the lack of accepted reference materials that simulate the lunar soil and regolith is a significant one that must be addressed now. The workshop organized today focuses on defining these materials for the use of NASA and other organizations involved in research and development efforts of lunar technologies. The range of these technologies is vast and so is the variety of lunar materials. The gathering of experts from a large number of disciplines is a pre-requisite to achieve a credible consensus on this issue and be able to make the right recommendations for NASA management as they consider solutions to the development and usage of these specific materials.

Sibille, L.↗
DOE OSTI2021

Surface prediction and measurement for modulated tool path (MTP) turning

Here, this paper describes a time-domain simulation for predicting surface finish in modulated tool path (MTP) turning, which uses sinusoidal axis motions in the feed direction to produce discontinuous chips for ductile workpiece materials. The simulation includes: the low frequency and low amplitude tool oscillation in the feed direction; the time-varying chip thickness, cutting force, and tool displacement; and the plastic side flow effect used to calibrate the effective tool nose radius. Comparisons between predicted and measured surface profiles are presented for turning a 6061-T6 aluminum cylinder as a function of the MTP oscillation frequency and amplitude with discontinuous chip formation.

42 ENGINEERING↗
NASA NTRS1976

Flight simulation testing equipment for composite material systems

A test program is discussed which aims at establishing the time-temperature-stress characteristics of several classes of high-temperature composite materials in order to determine their suitability for applications in supersonic cruise aircraft. Five advanced composite materials (a boron epoxy, a boron polyimide, a graphite epoxy, a graphite polyimide, and diffusion-bonded boron aluminum) are being evaluated using a flight-test simulator capable of long-term automatic testing based on random loading and realistic flight temperature profiles. The design, construction, and checkout of this simulator are described along with the digital load programmer, load magnitude controllers, the hydraulic pumping system, the heating and cooling systems, the control console, and the data recording system. Typical results for short-term tests performed at constant temperatures and accelerated load rates are presented in terms of a random-load spectrum and a wearout model.

Haskins, J. F.↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the Breakdown Behavior of Ultrawide-Bandgap Boron Nitride Power Diodes Using Device Modeling

Herein, a device study using technology computer-aided design simulation to theoretically analyze the electrical performance of ultrawide-bandgap boron nitride (BN)-based vertical junction devices is performed, including h-BN Schottky diode, h-BN pn diode, and h-BN/AlN pn diode; this is also the first demonstration of the BN power devices in simulation. The material properties of BN are defined with recently reported data, and the physical mechanisms of the device performance are systematically investigated. Additionally, the h-BN junctions in this simulation shows excellent performance, especially for breakdown behaviors. Schottky diode shows a turn-on voltage of 0.6 V for Pt Schottky contact and breakdown voltages over 450 V for 5 μm, 6 × 10 15 cm –3 p-type-doped drift layer; The h-BN pn diode shows a turn-on voltage of 6 V and breakdown voltages over 3 kV with a critical electric field of 13.6 MV cm –1 for 2.5 μm, 2 × 10 16 cm –3 p-type-doped drift layer. The h-BN/AlN heterojunction pn diode shows a turn-on voltage of 5.8 V and breakdown voltage over 2 kV for 2.5 μm, 2 × 10 16 cm –3 n-type-doped AlN drift layer. Herein, an understanding of the device principles of vertical BN junctions is provided, which can serve as a reference for the future development of robust BN power electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2021

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗
DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1977

Laboratory simulation of infrared astrophysical features

Intermediate resolution emission spectroscopy was used to study a group of 9 terrestrial silicates, 1 synthetic silicate, 6 meteorites and 2 lunar soils; comparisons were made with the intermediate resolution spectra of Comet Kohoutek in order to determine which materials best simulate the 10um astrophysical feature. Mixtures of silicates which would yield spectra matching the spectrum of the comet in the 10um region include: (1) A hydrous layer lattice silicate in combination with a high temperature condensate; (2) an amorphous magnesium silicate in combination with a high temperature condensate and (3) glassy olivine and glassy anorthite in approximately equal proportions.

Rose, L. A.↗
NASA NTRS1992

Current research in oxidation-resistant carbon-carbon composites at NASA. Langley Research Center

The significant potential of carbon-carbon composites for high-temperature structural applications is well established. For hypersonic vehicle applications, desirable properties include low density, high specific strength and stiffness, low coefficients of thermal expansion, and retention of mechanical properties above 3000 F. A significant problem associated with carbon materials, however, is that they oxidize rapidly in air at temperatures above about 800 F, and therefore must be protected from oxidation. Successful development of effective methods of oxidation protection is key to the eventual utilization of carbon-carbon composites on hypersonic vehicles such as NASP. In this presentation, the basic elements of an oxidation-protection system are described. Results from oxidation-performance evaluations of state-of-the-art ACC-4 type material in simulated airframe vehicle environments (temperature, pressure, and time) conducted at NASA Langley are also presented. NASA Langley has an active research effort to improve the oxidation resistance of carbon-carbon materials for airframe structural and vehicle thermal protection applications. Conversion coating and sealant development research is highlighted.

Ohlhorst, Craig W.↗
NASA NTRSDate not supplied

Determination of Interfacial Energy for Solid-State Precipitation

Precipitation of strengthening phases from solid state is a critical process for many materials. For example, the high-temperature creep resistance of Ni-based superalloys is governed by ’ precipitates. However, the interfacial energy which controls the precipitation is expensive to determine experimentally. For simulations, the atomic nature of precipitation events necessitates using Molecular Dynamics (MD) so relying on fitting to a-priori experimental behavior is not required. However, the solid-state precipitation takes orders of magnitude longer time than is accessible for traditional MD simulation. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation (and determination of the interfacial energy). Implementing these methods in LAMMPS, we calculated the interfacial energy for the ’ precipitation out of a super-saturated  phase in Ni-based superalloys. The resulting interfacial energy value for the Ni-Al system is well within the accuracy of experimental values, demonstrating the capability of this method. Enabling the calculation of the interfacial energy without a-priori experimental knowledge of the precipitation unlocks a more efficient ICME-based approach for guiding annealing protocols in Ni-based superalloys and other solid-state precipitation systems.

Precipitation↗
NASA NTRS1979

Laboratory simulation of infrared astrophysical features

Laboratory infrared emission and absorption spectra have been taken of terrestrial silicates, meteorites, and lunar soils in the form of micrometer and submicrometer grains. The emission spectra were taken in a way that imitates telescopic observations. The purpose was to see which materials best simulate the 10-micron astrophysical feature. The emission spectra of dunite, fayalite, and Allende give a good fit to the 10-micron broadband emission feature of comets Bennett and Kohoutek. A study of the effect of grain size on the presence of the 10-micron emission feature of dunite shows that for particles larger than 37 microns no feature is seen. The emission spectrum of the Murray meteorite, a Type 2 carbonaceous chrondrite, is quite similar to the intermediate-resolution spectrum of comet Kohoutek in the 10-micron region. Hydrous silicates or amorphous magnesium silicates in combination with high-temperature condensates, such as olivine or anorthite, would yield spectra that match the intermediate-resolution spectrum of comet Kohoutek in the 10-micron region. Glassy olivine and glassy anorthite in approximately equal proportions would also give a spectrum that is a good fit to the cometary 10-micron feature.

Rose, L. A.↗
NASA NTRS2009

An Overview of the State of the Art in Atomistic and Multiscale Simulation of Fracture

The emerging field of nanomechanics is providing a new focus in the study of the mechanics of materials, particularly in simulating fundamental atomic mechanisms involved in the initiation and evolution of damage. Simulating fundamental material processes using first principles in physics strongly motivates the formulation of computational multiscale methods to link macroscopic failure to the underlying atomic processes from which all material behavior originates. This report gives an overview of the state of the art in applying concurrent and sequential multiscale methods to analyze damage and failure mechanisms across length scales.

Saether, Erik↗
NASA NTRS1989

Modifications of comet materials by the sublimation process: Results from simulation experiments

An active comet like comet Halley loses by sublimation a surface layer of the order of 1 m thickness per perihelion passage. In situ measurements show that water ice is the main constituent which contributes to the gas emission although even more volatile species (CO, NH3, CH4, CO2 etc.) have been identified. Dust particles which were embedded in the ices are carried by the sublimating gases. Measurements of the chemical composition of cometary grains indicate that they are composed of silicates of approximate chondritic composition and refractory carbonaceous material. Comet simulation experiments show that significant modifications of cometary materials occur due to sublimation process in near surface layers which have to be taken into account in order to derive the original state of the material.

Gruen, E.↗
NASA NTRS2007

A Method to Perform Direct Oxygen Analysis on Lunar Simulants and Other Complex Oxide Materials

An essential requirement for making space travel and long term missions more efficient and affordable to NASA includes finding innovative ways to supply oxygen for life support and propulsion. In this experiment, carrier gas hot extraction was investigated as a possible method for measuring the oxygen from samples of lunar soil simulants before and after oxygen extraction. The determination of oxygen using the R0600 Oxygen Determinator is usually limited to oxides with low oxygen concentrations, but after the manipulation of certain furnace operating parameters such as analysis time and ramp rate, the R0600 was used to determine the oxygen content of high concentration oxides such as Fe 2O3 , Al2O3 , and SiO2.

Santiago-Maldonado, Edgardo↗
NASA NTRS1986

Simulation studies for surfaces and materials strength

During this reporting period three investigations were carried out. The first area of research concerned the analysis of the structure-energy relationship in small clusters. This study is very closely related to the improvement of the potential energy functions which are suitable and simple enough to be used in atomistic simulation studies. Parameters obtained from ab initio calculations for dimers and trimers of Al were used to estimate energetics and global minimum energy structures of clusters continuing up to 15 Al atoms. The second research topic addressed modeling of the collision process for atoms impinging on surfaces. In this simulation study qualitative aspects of the O atom collision with a graphite surface were analyzed. Four different O/graphite systems were considered and the aftermath of the impact was analyzed. The final area of investigation was related to the simulation of thin amorphous Si films on crystalline Si substrates. Parameters obtained in an earlier study were used to model an exposed amorphous Si surface and an a-Si/c-Si interface. Structural details for various film thicknesses were investigated at an atomistic level.

Halicioglu, T.↗
DOE OSTI2026

Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems

Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.

Chen, I-Chi [Ames Laboratory (AMES), Ames, IA (Uni↗