Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “configurational disorder”

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.

197 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An FDA-approved drug structurally and phenotypically corrects the K210del mutation in genetic cardiomyopathy models

Dilated cardiomyopathy (DCM) due to genetic disorders results in decreased myocardial contractility, leading to high morbidity and mortality rates. There are several therapeutic challenges in treating DCM, including poor understanding of the underlying mechanism of impaired myocardial contractility and the difficulty of developing targeted therapies to reverse mutation-specific pathologies. In this report, we focused on K210del, a DCM-causing mutation, due to 3-nucleotide deletion of sarcomeric troponin T (TnnT), resulting in loss of Lysine210. We resolved the crystal structure of the troponin complex carrying the K210del mutation. K210del induced an allosteric shift in the troponin complex resulting in distortion of activation Ca 2+ -binding domain of troponin C (TnnC) at S69, resulting in calcium discoordination. Next, we adopted a structure-based drug repurposing approach to identify bisphosphonate risedronate as a potential structural corrector for the mutant troponin complex. Cocrystallization of risedronate with the mutant troponin complex restored the normal configuration of S69 and calcium coordination. Risedronate normalized force generation in K210del patient-induced pluripotent stem cell–derived (iPSC-derived) cardiomyocytes and improved calcium sensitivity in skinned papillary muscles isolated from K210del mice. Systemic administration of risedronate to K210del mice normalized left ventricular ejection fraction. Collectively, these results identify the structural basis for decreased calcium sensitivity in K210del and highlight structural and phenotypic correction as a potential therapeutic strategy in genetic cardiomyopathies.

Research & Experimental Medicine↗
DOE OSTI2021

Thermodynamics of boron incorporation in BGaN

We study the thermodynamics of boron (B) incorporation into gallium nitride (GaN) using first-principles calculations. In the dilute limit, we have calculated the formation energies of different configurations of the B impurity in GaN and found that substitution on the cation site is favored over substitution on the anion site. Furthermore, under p-type conditions, interstitial boron can become the more favorable configuration and will ultimately limit the p-type conductivity. At higher B concentrations we use the generalized quasi-chemical approximation to elucidate the thermodynamic stability of boron gallium nitride (BGaN) alloys. We also investigate the effects of strain, which will be present if BGaN alloys are grown pseudomorphically on a GaN substrate.Without strain, B incorporation at typical growth conditions is limited to about 1.4% at 800 °C. Pseudomorphic strain raises the limit to 3.0% at the same temperature, close to experimentally observed levels of B incorporation.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
DOE OSTI2023

Range and strength of mechanical interactions of force dipoles in elastic fiber networks

Mechanical forces generated by myosin II molecular motors drive diverse cellular processes, most notably shape change, division and locomotion. These forces may be transmitted over long range through the cytoskeletal medium – a disordered, viscoelastic network of biopolymers. The resulting cell size scale force chains can in principle mediate mechanical interactions between distant actomyosin units, leading to self-organized structural order in the cell cytoskeleton. Inspired by such force transmission through elastic structures in the cytoskeleton, we consider a percolated fiber lattice network, where fibers are represented as linear elastic elements that can both bend and stretch, and the contractile activity of myosin motors is represented by force dipoles. Then, by using a variety of metrics, we show how two such contractile force dipoles interact with each other through their mutual mechanical deformations of the elastic fiber network. As a prelude to two-dipole interactions, we quantify how forces propagate through the network from a single anisotropic force dipole by analyzing clusters of nodes connected by highly strained bonds, as well as through the decay rate of strain energy with distance from a force dipole. We show that predominant fiber bending screens out force propagation, resulting in reduced and strongly network configuration-dependent dipole interactions. On the other hand, stretching-dominated networks support longer-ranged inter-dipole interactions that recapitulate the predictions of linear elasticity theory. By characterizing the differences between tensile and compressive force propagation in the fiber network, we show how inter-dipole interaction depends on the dipoles’ mutual separation and orientation. The resulting elastic interaction energy may mediate a force between multiple distant dipoles, leading to their self-organization into ordered configurations. In conclusion, this provides a potential pathway for active mechanical force-driven structural order in elastic biopolymer networks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hole doping in compositionally complex correlated oxide enables tunable exchange biasing

