Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The SIMRAND methodology: Theory and application for the simulation of research and development projects

A research and development (R&D) project often involves a number of decisions that must be made concerning which subset of systems or tasks are to be undertaken to achieve the goal of the R&D project. To help in this decision making, SIMRAND (SIMulation of Research ANd Development Projects) is a methodology for the selection of the optimal subset of systems or tasks to be undertaken on an R&D project. Using alternative networks, the SIMRAND methodology models the alternative subsets of systems or tasks under consideration. Each path through an alternative network represents one way of satisfying the project goals. Equations are developed that relate the system or task variables to the measure of reference. Uncertainty is incorporated by treating the variables of the equations probabilistically as random variables, with cumulative distribution functions assessed by technical experts. Analytical techniques of probability theory are used to reduce the complexity of the alternative networks. Cardinal utility functions over the measure of preference are assessed for the decision makers. A run of the SIMRAND Computer I Program combines, in a Monte Carlo simulation model, the network structure, the equations, the cumulative distribution functions, and the utility functions.

Miles, R. F., Jr.↗
DOE OSTI2020

Fundamental Phenomena and Applications of Swift Heavy Ion Irradiations

This review concentrates on the specific properties and characteristics of damage structures generated with high-energy ions in the electronic energy loss regime. Irradiation experiments with so-called swift heavy ions (SHIs) find applications in many different fields, with examples presented in ion-track nanotechnology, radiation hardness analysis of functional materials, and laboratory tests of cosmic radiation. The basics of the SHI-solid interaction are described with special attention to processes in the electronic subsystem. The broad spectrum of damage phenomena is exemplified for various materials and material classes, along with a description of typical characterization techniques. This review also presents state-of-the-art modeling efforts that try to account for the complexity of the coupled processes of the electronic and atomic subsystems. Lastly, a brief discussion at the end of this review on SHI phenomena induced by fission fragments in nuclear materials will highlight the relevance of energetic ion irradiation in nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

The structure of CaO–MgO–Al 2 O 3 –SiO 2 melts and glasses doped with FeO X –NiO

Neutron and x-ray diffraction measurements have been performed on CaO–MgO–Al 2 O 3 –SiO 2 (CMAS) glasses doped with NiO–Fe X O at room temperature, along with x-ray measurements on aerodynamically levitated liquids at ≥2000 K. The disordered structures have been modeled using empirical potential structure refinement to investigate the relation between the aluminosilicate network and the modifying cations. The SiO 4 and AlO 4 tetrahedra are found to have wider Si–O and Al–O bond distance distributions in the glass, and the first Ca–O n coordination shell is highly distorted, redistributing different populations of long and short bonds between the liquid and the glass. The addition of Fe and Ni at low aluminosilicate content increases the number of free oxygens not bonded to AlO 4 or SiO 4 . Mg–O and Fe–O are both found to be predominantly fourfold and fivefold in the liquid and glassy states. Despite these low coordination numbers, their bond angle distributions indicate that they are predominantly in nontetrahedral-type geometries, with ferrous and ferric iron possessing similar coordination environments. The Ca–O and Mg–O average coordination numbers and enthalpies of solution are consistent with their higher reactivity within relatively acidic aluminosilicate melts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modelling Development in Radical (Co)Polymerization of Multivinyl Monomers

Radical polymerization (RP) of multivinyl monomers (MVMs) provides a facile solution for manipulating polymer topology and has received increasing attention due to their industrial and academic significance. Continuous efforts have been made to understand their mechanism, which is the key to regulating materials structure. Modelling techniques have become a powerful tool that can provide detailed information on polymerization kinetics which is inaccessible by experiments. Many publications have reported the combination of experiments and modelling for free radical polymerization (FRP) and reversible-deactivation radical polymerizations (RDRP) of MVMs. Herein, a minireview is presented for the most important modelling techniques and their applications in FRP/RDRP of MVMs. This review hopes to illustrate that the combination of modelling and wet experiments can be a great asset to polymer researchers and inspire new thinking for the future MVMs experiment optimization and product design.

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

Scaling and statistics of bottom-up synthesized armchair graphene nanoribbon transistors

Bottom-up assembled nanomaterials and nanostructures allow for the studies of rich and unprecedented quantum-related and mesoscopic transport phenomena. However, it can be difficult to quantify the correlations between the geometrical or structural parameters obtained from advanced microscopy and measured electrical characteristics when they are made into macroscopic devices. Here, we propose a strategy to connect the nanomaterial morphologies and the device performance through a Monte Carlo device model and apply it to understand the scaling trends of bottom-up synthesized armchair graphene nanoribbon (GNR) transistors. A new nanofabrication process is developed for GNR transistors with channel length down to 7 nm. Further, the impacts of the GNR spatial distributions and the device geometries on the device performance are investigated systematically through comparison of experimental data with the model. Through this study, challenges and opportunities of transistor technologies based on bottom-up synthesized GNRs are pinpointed, paving the way to the further improvement of the GNR device performance for future transistor technology nodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sputter deposition of ultrathick Bi coatings onto rotating substrates for inertial confinement fusion

Sputter-deposited Bi is a potential candidate material for the production of inertial confinement fusion (ICF) hohlraums, which are sphero-cylindrical canisters that serve as the housing for hydrogen-fuel-filled capsules in ICF implosions. The fabrication of hohlraums involves deposition onto rotating sphero-cylindrical substrates and, hence, requires an understanding of oblique angle deposition phenomena. Here, we systematically study the effect of dynamic substrate tilt on properties of Bi films deposited onto rotating planar substrates. Results show that film morphology, porosity, and electrical resistivity depend on substrate tilt. In comparison with films deposited onto substrates adhered to a stationary sample holder, films deposited on rotating substrates have more uniform microstructure with fewer void defects. Based on Monte Carlo simulations, these changes in film morphology and microstructure are attributed to averaging and geometric effects related to sample position and rotation. In conclusion, the fabrication of an ∼ 50-µm-thick Bi film on a rotating sphero-cylindrical mandrel demonstrates a significant step toward the production of Bi hohlraums.

