Influence of particle size distribution on reflected and transmitted light from clouds
Particle size distribution effects on incident radiation from clouds calculated by Monte Carlo method
Kattawar, G. W.↗
SEARCH · Search NASA
Results for “Monte Carlo N Particle”
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.
Influence of particle size distribution on reflected and transmitted light from clouds.
Particle size distribution influence on reflected and transmitted light from clouds calculated by Monte Carlo technique
Kattawar, G. W.↗
Space Radiation Transport Code Development: 3DHZETRN
The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.
Wilson, John W.↗
Cosmic ray studies with a gas Cerenkov counter in association with an ionization spectrometer
The results from a balloon-borne gas Cerenkov counter (threshold 16.5 GeV/nuc) and an ionization spectrometer are presented. The gas Cerenkov counter provides an absolute energy calibration for the response of the calorimeter for the Z range of 5-26 nuclei of cosmic rays. The contribution of scintillation to the gas Cerenkov pulse height has been obtained by independently selecting particles below the gas Cerenkov threshold using the ionization spectrometer. Energy spectra were derived by minimizing the chi-squared between a Monte Carlo simulated data and flight data. Best fit power laws were determined for C, N, O, Ne, Mg, and Si. The power laws, all consistent with E exp-2.7, are not good fits to the data. A better fit is obtained using the spectrum derived from the spectrometer.
Balasubrahmanyan, V. K.↗
Influence of particle size distribution on reflected and transmitted light from clouds Scientific interim report
Light reflection and transmission from clouds with various drop size distributions calculated using Monte Carlo technique
CLOUD PHYSICS↗
Cosmic ray studies with a gas Cherenkov counter in association with an ionization spectrometer
The results from a balloon-borne gas Cherenkov counter (threshold 16.5 GeV/nucleon) and an ionization spectrometer are presented. The gas Cherenkov counter provides an absolute energy distribution for the response of the calorimeter for 5 or = Z 26 nuclei of cosmic rays. The contribution of scintillation to the gas Cherenkov pulse height was obtained by independently selecting particles below the gas Cherenkov threshold using the ionization spectrometer. Energy spectra were derived by minimizing the chi squared between Monte Carlo simulted data and flight data. Best fit power laws, dN/dE = AE-gamma, were determined for C, N, O, Ne, Mg, and Si. The power laws, all consistent with E (-2.7) are not good fits to the data. A better fit is obtained using the spectrum derived from the spectrometer. The data from the ionization calorimeter and the gas Cherenkov are thus completely self-consistent.
Balasubrahmanyan, V. K.↗
Calculations of neutron flux spectra induced in the earth's atmosphere by galactic cosmic rays.
Calculations have been carried out to determine the neutron flux induced in the earth's atmosphere by galactic protons and alpha particles at solar minimum for a geomagnetic latitude of 42 N. Neutron flux spectra in the energy range from about 10 to the minus 8th to about 100,000 MeV at various depths in the atmosphere were calculated by using Monte Carlo and discrete ordinates methods, and various comparisons with experimental data are presented. The magnitude and shape of the calculated neutron leakage spectrum at the particular latitude considered support the theory that the cosmic ray albedo neutron decay (Crand) mechanism is the source of the protons trapped in the inner radiation belt.
Armstrong, T. W.↗
Potential, velocity, and density fields from sparse and noisy redshift-distance samples - Method
A method for recovering the three-dimensional potential, velocity, and density fields from large-scale redshift-distance samples is described. Galaxies are taken as tracers of the velocity field, not of the mass. The density field and the initial conditions are calculated using an iterative procedure that applies the no-vorticity assumption at an initial time and uses the Zel'dovich approximation to relate initial and final positions of particles on a grid. The method is tested using a cosmological N-body simulation 'observed' at the positions of real galaxies in a redshift-distance sample, taking into account their distance measurement errors. Malmquist bias and other systematic and statistical errors are extensively explored using both analytical techniques and Monte Carlo simulations.
Dekel, Avishai↗
In vitro cell irradiation systems based on 210Po alpha source: construction and characterisation
One way of studying the risk to human health of low-level radiation exposure is to make biological experiments on living cell cultures. Two 210Po alpha-particle emitting devices, with 0.5 and 100 MBq activity, were designed and constructed to perform such experiments irradiating monolayers of cells. Estimates of dose rate at the cell surface were obtained from measurements by a PIPS alpha-particle spectrometer and from calculations by the SRIM 2000, Monte Carlo charged particle transport code. Particle fluence area distributions were measured by solid state nuclear track detectors. The design and dosimetric characterisation of the devices are discussed. c2002 Elsevier Science Ltd. All rights reserved.
Alpha Particles↗
Space Applications of the FLUKA Monte-Carlo Code: Lunar and Planetary Exploration
NASA has recognized the need for making additional heavy-ion collision measurements at the U.S. Brookhaven National Laboratory in order to support further improvement of several particle physics transport-code models for space exploration applications. FLUKA has been identified as one of these codes and we will review the nature and status of this investigation as it relates to high-energy heavy-ion physics.
Anderson, V.↗
Monte Carlo calculation of dynamical properties of the two-dimensional Hubbard model
A new method is introduced for analytically continuing imaginary-time data from quantum Monte Carlo calculations to the real-frequency axis. The method is based on a least-squares-fitting procedure with constraints of positivity and smoothness on the real-frequency quantities. Results are shown for the single-particle spectral-weight function and density of states for the half-filled, two-dimensional Hubbard model.
White, S. R.↗
Proton Straggling in Thick Silicon Detectors
Straggling functions for protons in thick silicon radiation detectors are computed by Monte Carlo simulation. Mean energy loss is constrained by the silicon stopping power, providing higher straggling at low energy and probabilities for stopping within the detector volume. By matching the first four moments of simulated energy-loss distributions, straggling functions are approximated by a log-normal distribution that is accurate for Vavilov k is greater than or equal to 0:3. They are verified by comparison to experimental proton data from a charged particle telescope.
Selesnick, R. S.↗
Study of Analytic Statistical Model for Decay of Light and Medium Mass Nuclei in Nuclear Fragmentation
The angular momentum independent statistical decay model is often applied using a Monte-Carlo simulation to describe the decay of prefragment nuclei in heavy ion reactions. This paper presents an analytical approach to the decay problem of nuclei with mass number less than 60, which is important for galactic cosmic ray (GCR) studies. This decay problem of nuclei with mass number less than 60 incorporates well-known levels of the lightest nuclei (A less than 11) to improve convergence and accuracy. A sensitivity study of the model level density function is used to determine the impact on mass and charge distributions in nuclear fragmentation. This angular momentum independent statistical decay model also describes the momentum and energy distribution of emitted particles (n, p, d, t, h, and a) from a prefragment nucleus.
Cucinotta, Francis A.↗
Multiple scattering in reflection nebulae. I - A Monte Carlo approach. II - Uniform plane-parallel nebulae with foreground stars. III - Nebulae with embedded illuminating stars
A method for calculating the surface brightness distribution on a plane-parallel reflection nebula of uniform density illuminated by a star located either in front of, behind, or arbitrarily inside the scattering medium is proposed. The Monte Carlo technique is used to find solutions to the radiative transfer problem. The scattering properties of the nebular particles are parameterized by the albedo for single scattering and a three-parameter analytic phase function. Calculations are then presented for the surface brightness distribution across the face of such nebulae with (1) a foreground star and (2) immersed stars. The calculations include the full effects of multiple scattering, are independent of a particular assumed grain material or size distribution, and are applicable to any wavelength region for which observations can be obtained.
Witt, A. N.↗
Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators
Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.
Ablators↗
Coherent-phase or random-phase acceleration of electron beams in solar flares
Time structures of electron beam signatures at radio wavelengths are investigated to probe correlated versus random behavior in solar flares. In particular we address the issue whether acceleration and injection of electron beams is coherently modulated by a single source, or whether the injection is driven by a stochastic (possibly spatially fragmented) process. We analyze a total of approximately = 6000 type III bursts observed by Ikarus (Zurich) in the frequency range of 100-500 MHz, during 359 solar flares with simultaneous greater than or = 25 keV hard X-ray emission, in the years 1890-1983. In 155 flares we find a total of 260 continuous type III groups, with an average number of 13 +/- 9 bursts per group, a mean duration of D = 12 +/- 14 s, a mean period of P = 2.0 +/- 1.2 s, with the highest burst rate at a frequency of nu = 310 +/- 120 MHz. Pulse periods have been measured between 0.5 and 10 s, and can be described by an exponential distribution, i.e., N(P) varies as e (exp -P/1.0s). The period shows a frequency dependence of P(nu)=46(exp-0.6)(sub MHz)s for different flares, but is invariant during a particular flare. We measure the mean period P and its standard deviation sigma (sub p) in each type III group, and quantify the degree of periodicity (or phase-coherence) by the dimensionless parameter sigma (sub p)P. The representative sample of 260 type III burst groups shows a mean periodicity of sigma (sub p/P) = 0.37 +/- 0.12, while Monte Carlo simulations of an equivalent set of truly random time series show a distinctly different value of sigma (sub p)P = 0.93 +/- 0.26. This result indicates that the injection of electron beams is coherently modulated by a particle acceleration source which is either compact or has a global organization on a timescale of seconds, in contrast to an incoherent acceleration source, which is stochastic either in time or space. We discuss the constraints on the size of the acceleration region resulting from electron beam propagation delays and from Alfvenic synchronization during a pulse period. We discuss two periodic preocesses in flares, which potentially control quasi-periodic particle acceleration: (1) MHD oscillations, and (2) current sheets with oscillatory dynamics.
Aschwanden, Markus J.↗
Issues in the reconstruction of environmental doses on the basis of thermoluminescence measurements in the Techa riverside
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.
NASA Discipline Radiation Health↗
Software for Processing Flight and Simulated Data of the ATIC Experiment
ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. A C++-class library for building different programs for processing flight and simulated data of the ATIC balloon experiment is described. This library is compatible with the ROOT-system and includes classes and methods for solving a number of problems as the following: Reading data files in different formats (raw-data format, ROOT-format, ASCII-format, different formats for simulated data); Transferring all these formats to the only inner format of the library; Reconstruction of trajectories of primary particles with BGO calorimeter only. The Monte-Carlo simulations with GEANT code were used to obtain the basic tables for computing error corridors and chi(sup 2)-values for the trajectories. Obtaining error corridors for searching for signal of primary particle in the Si-matrix; Searching for hit of primary particle in the Si-matrix with using of error corridor and other criteria (chi(sup 2)-values, agreement between signals in Si-matrix and in the upper layer of scintillator and others); Determination of charge of primary particle; Determination of energy deposit in BGO calorimeter.
Panov, A. D.↗