SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Energy conversion and transfer for plasmas in a magnetic expansion configuration

A two-dimensional axisymmetric particle-in-cell code with Monte Carlo collision conditions has been used to study particle energy transfer in plasmas and conversion in applied magnetic and electric fields appropriate to coaxial acceleration. The research incorporates a computation scheme with: a model of single particle magnetic interactions; a model of single particle interactions in electric and magnetic fields; and a model of multi-particle collisional interactions in order to understand the energy transfer processes and conversion mechanisms of charged plasma particles. This approach predicts electron and ion motions along with their energy variations for physical conditions that occur in the related models; the results allow comparison with experimental data for magnetic field strengths of 0.01–0.05 T and electrode voltages of 22.0–32.0 V. With the incorporation of magnetic and electric field effects on charged particles, the multi-particle model includes electron-neutral ionization collisions, ion-neutral charge exchange collisions, and electron-ion Coulomb collisions. This research presents a new approach to achieve an underlying understanding of the plasma energy transfer and conversion in the external electric and magnetic fields that is not possible using magnetohydrodynamics continuum representations. Results indicate the following innovative conclusions: (1) Radial and azimuthal energies of magnetized electrons are converted into an axial electron energy component in the diverging magnetic field, and the azimuthal kinetic energy of unmagnetized ions is converted into axial and radial components. (2) In electric and magnetic fields, electric field energy is primarily converted into axial kinetic energy of magnetized electrons by the energy transformation effects of magnetic fields, and for unmagnetized ions, the radial kinetic energy component dominates in the conversion of electric field energy. (3) For the collisional plasma, electron kinetic energy tends to increase (or decrease) to a terminal value since electrons lose energy in collisions then gain energy again from the field. Ions acquire most energy directly from the electric field, although part of the electric field energy arrives to the ions by collisions. Further, the ion axial energy component dominates the total ion energy. The collision processes are found to be integral and essential for the conversion of the plasma non-directed energy gain to be converted into the resultant axial energy, the magnitudes of which are found to be in agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Control rod calibration and reactivity effects at the IPEN/MB-01 reactor

Researches that aim to improve the performance of neutron transport codes and quality of nuclear cross section databases are very important to increase the accuracy of simulations and the quality of the analysis and prediction of phenomena in the nuclear field. In this context, relevant experimental data such as reactivity worth measurements are needed. Control rods may be made of several neutron absorbing materials that are used to adjust the reactivity of the core. For the reactor operation, these experimental data are also extremely important: with them it is possible to estimate the reactivity worth by the movement of the control rod, understand the reactor response at each rod position and to operate the reactor safely. This work presents a temperature correction approach for the control rod calibration problem. It is shown the control rod calibration data of the IPEN/MB-01 reactor, the integral and differential reactivity curves and a theoretical analysis, performed by the MCNP-5 reactor physics code, developed and maintained by Los Alamos National Laboratory, using the ENDF/B-VII.0 nuclear data library.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Design considerations for three-dimensional betavoltaics

Betavoltaic devices are suitable for delivering low-power over periods of years. Typically, their power density is on the order of nano to micro-Watts per cubic centimeter. In this work we evaluate the potential for using high-aspect ratio three-dimensional semiconductor structures to enhance the power and efficiency of these devices. The Monte Carlo transport code MCNP6 is used to provide realistic estimates of the theoretical levels of charge generation, which is in turn used to make predictions about the power output from three-dimensional betavoltaics. The focus of this work is on silicon and promethium-147, but other semiconductors and radioisotopes are considered as well. In the case of silicon diodes with three-dimensional features that are comparable to what is commercially available we estimate that power densities in the range of 20-25 mW/cm 3 can be achieved at efficiencies of 2.9-5.8% when coupled with promethium-147 oxide.

25 ENERGY STORAGE↗

DOE OSTI2024

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↗

DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗

DOE OSTI2019

Systematic Measurements and Analyses for Lead Void Reactivity Worth in a Plutonium Core and Two Uranium Cores with Different Enrichments

