SEARCH · Engineering Papers
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.
Vapor phase growth method for single crystalline preparation of gallium nitride, gallium arsenide alloy, and gallium antimonide alloy
Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.
Calcium, aluminum-rich inclusions (CAI's) from carbonaceous chondrite are believed to be among the first solids to have formed in the early solar system and thus to preserve a record of conditions and processes that prevailed early in the solar nebula. We have measured inter- and intra- crystalline oxygen isotope distribution of fassaite in Type B1 CAI of Allende meteorite. Secondary ion mass spectrometry (SIMS) using negative Au primary ion was applied to in situ oxygen isotope ratio analysis on a polished thin section of the CAI. Present data indicate that O isotope composition is homogeneous inside of fassaite grains as well as among grains in the HN3-1 CAI, and that the delta(exp 18) O value (-24.14 percent relative to SMOW) determined previously could be applied to most fassaite grains. The results suggest that most fassaite grains are preexisting solids like spinel grains.
Perovskite photovoltaics reached record efficiencies in the laboratory and, if sustainably commercialized, they would accelerate a green energy transition. Here, this article presents the development of life-cycle inventory material and energy databases of most promising four single-junction and three tandem scalable perovskite systems with assumptions regarding scalable production validated by industry experts. We conducted comprehensive “ex ante” Life Cycle Analysis (LCA) and Net Energy Analysis, analyzing their cumulative energy demand, global warming potential profiles, energy payback times and energy return on investment (EROI). LCA contribution analysis elucidates the most impactful material and process choices. It shows that solution-based perovskite manufacturing would have lower environmental impact than vapor-based methods, and that roll-to-roll (RtR) printing offers the lowest impact. Among material choices, MoOx/Al has lower impact than Ag, and fluorine-tin-oxide lower than indium-tin-oxide. Furthermore, we compare perovskites with commercial crystalline-silicon and thin-film PV, accounting for the most recent developments in crystalline-Si wafer production and differences in life expectancies and efficiencies. It is shown that perovskite systems produced with RtR manufacturing could reach in only 12 years of life the same EROI as that of single-crystalline-Si PV lasting 30-years. This work lays a foundation for sustainability investigations of perovskite large-scale deployment.
Groundwater contaminated with arsenic (As) must be treated prior to drinking, as human exposure to As at toxic levels can cause various diseases including cancer. Conventional aeration-filtration applied to anaerobic arsenite (As(III)) contaminated groundwater can remove As(III) by co-oxidizing native iron (Fe(II)) and As(III) with oxygen (O 2 ). However, the As(III) removal efficiency of conventional aeration can be low, in part, because of incomplete As(III) oxidation to readily-sorbed arsenate (As(V)). In this work, we investigated a new approach to enhance As(III) co-removal with native Fe(II) by the anaerobic addition of hydrogen peroxide (H 2 O 2 ) prior to aeration. Experiments were performed to co-oxidize Fe(II) and As(III) with H 2 O 2 (anaerobically), O 2 (aerobically), and by sequentially adding of H 2 O 2 and O 2 . Aqueous As(III) and As(V) measurements after the reaction were coupled with solid-phase speciation by Fe and As K-edge X-ray absorption spectroscopy (XAS). We found that complete anaerobic oxidation of 100 µM Fe(II) with 100 µM H 2 O 2 resulted in co-removal of 95% of 7 µM As(III) compared to 44% with 8.0-9.0 mg/L dissolved O 2 . Furthermore, we found that with 100 µM Fe(II), the initial Fe(II):H 2 O 2 ratio was a critical parameter to remove 7 µM As(III) to below the 10 µg/L (0.13 µM) WHO guideline, where ratios of 1:4 (mol:mol) Fe(II):H 2 O 2 led to As(III) removal matching that of 7 µM As(V). The improved As(III) removal with H 2 O 2 was found to occur partly because of the well-established enhanced efficiency of As(III) oxidation in Fe(II)+H 2 O 2 systems relatively to Fe(II)+O 2 systems. However, the XAS results unambiguously demonstrated that a large factor in the improved As(III) removal was also due to a systematic decrease in crystallinity, and thus increase in specific surface area, of the generated Fe(III) (oxyhydr)oxides from lepidocrocite in the Fe(II)+O 2 system to poorly-ordered Fe(III) precipitates in the Fe(II)+H 2 O 2 system. The combined roles of H 2 O 2 (enhanced As(III) oxidation and structural modification) can be easily overlooked when only aqueous species are measured, but this dual impact must be considered for accurate predictions of As removal in groundwater treatment.
Amorphous grain boundary complexions lack long-range crystalline order but are not featureless, as distinct gradients in structural short-range order have been reported through their thickness. In this work, we test the hypothesis that the distribution of short-range order is determined by the confining crystals using atomistic simulations of both Cu-Zr bicrystals and a random polycrystal. Voronoi polyhedra with structures similar to that of perfect face-centered cubic serve as signatures of high structural order and are only found at the amorphous-crystalline transition regions. The density of the ordered structural motifs within a specific amorphous-crystalline transition region is found to not be directly determined by the orientation and symmetry of the grain which touches it, but rather by the incompatibility between the two confining grains. Further, ordered polyhedra density is found to be inversely related to grain incompatibility, meaning that large incompatibilities between the confining crystals lead to less order in the amorphous-crystalline transition region. The finding that the entire grain-film-grain system must be considered to understand local structure unequivocally demonstrates that amorphous complexions are not simply a collection of independent phases which happen to nucleate at a grain boundary. Rather, an amorphous grain boundary complexion is a single entity that finds a local equilibrium configuration.
The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.
Supramolecular polymers can emulate some of the physical properties of covalent polymers but offer new opportunities given the possibility of designing monomers that will form highly ordered assemblies with defined shapes. Internally ordered supramolecular polymers formed through nucleation–elongation self-assembly are well-known but highly crystalline examples which exhibit important properties such as light harvesting, charge transport, and ferroelectricity are not common. We report here on a detailed study of supramolecular polymers formed in water by carboxylated naphtho-p-quinodimethane amphiphiles. We found that supramolecular polymerization of these amphiphiles in aqueous media yields crystalline assemblies with morphologies that included ribbons, helically rolled ribbons, and twisted filaments. This polymorphism was found to be controlled exclusively by repulsive electrostatic interactions controlled by the degree of protonation of the carboxylic head groups which also dictates the nature of supramolecular packing. Substoichiometric amounts of base lead to highly crystalline ribbons due to a decreased surface charge density and less electrostatic repulsion. Increasing deprotonation results in helically rolled ribbons with a different polymorph crystal lattice, whereas excessive deprotonation leads to twisted filaments with maximum surface charge density. Ribbons, helical rolled ribbons, and twisted filaments revealed an increasing red shift in their visible absorption maxima. These crystalline assemblies could be potential candidates for solar energy materials and photocatalytic systems.
Powder X-ray diffraction (XRD) is a foundational technique for characterizing crystalline materials. However, the reliable interpretation of XRD patterns, particularly in multiphase systems, remains a manual and expertise-demanding task. As a characterization method that only provides structural information, multiple reference phases can often be fit to a single pattern, leading to potential misinterpretation when alternative solutions are overlooked. To ease humans’ efforts and address the challenge, we introduce Dara (data-driven automated Rietveld analysis), a framework designed to automate the robust identification and refinement of multiple phases from powder XRD data. Dara performs an exhaustive tree search over all plausible phase combinations within a given chemical space and validates each hypothesis using the BGMN Rietveld refinement routine. Key features include structural database filtering, automatic clustering of isostructural phases during tree expansion, and peak-matching-based scoring to identify promising phases for refinement. When ambiguity exists, Dara generates multiple hypothesis which can then be decided between by human experts or with further characterization tools. By enhancing the reliability and accuracy of phase identification, Dara enables scalable analysis of realistic complex XRD patterns and provides a foundation for integration into multimodal characterization workflows, moving toward fully self-driving materials discovery.
Obtaining high-quality 3D reconstructions from electron tomography of crystalline particles embedded in lighter support elements is crucial for various material systems such as catalysts for fuel cell applications. However, significant challenges arise due to the limited tilt range, sparse and low signal-to-noise ratio of the measurements. In addition, small metal particles can cause strong streaking and shading artifacts in the 3D reconstructions when using conventional reconstruction algorithms due to the presence of Bragg diffraction and the large scattering cross-section difference between the materials of the particles and the background support regions. These artifacts lead to errors in the downstream characterization affecting extraction of critical features such as the size of the metal particles, their distribution and the volume of the lighter support regions. In this paper, we present a two-stage algorithm based on metal artifact reduction, utilizing model-based iterative reconstruction methods with adaptive adjustment of regularization parameters. Our approach yields high-quality 3D reconstructions compared to traditional algorithms, accurately capturing both the metal particles as well as the background support. We demonstrate the effectiveness of our algorithm through simulated and experimental bright-field electron tomography data, showing significant improvements in reconstruction quality compared to traditional methods.
Yielding behavior and dislocation motion in two phase crystalline solid
Composites were fabricated from unidirectional unsized AS-4 carbon fiber and two baseline polyimides: benzophenone dianhydride-3,3'-diaminodiphenylsulfone (PISO2) and benzophenone dianhydride-3,3'-diaminobenzophenone (LARC-TPI). In addition, each polymer solution prior to prepregging was doped with various amounts of crystalline LARC-TPI powder to enhance melt flow during press molding. The 2:1, 1:1, and 1:2 ratios of crystalline to amorphous resin, respectively, were studied in the LARC-TPI system and the 1:2 ratio in the PISO2 system. Matrix characterization, prepreg fabrication/characterization and composite fabrication and physical/mechanical properties are described. The latter include three point short beam shear and flexure, dry and wet, as well as fracture toughness properties in selected compositions.
The softening of glassy grain boundaries in ceramic matrix composites and Si3N4 at high temperatures reduces mechanical strength and the upper-use temperature. By crystallizing this glass to a more refractory crystalline phase, a material which performs at higher temperatures may result. Three systems were examined: a cordierite composition with ZrO2 as a nucleating agent; celsian compositions; and yttrium silicate glasses both in bulk and intergranular in Si3N4. For the cordierite compositions, a series of metastable phases was obtained. The crystallization of these compositions was summarized in terms of metastable ternary isothermal sections. Zircon formed at the expense of ZrO2 and spinel. In SiC composites, the transformations were slower. In celsian, two polymorphs were crystallized. One phase, hexacelsian, which always crystallized, even when metastable, had an undesirable volume change. The other phase, celsian, was very difficult to crystallize. In yttrium silicate bulk glasses, similar in composition to the intergranular glass in Si3N4, a number of polymorphs of Y2Si2O7 were crystallized. The conditions under which these polymorphs formed are compared with crystallization in Si3N4.
We examine the infrared emission of the Herbig Ae/Be stars and show that some possess characteristics indicative of partially crystalline grains similar to those seen in Beta Pictoris and some solar system comets.
Impact cratering is one of the most ubiquitous geologic processes shaping the surface of all solid bodies in our solar system. Impacts are also a major source of clay minerals and poorly crystalline, clay-like materials on Earth and Mars[1,2].These phyllosilicates and related clay-like phases comprise an incredibly complex group of materials, and their characterization, even in controlled laboratory settings, remains a challenging endeavor. The nature and origin of clay minerals and amorphous materials on Mars, which form a major component (~20-70 wt %) of rock and soil samples in Gale Crater as determined by the CheMin instrument on Curiosity, have remained ambiguous[3]. These amorphous phases likely fall on a spectrum between pristine volcanic and/or impact-produced primary materials and clay minerals and related phases formed from aqueous alteration. Amorphous materials are common weathering products in terrestrial sediments, soils, and paleosols [e.g., 4,5].Primary impact materials (glass, melt rocks)are comparable in some ways to those generated volcanically [6], and so it may be possible that an amorphous component is preserved within altered impactites. There are hundreds of thousands of impact craters on Mars, and Curiosity and Perseverance are currently exploring ancient impact craters. We hypothesize that the sediments and lithologies in Gale and Jezero Craters and elsewhere contain altered impact products; it is therefore important to better understand the composition and structure of these materials. Here we present the first results from a study characterizing materials produced from impact glass alteration–clay minerals, poorly crystalline/amorphous materials–in terrestrial craters using Transmission Electron Microscopy (TEM), powder X-ray diffraction (pXRD)and chemical (EDS) analysis.
The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.
The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.
This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.