Bismuth↗
DOE OSTI2024

Surface Dynamics of Selective Hydration of Rutile TiO 2 : A Kinetic Monte Carlo Approach

The utility of atomic layer deposition (ALD) continues to expand beyond conformal thin-film fabrication to include area- or site-selective ALD. We previously identified a strategy for site-selective ALD (SS-ALD) via the evaluation of thermodynamically preferred hydration of rutile TiO 2 surfaces, as deduced from electronic structure calculations. Here, we report a novel kinetic Monte Carlo (KMC) model that allows for the investigation of surface dynamics and kinetics that improves our understanding of and intuition for the selective hydration strategy. Here we demonstrate the validity of the strategy with respect to step-edge defects for the lowest energy (110) facet as well as report results for the other common facets which agree with experimental STM observations. The results here indicate that the selective hydration strategy is feasible both thermodynamically (evaluated in our previous publication) and kinetically (from the KMC model). Because diffusion has a slower rate than others, we find that any proximity effects between terrace and defect sites are unlikely to affect the selective hydration strategy for rutile TiO 2 . The KMC model further provides relevant timescales for achieving selectivity experimentally and establishes the kinetic viability of the selective hydration approach to SS-ALD.

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

Simultaneous Development of Antiferromagnetism and Local Symmetry Breaking in a Kagome Magnet (Co 0.45 Fe 0.55 )Sn

CoSn and FeSn, two kagome-lattice metals, have recently attracted significant attention as hosts of electronic flat bands and emergent physical properties. However, current understandings of their physical properties are limited to knowledge of the average crystal structure. Here, we report the Fe-doping induced coemergence of the antiferromagentic (AFM) order and local symmetry breaking in (Co 0.45 Fe 0.55 )Sn. Rietveld analysis on the neutron and synchrotron X-ray diffraction data indicates A-type antiferromagnetic order with the moment pointing perpendicular to the kagome layers, associated with the anomaly in the MSn(1) 2 Sn(2) 4 (M = Co/Fe) octahedral distortion and the lattice constant c. Reverse Monte Carlo (RMC) modeling of the synchrotron X-ray total scattering results captured the subtle local orthorhombic distortion involving off-axis displacements of Sn(2). Our results indicate that the stable hexagonal lattice above T N becomes unstable once the A-type AFM order is formed below T N . Here we argue that the local symmetry breaking has a magnetic origin, since the spatially varied M–Sn(2) bond lengths arise from out-of-plane magnetic exchange coupling J c via the exchange pathway M–Sn(2)–M. Our study provides comprehensive information on the crystal structure in both long-range scale and local scale, unveiling unique coupling between AFM order, octahedral distortion, and hidden local symmetry breaking.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

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

Kinetic simulations of ignited mode cesium vapor thermionic converters

Cesium vapor thermionic converters are an attractive method of converting high-temperature heat directly to electricity, but theoretical descriptions of the systems have been difficult due to the multi-step ionization of Cs through inelastic electron–neutral collisions. This work presents particle-in-cell simulations of these converters, using a direct simulation Monte Carlo collision model to track 52 excited states of Cs. Here, these simulations show the dominant role of multi-step ionization, which also varies significantly based on both the applied voltage bias and pressure. The electron energy distribution functions are shown to be highly non-Maxwellian in the cases analyzed here. A comparison with previous approaches is presented, and large differences are found in ionization rates due especially to the fact that previous approaches have assumed Maxwellian electron distributions. Finally, an open question regarding the nature of the plasma sheaths in the obstructed regime is discussed. The one-dimensional simulations did not produce stable obstructed regime operation and thereby do not support the double-sheath hypothesis.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2021

Ion heating in the PISCES-RF liquid-cooled high-power, steady-state, helicon plasma device

Radio frequency (RF) driven helicon plasma sources are commonly used for their ability to produce high-density argon plasmas ( n > 10 19 m -3 ) at relatively moderate powers (typical RF power < 2 kW). Typical electron temperatures are <10 eV and typical ion temperatures are <0.6 eV. A newly designed helicon antenna assembly (with concentric, double-layered, fully liquid-cooled RF-transparent windows) operates in steady-state at RF powers up to 10 kW. We report on the dependence of argon plasma density, electron temperature and ion temperature on RF power. At 10 kW, ion temperatures >2 eV in argon plasmas are measured with laser induced fluorescence, which is consistent with a simple volume averaged 0D power balance model. 1D Monte Carlo simulations of the neutral density profile for these plasma conditions show strong neutral depletion near the core and predict neutral temperatures well above room temperatures. The plasmas created in this high-power helicon source (when light ions are employed) are ideally suited for fusion divertor plasma-material interaction studies and negative ion production for neutral beams.

Physics↗
DOE OSTI2025

External Magnetic Field Suppression of Carbon Diffusion in Iron

External magnetic fields reduce diffusion of carbon in BCC iron, but the physical mechanism is not understood. Using DFT calculations with magnetic moments sampled from a Heisenberg model, we calculate diffusivities of carbon in iron at high temperatures and with field. Our model reproduces the measured suppression of diffusivity from field. Here, we find that increasing magnetic disorder flattens the electron density of states compared with the ferromagnetic case, which distorts the octahedral cages around carbon, lowering the activation barrier to diffusion; an applied field reverses these trends.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