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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2021

Monte carlo simulations and experiments of all-silica zeolite LTA assembly combining structure directing agents that match cage sizes

In this work, we investigated the influence of organic structure-directing agents (OSDAs) on the formation rates of all-silica zeolite LTA using both simulations and experiments, to shed light on the crystallization process. We compared syntheses using one OSDA with a diameter close to the size of the large cavity in LTA, and two OSDAs of diameters matching the sizes of both the small and large LTA cavities. Reaction-ensemble Monte Carlo (RxMC) simulations predict a speed up of LTA formation using two OSDAs matching the LTA pore sizes; this qualitative result is confirmed by experimental studies of crystallization kinetics, which find a speedup in all-silica LTA crystallization of a factor of 3. Analyses of simulated rings and their Si–O–Si angular energies during RxMC crystallizations show that all ring sizes in the faster crystallization exhibit lower angular energies, on average, than in the slower crystallization, explaining the origin of the speedup through packing effects.

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

Assessment of Grizzly Capabilities for Reactor Pressure Vessels and Reinforced Concrete Structures

Over the last several years, capabilities to simulate the progression and effects of degradation in critical structures in light water reactor (LWR) nuclear power plants have been under development in the Grizzly code. Age-related material degradation is important for a number of systems in LWRs, but the main focus for Grizzly development has been on reactor pressure vessels (RPVs) and concrete structures because of their central role and the difficulty of replacement of these structures if they are found to be degraded to an unacceptable degree. The capabilities for analyzing both RPVs and concrete structures in Grizzly have reached a point where the feature sets are sufficiently complete to perform credible analyses of those types of structures. Because of this, the emphasis has shifted from foundational development to assessing the accuracy of these modeling capabilities on representative problems of interest and on improving the physical basis of those models to improve their ability to predict the actual response of those structural systems. This report documents a first set of test cases that have been developed to assess the Grizzly code on real-world problems for these two types of structures. For RPVs, these test cases consist of a set of benchmark problems where Grizzly is compared to another code. For concrete structures, a set of models of experimental specimens designed to characterize the multidimensional swelling response of reinforced concrete members due to alkali-silica reaction (ASR) has been developed and compared with experimental results. Both the RPV and concrete test cases developed here are much more computationally intensive than the small regression tests that are used in the testing that Grizzly undergoes every time a proposed set of changes is made to ensure that those changes do not adversely affect previously established behavior. A system for regularly running these large problems to monitor the behavior of the code as it is developed has been instituted. While these test cases generally indicate good comparison with the benchmark results, they also indicate areas where further development is warranted.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

DYNAMICS OF TURBULENT CONVECTION AND CONVECTIVE OVERSHOOT IN A MODERATE-MASS STAR

We present results of realistic three-dimensional (3D) radiative hydrodynamic simulations of the outer layers of a moderate-mass star (1.47 M {sub ⊙}), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding one-dimensional (1D) standard stellar model shows an increase of the stellar radius by ∼800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km s{sup −1}, penetrate through the whole convection zone, hit the radiative zone, and form an 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes. This effect has to be taken into account in stellar modeling and the interpretation of asteroseismology data. In particular, we demonstrate that the deviations of the mean structure of the 3D model from the 1D standard model of the same mass and composition are qualitatively similar to the deviations for the Sun found by helioseismology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Novel Method to Train Classification Models for Structure Detection in In Situ Spacecraft Data

We present a method for creating spacecraft-like data which can be used to train Machine Learning (ML) models to detect and classify structures in in situ spacecraft data. First, we use the Grad-Shafranov equation to numerically solve for several magnetohydrostatic equilibria which are variations on a known analytic equilibrium. These equilibria are then used as the initial conditions for Particle-In-Cell simulations in which the structures of interest are observed and labeled. We then take one-dimensional slices through the simulations to replicate what a spacecraft collecting data from the simulation would observe. This sliced data then can be used as training data for the initial training of ML models intended for use on spacecraft data. We demonstrate the method applied to the problem of detecting small-scale plasmoids in the magnetotail, which is important for understanding complex magnetotail reconnection dynamics. The simple 1D classifier we train is able to detect more than 70% of the plasmoid points in the data set but also produces a large number of false positives. Our further work on this example problem is detailed, and further potential uses of the method are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Simultaneous coupling of surface plasmon resonance and photonic bandgap to InGaAs quantum well emission

A photonic bandgap structure was created on the 100 nm thick GaAs barrier layer with Au nanodisks deposited inside the holes. To mitigate the nonradiative surface recombination of GaAs, the Au nanodisks were formed on top of a 15 nm SiO{sub 2} deposited in the holes. A maximum 7.6-fold increase in photoluminescence intensity was obtained at the etch depth of 80 nm. In this configuration, the Au nanodisk is separated from the quantum well by 20 nm of GaAs and 15 nm of SiO{sub 2}. The experimental result was verified by the simulation based on this structure. There was a good agreement between experiments with simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Modeling the co-assembly of binary nanoparticles

