Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microscopic 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Berni Julian Alder, theoretical physicist and inventor of molecular dynamics, 1925–2020

Berni Julian Alder, one of the leading figures in the invention of molecular dynamics simulations used for a wide array of problems in physics and chemistry, died on September 7th, 2020. His career, spanning more than 65 years, transformed statistical mechanics, many body physics, the study of chemistry and the microscopic dynamics of fluids, by making atomistic computational simulation (in parallel with traditional theory and experiment) a new pathway to unexpected discoveries. Among his many honors, the CECAM prize, recognizing exceptional contributions to the simulation of the microscopic properties of matter is named for him. He was awarded the National Medal of Science by President Obama in 2008.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Reactor simulations with nuclear data uncertainties

This work demonstrates the influence of uncertainties in microscopic nuclear data on the results of reactor simulations, both for stationary and transient states. It gives an overview on the methods in use for uncertainty and sensitivity analyses with respect to nuclear data, and discusses the pros and cons. For full-scale reactor simulations, in particular with coupled neutron transport and thermo-hydraulics, random sampling provides a powerful means to propagate nuclear data uncertainties through the complete calculation sequence. Results of uncertainty analyses performed with the GRS XSUSA – “cross section (XS) Uncertainty and Sensitivity Analysis” methodology are shown for radial power distributions from steady-state PWR calculations and for the time evolution of the reactor power in the course of a control rod withdrawal from a PWR mini-core. In all cases, the output uncertainties are considerable. For the radial power distributions, relative 1σ uncertainties of up to 10% are observed, and for the power peak during the transient, the relative 1σ uncertainty reaches 20%. These large uncertainties strongly suggest to routinely accompany best-estimate simulations by uncertainty analyses with respect to nuclear data, in particular for systems beyond LWR for which much less operation experience is available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2013

Development of the Code RITRACKS

A document discusses the code RITRACKS (Relativistic Ion Tracks), which was developed to simulate heavy ion track structure at the microscopic and nanoscopic scales. It is a Monte-Carlo code that simulates the production of radiolytic species in water, event-by-event, and which may be used to simulate tracks and also to calculate dose in targets and voxels of different sizes. The dose deposited by the radiation can be calculated in nanovolumes (voxels). RITRACKS allows simulation of radiation tracks without the need of extensive knowledge of computer programming or Monte-Carlo simulations. It is installed as a regular application on Windows systems. The main input parameters entered by the user are the type and energy of the ion, the length and size of the irradiated volume, the number of ions impacting the volume, and the number of histories. The simulation can be started after the input parameters are entered in the GUI. The number of each kind of interactions for each track is shown in the result details window. The tracks can be visualized in 3D after the simulation is complete. It is also possible to see the time evolution of the tracks and zoom on specific parts of the tracks. The software RITRACKS can be very useful for radiation scientists to investigate various problems in the fields of radiation physics, radiation chemistry, and radiation biology. For example, it can be used to simulate electron ejection experiments (radiation physics).

Plante, Ianik↗
DOE OSTI2018

Radiocesium interaction with clay minerals: Theory and simulation advances Post–Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

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

Computational simulations of solvation force and squeezing out of dodecane chain molecules in an atomic force microscope

Understanding the squeeze out behaviors of liquid films at nanometer scale in an atomic force microscope (AFM) has been a significant interest since the 1990s. We carry out all-atom static-mode AFM simulations in a liquid-vapor molecular dynamics ensemble to investigate the solvation force oscillation and squeeze out mechanisms of a confined linear dodecane fluid between a gold AFM tip and a mica substrate. Solvation force oscillations are found to be associated with the layering transition of the liquid film and unstable jumps of the AFM tip. Detailed structural analyses and molecular animations show that the local permeation of chain molecules and the squeeze out of molecules near the edge of contact promote the layering transition under compression. The confinement-induced slow down dynamics is manifested by the decrease in diffusivity and increase in rotational relaxation times. Furthermore, the persistent diffusive behavior of dodecane chain molecules even in the single-monolayer film is attributed to the chain sliding motions in the film due to the substantial vacancy space and thermal fluctuations.

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

Start-to-end self-consistent simulation of the UEM beamline

Addressing the preliminary design review report of the ultrafast electron microscope (UEM), start-to-end self-consistent simulations were carried out to explore fundamental aspects affecting the resolution (space-charge effects, energy spread, shot-to-shot energy jitter, contrast formation, etc.). The electron beam energy spread in the single-shot mode and the effective energy spread in the accumulation mode, equivalent to the convolution of the single-shot energy spread and shot-to-shot energy jitter, have been identified as the ultimate limits of the instrument resolution. A full self-consistent simulation from the photocathode to the detector is being performed using the GPT code to find out what resolution can be achieved. Preliminary results of the start-to-end simulations are reported in this note and the plan of further simulations is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

