Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical Simulations of Flow-Induced Deflections in MITR LEU Fuel Plate Due to Channel Size Disparity

The hydromechanical stability of the fuel plates in parallel coolant channels of a Materials Testing Reactor (MTR) fuel element design is of great importance to the safety of research and test reactors. Previous analytical, experimental, and numerical efforts focused on parallel channels with the same or similar size; also, in the prior numerical simulations, the fuel plate was often assumed to be perfectly flat. This work presents the results of a fluid-structure interaction simulation performed to evaluate the flow-induced deflections of the fuel plates in the low-enriched uranium (LEU, <20 wt% 235 U) fuel element design for the conversion (from highly enriched uranium) of the Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR). Various manufacturing and assembly tolerances of the MITR LEU elements are considered in the analysis, and the effects of channel size disparity, nonideal plate shape, and flow rate uncertainty are investigated. Results show that, for all cases analyzed, the deflection occurs toward the larger channel, and the change in any channel stripe remains small (less than 0.021 mm) compared to fabrication tolerances. In addition to simulation work, a hydraulic performance test of the MITR LEU fuel element is currently planned to support conversion to the use of LEU fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Medium-Range Ordering in the Ionic Glass Electrolytes LiPON and LiSiPON

In this work, we provide an in-depth structural characterization of the amorphous ionic glasses LiPON and LiSiPON with high Li content. Based on ab initio molecular dynamics simulations, the structure of these materials is an inverted structure with either isolated polyanion tetrahedra or polyanion dimers in a Li+ matrix. Based on neutron scattering data, this type of inverted structure leads to a significant amount of medium-range ordering in the structure, as demonstrated by two sharp diffraction peaks and a periodic structural oscillation in the density function G(r). While this medium-range ordering is commonly observed in liquids and metallic glasses, it has not previously been observed in oxides. On a local scale, adding N and Si increases the number of anion bridges and polyanion dimer structures, leading to higher ionic conductivity. In the medium-range ordering, the addition of Si leads to more disorder in the polyanion substructure but a significant increase in the ordering of the O substructure. Finally, we demonstrate that this inverted structure with medium-range ordering results in a glassy material that is both mechanically stiff and ductile on the nanoscale.

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

PDADMAC/PSS Oligoelectrolyte Multilayers: Internal Structure and Hydration Properties at Early Growth Stages from Atomistic Simulations

We analyze the internal structure and hydration properties of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) oligoelectrolyte multilayers at early stages of their layer-by-layer growth process. Our study is based on large-scale molecular dynamics simulations with atomistic resolution that we presented recently [Sánchez et al., Soft Matter 2019, 15, 9437], in which we produced the first four deposition cycles of a multilayer obtained by alternate exposure of a flat silica substrate to aqueous electrolyte solutions of such polymers at 0.1M of NaCl. In contrast to any previous work, here we perform a local structural analysis that allows us to determine the dependence of the multilayer properties on the distance to the substrate. We prove that the large accumulation of water and ions next to the substrate observed in previous overall measurements actually decreases the degree of intrinsic charge compensation, but this remains as the main mechanism within the interface region. We show that the range of influence of the substrate reaches approximately 3 nm, whereas the structure of the outer region is rather independent from the position. This detailed characterization is essential for the development of accurate mesoscale models able to reach length and time scales of technological interest.

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

Electron tomography unravels new insights into fiber cell wall nanostructure; exploring 3D macromolecular biopolymeric nano-architecture of spruce fiber secondary walls

