Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo modeling”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Investigation of the significance of excited species in wave-heated dense non-ideal microplasmas using particle-in-cell Monte Carlo-collision modeling

A computational model for the interaction of a high-pressure microplasma with an electromagnetic wave is presented. A one-dimensional particle-in-cell Monte Carlo collision model is used to investigate the plasma non-ideality effects in a second-stage laser-heated xenon plasma with a comprehensive chemistry mechanism, including excited species. A xenon microplasma at a temperature of 300 K and pressure of 10 bar becomes non-ideal after the laser heating stage with lower ionization and electronic excitation barriers, which makes the presence of excited species important to the plasma generation process. For these conditions, two-step ionization via excited species becomes a major contributor to plasma ionization. We find that full ionization of the plasma is obtained in ~2 ps when excited species are included in the chemistry model, as opposed to ~2.5 ps for a chemistry mechanism that does not consider their effect. With excited species, the ionization mode transitions from direct ionization via electron-neutral collisions to direct and stepwise ionization as the plasma generation progresses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Studying strangeness and baryon production mechanisms through angular correlations between charged Ξ baryons and identified hadrons in pp collisions at $\sqrt{s}$ = 13 TeV

The angular correlations between charged Ξ baryons and associated identified hadrons (pions, kaons, protons, Λ baryons, and Ξ baryons) are measured in pp collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector to give insight into the particle production mechanisms and balancing of quantum numbers on the microscopic level. In particular, the distribution of strangeness is investigated in the correlations between the doubly-strange Ξ baryon and mesons and baryons that contain a single strange quark, K and Λ. As a reference, the results are compared to Ξπ and Ξp correlations, where the associated mesons and baryons do not contain a strange valence quark. These measurements are expected to be sensitive to whether strangeness is produced through string breaking or in a thermal production scenario. Furthermore, the multiplicity dependence of the correlation functions is measured to look for the turn-on of additional particle production mechanisms with event activity. The results are compared to predictions from the string-breaking model P YTHIA 8, including tunes with baryon junctions and rope hadronisation enabled, the cluster hadronisation model H ERWIG 7, and the core-corona model E POS-LHC . While some aspects of the experimental data are described quantitatively or qualitatively by the Monte Carlo models, no model can match all features of the data. These results provide stringent constraints on the strangeness and baryon number production mechanisms in pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Forward and inverse predictive transient models of TREAT using surrogate reactivity models

In this work, we present two novel approaches for predicting the power evolution and control rod height of the Transient Reactor Test Facility (TREAT) to support experiment modeling; specifically, transient analysis of the NASA-sponsored Sirius series of experiments. These approaches utilize steady-state Monte Carlo model, point kinetics model, and surrogate models to predict power evolution and control rod axial position during transient experiments. Both approaches were tested and validated against several Sirius experiments that were performed in the TREAT facility at different power levels. The validation test results show very good agreement with the experimental data, and the models were able to accurately predict the power evolution and the axial control rod position with an average error within 3.0%. Finally, this indicates that these approaches will help the reactor engineering team of the TREAT facility in preparing and predicting the power and temperature of the experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1996

High-Multiplicity Lead-Lead Interactions at 158 GeV/c per nucleon

The Krakow-Louisiana-Minnesota-Moscow Collaboration (KLMM) has exposed a set of emulsion chambers with lead targets to a 158 GeV/c per nucleon beam of Pb-208 nuclei, and we report the initial analysis of 40 high-multiplicity Pb-Pb collisions. To test the validity of the superposition model of nucleus-nucleus interactions in this new regime, we compare the shapes of the pseudorapidity distributions with FRITIOF Monte Carlo model calculations, and find close agreement for even the most central events. We characterize head-on collisions as having a mean multiplicity of 1550 +/- 120 and a peak pseudorapidity density of 390 +/- 30. These estimates are significantly lower than our FRITIOF calculations.

Deines-Jones, P.↗
DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Heating of the moon by heterogeneous accretion

