Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Lunar Surface Reactor Shielding Study

Nuclear reactor system could provide power to support a long term human exploration to the moon. Such a system would require shielding to protect astronauts from its emitted radiations. Shielding studies have been performed for a Gas Cooled Reactor (GCR) system because it is considered to be the most suitable nuclear reactor system available for lunar exploration, based on its tolerance of oxidizing lunar regolith and its good conversion efficiency (Wright, 2003). The goals of the shielding studies were to provide optimal material shielding configuration that reduces the dose (rem) to the required level in order to protect astronauts, and to estimate the mass of regolith that would provide an equivalent protective effect if it were used as the shielding material. All calculations were performed using MCNPX code, a Monte Carlo transport code.

Monte Carlo N-Particle eXtended (MCNPX)↗
NASA NTRS2008

Neutron Transport Models and Methods for HZETRN and Coupling to Low Energy Light Ion Transport

Exposure estimates inside space vehicles, surface habitats, and high altitude aircraft exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETCHEDS and FLUKA, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light ion (A<4) transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Blattnig, S.R.↗
NASA NTRS2009

Coupled Neutron Transport for HZETRN

Exposure estimates inside space vehicles, surface habitats, and high altitude aircrafts exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETC-HEDS, FLUKA, and MCNPX, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light particle transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Slaba, Tony C.↗
NASA NTRS1983

Computer simulation study of multiple germanium gamma-ray sensor arrays

Design studies of large high resolution gamma-ray spectrometer systems employing multiple sensors are greatly facilitated by the use of computer simulation techniques. A study has been made by interfacing a detector geometry code with the EGS-3 version of the Stanford Linear Accelerator Center Monte-Carlo Code which has been modified to run on our VAX computer. The detector responses have been simulated for input photons with energies between 0.1 and 10 MeV. Examples are given to illustrate the performance of complex sensor arrays with different configurations of the anticoincidence system. Of particular interest is a design composed of an array of 7 detector elements each consisting of a thin planar Ge (1 cm thick) detector mounted in front of and in tandem with a large coaxial Ge detector.

Nakano, G. H.↗
NASA NTRS1989

BRYNTRN: A baryon transport model

The development of an interaction data base and a numerical solution to the transport of baryons through an arbitrary shield material based on a straight ahead approximation of the Boltzmann equation are described. The code is most accurate for continuous energy boundary values, but gives reasonable results for discrete spectra at the boundary using even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O). The resulting computer code is self-contained, efficient and ready to use. The code requires only a very small fraction of the computer resources required for Monte Carlo codes.

Wilson, John W.↗
NASA NTRS1989

Minimum length Pb/SCIN detector for efficient cosmic ray identification

A study was made of the performance of a minimal length cosmic ray shower detector that would be light enough for space flight and would provide efficient identification of positions and protons. Cosmic ray positions are mainly produced in the decay chain of: Pion yields Muon yields Positron and they provide a measure of the matter density traversed by primary protons. Present positron flux measurements are consistent with the Leaky Box and Halo models for sources of cosmic rays. Abundant protons in the space environment are a significant source of background that would wash out the positron signal. Protons and positrons produced very distictive showers of particles when they enter matter; many studies have been published on their behavior on large calorimeter detectors. The challenge is to determine the minimal material necessary (minimal calorimeter depth) for positive particles identification. The primary instrument for the investigation is the Monte Carlo code GEANT, a library of programs from CERN that can be used to model experimental geometry, detector responses and particle interaction processes. The use of the Monte Carlo approach is crucial since statistical fluctuations in shower shape are significant. Studies conducted during the 1988 summer program showed that straightforward approaches to the problem achieved 85 to 90 percent correct identification, but left a residue of 10 to 15 percent misidentified particles. This percentage improved to a few percent when multiple shower-cut criteria were applied to the data. This summer, the same study was extended to employ several physical and statistical methods of identifying response of the calorimeter and the efficiency of the optimal shower cuts to off-normal incidence particle was determined.

Snyder, H. David↗
NASA NTRS2017

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗
NASA NTRS1991

ASRM radiation and flowfield prediction status

