SEARCH · Search NASA

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,711 records · Page 95

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Effects of alpha stopping power modelling on the ignition threshold in a directly-driven inertial confinement fusion capsule

The alpha-particle energy deposition mechanism modifies the ignition conditions of the thermonuclear Deuterium-Tritium fusion reactions, and constitutes a key issue in achieving high gain in Inertial Confinement Fusion implosions. One-dimensional hydrodynamic calculations have been performed with the code Multi-IFE to simulate the implosion of a capsule directly irradiated by a laser beam. The diffusion approximation for the alpha energy deposition has been used to optimize three laser profiles corresponding to different implosion velocities. A Monte-Carlo package has been included in Multi-IFE to calculate the alpha energy transport, and in this case the energy deposition uses both the LP and the BPS stopping power models. Homothetic transformations that maintain a constant implosion velocity have been used to map out the transition region between marginally-igniting and high-gain configurations. Furthermore, the results provided by the two models have been compared and it is found that – close to the ignition threshold – in order to produce the same fusion energy, the calculations performed with the BPS model require about 10% more invested energy with respect to the LP model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

BIhNNs

The code enables to perform Bayesian inference in an efficient manner through the use of Hamiltonian Neural Networks (HNNs), Deep Neural Networks (DNNs), Neural ODEs, and Symplectic Neural Networks (SympNets) used with state-of-the-art sampling schemes like Hamiltonian Monte Carlo (HMC) and the No-U-Turn-Sampler (NUTS).

Dhulipala, Som↗

NASA NTRS1982

Computer program for prediction of capture maneuver probability for an on-off reaction controlled upper stage

A FORTRAN coded computer program which computes the capture transient of a launch vehicle upper stage at the ignition and/or separation event is presented. It is for a single degree-of-freedom on-off reaction jet attitude control system. The Monte Carlo method is used to determine the statistical value of key parameters at the outcome of the event. Aerodynamic and booster induced disturbances, vehicle and control system characteristics, and initial conditions are treated as random variables. By appropriate selection of input data pitch, yaw and roll axes can be analyzed. Transient response of a single deterministic case can be computed. The program is currently set up on a CDC CYBER 175 computer system but is compatible with ANSI FORTRAN computer language. This routine has been used over the past fifteen (15) years for the SCOUT Launch Vehicle and has been run on RECOMP III, IBM 7090, IBM 360/370, CDC6600 and CDC CYBER 175 computers with little modification.

Knauber, R. N.↗

NASA NTRS1995

DS-CDMA satellite diversity reception for personal satellite communication: Downlink performance analysis

The downlink of a satellite-mobile personal communication system employing power-controlled Direct Sequence Code Division Multiple Access (DS-CDMA) and exploiting satellite-diversity is analyzed and its performance compared with a more traditional communication system utilizing single satellite reception. The analytical model developed has been thoroughly validated by means of extensive Monte Carlo computer simulations. It is shown how the capacity gain provided by diversity reception shrinks considerably in the presence of increasing traffic or in the case of light shadowing conditions. Moreover, the quantitative results tend to indicate that to combat system capacity reduction due to intra-system interference, no more than two satellites shall be active over the same region. To achieve higher system capacity, differently from terrestrial cellular systems, Multi-User Detection (MUD) techniques are likely to be required in the mobile user terminal, thus considerably increasing its complexity.

DeGaudenzi, Riccardo↗

NASA NTRS2007

Earth Global Reference Atmospheric Model (GRAM99): Short Course

Earth-GRAM is a FORTRAN software package that can run on a variety of platforms including PC's. For any time and location in the Earth's atmosphere, Earth-GRAM provides values of atmospheric quantities such as temperature, pressure, density, winds, constituents, etc.. Dispersions (perturbations) of these parameters are also provided and have realistic correlations, means, and variances - useful for Monte Carlo analysis. Earth-GRAM is driven by observations including a tropospheric database available from the National Climatic Data Center. Although Earth-GRAM can be run in a "stand-alone" mode, many users incorporate it into their trajectory codes. The source code is distributed free-of-charge to eligible recipients.