To validate lead (Pb) nuclear cross sections, a series of integral experiments to measure lead void reactivity worth was conducted systematically in three fast neutron spectra with different fuel compositions on the Comet critical assembly of the National Criticality Experiments Research Center. Previous experiments in high-enriched uranium (HEU)/Pb and low-enriched uranium (LEU)/Pb systems had been performed in 2016 and 2017, respectively. A follow-on experiment in a plutonium (Pu)/Pb system has been completed. The Pu/Pb system was constructed using lead plates and weapons-grade Pu plates that had been used in the Zero Power Physics Reactor (ZPPR) of Argonne National Laboratory until the 1990s. Furthermore, the HEU/Pb system was reexamined on the Comet critical assembly with a newly installed device that can measure the compression of the stack, improving reproducibility. Using the lead void reactivity worth measured in these three cores with different fuel compositions, the latest nuclear data libraries, JENDL-4.0 and ENDF/B-VIII.0, were tested with the Monte Carlo calculation code MCNP® version 6.1. As a result, the calculations by ENDF/B-VIII.0 were confirmed to agree with lead void reactivity worth measured in all the cores. It was furthermore found that the calculations by JENDL-4.0 overestimate by more than 20% for the Pu/Pb core while being in good agreement for the HEU/Pb and LEU/Pb cores.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2013

RBS/C, HRTEM and HRXRD study of damage accumulation in irradiated SrTiO3

Damage accumulation in argon-irradiated SrTiO3 single crystals has been studied by using combination of Rutherford Backscattering/Channeling (RBS/C), High Resolution Transmission Electron Microscopy (HRTEM) and High Resolution X-Ray Diffraction (HRXRD) techniques. The RBS/C spectra were fitted using McChasy, a Monte Carlo simulation code allowing the quantitative analysis of amorphous-like and dislocation-like types of defects. The results were interpreted by using a Multi-Step Damage Accumulation model which assumes, that the damage accumulation occurs in a series of structural transformations, the defect transformations are triggered by a stress caused by formation of a free volume in the irradiated crystal. This assumption has been confirmed by High Resolution Transmission Electron Microscopy and High Resolution X-Ray Diffraction analysis.

damage accumulation, radiation defects, oxides, ir↗

DOE OSTI2023

Comprehensive Analysis of Streaming and Shutdown Dose Rate Experiments at JET with ORNL Fusion Neutronics Workflows

Current experimental fusion systems and conceptual designs of fusion pilot plants (FPPs) are growing in complexity and size. Several radiation metrics are crucial to the safe operation of fusion machines, including neutron flux streaming through openings and the shutdown dose rate (SDDR). Most current designs of advanced experimental fusion systems—and the most probable candidates for FPPs—are based on the tokamak concept, which is prone to neutron streaming through the myriad openings needed for diagnostic and support systems. SDDR is caused by decay gamma rays from radionuclides that become activated by neutrons during the operation of a fusion system that use deuterium-deuterium (DD), tritium-tritium, or deuterium-tritium plasma. Because computational tools have become essential for determining these radiation metrics, they must be validated against reliable and applicable experimental data. Experiments at the Joint European Torus (JET) provide a unique source of experimental data for validating computational tools and nuclear data used to determine SDDR and neutron fluxes in streaming-dominated geometries. Here, this paper presents the comprehensive analysis of the high-performance DD JET SDDR, and streaming experiments performed using Oak Ridge National Laboratory (ORNL) fusion workflows. The computational results were compared with experimental results that consist of online SDDR measurements with ionization chambers and neutron fluence streaming measurements using thermoluminescent detectors. The ratio of calculated-to-experimental SDDR values ranges from 0.6 to 2.5, and the streaming results range from 0.5 to 8.0. Future work will include analyzing the JET 2021 DTE2 campaign alongside the integration of the Shift Monte Carlo transport code into all ORNL fusion neutronics workflows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

THE EVOLUTION OF SOLAR FLUX FROM 0.1 nm TO 160 {mu}m: QUANTITATIVE ESTIMATES FOR PLANETARY STUDIES