Existing and proposed methods for the prediction of plume radiation are discussed in terms of their application to the NASA Advanced Solid Rocket Motor (ASRM) and Space Shuttle Main Engine (SSME) projects. Extrapolations of the Solid Rocket Motor (SRM) are discussed with respect to preliminary predictions of the primary and secondary radiation environments. The methodology for radiation and initial plume property predictions are set forth, including a new code for scattering media and independent secondary source models based on flight data. The Monte Carlo code employs a reverse-evaluation approach which traces rays back to their point of absorption in the plume. The SRM sea-level plume model is modified to account for the increased radiation in the ASRM plume due to the ASRM's propellant chemistry. The ASRM cycle-1 environment predictions are shown to identify a potential reason for the shutdown spike identified with pre-SRM staging.

Reardon, J. E.↗
NASA NTRS1992

Model scattering envelopes of young stellar objects. I - Method and application to circumstellar disks

We describe a Monte Carlo code that accurately treats multiple scattering, absorption, and polarization by dust, and use this code to calculate images of dusty disks around young stellar objects. We present some approximate analytic results that describe the behavior of the Monte Carlo calculations. A geometrically thin disk illuminated by a central T Tauri star scatters very little light at distances of many AU from the star. Viewed at any inclination, the flux scattered by such a disk at the distance to the nearest star-forming region will be overwhelmed by the stellar image. An optically thick disk that has a flaring surface may be observable, especially if viewed nearly edge-on so that the stellar source becomes occulted. An optically thin disk with a finite opening angle, similar to the one surrounding beta Pictoris, is about as observable as the typical flared optically thick disk at a similar distance from the earth. The polarization position angle is perpendicular to the disk plane in all of the models, in contrast to observations of many young stellar objects which have the position angle oriented parallel to the presumed disk plane. We suggest that the scattered light structures observed around many premain-sequence objects are dusty envelopes rather than disks.

Whitney, Barbara A.↗
NASA NTRS1989

Monte Carlo simulation of particle acceleration at astrophysical shocks

A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.

Campbell, Roy K.↗
NASA NTRS2004

High Altitude Radiations Relevant to the High Speed Civil Transport (HSCT)

The Langley Research Center (LaRC) performed atmospheric radiation studies under the SST development program in which important ionizing radiation components were measured and extended by calculations to develop the existing atmospheric ionizing radiation (AIR) model. In that program the measured neutron spectrum was limited to less than 10 MeV by the available 1960-1970 instrumentation. Extension of the neutron spectrum to high energies was made using the LaRC PROPER-3C monte carlo code. It was found that the atmospheric neutrons contributed about half of the dose equivalent and approximately half of the neutron contribution was from high energy neutrons above 10 MeV. Furthermore, monte carlo calculations of solar particle events showed that potential exposures as large as 10-100 mSv/hr may occur on important high latitude routes but acceptable levels of exposure could be obtained if timely descent to subsonic altitudes could be made. The principal concern was for pregnant occupants onboard the aircraft. As a result of these studies the FAA Advisory Committee on the Radiobiological Aspects of the SST recommended: 1. Crew members will have to be informed of their exposure levels 2. Maximum exposures on any flight to be limited to 5 mSv 3. Airborne radiation detection devices for total exposure and exposure rates 4. Satellite monitoring system to provide SST aircraft real-time information on atmospheric radiation levels for exposure mitigation 5. A solar forecasting system to warn flight operations of an impending solar event for flight scheduling and alert status. These recommendations are a reasonable starting point to requirements for the HSCT with some modification reflecting new standards of protection as a result of changing risk coefficients.

Wilson, J. W.↗
NASA NTRS2013

Development of the Code RITRACKS

A document discusses the code RITRACKS (Relativistic Ion Tracks), which was developed to simulate heavy ion track structure at the microscopic and nanoscopic scales. It is a Monte-Carlo code that simulates the production of radiolytic species in water, event-by-event, and which may be used to simulate tracks and also to calculate dose in targets and voxels of different sizes. The dose deposited by the radiation can be calculated in nanovolumes (voxels). RITRACKS allows simulation of radiation tracks without the need of extensive knowledge of computer programming or Monte-Carlo simulations. It is installed as a regular application on Windows systems. The main input parameters entered by the user are the type and energy of the ion, the length and size of the irradiated volume, the number of ions impacting the volume, and the number of histories. The simulation can be started after the input parameters are entered in the GUI. The number of each kind of interactions for each track is shown in the result details window. The tracks can be visualized in 3D after the simulation is complete. It is also possible to see the time evolution of the tracks and zoom on specific parts of the tracks. The software RITRACKS can be very useful for radiation scientists to investigate various problems in the fields of radiation physics, radiation chemistry, and radiation biology. For example, it can be used to simulate electron ejection experiments (radiation physics).

