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 253 records · Page 14

Multi-dimensional modelling of electrostatic force distance curve over dielectric surface: Influence of tip geometry and correlation with experiment

Electric Force-Distance Curves (EFDC) is one of the ways whereby electrical charges trapped at the surface of dielectric materials can be probed. To reach a quantitative analysis of stored charge quantities, measurements using an Atomic Force Microscope (AFM) must go with an appropriate simulation of electrostatic forces at play in the method. This is the objective of this work, where simulation results for the electrostatic force between an AFM sensor and the dielectric surface are presented for different bias voltages on the tip. The aim is to analyse force-distance curves modification induced by electrostatic charges. The sensor is composed by a cantilever supporting a pyramidal tip terminated by a spherical apex. The contribution to force from cantilever is neglected here. A model of force curve has been developed using the Finite Volume Method. The scheme is based on the Polynomial Reconstruction Operator—PRO-scheme. First results of the computation of electrostatic force for different tip–sample distances (from 0 to 600 nm) and for different DC voltages applied to the tip (6 to 20 V) are shown and compared with experimental data in order to validate our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time evolution of electron waves in graphene superlattices

The time evolution of electron waves in graphene superlattices is studied using both microscopic and “effective medium” formalisms. The numerical simulations reveal that in a wide range of physical scenarios it is possible to neglect the granularity of the superlattice and characterize the electron transport using a simple effective Hamiltonian. It is verified that as general rule the continuum approximation is rather accurate when the initial state is less localized than the characteristic spatial period of the superlattice. This property holds even when the microsocopic electric potential has a strong spatial modulation or in presence of interfaces between different superlattices. Detailed examples are given both of the time evolution of initial electronic states and of the propagation of stationary states in the context of wave scattering. The theory also confirms that electrons propagating in tailored graphene superlattices with extreme anisotropy experience virtually no diffraction.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Atomistic simulation of frictional anisotropy on quasicrystal approximant surfaces

J. Y. Park et al. [Science 309, 1354 (2005)] have reported eight times greater atomic-scale friction in the periodic than in the quasiperiodic direction on the twofold face of a decagonal Al-Ni-Co quasicrystal. Here we present results of molecular-dynamics simulations intended to elucidate mechanisms behind this giant frictional anisotropy. Simulations of a bare atomic-force-microscope tip on several model substrates and under a variety of conditions failed to reproduce experimental results. On the other hand, including the experimental passivation of the tip with chains of hexadecane thiol, we reproduce qualitatively the experimental anisotropy in friction, finding evidence for entrainment of the organic chains in surface furrows parallel to the periodic direction.

36 MATERIALS SCIENCE↗

Effect of directed aging on nonlinear elasticity and memory formation in a material

Disordered solids often change their elastic response as they slowly age. Using experiments and simulations, we study how aging disordered planar networks under an applied stress affects their nonlinear elastic response. We are able to modify dramatically the elastic properties of our systems in the nonlinear regime. Using simulations, we study two models for the microscopic evolution of properties of such a material; the first considers changes in the material strength, while the second considers distortions in the microscopic geometry. Both models capture different aspects of the experiments including the encoding of memories of the aging history of the system and the dramatic effects on the material’s nonlinear elastic properties. Our results demonstrate how aging can be used to create complex elastic behavior in the nonlinear regime.

36 MATERIALS SCIENCE↗

Quantum Non-Equilibrium Dynamics Prediction of Electronic Transport Coefficients in Non-Linear Regimes

The goal of this work was to develop unique LLNL capabilities for the computational prediction of electrical, optical and thermal transport properties from quantum simulations. Using a recently developed theory of microscopic quantum transport, we were able to model compressed metallic hydrogen and it’s non-linear and out-of-equilibrium contributions in electrical and thermal transport. This theory was implemented in the open-source, real-space time-dependent density functional theory code Octopus. This model is general and widely applicable to high energy density systems that are important to LLNL and NNSA mission areas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of New Experimental Methods for Correlated Operando Surface/Gas Characterization