Lignocellulose biomass has a tremendous potential as renewable biomaterials for fostering the “bio-based society” and circular bioeconomy paradigm. It requires efficient use and breakdown of fiber cell walls containing mainly cellulose, hemicellulose and lignin biopolymers. Despite their great importance, there is an extensive debate on the true structure of fiber walls and knowledge on the macromolecular nano-organization is limited and remains elusive in 3D. We employed dual-axis electron tomography that allows visualization of previously unseen 3D macromolecular organization/biopolymeric nano-architecture of the secondary S2 layer of Norway spruce fiber wall. Unprecedented 3D nano-structural details with novel insights into cellulose microfibrils (~2 nm diameter), macrofibrils, nano-pore network and cell wall chemistry (volume %) across the S2 were explored and quantified including simulation of structure related permeability. Matrix polymer association with cellulose varied between microfibrils and macrofibrils with lignin directly associated with MFs. Simulated bio-nano-mechanical properties revealed stress distribution within the S2 and showed similar properties between the idealized 3D model and the native S2 (actual tomogram). Present work has great potential for significant advancements in lignocellulose research on nano-scale understanding of cell wall assembly/disassembly processes leading to more efficient industrial processes of functionalization, valorization and target modification technologies.

3-D reconstruction↗
DOE OSTI2021

Rate-Based Absorption Modeling for Postcombustion CO 2 Capture with Additively Manufactured Structured Packing

Carbon capture using amine-based solvents in an absorption process is a leading candidate for reducing greenhouse gas emissions in industrial flue gas streams. To reduce operating costs and associated parasitic energy of these processes, process intensification utilizing additively manufactured structured packing has emerged as a new technology to manage exothermic reactions during absorption while improving CO 2 capture. A rate-based model framework has been developed for these novel packings that incorporates the mass and heat transfer phenomena for amine-based absorption of CO 2 . The rate-based model framework is first benchmarked using available solubility data and pilot plant data for aqueous monoethanolamine. The model validation shows accurate prediction of both CO 2 equilibrium partial pressures and ion speciation for solubility data as well as CO 2 capture, temperature profile, and solvent composition for pilot plant data in the absorption column. The model framework is then applied toward predicting the CO 2 capture performance of additively manufactured structured packing. Simulations agree with experimental data in predicting the CO 2 capture and the capture performance increase due to cooling within the structured packing device. Finally, advantages of this rate-based model framework are the utilization of correlations that may predict mass transfer and heat transfer coefficients of the packing based on the geometric properties of the device and the implementation of the model framework in open-source programming.

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

X-ray Absorption Spectral Signature of Quantum Nuclei in Liquid Water

Based on electron-hole excitation theory, we investigate the X-ray absorption spectral signature of nuclear quantum effect in liquid water, whose molecular structure is simulated by path-integral molecular dynamics using the MB-pol model. Compared to spectra generated from classically modeled water structure, quantum nuclei has important effect on spectra in terms of both the spectral energies and their line shapes. At the short-range ordering of H-bond network, the delocalized protons increase the fluctuations on the intramolecular covalency and broaden the pre-edge of the spectra. For intermediate-range and long-range orderings, the observed red and blue shifts of the main-edge and post-edge are attributed to the so-called competing quantum effects, under which both the weak and well-formed H-bonds are promoted. The theoretical spectra are in nearly quantitative agreement with the available experimental data.

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

Construction of an Anion-Pillared MOF Database and the Screening of MOFs Suitable for Xe/Kr Separation

The separation of xenon/krypton (Xe/Kr) mixtures is a challenging process. Many porous materials allow the adsorption of both Xe and Kr but only with low selectivity. Anion-pillared metal–organic frameworks (MOFs), featuring the anion groups as structural pillars, show potential in gas separations, but only a limited number of them have been synthesized. Here, we describe a collection of 936 anion-pillared MOFs based on 22 experimentally available structures. We performed density functional theory (DFT) optimization and then assigned density-derived electrostatic and chemical (DDEC) charges for each MOF to make them well suited to molecular simulations. The structural properties of the MOFs vary more strongly with the choice of the organic ligand than with other aspects like fluorine groups and metal centers. We then screened the entire collection of MOFs in the context of Xe/Kr separation at room temperature. Compared with previously reported MOFs, the interpenetrated MOF SIFSIX-6-Cd-i is predicted to perform better for Xe/Kr separations, with a good balance between working capacity (1.62 mmol/g) and separation selectivity (16.4) at 298 K and 100 kPa. Here, we also found that the heterogeneity of fluorine groups within a MOF can help to enhance Xe working capacity without reducing the Xe/Kr selectivity, suggesting that synthesis of anion-pillared MOFs with mixed fluorine groups may lead to improved Xe/Kr separation performance.

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