Magnetic interfaces and the phenomena arising from them drive both the design of modern spintronics and fundamental research. Recently, it was revealed that through designing magnetic frustration in configurationally complex entropy stabilized oxides, exchange bias can occur in structurally single crystal films. This eliminates the need for complex heterostructures and nanocomposites in the design and control of magnetic response phenomena. In this work, we demonstrate through hole doping of a high entropy perovskite oxide that tuning of magnetic responses can be achieved. With detailed magnetometry, we show magnetic coupling exhibiting a variety of magnetic responses including exchange bias and antiferromagnetic spin reversal in the entropy stabilized ABO 3 perovskite oxide La 1-x Sr x (Cr 0.2 Mn 0.2 Fe 0.2 Co 0.2 Ni 0.2 )O 3 family. We find that manipulation of the A-site charge state can be used to balance magnetic phase compositions and coupling responses. This allows for the creation of highly tunable exchange bias responses. In the low Sr doping regime, a spin frustrated region arising at the antiferromagnetic phase boundary is shown to directly couple to the antiferromagnetic moments of the film and emerges as the dominant mechanism, leading to a vertical shift of magnetization loops in response to field biasing. At higher concentrations, direct coupling of antiferromagnetic and ferromagnetic regions is observed. This tunability of magnetic coupling is discussed within the context of these three competing magnetic phases, revealing critical features in designing exchange bias through exploiting spin frustration and disorder in high entropy oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
NASA NTRS2019

NASA’s Galactic Cosmic Ray Simulator at Brookhaven National Laboratory: Enabling Human Exploration Missions to the Moon and Mars

