Search NASASearch

SEARCH · Search NASA

Results for “Microstructural dynamics”

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 91 records · Page 5

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer

Influence of Microstructural Features on Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

Shape memory alloys (SMAs) exhibit several unique thermomechanical properties due to a reversible martensitic phase transformation. A key aspect of this transformation is the austenite-martensite interface. Here we present a molecular dynamics (MD) simulation methodology to study the atomic-scale features of austenite-martensite interface migration under near-equilibrium conditions. In single crystals, the interfaces migrate rapidly with only a small thermodynamic driving force. In polycrystals, however, interface migration is significantly impeded due to the change in orientation relationship at grain boundaries and the stored elastic energy resulting from microstructural constraints. This behavior can be linked to several mechanisms associated with transformation width and hysteresis in SMAs, properties of great importance for applications involving actuation. Additionally, we will present preliminary MD simulation results on the influence of precipitates on the formation and migration of austenite-martensite interfaces.

Gabriel Plummer

Special Grain Boundaries in NiTi Shape Memory Alloys as Sites for Preferential Martensite Nucleation

Shape memory alloys (SMAs) exhibit unique thermomechanical properties due to a reversible martensitic transformation which can be controlled via chemical composition and microstructural features. Among the latter, grain boundaries (GBs) are key in determining how the transformation nucleates and propagates. We performed molecular dynamics simulations to examine the roles of a few special GBs in the austenite phase of NiTi. The GBs can act as preferential host sites for martensite nuclei which can substantially lower the nucleation barrier for transformation and thereby result in reduced thermal hysteresis, an important SMA property for cyclic actuation applications. Free energy calculations show that a characteristic of these GBs is a negative entropy which drops sharply close to the transformation temperature, an anomalous behavior for GBs with fixed composition which is a direct result of the martensitic transformation. We discuss implications of these results with respect to SMA processing techniques for achieving improved properties.

Gabriel Plummer

Energy efficient engine, high pressure turbine thermal barrier coating. Support technology report

This report describes the work performed on a thermal barrier coating support technology task of the Energy Efficient Engine Component Development Program. A thermal barrier coating (TBC) system consisting of a Ni-Cr-Al-Y bond cost layer and ZrO2-Y2O3 ceramic layer was selected from eight candidate coating systems on the basis of laboratory tests. The selection was based on coating microstructure, crystallographic phase composition, tensile bond and bend test results, erosion and impact test results, furnace exposure, thermal cycle, and high velocity dynamic oxidation test results. Procedures were developed for applying the selected TBC to CF6-50, high pressure turbine blades and vanes. Coated HPT components were tested in three kinds of tests. Stage 1 blades were tested in a cascade cyclic test rig, Stage 2 blades were component high cycle fatigue tested to qualify thermal barrier coated blades for engine testing, and Stage 2 blades and Stage 1 and 2 vanes were run in factory engine tests. After completion of the 1000 cycle engine test, the TBC on the blades was in excellent condition over all of the platform and airfoil except at the leading edge above midspan on the suction side of the airfoil. The coating damage appeared to be caused by particle impingement; adjacent blades without TBC also showed evidence of particle impingement.

Duderstadt, E. C.

Temperature Dependent Cyclic Deformation Mechanisms in Haynes 188 Superalloy