Leslie, Fred W.↗

NASA NTRS2011

Trajectory Guidance for Mars Robotic Precursors: Aerocapture, Entry, Descent, and Landing

Future crewed missions to Mars require improvements in landed mass capability beyond that which is possible using state-of-the-art Mars Entry, Descent, and Landing (EDL) systems. Current systems are capable of an estimated maximum landed mass of 1-1.5 metric tons (MT), while human Mars studies require 20-40 MT. A set of technologies were investigated by the EDL Systems Analysis (SA) project to assess the performance of candidate EDL architectures. A single architecture was selected for the design of a robotic precursor mission, entitled Exploration Feed Forward (EFF), whose objective is to demonstrate these technologies. In particular, inflatable aerodynamic decelerators (IADs) and supersonic retro-propulsion (SRP) have been shown to have the greatest mass benefit and extensibility to future exploration missions. In order to evaluate these technologies and develop the mission, candidate guidance algorithms have been coded into the simulation for the purposes of studying system performance. These guidance algorithms include aerocapture, entry, and powered descent. The performance of the algorithms for each of these phases in the presence of dispersions has been assessed using a Monte Carlo technique.

Sostaric, Ronald R.↗

NASA NTRS2017

Cyclotron Resonant Scattering Feature Simulations II. Description of the CRSF Simulation Process

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons o_ quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems.Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time.Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua.Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗

DOE OSTI2017

Sensitivity and uncertainty analysis for a combined depletion and criticality calculation chain of a BWR FA

This paper examines the prerequisites and boundary conditions for extending burnup credit from pressurized water reactors to boiling water reactors with validated depletion codes by performing uncertainty and sensitivity analyses on a modern boiling water reactor fuel assembly design. We developed a computational chain consisting of an axial boiling water reactor axial irradiation model, followed by a series of depletion calculations and a three-dimensional criticality calculation. The irradiation model puts into relation the profiles of power, burnup and moderator density of the fuel assembly at 24 axial locations. In a simplified fashion it also includes the insertion of control blades. The depletion calculations of the fuel assembly were performed in axial slices with the 2D-depletion code HELIOS 1.12. The criticality calculations were done in 3D-geometry with KENO-Va from the SCALE 6.2.1 package, in which several fuel assemblies are placed into a generic transport cask. For the uncertainty and sensitivity analyses a Monte Carlo sampling based approach is used applying the GRS tool SUnCISTT. We performed several generic investigations testing the capabilities and the limitations of the whole computational chain. The impact of variations of the irradiation model on the resulting reactivity of the fuel assemblies inside a transport cask or storage pool was investigated using an uncertainty and sensitivity analysis on the whole computational chain. Geometrical uncertainties were also included in the evaluation and a parameter study on the use of control blades during irradiation was performed. The goal was to develop the procedures and tools to be able to obtain quantitative uncertainty data which can be used in the future to derive conservative irradiation profiles without unnecessary excess reactivity margins. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

SU-G-201-02: Application of RayStretch in Clinical Cases: A Calculation for Heterogeneity Corrections in LDR Permanent I-125 Prostate Brachytherapy

Purpose: Tissue heterogeneities and calcifications have significant impact on the dosimetry of low energy brachytherapy (BT). RayStretch is an analytical algorithm developed in our institution to incorporate heterogeneity corrections in LDR prostate brachytherapy. The aim of this work is to study its application in clinical cases by comparing its predictions with the results obtained with TG-43 and Monte Carlo (MC) simulations. Methods: A clinical implant (71 I-125 seeds, 15 needles) from a real patient was considered. On this patient, different volumes with calcifications were considered. Its properties were evaluated in three ways by i) the Treatment planning system (TPS) (TG-43), ii) a MC study using the Penelope2009 code, and iii) RayStretch. To analyse the performance of RayStretch, calcifications located in the prostate lobules covering 11% of the total prostate volume and larger calcifications located in the lobules and underneath the urethra for a total occupied volume of 30% were considered. Three mass densities (1.05, 1.20, and 1.35 g/cm3) were explored for the calcifications. Therefore, 6 different scenarios ranging from small low density calcifications to large high density ones have been discussed. Results: DVH and D90 results given by RayStretch agree within 1% with the full MC simulations. Although no effort has been done to improve RayStretch numerical performance, its present implementation is able to evaluate a clinical implant in a few seconds to the same level of accuracy as a detailed MC calculation. Conclusion: RayStretch is a robust method for heterogeneity corrections in prostate BT supported on TG-43 data. Its compatibility with commercial TPSs and its high calculation speed makes it feasible for use in clinical settings for improving treatment quality. It will allow in a second phase of this project, its use during intraoperative ultrasound planning. This study was partly supported by a fellowship grant from the Spanish Ministry of Education, by the Generalitat Valenciana under Project PROMETEOII/2013/010, by the Spanish Government under Project No. FIS2013-42156 and by the European Commission within the SeventhFramework Program through ENTERVISION (grant agreement number 264552).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Evaluating the X2 initial core zero power physics tests with Serpent-Ants

