Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Multi-mode evaluation of power-maximizing cross-flow turbine controllers

A general method for predicting and evaluating the performance of three candidate cross-flow turbine power-maximizing controllers is presented in this paper using low-order dynamic simulation, scaled laboratory experiments, and full-scale field testing. For each testing mode and candidate controller, performance metrics quantifying energy capture (ability of a controller to maximize power), variation in torque and rotation rate (related to drive train fatigue), and variation in thrust loads (related to structural fatigue) are quantified for two purposes. First, for metrics that could be evaluated across all testing modes, we considered the accuracy with which simulation or laboratory experiments could predict performance at full scale. Second, we explored the utility of these metrics to contrast candidate controller performance. For these turbines and set of candidate controllers, energy capture was found to only differentiate controller performance in simulation, while the other explored metrics were able to predict performance of the full-scale turbine in the field with various degrees of success. Finally, effects of scale between laboratory and full-scale testing are considered, along with recommendations for future improvements to dynamic simulations and controller evaluation.

13 HYDRO ENERGY↗
DOE OSTI2015

Precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections

The precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections in a quasistationary state is studied. The rates in question obtained analytically are compared with their numerical counterparts resulting from a dynamical simulation based on the Smoluchowski equation. It is shown that the integral Kramers formula is applicable even in the presence of shell corrections. As the excitation energy of the nucleus being considered decreases, the inaccuracy of this formula decreases monotonically, reaching 2% (inaccuracy of the dynamical simulation) in which case the dimensionless controlling parameter (entropy barrier) becomes greater than four. At small and intermediate values of this parameter, the accuracy of the integral Kramers formula is higher than the accuracy of the approximate Kramers formula obtained from it.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Cavitation and radicals drive the sonochemical synthesis of functional polymer spheres

Sonochemical synthesis can lead to a dramatic increase in the kinetics of formation of polymer spheres (templates for carbon spheres) compared to the modified Stöber silica method applied to produce analogous polymer spheres. Reactive molecular dynamics simulations of the sonochemical process indicate a significantly enhanced rate of polymer sphere formation starting from resorcinol and formaldehyde precursors. The associated chemical reaction kinetics enhancement due to sonication is postulated to arise from the localized lowering of atomic densities, localized heating, and generation of radicals due to cavitation collapse in aqueous systems. This dramatic increase in reaction rates translates into enhanced nucleation and growth of the polymer spheres. The results are of broad significance to understanding mechanisms of sonication induced synthesis as well as technologies utilizing polymers spheres.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A kinetic model of friction in strongly coupled strongly magnetized plasmas

Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge, in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that the “gyro” component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyrofriction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗
DOE OSTI2019

Radiation-induced Swelling and Microcracking in SiC Cladding for LWRs

After the Fukushima accident in 2011, SiC has attracted considerable attention as a potential material for fuel cladding that can offer increased accident tolerance in light water reactors (LWRs) compared to zircaloy cladding. Existing measurements on nuclear grade SiC/SiC composites made of high-purity, stoichiometric SiC have already shown very good radiation tolerance in various harsh environments. However, there are a number of challenges that need to be addressed before SiC composites can be successfully used in LWR technologies. One of them is irradiation-induced swelling, which can lead to development of significant stresses in the matrix, followed by microcracking and subsequent release of fission products. Models of irradiation-induced swelling had been largely empirical and therefore they cannot be used to predict the microstructural dependence of microcracking. Existing models are also limited by the lack of understanding of defects that form due to irradiation in the range of temperatures relevant to fuel cladding in LWRs (<1000°C). Many of the defects in this regime of temperature are too small to be detected with traditional transmission electron microscopy (TEM) techniques. In fact the defects observed in traditional TEM account only for 10-45% of the swelling measured in irradiated SiC. Here, we developed a multi-scale simulation methodology combined with state of the art experimental imaging techniques based on aberration corrected scanning transmission electron microscopy (STEM) to provide a scientific basis for predictive models of radiation swelling and microcracking in SiC. Specifically, we have developed a phase field fracture model that predicts microcrack initiation and growth in polycrystalline SiC. The model accounts for single crystal anisotropy, including anisotropy in the elastic constants and anisotropic surface energy. It also considers the reduced fracture strength of grain boundaries. The model was implemented using the MARMOT mesoscale nuclear materials tool. We have parameterized the model for SiC using molecular dynamics simulation results. The model has been verified using single crystal, bicrystal, and polycrystal examples. Swelling model has also been developed based on cluster dynamics simulations of defect evolution in SiC. We found that point defects and small defect clusters have a significant contribution to swelling. This finding is important because point defects and very small defect clusters are not visible in experiments (e.g., in TEM). Combined modeling and experiments can bring critical insights into how swelling depends on irradiation conditions and the microstructure. In situ nanoindentation tests and ex situ nanoindentation tests were performed on assynthesized 3C-SiC and 3C-SiC irradiated under different radiation conditions. Without radiation, bend contour movement and residual bend contours after indentation demonstrate some room temperature plastic flow and residual plastic deformation. Plastic deformation recovery is also observed on unirradiated samples during unloading. Radiation-induced embrittlement is demonstrated by differences in crack geometry in in situ tests, and is consistent with ex situ test results. Failure by single, straight, cracks indicates decreased fracture toughness due to irradiation in 3C-SiC. Future work includes coupling of the swelling and fracture models and development of a master equation that can be used in the BISON code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Visualizing interfacial collective reaction behaviour of Li–S batteries

