SEARCH · Engineering Papers
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.
A total dose analysis for exposure of electronic parts at the box level is presented for the Hubble Space Telescope. This was done using solid angle sectoring/3-dimensional ray trace and Monte Carlo radiation transport simulations. Results are discussed in terms of parts that are potential total dose concerns.
GEANT 4 is a Monte Carlo tool set developed by the High Energy Physics Community (CERN, SLAC, etc) to perform simulations of complex particle detectors. GEANT4 is the ideal tool to study radiation transport and should be applied to space environments and the complex geometries of modern day spacecraft.
Explore the source record for details and available documents.
Computations of cosmic ray transport based upon finite difference methods are afflicted by instabilities, inaccuracies, and artifacts. To avoid these problems, a Monte Carlo formulation was developed which is closely related not only to the finite difference formulation, but also to the underlying physics of transport phenomena. Implementations of this approach are currently running on the Massively Parallel Processor at Goddard, whose enormous computing power overwhelms the poor statistical accuracy that usually limits the use of stochastic methods. In a Monte Carlo simulation of rectilinear transport, the coherent and diffusive effects that appeared are in good quantitative agreement with both finite difference and analytical calculations.
The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.
A Monte Carlo simulation of ultrarelativistic electron transport in solar flare magnetic loops has been developed. It includes Coulomb, synchrotron, and bremsstrahlung energy losses; pitch-angle scattering by Alfven and whistler turbulence in the coronal region of the loop; and magnetic mirroring in the converging magnetic flux tubes beneath the transition region. Depth distributions, time profiles, energy spectra, and angular distributions of the resulting bremsstrahlung emission are calculated. It is found that both the preferential detection of solar flares with greater than 10 MeV emission near the limb of the sun and the observation of ultrarelativistic electron bremsstrahlung from flares on the disk are consequences of the loop transport model. The declining portions of the observed time profiles of greater than 10 MeV emission from solar flares can also be accounted for, and it is proposed that these portions are determined by transport and not acceleration.
Recent models for the ion and relativistic electron transport in solar flare magnetic loops are reviewed, focusing on the ways in which the models handle loop structure, particle acceleration, and particle transport. Results are presented from Monte Carlo simulations of ion and relativistic ion transport. Consideration is given to results concerning the depth distributions of gamma-ray production, the attenuation of gamma-ray lines, the time dependences of the various emissions and the angular distribution of the bremsstrahlung.
Monte Carlo radiation transport codes such as RITRACKS or Geant4 are used to simulate the interaction of ions with matter. These codes rely on sampling algorithms to determine interactions and various physical properties of particles involved in the simulations. It is crucial to develop efficient sampling algorithms since Monte Carlo radiation transport simulations can be time consuming. This work presents an efficient sampling algorithm to determine the energy of secondary electrons following ion-water interactions. The applicability and intended use of the algorithm are discussed in detail, and it is shown that the new algorithm is up to 6X10 4 times faster than the method currently used in Geant4-DNA.
When charged particles spiral along a large constant magnetic field, their trajectories are scattered by any random field components that are superposed on the guiding field. If the random field configuration embodies helicity, the scattering is asymmetrical with respect to a plane perpendicular to the guiding field, for particles moving into the forward hemisphere are scattered at different rates from those moving into the backward hemisphere. This asymmetry gives rise to new terms in the transport equations that describe propagation of charged particles. Helicity has virtually no impact on qualitative features of the diffusive mode of propagation. However, characteristic velocities of the coherent modes that appear after a highly anisotropic injection exhibit an asymmetry related to helicity. Explicit formulas, which embody the effects of helicity, are given for the anisotropies, the coefficient diffusion, and the coherent velocities. Predictions derived from these expressions are in good agreement with Monte Carlo simulations of particle transport, but the simulations reveal certain phenomena whose explanation calls for further analytical work.
A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.
Forest fires can burn large areas, but can also inject smoke into the upper troposphere/lower stratosphere (UT/LS), where stakes are even higher for climate, because emissions tend to have a longer lifetime, and can produce significant regional and even global climate effects, as is the case with some volcanoes. Large forest fires are now believed to be more common in summer, especially in the boreal regions, where pyrocumulus (pyroCu), and occasionally pyrocumuionimbus (pyroCb) clouds are formed, which can transport emissions into the UT/LS. A major difficulty in developing realistic fire plume models is the lack of observational data within fire plumes that resolves structure at a few 100 m scales, which can be used to validate these models. Here, we report detailed airborne radiation measurements within strong pyroCu taken over boreal forest fires in Saskatchewan, Canada during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) summer field campaign in 2008. We find that the angular distribution of radiance within the pyroCu is closely related to the diffusion domain in water clouds and can be described by very similar simple cosine functions. We demonstrate with Monte Carlo simulations that radiation transport in pyroCu is inherently a 3D phenomenon and must account for particle absorption. However, the simple cosine function promises to offer an easy solution for climate models. The presence of a prominent smoke core, defined by strong extinction in the UV, VIS and NIR, suggests that the core might be an important pathway for emission transport to the upper troposphere and lower stratosphere. We speculate that this plume injection core is generated and sustained by complex processes not yet well understood, but not necessarily related directly to the intense fires that originally initiated the plume rise.
The potential of thermoluminescence measurements of bricks from the contaminated area of the Techa river valley, Southern Urals, Russia, for reconstructing external exposures of affected population groups has been studied. Thermoluminescence dating of background samples was used to evaluate the age of old buildings available on the river banks. The anthropogenic gamma dose accrued in exposed samples is determined by subtracting the natural radiation background dose for the corresponding age from the accumulated dose measured by thermoluminescence. For a site in the upper Techa river region, where the levels of external exposures were extremely high, the depth-dose distribution in bricks and the dependence of accidental dose on the height of the sampling position were determined. For the same site, Monte Carlo simulations of radiation transport were performed for different source configurations corresponding to the situation before and after the construction of a reservoir on the river and evacuation of the population in 1956. A comparison of the results provides an understanding of the features of the measured depth-dose distributions and height dependencies in terms of the source configurations and shows that bricks from the higher sampling positions are likely to have accrued a larger fraction of anthropogenic dose from the time before the construction of the reservoir. The applicability of the thermoluminescent dosimetry method to environmental dose reconstruction in the middle Techa region, where the external exposure was relatively low, was also investigated.
Monte Carlo transport program for simulating penetration and interaction of nucleons with absorbing matter
Various papers on scatter from optical components are presented. Individual topics addressed include: BRDF of SiC and Al foam compared to black paint at 3.39 microns, characterization of optical baffle materials, bidirectional transmittance distribution function of several IR materials at 3.39 microns, thermal cycling effects on the BRDF of beryllium mirrors, BTDV of ZnSe with multilayer coatings at 3.39 microns, scattering from contaminated surfaces, cleanliness correlation by BRDF and PFO instruments, contamination effects on optical surfaces, means of eliminating the effects of particulate contamination on scatter measurements of superfine optical surfaces, vacuum BRDF measurement of cryogenic optical surfaces, Monte Carlo simulation of contaminant transport to and deposition on complex spacecraft surfaces, surface particle observation and BRDF predictions, satellite material contaminant optical properties, dark field photographic techniques for documenting optical surface contamination, design of a laboratory study of contaminant film darkening in space, contamination monitoring approaches for EUV space optics.
The Slope Imaging Multi-polarization Photon-counting Lidar (SIMPL) is a multi-beam, micropulse airborne laser altimeter that acquires active and passive polarimetric optical remote sensing measurements at visible and near-infrared wavelengths. SIMPL was developed to demonstrate advanced measurement approaches of potential benefit for improved, more efficient spaceflight laser altimeter missions. SIMPL data have been acquired for wide diversity of forest types in the summers of 2010 and 2011 in order to assess the potential of its novel capabilities for characterization of vegetation structure and composition. On each of its four beams SIMPL provides highly-resolved measurements of forest canopy structure by detecting single-photons with 15 cm ranging precision using a narrow-beam system operating at a laser repetition rate of 11 kHz. Associated with that ranging data SIMPL provides eight amplitude parameters per beam unlike the single amplitude provided by typical laser altimeters. Those eight parameters are received energy that is parallel and perpendicular to that of the plane-polarized transmit pulse at 532 nm (green) and 1064 nm (near IR), for both the active laser backscatter retro-reflectance and the passive solar bi-directional reflectance. This poster presentation will cover the instrument architecture and highlight the performance of the SIMPL instrument with examples taken from measurements for several sites with distinct canopy structures and compositions. Specific performance areas such as probability of detection, after pulsing, and dead time, will be highlighted and addressed, along with examples of their impact on the measurements and how they limit the ability to accurately model and recover the canopy properties. To assess the sensitivity of SIMPL's measurements to canopy properties an instrument model has been implemented in the FLIGHT radiative transfer code, based on Monte Carlo simulation of photon transport. SIMPL data collected in 2010 over the Smithsonian Environmental Research Center, MD are currently being modelled and compared to other remote sensing and in situ data sets. Results on the adaptation of FLIGHT to model micropulse, single'photon ranging measurements are presented elsewhere at this conference. NASA's ICESat-2 spaceflight mission, scheduled for launch in 2016, will utilize a multi-beam, micropulse, single-photon ranging measurement approach (although non-polarimetric and only at 532 nm). Insights gained from the analysis and modelling of SIMPL data will help guide preparations for that mission, including development of calibration/validation plans and algorithms for the estimation of forest biophysical parameters.
Three dimensional unsteady and steady sediment distribution obtained by Monte Carlo simulation method, using stochastic hydrodynamics model for transport
Humans embarking on deep space exploration missions will encounter persistent exposure to galactic cosmic rays (GCR) - an energetic and highly complex radiation field that is unlike anything found on Earth. Exposure to such radiation fields is attributed to various adverse health effects, including cancer, cardiovascular disease, and cognitive impairment and is identified by NASA as one of the five main hazards of human spaceflight. It is therefore critical to be able to fully characterize the exposure received by humans behind shielding in space and project consequent health risks. A wide variety of computational models have been developed over the years to help meet this requirement. In this talk, an overview of the GCR environment in deep space is provided. Methods of propagating GCR fields through the shielding that protects humans in space are described along with simulation tools used to assess biological damage at the cellular scale. Finally, the NASA cancer risk model is briefly described, and risk projections are provided for various mission scenarios. Radiation transport models and solution methods pervade many aspects of this talk. For example, the GCR spectrum impinging on spacecraft is determined by solving the Fokker-Planck transport equation to propagate cosmic rays (believed to originate from supernovae) from the edge of the heliosphere to the vicinity of Earth. The Boltzmann transport equation is solved to transport this GCR spectrum through shielding and human tissue. Monte Carlo methods are used to simulate the transport of low energy electrons that dominate biological damage at the cellular scale. Progress and challenges in each of these areas will be highlighted.
A hybrid method has been developed for the numerical prediction of turbulent mixing in a spatially-developing, free shear layer. Most significantly, the computation incorporates the effects of large-scale structures, Schmidt number and Reynolds number on mixing, which have been overlooked in the past. In flow field prediction, large-eddy simulation was conducted by a modified 2-D vortex method with subgrid-scale modeling. The predicted mean velocities, shear layer growth rates, Reynolds stresses, and the RMS of longitudinal velocity fluctuations were found to be in good agreement with experiments, although the lateral velocity fluctuations were overpredicted. In scalar transport, the Monte Carlo method was extended to the simulation of the time-dependent pdf transport equation. For the first time, the mixing frequency in Curl's coalescence/dispersion model was estimated by using Broadwell and Breidenthal's theory of micromixing, which involves Schmidt number, Reynolds number and the local vorticity. Numerical tests were performed for a gaseous case and an aqueous case. Evidence that pure freestream fluids are entrained into the layer by large-scale motions was found in the predicted pdf. Mean concentration profiles were found to be insensitive to Schmidt number, while the unmixedness was higher for higher Schmidt number. Applications were made to mixing layers with isothermal, fast reactions. The predicted difference in product thickness of the two cases was in reasonable quantitative agreement with experimental measurements.