A two-dimensional type I superionic conductor

Superionic conductors (SICs) possess liquid-like ionic diffusivity in the solid state, finding wide applicability from electrolytes in energy storage to materials for thermoelectric energy conversion. Type I SICs (e.g., AgI, Ag 2 Se, etc.) are defined by a first-order transition to the superionic state and have so far been found exclusively in three-dimensional crystal structures. Here, we reveal a two-dimensional type I SIC, α-KAg 3 Se 2 by scattering techniques and complementary simulations. Quasi-elastic neutron scattering and ab initio molecular dynamics simulations confirm that the superionic Ag + ions are confined to sub-nanometre sheets, with the simulated local structure validated by experimental X-ray powder pair-distribution-function analysis. Finally, we demonstrate that the phase transition temperature can be controlled by chemical substitution of the alkali metal ions that comprise the immobile charge-balancing layers. Our work thus extends the known classes of SICs and will facilitate the design of new materials with tailored ionic conductivities and phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

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

Dependence of simulated radiation damage on crystal structure and atomic misfit in metals

This study investigates the evolution of radiation damage in three metals in the low temperature and high radiant flux regime using molecular statics and a Frenkel pair accumulation method to simulate up to 2.0 displacements per atom. The metals considered include Fe, equiatomic CrCoNi, and a fictitious metal with similar bulk properties to the CrCoNi composed of a single atom type referred to as an A-atom. CrCoNi is found to sustain higher concentrations of dislocations than either the Fe or A-atom systems and more stacking faults than the A-atom system. The results suggest that the difference between the concentrations of vacancies and interstitials is substantially smaller for CrCoNi than the A-atom system, perhaps reflecting that the sink capture radius is smaller in CrCoNi due to the roughened potential energy landscape. A model that partitions the major contributions from defects to the stored energy is described, and serves to highlight a general need for higher fidelity approaches to point defect identification.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulation and characterization of cation disorder in ZnGeP 2

New optoelectronic materials are needed for improving the efficiency and reliability of devices such as solar cells. Cation ordering presents one means of controlling optoelectronic properties while introducing potential to also diversify the mineral constituents of electronic devices; however, the mechanisms of ordering are not yet well understood. To better understand cation ordering in a system integratable with current devices, we assess short- and long-range order parameters of ZnGeP 2 , a material closely lattice matched to Si. Structures are simulated using cluster-based Monte Carlo and first-principles calculations to compare structural distortions, periodicity, and local coordination environments in ZnGeP 2 to experimental data both from the literature and presented here. Comparing order parameters, we relate the transition in order parameters of ZnGeP 2 to that of ZnGeP 2 , discuss the reduction of band gaps with disorder, and show that traditional structural characterization alone is insufficient for understanding order in ZnGeP 2 .

14 SOLAR ENERGY↗
DOE OSTI2021

A deep learning interatomic potential developed for atomistic simulation of carbon materials

Interatomic potentials based on neural-network machine learning method have attracted considerable attention in recent years owing to their outstanding ability to balance the accuracy and efficiency in atomistic simulations. In this work, a neural-network potential (NNP) for carbon is generated to simulate the structural properties of various carbon structures. The potential is trained using a database consisting of crystalline and liquid structures obtained by the first-principles density functional theory (DFT) calculations. The developed potential accurately predicts the energies and forces in crystalline and liquid carbon structures, the energetic stability of defected graphene, and the structures of amorphous carbon as the function of density. As a result, the excellent accuracy and transferability of the NNP provide a promising tool for accurate atomistic simulations of various carbon materials with faster speed and much lower cost.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simulation of the Inelastic Alpha-Particle Scattering Reaction for Investigating the Cluster Structure of Excited States of {sup 9}Be