The validation of the Ants nodal neutronics code for VVER applications is started by modelling the zero power physics tests for the initial core of the Khmelnitsky 2 nuclear power plant as described in the X2 benchmark using the Serpent-Ants two step neutronics calculation chain. The Ants prediction compare favorably against an earlier continuous energy Monte-Carlo reference solution as well as the measured data, except for discrepancies in the SCRAM worth between predicted and measured values. Similar discrepancies have been previously reported for VVER-1000 reactors, and the Ants predicted SCRAM worths match well with the Serpent results. Additionally, the assembly power distribution and axial power distribution predicted by Ants for the critical HZP state of the reactor are also compared to a reference Serpent solution with very good results. Lastly, an investigation into the effects of the few-group structure used in the nodal calculations and the use or lack of leakage correction for group constants on the results shows that the best accuracy is reached with seven or eight energy groups without leakage correction, although a reasonable accuracy can also be obtained with two or three energy groups using fundamental mode leakage correction. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1999

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.

Cahalan, Robert F.↗

NASA NTRS2001

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

NASA NTRS2001

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

NASA NTRS2016

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Modeling↗

DOE OSTI2016

SU-F-T-364: Monte Carlo-Dose Verification of Volumetric Modulated Arc Therapy Plans Using AAPM TG-119 Test Patterns

Purpose: To investigate the Monte Carlo (MC)-based dose verification for VMAT plans by a treatment planning system (TPS). Methods: The AAPM TG-119 test structure set was used for VMAT plans by the Pinnacle3 (convolution/superposition), using a Synergy radiation head of a 6 MV beam with the Agility MLC. The Synergy was simulated with the EGSnrc/BEAMnrc code, and VMAT dose distributions were calculated with the EGSnrc/DOSXYZnrc code by the same irradiation conditions as TPS. VMAT dose distributions of TPS and MC were compared with those of EBT3 film, by 2-D gamma analysis of ±3%/3 mm criteria with a threshold of 30% of prescribed doses. VMAT dose distributions between TPS and MC were also compared by DVHs and 3-D gamma analysis of ±3%/3 mm criteria with a threshold of 10%, and 3-D passing rates for PTVs and OARs were analyzed. Results: TPS dose distributions differed from those of film, especially for Head & neck. The dose difference between TPS and film results from calculation accuracy for complex motion of MLCs like tongue and groove effect. In contrast, MC dose distributions were in good agreement with those of film. This is because MC can model fully the MLC configuration and accurately reproduce the MLC motion between control points in VMAT plans. D95 of PTV for Prostate, Head & neck, C-shaped, and Multi Target was 97.2%, 98.1%, 101.6%, and 99.7% for TPS and 95.7%, 96.0%, 100.6%, and 99.1% for MC, respectively. Similarly, 3-D gamma passing rates of each PTV for TPS vs. MC were 100%, 89.5%, 99.7%, and 100%, respectively. 3-D passing rates of TPS reduced for complex VMAT fields like Head & neck because MLCs are not modeled completely for TPS. Conclusion: MC-calculated VMAT dose distributions is useful for the 3-D dose verification of VMAT plans by TPS.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction

Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