The predictive understanding of catalytic surface reactions requires accurate microkinetic models, and while decades of work has been devoted to the elucidation of the reaction steps in these models, many open questions remain. One key issue is a lack of approaches enabling the local spatially resolved assessment of catalytic activity over a surface. In this report, we detail efforts to develop a new diagnostic approach to solve this problem. The approach is based upon laser resonance enhanced multiphoton ionization of reaction products emitted into the gas phase followed by spatially resolved imaging of the resultant ions or electrons. Ion imaging is pursued with a velocity-selected spatially resolved ion imaging microscope, while electron imaging was attempted in a low energy electron microscope. Successful demonstration of the ion imaging microscope coupled with the development of transport simulations shows promise for a revolutionary new tool to assess local catalytic activity

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiconductor‐Electrocatalyst Contacts: Theory, Experiment, and Applications to Solar Water Photoelectrolysis

The goal of the work supported by this program was to use new types of photoelectrochemical experiments, supported by simulation and theory, to develop an improved microscopic picture of electron transfer in combined semiconductor|catalyst|solution (sem|cat|sol) systems. An emphasis was placed on developing techniques which can directly characterize electronic interfacial phenomena without first assuming a junction model. We then aimed to illustrate the design principles that enhance interfacial processes and enable high-efficiency devices. Our work thus provides deep insight into the processes that, in general, use light to drive chemical reactions.

14 SOLAR ENERGY↗

WE-D-BRF-01: FEATURED PRESENTATION - Investigating Particle Track Structures Using Fluorescent Nuclear Track Detectors and Monte Carlo Simulations

Purpose: To report on the efforts funded by the AAPM seed funding grant to develop the basis for fluorescent nuclear track detector (FNTD) based radiobiological experiments in combination with dedicated Monte Carlo simulations (MCS) on the nanometer scale. Methods: Two confocal microscopes were utilized in this study. Two FNTD samples were used to find the optimal microscope settings, one FNTD irradiated with 11.1 MeV/u Gold ions and one irradiated with 428.77 MeV/u Carbon ions. The first sample provided a brightly luminescent central track while the latter is used to test the capabilities to observe secondary electrons. MCS were performed using TOPAS beta9 version, layered on top of Geant4.9.6p02. Two sets of simulations were performed, one with the Geant4-DNA physics list and approximating the FNTDs by water, a second set using the Penelope physics list in a water-approximated FNTD and a aluminum-oxide FNTD. Results: Within the first half of the funding period, we have successfully established readout capabilities of FNTDs at our institute. Due to technical limitations, our microscope setup is significantly different from the approach implemented at the DKFZ, Germany. However, we can clearly reconstruct Carbon tracks in 3D with electron track resolution of 200 nm. A second microscope with superior readout capabilities will be tested in the second half of the funding period, we expect an improvement in signal to background ratio with the same the resolution.We have successfully simulated tracks in FNTDs. The more accurate Geant4-DNA track simulations can be used to reconstruct the track energy from the size and brightness of the observed tracks. Conclusion: We have achieved the goals set in the seed funding proposal: the setup of FNTD readout and simulation capabilities. We will work on improving the readout resolution to validate our MCS track structures down to the nanometer scales.

60 APPLIED LIFE SCIENCES↗

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗

Ptychographic atomic electron tomography: Towards three-dimensional imaging of individual light atoms in materials