Monte Carlo models of lunar accretion are developed by randomly sampling mass and velocity distributions of objects in the solar nebula in circumsolar and circumterrestrial orbits. The thermal effect of an impact is expressed as the volume of melt produced; the melt volumes produced by each Monte Carlo case are sorted into regions to produce percentage melt estimates as a function of radius. One-hundred cases are calculated to obtain a range of plausible thermal histories for lunar accretion. The average melt percentage estimates show substantial deep melting, but the spread of melt percentage values about the mean from the random sampling procedure is wide enough to include 0% melting for radii below 100 km and 100% melting for radii below 500 km.

Ransford, G. A.↗
NASA NTRS2003

Modeling Charge Collection in Detector Arrays

A detector array charge collection model has been developed for use as an engineering tool to aid in the design of optical sensor missions for operation in the space radiation environment. This model is an enhancement of the prototype array charge collection model that was developed for the Next Generation Space Telescope (NGST) program. The primary enhancements were accounting for drift-assisted diffusion by Monte Carlo modeling techniques and implementing the modeling approaches in a windows-based code. The modeling is concerned with integrated charge collection within discrete pixels in the focal plane array (FPA), with high fidelity spatial resolution. It is applicable to all detector geometries including monolithc charge coupled devices (CCDs), Active Pixel Sensors (APS) and hybrid FPA geometries based on a detector array bump-bonded to a readout integrated circuit (ROIC).

Hardage, Donna↗
DOE OSTI2021

The Initial–Final Mass Relation for Hydrogen-deficient White Dwarfs

The initial–final mass relation represents the total mass lost by a star during the entirety of its evolution from the zero age main sequence to the white-dwarf cooling track. The semiempirical initial–final mass relation (IFMR) is largely based on observations of DA white dwarfs, the most common spectral type of white dwarf and the simplest atmosphere to model. We present a first derivation of the semiempirical IFMR for hydrogen-deficient (non-DA) white dwarfs in open star clusters. We identify a possible discrepancy between the DA and non-DA IFMRs, with non-DA white dwarfs ≈0.07 M {sub ⊙} less massive at a given initial mass. Such a discrepancy is unexpected based on theoretical models of non-DA formation and observations of field white dwarf mass distributions. If real, the discrepancy is likely due to enhanced mass loss during the final thermal pulse and renewed post-AGB evolution of the star. However, we are dubious that the mass discrepancy is physical and instead is due to the small sample size, to systematic issues in model atmospheres of non-DAs, and to the uncertain evolutionary history of Procyon B (spectral type DQZ). A significantly larger sample size is needed to test these assertions. In addition, we also present Monte Carlo models of the correlated errors for DA and non-DA white dwarfs in the initial–final mass plane. We find the uncertainties in initial–final mass determinations for individual white dwarfs can be significantly asymmetric, but the recovered functional form of the IFMR is grossly unaffected by the correlated errors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Innovative dead-time correction and background subtraction for neutron multiplicity measurements using neural networks

Abstract The number of neutrons emitted from a nuclear reaction plays a crucial role in various fields, including nuclear theory, nuclear nonproliferation, nuclear energy and nuclear criticality safety. Accurate determination of neutron multiplicities requires the application of several corrections, with dead-time correction and background subtraction being particularly significant. These corrections become more challenging for neutron detectors with time-dependent neutron capture. In this work, we perform a comprehensive study of three existing methods used for dead-time correction and background subtraction in neutron detectors with time-dependent neutron capture. The methods were tested for dead-times in the range from 0 to 1 μs using a Monte Carlo model simulating the dead-time and background effects in the standard neutron multiplicity probability distribution of $$^{252}$$ 252 Cf. The previous methods showed larger than desired uncertainty or systematic trade off. Those uncertainties prompted the development of a novel approach using neural networks trained with data from Monte Carlo simulations. The Neural Network method enabled the correction of neutron multiplicity probabilities more accurately than the other methods with fractional errors smaller than 3% for multiplicities around the peak of $$^{252}$$ 252 Cf. A similar approach using neural networks could be applied to problems where the system being studied can be accurately simulated without having an accurate analytical description available. The neural network method presented in this paper can be easily expanded if multiplicities greater than 10 are expected.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Understanding Temporal and Spatial Variability of the Lunar Helium Atmosphere Using Simultaneous Observations from LRO, LADEE, and ARTEMIS