Abstract In this work, we present a binary assembly model that can predict the co-assembly structure and spatial frequency spectra of monodispersed nanoparticles with two different particle sizes. The approach relies on an iterative algorithm based on geometric constraints, which can simulate the assembly patterns of particles with two distinct diameters, size distributions, and at various mixture ratios on a planar surface. The two-dimensional spatial-frequency spectra of the modeled assembles can be analyzed using fast Fourier transform analysis to examine their frequency content. The simulated co-assembly structures and spectra are compared with assembled nanoparticles fabricated using transfer coating method are in qualitative agreement with the experimental results. The co-assembly model can also be used to predict the peak spatial frequency and the full-width at half-maximum bandwidth, which can lead to the design of the structure spectra by selection of different monodispersed particles. This work can find applications in fabrication of non-periodic nanostructures for functional surfaces, light extraction structures, and broadband nanophotonics.

Mohanty, Saurav (ORCID:0009000302592855)↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2021

First Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the sensor was initially investigated with analytical derivations and with COMSOL RF Module computer simulations of resonant frequency spectral shift due to uniform load. We also investigated the mechanical integrity of the flowmeter’s membrane through analytical modelling and COMSOL Structural Mechanics Module computer simulations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the material’s yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure. Using results from models, we have developed an initial design for a microwave K-band sensor. A cylindrical resonator prototype was fabricated from brass for the initial tests. The external dimensions of the cavity are matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed consisting of a piping Tee with a bulkhead WR-42 microwave waveguide installed in a leak-proof assembly. A microwave waveguide circulator was installed in the setup to suppress the effect of reflections at the cavity entrance by increasing the isolation between the input and the output port. Preliminary spectral characterization of cavity spectral response was performed with a portable PXIe chassis microwave VNA with a custom GUI. Preliminary dry tests of the transducer response were conducted with a set of calibrated weights. Transducer frequency shift was shown to be monotonically increasing with increasing pressure. The next steps will involve investigation of the transducer performance for water flow sensing.

42 ENGINEERING↗
DOE OSTI2015

Ligand structure and mechanical properties of single-nanoparticle thick membranes

We believe that the high mechanical stiffness of single-nanoparticle-thick membranes is the result of the local structure of ligand coatings that mediate interactions between nanoparticles. These ligand structures are not directly observable experimentally. We use molecular dynamics simulations to observe variations in ligand structure and simultaneously measure variations in membrane mechanical properties. We have shown previously that ligand end group has a large impact on ligand structure and membrane mechanical properties. Here we introduce and apply quantitative molecular structure measures to these membranes and extend analysis to multiple nanoparticle core sizes and ligand lengths. Simulations of nanoparticle membranes with a nanoparticle core diameter of 4 or 6 nm, a ligand length of 11 or 17 methylenes, and either carboxyl (COOH) or methyl (CH 3 ) ligand end groups are presented. In carboxyl-terminated ligand systems, structure and interactions are dominated by an end-to-end orientation of ligands. In methyl-terminated ligand systems large ordered ligand structures form, but nanoparticle interactions are dominated by disordered, partially interdigitated ligands. Core size and ligand length also affect both ligand arrangement within the membrane and the membrane's macroscopic mechanical response, but are secondary to the role of the ligand end group. Additionally, the particular end group (COOH or CH 3 ) alters the nature of how ligand length, in turn, affects the membrane properties. The effect of core size does not depend on the ligand end group, with larger cores always leading to stiffer membranes. Asymmetry in the stress and ligand density is observed in membranes during preparation at a water-vapor interface, with the stress asymmetry persisting in all membranes after drying.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Cooling rate dependence of structural order in Ni 62 Nb 38 metallic glass

In this article, molecular dynamics (MD) simulations are performed to study the structure of Ni 62 Nb 38 bulk metallic glass at the atomistic level. Structural analysis based on the cluster alignment method is carried out and a new Ni-centered distorted-icosahedra (DISICO) motif is excavated. We show that the short-range order and medium-range order in the glass are enhanced with lower cooling rate. Almost 50% of the clusters around the Ni atoms in the well-annealed Ni 62 Nb 38 glass sample from our MD simulations can be classified as DISICO. It is revealed that the structural distortion with respect to the perfect icosahedra is driven by chemical ordering in the distorted region of the DISICO motif. The relationship between the structure, energy, and dynamics in this glass-forming alloy during the cooling and annealing processes is also established.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Influence of heterogeneous sulfur atoms on the negative differential resistance effect in polythiophene

