CAVEAT - A revised version of the general purpose Monte Carlo program, COHORT Volume 2 - Users' manual
Revised version of general purpose Monte Carlo computer code for radiation transport calculations - users manual
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.
Revised version of general purpose Monte Carlo computer code for radiation transport calculations - users manual
The mechanism of merging of like-signed aircraft vortices leading to a rapid redistribution of trailed vorticity in a wake through both convective and turbulent processes was investigated. Research was done experimentally in a small wind tunnel and analytically through the use of a code which computes turbulent transport using a second-order closure turbulent model. Computations are reported which demonstrate the merging phenomenon, and comparisons are made with experimental results. The usefulness of point vortex computations in predicting merging was explored. Limited computations showed that jet exhaust does not appreciably alter the merging phenomenon. The effect of ambient atmospheric turbulence on the aging of an aircraft wake was investigated at a constant turbulent dissipation rate. It was shown that under stable atmospheric conditions, when atmospheric macroscales are less than or equal to the vortex spacing, misleading results may be obtained.
Excess cancer risks from exposures to space radiation are estimated for various orbits of the International Space Station (ISS). Organ exposures are computed with the transport codes, BRYNTRN and HZETRN, and the computerized anatomical male and computerized anatomical female models. Cancer risk coefficients in the National Council on Radiation Protection and Measurements report No. 98 are used to generate lifetime excess cancer incidence and cancer mortality after a one-month mission to ISS. The generated data are tabulated to serve as a quick reference for assessment of radiation risk to astronauts on ISS missions.
To estimate astronaut health risk due to space radiation, one must have the ability to calculate exposure-related quantities averaged over specific organs and tissue types. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various tissues to the reference values specified by the International Commission on Radiological Protection (ICRP). Major discrepancies are found between the CAM and CAF tissue masses and the ICRP reference data for almost all of the tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN to compute mass averaged exposure quantities. A numerical algorithm is used to generate multiple point distributions for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.
A worst case assessment was performed on the Command Data System (CDS) multilayer printed circuit board and an output power transformer module in the power subsystem. An estimate of the Jovian environment during the 35 hour orbit insertion was supplied by JPL and used as an input to calculate the electron transport into the Galileo components. A radiation shielding analysis computer code, CHARGE, calculated the electron transport deposition trapped in the anticipated sensitive areas of the multilayer board and transformer module. Based on these trapped charge calculations electric fields were calculated between the identified isolated areas and the spacecraft ground. The results of the assessment of electrostatic discharge (DSD) in the CDS multilayer printed circuit board indicate that the probability of ESD in the FR4 is low. The probability of ESD in the components attached to the multilayer board, however, is uncertain based on a lack of prior experimental data.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multiple-layered materials. The code is established to operate on the Langley nuclear fragmentation model used in space engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurement. Comparison with Fe-56 ion with Pb-Al and Pb-(CH2)(x) targets shows reasonable agreement.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multiple-layered materials. The code is established to operate on the Langley nuclear fragmentation model used in space engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurement. Comparison with 56Fe ion with Pb-Al and Pb-(CH2)x targets shows reasonable agreement.
A thorough understanding of dynamic interactions between inlets and compressors is extremely important to the design and development of propulsion control systems, particularly for supersonic aircraft such as the High-Speed Civil Transport (HSCT). Computational fluid dynamics (CFD) codes are routinely used to analyze individual propulsion components. By coupling the appropriate CFD component codes, it is possible to investigate inlet-compressor interactions. The objectives of this work were to gain a better understanding of inlet-compressor interaction physics, formulate a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets, and to take a first step toward the CFD simulation of an entire engine by coupling multidimensional component codes. This work was conducted at the NASA Lewis Research Center by a team of civil servants and support service contractors as part of the High Performance Computing and Communications Program (HPCCP).
A general computer code was developed to allow calculation of atom-atom and ion-atom transport collision integrals from accurate potential energy curves described by a set of discrete data points for a broad range of scattering conditions. This code is based upon semiclassical approximations that properly account for quantum mechanical behavior such as tunneling effects near a barrier maximum, resonance charge exchange, and nuclear symmetry effects. Transport collision integrals were determined for N-N, O-O, N(+)-N, and O(+)-O interactions from complete sets of accurate potential functions derived from combined experimental and ab initio structure calculations. For the O-O case, this includes results for excited states. The calculated values of the N(+)-N and O(+)-O resonance charge exchange cross section Q(ex) agree well with measurements from beam experiment that are available at high energies where the diffusion cross section Q(d) satisfies Q(d) approximately equal to 2Q(ex).
One of the objectives of the AO 138-7 experiment on board the Long Duration Exposure Facility (LDEF) was a total dose measurement with Thermo Luminescent Detectors (TLD 100). Two identical packages, both of them including five TLD's inside various aluminum shields, are exposed to the space environment in order to obtain the absorbed dose profile. Radiation fluence received during the total mission length was computed, taking into account the trapped particles (AE8 and AP8 models during solar maximum and minimum periods) and the cosmic rays; due to the magnetospheric shielding the solar proton fluences are negligible on the LDEF orbit. The total dose induced by these radiations inside a semi infinite plane shield of aluminum are computed with the radiation transport codes available at DERTS. The dose profile obtained is in good agreement with the evaluation by E.V. Benton. TLD readings are performed after flight; due to the mission duration increase a post flight calibration was necessary in order to cover the range of the in flight induced dose. The results obtained, similar (plus or minus 30 percent) for both packages, are compared with the dose profile computation. For thick shields it seems that the measurements exceed the forecast (about 40 percent). That can be due to a cosmic ray and trapped proton contributions coming from the backside (assumed as perfectly shielded by the LDEF structure in the computation), or to an underestimate of the proton or cosmic ray fluences. A fine structural shielding analysis should be necessary in order to determine the origin of this slight discrepancy between forecast and in flight measurements. For the less shielded dosimeters, mainly exposed to the trapped electron flux, a slight overestimation of the dose (less than 40 percent) appears. Due to the dispersion of the TLD's response, this cannot be confirmed. In practice these results obtained on board LDEF, with less than a factor 1.4 between measurements and forecast, reinforce the validity of the computation methods and models used for the long term evaluation of the radiation levels (flux and dose) encountered in space on low inclination and altitude Earth orbits.
Information analysis center concept and computer codes for calculating radiation transport, and shield designs
The evolution of the protoplanet Jupiter is followed, using a hydrodynamic computer code with radiative energy transport. Jupiter is assumed to have formed as a subcondensation in the primitive solar nebula at a density just high enough for gravitational collapse to occur. The initial state has a density of 0.0015 nanograms per cu cm and a temperature of 43 K; the calculations are carried to an equilibrium state where the central density reaches 0.5 g per cu cm and the central temperature reaches 25,000 K. During the early part of the evolution the object contracts in quasi-hydrostatic equilibrium; later on hydrodynamic collapse occurs, induced by the dissociation of hydrogen molecules. After dissociation is complete, the planet regains hydrostatic equilibrium with a radius of a few times the present value. Further evolution beyond this point is not treated here; however the results are consistent with the existence of a high-luminosity phase shortly after the planet settles into its final quasi-static contraction.-
A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.
Radiation shielding analyses are performed for candidate lunar base habitation modules. The study primarily addresses potential hazards due to contributions from the galactic cosmic rays. The NASA Langley Research Center's high energy nucleon and heavy ion transport codes are used to compute propagation of radiation through conventional and regolith shield materials. Computed values of linear energy transfer are converted to biological dose-equivalent using quality factors established by the International Commision of Radiological Protection. Special fluxes of heavy charged particles and corresponding dosimetric quantities are computed for a series of thicknesses in various shield media and are used as an input data base for algorithms pertaining to specific shielded geometries. Dosimetric results are presented as isodose contour maps of shielded configuration interiors. The dose predictions indicate that shielding requirements are substantial, and an abbreviated uncertainty analysis shows that better definition of the space radiation environment as well as improvement in nuclear interaction cross-section data can greatly increase the accuracy of shield requirement predictions.
The goal of the National Aeronautics and Space Agency and the Space Radiation Project is to ensure that astronauts can safely live and work in the space radiation environment. The space radiation environment poses both acute and chronic risks to crew health and safety, but unlike some other aspects of space travel, space radiation exposure has clinically relevant implications for the lifetime of the crew. The term safely means that risks are sufficiently understood such that acceptable limits on mission, post-mission and multi-mission consequences (for example, excess lifetime fatal cancer risk) can be defined. The Space Radiation Project strategy has several elements. The first element is to use a peer-reviewed research program to increase our mechanistic knowledge and genetic capabilities to develop tools for individual risk projection, thereby reducing our dependency on epidemiological data and population-based risk assessment. The second element is to use the NASA Space Radiation Laboratory to provide a ground-based facility to study the understanding of health effects/mechanisms of damage from space radiation exposure and the development and validation of biological models of risk, as well as methods for extrapolation to human risk. The third element is a risk modeling effort that integrates the results from research efforts into models of human risk to reduce uncertainties in predicting risk of carcinogenesis, central nervous system damage, degenerative tissue disease, and acute radiation effects. To understand the biological basis for risk, we must also understand the physical aspects of the crew environment. Thus the fourth element develops computer codes to predict radiation transport properties, evaluate integrated shielding technologies and provide design optimization recommendations for the design of human space systems. Understanding the risks and determining methods to mitigate the risks are keys to a successful radiation protection strategy.
Potential curves for the higher-lying states of NO are used together with a general scattering code to compute accurate values of transport collision integrals for N(4S/0/)-O(3P) interactions. It is found that the major contribution to the collision integrals comes from these higher-lying states of NO. New values of the collision integrals for N(4S/0/)-N(4S/0/) and O(3P)-O(3P) interactions are computed as well. The tabulations of the collision integrals cover a broad range of temperatures from 250 to 100,000 K and can be used to determine transport properties such as viscosity, thermal conductivity, and the diffusion coefficient to the second order.