With exciting new Agency plans for a sustainable return to the moon, astronauts will once again leave earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic rays (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPE’s with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever-penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include epithelial carcinogenesis and leukemia, central nervous system effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including cataracts, circulatory and heart disease, as well as, potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy ion exposures and the validation of countermeasures in a relevant space environment. NASA has developed the “GCR Simulator” at Brookhaven National Laboratory to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (~65-75%) and alpha particles (~10-20%) with heavier ions (Z≤3) contributing the remainder. The “GCR Simulator” exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the five heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires 75-90 minutes. To more closely simulate the low dose rates found in space, sequential field exposures can be divided into daily fractions over 2-4 weeks, with individual fractions as low as 0.1-0.2 mGy. In the large beam configuration (60 x 60 cm(exp 2)), 54 special housing cages can accommodate 2-3 mice each for a 70-75 min duration or ~15 individually housed rats. Emerging research results from our 2018 runs utilizing mixed heavy ion fields and protracted space exposures are forthcoming and deepen our understanding of the numerous health risks faced by our astronauts. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NASA Space Radiation Laboratory to enable future exploration missions.

Lisa C Simonsen↗
NASA NTRSDate not supplied

NASA Space Radiation Laboratory Galactic Cosmic Ray Simulator

With exciting new Agency plans for a sustainable return to the moon, astronauts will once again leave earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic rays (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever-penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including cataracts, circulatory and heart disease, as well as, potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy ion exposures and the validation of countermeasures in a relevant space environment. NASA has developed the “GCR Simulator” at Brookhaven National Laboratory to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (~65-75%) and helium ions (~10-20%) with heavier ions (Z>3) contributing the remainder. The “GCR Simulator” exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the five heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires 75-90 minutes. To more closely simulate the low dose rates found in space, sequential field exposures can be divided into daily fractions over 2-6 weeks, with individual beam fractions as low as 0.1-0.2 mGy. In the large beam configuration (60 x 60 cm2), 54 special housing cages can accommodate 2-3 mice each for a 70-75 min duration or ~15 individually housed rats. Emerging research results from our 2018 runs utilizing mixed heavy ion fields and protracted space exposures are forthcoming and will deepen our understanding of the numerous health risks faced by astronauts. This talk discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NASA Space Radiation Laboratory to enable future exploration missions.

galactic cosmic ray simulator↗
DOE OSTI2023

Correlation between cation order/disorder and the electrocaloric effect in the MLCCs of complex perovskite ferroelectrics

The physical properties (dielectric, ferroelectric, piezoelectric, etc.) of complex perovskite ferroelectrics depend on the degree of order/disorder and the scale of the ordered domains. In this study, the electrocaloric (EC) properties of three representative complex perovskite ferroelectrics, Pb(Mg 1/3 Nb 2/3 )O 3 -8PbTiO 3 (PMN-8PT), 1mol% Sm-doped Pb(Mg 1/3 Nb 2/3 )O 3 -8PbTiO 3 (1S-PMN-8PT) and Pb(Sc 1/2 Ta 1/2 )O 3 (PST) are evaluated. Multi-layer ceramic capacitors (MLCCs) with identical structural configurations were fabricated for these three compounds, and their EC properties were characterized by direct measurement using a thermocouple. The EC temperature change of PMN-8PT, 1S-PMN-8PT and PST MLCCs under 20 V μm –1 at room temperature were found to be 1.42 K, 1.54 K, and 3.10 K, respectively. X-ray diffraction and high-resolution transmission electron microscopy data suggests that the high EC performance of PST is related to the ordering of B-site cations (Sc 3+ and Ta 5+ ) with the ordering parameter S = 0.82 and a long coherence length of ~100 nm, such that the sample transitioned from a relaxor ferroelectric to a normal ferroelectric. Furthermore, these results provide pathway towards design of high performance EC materials required for solid state refrigeration and air-conditioning technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oxygen Electrocatalysis on Mixed-Metal Oxides/Oxyhydroxides: From Fundamentals to Membrane Electrolyzer Technology

Catalyzing the oxygen evolution reaction (OER) is important for key energy-storage technologies, particularly water electrolysis and photoelectrolysis for hydrogen fuel production. Under neutral-to-alkaline conditions, first-row transition-metal oxides/(oxy)hydroxides are the fastest-known OER catalysts and have been the subject of intense study for the past decade. Critical to their high performance is the intentional or accidental addition of Fe to Ni/Co oxides that convert to layered (oxy)hydroxide structures during the OER. Unraveling the role that Fe plays in the catalysis and the molecular identity of the true “active site” has proved challenging, however, due to the dynamics of the host structure and absorbed Fe sites as well as the diversity of local structures in these disordered active phases. In this Account, we highlight our work to understand the role of Fe in Ni/Co (oxy)hydroxide OER catalysts. We first discuss how we characterize the intrinsic activity of the first-row transition-metal (oxy)hydroxide catalysts as thin films by accounting for the contributions of the catalyst-layer thickness (mass loading) and electrical conductivity as well as the underlying substrate’s chemical interactions with the catalyst and the presence of Fe species in the electrolyte. We show how Fe-doped Ni/Co (oxy)hydroxides restructure during catalysis, absorb/desorb Fe, and in some cases degrade or regenerate their activity during electrochemical testing. We highlight the relevant techniques and procedures that allowed us to better understand the role of Fe in activating other first-row transition metals for OER. We find several modes of Fe incorporation in Ni/Co (oxy)hydroxides and show how those modes correlate with activity and durability. We also discuss how this understanding informs the incorporation of earth-abundant transition-metal OER catalysts in anion-exchange-membrane water electrolyzers (AEMWE) that provide a locally basic anode environment but run on pure water and have advantages over the more-developed proton-exchange-membrane water electrolyzers (PEMWE) that use platinum-group-metal (PGM) catalysts. We outline the key issues of introducing Fe-doped Ni/Co (oxy)hydroxide catalysts at the anode of the AEMWE, such as the oxidative processes triggered by Fe species traveling through the polymer membrane, pH-gradient effects on the catalyst stability, and possibly limited catalyst utilization in the compressed stack configuration. We also suggest possible mitigation strategies for these issues. Lastly, we summarize remaining challenges including the long-term stability of Fe-doped Ni/Co (oxy)hydroxides under OER conditions and the lack of accurate models of the dynamic active surface that hinder our understanding of, and thus ability to design, these catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bottlebrush Block Copolymer Thin Films

The self-assembly of block copolymers (BCPs) in thin films is governed by interfacial interactions (enthalpy) between the blocks and interfaces (at both the substrate and the surface). Advantage can also be taken of entropy to control the polymer orientation. Here, we synthesized a series of bottlebrush block copolymers (BBCPs), with deuterated polystyrene (DPS) as the core block and poly­(solketal acrylate) (PSA) as the corona block, where the backbone length (N BB) and grafting densities (GDs) were varied. The hydrophobic PSA block was converted to a hydrophilic poly­(glyceryl acrylate) (PGA) block by solid-state hydrolysis, bringing the BBCP from a disordered state into an ordered state with a lamellar microdomain morphology. The orientation of the morphology as a function of distance from the interfaces was systematically investigated by atomic force microscopy (AFM), interference microscopy, grazing-incidence small-angle X-ray scattering (GISAXS), and grazing-incidence small angle neutron scattering (GISANS). For N BB = 1 (diblock BCP), a mixed lamellar orientation was found. For N BB = 5 (star-like BBCP), a vertical orientation of the lamellar microdomains was found dominantly at the interfaces that propagated into the bulk of the film due to the entropic penalty associated with constraints on the junction points. Increasing the N BB to 50 (rod-like BBCP) and 100 (worm-like BBCP) yielded a random orientation, as the side chains became overcrowded when aligned in the same direction. By reducing the grafting density (GD) of a N BB = 50 BBCP from 100% to 75% and 50% by the incorporation of spacers into the backbone, the steric hindrance of side chains is reduced markedly, facilitating a vertical orientation of the lamellar microdomains. However, the vertical orientation becomes more random within the bulk of the thin film for different BBCP architectures. By modifying the substrate with a layer of PS, the preferential interactions between the core block of the BBCP and the interfaces induce a reconfiguration of the BBCP, drawing the core block out to contact substrate directly driven by enthalpy. Under these conditions, N BB = 5 showed a predominantly parallel orientation of the lamellar microdomains with the corona block minimizing contact with the substrate. Higher N BB values of 50 and 100 showed a more random orientation of the lamellar microdomains at different depths on the PS-modified substrate. Overall, thin-film morphology was regulated by tuning the configuration of the BBCP, varying the entropic contribution to the microdomain orientation.

Chen, Zhan↗
DOE OSTI2020

NASA’s first ground-based Galactic Cosmic Ray Simulator: Enabling a new era in space radiobiology research

With exciting new NASA plans for a sustainable return to the moon, astronauts will once again leave Earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic radiation (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system (CNS) effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including circulatory and heart disease, as well as potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy-ion exposures and the validation of countermeasures in a relevant space environment. Historically, most research on understanding space radiation-induced health risks has been performed using acute exposures of monoenergetic single-ion beams. However, the space radiation environment consists of a wide variety of ion species over a broad energy range. Using the fast beam switching and controls systems technology recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory, a new era in radiobiological research is possible. NASA has developed the “GCR Simulator” to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (approximately 65%–75%) and helium ions (approximately 10%–20%) with heavier ions (Z ≥ 3) contributing the remainder. The GCR simulator exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the 5 heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader system is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low-energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires approximately 75 minutes. To more closely simulate the low-dose rates found in space, sequential field exposures can be divided into daily fractions over 2 to 6 weeks, with individual beam fractions as low as 0.1 to 0.2 mGy. In the large beam configuration (60 × 60 cm 2 ), 54 special housing cages can accommodate 2 to 3 mice each for an approximately 75 min duration or 15 individually housed rats. On June 15, 2018, the NSRL made a significant achievement by completing the first operational run using the new GCR simulator. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NSRL to accelerate our understanding and mitigation of health risks faced by astronauts. Ultimately, the GCR simulator will require validation across multiple radiogenic risks, endpoints, doses, and dose rates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

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↗
DOE OSTI2025

Minimal cyclic behavior in sheared amorphous solids

Although jammed packings of soft spheres exist in potential energy landscapes with a vast number of minima, when subjected to cyclic shear they may revisit the same configurations repeatedly. Simple hysteretic spin models, in which particle rearrangements are represented by interacting spin flips called hysterons, capture many features of this periodic behavior. Yet it has been unclear to what extent individual rearrangements can be described by such binary objects and how such objects interact with one another. Using a particularly sensitive algorithm, we identify rearrangements in simulated jammed packings and select pairs of rearrangements that undo one another to create periodic cyclic behavior. We find that the rearrangement pairs surprisingly persist down to the smallest increments in strain, even in the smallest systems we can study. We explore the statistics of these rearrangement pairs and find that there is a relation between the amount of hysteresis and the energy drop and mean-square displacement of the particles; these results are inconsistent with the scaling found in models that treat rearrangements as localized buckling events. Finally, our analysis shows that there is no clean distinction between the particle motions that represent the identity of a single, individual rearrangement and the particle motions that lead to interactions between separated rearrangements or hysterons. These results offer insight into how complex systems such as amorphous solids can reach a limit cycle.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Advanced sorbents for modular oxygen production for REMS gasifiers

Under a DOE funded effort, Thermosolv LLC has been developing a sorbent-based oxygen production technology for gasification and oxy-combustion applications. The sorbents utilize oxygen-storage properties of certain perovskites to (1) selectively adsorb oxygen at moderate temperature from compressed air and (2) release the adsorbed oxygen into a vacuum or a sweep gas such as CO2 and/or steam. Through cyclic operations of multiple sorbent beds such a process can be made continuous. Pressure drop across the sorbent bed and other similar design considerations dictate that the sorbent be in the form of a high surface-area-to-volume pellet and yet still possess sufficient crush strength and integrity. Process efficiency is determined by the sorption and desorption kinetics and sorbent capacity. Typical process involves cycle times in the order of 100 seconds, a time too short to fully utilize the full volume of the sorbent and the sorption capacity of the sorbent pellets. As a part of this project, Thermosolv LLC undertook the development of oxygen sorbents as a high-surface area supported sorbent on a light-weight inert support to (1) reduce the cost, (2) increase the productivity, and (3) reduce the overall weight of the reactor. Using a previously developed perovskite (LSCF 1991) with its well-characterized performance, composite sorbent pellets each consisting of an inert core coated by a thin layer of the functional perovskite material were produced. Several low-density, inexpensive inert solids supports in the 1/8”-1/4” size were selected to keep the pressure drop across the sorbent bed in the acceptable range. A number of commercial technologies including spray coating, dip coating, incipient impregnation of porous substrate followed by thermal annealing at various temperatures and duration were employed and tested to produce robust composite supported sorbent pellets. A sense of optimum coating thickness was developed by testing sorbent pellets of various increasing diameter pellets. LSCF-1991extrudates ranging in size from 1/32” to 3/16” were tested in a TGA and in a fixed-bed pressure swing test set-up for sorption/desorption cycles of interest. The data show that for the operational conditions of interest the coating thickness for a supported sorbent pellet needs to be approximately 1/32” (about 0.8 mm). Subsequent work thereby concentrated on developing composite pellets of various substrates, shapes and sizes with coating layer of about this size range. Candidate support materials were chosen based on cost, inertness, mechanical strength, thermal expansion and chemical stability with sorbent material, and in a size range to give an acceptable pressure drop in fixed-bad reactor configurations. Coating application methods used for application of the sorbent onto the support included precipitation, spray coating and dip coating from sorbent slurries and sol gels. For all coated supports where we could successfully apply a uniform coating of desired thickness with an acceptable handling performance in terms of exfoliation, TGA-based cyclic sorption/desorption testing was performed to determine cycling capacity of oxygen. Among nearly twenty different support materials tested, best adhesion performance was obtained from stainless saddle supports. In the bench-scale fixed-bad tests stainless steel saddle support composite pellets with a coating thickness in the 0.6 mm or so range, the composite sorbent pellet performance approached up to 95% of that of the 2 mm parent material pellets. In a parallel approach to reducing the cost of the sorbent, alternate sorbent formulations replacing/reducing the amount of cobalt in the LSCF family were also investigated. Successful formulations that could match the performance of LSCF 1991 were identified based on TGA cyclic tests as LSCF 1919 and LSF 1910. In the range of operational envelop of cycle times and other relevant process operating conditions, the reduced Co formulations showed comparable performance in the bench-scale fixed-bad cyclic operations. As a part of this project, we also attempted substitution of La and Sr with Ca and Ba, and Co with Mn, Ni and Cu with little success, but the overall project goal of reducing sorbent cost in terms of raw material and manufacturing expenses was successful.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