It is proposed to investigate inelastic scattering of alpha particles on the clusters of the {sup 9}Be nucleus for determining the contribution from the cluster configurations ({sup 8}Be+n, α+α+n, {sup 5}He+α ) to the {sup 9}Be structure. The kinematic simulation of the {sup 9}Be(α,α {sup ′}n)X and {sup 9}Be(α,α {sup ′}α)X reactions shows that detection of the scattered alpha particle in coincidence with the particle from the breakup of a certain excited state of {sup 9}Be allows determining the contribution from different cluster configurations to the structure of those states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Microbial Community Structure and Ecological Networks during Simulation of Diatom Sinking

Microbial-mediated utilization of particulate organic matter (POM) during its downward transport from the surface to the deep ocean constitutes a critical component of the global ocean carbon cycle. However, it remains unclear as to how high hydrostatic pressure (HHP) and low temperature (LT) with the sinking particles affects community structure and network interactions of the particle-attached microorganisms (PAM) and those free-living microorganisms (FLM) in the surrounding water. In this study, we investigated microbial succession and network interactions in experiments simulating POM sinking in the ocean. Diatom-derived 13C- and 12C-labeled POM were used to incubate surface water microbial communities from the East China Sea (ECS) under pressure (temperature) of 0.1 (25 °C), 20 (4 °C), and 40 (4 °C) MPa (megapascal). Our results show that the diversity and species richness of the PAM and FLM communities decreased significantly with HHP and LT. Microbial community analysis indicated an increase in the relative abundance of Bacteroidetes at high pressure (40 MPa), mostly at the expense of Gammaproteobacteria, Alphaproteobacteria, and Gracilibacteria at atmospheric pressure. Hydrostatic pressure and temperature affected lifestyle preferences between particle-attached (PA) and free-living (FL) microbes. Ecological network analysis showed that HHP and LT enhanced microbial network interactions and resulted in higher vulnerability to networks of the PAM communities and more resilience of those of the FLM communities. Most interestingly, the PAM communities occupied most of the module hubs of the networks, whereas the FLM communities mainly served as connectors of the modules, suggesting their different ecological roles of the two groups of microbes. These results provided novel insights into how HHP and LT affected microbial community dynamics, ecological networks during POM sinking, and the implications for carbon cycling in the ocean.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Insights on the Alumina–Water Interface Structure by Direct Comparison of Density Functional Simulations with X-ray Reflectivity

Density functional theory molecular dynamics (DFT-MD) simulations are frequently used to predict the interfacial structures and dynamical processes at solid-water interfaces in efforts to gain a deeper understanding of these systems. However, the accuracy of these predictions of interfacial structure has not been rigorously quantified. Here, direct comparisons between large-scale DFT-MD simulations and high-resolution X-ray reflectivity (XR) measurements of the well-defined Al2O3(001)-water interface reveal the relative accuracy of these two methods to describe interfacial structure, a comparison that is enabled by XR’s high sensitivity to atomic-scale displacements. The DFT-MD simulated and XR model-fit structures are qualitatively similar, but XR signals calculated directly from the DFT-MD predictions deviate significantly from the experimental data, revealing discrepancies in these two approaches. Differences in the derived interfacial Al2O3 relaxation profiles of ~0.02 A within the top 5 layers are significant to XR but at the limit of the accuracy of DFT. Further differences are found in the surface hydration layer with a simulated average water layer height ~0.2 A higher than that observed experimentally. This is outside the accuracy of both XR and DFT and is not improved by the inclusion of a phenomenological correction for hydrogen bonding (e.g., Grimme).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