The Role of Organic Aerosol in Atmospheric Ice Nucleation: A Review

Predicting the formation of ice crystals in the atmosphere presents one of the great challenges in physical sciences with important implications for the chemistry and composition of the Earth’s atmosphere, the hydrological cycle, and climate. Among atmospheric ice formation processes, heterogeneous ice nucleation proceeds on aerosol particles ranging from a few nanometers to micrometers in size, commonly referred to as ice nucleating particles (INPs). Research over the last two decades has demonstrated that organic matter (OM) is ubiquitous in the atmosphere, present as organic aerosol (OA) particles or as coatings on other particle types. While it is well understood that mineral dust particles can induce ice nucleation under atmospheric conditions, it is currently not clear how much OA particles including pollution aerosol and secondary organic aerosol (SOA) contributes to the ambient INP population. The physicochemical properties of OM make predicting how OM interacts with water vapor leading ultimately to ice crystals challenging. This review focuses on the role of OM in INPs, summarizing and highlighting recent advances in our understanding of the ice nucleation process gained from theoretical, laboratory, and field studies. Examination of ice residuals and INPs with analytical techniques demonstrates that OM participates in atmospheric ice crystal formation. Molecular dynamic simulations provide insight into the microscopic processes that initiate ice nucleation. The amorphous phase state of OM in the supercooled and metastable regime is identified as a key factor in assessing the particles’ nucleation pathways and rates. A theoretical model is advanced, based on particle water activity, to holistically predict amorphous phase changes and ice nucleation rates of particles coated by OM. The goal of this review is to synthesize our current understanding and propose future research directions needed to fully evaluate how OA particles contribute to INPs in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Multi-pass transmission electron microscopy

Feynman once asked physicists to build better electron microscopes to be able to watch biology at work. While electron microscopes can now provide atomic resolution, electron beam induced specimen damage precludes high resolution imaging of sensitive materials, such as single proteins or polymers. Here, we use simulations to show that an electron microscope based on a multi-pass measurement protocol enables imaging of single proteins, without averaging structures over multiple images. While we demonstrate the method for particular imaging targets, the approach is broadly applicable and is expected to improve resolution and sensitivity for a range of electron microscopy imaging modalities, including, for example, scanning and spectroscopic techniques. The approach implements a quantum mechanically optimal strategy which under idealized conditions can be considered interaction-free.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Study of Ferrite During Refinement of Prior Austenite Grains in Microalloyed Steel Continuous Casting

The formation of coarse prior austenite grain is a key factor to promote transverse crack, and the susceptibility to the transverse crack can be reduced by refining the austenite grain size. In the present study, the high-temperature confocal laser scanning microscope (CLSM) was used to simulate two types of double phase-transformation technologies. The distribution and morphology of ferrites under different cooling conditions were analyzed, and the effects of ferrite distribution and morphology on the double phase-transformation technologies were explored to obtain the suitable double phase-change technology for the continuous casting process. The results indicate that, under the thermal cycle TH0 [the specimens were cooled down to 913 K (640 °C) at a cooling rate of 5.0 K/s (5.0 °C/s)], the width of prior austenite grain boundaries was thick, and the dislocation density at grain boundaries was high. It had strong inhibition effect on crack propagation; under the thermal cycle TH1 [the specimens were cooled down to 1073 K (800 °C) at a cooling rate of 5.0 K/s (5.0 °C/s) and then to 913 K (640 °C) at a cooling rate of 1.0 K/s (1.0 °C/s)], the width of prior austenite grain boundary was thin, and the dislocation density at grain boundaries was low. It was beneficial to crack propagation. After the first phase change, the developed film-like ferrite along the austenite grain boundaries improved the nucleation conditions of new austenitic grains and removed the inhibition effect of the prior austenite grain boundaries on the austenite grain size.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Research on the Corrosion Behavior of X70 Pipeline Steel Under Coupling Effect of AC + DC and Stress

The coupling effect of stray AC (alternating current) + DC (direct current) and tensile stress-induced X70 pipeline corrosion is complicated with its mechanism uncovered. In this work, the corrosion behavior of X70 pipeline steel under three different corrosion (including an individual AC corrosion, mixed AC + 150 MPa tensile stress corrosion, coupling effect of AC + a 30 A/m{sup 2} DC interference and 150 MPa tensile stress corrosion) conditions in Dagang soil simulated solution was investigated by COMSOL Multiphysics software simulation, electrochemical tests, scanning electron microscope and other techniques. The results indicated that with the increase in AC density, the open-circuit potential of the X70 steel samples was negatively shifted and corrosion rate increased in three different corrosion conditions. The corrosion rates of X70 steel samples under individual AC interference were almost unchanged when 150.9 MPa tensile stress was applied. When applying 30 A/m{sup 2} DC interference, the corrosion rate of X70 steel increases significantly. Therefore, the factors affecting corrosion performance of X70 pipeline steel were DC > AC > stress. Under individual AC interference, uniform corrosion occurred on the X70 steel at an AC current density of less than 30 A/m{sup 2}, while the pitting corrosion occurred at an AC current density of greater than 100 A/m{sup 2}. The application of tensile stress and DC had a certain inhibitory effect on pitting corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mechanical properties of metal-ceramic nanolaminates: Effect of constraint and temperature