Understanding changes in the solar flux over geologic time is vital for understanding the evolution of planetary atmospheres because it affects atmospheric escape and chemistry, as well as climate. We describe a numerical parameterization for wavelength-dependent changes to the non-attenuated solar flux appropriate for most times and places in the solar system. We combine data from the Sun and solar analogs to estimate enhanced UV and X-ray fluxes for the young Sun and use standard solar models to estimate changing visible and infrared fluxes. The parameterization, a series of multipliers relative to the modern top of the atmosphere flux at Earth, is valid from 0.1 nm through the infrared, and from 0.6 Gyr through 6.7 Gyr, and is extended from the solar zero-age main sequence to 8.0 Gyr subject to additional uncertainties. The parameterization is applied to a representative modern day flux, providing quantitative estimates of the wavelength dependence of solar flux for paleodates relevant to the evolution of atmospheres in the solar system (or around other G-type stars). We validate the code by Monte Carlo analysis of uncertainties in stellar age and flux, and with comparisons to the solar proxies {kappa}{sup 1} Cet and EK Dra. The model is applied to the computation of photolysis rates on the Archean Earth.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

A simplified lithium vapor box divertor

A detached divertor plasma is predicted to be necessary to mitigate the heat fluxes and ion energy, and so sputtering, at the divertor plate of a fusion power plant. It has proven difficult in current experiments, however, to prevent the detachment front, once formed, from running up to the main plasma resulting in deterioration of pedestal performance. The lithium vapor box divertor localizes a dense cloud of lithium vapor away from the main plasma, in order to induce detachment at a stable, distant location. The vapor localization is created by local evaporation and nearby condensation, a configuration inaccessible for non-condensing seed impurities. This paper provides simulations of lithium vapor flow using the SPARTA direct simulation Monte Carlo (DSMC) code, in which the lithium vapor dynamics are governed by Li–Li collisions that should dominate in regions with little plasma. Inertia is included in this code, permitting the formation of shocks. The plasma model implemented in SPARTA is based on UEDGE results for the location of ionization and recombination. We find that a simplified lithium vapor box configuration, without baffles, provides robust stabilization of the detachment front and acceptable lithium vapor flow to the main chamber. Lithium is evaporated from a capillary porous surface on the private-flux side of the divertor, distant from the main plasma, and is condensed on the other side walls of the divertor chamber. The condensed lithium flows back to the evaporator through 2 cm internal diameter pipes. The capillary pressure differential at the surface of the evaporator is capable of overcoming MHD back-force in the piping, with a great deal of margin if a sandwich flow channel insert is implemented. Very little lithium should be accumulated on the first wall of the main chamber if it is maintained at 600 °C, as expected for a power-producing fusion system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Advanced assessment of beam-induced background at a muon collider

Renewed international interest in muon colliders motivates the continued investigation of the impacts of beam-induced background on detector performance. This continues the effort initiated by the Muon Accelerator Program and carried out until 2017. The beam-induced background from muon decays directly impacts detector performance and must be mitigated by optimizing the overall machine design, with particular attention paid to the machine detector interface region. In order to produce beam-induced background events and to study their characteristics in coordination with the collider optimization, a flexible simulation approach is needed. To achieve this goal we have chosen to utilize the combination of LineBuilder and Monte Carlo FLUKA codes. We report the results of beam-induced background studies with these tools obtained for a 1.5 TeV center of mass energy collider configuration. Good agreement with previous simulations using the MARS15 code demonstrates that our choice of tools meet the accuracy and performance requirements to perform future optimization studies on muon collider designs.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

ASAS J071404+7004.3 – a close, bright nova-like cataclysmic variable with gusty winds