Simultaneous measurements of helium in the exosphere of the Moon are made from the Lunar Reconnaissance Orbiter (LRO) Lyman Alpha Mapping Project (LAMP) and the Lunar Atmosphere and Dust Environment Explorer (LADEE) Neutral Mass Spectrometer (NMS) through the entire 5-month span of the LADEE mission. In addition, the ARTEMIS mission monitored the solar wind alpha particle flux to the Moon. Modeling the lunar helium exosphere, we relate the LAMP polar observations to the LADEE equatorial observations. Further, using the ARTEMIS alpha flux in the Monte Carlo model reproduces the temporal variations in helium density. Comparing the LAMP data to the LADEE data shows excellent agreement. Comparing those with the ARTEMIS data reveals that the solar wind alpha flux is the primary driver to variability in the helium exosphere throughout the LADEE mission. Using a decay time for exospheric helium of 5 days, we determine that the solar wind contributes 64 +/- 5% of the helium to the lunar exosphere. The remaining 36 +/- 5% is presumed to come from outgassing of radiogenic helium from the interior of the Moon. Furthermore, the model reproduces the measurements if 63 +/- 6% of the incident alpha particles are converted to thermalized helium atoms through the interaction between the alphas and the lunar surface. However, these values are dependent on both inferred source rates from LAMP and LADEE observations and on the assumed time constant of the exospheric decay rate.

Moon↗
NASA NTRS1979

Visible flux variations across finite clouds

A radiative transfer model developed to computer visible flux variations across and within finite clouds of varying shapes is described. The technique has much of the versatility of Monte Carlo models, as well as the speed of the analytical finite cloud models.

Mosher, F. R.↗
DOE OSTI2026

ACRRF High-Bay Dose Calculations using MCNP (Part A)

Analytical tools and models have been developed as a starting point for directly assessing dose in the Annular Core Research Reactor Facility (ACRRF) due to reactor operation. Key results include peak dose along the Central Cavity (CC) centerline (beamline) at the cavity level, dose throughout the High-Bay (HB), and dose on the facility roof for partially-shielded reactor operation where the 4” insert is removed from the CC Shield Plug (SP). Model results in the beamline are benchmarked against measured doses from passive dosimetry evaluations. Personnel total (neutron and gamma) dose in the ACRRF HB is calculated using Monte Carlo N-Particle (MCNP). Various CC and SP configurations are analyzed, including unshielded (no SP) and partially shielded (SP installed but 4” insert removed). Novel application of Variance Reduction (VR) techniques, namely the Surface Source Write (SSW) and Surface Source Read (SSR) capabilities in MCNP, enable impressive resolution (in a Monte Carlo modeling sense) of dose throughout much the facility. The VR techniques reduce stochastic error for challenging tallies, with more advanced techniques explored in the companion to this report (Part B) [1]. Supplementary studies (including a verification analysis) and pedagogic evaluations in Part B involve neutron spectra, angular distributions, and the dose impact of facility characteristics. With the SP 4” insert removed and the Lead-Boron (44”) Bucket (LB–44) in the reactor cavity, Total Effective Dose (TED) within the CC beamline is ≈140 rem per 300 MJ of reactor yield (or 3900 rem per hour at 100% Steady-State (SS) power). With no SP (unshielded) and a Free-Field (FF) cavity, TED within the beamline is ≈610 rem per 300 MJ (or 17000 rem per hour at 100% SS power). Due to the predicted collimation of radiation by the reactor pool (and partial SP, if present), beamline dose is much greater than the scattered radiation field surrounding the cavity and reactor tank. Comparisons are made to beamline dosimetry measurements to validate the model. Model predictions agree reasonably well (⪅10%) with measured quantities of neutron fluence, gamma fluence, and spectral metrics. Away from the beamline, comparisons made to previous dose measurements in the HB agreement within an order of magnitude.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1974

A three-dimensional study of 30- to 300-MeV atmospheric gamma rays

