SEARCH · Search NASA
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.
MCFET computer program, named for combined Monte Carlo Finite Element Technique, implements specialized microstructural-lattice simulation technique. Developed to simulate microstructural evolution in material systems in which modulated phases occur and directionality of modulation influenced by internal and external stresses. Written in FORTRAN.
Thermal and compositional gradients induced during the growth process contribute significantly to the development of defects in the solidified boule. Thermal gradients and the solid-liquid interface shape can be greatly effected by ampoule material. Compositional gradients are strongly influenced by interface curvature and convective flow in the liquid. Results of this investigation illustrate the combined influences of interface shape and convective fluid flow. An applied magnetic field was used to reduce the effects of convective fluid flow in the electrically conductive melt during directional solidification. Several 8 mm diameter boules of Ga-doped Ge were grown at different field strengths, up to 5 Tesla, in four different ampoule materials. Compositional profiles indicate mass transfer conditions ranged from completely mixed to diffusion controlled. The influence of convection in the melt on the developing crystal microstructure and defect density was investigated as a function of field strength and ampoule material. Chemical etching and electron backscattered electron diffraction were used to map the crystal structure of each boule along the center plane. Dislocation etch pit densities were measured for each boule. Results show the influence of magnetic field strength and ampoule material on overall crystal quality.
Austenitic stainless steels, such as types 316 and 304, are widely used as pressure vessel materials in the temperature range of 425 to 650 C. Stainless steel specimens were tested to rupture at two different stress levels sigma and sigma 2 sigma 1 sigma 2) to establish the normal stain-time behavior. A subsequent test was performed in which the specimen was crept at the higher stress (sigma 1) to the beginning of the secondary stage of creep, presumed to be the strain/time conditions at which a steady state microstructure is developed, and then the stress was reduced to the lower level (sigma 2). The associated microstructure, and significance of this microstructure on the creep strain-hardening model for variable uniaxial loads were assesed and found to be consistent with the use of creep-recovery models at high stresses and temperatures and strain-hardening models at low stresses and tempertures.
This project made advancements in developing lower cost steel alloys with novel microstructural design for use in hydrogen refueling infrastructure such as storage, compressors, and dispensing components, and utilizing accelerated test methods to efficiently evaluate variations in alloy and microstructure design. The project specifically sought to design alloys with lower nickel contents to reduce alloy cost, which was accomplished through substituting manganese for nickel along with other alloy additions to control deformation characteristics known to be important for hydrogen embrittlement resistance. Through Mn substitutions for nickel, austenitic and duplex austenite-ferrite steels were successfully developed with lower cost than currently available commercial stainless steel products that are employed for hydrogen refueling infrastructure. The steels were processed to achieve comparable strength and toughness in hydrogen environments as the commercially available steels containing high Ni contents. To evaluate mechanical performance in hydrogen, a testing methodology was employed to compare ubiquitous laboratory testing using electrochemical hydrogen charging in a liquid electrolyte to less accessible high pressure gaseous testing. While the application of these steels is in high pressure gaseous environments, the electrochemical hydrogen charging tests produced comparable results. Additionally, the Los Alamos Neutron Scatting Center enabled characterization of deformation mechanisms of the steel alloys in the presence of hydrogen, which has been associated with steel alloy characteristics associated with hydrogen embrittlement. Finally, a fracture mechanics based test bed model was developed to predict the influence of hydrogen gas pressure and fatigue conditions on fatigue lifetimes of pressure vessel steels. Together, these developments can be employed to enable lower cost hydrogen fueling infrastructure and more reliable prediction of steel alloy components in hydrogen service conditions. In particular, the newly alloys have the potential to replace stainless steels, having demonstrated comparable performance at substantially reduced cost.
In order to identify the possible hierarchy of microstructural pathways and transitions, a systematic evaluation of the microstructural evolution in undercooled Fe-Ni alloys was performed on uniformly processed samples covering seven orders of magnitude in volume. At appropriate undercooling levels, alternate solidification pathways become thermodynamically possible and metastable product structures can result from the operation of competitive solidification kinetics. For thermal history evaluation, a heat flow analysis was applied and tested with large Fe-Ni alloy particles (1 to 3 mm) to assess undercooling potential. Alloy powders (10 to 150 microns), with large liquid undercoolings, were studied under the same composition and processing conditions to evaluate the solidification kinetics and microstructural evolution, including fcc/bcc phase selection and the thermal stability of a retained metastable bcc phase. The identification of microstructural transitions with controlled variations in sample size and composition during containerless solidification processing was used to develop a microstructure map which delineates regimes of structural evolutions and provides a unified analysis of experimental observations in the Fe-Ni system.
A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.
Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use the three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15 percent and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50 to 57 percent of the cell faces were pentagonal while 24 to 28 percent were quadrilateral and 15 to 22 percent were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with 3 quadrilateral, 6 pentagonal faces and 2 hexagonal faces consistent with the 3-6-2 cell.
Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use an idealized three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15 percent and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50 to 57 percent of the cell faces were pentagonal while 24 to 28 percent were quadrilateral and 15 to 22 percent were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with three quadrilateral, six pentagonal faces and two hexagonal faces consistent with the 3-6-2 Matzke cell. A compilation of 90 years of experimental data reveals that the average number of cell faces decreases linearly with the increasing ratio of quadrilateral to pentagonal faces. It is concluded that the Kelvin model is not supported by these experimental data.
Establishing the geometry of foam cells is useful in developing microstructure-based acoustic and structural models. Since experimental data on the geometry of the foam cells are limited, most modeling efforts use an idealized three-dimensional, space-filling Kelvin tetrakaidecahedron. The validity of this assumption is investigated in the present paper. Several FeCrAlY foams with relative densities varying between 3 and 15% and cells per mm (c.p.mm.) varying between 0.2 and 3.9 c.p.mm. were microstructurally evaluated. The number of edges per face for each foam specimen was counted by approximating the cell faces by regular polygons, where the number of cell faces measured varied between 207 and 745. The present observations revealed that 50-57% of the cell faces were pentagonal while 24-28% were quadrilateral and 15-22% were hexagonal. The present measurements are shown to be in excellent agreement with literature data. It is demonstrated that the Kelvin model, as well as other proposed theoretical models, cannot accurately describe the FeCrAlY foam cell structure. Instead, it is suggested that the ideal foam cell geometry consists of 11 faces with 3 quadrilateral, 6 pentagonal faces and 2 hexagonal faces consistent with the 3-6-2 Matzke cell
Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.
Understanding and controlling the development of deformation twins is paramount for engineering strong and stable hexagonal close-packed (HCP) Mg alloys. Actual twins are often irregular in boundary morphology and twin crystallography, deviating from the classical picture commonly used in theory and simulation. In this work, the elastic strains and stresses around irregular twins are examined both experimentally and computationally to gain insight into how twins develop and the microstructural features that influence their development. A nanoprecession electron diffraction (N-PED) technique is used to measure the elastic strains within and around a $\left\{ {10\overline{1}2} \right\}$ tensile twin in AZ31B Mg alloy with nm scale resolution. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) crystal plasticity model of the same sub-grain and irregular twin structure is employed to understand and interpret the measured elastic strain fields. The calculations predict spatially resolved elastic strain fields in good agreement with the measurement, as well as all the stress components and the dislocation density fields generated by the twin, which are not easily obtainable from the experiment. The model calculations find that neighboring twins, several twin thicknesses apart, have little influence on the twin-tip micromechanical fields. Furthermore, this work reveals that irregularity in the twin-tip shape has a negligible effect on the development of the elastic strains around and inside the twin. Importantly, the major contributor to these micromechanical fields is the alignment of the twinning shear direction with the twin boundary.
Measurements of convective flow fields and solid-liquid interface shapes during the solidification of a pure and a slightly alloyed transparent material reveal that the convective transport of solute can cause a macroscopic depression to develop in the solid-liquid interface. This effect occurs under conditions close to those which are predicted to produce morphological instability of a planar interface. A cellular or dendritic microstructure later develops within the interface depression. The convection is attributed to the effect of radial temperature gradients in the crystal growth apparatus.
This report describes the development of a microstructural model that can quantify the uncertainty in the observed rupture life of Grade 91 steel. The model is microstructural, meaning it relates microstructural characteristics of the material to the resulting material response. As such, one of the uses of this model is to identify the key microstructural parameters controlling the development of damage in Grade 91 operating at elevated temperatures. The report describes two veins of work: improvements to the crystal plasticity model required to run the uncertainty quantification analysis and the results of that UQ analysis. For creep, the model identifies the grain boundary diffusivity as the critical parameter controlling the rupture life of the material. The report demonstrates that a reasonable microstructural distribution of grain boundary diffusivity can account for the observed macroscale variation in rupture life at fixed temperature and load.
In third generation advanced high strength steels (AHSS) there is considerable potential to alter the transformation characteristics through variation of process parameters and through changes in chemical composition. This proposal aims to utilize this known phenomenon to develop microstructure-based transformation models to predict austenite stability as a function of chemical composition, strain level, deformation mode. The ultimate intent is to predict microstructure and performance during room temperature forming to produce automotive structural components with superior strength and failure resistance.
The poster summarizes the main results for a LDRD focused on developing a arc-based additive manufacturing (AM) system capable of fabricating large components with a composition gradient otherwise known as functionally graded materials (FGMs). FGMs are capable of addressing local extreme conditions in an optimal and cost-effective manner. The AM system (referred to as multi-Wire Arc Additive Manufacturing or m-WAAM) consisted of modifying a standard automated gas tungsten arc welding system to deliver multiple welding wires to the weld pool at controlled rates to deposit desired alloy compositions at critical locations. Numerous shapes and composition gradients were fabricated with the system and the effect of various deposition parameters on the nature of the composition profile were evaluated. Additionally, heat treatment of the FGM was performed to develop graded microstructures which should enhance dynamic impact performance of inexpensive steels. The project highlights the potential of m-WAAM to efficiently fabricate FGMs. However, this project only explored the science behind the technique and additional development work will be required to develop a fully-functional system that can build more complex components.
A model is presented which predicts the stress-strain behavior of continuous fiber reinforced laminated composites in the presence of microstructural damage. The model is based on the concept of continuum damage mechanics and uses internal state variables to characterize the various damage modes. The associated internal state variable growth laws are mathematical models of the loading history induced development of microstructural damage. The model is demonstrated by using it to predict the response of damaged AS-4/3502 graphite/epoxy laminate panels.
Microstructure formation of a hypo-eutectic Al-11Cu (atom percent) alloy during solidification after laser melting has been studied by transmission electron microscopy (TEM). The evolution of the solid-liquid-interface velocity, V SL , during the solidification process has been determined from direct observation by in-situ TEM. This enabled correlating V SL with development of four distinct microstructure zones. Crystal growth mode transitions from planar to cellular, cellular to dendritic, dendritic to cellular, and cellular to planar, have been observed for the accelerating solid-liquid-interface. The transition from coupled two-phase growth to the single-phase growth occurred for V SL = V a = (0.80±0.05) m/s, where V a is the velocity of absolute stability, at the onset of banded morphology grain formation. The in-situ and post-mortem TEM uniquely permitted determination of the non-equilibrium solidus for the rapidly solidifying Al-11Cu alloy. Experimental evidence for solute clustering and chemical ordering tendencies at length scales on the order of ≤ 5nm has been detected in the single-phase regions of the banded grains. The structural features of the single-phase bands have been interpreted as signatures of ‘frozen in’ configurations present in the liquid boundary layer adjacent to the growing crystal, which has a width of about 3nm. Finally, the nano-scale spatiotemporal resolution experimental TEM studies performed here provided quantitative metrics, e.g. the solidification interface velocity dependent solute concentration of the α-Al phase and the near-atomic scale structure in the single-phase bands, that are uniquely suitable for comparison with theory and model predictions for solidification microstructure development in multicomponent alloys after laser melting.
Microstructure-aware models are necessary to predict the behavior of material based on process knowledge or to extrapolate mechanical properties of materials to environmental conditions which are not easily reproduced in the laboratory, e.g., nuclear reactor environments. Elemental Ta provides a relatively simple BCC system in which to develop a microstructural understanding of deformation processes which can then be applied to more complicated BCC alloys. In situ neutron diffraction during compressive deformation and subsequent heat treatment have been used to monitor the evolution of microstructural features in Ta throughout simulated processing steps. Crystallographic texture and dislocation density are determined as a function of first plastic strain, then temperature. Lattice strains are determined and attributed to stresses at macroscopic, grain and dislocation length scales. The increase of the dislocation density through deformation and subsequent recovery during heat treatment is monitored through the changing diffraction line profile. Also, randomization of the texture is used as a signature of recrystallization. The recovery of dislocations through annihilation is not observed to depend on the initial dislocation density in the range studied here. In contrast, recrystallization is observed to depend strongly on the initially dislocation density.