The cyclic deformation behavior of a wrought cobalt-base superalloy, Haynes 188, has been investigated over a range of temperatures between 25 and 1000 C under isothermal and in-phase thermomechanical fatigue (TMF) conditions. Constant mechanical strain rates (epsilon-dot) of 10(exp -3)/s and 10(exp -4)/s were examined with a fully reversed strain range of 0.8%. Particular attention was given to the effects of dynamic strain aging (DSA) on the stress-strain response and low cycle fatigue life. A correlation between cyclic deformation behavior and microstructural substructure was made through detailed transmission electron microscopy. Although DSA was found to occur over a wide temperature range between approximately 300 and 750 C the microstructural characteristics and the deformation mechanisms responsible for DSA varied considerably and were dependent upon temperature. In general, the operation of DSA processes led to a maximum of the cyclic stress amplitude at 650 C and was accompanied by pronounced planar slip, relatively high dislocation density, and the generation of stacking faults. DSA was evidenced through a combination of phenomena, including serrated yielding, an inverse dependence of the maximum cyclic hardening with epsilon-dot, and an instantaneous inverse epsilon-dot sensitivity verified by specialized epsilon-dot -change tests. The TMF cyclic hardening behavior of the alloy appeared to be dictated by the substructural changes occuring at the maximum temperature in the TMF cycle.

Rao, K. Bhanu Sankara

Tensile behavior of the L(1)2 compound Al67Ti25Cr8

Temperature-related variations in tensile yield strength and ductility were studied on cast, homogenized and isothermally forged Al67Ti25Cr8. Yield strength dropped discontinuously between 623 K and 773 K and then decreased gradually with increasing temperature. Below 623 K, fracture occurred prior to macroscopic yielding. Ductility decreased from 0.2 percent at 623 K to zero at 773 K, but increased again at higher temperatures. At 1073 K, an elongation of 19 percent was obtainable, and ultimate tensile strength and localized necking were observed. Fracture surfaces and deformed microstructures were examined. The 1073 K tensile specimen that exhibited 19 percent elongation showed grain boundary serrations and some evidence of recrystallization (likely dynamic) although fracture occurred predominantly via an intergranular mode.

Kumar, K. S.

Toward an In-Depth Material Model for Cermet Nuclear Thermal Rocket Fuel Elements

The development and qualification of nuclear thermal propulsion (NTP) fuel element technologies would be aided by an in-depth model of material response and failure modes at operating conditions. Integrated computational materials engineering techniques have the potential to provide such a model, as demonstrated here through three case studies focused on a tungsten-uranium mononitride cermet fuel. The first case focuses on the erosion of tungsten (W, also named wolfram), a nominal coating/cladding material, in hot hydrogen. Ab initio techniques are used to calculate erosion rates and thermal expansion at NTP operating conditions. The second focuses on the stability of uranium mononitride (UN) fuels at high temperature and in the presence of hydrogen. Phase diagram techniques reveal potential instabilities and decomposition pathways at high hydrogen concentrations. The third focuses on using microstructure information to predict high temperature mechanical response and failure of tungsten, used in refractory cermet materials. Combined finite element and discrete dislocation dynamics techniques provide mechanical properties in agreement with experimental methods. The integration of these techniques for an all-encompassing material model is discussed.

cermet

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics

MRI Measurements and Granular Dynamics Simulation of Segregation of Granular Mixture

