Neutron Transport Modeling using Monte Carlo Codes
What is Monte Carlo: Monte Carlo Versus Deterministic
36 MATERIALS SCIENCE↗
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.
The origin of meteorites as small bodies. ii - the model.
Monte Carlo technique to determine final orbital parameters and bombardment age distribution of meteoritic fragments ejected into inner solar system
ASTEROID↗
Monte Carlo lattice models for adsorbed polymer conformation
The adhesion between a polymer film and a metal surface is of great technological interest. However, the prediction of adhesion and wear properties of polymer coated metals is quite difficult because a fundamental understanding of the polymer surface interaction does not yet exist. A computer model for the conformation of a polymer molecule adsorbed on a surface is discussed. The chain conformation is assumed to be described by a partially directed random walk on a three dimensional simple cubic lattice. An attractive surface potential is incorporated into the model through the use of a random walk step probability distribution that is anisotropic in the direction normal to the attractive surface. The effects of variations in potential characteristics are qualitatively included by varying both the degree of anisotropy of the step distribution and the range of the anisotropy. Polymer conformation is characterized by the average end to end distance, average radius of gyration, and average number of chain segments adsorbed on the surface.
Good, B. S.↗
Observations and Models of the Lunar Sodium Exosphere 1988 - 1999
Sodium in the lunar exosphere is easily observed from the Earth's surface due to its strong resonance emission lines in the visible region of the spectrum. Although sodium is a trace element, it is easily ejected from the surface by a number of processes. The variation of this exospheric constituent both spatially and temporally can help to constrain these sources and the loss processes and their timescales. Due to a revival of interest in the Moon and its volatiles, observations of the lunar exosphere obtained at the McMath-Pierce solar telescope in 1998 and 1999 have recently been reduced and analyzed. In addition, observations of the lunar sodium exosphere obtained with the Mt. Lemmon Lunar Coronagraph on Mt. Lemmon, Arizona, have also been published. We combine these new data with data previously published and reanalyzed by Sarantos et al. This comprehensive data set will be modeled using both a simple Chamberlain exosphere model and a comprehensive Monte Carlo model.
Killen, Rosemary↗
Atomistic Multi-Lattice Kinetic Monte Carlo (Kmc) Modeling of Hyperthermal Oxidation of Multi-Layer Graphene
Explore the source record for details and available documents.
Sharon Edward↗
Monte Carlo Radiative Transfer Modeling of Lightning Observed in Galileo Images of Jupiter
We study lightning on Jupiter and the clouds illuminated by the lightning using images taken by the Galileo orbiter. The Galileo images have a resolution of 25 km/pixel and axe able to resolve the shape of the single lightning spots in the images, which have full widths at half the maximum intensity in the range of 90-160 km. We compare the measured lightning flash images with simulated images produced by our ED Monte Carlo light-scattering model. The model calculates Monte Carlo scattering of photons in a ED opacity distribution. During each scattering event, light is partially absorbed. The new direction of the photon after scattering is chosen according to a Henyey-Greenstein phase function. An image from each direction is produced by accumulating photons emerging from the cloud in a small range (bins) of emission angles. Lightning bolts are modeled either as points or vertical lines. Our results suggest that some of the observed scattering patterns axe produced in a 3-D cloud rather than in a plane-parallel cloud layer. Lightning is estimated to occur at least as deep as the bottom of the expected water cloud. For the six cases studied, we find that the clouds above the lightning are optically thick (tau > 5). Jovian flashes are more regular and circular than the largest terrestrial flashes observed from space. On Jupiter there is nothing equivalent to the 30-40-km horizontal flashes which axe seen on Earth.
Dyudine, U. A.↗
Linear Boltzmann transport for jet propagation in the quark-gluon plasma: Inelastic processes and jet modification
A linear Boltzmann transport (LBT) Monte Carlo model has been developed to describe jet propagation and interaction with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions. A complete set of elastic-scattering processes and medium-induced gluon emissions based on the higher-twist formalism are incorporated for both jet shower and medium recoil partons. It has been employed to describe experimental data on large transverse momentum hadron and jet spectra, correlation and jet substructures in high-energy heavy-ion collisions. Here, we document in detail the structure of the model and validation of the Monte Carlo implementations of the physics processes in LBT, in particular, the inelastic process of medium-induced gluon radiation. We carry out a comprehensive examination of the jet-medium interaction as implemented in LBT through energy loss and momentum broadening of a single hard parton, the energy and transverse momentum transfer from leading partons to medium-induced gluons and jet-induced medium excitation, and medium modification of reconstructed jets in a static and uniform medium. With realistic and event-by-event hydrodynamic medium in heavy-ion collisions, we compute and compare with experimental data on the jet cone-size dependence of the single inclusive jet suppression at both the BNL Relativistic Heavy-Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), the dijet asymmetry at the LHC and γ-jet correlation at RHIC. Effects of medium-induced gluon emissions and jet-induced medium excitation on jet observables are systematically examined. Rescatterings of the radiated gluons and recoil partons with the QGP are found essential to account for the enhancement of soft particle yield toward the edge of the jet cone.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
An analysis of the OI 1304 A dayglow using a Monte Carlo resonant scattering model with partial frequency redistribution
The transport of resonance radiation under optically thick conditions is shown to be accurately described by a Monte Carlo model of the atomic oxygen 1304 A airglow triplet in which partial frequency redistribution, temperature gradients, pure absorption and multilevel scattering are accounted for. All features of the data can be explained by photoelectron impact excitation and the resonant scattering of sunlight, where the latter source dominates below 100 and above 500 km and is stronger at intermediate altitudes than previously thought. It is concluded that the OI 1304 A emission can be used in studies of excitation processes and atomic oxygen densities in planetary atmospheres.
Meier, R. R.↗
A Monte Carlo-finite element model for strain energy controlled microstructural evolution - 'Rafting' in superalloys
This paper presents a specialized microstructural lattice model, MCFET (Monte Carlo finite element technique), which simulates microstructural evolution in materials in which strain energy has an important role in determining morphology. The model is capable of accounting for externally applied stress, surface tension, misfit, elastic inhomogeneity, elastic anisotropy, and arbitrary temperatures. The MCFET analysis was found to compare well with the results of analytical calculations of the equilibrium morphologies of isolated particles in an infinite matrix.
Gayda, J.↗
Improving Multi-Model Trajectory Simulation Estimators using Model Selection and Tuning
Multi-model Monte Carlo methods have been demonstrated to be an efficient and accurate alternative to standard Monte Carlo (MC) in the model-based propagation of uncertainty in entry, descent, and landing (EDL) applications. These multi-model MC methods fuse predictions from low-fidelity models with the high-fidelity EDL model of interest to produce unbiased statistics with a fraction of the computational cost. The accuracy and efficiency of the multi-model MC methods are dependent upon the magnitude of correlations of the low-fidelity models with the high-fidelity model, but also upon the correlation among the low-fidelity models, and their relative computational cost. Because of this layer of complexity, the question of how to optimally select the set of low-fidelity models has remained open. In this work, methods for optimal model construction and tuning are investigated as a means to increase the speed and precision of trajectory simulation for EDL. Specifically, the focus is on the inclusion of low-fidelity model tuning within the sample allocation optimization that accompanies multi-model MC methods. Preliminary results indicate that low-fidelity model tuning can significantly improve efficiency and precision of trajectory simulations and provide an increased edge to multi-model MC methods when compared to standard MC. The challenges and potential benefits to exploring a fully iterative and comprehensive optimization strategy in future work are highlighted.
uncertainty quantification↗
Resolving Atomic-Scale Structure and Chemical Coordination in High-Entropy Alloy Electrocatalysts for Structure–Function Relationship Elucidation
The recent breakthrough in confining five or more atomic species in nanocatalysts, referred to as high-entropy alloy nanocatalysts (HEAs), has revealed the possibilities of multielemental interactions that can surpass the limitations of binary and ternary electrocatalysts. The wide range of potential surface configurations in HEAs, however, presents a significant challenge in resolving active structural motifs, preventing the establishment of structure-function relationships for rational catalyst design and optimization. Here, we present a methodology for creating sub-5 nm HEAs using an aqueous-based peptide-directed route. Using a combination of pair distribution function and X-ray absorption spectroscopy, HEA structure models are constructed from reverse Monte Carlo modeling of experimental data sets and showcase a clear peptide-induced influence on atomic-structure and chemical miscibility. Coordination analysis of our structure models facilitated the construction of structure-function correlations applied to electrochemical methanol oxidation reactions, revealing the complex interplay between multiple metals that leads to improved catalytic properties. Our results showcase a viable strategy for elucidating structure-function relationships in HEAs, prospectively providing a pathway for future materials design.
36 MATERIALS SCIENCE↗
Feature Profile Evolution of SiO2 Trenches In Fluorocarbon Plasmas
Etching of silicon microstructures for semiconductor manufacturing in chlorine plasmas has been well characterized. The etching proceeds in a two-part process, where the chlorine neutrals passivate the Si surface and then the ions etch away SiClx. However, etching in more complicated gas mixtures and materials, such as etching of SiO2 in Ar/C4F8, requires knowledge of the ion and neutral distribution functions as a function of angle and velocity, in addition to modeling the gas surface reactions. In order to address these needs, we have developed and integrated a suite of models to simulate the etching process from the plasma reactor level to the feature profile evolution level. This arrangement allows for a better understanding, control, and prediction of the influence of equipment level process parameters on feature profile evolution. We are currently using the HPEM (Hybrid Plasma Equipment Model) and PCMCM (Plasma Chemistry Monte Carlo Model) to generate plasma properties and ion and neutral distribution functions for argon/fluorocarbon discharges in a GEC Reference Cell. These quantities are then input to the feature scale model, Simulation of Profile Evolution by Level Sets (SPELS). A surface chemistry model is used to determine the interaction of the incoming species with the substrate material and simulate the evolution of the trench profile. The impact of change of gas pressure and inductive power on the relative flux of CFx and F to the wafer, the etch and polymerization rates, and feature profiles will be examined. Comparisons to experimental profiles will also be presented.
Hwang, Helen↗
The 2004 NASA Aerospace Battery Workshop
Topics covered include: Super NiCd(TradeMark) Energy Storage for Gravity Probe-B Relativity Mission; Hubble Space Telescope 2004 Battery Update; The Development of Hermetically Sealed Aerospace Nickel-Metal Hydride Cell; Serial Charging Test on High Capacity Li-Ion Cells for the Orbiter Advanced Hydraulic Power System; Cell Equalization of Lithium-Ion Cells; The Long-Term Performance of Small-Cell Batteries Without Cell-Balancing Electronics; Identification and Treatment of Lithium Battery Cell Imbalance under Flight Conditions; Battery Control Boards for Li-Ion Batteries on Mars Exploration Rovers; Cell Over Voltage Protection and Balancing Circuit of the Lithium-Ion Battery; Lithium-Ion Battery Electronics for Aerospace Applications; Lithium-Ion Cell Charge Control Unit; Lithium Ion Battery Cell Bypass Circuit Test Results at the U.S. Naval Research Laboratory; High Capacity Battery Cell By-Pass Switches: High Current Pulse Testing of Lithium-Ion; Battery By-Pass Switches to Verify Their Ability to Withstand Short-Circuits; Incorporation of Physics-Based, Spatially-Resolved Battery Models into System Simulations; A Monte Carlo Model for Li-Ion Battery Life Projections; Thermal Behavior of Large Lithium-Ion Cells; Thermal Imaging of Aerospace Battery Cells; High Rate Designed 50 Ah Li-Ion Cell for LEO Applications; Evaluation of Corrosion Behavior in Aerospace Lithium-Ion Cells; Performance of AEA 80 Ah Battery Under GEO Profile; LEO Li-Ion Battery Testing; A Review of the Feasibility Investigation of Commercial Laminated Lithium-Ion Polymer Cells for Space Applications; Lithium-Ion Verification Test Program; Panasonic Small Cell Testing for AHPS; Lithium-Ion Small Cell Battery Shorting Study; Low-Earth-Orbit and Geosynchronous-Earth-Orbit Testing of 80 Ah Batteries under Real-Time Profiles; Update on Development of Lithium-Ion Cells for Space Applications at JAXA; Foreign Comparative Technology: Launch Vehicle Battery Cell Testing; 20V, 40 Ah Lithium Ion Polymer Battery for the Spacesuit; Low Temperature Life-Cycle Testing of a Lithium-Ion Battery for Low-Earth-Orbiting Spacecraft; and Evaluation of the Effects of DoD and Charge Rate on a LEO Optimized 50 Ah Li-Ion Aerospace Cell.
Source record↗
Improving Multi-Model Trajectory Simulation Estimators using Model Selection and Tuning
Multi-model Monte Carlo methods have been demonstrated to be an efficient and accurate alternative to standard Monte Carlo (MC) in the model-based propagation of uncertainty in entry, descent, and landing (EDL) applications. These multi-model MC methods fuse predictions from low-fidelity models with the high-fidelity EDL model of interest to produce unbiased statistics with a fraction of the computational cost. The accuracy and efficiency of the multi-model MC methods are dependent upon the magnitude of correlations of the low-fidelity models with the high-fidelity model, but also upon the correlation amongst the low-fidelity models, and their relative computational cost. Because of this layer of complexity, the question of how to optimally select the set of low-fidelity models has remained open. In this work, methods for optimal model construction and tuning are investigated as a means to increase the speed and precision of trajectory simulation for EDL. Specifically, the focus is on the inclusion of low-fidelity model tuning within the sample allocation optimization that accompanies multi-model MC methods. Preliminary results indicate that low-fidelity model tuning can significantly improve efficiency and precision of trajectory simulations and provide an increased edge to multi-model MC methods when compared to standard MC. The challenges and potential benefits to exploring a fully iterative and comprehensive optimization strategy in future work are highlighted.
uncertainty quantification↗
The importance of delayed neutron precursors in gamma dose calculations for activated primary heat exchanger components in the Molten Salt Breeder Reactor
The Molten Salt Reactor (MSR) Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation program has supported the development of 3D Monte Carlo models of the Molten Salt Breeder Reactor (MSBR) over the last couple of years. This MSBR model was previously run with the Shift Monte Carlo code to perform radiation shielding calculations in the reactor cell area. The MSBR is a 2250 MWth (1000 MWe) liquid-fueled molten salt reactor design developed at Oak Ridge National Laboratory in the 1970s. Determining the source terms from activated primary heat exchanger (HX) components is important because delayed and prompt neutron fluxes incident on these components affect the dose rate in the primary HX maintenance areas. This information can be used in the development of remote handling procedures required during shutdown for maintenance. A methodology has been developed and is proposed in this paper to quantify the activated source term from the primary HX components as a result of the movement of the delayed neutron precursors in flowing primary fuel salt through the primary HXs in the MSBR. The goal of this research is to evaluate the gamma dose rates in the maintenance hatches above the primary HX using the activated HX source terms. The study showed that the gamma dose rates are approximately two orders of magnitudes higher when accounting for the neutron activation from the movement of delayed neutron precursors through the HXs than when flowing fuel is not considered. Thus, the movement of delayed neutron precursors must be taken into account for accurately predicting the neutron activation of primary loop components.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A model for astrophysical spallation reactions
A Monte-Carlo model (RENO) for spallation reactions is described which can treat both the spallations induced by a free nucleon and those induced by a complex nucleus. It differs from other such models in that it employs a discrete-nucleon representation of the nucleus and allows clusters of nucleons to form and to participate in the reaction. The RENO model is particularly suited for spallations involving the relatively light nuclei of astrophysical and cosmic-ray interest.
Schmitt, W. F.↗
A triple junction energy study using an inclination-dependent anisotropic Monte Carlo Potts grain growth model
This work presents a Monte Carlo Potts grain growth model in which the grain boundary (GB) energies depend on the GB inclination. The inclination is calculated using a linear smoothing approach developed by the authors. In bicrystal simulations with a shrinking grain, the grain changes shape to prefer low energy GB inclinations. However, in polycrystal simulations the preferred inclinations depend on the approach used to assign the triple junction (TJ) energies. Approaches that produce unimodal TJ energy distributions result in the expected behavior of preferring low energy GB inclinations. However, approaches that produce bimodal TJ energy distributions result in medium energy or even high energy inclinations being preferred. Overall, this study underscores the importance of TJs in anisotropic grain growth.
Grain boundary inclination↗
Radiation modeling with direct simulation Monte Carlo
Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available experimental data. The new method is compared with the current modeling technique and both techniques are compared with available data. The differences in the results are evaluated. The test case is based on an AVCO-Everett shock tube experiment, a 10-km/s standing shock wave in air at 0.1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to different flow conditions.
Carlson, Ann B.↗