Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Single crystal model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Novel Electric Nucleation Technique for Growing Large Single Crystal in Space

We present, herein, an electrical model for growing crystals without a seed which might not be free of defects and thereby still hinder the growth of a perfect crystal in space. The system is designed to confine nucleation to a single site automatically in an under saturated solution to avoid multiple nucleation. The technique is based on the effect of electrostriction, which is the tendency of a material to become more compressed in the presence of an electric field. The system is designed to create an electrical potential well between two hyperboloid electrodes with applied voltage at low frequency. The induced potential well between the electrodes oscillates at low frequency and attracts the solute and condenses it into the region of maximum field intensity. The alternating voltage prevents molecules with intrinsic charge from being attracted to the electrodes. The continuous presence of the electric field during the duration of the experiment, provides a continuous migration of the molecules toward the trapping site. This will eliminate the creation of a depletion region around the nucleation center and will enhance the crystal growth rate. Aside from the above mentioned advantages, the system is compact, safe to operate, and inexpensive to build.

Abdeldayem, Hossin↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗
DOE OSTI2022

Dendrite-resolved, full-melt-pool phase-field simulations to reveal non-steady-state effects and to test an approximate model

In this report we study the epitaxial, columnar growth of (multiply oriented) dendrites/cells for a spot melt in a polycrystalline Al–Cu substrate using two-dimensional, phase-field, direct numerical simulations (DNS) at the full-melt-pool scale. Our main objective is to compare the expensive DNS model to a much cheaper but approximate “line” model in which a single-crystal phase-field simulation is confined to a narrow rectangular geometry. To perform this comparison, we develop algorithms that automatically extract quantities of interest (QoIs) from both DNS and line models. These QoIs allow us to quantitatively assess the assumptions in the line model and help us analyze its discrepancy with the DNS model. We consider four sets of heat source parameters, mimicking welding and additive manufacturing conditions, that create a combination shallow and deep melt pools. Our largest DNS simulation used 16K x 14K grid points in space. Our main findings can be summarized as follows. Under AM conditions, the QoIs of line models are in excellent agreement with the full DNS results for both shallow and deep melt pools. Under welding conditions, the primary spacing of the DNS model is smaller than the prediction of line model. We identify a geometric crowding effect that accounts for the discrepancies between the DNS and line models. We propose two potential mechanisms that determine the response of the microstructure to geometric crowding.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synthesis and characterization of the novel breathing pyrochlore compound Ba 3 ⁢Tm 2⁢ Zn 5⁢ O 11

In this study, a novel material from the rare-earth based breathing pyrochlore family, Ba 3 ⁢Tm 2 ⁢Zn 5 ⁢O 11 , was successfully synthesized. Powder x-ray diffraction and high-resolution powder neutron diffraction confirmed phase purity and the $F\bar4$3⁢𝑚 breathing pyrochlore crystal structure, while thermogravimetric analysis revealed incongruent melting behavior compared to its counterpart, Ba 3 ⁢Yb 2⁢ Zn 5⁢ O 11 . High-quality single crystals of Ba 3 ⁢Tm 2 ⁢Zn 5 ⁢O 11 were grown using the traveling solvent floating zone technique and assessed using Laue x-ray diffraction and single crystal x-ray diffraction. Further, thermodynamic characterization indicated paramagnetic behavior down to 0.05 K, and inelastic neutron scattering measurements identified distinct dispersionless crystal electric field energy bands, with the fitted crystal electric field model predicting a single-ion singlet ground state and an energy gap of ∼9 meV separating it from the first excited (singlet) state. Additional low-energy excitation studies on single crystals revealed dispersionless bands at 0.8 and 1 meV. Computed phonon dispersions from first-principles calculations ruled out phonons as the origin of these modes, further illustrating the puzzling and unique properties of Ba 3 ⁢Tm⁢ 2 Zn 5 ⁢O 11 .

36 MATERIALS SCIENCE↗
NASA NTRS1988

Modeling and simulation of InP homojunction solar cells

Modeling and simulation of single-crystal InP homojunction solar cells has been performed using the PC-1D code. Cell design and performance have been optimized using the best available estimates of the various materials parameters. A comparison has been made of the predictions of the PC-1D model to those of other models. The optimum performance is predicted to give an efficiency approaching 21 percent at AM0. It is shown that in order to describe the performance of actual cells it is necessary to use larger values of the intrinsic carrier concentration and the surface recombination velocity (SRV) than have been reported in the literature. However, even with a near-maximum value of SRV (e.g., 10 to the 7th cm/s) it is necessary to reduce the minority carrier diffusion length in the emitter to only 0.01 micron in order to account for the relatively low quantum efficiency in the blue part of the spectrum. This indicates that improvement in the emitter bulk properties could be much more important than the SRV. Other loss mechanisms are also discussed; in particular, it is shown that recombination in the base, for good quality material, is relatively insignificant.

