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 739 records · Page 41

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

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↗

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

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

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

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↗

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

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↗

Probabalistic Risk Analysis and Thermal Margin Process for an Inflatable Aeroshell

Uncertainties always exist in atmospheric entry aeroheating environments and the thermal response of thermal protection system (TPS) material. These uncertainties are mitigated in the design by ap-plying margin and factors of safety to the TPS. Entry vehicle TPS is often conservatively over-sized for the heat loads that are experienced along the entry trajectory by stacking worst-case scenarios together. Additionally, the current TPS design and margin process used by NASA offers very little insight into the risk of over-temperature during flight and the reliability of the heat shield performance [1,3]. A probabilistic margin process can be used to calculate the amount of TPS margin necessary to survive a given entry heat load at a specified level of risk [2,3,4]. The vehicle’s initial entry state (entry velocity, flight path angle, and entry mass) determines the expected atmospheric entry environmental conditions and resulting heat load that the entry vehicle will experience. If there is flexibility in the entry state, then this process can be used to select an appropriate combination of entry state parameters and TPS size to target a desired reentry reliability. This probabilistic margin process allows engineers to make informed aeroshell design, entry-trajectory design, and TPS performance risk trades while preventing excessive TPS margin from being applied. The probabilistic TPS margin process has been performed to determine TPS thickness and entry heating constraints given an acceptable risk level for the Low Earth Orbit Flight Experiment of an Inflatable Decelerator (LOFTID) flight project. The process is used in a manner to size the entry heat load for a given flexible TPS (FTPS) thickness so that it meets project reliability standards while allowing the FTPS and the underlying inflatable structure (IS) to be pushed to adequately high temperatures. Since the LOFTID project is an experimental flight demonstration, it is de-sired to drive the FTPS and IS to temperatures that cover a large range of their thermal response models’ applicability. This will allow the thermal response models to be better improved and validated post-flight using LOFTID’s extensive instrumentation embedded within the aeroshell. The presentation demonstrates how uncertainty analysis is carried out using an end-to-end Monte Carlo process where three separate Monte Carlo simulations are run in sequence. The first Monte Carlo simulation operates on the entry trajectory model to generate trajectory parameter dispersions that are fed into the second Monte Carlo simulation. The second Monte Carlo simulation operates on the aerothermodynamics model to generate aeroheating parameter dispersions that are fed into the third Monte Carlo simulation. The third Monte Carlo simulation operates on the FTPS material thermal response model to generate the final FTPS/IS thermal response dispersions. The end-to-end Monte Carlo simulation propagates the uncertainties of each model into the next to quantify the resulting uncertainty of the FTPS/IS thermal response. The fractional contributions of the uncertain parameters in the trajectory, aerothermal, and thermal response models to the variance in the FTPS/IS thermal response is determined as a byproduct of the Monte Carlo analysis. The structural uncertainty of the FTPS thermal response model is evaluated by flight relevant ground testing and model error analysis using test measurements. This probabilistic TPS margin process had never been applied to an entry vehicle and it is one of the LOFTID project’s goals to demonstrate its merits.

Steven A. Tobin↗

High gas velocity burner tests on silicon carbide and silicon nitride at 1200 C

Ten SiC materials and five Si3N4 materials were exposed in a Mach 1-gas-velocity burner simulating a gas-turbine engine environment. All materials studied are commercially available. Cyclic tests up to 100 hours' duration were conducted at specimen temperatures of 1200 C. A specimen geometry was used that develops thermal stresses during thermal cycling in a manner similar to blades and vanes of a gas turbine engine. Materials were compared on a basis of weight change, dimensional reductions, metallography, fluorescent-penetrant inspection, X-ray diffraction analyses, failure mode, and general appearance. One hot-pressed SiC, one reaction-sintered SiC, and three hot-pressed Si3N4 materials survived the program goal of 100 one-hour cycles of exposures. Of the materials that failed to meet the program goal, thermal fatigue was identified as the exclusive failure mode.

Sanders, W. A.↗

Initial Investigation of Cryogenic Wind Tunnel Model Filler Materials

Filler materials are used for surface flaws, instrumentation grooves, and fastener holes in wind tunnel models. More stringent surface quality requirements and the more demanding test environment encountered by cryogenic wind tunnels eliminate filler materials such as polyester resins, plaster, and waxes used on conventional wind tunnel models. To provide a material data base for cryogenic models, various filler materials are investigated. Surface quality requirements and test temperature extremes require matching of coefficients of thermal expansion or interfacing materials. Microstrain versus temperature curves are generated for several candidate filler materials for comparison with cryogenically acceptable materials. Matches have been achieved for aluminum alloys and austenitic steels. Simulated model surfaces are filled with candidate filler materials to determine finishing characteristics, adhesion and stability when subjected to cryogenic cycling. Filler material systems are identified which meet requirements for usage with aluminum model components.