A three-dimensional study of atmospheric gamma rays with energy greater than 30 MeV has been carried out. A knowledge of these atmospheric secondaries has significant applications to the study of cosmic gamma rays. For detectors carried on balloons, atmospherically produced gamma rays are the major source of background. For satellite detectors, atmospheric secondaries provide a calibration source. Experimental results were obtained from four balloon flights from Palestine, Texas, with a 15 cm by 15 cm digitized wire grid spark chamber. The energy spectrum for downward-moving gamma rays steepens with increasing atmospheric depth. Near the top of the atmosphere, the spectrum steepens with increasing zenith angle. A new model of atmospheric secondary production has calculated the depth, the energy, and the zenith angle dependence of gamma rays above 30 MeV, using a comprehensive three-dimensional Monte Carlo model of the nucleon-meson-electromagnetic cascade.

Thompson, D. J.↗
NASA NTRS2001

Development of Numerical Models for Performance Predictions of Single-Photon Avalanche Photodetectors (SPAP) for the 2-Micron Regime

Field dependent drift velocity results are presented for electron transport in bulk Indium Arsenide (InAs) material based on a Monte Carlo model, which includes an analytical treatment of band-to-band impact ionization. Avalanche multiplication and related excess noise factor (F) are computed as a function of device length and applied voltage. A decrease in F with increases in device length is obtained. The results suggest an inherent utility for InAs-based single-photon avalanche detectors, particularly around the 2 microns region of interest for atmospheric remote sensing applications. The dark current response was also evaluated. The role of the various components has been analyzed. For shorter devices, the tunneling component is shown to dominate at low temperatures. Finally, possible structures for enhanced photodetection are proposed for future research.

Joshi, Ravindra P.↗
DOE OSTI2020

Multivariate discrimination in quantum target detection

In this work, we describe a simple multivariate technique of likelihood ratios for improved discrimination of signal and background in multi-dimensional quantum target detection. Furthermore, the technique combines two independent variables, time difference and summed energy, of a photon pair from the spontaneous parametric down-conversion source into an optimal discriminant. The discriminant performance was studied in experimental data and in Monte-Carlo modelling with clear improvement shown compared to previous techniques. As novel detectors become available, we expect this type of multivariate analysis to become increasingly important in multi-dimensional quantum optics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2012

Modeling MESSENGER Observations of Calcium in Mercury's Exosphere

The Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MESSENGER spacecraft has made the first high-spatial-resolution observations of exospheric calcium at Mercury. We use a Monte Carlo model of the exosphere to track the trajectories of calcium atoms ejected from the surface until they are photoionized, escape from the system, or stick to the surface. This model permits an exploration of exospheric source processes and interactions among neutral atoms, solar radiation, and the planetary surface. The MASCS data have suggested that a persistent, high-energy source of calcium that was enhanced in the dawn, equatorial region of Mercury was active during MESSENGER's three flybys of Mercury and during the first seven orbits for which MASCS obtained data. The total Ca source rate from the surface varied between 1.2x10(exp 23) and 2.6x10(exp 23) Ca atoms/s, if its temperature was 50,000 K. The origin of this high-energy, asymmetric source is unknown, although from this limited data set it does not appear to be consistent with micrometeoroid impact vaporization, ion sputtering, electron-stimulated desorption, or vaporization at dawn of material trapped on the cold nightside.

Messenger↗
DOE OSTI2022

XE-100 modeling and simulation for neutronic analysis in MCNP6.2

XE-100 is a generation IV helium-cooled, graphite-moderated, pebble-bed reactor (HTGR). As part of the pathway toward a conceptually designing and licensing this reactor, an independent Monte Carlo model was created in MCNP6.2, and several distinct neutronic analyses were then performed. The double heterogeneity of TRISO fuel within graphite pebbles introduces unique modeling challenges related to particle and pebble clipping. The results show that for neutron and photon heating of ex-core components such as the reflector, RCSS, core barrel (CB), the model that contains clipping produces higher heating values. It is therefore concluded that removing clipping via compression of the particles and pebbles within the model distributes the neutrons and gammas preferentially toward the core center, and reduces the heating that is experienced toward the reactor periphery. Thus, the most conservative model for ex-core heating contains particle and pebble clipping. Also presented are results on the impact of chamfers that exist on the corners of graphite reflector blocks. As these chamfers could potentially create streaming paths, the neutron and gamma flux from the core to the CB were analyzed. It was determined that the chamfers do not significantly impact the neutron or gamma signatures on the CB, in that the shape of the neutron and photon flux on a detector imposed on the CB shows no preferential streaming path. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