Al/SiC nanolaminates with equal nominal thicknesses of the Al and SiC layers (10, 25, 50 and 100 nm) were manufactured by magnetron sputtering. The mechanical properties were measured at 25 °C and 100 °C by means of nanoindentation and micropillar compression tests and the deformation mechanisms were analyzed by in situ micropillar compression tests in the transmission electron microscope. In addition, finite element simulations of both tests were carried out to ascertain the role played by the strength of the Al layers and by the elastic constraint of the ceramic layers on the plastic flow of Al in the mechanical response. It was found that the mechanical response was mainly controlled by the constraint during nanoindentation or micropillar compression tests of very thin layered (≈10 nm) laminates, while the influence of the strength of Al layers was not as critical. This behavior was reversed, however, for thick layered laminates (100 nm). Here, these mechanisms point to the different effects of layer thickness during nanoindentation and micropillar compression, at both temperatures, and showed the critical role played by constraint on the mechanical response of nanolaminates made of materials with a very large difference in the elasto-plastic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Multi-scale computational analysis of unidirectional carbon fiber reinforced polymer composites under various loading conditions

In this study, a multi-scale computational analysis based on representative volume element (RVE) modeling and molecular dynamics (MD) simulations is developed to investigate the microscopic failure mechanisms of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The average properties of the 200-nm thickness interphase region between fiber and matrix are characterized through MD simulations and an analytical gradient model. The results demonstrate that the interphase region has higher Young’s modulus and strength, compared to the bulk matrix. This stiffened interphase region influences the composite response significantly. Specifically, the traditional two-phase model with zero-thickness interface fails to capture the stress-strain behavior compared to the experimental data. However, by adding the interphase region to a modified RVE model, the accuracy of simulation results will be improved significantly. Furthermore, a coupled experimental-computational micromechanics approach is adopted to calibrate and validate the cohesive parameters of the interface. By including the cohesive interface, our modified RVE model accurately captures the failure strength of the composites. Finally, different failure mechanisms for specimens are investigated using our multi-scale computational framework. In conclusion, the results show that the failure modes of UD CFRP composites are very complex and multiple failure mechanisms co-exist depending on the loading conditions, agreeing well with our experimental analyses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mobility of a spatially modulated electron liquid on the helium surface

Here, we present the results of large-scale numerical simulations of the mobility of a two-dimensional electron liquid on the helium surface in the presence of a one-dimensional periodic potential. Even where the potential is much weaker than the electron-electron interaction, it can strongly change the mobility. The effect depends on the interrelation between the potential period and the mean interelectron distance. It is most pronounced where the period is close to the period of the Wigner crystal that would form if the liquid were cooled to a lower temperature. The results suggest, in particular, that the correlation length in the electron liquid can be found by measuring the mobility in a weak periodic potential. The simulations are based on the microscopic model of the electron scattering by the excitations in helium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Benchmarking Quantum Simulators Using Ergodic Quantum Dynamics

Here we propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced spatiotemporal control. Our protocol relies on universal fluctuations emerging from generic Hamiltonian dynamics, which we discover in the present work. It does not require fine-tuned control over state preparation, quantum evolution, or readout capability, while achieving near optimal sample complexity: a percent-level precision is obtained with ~103 measurements, independent of system size. Furthermore, the accuracy of our fidelity estimation improves exponentially with increasing system size. We numerically demonstrate our protocol in a variety of quantum simulator platforms, including quantum gas microscopes, trapped ions, and Rydberg atom arrays. We discuss applications of our method for tasks such as multiparameter estimation of quantum states and processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Directional dependency of electronic stopping in nickel, projectile’s excited charge state and momentum transfer

Abstract We studied the directional dependency of electronic stopping power of swift light ions in nickel using real-time time-dependent density functional theory. We report a variation of electronic stopping for moving ions as the projectile probes different electronic densities of the host material. These results show that while the predicted magnitude stays in reasonable agreement with experiment, for $$v > 2$$ v > 2 . a.u. simulating only low index crystallographic directions is not enough to sample the experimental average values. The ab initio simulations give us access to microscopic quantities, such as non-adiabatic forces, momentum transfer and transient excited state charges of the projectile and host ions, which are not available through other methods. We report these quantities for the first time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