Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Transport Simulations”

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 55 records · Page 3

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monte Carlo simulations of particle acceleration at oblique shocks

The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.

Baring, Matthew G.↗

Relativistic electron transport and bremsstrahlung production in solar flares

A Monte Carlo simulation of ultrarelativistic electron transport in solar flare magnetic loops has been developed. It includes Coulomb, synchrotron, and bremsstrahlung energy losses; pitch-angle scattering by Alfven and whistler turbulence in the coronal region of the loop; and magnetic mirroring in the converging magnetic flux tubes beneath the transition region. Depth distributions, time profiles, energy spectra, and angular distributions of the resulting bremsstrahlung emission are calculated. It is found that both the preferential detection of solar flares with greater than 10 MeV emission near the limb of the sun and the observation of ultrarelativistic electron bremsstrahlung from flares on the disk are consequences of the loop transport model. The declining portions of the observed time profiles of greater than 10 MeV emission from solar flares can also be accounted for, and it is proposed that these portions are determined by transport and not acceleration.

Miller, James A.↗

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗

COG User's Manual: A Multiparticle Monte Carlo Transport Code (Sixth Edition)

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary 3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to make it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems while retains all of the standard capabilities found in other Monte Carlo transport codes. COG uses high-resolution pointwise cross-section databases and makes no compromises in the transport physics, so that the results of a COG run are limited only by the accuracy of the databases used. COG runs primarily on Linux Operating System workstations with MPICH software installed – currently, Red Hat 7 & 8, Windows 10 (Windows Subsystem for Linux –WSL), Ubuntu 16, 18 & 20, OpenSUSE Leap 15.2, Fedora 32, Apple Power Mac with Intel CPU (with MacPorts installed) workstations, and LLNL LC supercomputer CTS-1 cluster with TOSS 3 are supported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

General relativistic simulations of collapsing binary neutron star mergers with Monte Carlo neutrino transport

Recent gravitational wave observations of neutron-star-neutron-star and neutron-star-black-hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. These discoveries have led to an increased interest in the simulation of merging compact binaries involving massive stars. In this paper, we present a first set of evolution of massive neutron star binaries using Monte Carlo radiation transport for the evolution of neutrinos. Here, we study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting postmerger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the whiskyTHC code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the postmerger remnant. One of the main potential uses of our simulations is, thus, as improved initial conditions for longer evolutions of such remnants.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion and relativistic electron transport in solar flares

Recent models for the ion and relativistic electron transport in solar flare magnetic loops are reviewed, focusing on the ways in which the models handle loop structure, particle acceleration, and particle transport. Results are presented from Monte Carlo simulations of ion and relativistic ion transport. Consideration is given to results concerning the depth distributions of gamma-ray production, the attenuation of gamma-ray lines, the time dependences of the various emissions and the angular distribution of the bremsstrahlung.

Ramaty, R.↗

Simulation of Particle-Material Interactions

This paper gives an overview of the particle transport theory essentials, the basics of particle-material interaction simulation, physical quantities needed to simulate particle transport and interactions in materials, Monte Carlo simulation flow, response of additive detectors, statistical weights and other techniques to minimize statistical errors. Effects in materials under irradiation, materials response related to component lifetime and performance are considered with a focus on high-energy and high-power accelerator applications. Implementation of simulation of particle-material interactions in the modern Monte Carlo codes along with the code s main features and results of recent benchmarking are described.

43 PARTICLE ACCELERATORS↗

A Computationally Efficient Algorithm for Sampling the Rudd Differential Cross Section

Monte Carlo radiation transport codes such as RITRACKS or Geant4 are used to simulate the interaction of ions with matter. These codes rely on sampling algorithms to determine interactions and various physical properties of particles involved in the simulations. It is crucial to develop efficient sampling algorithms since Monte Carlo radiation transport simulations can be time consuming. This work presents an efficient sampling algorithm to determine the energy of secondary electrons following ion-water interactions. The applicability and intended use of the algorithm are discussed in detail, and it is shown that the new algorithm is up to 6X10 4 times faster than the method currently used in Geant4-DNA.

Floriane Poignant↗

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of magnetic helicity upon rectilinear propagation of charged particles in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by any random field components that are superposed on the guiding field. If the random field configuration embodies helicity, the scattering is asymmetrical with respect to a plane perpendicular to the guiding field, for particles moving into the forward hemisphere are scattered at different rates from those moving into the backward hemisphere. This asymmetry gives rise to new terms in the transport equations that describe propagation of charged particles. Helicity has virtually no impact on qualitative features of the diffusive mode of propagation. However, characteristic velocities of the coherent modes that appear after a highly anisotropic injection exhibit an asymmetry related to helicity. Explicit formulas, which embody the effects of helicity, are given for the anisotropies, the coefficient diffusion, and the coherent velocities. Predictions derived from these expressions are in good agreement with Monte Carlo simulations of particle transport, but the simulations reveal certain phenomena whose explanation calls for further analytical work.

Earl, James A.↗

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗

A composite electrodynamic mechanism to reconcile spatiotemporally resolved exciton transport in quantum dot superlattices