ABSTRACT Despite being bright (V ≃ 11.8) and nearby (d = 212 pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton and the Hiltner and McGraw-Hill Telescopes, together with Swift X-ray and ultraviolet observations. We combined these with TESS photometry and find a period of 3.28 h and a mass transfer rate of $4\!-\!9\times 10^{-9}\, {\mathrm{M_{\odot }\, yr}^{-1}}$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo python code and shown that all the emission lines could emanate from the wind – which would explain the lack of double-peaked lines in such systems. In passing, we discuss the effect of variability on the position of cataclysmic variables in the Gaia Hertzsprung–Russell diagram.

Inight, K. (ORCID:0000000222002416)↗

DOE OSTI2022

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

ABSTRACT We implement euclidemulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the Markov chain Monte Carlo forecasting code montepython. We compare the performance of halofit, hmcode, and euclidemulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wavenumber cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wavenumber cut-off of $0.2\, h\, {\rm Mpc}^{-1}$. We find that euclidemulator1 is on average up to 17 per cent more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to 42 per cent and the equation of state of dark energy of up to 26 per cent, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5σ and 2σ (from below 1σ up to 6σ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of ≤1.4 in all cases.

Astronomy & Astrophysics↗

DOE OSTI2017

Total cross sections for ultracold neutrons scattered from gases

Here, we have followed up on our previous measurements of upscattering of ultracold neutrons (UCNs) from a series of gases by making measurements of total cross sections on the following gases hydrogen, ethane, methane, isobutene, n-butane, ethylene, water vapor, propane, neopentane, isopropyl alcohol, and 3 He. The values of these cross sections are important for estimating the loss rate of trapped neutrons due to residual gas and are relevant to neutron lifetime measurements using UCNs. The effects of the UCN velocity and path-length distributions were accounted for in the analysis using a Monte Carlo transport code. Results are compared to our previous measurements and with the known absorption cross section for 3 He scaled to our UCN energy. We find that the total cross sections for the hydrocarbon gases are reasonably described by a function linear in the number of hydrogen atoms in the molecule.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2012

Boron-10 Lined Proportional Counter Model Validation

The decreasing supply of 3He is stimulating a search for alternative neutron detectors; one potential 3He replacement is 10B-lined proportional counters. Simulations are being performed to predict the performance of systems designed with 10B-lined tubes. Boron-10-lined tubes are challenging to model accurately because the neutron capture material is not the same as the signal generating material. Thus, to simulate the efficiency, the neutron capture reaction products that escape the lining and enter the signal generating fill gas must be tracked. The tube lining thickness and composition are typically proprietary vendor information, and therefore add additional variables to the system simulation. The modeling methodologies used to predict the neutron detection efficiency of 10B-lined proportional counters were validated by comparing simulated to measured results. The measurements were made with a 252Cf source positioned at several distances from a moderated 2.54-cm diameter 10B-lined tube. Models were constructed of the experimental configurations using the Monte Carlo transport code MCNPX, which is capable of tracking the reaction products from the (n,10B) reaction. Several different lining thicknesses and compositions were simulated for comparison with the measured data. This paper presents the results of the evaluation of the experimental and simulated data, and a summary of how the different linings affect the performance of a coincidence counter configuration designed with 10B-lined proportional counters.

Lintereur, Azaree T.↗

DOE OSTI2014

SU-E-T-274: Radiation Therapy with Very High-Energy Electron (VHEE) Beams in the Presence of Metal Implants

Purpose: To evaluate the effect of metal implants on treatment plans for radiation therapy with very high-energy electron (VHEE) beams. Methods: The DOSXYZnrc/BEAMnrc Monte Carlo (MC) codes were used to simulate 50–150MeV VHEE beam dose deposition and its effects on steel and titanium (Ti) heterogeneities in a water phantom. Heterogeneities of thicknesses ranging from 0.5cm to 2cm were placed at 10cm depth. MC was also used to calculate electron and photon spectra generated by the VHEE beams' interaction with metal heterogeneities. The original VMAT patient dose calculation was planned in Eclipse. Patient dose calculations with MC-generated beamlets were planned using a Matlab GUI and research version of RayStation. VHEE MC treatment planning was performed on water-only geometry and water with segmented prostheses (steel and Ti) geometries with 100MeV and 150MeV beams. Results: 100MeV PDD 5cm behind steel/Ti heterogeneity was 51% less than in the water-only phantom. For some cases, dose enhancement lateral to the borders of the phantom increased the dose by up to 22% in steel and 18% in Ti heterogeneities. The dose immediately behind steel heterogeneity decreased by an average of 6%, although for 150MeV, the steel heterogeneity created a 23% increase in dose directly behind it. The average dose immediately behind Ti heterogeneities increased 10%. The prostate VHEE plans resulted in mean dose decrease to the bowel (20%), bladder (7%), and the urethra (5%) compared to the 15MV VMAT plan. The average dose to the body with prosthetic implants was 5% higher than to the body without implants. Conclusion: Based on MC simulations, metallic implants introduce dose perturbations to VHEE beams from lateral scatter and backscatter. However, when performing clinical planning on a prostate case, the use of multiple beams and inverse planning still produces VHEE plans that are dosimetrically superior to photon VMAT plans. BW Loo and P Maxim received research support from RaySearch laboratories; B Hardemark and E Hynning are employees of RaySearch.

60 APPLIED LIFE SCIENCES↗