Through numerical simulations, we demonstrate the combination of ptychography and atomic electron tomography as an effective method for low dose imaging of individual low-Z atoms in three dimensions. After generating noisy diffraction patterns with multislice simulations of an aberration-corrected scanning transmission electron microscope through a 5-nm zinc-oxide nanoparticle, we have achieved three-dimensional (3D) imaging of individual zinc and oxygen atoms and their defects by performing tomography on ptychographic projections. The methodology has also been simulated in 2D materials, resolving individual sulfur atoms in vertical WS 2 /WSe 2 van der Waals heterostructure with a low total electron dose where annular-dark-field images fail to resolve. Overall, we envision that the development of this method could be instrumental in studying the precise 3D atomic structures of radiation sensitive systems and low-Z atomic structures such as 2D heterostructures, catalysts, functional oxides, and glasses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Final Technical Report "Multiscale Simulation Algorithms for Biochemical Systems"

Biochemical systems are inherently multiscale and stochastic. In microscopic systems formed by living cells, the small numbers of reactant molecules can result in dynamical behavior that is discrete and stochastic rather than continuous and deterministic. An analysis tool that respects these dynamical characteristics is the stochastic simulation algorithm (SSA, Gillespie, 1976), a numerical simulation procedure that is essentially exact for chemical systems that are spatially homogeneous or well stirred. Despite recent improvements, as a procedure that simulates every reaction event, the SSA is necessarily inefficient for most realistic problems. There are two main reasons for this, both arising from the multiscale nature of the underlying problem: (1) stiffness, i.e. the presence of multiple timescales, the fastest of which are stable; and (2) the need to include in the simulation both species that are present in relatively small quantities and should be modeled by a discrete stochastic process, and species that are present in larger quantities and are more efficiently modeled by a deterministic differential equation (or at some scale in between). This project has focused on the development of fast and adaptive algorithms, and the fun- damental theory upon which they must be based, for the multiscale simulation of biochemical systems. Areas addressed by this project include: (1) Theoretical and practical foundations for ac- celerated discrete stochastic simulation (tau-leaping); (2) Dealing with stiffness (fast reactions) in an efficient and well-justified manner in discrete stochastic simulation; (3) Development of adaptive multiscale algorithms for spatially homogeneous discrete stochastic simulation; (4) Development of high-performance SSA algorithms.

97 MATHEMATICS AND COMPUTING↗

Dimensionality Reduction Algorithms for Big Nuclear Data Streams with Application to SCALE Cross Sections

A central problem in performing accurate nuclear reactor simulations is the ability to capture the changes in the few group microscopic cross sections as functions of system state, e.g. moderator density, fuel temperature, or neutron spectrum, and time and/or burnup of the fuel. For most reactor simulation scenarios, this results in huge data stream, which must be interpolated as needed by downstream calculations (e.g., core-wide simulation). This is typically accomplished by using sophisticated interpolating functions to capture the few-group variations over the wide range of conditions expected during operation. To address handling this big data, this paper explores the use of dimensionality reduction (DR) techniques to identify the true dimensionality of the data prior to its use in the calculations. This could significantly reduce the computational cost required for data processing/manipulation, thus reducing a significant portion of the total calculations time. By way of an example, in the current version of the ORIGEN-ARP module, the assembly-averaged microscopic cross sections at different conditions, e.g., fuel enrichments, irradiation history and water density, which are calculated using the TRITON sequence and interpolated on-demand by the ARP module to model assembly depletion for a wide range of assembly types (e.g. Westinghouse 17x17 and GE 9x9) over a wide range of conditions. This approach typically results in ∼1 GB per assembly type, which although not insurmountable, can be impractical for distribution. and it is the premise of this paper that most of this data is heavily correlated and hence can take great benefit of condensation via DR. To demonstrate DR, the proposed algorithm is applied to the nuclear data streams generated by SCALE 6.1. Using the depletion sequence (TRITON), microscopic cross sections are calculated for a mixed oxide (MOX) fuel assembly design at a representative range of fuel enrichments and fuel burnup. In conventional core simulators, multi dimensional interpolation of these cross sections is used to provides neutronics feedback due to local conditions without returning to expensive TRITON calculations for every change. This paper shows that dimensionality reduction can help compress the size of the cross section data without sacrificing the accuracy of their reconstruction at all conditions of interest. This work explores the application of dimensionality reduction algorithms to nuclear data streams typically employed in support of conventional full-core reactor calculations which rely on interpolation of pre-calculated microscopic cross sections. In particular, analysis of the few group MOX cross sections revealed that a large degree of dependencies exist between the few-group cross sections, which can be revealed via range-finding algorithms. Specifically, correlations between different reaction types for the same isotope have been identified, and correlations between different enrichments and correlations across different isotopes have been identified. This paper demonstrates that the reduction errors can be quantified, and they are uncorrelated with the range of variations for the various cross sections. In future work, we will focus on identifying the smallest rank that takes into account all types of correlations across the range of conditions used to model cross section variations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Geometric state function for two-fluid flow in porous media