Benefiting from high energy density (2,600 Wh kg -1 )) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing the liquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured on the nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li 2 S 2 and Li 2 S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of a dense phase consisting of Li + and S n 2- (2 < n ≤ 6), and the collective charge transfer in the dense phase was further verified by ab initio molecular dynamics simulations. Finally, the collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

25 ENERGY STORAGE↗
DOE OSTI2015

Influence of changes in initial conditions for the simulation of dynamic systems

Chaos theory is a field of study in mathematics, with applications in several disciplines including meteorology, sociology, physics, engineering, economics, biology, and philosophy. Chaos theory studies the behavior of dynamical systems that are highly sensitive to initial conditions—a paradigm popularly referred to as the butterfly effect. Small differences in initial conditions field widely diverging outcomes for such dynamical systems, rendering long-term prediction impossible in general. This happens even though these systems are deterministic, meaning that their future behavior is fully determined by their initial conditions, with no random elements involved. In this paperinfluence of changes in initial conditions will be presented for the simulation of Lorenz system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

Tensile Strength of Carbon Nanotubes Under Realistic Temperature and Strain Rate

Strain rate and temperature dependence of the tensile strength of single-wall carbon nanotubes has been investigated with molecular dynamics simulations. The tensile failure or yield strain is found to be strongly dependent on the temperature and strain rate. A transition state theory based predictive model is developed for the tensile failure of nanotubes. Based on the parameters fitted from high-strain rate and temperature dependent molecular dynamics simulations, the model predicts that a defect free micrometer long single-wall nanotube at 300 K, stretched with a strain rate of 1%/hour, fails at about 9 plus or minus 1% tensile strain. This is in good agreement with recent experimental findings.

Wei, Chen-Yu↗
DOE OSTI2016

Vacancy enhanced formation and phase transition of Cu-rich precipitates in α - iron under neutron irradiation

In this paper, we employed both molecular statics and molecular dynamics simulation methods to investigate the role of vacancies in the formation and phase transition of Cu-rich precipitates in α-iron. The results indicated that vacancies promoted the diffusion of Cu atoms to form Cu-rich precipitates. After Cu-rich precipitates formed, they further trapped vacancies. The supersaturated vacancy concentration in the Cu-rich precipitate induced a shear strain, which triggered the phase transition from bcc to fcc structure by transforming the initial bcc (110) plane into fcc (111) plane. In addition, the formation of the fcc-twin structure and the stacking fault structure in the Cu-rich precipitates was observed in dynamics simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

New insights into the mechanism governing the elasticity of calcium silicate hydrate gels exposed to high temperature: A molecular dynamics study

When exposed to fire, the integrity of cement-based materials is governed by thermally-induced changes in the mechanical properties of their binding phase, i.e., the calcium–silicate–hydrate (C–S–H) gel. However, the effect of temperature on the structure, density, and mechanical properties of C–S–H remains only partially known. Here, based on reactive molecular dynamics simulations, we reveal the nature of thermally-induced damage in C–S–H gels. In general, we show that, at the atomic scale, exposure to high temperature results in partial dehydration, volumetric shrinkage, disordering, and stiffening in the C–S–H grains. However, we show that the thermal response of C–S–H strongly depends on its chemical composition, wherein C–S–H systems associated with lower Ca/Si molar ratios are able to undergo higher temperatures before amorphization. Based on these results, we demonstrate that the stiffness of C–S–H gels (i.e., including porosity—as probed by nanoindentation) is governed by a competition between the stiffening of the grains and the decrease in packing density—wherein the latter eventually become predominant.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. In conclusion, a proposed order parameter distinguishing glass forming ability, 1/(S(T,q 1 )−1, shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .

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

Free Surface Relaxations of Star-Shaped Polymer Films

Here, the surface relaxation dynamics of supported star-shaped polymer thin films are shown to be slower than the bulk, persisting up to temperatures at least 50 degrees above the bulk glass transition temperature T$bulk\atop{g}$. This behavior, exhibited by star-shaped polystyrenes (SPSs) with functionality f = 8 arms and molecular weights per arm M arm < M e (M e is the entanglement molecular weight), is shown by molecular dynamics simulations to be associated with a preferential localization of these macromolecules at the free surface. This new phenomenon is in notable contrast to that of linear chain polymer thin film systems where the surface relaxations are enhanced in relation to the bulk; this enhancement persists only for a limited temperature range above the bulk T$bulk\atop{g}$. Finally, evidence of the slow surface dynamics, compared to the bulk, for temperatures well above T g and at length and time scales not associated with the glass transition has not previously been reported for polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Effects of electronic excitation on cascade dynamics in nickel–iron and nickel–palladium systems

Using molecular dynamics simulations and the two-temperature model, we provide in this paper a comparison of the surviving damage from single ion irradiation events in nickel-based alloys, for cascades with and without taking into account the effects of the electronic excitations. We find that including the electronic effects impacts the amount of the resulting damage and the production of isolated defects. Finally, irradiation of nickel–palladium systems results in larger numbers of defects compared to nickel–iron systems, with similar numbers of isolated defects. We additionally investigate the mass effect on the two-temperature model in molecular dynamics simulations of cascades.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