A counter intuitive axial segregation phenomenon in a rotating horizontal cylinder has recently captured attention of many researchers in different disciplines. There is a growing consensus that the interplay between the particle dynamics and the evolution of the internal structure during the segregation process must be carefully investigated. Magnetic resonance imaging (MRI) has been used to non-invasively obtain much needed dynamic/static information such as velocity and concentration profiles, and it has proven to be capable of depicting the evolution of segregation processes. Segregation in a rotating cylinder involves two processes: the first is to transport small particles in the radial direction to form a radial core, and the second is to transform the radial core into axially segregated bands. Percolation and/or "stopping" have been proposed as mechanisms for the radial segregation. As to mechanisms for axial band formation, much less is known. The difference in the dynamic angle of repose has been proposed to segregate different components in the axial direction. Recently, Hill and Kakalios have reported that particles mix or demix depending upon the competition between diffusion and preferential drift whose order can be determined by the dynamic angle of repose through the adjustment of the rotation rate. We claim that the dynamic angle of repose could be one of the causes, however, it fails to offer reasonable explanations for certain aspects of the axial migration. For example, we always observe that the radial segregation precedes the axial segregation and small particles migrate in the radial direction to form an axially extended radial core. It then transforms into axially segregated bands. By definition, the effects of the dynamic angle of repose are restricted near the free surface where the flowing layer is present. However, during the process of transforming from the radially segregated core to axially segregated bands, small particles located in the deep core region, which is untouched by the flowing layer, also completely disappear. Usually, the dynamics angle of repose are uniquely defined for individual species to characterize particle properties, and the dynamic angle of repose thus defined provides little information for the dynamic angle of repose of the mixture since the concentration ratio and the internal packing structure do not remain the same during the segregation processes. Under microgravity environment, the dynamics angle of repose argument does not hold since there is simply no flowing layer to influence/determine the preferred directions of segregation. We have thus designed an experiment so that the effects of the dynamic angle of repose can be minimized by filling the cylinder almost completely full. Small particles still formed a radial core and also migrated to form axial bands. As ground based experiments we have designed and conducted both 2D and 3D segregation experiments. The 2D experiments are performed using a thin cylinder (the gap between two end caps is about 5 mm) filled with different combinations of particles. The 3D experiments are conducted with a long cylinder of its length and diameter of 27cm and 7cm, respectively. Results of 2D experiments indicate that different mechanisms govern particle motion in regions near and far from the axis of rotation. Results of 3D experiments indicate that a series of collapses of microstructures of particle packing (micro-collapses) may be responsible for the creation of voids for small particles to migrate through in the axial direction. We have successfully eliminated the dynamic angle of repose as a cause for segregation, however, by almost completely filling the cylinder with the particles, we have lost an opportunity to investigate a possibility of particle "mobility" being a cause for segregation which requires a flowing surface but not the difference in the angle of repose. This is currently being investigated.

Nakagawa, M.

Fast Assessment of Metal Performance through Dislocation Physics and Machine Learning

The microstructure of metals is key to their mechanical properties. The types, density, composition and morphology of crystal defects all have pronounced impact on the properties. Changes to the microstructure occurring during processing and use can be very striking. The emerging technology additive manufacturing (AM) has the potential to improve performance by allowing optimized designs, but the process and environments can lead to unusual microscale features whose properties must be understood and characterized to enable higher technological readiness levels and application. Experimentally, an extensive evaluation of mechanical properties of 3D printed metals is a challenge, and anomalous effects related to the AM process add complexity. We present a new machine learning (ML) model predicting mechanical response based on dislocation mediated plasticity simulations. A large set of 3D discrete dislocation dynamics simulations with wide ranges of loading conditions is transformed to preprocessed data ready for training with the ML model. The trained model can predict the mechanical response of Mo30W for a given microstructure evolution, providing key information essential for optimization of AM processing.

Jaehyun Cho

Fast Assessment of Metal Performance through Dislocation Physics and Machine Learning

The microstructure of metals is key to their mechanical properties. The types, density, composition and morphology of crystal defects all have pronounced impact on the properties. Changes to the microstructure occurring during processing and use can be very striking. The emerging technology additive manufacturing (AM) has the potential to improve performance by allowing optimized designs, but the process and environments can lead to unusual microscale features whose properties must be understood and characterized to enable higher technological readiness levels and application. Experimentally, an extensive evaluation of mechanical properties of 3D printed metals is a challenge, and anomalous effects related to the AM process add complexity. We present a new machine learning (ML) model predicting mechanical response based on dislocation mediated plasticity simulations. A large set of 3D discrete dislocation dynamics simulations with wide ranges of loading conditions is transformed to preprocessed data ready for training with the ML model. The trained model can predict the mechanical response of Mo30W for a given microstructure evolution, providing key information essential for optimization of AM processing.

Jaehyun Cho

The Isothermal Dendritic Growth Experiment (IDGE)