Firth, G. C.↗

Coupled Aerosol-Chemistry Simulations of the January 2022 Eruptions of Hunga Tonga-Hunga Ha’apai in the NASA GEOS Earth System Model

The January 2022 eruptions of the underwater Hunga Tonga-Hunga Ha’apai volcano injected more than 100 Tg of water vapor and approximately 0.5 Tg of sulfur dioxide into the stratosphere. Injected materials reached as high as ~55 km in altitude, while the main plume from the eruption travelled west over Australia and the Indian Ocean between about 20 – 30 km altitude. We investigate the transport and the impact of the erupted materials in coupled aerosol-chemistry simulations performed with the NASA Goddard Earth Observing System (GEOS) Earth system model. A sulfur mechanism introduced into the Global Modeling Initiative (GMI) stratospheric-tropospheric chemistry package allows for an interactive simulation of the water vapor-chemistry impacts, and the large water vapor perturbation results in rapid conversion of sulfur dioxide to sulfate aerosol in our aerosol mechanism. We report on the results of a multi-year ensemble of simulations performed with this system that include a control ensemble (no eruption), a water vapor-only injection ensemble, and a water vapor and sulfur dioxide injection ensemble. Simulations of the near-field aerosol and chemistry transport are compared to available measurements from OMPS-LP, OMPS-NM, CALIOP, MLS, and other sensors. The extended impact of the volcanic materials on the stratospheric composition and chemistry over the next several years is further investigated in forecast simulations.

Peter Colarco↗

Experimental performance of an ablative material as an external insulator for a hypersonic research aircraft

An ablative material composed of silica-filled elastomeric silicone was tested to evaluate its thermal and structural performance as an external insulator, or heat shield, for a hypersonic research aircraft. The material was also tested to determine whether it would form a durable char layer when initially heated and thereafter function primarily as an insulator with little further pyrolysis or char removal. Aerothermal tests were representative of nominal Mach 6 cruise conditions of the aircraft, and additional tests were representative of Mach 8 cruise and interference heating conditions. Radiant heating tests were used to simulate the complete nominal Mach 6 surface-temperature history. The silica char that formed during aerothermal tests was not durable. The char experienced a general and preferential surface recession, with the primary mechanism for char removal being erosion. Tests revealed that radiant heating is not a valid technique for simulating aerodynamic heating of the material.

Puster, R. L.↗

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

Fatigue-test acceleration with flight-by-flight loading and heating to simulate supersonic-transport operation

Possibilities for reducing fatigue-test time for supersonic-transport materials and structures were studied in tests with simulated flight-by-flight loading. In order to determine whether short-time tests were feasible, the results of accelerated tests (2 sec per flight) were compared with the results of real-time tests (96 min per flight). The effects of design mean stress, the stress range for ground-air-ground cycles, simulated thermal stress, the number of stress cycles in each flight, and salt corrosion were studied. The flight-by-flight stress sequences were applied to notched sheet specimens of Ti-8Al-1Mo-1V and Ti-6Al-4V titanium alloys. A linear cumulative-damage analysis accounted for large changes in stress range of the simulated flights but did not account for the differences between real-time and accelerated tests. The fatigue lives from accelerated tests were generally within a factor of two of the lives from real-time tests; thus, within the scope of the investigation, accelerated testing seems feasible.

Imig, L. A.↗

Computational Methods For Composite Structures

Selected methods of computation for simulation of mechanical behavior of fiber/matrix composite materials described in report. For each method, report describes significance of behavior to be simulated, procedure for simulation, and representative results. Following applications discussed: effects of progressive degradation of interply layers on responses of composite structures, dynamic responses of notched and unnotched specimens, interlaminar fracture toughness, progressive fracture, thermal distortions of sandwich composite structure, and metal-matrix composite structures for use at high temperatures. Methods demonstrate effectiveness of computational simulation as applied to complex composite structures in general and aerospace-propulsion structural components in particular.

Chamis, Christos C.↗

Progressive failure analysis of fibrous composite materials and structures

A brief description is given of the modifications implemented in the PAFAC finite element program for the simulation of progressive failure in fibrous composite materials and structures. Details of the memory allocation, input data, and the new subroutines are given. Also, built-in failure criteria for homogeneous and fibrous composite materials are described.

Bahei-El-din, Yehia A.↗