Next Generation Uncertainty Quantification and Stochastic Media Monte Carlo Transport Methods.
Abstract not provided.
Olson, Aaron Jeffrey↗
SEARCH · Engineering Papers
Results for “Monte Carlo Methods”
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.
Next Generation Uncertainty Quantification and Stochastic Media Monte Carlo Transport Methods.
Abstract not provided.
Olson, Aaron Jeffrey↗
Uncertainty Quantification and Global Sensitivity Analysis Methods for Monte Carlo Radiation Transport Solvers
Explore the source record for details and available documents.
Clements, Kayla [Oregon State University]↗
Monte carlo solution of thermal transfer through radiant media between gray walls.
Radiative heat transfer to a nonisothermal absorbing and emitting gray gas between gray walls analyzed, using monte carlo method
GRAY GAS↗
Assessing the social sustainability contribution of an infrastructure project under conditions of uncertainty
Assessing the viability of a public infrastructure includes economic, technical and environmental aspects; however, on many occasions, the social aspects are not always adequately considered. This article proposes a procedure to estimate the social sustainability of infrastructure projects under conditions of uncertainty, based on a multicriteria deterministic method. The variability of the method inputs is contributed by the decision-makers. Uncertain inputs are treated through uniform and beta PERT distributions. The Monte Carlo method is used to propagate uncertainty in the method. A case study of a road infrastructure improvement in El Salvador is used to illustrate this treatment. The main results determine the variability of the short and long-term social improvement indices by infrastructure and the probability of the position in the prioritization of the alternatives. The proposed mechanism improves the reliability of the decision making early in infrastructure projects, taking their social contribution into account. The results can complement environmental and economic sustainability assessments. - Highlights: •Estimate the social sustainability of infrastructure projects under conditions of uncertainty •The method uses multicriteria and Monte Carlo techniques and beta PERT distributions •Determines variability of the short and long term social improvement •Determines probability in the prioritization of alternatives •Improves reliability of decision making considering the social contribution.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
GRID-BASED EXPLORATION OF COSMOLOGICAL PARAMETER SPACE WITH SNAKE
We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the 'curse of dimensionality' problem plaguing standard grid-based parameter estimation simply by disregarding grid cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings Markov Chain Monte Carlo methods include (1) trivial extraction of arbitrary conditional distributions; (2) direct access to Bayesian evidences; (3) better sampling of the tails of the distribution; and (4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N{sub par}. One of the main goals of the present paper is to determine how large N{sub par} can be, while still maintaining reasonable computational efficiency; we find that N{sub par} = 12 is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme.
79 ASTRONOMY AND ASTROPHYSICS↗
Numerical method for determination of resonance electron velocities in periodic electric fields
A Monte Carlo method is applied for modeling of a plasma source with a periodic structure of a radiofrequency electric field. The method is modified in order to detect resonance increase of the electron energy at given velocities of the electrons as expected from theoretical considerations. The numerical model is validated by comparison with known analytical results. The applicability of the method for detection of resonance velocities in plasma sources with electric field variation parallel or perpendicular to the electric field vector is proven.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A space radiation transport method development
Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.
NASA Center JSC↗
SU-F-207-03: Dosimetric Effect of the Position of Arms in Torso CT Scan with Tube Current Modulation
Purpose: To evaluate the patient organ dose differences between the arms-raised and arms-lowered postures in Torso multidetector computed tomography (MDCT) scan protocols with tube current modulation (TCM). Methods: Patient CT organ doses were simulated using the Monte Carlo method with human phantoms and a validated CT scanner model. A set of adult human phantoms with arms raised and arms lowered postures were developed using advanced BREP-based mesh surface geometries. Organ doses from routine Torso scan protocols such as chest, abdomen-pelvis, and CAP scans were simulated. The organ doses differences caused by two different posutres were investigated when tube current modulation (TCM) were applied during the CT scan. Results: With TCM applied, organ doses of all the listed organs of arms-lowered posture phantom are larger than those of arms raised phantom. The dose difference for most of the organs or tissues are larger than 50%, and the skin doses difference for abdomen-pelvis scan even reaches 112.03%. This is due to the fact that the tube current for patient with arms-lowered is much higher than for the arms raised posture. Conclusion: Considering CT scan with TCM, which is commonly applied clinically, patients who could not raise their arms will receive higher radiation dose than the arms raised patient, with dose differences for some tissues such as the skin being larger than 100%. This is due to the additional tube current necessary to penetrate the arms while maintaining consistent image quality. National Nature Science Foundation of China(No.11475047)
60 APPLIED LIFE SCIENCES↗
Testing of a Code for the Calculation of Spectra of Neutrons Produced in a Target of a Neutron Generator
The correctness of calculations performed with the SRIANG code for modeling the spectra of DT neutrons is estimated by comparing the obtained spectra to the results of calculations carried out with five different codes based on the Monte Carlo method.
43 PARTICLE ACCELERATORS↗
Monte Carlo solution for the characteristics of a highly rarefied ionized gas flowing through a channel with a transverse magnetic field.
Steady flow of highly rarefied ionized gas through channel with magnetic field solved by Monte Carlo method obtaining density, energies, wall shear stress, etc
CHANNEL FLOW↗
Analysis of Couette flow and heat transfer between parallel plates enclosing rarefied gas by Monte Carlo.
Couette flow and heat transfer of rarefied gas between parallel plates analyzed by Monte Carlo method
HEAT TRANSFER↗
Fe and Co nanostructures embedded into the Cu(100) surface: Self-Organization and magnetic properties
The self-organization and magnetic properties of small iron and cobalt nanostructures embedded into the first layer of a Cu(100) surface are investigated using the self-learning kinetic Monte Carlo method and density functional theory. The similarities and differences between the Fe/Cu(100) and the Co/Cu(100) are underlined. The time evolution of magnetic properties of a copper monolayer with embedded magnetic atoms at 380 K is discussed.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Calculating Control Rod Differential Worth Using Adjoint-Weighted Method
Recently, there are more and more research on using Monte Carlo method to calculate reactor important parameters. Control rod worth is one of the important parameters as it can determine the critical control rod position. There are two presentations of the control rod worth for utilization purpose, integral and differential. The former is the total reactivity change, while the latter is the reactivity change caused by per unit movement. In Monte Carlo calculation, differential control rod worth of a particular position is estimated by finite difference of reactivity between a small gap divided by gap distance. Though easy to implemented, the problem of this approach lies in that too large gap will enlarge error between 'differential' and 'difference', while too small gap will be disturbed by error caused by random fluctuation. Monte Carlo perturbation theory has been used to calculate the k{sub eff} sensitivity to certain parameters for a long time. By calculating the differential coefficient along with the k{sub eff}, this method is quite time-consuming and relatively accurate compared to finite difference method. Nevertheless, most of perturbation research focus on the perturbation of nuclide density or nuclear data. For geometrical perturbation, limited research was done before. Takahashi estimated the k{sub eff} sensitivity to reflector movement of pulsed fast reactor. MCNP6 used adjoint-weighted method for Godiva, Zeus and uranium oxyfluoride solution container benchmarks. This paper tries to apply the perturbation method to the determination of control rod differential worth. In this paper, based on RMC, a continuous-energy Reactor Monte Carlo transport code, Iterated Fission Probability method (IFP) is researched to calculate control rod differential worth on assembly level. IFP was originally developed by Hurwitz, and has been adopted by most of Monte Carlo code to compute adjoint-weighted parameters. This paper is organized with four parts. The first is introduction. Adjoint-weighted method is introduced in the second part. Problem 4 of benchmark from CASL is built to verify the method. Conclusion is given in final part.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Transport error estimation using residual Monte Carlo
The residual Monte Carlo (RMC) method is also known in the literature as sequential Monte Carlo and reduced-source Monte Carlo. Given a Monte Carlo method for solving a linear equation and an approximate solution to that system, the residual method enables use of essentially the same Monte Carlo algorithm to directly compute the additive error or “defect” associated with the approximate solution. As the size of the defect decreases relative to the size of the solution, the residual Monte Carlo method becomes increasingly efficient relative to the standard Monte Carlo (SMC) method. Here we present a new RMC algorithm for evaluating the space-angle error in S n radiation transport solutions, and provide computational examples demonstrating that it can be far more efficient than SMC for this purpose. Herein we also describe a particular pitfall that must be avoided if RMC is to be efficient, and explain why the performance of RMC can significantly differ between different transport problems and different quantities of interest for the same problem.
97 MATHEMATICS AND COMPUTING↗
Tungsten anode spectral model using interpolating cubic splines: Unfiltered x-ray spectra from 20 kV to 640 kV
Purpose: Monte Carlo methods were used to generate lightly filtered high resolution x-ray spectra spanning from 20 kV to 640 kV. Methods: X-ray spectra were simulated for a conventional tungsten anode. The Monte Carlo N-Particle eXtended radiation transport code (MCNPX 2.6.0) was used to produce 35 spectra over the tube potential range from 20 kV to 640 kV, and cubic spline interpolation procedures were used to create piecewise polynomials characterizing the photon fluence per energy bin as a function of x-ray tube potential. Using these basis spectra and the cubic spline interpolation, 621 spectra were generated at 1 kV intervals from 20 to 640 kV. The tungsten anode spectral model using interpolating cubic splines (TASMICS) produces minimally filtered (0.8 mm Be) x-ray spectra with 1 keV energy resolution. The TASMICS spectra were compared mathematically with other, previously reported spectra. Results: Using pairedt-test analyses, no statistically significant difference (i.e., p > 0.05) was observed between compared spectra over energy bins above 1% of peak bremsstrahlung fluence. For all energy bins, the correlation of determination (R{sup 2}) demonstrated good correlation for all spectral comparisons. The mean overall difference (MOD) and mean absolute difference (MAD) were computed over energy bins (above 1% of peak bremsstrahlung fluence) and over all the kV permutations compared. MOD and MAD comparisons with previously reported spectra were 2.7% and 9.7%, respectively (TASMIP), 0.1% and 12.0%, respectively [R. Birch and M. Marshall, “Computation of bremsstrahlung x-ray spectra and comparison with spectra measured with a Ge(Li) detector,” Phys. Med. Biol. 24, 505–517 (1979)], 0.4% and 8.1%, respectively (Poludniowski), and 0.4% and 8.1%, respectively (AAPM TG 195). The effective energy of TASMICS spectra with 2.5 mm of added Al filtration ranged from 17 keV (at 20 kV) to 138 keV (at 640 kV); with 0.2 mm of added Cu filtration the effective energy was 9 keV at 20 kV and 169 keV at 640 kV. Conclusions: Ranging from 20 kV to 640 kV, 621 x-ray spectra were produced and are available at 1 kV tube potential intervals. The spectra are tabulated at 1 keV intervals. TASMICS spectra were shown to be largely equivalent to published spectral models and are available in spreadsheet format for interested users by emailing the corresponding author (JMB)
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Correction of Fluorescent Images of Biological Tissues for the Effects of Light Absorption and Scattering
A method for correcting the influence of light attenuation processes in biological tissues on their fluorescent images is proposed. The transfer function that takes into account the radiation losses in the medium at the wavelengths of excitation and emission of fluorescence is calculated analytically, depending on the transport scattering index and the hemoglobin tissue index. The latter is determined on the basis of a color image of the tissue in reflected visible light. Verification of the developed method was carried out on the basis of computer simulation of fluorescent images of biological tissue by the Monte Carlo method.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Numerical thermalization in particle-in-cell simulations with Monte-Carlo collisions
Numerical thermalization in collisional one-dimensional (1D) electrostatic (ES) particle-in-cell (PIC) simulations was investigated. Two collision models, the pitch-angle scattering of electrons by the stationary ion background and large-angle collisions between the electrons and the neutral background, were included in the PIC simulation using Monte-Carlo methods. The numerical results show that the thermalization times in both models were considerably reduced by the additional Monte-Carlo collisions as demonstrated by comparisons with Turner's previous simulation results based on a head-on collision model [M. M. Turner, Phys. Plasmas 13, 033506 (2006)]. However, the breakdown of Dawson's scaling law in the collisional 1D ES PIC simulation is more complicated than that was observed by Turner, and the revised scaling law of the numerical thermalization time with numerical parameters are derived on the basis of the simulation results obtained in this study.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