In this work, we have carried out theoretical investigations aiming to clarify the effects of sulfur heteroatoms on the transport characteristics in polythiophene. Sulfur atoms in polythiophene are demonstrated to influence the structure and transport process by two aspects: the electron hopping between carbon atoms on both sides of the sulfur atom as well as the effective confinement of π electrons from the sulfur atom. Based on the static Su-Schrieffer-Heeger model and the nonequilibrium Green's function formalism, we simulate the electron transportation in a metal/polythiophene/metal structure. The simulation results show that the electron hopping via sulfur atoms is responsible for the observed negative differential resistance (NDR) behavior in the I-V curves. The NDR disappears if the electron transport channels from carbon to carbon via sulfur atoms are forbidden. The weaker the effective confinement of π electrons and the electron hopping between carbon atoms on both sides of the sulfur atom are, the higher is the peak-to-valley ratio of the NDR and the wider the voltage range where the current remains at low levels. These results can help in understanding the NDR effect in polythiophene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗
DOE OSTI2014

Stoichiometry determined exchange interactions in amorphous ternary transition metal oxides: Theory and experiment

Amorphous transition metal oxides exhibit exotic transport and magnetic properties, while the absence of periodic structure has long been a major obstacle for the understanding of their electronic structure and exchange interaction. In this paper, we have formulated a theoretical approach, which combines the melt-quench approach and the spin dynamic Monte-Carlo simulations, and based on it, we explored amorphous Co{sub 0.5}Zn{sub 0.5}O{sub 1−y} ternary transition metal oxides. Our theoretical results reveal that the microstructure, the magnetic properties, and the exchange interactions of Co{sub 0.5}Zn{sub 0.5}O{sub 1−y} are strongly determined by the oxygen stoichiometry. In the oxygen-deficient sample (y > 0), we have observed the long-range ferromagnetic spin ordering which is associated with the non-stoichiometric cobalt-rich region rather than metallic clusters. On the other hand, the microstructure of stoichiometric sample takes the form of continuous random networks, and no long-range ferromagnetism has been observed in it. Magnetization characterization of experimental synthesized Co{sub 0.61}Zn{sub 0.39}O{sub 1−y} films verifies the relation between the spin ordering and the oxygen stoichiometry. Furthermore, the temperature dependence of electrical transport shows a typical feature of semiconductors, in agreement with our theoretical results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Ideal plasma response to vacuum magnetic fields with resonant magnetic perturbations in non-axisymmetric tokamaks

Ideal plasma shielding and amplification of resonant magnetic perturbations in non-axisymmetric tokamak is presented by field line tracing simulation with full ideal plasma response, compared to measurements of divertor lobe structures. Magnetic field line tracing simulations in NSTX with toroidal non-axisymmetry indicate the ideal plasma response can significantly shield/amplify and phase shift the vacuum resonant magnetic perturbations. Ideal plasma shielding for n = 3 mode is found to prevent magnetic islands from opening as consistently shown in the field line connection length profile and magnetic footprints on the divertor target. It is also found that the ideal plasma shielding modifies the degree of stochasticity but does not change the overall helical lobe structures of the vacuum field for n = 3. Furthermore, amplification of vacuum fields by the ideal plasma response is predicted for low toroidal mode n = 1, better reproducing measurements of strong striation of the field lines on the divertor plate in NSTX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

V67L Mutation Fills an Internal Cavity To Stabilize RecA Mtu Intein

Inteins mediate protein splicing, which has found extensive applications in protein science and biotechnology. In the Mycobacterium tuberculosis RecA mini–mini intein (ΔΔIhh), a single valine to leucine substitution at position 67 (V67L) dramatically increases intein stability and activity. However, crystal structures show that the V67L mutation causes minimal structural rearrangements, with a root-mean-square deviation of 0.2 Å between ΔΔIhh-V67 and ΔΔIhh-L67. Thus, the structural mechanisms for V67L stabilization and activation remain poorly understood. In this paper, we used intrinsic tryptophan fluorescence, high-pressure nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations to probe the structural basis of V67L stabilization of the intein fold. Guanidine hydrochloride denaturation monitored by fluorescence yielded free energy changes (ΔG f °) of -4.4 and -6.9 kcal mol –1 for ΔΔIhh-V67 and ΔΔIhh-L67, respectively. High-pressure NMR showed that ΔΔIhh-L67 is more resistant to pressure-induced unfolding than ΔΔIhh-V67 is. The change in the volume of folding (ΔV f ) was significantly larger for V67 (71 ± 2 mL mol –1 ) than for L67 (58 ± 3 mL mol –1 ) inteins. The measured difference in ΔV f (13 ± 3 mL mol –1 ) roughly corresponds to the volume of the additional methylene group for Leu, supporting the notion that the V67L mutation fills a nearby cavity to enhance intein stability. In addition, we performed MD simulations to show that V67L decreases side chain dynamics and conformational entropy at the active site. Finally, it is plausible that changes in cavities in V67L can also mediate allosteric effects to change active site dynamics and enhance intein activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