Quantum dot (QD) solids are promising optoelectronic materials; further advancing their device functionality requires understanding their energy transport mechanisms. The commonly invoked near-field Förster resonance energy transfer (FRET) theory often underestimates the exciton hopping rate in QD solids, yet no consensus exists on the underlying cause. In response, we use time-resolved ultrafast stimulated emission depletion (STED) microscopy, an ultrafast transformation of STED to spatiotemporally resolve exciton diffusion in tellurium-doped cadmium selenide–core/cadmium sulfide–shell QD superlattices. We measure the concomitant time-resolved exciton energy decay due to excitons sampling a heterogeneous energetic landscape within the superlattice. The heterogeneity is quantified by single-particle emission spectroscopy. This powerful multimodal set of observables provides sufficient constraints on a kinetic Monte Carlo simulation of exciton transport to elucidate a composite transport mechanism that includes both near-field FRET and previously neglected far-field emission/reabsorption contributions. Uncovering this mechanism offers a much-needed unified framework in which to characterize transport in QD solids and additional principles for device design.

36 MATERIALS SCIENCE↗

Taking the Pulse of PyroCumulus Clouds

Forest fires can burn large areas, but can also inject smoke into the upper troposphere/lower stratosphere (UT/LS), where stakes are even higher for climate, because emissions tend to have a longer lifetime, and can produce significant regional and even global climate effects, as is the case with some volcanoes. Large forest fires are now believed to be more common in summer, especially in the boreal regions, where pyrocumulus (pyroCu), and occasionally pyrocumuionimbus (pyroCb) clouds are formed, which can transport emissions into the UT/LS. A major difficulty in developing realistic fire plume models is the lack of observational data within fire plumes that resolves structure at a few 100 m scales, which can be used to validate these models. Here, we report detailed airborne radiation measurements within strong pyroCu taken over boreal forest fires in Saskatchewan, Canada during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) summer field campaign in 2008. We find that the angular distribution of radiance within the pyroCu is closely related to the diffusion domain in water clouds and can be described by very similar simple cosine functions. We demonstrate with Monte Carlo simulations that radiation transport in pyroCu is inherently a 3D phenomenon and must account for particle absorption. However, the simple cosine function promises to offer an easy solution for climate models. The presence of a prominent smoke core, defined by strong extinction in the UV, VIS and NIR, suggests that the core might be an important pathway for emission transport to the upper troposphere and lower stratosphere. We speculate that this plume injection core is generated and sustained by complex processes not yet well understood, but not necessarily related directly to the intense fires that originally initiated the plume rise.

Gatebe, C. K.↗

Issues in the reconstruction of environmental doses on the basis of thermoluminescence measurements in the Techa riverside

The potential of thermoluminescence measurements of bricks from the contaminated area of the Techa river valley, Southern Urals, Russia, for reconstructing external exposures of affected population groups has been studied. Thermoluminescence dating of background samples was used to evaluate the age of old buildings available on the river banks. The anthropogenic gamma dose accrued in exposed samples is determined by subtracting the natural radiation background dose for the corresponding age from the accumulated dose measured by thermoluminescence. For a site in the upper Techa river region, where the levels of external exposures were extremely high, the depth-dose distribution in bricks and the dependence of accidental dose on the height of the sampling position were determined. For the same site, Monte Carlo simulations of radiation transport were performed for different source configurations corresponding to the situation before and after the construction of a reservoir on the river and evacuation of the population in 1956. A comparison of the results provides an understanding of the features of the measured depth-dose distributions and height dependencies in terms of the source configurations and shows that bricks from the higher sampling positions are likely to have accrued a larger fraction of anthropogenic dose from the time before the construction of the reservoir. The applicability of the thermoluminescent dosimetry method to environmental dose reconstruction in the middle Techa region, where the external exposure was relatively low, was also investigated.

NASA Discipline Radiation Health↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Impacts of Focused Ion Beam Processing on the Fabrication of Nanoscale Functionalized Probes

Herein, we examine the impact of Ga + ion kinetic energy and the target material type on the extent of ion implantation and structural damage in atomic force microscopy probes made of Al 2 O 3 and ZnO manufactured by focused ion beam (FIB) using scanning transmission electron microscopy and energy-dispersive X-ray mapping. Penetration of Ga into the Al 2 O 3 lattice induced structural distortions and amorphization. For ZnO probes, Ga is uniformly dispersed across the surface, resulting in the formation of distinct clusters. Atom probe tomography further validates the Ga distributions in Al 2 O 3 and ZnO nanoprobes. Complementary Monte Carlo simulations with the transport of ions in the matter program indicated that the introduction of Ga + prompts the generation of cation and anion vacancies, an occurrence more pronounced in Al 2 O 3 compared to ZnO. In conclusion, this study not only enriches the knowledge of ion-matter interactions, but also serves as a practical guide for the fabrication of nanoscale functionalized AFM probes.

Al2O3↗

Dependence of high- Z redeposition on the field-to-surface pitch angle and other sheath parameters in tokamaks

Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