Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HEXAGONAL CONFIGURATION”

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.

134 records · Page 8

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical and experimental characterization of twin transmission across grain boundaries along the forward and lateral directions

Pervasive deformation twinning and transmission events across grain boundaries (GBs) affect the strength and failure of hexagonal close-packed (HCP) magnesium. A three-dimensional twin can transmit along the twinning shear direction, (forward), and along the direction perpendicular to both the twinning plane normal and the shear direction, (lateral). For the first time, phase-field calculations and electron backscatter diffraction (EBSD)-based statistical analysis are combined to investigate the effect of the twinned grain boundary characteristics on twin transmission (TT) along the forward and lateral directions. This combined analysis reveals that TT propensity decreases with increasing misorientation angle across the GB for both forward and lateral directions. Also, the TT is more favorable along the lateral than along the forward direction. Twin transmission seems harder across GBs with a misorientation axis closer to the twin -direction than the other directions ( and ). Further, the PF calculations reveal that, at the onset of a transmission process, the crystallography tends to be preserved in the case of lateral transmission, whereas, in the forward transmission case, the transmitted twin punches straight through the GBs and its morphology prevails. Additionally, the EBSD analysis finds that pure forward and lateral transmissions do not occur often, yet reveals a preference for lateral propagation consistent with PF simulations. Further, the local twin transmission configurations observed in the actual material do not correspond to pure tilt or twist GBs, which are most commonly considered as model cases.

36 MATERIALS SCIENCE↗

Determination of Ice Cloud Models Using MODIS and MISR Data

Representation of ice clouds in radiative transfer simulations is subject to uncertainties associated with the shapes and sizes of ice crystals within cirrus clouds. In this study, we examined several ice cloud models consisting of smooth, roughened, homogeneous and inhomogeneous hexagonal ice crystals with various aspect ratios. The sensitivity of the bulk scattering properties and solar reflectances of cirrus clouds to specific ice cloud models is investigated using the improved geometric optics method (IGOM) and the discrete ordinates radiative transfer (DISORT) model. The ice crystal habit fractions in the ice cloud model may significantly affect the simulations of cloud reflectances. A new algorithm was developed to help determine an appropriate ice cloud model for application to the satellite-based retrieval of ice cloud properties. The ice cloud particle size retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) data, collocated with Multi-angle Imaging Spectroradiometer (MISR) observations, is used to infer the optical thicknesses of ice clouds for nine MISR viewing angles. The relative differences between view-dependent cloud optical thickness and the averaged value over the nine MISR viewing angles can vary from −0.5 to 0.5 and are used to evaluate the ice cloud models. In the case for 2 July 2009, the ice cloud model with mixed ice crystal habits is the best fit to the observations (the root mean square (RMS) error of cloud optical thickness reaches 0.365). This ice cloud model also produces consistent cloud property retrievals for the nine MISR viewing configurations within the measurement uncertainties.

MODIS↗

Interplay between charge ordering and geometric ferroelectricity in LuFe 2 O 4 /LuFeO 3 superlattices

Oxide superlattices have drawn great attentions owing to the intriguing coupling among elastic, electrical, and magnetic orderings at the interfaces and the emerging of improper ferroelectricity. Here, superlattices composed of hexagonal LuFeO 3 (h-LuFeO 3 ) and LuFe 2 O 4 are investigated via density functional theory calculations. h-LuFeO 3 is a well-known multiferroic material that is stable only in thin-film or doped bulk state, while LuFe 2 O 4 is a charge ordered material where the existence of ferroelectricity is still in controversy. We have found that the charge ordering (CO) induced polarizations in LuFe 2 O 4 layers coexist with the geometric polarizations in h-LuFeO 3 layers in the (LuFe 2 O 4 ) m /(LuFeO 3 ) n superlattices with different periodicities, and the ferroelectric states are generally preferred over the antiferroelectric states for LuFe 2 O 4 in superlattices. The out-of-plane polarizations in h-LuFeO 3 and LuFe 2 O 4 layers tend to be aligned in parallel, and the overall polarization increases with the ratio of h-LuFeO 3 . The influence of layered polarizations on the local electrostatic potential is not significant except the detected small trend caused by the CO-induced polarization within a FeO bilayer. Besides, the local electronic structures show that the Fermi level position in a certain layer can be tuned by the valences of Fe in this layer and the polarization distributions in neighboring layers. LuFe 2 O 4 layers sandwiched between thick h-LuFeO 3 layers are more susceptible. The calculated configurations of the superlattices are supported by atomic-resolution transmission electron microscopy experiments. Our results pave a way for tunable ferroelectricity in superlattice systems and create a new playground for manipulating the coupling between various degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗

A Novel High-Performance Mission-Enabling Multi-Purpose Radioisotope Heat Source

Recent studies indicate science mission concepts targeting access to the sub-surface oceans of icy moons require ice-penetrating cryobots powered by advanced Radioisotope Power Systems (RPS). These systems would deliver waste heat for ice-melting in the range of 10 kW. Minimizing the transit time through the kilometers-thick ice shells to just a few years requires these RPS to utilize heat sources having a higher thermal energy volumetric density than the existing flight-qualified General-Purpose Heat Source (GPHS). A Compact Heat Source (CPHS)has been conceptualized in which the graphite impact shells(GIS) of the existing GPHS are rearranged in a hexagonal aeroshell containing seven GIS per module, as opposed to the standard two per module; offering a thermal energy density of 0.57 W/cm3versus 0.29 W/cm3 offered by the GPHS simply from the repackaging of Technology Readiness Level (TRL)9 subassemblies. Preliminary thermal modeling of the CPHS integrated into a notional radioisotope thermoelectric generator(RTG) structure further suggests that centerline temperatures are well within allowable limits during nominal operation. Given the need for the CPHS for a subset of missions, it is worth exploring the applicability of the CPHS for more general RTG purposes. We discuss herein how the CPHS may be implemented with either heritage or in-development thermoelectric converter technologies into a Next-Generation RTG concept. Due to a higher energy density, the legacy heat rejection fin arrangement must be modified to permit a sufficiently low cold-side temperature. Preliminary finite element analysis suggests fin-root temperatures can be kept as low as 520 K while allowing the generator to fit within the usable dimensions of currently available United States Department of Energy shipping containers. Such temperatures would certainly be compatible with the use of high temperature thermoelectric converter technologies.. A prime candidate is the heritage silicon-germanium (SiGe) unicouple, whose design could be adapted by approximately halving the leg-length, but without changes in hot and cold junction interfaces, which are features critical to the proven performance and reliability of these devices. The estimated Beginning of Life power for a SiGe-based CPHS-RTG using 12 CPHS for a thermal inventory of 10.5 kW is greater than 600 W under deep space operating conditions. Using higher performance segmented couples currently in development that are based on skutterudite,La3−xTe4and 14-1-11 Zintl thermoelectric materials in lieu of the SiGe unicouples would increase the power level to more than1 kW. The high specific power (We/kg) attribute of CPHS-RTGs found in this study could potentially enable Radioisotope Electric Propulsion (REP) mission concepts. Past NASA REP mission concept studies identified specific power needs in excess of 6to 8 We/kg. Based on a GPHS-RTG-like system configuration, we show that at fin root temperatures between 530 K and 570K (deep space environment), specific powers exceeding 10 We/kg are achievable using high performance segmented thermoelectric converters. The compact sizing and power density of the CPHS-RTG would constitute a significant step upgrade in specific power when compared to heritage GPHS-RTG (approximately5.1 We/kg) and off-the-shelf Multi-Mission RTG (approximately 2.6 We/kg).

Nesmith, Bill J.↗

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