Plante, Ianik↗
NASA NTRS2011

Comparison of DSMC and CFD Solutions of Fire II Including Radiative Heating

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. These flows may also contain significant radiative heating. To prepare for these missions, NASA is developing the capability to simulate rarefied, ionized flows and to then calculate the resulting radiative heating to the vehicle's surface. In this study, the DSMC codes DAC and DS2V are used to obtain charge-neutral ionization solutions. NASA s direct simulation Monte Carlo code DAC is currently being updated to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced Quantum-Kinetic chemistry model, and to include electronic energy levels as an additional internal energy mode. The Fire II flight test is used in this study to assess these new capabilities. The 1634 second data point was chosen for comparisons to be made in order to include comparisons to computational fluid dynamics solutions. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid. It is shown that there can be quite a bit of variability in the vibrational temperature inferred from DSMC solutions and that, from how radiative heating is computed, the electronic temperature is much better suited for radiative calculations. To include the radiative portion of heating, the flow-field solutions are post-processed by the non-equilibrium radiation code HARA. Acceptable agreement between CFD and DSMC flow field solutions is demonstrated and the progress of the updates to DAC, along with an appropriate radiative heating solution, are discussed. In addition, future plans to generate more high fidelity radiative heat transfer solutions are discussed.

Liechty, Derek S.↗
NASA NTRS2009

Comparison of Radiation Transport Codes, HZETRN, HETC and FLUKA, Using the 1956 Webber SPE Spectrum

Protection of astronauts and instrumentation from galactic cosmic rays (GCR) and solar particle events (SPE) in the harsh environment of space is of prime importance in the design of personal shielding, spacec raft, and mission planning. Early entry of radiation constraints into the design process enables optimal shielding strategies, but demands efficient and accurate tools that can be used by design engineers in every phase of an evolving space project. The radiation transport code , HZETRN, is an efficient tool for analyzing the shielding effectiveness of materials exposed to space radiation. In this paper, HZETRN is compared to the Monte Carlo codes HETC-HEDS and FLUKA, for a shield/target configuration comprised of a 20 g/sq cm Aluminum slab in front of a 30 g/cm^2 slab of water exposed to the February 1956 SPE, as mode led by the Webber spectrum. Neutron and proton fluence spectra, as well as dose and dose equivalent values, are compared at various depths in the water target. This study shows that there are many regions where HZETRN agrees with both HETC-HEDS and FLUKA for this shield/target configuration and the SPE environment. However, there are also regions where there are appreciable differences between the three computer c odes.

Heinbockel, John H.↗
NASA NTRS2000

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗
NASA NTRS1993

Gravitational microlensing of gamma-ray bursts

A Monte Carlo code is developed to calculate gravitational microlensing in three dimensions when the lensing optical depth is low or moderate (not greater than 0.25). The code calculates positions of microimages and time delays between the microimages. The majority of lensed gamma-ray bursts should show a simple double-burst structure, as predicted by a single point mass lens model. A small fraction should show complicated multiple events due to the collective effects of several point masses (black holes). Cosmological models with a significant fraction of mass density in massive compact objects can be tested by searching for microlensing events in the current BATSE data. Our catalog generated by 10,000 Monte Carlo models is accessible through the computer network. The catalog can be used to take realistic selection effects into account.

Mao, Shude↗
NASA NTRS2017

3DHZETRN: Inhomogeneous Geometry Issues

Historical methods for assessing radiation exposure inside complicated geometries for space applications were limited by computational constraints and lack of knowledge associated with nuclear processes occurring over a broad range of particles and energies. Various methods were developed and utilized to simplify geometric representations and enable coupling with simplified but efficient particle transport codes. Recent transport code development efforts, leading to 3DHZETRN, now enable such approximate methods to be carefully assessed to determine if past exposure analyses and validation efforts based on those approximate methods need to be revisited. In this work, historical methods of representing inhomogeneous spacecraft geometry for radiation protection analysis are first reviewed. Two inhomogeneous geometry cases, previously studied with 3DHZETRN and Monte Carlo codes, are considered with various levels of geometric approximation. Fluence, dose, and dose equivalent values are computed in all cases and compared. It is found that although these historical geometry approximations can induce large errors in neutron fluences up to 100 MeV, errors on dose and dose equivalent are modest (<10%) for the cases studied here.

Wilson, John W.↗