Yahia, A. H.↗
NASA NTRS1989

The evolution of complex type B Allende inclusion - An ion microprobe trace element study

Results are presented of a detailed trace-element and isotopic analyses of the constituent phases in each of the major textural parts (mantle, core, and islands) of a Type B refractory inclusion, the USNM 5241 inclusion from Allende, first described by El Goresy et al. (1985). The REE data on 5241 were found to be largely consistent with a model in which the mantle and the core of 5241 formed sequentially out of a single melt by fractional crystallization. The numerical models of REE evolution in the 5241 melt, especially that of Eu, require that a significant mass of spinel-free island material was assimilated into the evolving melt during the last half of the solidification history of 5241. The trace element results pbtained thus strongly support the interpretation of El Goresy et al. (1985) that the spinel-free islands in the 5241 are trapped xenoliths.

Macpherson, Glenn J.↗
NASA NTRS2020

A Flexible Parameterization for Shortwave and Longwave Optical Properties of Ice Crystals and Derived Bulk Optical Properties for Climate Models

We provide a parameterization of the extinction efficiency, single scattering albedo and asymmetry parameter of single ice crystals with any combination of particle volume, projected area, component aspect ratio and crystal distortion at any wavelength between 0.2 and 100 mm. The paramerization is an extension of the one previously published by van Diedenhoven et al. (2014, https://doi.org/10.1175/JAS-D-13-0205.1). In addition, the parameterized optical properties are integrated over size distributions yielding bulk extinction efficiencies, single scattering albedos and asymmetry parameters for large ranges of effective radii, particle component aspect ratios and crystal distortion values. The parameterization of single particle optical properties is evaluated with a reference database. The bulk optical properties are evaluated against the ice model selected for the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 products, for which accurate optical properties are available. Mean absolute errors in parameterized extinction efficiency, asymmetry parameter and single scattering albedo are shown to be 0.0272, 0.00890 and 0.00468, respectively for shortwave wavelengths, while they are 0.0641, 0.0368 and 0.0200 in the longwave. Shortwave and longwave asymmetry parameter and single scattering albedos are shown to vary strongly with particle component aspect ratio and distortion, resulting in substantial variation in shortwave fluxes, but relatively small variations in longwave cloud emissivity. The parameterization and bulk optical properties are made publicly available.

Bastiaan Van Diedenhoven↗
NASA NTRS2012

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

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.

Saether, Erik↗
DOE OSTI2024

Deciphering the morphology of transition metal carbonate cathode precursors

The performance and life of Li-ion battery cathode materials are determined by both the composition (crystal structure and transition metal ratio) and the morphology (particle size, size distribution, and surface area). Careful control of these two aspects is the key to long lasting, high-energy batteries that can undergo fast charge. Developing such cathodes requires manipulation of the synthesis conditions, namely the coprecipitation process to develop the precursor and a calcination step to lithiate and convert it to a transition metal oxide. In this paper, we utilize a combination of controlled synthesis, microscopic and spectroscopic characterization, and multi-scale mathematical modeling to shed light on the synthesis of cathode precursors. The complex interplay between the various chemical reactions in the co-precipitation process is studied to provide experimentalists with guidance on achieving composition control during synthesis. Further, the formation of a variety of morphologies of the primary particles and the driving force for agglomeration is mathematically described, for the first time, based on an energy minimization approach. Results suggest that the presence of Ni and/or Co significantly lowers the reaction rate constant compared to Mn, resulting in agglomerated growth in the former and single crystal growth in the latter. Modeling studies are used to provide a phase map describing the synthesis conditions needed to control the secondary particle size and corresponding size distribution. Finally, this paper represents an important step in developing a computationally guided approach to the synthesis of battery cathode materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Sensitivity of Cirrus Bidirectional Reflectance at MODIS Bands to Vertical Inhomogeneity of Ice Crystal Habits and Size Distribution

A common assumption in satellite imager-based cirrus retrieval algorithms is that the radiative properties of a cirrus cloud may be represented by those associated with a specific ice crystal shape (or habit) and a single particle size distribution. However, observations of cirrus clouds have shown that the shapes and sizes of ice crystals may vary substantially with height within the clouds. In this study we investigate the sensitivity of the top-of-atmosphere bidirectional reflectances at two MODIS bands centered at 0.65 micron and 2.11 micron to the cirrus models assumed to be either a single homogeneous layer or three distinct but contiguous, layers. First, we define the single- and three-layer cirrus cloud models with respect to ice crystal habit and size distribution on the basis of in situ replicator data acquired during the First ISCCP Regional Experiment (FIRE-II), held in Kansas during the fall of 1991. Subsequently, fundamental light scattering and radiative transfer theory is employed to determine the single scattering and the bulk radiative properties of the cirrus cloud. Regarding the radiative transfer computations, we present a discrete form of the adding/doubling principle by introducing a direct transmission function, which is computationally straightforward and efficient an improvement over previous methods. For the 0.65 micron band, at which absorption by ice is negligible, there is little difference between the bidirectional reflectances calculated for the one- and three-layer cirrus models, suggesting that the vertical inhomogeneity effect is relatively unimportant. At the 2.11 micron band, the bidirectional reflectances computed for both optically thin (tau = 1) and thick (tau = 10) cirrus clouds show significant differences between the results for the one- and three-layer models. The reflectances computed for the three-layer cirrus model are substantially larger than those computed for the single-layer cirrus. Finally, we find that cloud reflectance is very sensitive to the optical properties of the small crystals that predominate in the top layer of the three-layer cirrus model. It is critical to define the most realistic geometric shape for the small "quasi-spherical" ice crystals in the top layer for obtaining reliable single-scattering parameters and bulk radiative properties of cirrus.

Yang, P.↗
DOE OSTI2022

Pulling thin single crystal silicon wafers from a melt: The new leading-edge solar substrate

The Floating Silicon Method (FSM) has been established as a viable, stable method for growing single crystal ribbons directly from a silicon melt. With intense helium jet cooling to drive the linear progress of a [111] facet, pulled in the <110> direction, ribbons in the 0.6 – 3.0 mm thickness range can be grown at linear growth rates from 0.3 mm/s to >6 mm/s as reported in the literature. We report on recent progress towards growing (100) oriented ribbons with a net thickness of less than 200 microns and a ribbon width up to 18 cm using a stable, continuous process in the Leading Edge prototype furnace. Here, the 3D details of the single crystal growth are explained by extending the mechanics of the Limit Cycle Theory, with novel Internal Side Effect morphology described by a proposed Facet Flow Theory. Grown-in crystalline defect distributions are described as well as values of critical impurities like oxygen, carbon, dopants, and metals that are relevant for use as wafers for solar cells.

14 SOLAR ENERGY↗
NASA NTRS1992

The effects of small ice crystals on cirrus infrared radiative properties

Using a model that combines single-scattering properties for spheroidal and hexagonal ice crystals, the thermal infrared radiative properties of cirrus clouds have been investigated. Infrared scattering and absorption properties for randomly oriented spheroids and hexagons are parameterized based on the anomalous diffraction theory and a geometric ray-tracing method, respectively. Using observed ice crystal size distributions, upwelling radiances at the top of cirrus cloudy atmospheres have been computed. Results show that the presence of small ice crystals can produce significant brightness temperature differences between two infrared wavelengths in the 10-micron window. Theoretical results have been compared with observed brightness temperature differences between 8.35 and 11.16 microns and between 11.16 and 12 microns. The observed values were obtained from the high-spectral resolution interferometer sounder. It is shown that the use of the present nonspherical model for ice crystals in radiative transfer calculations leads to a significantly better interpretation of the observed data than does the use of the spherical model.

Takano, Y.↗
NASA NTRS1990

A crystallographic model for the tensile and fatigue response for Rene N4 at 982 C

An anisotropic constitutive model based on crystallographic slip theory was formulated for nickel-base single-crystal superalloys. The current equations include both drag stress and back stress state variables to model the local inelastic flow. Specially designed experiments have been conducted to evaluate the existence of back stress in single crystals. The results showed that the back stress effect of reverse inelastic flow on the unloading stress is orientation-dependent, and a back stress state variable in the inelastic flow equation is necessary for predicting inelastic behavior. Model correlations and predictions of experimental data are presented for the single crystal superalloy Rene N4 at 982 C.

Sheh, M. Y.↗
DOE OSTI2020

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