Models that describe two-fluid flow in porous media suffer from a widely recognized problem that the constitutive relationships used to predict capillary pressure as a function of the fluid saturation are nonunique, thus requiring a hysteretic description. As an alternative to the traditional perspective, we consider a geometric description of the capillary pressure, which relates the average mean curvature, the fluid saturation, the interfacial area between fluids, and the Euler characteristic. The state equation is formulated using notions from algebraic topology and cast in terms of measures of the macroscale state. Synchrotron-based x-ray microcomputed tomography and high-resolution pore-scale simulation is applied to examine the uniqueness of the proposed relationship for six different porous media. We show that the geometric state function is able to characterize the microscopic fluid configurations that result from a wide range of simulated flow conditions in an averaged sense. The geometric state function can serve as a closure relationship within macroscale models to effectively remove hysteretic behavior attributed to the arrangement of fluids within a porous medium. As a result, this provides a critical missing component needed to enable a new generation of higher fidelity models to describe two-fluid flow in porous media.

42 ENGINEERING↗

Soft collective fluctuations governing hydrophobic association

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. Furthermore, the microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.

74 ATOMIC AND MOLECULAR PHYSICS↗

Performance Prediction of Irradiated TREAT Fuel

Poster Title: Performance Prediction of Irradiated TREAT Fuel The poster presents findings about the irradiation damage layer effects on the microscopic fuel temperature profile in TREAT. Microscopic fuel particles in TREAT are small enough that fission fragments travel out of the particles causing thermal property reduction in a shell around the fuel particle. The use of Magpie and MAMMOTH code is illustrated. The thermal temperature solution is simulated using EOS pulse power data. The reduction of microscopic fuel particle size reduced the irradiation damage due to fission fragments. The reduction of fission events is shown to be greater than the increase in fission fragments leaving the particle as size decreases. The temperature difference between the microscopic fuel particle and bulk graphite is shown to increase with irradiation as the damage layer grows.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Global System for Transportation Simulation and Visualization in Emergency Evacuation Scenarios

Simulation-based studies are frequently used for evacuation planning and decision making processes. Given the transportation systems complexity and data availability, most evacuation simulation models focus on certain geographic areas. With routine improvement of OpenStreetMap road networks and LandScanTM global population distribution data, we present WWEE, a uniform system for world-wide emergency evacuation simulations. WWEE uses unified data structure for simulation inputs. It also integrates a super-node trip distribution model as the default simulation parameter to improve the system computational performance. Two levels of visualization tools are implemented for evacuation performance analysis, including link-based macroscopic visualization and vehicle-based microscopic visualization. For left-hand and right-hand traffic patterns in different countries, the authors propose a mirror technique to experiment with both scenarios without significantly changing traffic simulation models. Ten cities in US, Europe, Middle East, and Asia are modeled for demonstration. With default traffic simulation models for fast and easy-to-use evacuation estimation and visualization, WWEE also retains the capability of interactive operation for users to adopt customized traffic simulation models. For the first time, WWEE provides a unified platform for global evacuation researchers to estimate and visualize their strategies performance of transportation systems under evacuation scenarios.

Emergency Evacuation↗