Dendritic solidification is one of the simplest examples of pattern formation where a structureless melt evolves into a ramified crystalline microstructure; it is a common mode of solidification in many materials, but especially so in metals and alloys. There is considerable engineering interest in dendrites because of the role dendrites play in the determination of microstructure, and thereby in influencing the physical properties of cast metals and alloys. Dendritic solidification provides important examples of non-equilibrium physics, pattern formation dynamics, and models for computational condensed matter and material physics. Current theories of dendritic growth generally couple diffusion effects in the melt with the physics introduced by the interface. Unfortunately, in terrestrial based experiments, convective effects in the melt alter the growth process in such a manner as to prevent definitive analysis of convective, diffusive or interfacial effects. Thus, the effective elimination of convection in the melt by operating experiments on orbit were required to produce high-fidelity data needed for achieving further progress. This simple fact comprised the scientific justification for the IDGE.

Glicksman, M. E.

Microstructure Modeling of 3rd Generation Disk Alloy

The objective of this initiative, funded by NASA's Aviation Safety Program, is to model, validate, and predict, with high fidelity, the microstructural evolution of third-generation high-refractory Ni-based disc superalloys during heat treating and service conditions. This initiative is a natural extension of the DARPA-AIM (Accelerated Insertion of Materials) initiative with GE/Pratt-Whitney and with other process simulation tools. Strong collaboration with the NASA Glenn Research Center (GRC) is a key component of this initiative and the focus of this program is on industrially relevant disk alloys and heat treatment processes identified by GRC. Employing QuesTek s Computational Materials Dynamics technology and PrecipiCalc precipitation simulator, physics-based models are being used to achieve high predictive accuracy and precision. Combining these models with experimental data and probabilistic analysis, "virtual alloy design" can be performed. The predicted microstructures can be optimized to promote desirable features and concurrently eliminate nondesirable phases that can limit the reliability and durability of the alloys. The well-calibrated and well-integrated software tools that are being applied under the proposed program will help gas turbine disk alloy manufacturers, processing facilities, and NASA, to efficiently and effectively improve the performance of current and future disk materials.

Jou, Herng-Jeng

Measurement of the time-temperature dependent dynamic mechanical properties of boron/aluminum composites

A relatively simple flexural vibration test is developed for accurate measurement of the low-strain dynamic modulus and damping capacity of B/Al composite bar specimens from -200 C to over 500 C. The specimens are prepared from 8-ply unidirectional panels containing 50 volume percent fibers composed of 203-micron commercial boron-on-tungsten fibers. The basic test technique consists of the forced flexural vibration of the composite bar specimens at their two lowest free-free symmetrical resonant modes in a high-vacuum cryostat furnace. Specimen damping is determined from oscilloscope photographs of the free decay obtained after simultaneously removing the resonant drive signal and grounding the vibration-drive electrode. The availability of time-temperature dynamic data coupled with the predictive accuracy of composite theory suggests a future potential for using such data in examining environmental effects on composite macrostructure and microstructure.

Dicarlo, J. A.

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan

Estimation of high temperature low cycle fatigue on the basis of inelastic strain and strainrate

Fatigue life at elevated temperature can be predicted by introducing parametric values obtained from monotonic constitutive behavior into the Universal-Slopes Equation. For directionally solidified MAR-M200+HF at 975 C, these parameters are the maximum stress achievable under entirely plastic (time-independent) and purely creep (time-dependent) conditions and the corresponding inelastic strains, as well as the elastic modulus. For materials which exhibit plasticity/creep interaction, two more pairs of monotonic parameters must be evaluated for fatigue life prediction. This life-prediction method based on the Universal-Slopes Equation, resulted from a constitutive model characterizing monotonic and cyclic data as inelastic strainrate as a function of inelastic strain. Characterizing monotonic data is this way, permitted distinction between different material responses such as strain-hardening, strain-softening, and dynamic recovery effects. Understanding and defining the region of influence of each of these effects facilitated formulation of the constitutive model in relation to the mechanical and microstructural processes occurring in the material under cyclic loading.

Berkovits, A.

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics