Los Alamos radiation transport code system on desktop computing platforms
Briesmeister, J.F.; Brinkley, F.W.; Clark, B.A.; West, J.T. )
1990-01-01
The Los Alamos Radiation Transport Code System (LARTCS) consists of state-of-the-art Monte Carlo and discrete ordinates transport codes and data libraries. These codes were originally developed many years ago and have undergone continual improvement. With a large initial effort and continued vigilance, the codes are easily portable from one type of hardware to another. The performance of scientific work-stations (SWS) has evolved to the point that such platforms can be used routinely to perform sophisticated radiation transport calculations. As the personal computer (PC) performance approaches that of the SWS, the hardware options for desk-top radiation transport calculations expands considerably. The current status of the radiation transport codes within the LARTCS is described: MCNP, SABRINA, LAHET, ONEDANT, TWODANT, TWOHEX, and ONELD. Specifically, the authors discuss hardware systems on which the codes run and present code performance comparisons for various machines.
Recent developments in the Los Alamos radiation transport code system
Forster, R.A.; Parsons, K.
1997-06-01
A brief progress report on updates to the Los Alamos Radiation Transport Code System (LARTCS) for solving criticality and fixed-source problems is provided. LARTCS integrates the Diffusion Accelerated Neutral Transport (DANT) discrete ordinates codes with the Monte Carlo N-Particle (MCNP) code. The LARCTS code is being developed with a graphical user interface for problem setup and analysis. Progress in the DANT system for criticality applications include a two-dimensional module which can be linked to a mesh-generation code and a faster iteration scheme. Updates to MCNP Version 4A allow statistical checks of calculated Monte Carlo results.
Urbatsch, Todd James
2015-06-15
We present an overview of radiation transport, covering terminology, blackbody raditation, opacities, Boltzmann transport theory, approximations to the transport equation. Next we introduce several transport methods. We present a section on Caseology, observing transport boundary layers. We briefly broach topics of software development, including verification and validation, and we close with a section on high energy-density experiments that highlight and support radiation transport.
DETERMINISTIC TRANSPORT METHODS AND CODES AT LOS ALAMOS
J. E. MOREL
1999-06-01
The purposes of this paper are to: Present a brief history of deterministic transport methods development at Los Alamos National Laboratory from the 1950's to the present; Discuss the current status and capabilities of deterministic transport codes at Los Alamos; and Discuss future transport needs and possible future research directions. Our discussion of methods research necessarily includes only a small fraction of the total research actually done. The works that have been included represent a very subjective choice on the part of the author that was strongly influenced by his personal knowledge and experience. The remainder of this paper is organized in four sections: the first relates to deterministic methods research performed at Los Alamos, the second relates to production codes developed at Los Alamos, the third relates to the current status of transport codes at Los Alamos, and the fourth relates to future research directions at Los Alamos.
Radiation transport. Progress report, October 1, 1982-March 31, 1983
O'Dell, R.D.
1984-05-01
Research and development progress in radiation transport by the Los Alamos National Laboratory's Group X-6 for the first half of FY 83 is reported. Included are tasks in the areas of Fission Reactor Neutronics, Deterministic Transport Methods, and Monte Carlo Radiation Transport.
Radiation transport. Progress report, April 1-December 31, 1983
O'Dell, R.D.
1984-10-01
Research and development progress in radiation transport by the Los Alamos National Laboratory's Group X-6 for the last nine months of CY 83 is reported. Included are unclassified tasks in the areas of Fission Reactor Neutronics, Deterministic Transport Methods, Monte Carlo Radiation Transport, and Cross Sections and Physics.
NASA Technical Reports Server (NTRS)
Friedel, R. H. W.; Bourdarie, S.; Fennell, J.; Kanekal, S.; Cayton, T. E.
2004-01-01
The highly energetic electron environment in the inner magnetosphere (GEO inward) has received a lot of research attention in resent years, as the dynamics of relativistic electron acceleration and transport are not yet fully understood. These electrons can cause deep dielectric charging in any space hardware in the MEO to GEO region. We use a new and novel approach to obtain a global representation of the inner magnetospheric energetic electron environment, which can reproduce the absolute environment (flux) for any spacecraft orbit in that region to within a factor of 2 for the energy range of 100 KeV to 5 MeV electrons, for any levels of magnetospheric activity. We combine the extensive set of inner magnetospheric energetic electron observations available at Los Alamos with the physics based Salammbo transport code, using the data assimilation technique of "nudging". This in effect input in-situ data into the code and allows the diffusion mechanisms in the code to interpolate the data into regions and times of no data availability. We present here details of the methods used, both in the data assimilation process and in the necessary inter-calibration of the input data used. We will present sample runs of the model/data code and compare the results to test spacecraft data not used in the data assimilation process.
Introduction to radiation transport
Olson, G.L.
1998-12-31
This lecture will present time-dependent radiation transport where the radiation is coupled to a static medium, i.e., the material is not in motion. In reality, radiation exerts a pressure on the materials it propagates through and will accelerate the material in the direction of the radiation flow. This fully coupled problem with radiation transport and materials in motion is referred to as radiation-hydrodynamics (or in a shorthand notation: rad-hydro) and is beyond the scope of this lecture.
Transportation of pyrochemical salts from Rocky Flats to Los Alamos
Schreiber, S.B.
1997-02-01
Radioactive legacy wastes or residues are currently being stored on numerous Sites around the former Department of Energy`s (DOE) Nuclear Weapons Complex. Since most of the operating facilities were shut down and have not operated since before the declared end to the Cold War in 1993, the historical method for treating these residues no longer exists. The risk associated with continued storage of these residues will dramatically increase with time. Thus, the DOE was directed by the Defense Nuclear Facility Safety Board in its Recommendation 94-1 to address and stabilize these residues and established an eight year time frame for doing so. There are only two options available to respond to this requirement: (1) restart existing facilities to treat and package the residues for disposal or (2) transport the residues to another operating facility within the Complex where they can be treated and packaged for disposal. This paper focuses on one such residue type, pyrochemical salts, produced at one Complex site, the Rocky Flats Plant located northwest of Denver, Colorado. One option for treating the salts is their shipment to Los Alamos, New Mexico, for handling at the Plutonium Facility. The safe transportation of these salts can be accomplished at present with several shipping containers including a DOT 6M, a DOE 9968, Type A or Type B quantity 55-gallon drum overpacks, or even the TRUPACT II. The tradeoffs between each container is examined with the conclusion that none of the available shipping containers is fully satisfactory. Thus, the advantageous aspects of each container must be utilized in an integrated and efficient way to effectively manage the risk involved. 1 fig.
Penetrating radiation: applications at Los Alamos National Laboratory
NASA Astrophysics Data System (ADS)
Watson, Scott; Hunter, James; Morris, Christopher
2013-09-01
Los Alamos has used penetrating radiography extensively throughout its history dating back to the Manhattan Project where imaging dense, imploding objects was the subject of intense interest. This interest continues today as major facilities like DARHT1 have become the mainstay of the US Stockpile Stewardship Program2 and the cornerstone of nuclear weapons certification. Meanwhile, emerging threats to national security from cargo containers and improvised explosive devices (IEDs) have invigorated inspection efforts using muon tomography, and compact x-ray radiography. Additionally, unusual environmental threats, like those from underwater oil spills and nuclear power plant accidents, have caused renewed interest in fielding radiography in severe operating conditions. We review the history of penetrating radiography at Los Alamos and survey technologies as presently applied to these important problems.
Radiation transport in inhomogeneous media
Keiter, Paul; Gunderson, Mark; Foster, John; Rosen, Paula; Comley, Andrew; Taylor, Mark; Perry, Ted
2008-05-15
Calculations of radiation transport in heated materials are greatly complicated by the presence of regions in which two or more materials are inhomogeneously mixed. This phenomenon is important in many systems, such as astrophysical systems where density clumps can be found in star-forming regions and molecular clouds. Laboratory experiments have been designed to test the modeling of radiation transport through inhomogeneous plasmas. A laser-heated hohlraum is used as a thermal source to drive radiation through polymer foam containing randomly distributed gold particles. Experimental measurements of radiation transport in foams with gold particle sizes ranging from 5-9 {mu}m to submicrometer diameters as well as the homogeneous foam case are presented. The simulation results of the radiation transport are compared to the experiment and show that an inhomogeneous transport model must be applied to explain radiation transport in foams loaded with 5 {mu}m diameter gold particles.
THE LOS ALAMOS NATIONAL LABORATORY ATMOSPHERIC TRANSPORT AND DIFFUSION MODELS
M. WILLIAMS
1999-08-01
The LANL atmospheric transport and diffusion models are composed of two state-of-the-art computer codes. The first is an atmospheric wind model called HOThlAC, Higher Order Turbulence Model for Atmospheric circulations. HOTMAC generates wind and turbulence fields by solving a set of atmospheric dynamic equations. The second is an atmospheric diffusion model called RAPTAD, Random Particle Transport And Diffusion. RAPTAD uses the wind and turbulence output from HOTMAC to compute particle trajectories and concentration at any location downwind from a source. Both of these models, originally developed as research codes on supercomputers, have been modified to run on microcomputers. Because the capability of microcomputers is advancing so rapidly, the expectation is that they will eventually become as good as today's supercomputers. Now both models are run on desktop or deskside computers, such as an IBM PC/AT with an Opus Pm 350-32 bit coprocessor board and a SUN workstation. Codes have also been modified so that high level graphics, NCAR Graphics, of the output from both models are displayed on the desktop computer monitors and plotted on a laser printer. Two programs, HOTPLT and RAPLOT, produce wind vector plots of the output from HOTMAC and particle trajectory plots of the output from RAPTAD, respectively. A third CONPLT provides concentration contour plots. Section II describes step-by-step operational procedures, specifically for a SUN-4 desk side computer, on how to run main programs HOTMAC and RAPTAD, and graphics programs to display the results. Governing equations, boundary conditions and initial values of HOTMAC and RAPTAD are discussed in Section III. Finite-difference representations of the governing equations, numerical solution procedures, and a grid system are given in Section IV.
Automatic beam position control at Los Alamos Spallation Radiation Effects Facility (LASREF)
Oothoudt, M.; Pillai, C.; Zumbro, M.
1997-08-01
Historically the Los Alamos Spallation Radiation Effects Facility (LASREF) has used manual methods to control the position of the 800 kW, 800 MeV proton beam on targets. New experiments, however, require more stringent position control more frequently than can be done manually for long periods of time. Data from an existing harp is used to automatically adjust steering magnets to maintain beam position to required tolerances.
Preliminary results from the Los Alamos TA54 complex terrain Atmospheric Transport Study (ATS)
Vold, E.; Chan, M.; Sanders, L.
1995-09-01
The Los Alamos National Laboratory (LANL) Low-Level Radioactive Waste (LLRW) disposal site at TA54, Area G la located on a mesa top amidst a complex terrain of finger like mesas typically 30 motors or more In height above canyons of widths varying from 100 to 300 motors. Atmospheric dispersion from this site is of concern for routine operations and for potential Incidents during waste retrieval operations. Indian lands are located In the dominant downwind direction within 500 m from the site and provide further incentive to understand the potential and actual impacts of waste disposal operations. The permanent network of meteorological towers at LANL have been located primarily at mesa-top locations to coincide with most laboratory facilities and as such do not resolve the effects of channeling in the canyons and the influence this has on potential surface releases. An Atmospheric Transport Study (ATS) was initiated to better understand the wind flow fields and dispersion from the LANL Waste Storage and Disposal facilities at TA-54, Area G. As part of this effort, a series of six portable meteorological towers were sited in the vicinity of Area G, two at mesa top locations, one just east of the site where the mesas have dissipated to mild ridges, and three in the canyons adjacent to the disposal site mesa as indicated on the topographic representation of the local terrain. Since 1994, the towers have collected horizontal wind velocities, pressure, temperature, relative humidity and a radiation gamma reading every fifteen minutes. The data bass is being analyzed for trends and to provide a basis for comparison to computational modeling efforts to predict the flow fields.
Deterministic methods in radiation transport
Rice, A.F.; Roussin, R.W.
1992-06-01
The Seminar on Deterministic Methods in Radiation Transport was held February 4--5, 1992, in Oak Ridge, Tennessee. Eleven presentations were made and the full papers are published in this report, along with three that were submitted but not given orally. These papers represent a good overview of the state of the art in the deterministic solution of radiation transport problems for a variety of applications of current interest to the Radiation Shielding Information Center user community.
Radiation transport calculations for cosmic radiation.
Endo, A; Sato, T
2012-01-01
The radiation environment inside and near spacecraft consists of various components of primary radiation in space and secondary radiation produced by the interaction of the primary radiation with the walls and equipment of the spacecraft. Radiation fields inside astronauts are different from those outside them, because of the body's self-shielding as well as the nuclear fragmentation reactions occurring in the human body. Several computer codes have been developed to simulate the physical processes of the coupled transport of protons, high-charge and high-energy nuclei, and the secondary radiation produced in atomic and nuclear collision processes in matter. These computer codes have been used in various space radiation protection applications: shielding design for spacecraft and planetary habitats, simulation of instrument and detector responses, analysis of absorbed doses and quality factors in organs and tissues, and study of biological effects. This paper focuses on the methods and computer codes used for radiation transport calculations on cosmic radiation, and their application to the analysis of radiation fields inside spacecraft, evaluation of organ doses in the human body, and calculation of dose conversion coefficients using the reference phantoms defined in ICRP Publication 110. PMID:23089013
Becker, N.M.; Van Eeckhout, E.; David, N.A.; Irvine, J.M.
1995-10-01
During World War II, Los Alamos, New Mexico was chosen as the site for the secret development of the first atomic bomb. The remote location in the southwestern United States was ideal for such a project. After the war, research activities continued at the Los Alamos installation, focusing on new nuclear weapons models as well as greater effectiveness and reliability of existing weapons. Due to the emphasis on nuclear and non-nuclear weapons development as well as associated nuclear research, a large inventory of radionuclides and heavy metals have been tested, expended, and disposed of in the local environment, a high plateau of tuffaceous volcanic rocks incised by deep canyons in a semi-arid climate. In recent years an intensive evaluation of the environmental, impact of weapons testing at Los Alamos and elsewhere has been undertaken. GIS system utilization and image processing of past and current data has been an important part of this evaluation. Important problems can be more easily displayed and understood using this methodology. The main objective in this paper is to illustrate how transport of depleted uranium and associated heavy metals (copper in this case) used in dynamic testing of weapons components at open air firing sites can be evaluated and visualized. In our studies, surface water has been found to be the predominant transport mechanism. We have sampled soils, sediments, fallout, runoff water and snowmelt over a number of years in order to understand contaminant transport on- and offsite. Statistical analyses of these data have assisted in our characterization of issues such as contaminant variability, spatially and temporally, as well as in development of transport rates.
NASA Astrophysics Data System (ADS)
Longmire, P.; Ding, M.; Rearick, M.; Vaniman, D.; Katzman, D.
2008-12-01
The environmental aqueous geochemistry of Cr is of considerable interest to physical scientists and toxicologists in quantifying the fate and transport of this metal in surface and subsurface environments. Chromium(VI) solutions were released from cooling towers to a stream channel within Sandia Canyon at Los Alamos National Laboratory, NM from 1956 to 1971. These solutions have migrated 293 m depth through the vadose zone, containing several saturated zones, to the regional water table. Concentrations of total dissolved Cr, mainly as Cr(VI), in the regional aquifer range between 0.17 to 8.46 mM. The regional aquifer is characterized by calcium-sodium-bicarbonate solution, contains dissolved oxygen (0.09 to 0.22 mM), and has a circumneutral pH (6.8 to 8.3). Geochemical processes controlling the fate and transport of Cr in groundwater at Los Alamos include a combination of adsorption and precipitation reactions within aquifer systems. Vadose zone material containing hydrous ferric oxide, smectite, silica glass, and calcite widely range in their ability to adsorb Cr(VI) under basic pH conditions. Overall, the vadose zone at Los Alamos is relatively oxidizing, however, basalt flows are locally reducing with respect to Fe. Ferrous iron concentrated within the Cerros del Rio basalt has been shown through batch experiments to reduce Cr(VI) to Cr(III) resulting in precipitation of chromium(III) hydroxide. Regional aquifer material, consisting of silicates, oxides, and calcite, vary in the amount of Fe(II) available in reactive minerals to effectively reduce Cr(VI) to Cr(III). The results of our studies (1) directly assess the relationship between mineralogical characterization and transport behavior of Cr using site-specific hydrogeologic material and (2) provide site-specific adsorption and precipitation parameters obtained through the experiments to refine the fate and transport modeling of Cr within the vadose zone and regional aquifer. Natural attenuation of Cr at Los
Radiation transport in numerical astrophysics
Lund, C.M.
1983-02-01
In this article, we discuss some of the numerical techniques developed by Jim Wilson and co-workers for the calculation of time-dependent radiation flow. Difference equations for multifrequency transport are given for both a discrete-angle representation of radiation transport and a Fick's law-like representation. These methods have the important property that they correctly describe both the streaming and diffusion limits of transport theory in problems where the mean free path divided by characteristic distances varies from much less than one to much greater than one. They are also stable for timesteps comparable to the changes in physical variables, rather than being limited by stability requirements.
Graf, W.L.
1996-10-01
Between 1945 and 1952 the development of nuclear weapons at Los Alamos National Laboratory, New Mexico, resulted in the disposal of plutonium into the alluvium of nearby Acid and (to a lesser degree) DP Canyons. The purpose of this paper is to explore the connection between the disposal sites and the main river, a 20 km link formed by the fluvial system of Acid, Pueblo, DP, and Los Alamos Canyons. Empirical data from 15 yr of annual sediment sampling throughout the canyon system has produced 458 observations of plutonium concentration in fluvial sediments. These data show that, overall, mean plutonium concentrations in fluvial sediment decline from 10,000 fCi/g near the disposal area to 100 fCi/g at the confluence of the canyon system and the Rio Grande. Simulations using a computer model for water, sediment, and plutonium routing in the canyon system show that discharges as large as the 25 yr event would fail to develop enough transport capacity to completely remove the contaminated sediments from Pueblo Canyon. Lesser flows would move some materials to the Rio Grande by remobilization of stored sediments. The simulations also show that the deposits and their contaminants have a predictable geography because they occur where stream power is low, hydraulic resistance is high, and the geologic and/or geomorphic conditions provide enough space for storage. 38 refs., 13 figs., 1 tab.
Reeves, Geoffrey D; Friedel, Reiner H W; Chen, Yue; Koller, Josef; Henderson, Michael G
2008-01-01
The Dynamic Radiation Environment Assimilation Model (DREAM) was developed at Los Alamos National Laboratory to assess, quantify, and predict the hazards from the natural space environment and the anthropogenic environment produced by high altitude nuclear explosions (HANE). DREAM was initially developed as a basic research activity to understand and predict the dynamics of the Earth's Van Allen radiation belts. It uses Kalman filter techniques to assimilate data from space environment instruments with a physics-based model of the radiation belts. DREAM can assimilate data from a variety of types of instruments and data with various levels of resolution and fidelity by assigning appropriate uncertainties to the observations. Data from any spacecraft orbit can be assimilated but DREAM was designed to function with as few as two spacecraft inputs: one from geosynchronous orbit and one from GPS orbit. With those inputs, DREAM can be used to predict the environment at any satellite in any orbit whether space environment data are available in those orbits or not. Even with very limited data input and relatively simple physics models, DREAM specifies the space environment in the radiation belts to a high level of accuracy. DREAM has been extensively tested and evaluated as we transition from research to operations. We report here on one set of test results in which we predict the environment in a highly-elliptical polar orbit. We also discuss long-duration reanalysis for spacecraft design, using DREAM for real-time operations, and prospects for 1-week forecasts of the radiation belt environment.
FY2008 Report on GADRAS Radiation Transport Methods.
Mattingly, John K.; Mitchell, Dean James; Harding, Lee; Varley, Eric S.; Hilton, Nathan R.
2008-10-01
The primary function of the Gamma Detector Response and Analysis Software (GADRAS) is the solution of inverse radiation transport problems, by which the con-figuration of an unknown radiation source is inferred from one or more measured radia-tion signatures. GADRAS was originally developed for the analysis of gamma spec-trometry measurements. During fiscal years 2007 and 2008, GADRAS was augmented to implement the simultaneous analysis of neutron multiplicity measurements. This report describes the radiation transport methods developed to implement this new capability. This work was performed at the direction of the National Nuclear Security Administration's Office of Nonproliferation Research and Development. It was executed as an element of the Proliferation Detection Program's Simulation, Algorithm, and Modeling element. Acronyms BNL Brookhaven National Laboratory CSD Continuous Slowing-Down DU depleted uranium ENSDF Evaluated Nuclear Structure Data Files GADRAS Gamma Detector Response and Analysis Software HEU highly enriched uranium LANL Los Alamos National Laboratory LLNL Lawrence Livermore National Laboratory NA-22 Office of Nonproliferation Research and Development NNDC National Nuclear Data Center NNSA National Nuclear Security Administration ODE ordinary differential equation ONEDANT One-dimensional diffusion accelerated neutral particle transport ORNL Oak Ridge National Laboratory PARTISN Parallel time-dependent SN PDP Proliferation Detection Program RADSAT Radiation Scenario Analysis Toolkit RSICC Radiation Safety Information Computational Center SAM Simulation, Algorithms, and Modeling SNL Sandia National Laboratories SNM special nuclear material ToRI Table of Radioactive Isotopes URI uniform resource identifier XML Extensible Markup Language
Biomedical program leading to therapeutic trials on pion radiation at Los Alamos.
Kligerman, M M; Knapp, E A; Petersen, D F
1975-11-01
Hypoxia and variations in cell cycle phase protect tumor cells being treated with x rays or gamma rays (cobalt). Heavy particles can overcome these protective effects, because of the dense ionization they deposit in tissues. Pions (negative pi mesons) can be directed to and stopped in a specific area, where they are captured by the nuclei of atoms, rendering the nuclei unstable. The nuclei disintegrate, releasing densely ionizing radiation. By confining the dense ionization to the tumor-bearing volume, pions have the potential of increasing the tolerance of the area under treatment to radiation, thus increasing the probability of destroying the tumor. A special channel at the proton factory at the Los Alamos Scientific Laboratory is producing pions for biomedical research. Considerable physical dosimetry has been completed. Cellular studies are underway to provide depth-dose-biological-effect curves. Animal studies will provide information on acute and late effects, which will permit the safe application of pions to a series of anatomical sites established by protocols for radiotherapy clinical trials. PMID:1192358
Becker, N.M.; Vanta, E.B.
1995-05-01
Hydrologic investigations on depleted uranium fate and transport associated with dynamic testing activities were instituted in the 1980`s at Los Alamos National Laboratory and Eglin Air Force Base. At Los Alamos, extensive field watershed investigations of soil, sediment, and especially runoff water were conducted. Eglin conducted field investigations and runoff studies similar to those at Los Alamos at former and active test ranges. Laboratory experiments complemented the field investigations at both installations. Mass balance calculations were performed to quantify the mass of expended uranium which had transported away from firing sites. At Los Alamos, it is estimated that more than 90 percent of the uranium still remains in close proximity to firing sites, which has been corroborated by independent calculations. At Eglin, we estimate that 90 to 95 percent of the uranium remains at test ranges. These data demonstrate that uranium moves slowly via surface water, in both semi-arid (Los Alamos) and humid (Eglin) environments.
Los Alamos Science: Number 23, 1995. Radiation protection and the human radiation experiments
Cooper, N.G.
1995-12-31
There are a variety of myths and misconceptions about the ionizing radiation that surrounds and penetrates us all. Dispel a few of these by taking a leisurely tour of radiation and its properties, of the natural and man-made sources of ionizing radiation, and of the way doses are calculated. By damaging DNA and inducing genetic mutations, ionizing radiation can potentially initiate a cell on the road to cancer. The authors review what is currently known about regulation of cellular reproduction, DNA damage and repair, cellular defense mechanisms, and the specific cancer-causing genes that are susceptible to ionizing radiation. A rapid survey of the data on radiation effects in humans shows that high radiation doses increase the risk of cancer, whereas the effects of low doses are very difficult to detect. The hypothetical risks at low doses, which are estimated from the atomic-bomb survivors, are compared to the low-dose data so that the reader can assess the present level of uncertainty. As part of the openness initiative, ten individuals who have worked with plutonium during various periods in the Laboratory`s history were asked to share their experiences including their accidental intakes. The history and prognosis of people who have had plutonium exposures is discussed by the Laboratory`s leading epidemiologist.
Boenig, H.J.; Rogers, J.D.; McLelland, G.R.; Pelts, C.T.; McLelland Engineering, Dallas, TX; Reliance Crane and Rigging, Inc., Phoenix, AZ )
1989-01-01
A 1430 MVA steam turbine generator was acquired from a cancelled nuclear power plant in Tennessee to be used as the pulsed power and energy storage unit for the Confinement Physics Research Facility being built at Los Alamos, NM. The transportation from Hartsville, near Nashville, TN, to Los Alamos, NM, of the two largest single pieces of the generator, a 451 t stator and a 234 t rotor presented a special challenge. Details of the move, by barge from Hartsville to Catoosa, near Tulsa, OK, by rail from Catoosa to Lamy, near Santa Fe, NM, and by road from Lamy to Los Alamos are described. The greatest difficulty of the successful move was the crossing of the Rio Grande river on an existing reinforced concrete bridge. The two-lane wide road transporters for the stator and rotor were fitted with outriggers to provide a four-lane wide vehicle, thus spreading the load over the entire bridge width and meeting acceptable load distribution and bridge safety factors. 2 refs., 6 figs.
Space Radiation Transport Methods Development
NASA Astrophysics Data System (ADS)
Wilson, J.; Tripathi, R.; Qualls, G.; Cucinotta, F.; Prael, R.; Norbury, J.
Early space radiation shield code development relied on Monte Carlo methods for proton, neutron and pion transport and made important contributions to the space program. More recently Monte Carlo code LAHET has been upgraded to include high-energy multiple-charged light ions for GCR simulations and continues to be expanded in capability. To compensate for low computational efficiency, Monte Carlo methods have resorted to restricted one-dimensional problems leading to imperfect representations of appropriate boundary conditions. Even so, intensive computational requirements resulted and shield evaluation was made near the end of the design process and resolving shielding issues usually had a negative impact on the design. We evaluate the implications of these common one-dimensional assumptions on the evaluation of the Shuttle internal radiation field. Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be
The Radiation Transport Conundrum in Radiation Hydrodynamics
Castor, J I
2005-03-18
The summary of this paper is: (1) The conundrum in the title is whether to treat radiation in the lab frame or the comoving frame in a radiation-hydrodynamic problem; (2) Several of the difficulties are associated with combining a somewhat relativistic treatment of radiation with a non-relativistic treatment of hydrodynamics; (3) The principal problem is a tradeoff between easily obtaining the correct diffusion limit and describing free-streaming radiation with the correct wave speed; (4) The computational problems of the comoving-frame formulation in more than one dimension, and the difficulty of obtaining both exact conservation and full u/c accuracy argue against this method; (5) As the interest in multi-D increases, as well as the power of computers, the lab-frame method is becoming more attractive; and (6) The Monte Carlo method combines the advantages of both lab-frame and comoving-frame approaches, its only disadvantage being cost.
Coupled electron-photon radiation transport
Lorence, L.; Kensek, R.P.; Valdez, G.D.; Drumm, C.R.; Fan, W.C.; Powell, J.L.
2000-01-17
Massively-parallel computers allow detailed 3D radiation transport simulations to be performed to analyze the response of complex systems to radiation. This has been recently been demonstrated with the coupled electron-photon Monte Carlo code, ITS. To enable such calculations, the combinatorial geometry capability of ITS was improved. For greater geometrical flexibility, a version of ITS is under development that can track particles in CAD geometries. Deterministic radiation transport codes that utilize an unstructured spatial mesh are also being devised. For electron transport, the authors are investigating second-order forms of the transport equations which, when discretized, yield symmetric positive definite matrices. A novel parallelization strategy, simultaneously solving for spatial and angular unknowns, has been applied to the even- and odd-parity forms of the transport equation on a 2D unstructured spatial mesh. Another second-order form, the self-adjoint angular flux transport equation, also shows promise for electron transport.
Status of Monte Carlo at Los Alamos
Thompson, W.L.; Cashwell, E.D.
1980-01-01
At Los Alamos the early work of Fermi, von Neumann, and Ulam has been developed and supplemented by many followers, notably Cashwell and Everett, and the main product today is the continuous-energy, general-purpose, generalized-geometry, time-dependent, coupled neutron-photon transport code called MCNP. The Los Alamos Monte Carlo research and development effort is concentrated in Group X-6. MCNP treats an arbitrary three-dimensional configuration of arbitrary materials in geometric cells bounded by first- and second-degree surfaces and some fourth-degree surfaces (elliptical tori). Monte Carlo has evolved into perhaps the main method for radiation transport calculations at Los Alamos. MCNP is used in every technical division at the Laboratory by over 130 users about 600 times a month accounting for nearly 200 hours of CDC-7600 time.
Comparison of beam transport simulations to measurements at the Los Alamos Proton Storage Ring
Wilkinson, C.; Neri, F.; Fitzgerald, D.H.; Blind, B.; Macek, R.; Plum, M.; Sander, O.; Thiessen, H.A.
1997-10-01
The ability to model and simulate beam behavior in the Proton Storage Ring (PSR) of the Los Alamos Neutron Science Center (LANSCE) is an important diagnostic and predictive tool. This paper gives the results of an effort to model the ring apertures and lattice and use beam simulation programs to track the beam. The results are then compared to measured activation levels from beam loss in the ring. The success of the method determines its usefulness in evaluating the effects of planned upgrades to the Proton Storage Ring.
Radiation transport in dust in disk geometry
NASA Technical Reports Server (NTRS)
Chun, Ming Leung
1986-01-01
The main objective of the research program is twofold: (1) to develop a computer code to solve the problem of scattering, absorption and emission of photons by dust grains in a dusty medium with 2 dimensional disk geometry, and (2) to study the various physical and geometrical effects of 2 dimensional radiation transport on the thermal structure and radiation field. These tasks were accomplished and are briefly summarized. The method for solving the radiation transport problem in disk geometry is a generalization of the quasi-diffusion method (QDM) previously developed by the author.
Radiation Transport Calculations and Simulations
Fasso, Alberto; Ferrari, A.; /CERN
2011-06-30
This article is an introduction to the Monte Carlo method as used in particle transport. After a description at an elementary level of the mathematical basis of the method, the Boltzmann equation and its physical meaning are presented, followed by Monte Carlo integration and random sampling, and by a general description of the main aspects and components of a typical Monte Carlo particle transport code. In particular, the most common biasing techniques are described, as well as the concepts of estimator and detector. After a discussion of the different types of errors, the issue of Quality Assurance is briefly considered.
Radiative Transport Based Frequency Domain Fluorescence Tomography
Joshi, Amit; Rasmussen, John C.; Sevick-Muraca, Eva M.; Wareing, Todd A.; McGhee, John
2011-01-01
We report the development of radiative transport model based fluorescence optical tomography from frequency domain boundary measurements. The coupled radiative transport model for describing NIR fluorescence propagation in tissue is solved by a novel software based on the established Attila™ particle transport simulation platform. The proposed scheme enables the prediction of fluorescence measurements with non-contact sources and detectors at minimal computational cost. An adjoint transport solution based fluorescence tomography algorithm is implemented on dual grids to efficiently assemble the measurement sensitivity Jacobian matrix. Finally, we demonstrate fluorescence tomography on a realistic computational mouse model to locate nM to μM fluorophore concentration distributions in simulated mouse organs. PMID:18364555
3D unstructured-mesh radiation transport codes
Morel, J.
1997-12-31
Three unstructured-mesh radiation transport codes are currently being developed at Los Alamos National Laboratory. The first code is ATTILA, which uses an unstructured tetrahedral mesh in conjunction with standard Sn (discrete-ordinates) angular discretization, standard multigroup energy discretization, and linear-discontinuous spatial differencing. ATTILA solves the standard first-order form of the transport equation using source iteration in conjunction with diffusion-synthetic acceleration of the within-group source iterations. DANTE is designed to run primarily on workstations. The second code is DANTE, which uses a hybrid finite-element mesh consisting of arbitrary combinations of hexahedra, wedges, pyramids, and tetrahedra. DANTE solves several second-order self-adjoint forms of the transport equation including the even-parity equation, the odd-parity equation, and a new equation called the self-adjoint angular flux equation. DANTE also offers three angular discretization options: $S{_}n$ (discrete-ordinates), $P{_}n$ (spherical harmonics), and $SP{_}n$ (simplified spherical harmonics). DANTE is designed to run primarily on massively parallel message-passing machines, such as the ASCI-Blue machines at LANL and LLNL. The third code is PERICLES, which uses the same hybrid finite-element mesh as DANTE, but solves the standard first-order form of the transport equation rather than a second-order self-adjoint form. DANTE uses a standard $S{_}n$ discretization in angle in conjunction with trilinear-discontinuous spatial differencing, and diffusion-synthetic acceleration of the within-group source iterations. PERICLES was initially designed to run on workstations, but a version for massively parallel message-passing machines will be built. The three codes will be described in detail and computational results will be presented.
Finite element radiation transport in one dimension
Painter, J.F.
1997-05-09
A new physics package solves radiation transport equations in one space dimension, multiple energy groups and directions. A discontinuous finite element method discretizes radiation intensity with respect to space and angle, and a continuous finite element method discretizes electron temperature `in space. A splitting method solves the resulting linear equations. This is a one-dimensional analog of Kershaw and Harte`s two-dimensional package. This package has been installed in a two-dimensional inertial confinement fusion code, and has given excellent results for both thermal waves and highly directional radiation. In contrast, the traditional discrete ordinate and spherical harmonic methods show less accurate results in both cases.
Coupled Deterministic-Monte Carlo Transport for Radiation Portal Modeling
Smith, Leon E.; Miller, Erin A.; Wittman, Richard S.; Shaver, Mark W.
2008-01-14
Radiation portal monitors are being deployed, both domestically and internationally, to detect illicit movement of radiological materials concealed in cargo. Evaluation of the current and next generations of these radiation portal monitor (RPM) technologies is an ongoing process. 'Injection studies' that superimpose, computationally, the signature from threat materials onto empirical vehicle profiles collected at ports of entry, are often a component of the RPM evaluation process. However, measurement of realistic threat devices can be both expensive and time-consuming. Radiation transport methods that can predict the response of radiation detection sensors with high fidelity, and do so rapidly enough to allow the modeling of many different threat-source configurations, are a cornerstone of reliable evaluation results. Monte Carlo methods have been the primary tool of the detection community for these kinds of calculations, in no small part because they are particularly effective for calculating pulse-height spectra in gamma-ray spectrometers. However, computational times for problems with a high degree of scattering and absorption can be extremely long. Deterministic codes that discretize the transport in space, angle, and energy offer potential advantages in computational efficiency for these same kinds of problems, but the pulse-height calculations needed to predict gamma-ray spectrometer response are not readily accessible. These complementary strengths for radiation detection scenarios suggest that coupling Monte Carlo and deterministic methods could be beneficial in terms of computational efficiency. Pacific Northwest National Laboratory and its collaborators are developing a RAdiation Detection Scenario Analysis Toolbox (RADSAT) founded on this coupling approach. The deterministic core of RADSAT is Attila, a three-dimensional, tetrahedral-mesh code originally developed by Los Alamos National Laboratory, and since expanded and refined by Transpire, Inc. [1
Radiative heat transport instability in ICF plasmas
NASA Astrophysics Data System (ADS)
Rozmus, W.; Bychenkov, V. Yu.
2015-11-01
A laser produced high-Z plasma in which an energy balance is achieved due to radiation losses and radiative heat transfer supports ion acoustic wave instability. A linear dispersion relation is derived and instability is compared to the radiation cooling instability. This instability develops in the wide range of angles and wavenumbers with the typical growth rate on the order of cs/LT (cs is the sound speed, LT is the temperature scale length). In addition to radiation dominated systems, a similar thermal transport driven ion acoustic instability was found before in plasmas where the thermal transport coefficient depends on electron density. However, under conditions of indirect drive ICF experiments the driving term for the instability is the radiative heat flux and in particular, the density dependence of the radiative heat conductivity. A specific example of thermal Bremsstrahlung radiation source has been considered corresponding to a thermal conductivity coefficient that is inversely proportional to the square of local particle density. In the nonlinear regime this instability may lead to plasma jet formation and anisotropic x-ray generation.
George, G. L.; Olsher, R. H.; Seagraves, D. T.
2002-01-01
MCNP-4C1 was used to perform the shielding design for the new Central Health Physics Calibration Facility (CHPCF) at Los Alamos National Laboratory (LANL). The problem of shielding the facility was subdivided into three separate components: (1) Transmission; (2) Skyshine; and (3) Maze Streaming/ Transmission. When possible, actual measurements were taken to verify calculation results. The comparison of calculation versus measurement results shows excellent agreement for neutron calculations. For photon comparisons, calculations resulted in conservative estimates of the Effective Dose Equivalent (EDE) compared to measured results. This disagreement in the photon measurements versus calculations is most likely due to several conservative assumptions regarding shield density and composition. For example, reinforcing steel bars (Rebar) in the concrete shield walls were not included in the shield model.
Space Radiation Transport Methods Development
NASA Technical Reports Server (NTRS)
Wilson, J. W.; Tripathi, R. K.; Qualls, G. D.; Cucinotta, F. A.; Prael, R. E.; Norbury, J. W.; Heinbockel, J. H.; Tweed, J.
2002-01-01
Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be utilized in the final design as verification of the deterministic method optimized design.
Radiation Transport in Type IA Supernovae
Eastman, R
1999-11-16
It has been said more than once that the critical link between explosion models and observations is the ability to accurately simulate cooling and radiation transport in the expanding ejecta of Type Ia supernovae. It is perhaps frustrating to some of the theorists who study explosion mechanisms, and to some of the observers too, that more definitive conclusions have not been reached about the agreement, or lack thereof, between various Type Ia supernova models and the data. Although claims of superlative accuracy in transport simulations are sometimes made, I will argue here that there are outstanding issues of critical importance and in need of addressing before radiation transport calculations are accurate enough to discriminate between subtly different explosion models.
Monte Carlo radiation transport: A revolution in science
Hendricks, J.
1993-04-01
When Enrico Fermi, Stan Ulam, Nicholas Metropolis, John von Neuman, and Robert Richtmyer invented the Monte Carlo method fifty years ago, little could they imagine the far-flung consequences, the international applications, and the revolution in science epitomized by their abstract mathematical method. The Monte Carlo method is used in a wide variety of fields to solve exact computational models approximately by statistical sampling. It is an alternative to traditional physics modeling methods which solve approximate computational models exactly by deterministic methods. Modern computers and improved methods, such as variance reduction, have enhanced the method to the point of enabling a true predictive capability in areas such as radiation or particle transport. This predictive capability has contributed to a radical change in the way science is done: design and understanding come from computations built upon experiments rather than being limited to experiments, and the computer codes doing the computations have become the repository for physics knowledge. The MCNP Monte Carlo computer code effort at Los Alamos is an example of this revolution. Physicians unfamiliar with physics details can design cancer treatments using physics buried in the MCNP computer code. Hazardous environments and hypothetical accidents can be explored. Many other fields, from underground oil well exploration to aerospace, from physics research to energy production, from safety to bulk materials processing, benefit from MCNP, the Monte Carlo method, and the revolution in science.
RADIATION TRANSPORT FOR EXPLOSIVE OUTFLOWS: OPACITY REGROUPING
Wollaeger, Ryan T.; Van Rossum, Daniel R. E-mail: daan@flash.uchicago.edu
2014-10-01
Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure ''opacity regrouping''. Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ∼10% less than that calculated using PHOENIX.
Radiation Transport for Explosive Outflows: Opacity Regrouping
NASA Astrophysics Data System (ADS)
Wollaeger, Ryan T.; van Rossum, Daniel R.
2014-10-01
Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure "opacity regrouping." Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ~10% less than that calculated using PHOENIX.
Radiation-shielding transport and storage container
Baatz, H.; Ritsscher, D.; Wriegt, J.
1985-07-09
A container for the storage and transportation of radioactive material comprises an elongated vessel adapted to receive the material and has a wall thickness and composition attenuating radioactive transmission therefrom. The vessel has an open end formed with an annular thickened portion defining a mouth communicating with the interior of the vessel; a radiation-shielding cover received in the mouth and having a plug-forming portion juxtaposed with a complimentary seat-forming portion of the vessel at the mouth, and a flange extending outwardly from the plug-forming portion, the vessel is provided with a wall bore communicating at one end with the interior of the vessel and terminating at its opposite end within the outline of the radiation-shielding cover, the radiation-shielding cover is provided with a connecting bore registering with the wall bore; an obturating element received in the connecting bore and adapted to block the wall bore; and a further cover secured directly to the vessel outwardly of the radiation-shielding cover and overlying the wall bore and the radiation-shielding cover.
Monte Carlo radiation transport¶llelism
Cox, L. J.; Post, S. E.
2002-01-01
This talk summarizes the main aspects of the LANL ASCI Eolus project and its major unclassified code project, MCNP. The MCNP code provide a state-of-the-art Monte Carlo radiation transport to approximately 3000 users world-wide. Almost all hardware platforms are supported because we strictly adhere to the FORTRAN-90/95 standard. For parallel processing, MCNP uses a mixture of OpenMp combined with either MPI or PVM (shared and distributed memory). This talk summarizes our experiences on various platforms using MPI with and without OpenMP. These platforms include PC-Windows, Intel-LINUX, BlueMountain, Frost, ASCI-Q and others.
Transport methods and interactions for space radiations
NASA Technical Reports Server (NTRS)
Wilson, John W.; Townsend, Lawrence W.; Schimmerling, Walter S.; Khandelwal, Govind S.; Khan, Ferdous S.; Nealy, John E.; Cucinotta, Francis A.; Simonsen, Lisa C.; Shinn, Judy L.; Norbury, John W.
1991-01-01
A review of the program in space radiation protection at the Langley Research Center is given. The relevant Boltzmann equations are given with a discussion of approximation procedures for space applications. The interaction coefficients are related to solution of the many-body Schroedinger equation with nuclear and electromagnetic forces. Various solution techniques are discussed to obtain relevant interaction cross sections with extensive comparison with experiments. Solution techniques for the Boltzmann equations are discussed in detail. Transport computer code validation is discussed through analytical benchmarking, comparison with other codes, comparison with laboratory experiments and measurements in space. Applications to lunar and Mars missions are discussed.
Effects of Nuclear Interactions in Space Radiation Transport
NASA Technical Reports Server (NTRS)
Lin, Zi-Wei; Barghouty, A. F.
2004-01-01
Space radiation transport codes have been developed to calculate radiation effects behind materials in human missions to the Moon, Mars or beyond. We study how nuclear fragmentation processes affect predictions from such radiation transport codes. In particular, we investigate the effects of fragmentation cross sections at different energies on fluxes, dose and dose-equivalent from galactic cosmic rays behind typical shielding materials.
Effects of Nuclear Interactions in Space Radiation Transport
NASA Technical Reports Server (NTRS)
Lin, Zi-Wei; Barghouty, A. F.
2005-01-01
Space radiation transport codes have been developed to calculate radiation effects behind materials in human mission to the Moon, Mars or beyond. We study how nuclear fragmentation processes affect predictions from such radiation transport codes. In particular, we investigate the effects of fragmentation cross sections at different energies on fluxes, dose and dose-equivalent from galactic cosmic rays behind typical shielding materials.
High Performance Radiation Transport Simulations on TITAN
Baker, Christopher G; Davidson, Gregory G; Evans, Thomas M; Hamilton, Steven P; Jarrell, Joshua J; Joubert, Wayne
2012-01-01
In this paper we describe the Denovo code system. Denovo solves the six-dimensional, steady-state, linear Boltzmann transport equation, of central importance to nuclear technology applications such as reactor core analysis (neutronics), radiation shielding, nuclear forensics and radiation detection. The code features multiple spatial differencing schemes, state-of-the-art linear solvers, the Koch-Baker-Alcouffe (KBA) parallel-wavefront sweep algorithm for inverting the transport operator, a new multilevel energy decomposition method scaling to hundreds of thousands of processing cores, and a modern, novel code architecture that supports straightforward integration of new features. In this paper we discuss the performance of Denovo on the 10--20 petaflop ORNL GPU-based system, Titan. We describe algorithms and techniques used to exploit the capabilities of Titan's heterogeneous compute node architecture and the challenges of obtaining good parallel performance for this sparse hyperbolic PDE solver containing inherently sequential computations. Numerical results demonstrating Denovo performance on early Titan hardware are presented.
Radiation heat transport in disordered media
NASA Astrophysics Data System (ADS)
Strieder, William
Radiation heat transport through the internal void spaces of particle beds, fiber beds, packed beds, reactors and porous media with opaque, diffusely reflecting, gray body surfaces and large solid dimensions ( πd/ λth > 100) is considered. A the local particle surface radiosity, differential view factor and solid temperature, permits a rigorous solution of the dependent, long range multiple scattering problem. The conductivity results, applied to a bed of randomly overlapping spheres, agree exactly with pseudohomogeneous results in both the isotropic and anisotropic scattering limits and shed rigorous light on the anisotropic phase function expansion theory. Explicit calculations, performed for several other standard packings, e.g. fiber beds, exhibit a parallll upper and series lower bound over the various particle shapes and dispersion structures. Results show that an empirical equation first suggested by Vortmeyer (German Chem. Engng, 3, (1980) 124-137), but generalized herein from one P to four P0, P1, P2, P3 coefficients, which vary substantially with the various industrial packings, will provide a suitable generalization of the emissivity factor of krad for engineering conductivity modeling of radiation heat transport.
Rare event simulation in radiation transport
Kollman, C.
1993-10-01
This dissertation studies methods for estimating extremely small probabilities by Monte Carlo simulation. Problems in radiation transport typically involve estimating very rare events or the expected value of a random variable which is with overwhelming probability equal to zero. These problems often have high dimensional state spaces and irregular geometries so that analytic solutions are not possible. Monte Carlo simulation must be used to estimate the radiation dosage being transported to a particular location. If the area is well shielded the probability of any one particular particle getting through is very small. Because of the large number of particles involved, even a tiny fraction penetrating the shield may represent an unacceptable level of radiation. It therefore becomes critical to be able to accurately estimate this extremely small probability. Importance sampling is a well known technique for improving the efficiency of rare event calculations. Here, a new set of probabilities is used in the simulation runs. The results are multiple by the likelihood ratio between the true and simulated probabilities so as to keep the estimator unbiased. The variance of the resulting estimator is very sensitive to which new set of transition probabilities are chosen. It is shown that a zero variance estimator does exist, but that its computation requires exact knowledge of the solution. A simple random walk with an associated killing model for the scatter of neutrons is introduced. Large deviation results for optimal importance sampling in random walks are extended to the case where killing is present. An adaptive ``learning`` algorithm for implementing importance sampling is given for more general Markov chain models of neutron scatter. For finite state spaces this algorithm is shown to give with probability one, a sequence of estimates converging exponentially fast to the true solution.
D. Kraig; W. A. Treadaway; R. J. Wechsler
1999-10-01
Recent advances and commercialization of electret-ion-chamber (EIC) technology for photon measurements prompted us to consider EKs as a replacement for our TLD system. After laboratory tests indicated that both systems gave adequate results for controlled exposures, throughout 1998 we conducted field tests with paired TLDs and EICS, in LANL technical areas and in public areas. We had approximately 30 paired sampling sites at Area G. At each sampling site, we deployed three TLDs and three EICS. The EICS were contained in air-tight jars, each of which was placed in a Tyvek envelope and hung about 1 m above the ground. The dosimeters were read (and, if necessary, replaced) every three months. At the sites outside Area G, the TLD readings for the first two quarters were statistically significantly higher than those of the EICS: group average exposures were 38 and 36, compared with 33 mR (both quarters) for the EICS; during quarter 3, the EIC average (40 mR) was higher than the TLD average (34 mR); and during quarter 4, the two systems were statistically the same: EIC = 42, TLD = 41 with a p-value of 0.61. We are still evaluating these differences and performing additional laboratory studies to determine causes. At the Area G sites,we noticed that several of the TLDs gave much higher readings than their co-located EICS; we believe that the TLDs were over-responding by {approx}50% to the low-energy (60-keV) gamma radiation associated with {sup 241}Am, whereas the EICS were responding accurately. We conclude that EICS are more accurate at a wide range of gamma energies and are preferable to TLDs in environments where a significant fraction of the photons are low energy.
THE MCNPX MONTE CARLO RADIATION TRANSPORT CODE
WATERS, LAURIE S.; MCKINNEY, GREGG W.; DURKEE, JOE W.; FENSIN, MICHAEL L.; JAMES, MICHAEL R.; JOHNS, RUSSELL C.; PELOWITZ, DENISE B.
2007-01-10
MCNPX (Monte Carlo N-Particle eXtended) is a general-purpose Monte Carlo radiation transport code with three-dimensional geometry and continuous-energy transport of 34 particles and light ions. It contains flexible source and tally options, interactive graphics, and support for both sequential and multi-processing computer platforms. MCNPX is based on MCNP4B, and has been upgraded to most MCNP5 capabilities. MCNP is a highly stable code tracking neutrons, photons and electrons, and using evaluated nuclear data libraries for low-energy interaction probabilities. MCNPX has extended this base to a comprehensive set of particles and light ions, with heavy ion transport in development. Models have been included to calculate interaction probabilities when libraries are not available. Recent additions focus on the time evolution of residual nuclei decay, allowing calculation of transmutation and delayed particle emission. MCNPX is now a code of great dynamic range, and the excellent neutronics capabilities allow new opportunities to simulate devices of interest to experimental particle physics; particularly calorimetry. This paper describes the capabilities of the current MCNPX version 2.6.C, and also discusses ongoing code development.
Validation of comprehensive space radiation transport code
Shinn, J.L.; Simonsen, L.C.; Cucinotta, F.A.
1998-12-01
The HZETRN code has been developed over the past decade to evaluate the local radiation fields within sensitive materials on spacecraft in the space environment. Most of the more important nuclear and atomic processes are now modeled and evaluation within a complex spacecraft geometry with differing material components, including transition effects across boundaries of dissimilar materials, are included. The atomic/nuclear database and transport procedures have received limited validation in laboratory testing with high energy ion beams. The codes have been applied in design of the SAGE-III instrument resulting in material changes to control injurious neutron production, in the study of the Space Shuttle single event upsets, and in validation with space measurements (particle telescopes, tissue equivalent proportional counters, CR-39) on Shuttle and Mir. The present paper reviews the code development and presents recent results in laboratory and space flight validation.
Transport of infrared radiation in cuboidal clouds
NASA Technical Reports Server (NTRS)
HARSHVARDHAN; Weinman, J. A.; Davies, R.
1981-01-01
The transport of infrared radiation in a single cuboidal cloud using a vertical two steam approximation was modeled. The emittance of the top face of the model cloud is always less than that for a plane parallel cloud of the same optical depth. The hemisphere flux escaping from the cloud top has a gradient from the center to the edges which brighten when the cloud is over warmer ground. Cooling rate calculations in the 8 to 13.6 micrometer region show that there is cooling from the sides of the cloud at all levels even when there is heating of the core from the ground below. The radiances exiting from model cuboidal clouds were computed by path integration over the source function obtained with the two stream approximation. It is suggested that the brightness temperature measured from finite clouds will overestimate the cloud top temperature.
Space Radiation Transport Code Development: 3DHZETRN
NASA Technical Reports Server (NTRS)
Wilson, John W.; Slaba, Tony C.; Badavi, Francis F.; Reddell, Brandon D.; Bahadori, Amir A.
2015-01-01
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
Simulating Radiation Transport in Curved Spacetimes
NASA Astrophysics Data System (ADS)
Endeve, Eirik; Hauck, Cory; Xing, Yulong; Cardall, Christian; Mezzacappa, Anthony
2014-03-01
We are developing methods for simulation of radiation transport in systems governed by strong gravity (e.g., neutrino transport in core-collapse supernovae). By employing conservative formulations of the general relativistic Boltzmann equation, we aim to develop methods that are (i) high-order accurate for computational efficiency; (ii) robust in the sense that the phase space density f preserves the maximum principle of the physical model (f ∈ [ 0 , 1 ] for fermions); and (iii) applicable to curvilinear coordinate systems to accommodate curved spacetimes, which result in gravity-induced frequency shift and angular aberration. Our approach is based on the Runge-Kutta discontinuous Galerkin method, which has many attractive properties, including high-order accuracy on a compact stencil. We present the physical model, describe our numerical methods, and show results from implementations in spherical and axial symmetry. Our tests show that the method is high-order accurate and strictly preserves the maximum principle on f. We also demonstrate the ability of our method to accurately include effects of a strong gravitational field.
Transport and radiation in complex LTE mixtures
NASA Astrophysics Data System (ADS)
Janssen, Jesper; Peerenboom, Kim; Suijker, Jos; Gnybida, Mykhailo; van Dijk, Jan
2014-10-01
Complex LTE mixtures are for example encountered in re-entry, welding, spraying and lighting. These mixtures typically contain a rich chemistry in combination with large temperature gradients. LTE conditions are also interesting because they can aid in the validation of NLTE algorithms. An example is the calculation of transport properties. In this work a mercury free high intensity discharge lamp is considered. The investigation focusses on using salts like InI or SnI as a buffer species. By using these species a dominant background gas like mercury is no longer present. As a consequence the diffusion algorithms based on Fick's law are no longer applicable and the Stefan-Maxwell equations must be solved. This system of equations is modified with conservation rules to set a coldspot pressure for saturated species and enforce the mass dosage for unsaturated species. The radiative energy transport is taken into account by raytracing. Quantum mechanical simulations have been used to calculate the potential curves and the transition dipole moments for indium with iodine and tin with iodine. The results of these calculations have been used to predict the quasistatic broadening by iodine. The work was supported by the project SCHELP from the Belgium IWT (Project Number 110003) and the CATRENE SEEL Project (CA502).
An Integrated Radiation Transport Particle-in-Cell Method
NASA Astrophysics Data System (ADS)
Lee, H. J.; Verboncoeur, J. P.; Smith, H. B.; Parker, G. J.; Birdsall, C. K.
2000-10-01
The study of radiation transport is important to understand the basic physics and to calculate the efficiency in a lamp discharge or laser induced plasma. Many models neglect radiation transport effects in evolving the steady state. In this study, we established a basic model to calculate radiation transport, including the effects of nonuniform ground state density and atomic collisions in one dimensional cylindrical and planar geometries. We coupled radiation transport with the self-consistent kinetic particle-in-cell codes, XPDP1 and XPDC1[1]. We treat electrons and ions with a particle-in-cell method, and the neutral ground and excited states with a fluid model to calculate radiation transport and atomic collisions. The steady state result of this model compares well with the solution of Holstein equation[2]. [1] J. P. Verboncoeur, M. V. Alves, V. Vahedi, and C. K. Birdsall, Journal of Computational Physics 104, 321 (1993). [2] T. Holstein, Phys. Rev. 72, 1213 (1947).
Description of Transport Codes for Space Radiation Shielding
NASA Technical Reports Server (NTRS)
Kim, Myung-Hee Y.; Wilson, John W.; Cucinotta, Francis A.
2011-01-01
This slide presentation describes transport codes and their use for studying and designing space radiation shielding. When combined with risk projection models radiation transport codes serve as the main tool for study radiation and designing shielding. There are three criteria for assessing the accuracy of transport codes: (1) Ground-based studies with defined beams and material layouts, (2) Inter-comparison of transport code results for matched boundary conditions and (3) Comparisons to flight measurements. These three criteria have a very high degree with NASA's HZETRN/QMSFRG.
Transport of infrared radiation in cuboidal clouds
NASA Technical Reports Server (NTRS)
Harshvardhan, MR.; Weinman, J. A.; Davies, R.
1981-01-01
The transport of infrared radiation in a single cuboidal cloud is modeled using a variable azimuth two-stream approximation. Computations are made at 10 microns for a Deirmendjian (1969) C-1 water cloud where the single scattering albedo is equal to 0.638 and the asymmetry parameter is 0.865. The results indicate that the emittance of the top face of the model cloud is always less than that for a plane parallel cloud of the same optical depth. The hemispheric flux escaping from the cloud top possesses a gradient from the center to the edges which are warmer when the cloud is over warmer ground. Cooling rate calculations in the 8-13.6 micron region demonstrate that there is cooling out of the sides of the cloud at all levels even when there is heating of the core from the ground below. The radiances exiting from model cuboidal clouds are computed by path integration over the source function obtained with the two-stream approximation. Results indicate that the brightness temperature measured from finite clouds will overestimate the cloud-top temperature.
Radiation energy transport through hydrodynamically evolving slits
NASA Astrophysics Data System (ADS)
Foster, J. M.; Graham, P.; Taylor, M.; Moore, A.; Sorce, C.; Reighard, A.; MacLaren, S.; Young, P.; Glendinning, G.; Blue, B.; Back, C.; Hund, J.
2008-11-01
Radiation transport through enclosed spaces with inwardly moving walls is a key component of the physics of laser-heated hohlraums. It arises in the cavity itself (where inward motion of the wall results in late-time stagnation of dense plasma on the hohlraum axis), and also in the laser-entry and diagnostic holes (where an understanding of hole-closure is important to hohlraum design and the interpretation of diagnostic data). To understand these phenomena better, we have carried out a series of experiments at the Omega laser facility. A laser-heated hohlraum is used to illuminate linear and annular slits machined in samples of solid-density tantalum and low-density, tantalum-oxide foam. Measurements of the transmitted energy are made indirectly (by measuring the temperature rise of a ``calorimeter'' hohlraum) and directly (by measuring the emission from the slit component, using a target in which the calorimeter hohlraum was omitted). The hydrodynamics is investigated by self-emission and absorption (backlighting) x-ray imaging of the closing slits. Simulations (using a 2-D Eulerian hydrocode) reproduce the overall energetics, the detail of the deceleration shock and axial stagnation region at the centre of the slit, and the complex shock interactions that occur at corners of the slits.
Environmental analysis of Lower Pueblo/Lower Los Alamos Canyon, Los Alamos, New Mexico
Ferenbaugh, R.W.; Buhl, T.E.; Stoker, A.K.; Becker, N.M.; Rodgers, J.C.; Hansen, W.R.
1994-12-01
The radiological survey of the former radioactive waste treatment plant site (TA-45), Acid Canyon, Pueblo Canyon, and Los Alamos Canyon found residual contamination at the site itself and in the channel and banks of Acid, Pueblo, and lower Los Alamos Canyons all the way to the Rio Grande. The largest reservoir of residual radioactivity is in lower Pueblo Canyon, which is on DOE property. However, residual radioactivity does not exceed proposed cleanup criteria in either lower Pueblo or lower Los Alamos Canyons. The three alternatives proposed are (1) to take no action, (2) to construct a sediment trap in lower Pueblo Canyon to prevent further transport of residual radioactivity onto San Ildefonso Indian Pueblo land, and (3) to clean the residual radioactivity from the canyon system. Alternative 2, to cleanup the canyon system, is rejected as a viable alternative. Thousands of truckloads of sediment would have to be removed and disposed of, and this effort is unwarranted by the low levels of contamination present. Residual radioactivity levels, under either present conditions or projected future conditions, will not result in significant radiation doses to persons exposed. Modeling efforts show that future transport activity will not result in any residual radioactivity concentrations higher than those already existing. Thus, although construction of a sediment trap in lower Pueblo Canyon is a viable alternative, this effort also is unwarranted, and the no-action alternative is the preferred alternative.
Radiation transport phenomena and modeling - part A: Codes
Lorence, L.J.
1997-06-01
The need to understand how particle radiation (high-energy photons and electrons) from a variety of sources affects materials and electronics has motivated the development of sophisticated computer codes that describe how radiation with energies from 1.0 keV to 100.0 GeV propagates through matter. Predicting radiation transport is the necessary first step in predicting radiation effects. The radiation transport codes that are described here are general-purpose codes capable of analyzing a variety of radiation environments including those produced by nuclear weapons (x-rays, gamma rays, and neutrons), by sources in space (electrons and ions) and by accelerators (x-rays, gamma rays, and electrons). Applications of these codes include the study of radiation effects on electronics, nuclear medicine (imaging and cancer treatment), and industrial processes (food disinfestation, waste sterilization, manufacturing.) The primary focus will be on coupled electron-photon transport codes, with some brief discussion of proton transport. These codes model a radiation cascade in which electrons produce photons and vice versa. This coupling between particles of different types is important for radiation effects. For instance, in an x-ray environment, electrons are produced that drive the response in electronics. In an electron environment, dose due to bremsstrahlung photons can be significant once the source electrons have been stopped.
NASA Astrophysics Data System (ADS)
Ding, M.; Longmire, P.; Vaniman, D.; Broxton, D.; Rearick, M.; Ware, S.; Katzman, D.; Bargar, J.; Webb, S.
2009-12-01
Routine groundwater monitoring conducted in 2005 identified significant chromium(Cr) contamination within the regional aquifer at Los Alamos, New Mexico. Site-specific Cr(III)/Cr(VI) transport characterization is needed (1) to refine conceptual and numerical models; (2) to evaluate long-term stability of Cr(III)/Cr(VI) in the Sandia Canyon groundwater system along the flow path from the source to the regional aquifer; and (3) to forecast natural attenuation and assess remediation possibilities. We have conducted both static (batch) and flowing (column) experiments using water chemistries and geological materials representative of potential flow pathways in Sandia Canyon. We also conducted synchrotron-based X-ray absorption spectroscopy measurements for determining Cr(III)/Cr(VI) speciation and oxidation state of Cr surface complexes in vadose zone and regional aquifer materials. In addition, characterization of materials used in experiments was performed using SEM, XRD, and BET Surface Area Measurement. Results of adsorption and flow-through experiments suggest that vadose zone and regional aquifer materials have different retention capacities for Cr(VI), reflecting large variations in measured adsorption coefficients (Kd values). Additionally, higher Kd values tend to be associated with crushed rocks and reduced particle sizes, where increased surface area and exposure of fresh mineral provides surfaces containing ferrous iron that acts as a chemical reductant. Nevertheless, X-ray Absorption Near Edge Spectroscopy (XANES) measurements suggest that reduction of Cr(VI) to Cr(III) occurred after Cr(VI) containing solution was exposed to alluvial and rock materials, both uncrushed and crushed. The results indicate that natural attenuation of anthropogenic chromium occurs within the Sandia Canyon vadose zone and in the regional aquifer. Findings from this study provide data necessary to reveal the nature and extent of site-specific long-term fate and transport of
Path Toward a Unifid Geometry for Radiation Transport
NASA Technical Reports Server (NTRS)
Lee, Kerry; Barzilla, Janet; Davis, Andrew; Zachmann
2014-01-01
The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their
Path Toward a Unified Geometry for Radiation Transport
NASA Astrophysics Data System (ADS)
Lee, Kerry
The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex CAD models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN (high charge and energy transport code developed by NASA LaRC), are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading
Effects of Nuclear Interactions on Accuracy of Space Radiation Transport
NASA Technical Reports Server (NTRS)
Lin, Zi-Wei; Barghouty, A. F.
2005-01-01
Space radiation risk to astronauts and electronic equipments is one major obstacle in long term human space explorations. Space radiation transport codes have been developed to calculate radiation effects behind materials in human missions to the Moon, Mars or beyond. We study how nuclear fragmentation processes affect the accuracy of predictions from such radiation transport. In particular, we investigate the effects of fragmentation cross sections at different energies on fluxes, dose and dose-equivalent from galactic cosmic rays behind typical shielding materials. These results tell us at what energies nuclear cross sections are the most important for radiation risk evaluations, and how uncertainties in our knowledge about nuclear fragmentations relate to uncertainties in space transport predictions.
Near-field radiative thermal transport: From theory to experiment
Song, Bai Fiorino, Anthony; Meyhofer, Edgar; Reddy, Pramod
2015-05-15
Radiative thermal transport via the fluctuating electromagnetic near-field has recently attracted increasing attention due to its fundamental importance and its impact on a range of applications from data storage to thermal management and energy conversion. After a brief historical account of radiative thermal transport, we summarize the basics of fluctuational electrodynamics, a theoretical framework for the study of radiative heat transfer in terms of thermally excited propagating and evanescent electromagnetic waves. Various approaches to modeling near-field thermal transport are briefly discussed, together with key results and proposals for manipulation and utilization of radiative heat flow. Subsequently, we review the experimental advances in the characterization of both near-field heat flow and energy density. We conclude with remarks on the opportunities and challenges for future explorations of radiative heat transfer at the nanoscale.
Description of transport codes for space radiation shielding.
Kim, Myung-Hee Y; Wilson, John W; Cucinotta, Francis A
2012-11-01
Exposure to ionizing radiation in the space environment is one of the hazards faced by crews in space missions. As space radiations traverse spacecraft, habitat shielding, or tissues, their energies and compositions are altered by interactions with the shielding. Modifications to the radiation fields arise from atomic interactions of charged particles with orbital electrons and nuclear interactions leading to projectile and target fragmentation, including secondary particles such as neutrons, protons, mesons, and nuclear recoils. The transport of space radiation through shielding can be simulated using Monte Carlo techniques or deterministic solutions of the Boltzmann equation. To determine shielding requirements and to resolve radiation constraints for future human missions, the shielding evaluation of a spacecraft concept is required as an early step in the design process. To do this requires (1) accurate knowledge of space environmental models to define the boundary condition for transport calculations, (2) transport codes with detailed shielding and body geometry models to determine particle transmission into areas of internal shielding and at each critical body organ, and (3) the assessment of organ dosimetric quantities and biological risks by applying the corresponding response models for space radiation against the particle spectra that have been accurately determined from the transport code. This paper reviews current transport codes and analyzes their accuracy through comparison to laboratory and spaceflight data. This paper also introduces a probabilistic risk assessment approach for the evaluation of radiation shielding. PMID:23032892
NASA Astrophysics Data System (ADS)
Zhao, J. M.; Tan, J. Y.; Liu, L. H.
2012-02-01
Light transport in graded index media follows a curved trajectory determined by Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.
Peterson, D.L.; Bowers, R.L.; Lebeda, C.F.; Matuska, W.; Benage, J.; Idzorek, G.; Oona, H.; Stokes, J.; Roderick, N.F.
1995-09-01
Two experiments, PegI-41, conducted on the Los Alamos Pegasus I capacitor bank, and PegII-25, on the Pegasus II bank, consisted of the implosions of 13 mg (nominal), 5 cm radius, 2 cm high thin cylindrical aluminum foils resulting in soft x-ray radiation pulses from the plasma thermalization on axis. The implosions were conducted in direct-drive (no intermediate switching) mode with peak currents of about 4 MA and 5 MA respectively, and implosion times of about 2.5 {micro}s and 2.0 {micro}s. A radiation yield of about 250 kJ was measured for PegII-25. The purpose of these experiments was to examine the physics of the implosion and relate this physics to the production of the radiation pulse and to provide detailed experimental data which could be compared with 2-D radiation-magnetohydrodynamic (RMHD) simulations. Included in the experimental diagnostic suites were faraday rotation and dB/dt current measurements, a visible framing camera, an x-ray stripline camera, time-dependent spectroscopy, bolometers and XRD`S. A comparison of the results from these experiments shows agreement with 2-D simulation results in the instability development, current, and radiation pulse data, including the pulsewidth, shape, peak power and total radiation yield as measured by bolometry. Instabilities dominate the behavior of the implosion and largely determine the properties of the resulting radiation pulse. The 2-D simulations can be seen to be an important tool in understanding the implosion physics.
Improved geometry representations for Monte Carlo radiation transport.
Martin, Matthew Ryan
2004-08-01
ITS (Integrated Tiger Series) permits a state-of-the-art Monte Carlo solution of linear time-integrated coupled electron/photon radiation transport problems with or without the presence of macroscopic electric and magnetic fields of arbitrary spatial dependence. ITS allows designers to predict product performance in radiation environments.
TORUS: Radiation transport and hydrodynamics code
NASA Astrophysics Data System (ADS)
Harries, Tim
2014-04-01
TORUS is a flexible radiation transfer and radiation-hydrodynamics code. The code has a basic infrastructure that includes the AMR mesh scheme that is used by several physics modules including atomic line transfer in a moving medium, molecular line transfer, photoionization, radiation hydrodynamics and radiative equilibrium. TORUS is useful for a variety of problems, including magnetospheric accretion onto T Tauri stars, spiral nebulae around Wolf-Rayet stars, discs around Herbig AeBe stars, structured winds of O supergiants and Raman-scattered line formation in symbiotic binaries, and dust emission and molecular line formation in star forming clusters. The code is written in Fortran 2003 and is compiled using a standard Gnu makefile. The code is parallelized using both MPI and OMP, and can use these parallel sections either separately or in a hybrid mode.
Relativistic radiation transport in dispersive media
Kichenassamy, S.; Krikorian, R.A.
1985-10-15
A general-relativistic radiative transfer equation in an isotropic, weakly absorbing, nonmagnetized dispersive medium is derived using the kinetic-theoretical approach and the relativistic Hamiltonian theory of geometrical optics in those media. It yields the generally accepted classical equation in the special-relativistic approximation and in stationary conditions. The influence of the gravitational field and of space-time variations of the refractive index n on the radiation distribution is made explicit in the case of spherical symmetry.
Radiation transport Part B: Applications with examples
Beutler, D.E.
1997-06-01
In the previous sections Len Lorence has described the need, theory, and types of radiation codes that can be applied to model the results of radiation effects tests or working environments for electronics. For the rest of this segment, the author will concentrate on the specific ways the codes can be used to predict device response or analyze radiation test results. Regardless of whether one is predicting responses in a working or test environment, the procedures are virtually the same. The same can be said for the use of 1-, 2-, or 3-dimensional codes and Monte Carlo or discrete ordinates codes. No attempt is made to instruct the student on the specifics of the code. For example, the author will not discuss the details, such as the number of meshes, energy groups, etc. that are appropriate for a discrete ordinates code. For the sake of simplicity, he will restrict himself to the 1-dimensional code CEPXS/ONELD. This code along with a wide variety of other radiation codes can be obtained form the Radiation Safety Information Computational Center (RSICC) for a nominal handling fee.
Environmental surveillance at Los Alamos during 1994
1996-07-01
This report describes environmental monitoring activities at Los Alamos National Laboratory for 1994. Data were collected to assess external penetrating radiation, airborne emissions, liquid effluents, radioactivity of environmental materials and food stuffs, and environmental compliance.
Radiation transport in ultrafast heated high Z solid targets
NASA Astrophysics Data System (ADS)
Paraschiv, Ioana; Sentoku, Yasuhiko; Mancini, Roberto; Johzaki, Tomoyuki
2013-10-01
Ultra-intense laser-target interactions generate hot, dense, and radiating plasmas, especially in the case of high-Z target materials. In order to evaluate the effect of radiation and its transport on the laser-produced plasmas we have developed a radiation transport (RT) code and implemented it in a collisional particle-in-cell code, PICLS. The code uses a database of emissivities and opacities as functions of photon frequency, created for given densities and temperatures by the non-equilibrium, collisional-radiative atomic kinetics 0-D code FLYCHK together with its postprocessor FLYSPECTRA. Using the two-dimensional RT-PICLS code we have studied the X-ray transport in an ultrafast heated copper target, the X-ray conversion efficiency, and the exchange of energy between the radiation field and the target. The details of these results obtained from the implementation of the radiation transport model into the PICLS calculations will be reported in this presentation. Work supported by the DOE Office of Science grant no. DE-SC0008827 and by the NNSA/DOE grant no. DE-FC52-06NA27616.
A new, coupled transport-diffusion method for radiative transfer calculations
Wollaber, A. B.; Warsa, J. S.
2013-07-01
We derive and present a new frequency- and angle-integrated low-order system of equations designed to enhance the accuracy of a coupled, high-order (transport) solution of the thermal radiative transfer equations. In particular, our new low-order system is designed to use intensity-weighted opacities and anisotropic diffusion coefficients generated by a solution of the Implicit Monte Carlo (IMC) equations in order to predict the spatial dependence of the material temperature and radiation energies in the ensuing time cycle. The predicted temperature solution can then be exploited to generate appropriately time-centered opacities, specific heats, and Planck emission spectra for the upcoming IMC solution. Additionally, the relatively inexpensive solution of the low-order system can be iteratively solved to recommend an adaptive time step size before the IMC solution is computed. A test implementation has been implemented using existing software available from the Jayenne and Capsaicin projects at Los Alamos National Laboratory. We present initial results from a new driver code that has integrated these stochastic and deterministic software packages. (authors)
Comparison of Space Radiation Calculations from Deterministic and Monte Carlo Transport Codes
NASA Technical Reports Server (NTRS)
Adams, J. H.; Lin, Z. W.; Nasser, A. F.; Randeniya, S.; Tripathi, r. K.; Watts, J. W.; Yepes, P.
2010-01-01
The presentation outline includes motivation, radiation transport codes being considered, space radiation cases being considered, results for slab geometry, results from spherical geometry, and summary. ///////// main physics in radiation transport codes hzetrn uprop fluka geant4, slab geometry, spe, gcr,
Method for modeling radiative transport in luminescent particulate media.
Hughes, Michael D; Borca-Tasciuc, Diana-Andra; Kaminski, Deborah A
2016-04-20
Modeling radiative transport in luminescent particulate media is important to a variety of applications, from biomedical imaging to solar power harvesting. When absorption and scattering from individual particles must be considered, the description of radiative transport is not straightforward. For large particles and interparticle spacing, geometrical optics can be employed. However, this approach requires accurate knowledge of several particle properties, such as index of refraction and absorption coefficient, along with particle geometry and positioning. Because the determination of these variables is often nontrivial, we developed an approach for modeling radiative transport in such media, which combines two simple experiments with Monte Carlo simulations to determine the particle extinction coefficient (Γ) and the probability of absorption of light by a particle (P_{A}). The method is validated on samples consisting of luminescent phosphor powder dispersed in a silicone matrix. PMID:27140095
The essential linearity of thermal radiation net transport
NASA Astrophysics Data System (ADS)
Barker, J. J.
1991-06-01
Thermal radiation net transport is proportional to the difference in the fourth powers of the absolute temperatures of the radiator and the receiver. During a study of heat transfer mechanisms in the underground storage tanks for the radioactive wastes at the Hanford Site in Washington, it was discovered that the fourth power difference is approximated by a linear function of the temperature difference between the radiator and the receiver, with astounding fidelity. This linearity makes satisfying boundary conditions for analytical and numerical calculations much simpler and simplifies comparisons with other modes of heat transfer. The paper derives the relationship and illustrates its use by demonstrating that radiation transports more heat through the air space in the tank than does natural convection.
Radiative transport-based frequency-domain fluorescence tomography
NASA Astrophysics Data System (ADS)
Joshi, Amit; Rasmussen, John C.; Sevick-Muraca, Eva M.; Wareing, Todd A.; McGhee, John
2008-04-01
We report the development of radiative transport model-based fluorescence optical tomography from frequency-domain boundary measurements. The coupled radiative transport model for describing NIR fluorescence propagation in tissue is solved by a novel software based on the established Attila™ particle transport simulation platform. The proposed scheme enables the prediction of fluorescence measurements with non-contact sources and detectors at a minimal computational cost. An adjoint transport solution-based fluorescence tomography algorithm is implemented on dual grids to efficiently assemble the measurement sensitivity Jacobian matrix. Finally, we demonstrate fluorescence tomography on a realistic computational mouse model to locate nM to µM fluorophore concentration distributions in simulated mouse organs.
Discontinuous Galerkin Methods for Neutrino Radiation Transport
NASA Astrophysics Data System (ADS)
Endeve, Eirik; Hauck, Cory; Xing, Yulong; Mezzacappa, Anthony
2015-04-01
We are developing new computational methods for simulation of neutrino transport in core-collapse supernovae, which is challenging since neutrinos evolve from being diffusive in the proto-neutron star to nearly free streaming in the critical neutrino heating region. To this end, we consider conservative formulations of the Boltzmann equation, and aim to develop robust, high-order accurate methods. Runge-Kutta discontinuous Galerkin (DG) methods, offer several attractive properties, including (i) high-order accuracy on a compact stencil and (ii) correct asymptotic behavior in the diffusion limit. We have recently developed a new DG method for the advection part for the transport solve, which is high-order accurate and strictly preserves the physical bounds of the distribution function; i.e., f ∈ [ 0 , 1 ] . We summarize the main ingredients of our bound-preserving DG method and discuss ongoing work to include neutrino-matter interactions in the scheme. Research sponsored in part by Oak Ridge National Laboratory, managed by UT-Battelle, LLC for the U. S. Department of Energy
New Algorithms for Large-scale 3D Radiation Transport
NASA Astrophysics Data System (ADS)
Lentz, Eric J.
2009-05-01
Radiation transport is critical not only for analysis of astrophysical objects but also for the dynamical transport of energy within. Increased fidelity and dimensionality of the other components of such models requires a similar improvement in the radiation transport. Modern astrophysical simulations can be large enough that the values for a single variable for the entire computational domain cannot be stored on a single compute node. The natural solution is to decompose the physical domain into pieces with each node responsible for a single sub-domain. Using localized plus "ghost" zone data works well for problems like explicit hydrodynamics or nuclear reaction networks with modest impact from inter-process communication. Unfortunately, radiation transport is an inherently non-local process that couples the entire model domain together and efficient algorithms are needed to conquer this problem. In this poster, I present the early development of a new parallel, 3-D transport code using ray tracing to formally solve the transport equation across numerically decomposed domains. The algorithm model takes advantage of one-sided communication to develop a scalable, parallel formal solver. Other aspects and future direction of the parallel code development such as scalability and the inclusion of scattering will also be discussed.
Diagnosing ocean energy transports from earth radiation budget measurements
NASA Technical Reports Server (NTRS)
Sohn, Byung-Ju; Smith, Eric A.
1992-01-01
The maximum energy production (MEP) principle suggested by Paltridge (1975) is applied to separate the satellite-inferred required total transports into the atmospheric and the oceanic components within a two-dimensional (2D) framework. For this purpose, the required 2D energy transports (Sohn and Smith, 1991) are imposed on Paltridge's energy balance model which is then solved as a variational problem. The results provide separated atmospheric and oceanic transports on a 2D basis such that the total divergence is equal to the net radiation measured from a satellite.
Radiation transport on unstructured mesh with parallel computers
Fan, W.C.; Drumm, C.R.
2000-07-01
This paper summarizes the developmental work on a deterministic transport code that provides multidimensional radiation transport capabilities on an unstructured mesh. The second-order form of the Boltzmann transport equation is solved utilizing the discrete ordinates angular differencing and the Galerkin finite element spatial differencing. The discretized system, which couples the spatial-angular dependence, is solved simultaneously using a parallel conjugate-gradient (CG) iterative solver. This approach eliminates the need for the conventional inner iterations over the discrete directions and is well-suited for massively parallel computers.
RADIATION EFFECTS ON TRANSPORT AND BUBBLE FORMATION IN SILICATE GLASSES
The objective of the research is to discover the molecular details of chemistry induced by -, γ-, and neutron-irradiation of silicate glasses. The ionization and ballistic effects of radiation will be studied from the viewpoint of defect formation and transport properties. D...
A space radiation transport method development
NASA Technical Reports Server (NTRS)
Wilson, J. W.; Tripathi, R. K.; Qualls, G. D.; Cucinotta, F. A.; Prael, R. E.; Norbury, J. W.; Heinbockel, J. H.; Tweed, J.
2004-01-01
Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.
A space radiation transport method development.
Wilson, J W; Tripathi, R K; Qualls, G D; Cucinotta, F A; Prael, R E; Norbury, J W; Heinbockel, J H; Tweed, J
2004-01-01
Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. PMID:15880919
Numerical solution of the radiation transport equation in disk geometry
NASA Technical Reports Server (NTRS)
Spagna, George F., Jr.; Leung, Chun Ming
1987-01-01
An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.
Radiation Belt Transport Driven by Solar Wind Dynamic Pressure Fluctuations
NASA Astrophysics Data System (ADS)
Kress, B. T.; Hudson, M. K.; Ukhorskiy, A. Y.; Mueller, H.
2012-12-01
The creation of the Earth's outer zone radiation belts is attributed to earthward transport and adiabatic acceleration of electrons by drift-resonant interactions with electromagnetic fluctuations in the magnetosphere. Three types of radial transport driven by solar wind dynamic pressure fluctuations that have been identified are: (1) radial diffusion [Falthammer, 1965], (2) significant changes in the phase space density radial profile due to a single or few ULF drift-resonant interactions [Ukhorskiy et al., 2006; Degeling et al., 2008], and (3) shock associated injections of radiation belt electrons occurring in less than a drift period [Li et al., 1993]. A progress report will be given on work to fully characterize different forms of radial transport and their effect on the Earth's radiation belts. The work is being carried out by computing test-particle trajectories in electric and magnetic fields from a simple analytic ULF field model and from global MHD simulations of the magnetosphere. Degeling, A. W., L. G. Ozeke, R. Rankin, I. R. Mann, and K. Kabin (2008), Drift resonant generation of peaked relativistic electron distributions by Pc 5 ULF waves, textit{J. Geophys. Res., 113}, A02208, doi:10.1029/2007JA012411. Fälthammar, C.-G. (1965), Effects of Time-Dependent Electric Fields on Geomagnetically Trapped Radiation, J. Geophys. Res., 70(11), 2503-2516, doi:10.1029/JZ070i011p02503. Li, X., I. Roth, M. Temerin, J. R. Wygant, M. K. Hudson, and J. B. Blake (1993), Simulation of the prompt energization and transport of radiation belt particles during the March 24, 1991 SSC, textit{Geophys. Res. Lett., 20}(22), 2423-2426, doi:10.1029/93GL02701. Ukhorskiy, A. Y., B. J. Anderson, K. Takahashi, and N. A. Tsyganenko (2006), Impact of ULF oscillations in solar wind dynamic pressure on the outer radiation belt electrons, textit{Geophys. Res. Lett., 33}(6), L06111, doi:10.1029/2005GL024380.
Instabilities and the transport of polarized astrophysical maser radiation
NASA Technical Reports Server (NTRS)
Wallin, Bradley K.; Watson, William D.
1995-01-01
Time-dependent, radiative instabilities in the creation and transport of polarized astrophysical maser radiation in the presence of a magnetic field are calculated. The instabilities are similar to and occur under the same conditions as those found previously by Scappaticci & Watson for unpolarized maser radiation. The common limits in which the Zeeman splitting is much greater than, and much less than, the spectral line breadths are both considered in the current investigation. The resulting fluctuations in the emergent radiation are potentially relevant for the OH 1665 MHz masers which have been reported to fluctuate on timescales of 1000 s and which tend to be polarized. Instabilities occur and alter the transport of maser radiation only under a quite limited range of conditions. In particular, we find here that the instabilities do not occur for conditions that are appropriate for astrophysical masers with small Zeeman splittings such as the SiO and H2O masers. The time-independent, numerical solutions to the GKK equations of radiative transfer that have been obtained in previous investigations are thus essentially always valid except within a narrow range of conditions relevant for the OH masers.
Thermal radiation transport in a fluidized dry water system
NASA Astrophysics Data System (ADS)
Taylan, Onur; Berberoglu, Halil
2013-05-01
This paper reports a numerical study on thermal radiation transport in a novel material called dry water. Dry water is a water-in-air inverse foam which consists of micrometer-sized water droplets encapsulated by hydrophobic fumed-silica nanoparticles. First, the size distribution of dry water was measured using a particle size analyzer. Then, the radiation characteristics of dry water were obtained using the Mie theory for coated spheres. One-dimensional, steady radiative transport in fluidized dry water system was modeled using the radiative transport equation (RTE) and was solved spectrally with the discrete ordinates method. The effects of silica coating and water droplet size as well as the volume fraction of dry water particles on reducing radiative heat transfer were studied parametrically. The results obtained using the size distributions from experimental measurements at a volume fraction of 10-4 showed that dry water reduced the local radiative heat flux by more than 60% with respect to that by silica particles alone whereas its performance was comparable to that of fine water mists. Moreover, reduction of the diameter of dry water particles from 150 to 50 μm and increasing their volume fraction from 10-4 to 10-3 decreased the radiative heat flux by 45% and 67%, respectively. Dry water is a novel and unique material that does not require high pressure fluid lines for producing fine mists and features a silica shell that can serve to encapsulate water soluble compounds, retard water evaporation from the core as well as increase scattering. With these unique features, dry water finds diverse engineering applications serving as a base for photo-catalytic nanoreactors, gas and chemical storage and delivery systems, as well as alternative mist systems in firefighting.
Kay H. Birdsell; Kathleen M. Bower; Andrew V. Wolfsberg; Wendy E. Soll; Terry A. Cherry; Tade W. Orr
1999-07-01
Numerical simulations are used to predict the migration of radionuclides from the disposal units at Material Disposal Area G through the vadose zone and into the main aquifer in support of a radiological performance assessment and composite analysis for the site. The calculations are performed with the finite element code, FEHM. The transport of nuclides through the vadose zone is computed using a three-dimensional model that describes the complex mesa top geology of the site. The model incorporates the positions and inventories of thirty-four disposal pits and four shaft fields located at Area G as well as those of proposed future pits and shafts. Only three nuclides, C-14, Tc-99, and I-129, proved to be of concern for the groundwater pathway over a 10,000-year period. The spatial and temporal flux of these three nuclides from the vadose zone is applied as a source term for the three-dimensional saturated zone model of the main aquifer that underlies the site. The movement of these nuclides in the aquifer to a downstream location is calculated, and aquifer concentrations are converted to doses. Doses related to aquifer concentrations are six or more orders of magnitude lower than allowable Department of Energy performance objectives for low-level radioactive waste sites. Numerical studies were used to better understand vadose-zone flow through the dry mesa-top environment at Area G. These studies helped define the final model used to model flow and transport through the vadose zone. The study of transient percolation indicates that a steady flow vadose-zone model is adequate for computing contaminant flux to the aquifer. The fracture flow studies and the investigation of the effect of basalt and pumice properties helped us define appropriate hydrologic properties for the modeling. Finally, the evaporation study helped to justify low infiltration rates.
Radiosity Modeling of Radiation Transport in Z Hohlraum Configurations
NASA Astrophysics Data System (ADS)
Vesey, R. A.; Mehlhorn, T. A.
1998-11-01
Vacuum radiation transport in three-dimensional geometry has been modeled, with simple physics assumptions, using the Lightscape^TM commercial radiosity code. This code utilizes progressive, hierarchical radiosity techniques and adaptive mesh refinement to allow greater spatial resolution and unlimited geometric model complexity compared to traditional radiosity methods. Applications to current Z experiments at Sandia will be presented, including (a) uniformity of the radiation flux driving shock physics samples in offset and direct-viewing secondary hohlraums, (b) radiation temperature gradients and capsule illumination in on-axis secondaries driven by static-walled primary hohlraums, and (c) radiation loss from primary hohlraums into the Z anode-cathode gap and feed hardware.
Radiation Transport through cylindrical foams with heated walls
NASA Astrophysics Data System (ADS)
Baker, Kevin; MacLaren, Steve; Kallman, Joshua; Heinz, Ken; Hsing, Warren
2012-10-01
Radiation transport through low density SiO2 foams has been experimentally studied on the Omega laser. In particular these experiments examined the effects on radiation transport when the boundaries of the SiO2 foam are heated such that energy loss to the boundaries is minimized. The initial density of the SiO2 foams was determined by taking an x-ray radiograph of the foams using a monochromatic Henke source at multiple x-ray energies. The radiation drive used to both study the transport in the SiO2 foam as well as to heat the higher density CRF wall was generated in a laser-heated gold hohlraum using ˜7.5 kJ of the laser energy. The time-dependent spatial profile of the heat wave breaking out of the SiO2 foam was detected with an x-ray streak camera coupled with a soft x-ray transmission grating. The Omega DANTE diagnostic measured the radiation drive in the hohlraum and the Omega VISAR diagnostic monitored the spatial temperature gradient in the foam section of the hohlraum.
Radiation inactivation target size of rat adipocyte glucose transporter
Jung, C.Y.; Jacobs, D.B.; Berenski, C.J.; Spangler, R.A.
1987-05-01
In situ assembly states of rat adipocyte glucose transport protein in plasma membrane (PM) and in microsomal pool (MM) were assessed by measuring target size (TS) of D glucose-sensitive, cytochalasin B binding activity. High energy radiation inactivated the binding in both PM and MM by reducing the total capacity of the binding (B/sub T/) without affecting the dissociation constant (K/sub D/). The reduction in B/sub T/ as a function of radiation dose was analyzed based on classical target theory, from which TS was calculated. TS in the PM of insulin-treated adipocytes was 58 KDa. TS in the MM of noninsulin-treated and insulin-treated adipocytes were 112 and 109 KDa, respectively. With MM, however, inactivation data showed anomalously low radiation sensitivities at low radiation doses showing a shoulder in the semilog plots, which may be due to an interaction with a radiation sensitive inhibitor. With these results, they propose the following model: Adipocyte glucose transporter, while exists as a monomer (T) in PM, occurs in MM either as a homodimer (T/sub 2/) or as a heterodimer (TX) with a protein X of a similar size. These dimers (T/sub 2/ or TX) in MM, furthermore, may form a multi-molecular assembly with another, large (300-400 KDa) protein Y, and insulin increases this assembly formation. These putative, transporter-associated proteins X and Y may play an important role in control of transporter distribution between PM and MM, particularly in response to insulin.
Radiation Transport and Shielding for Space Exploration and High Speed Flight Transportation
NASA Technical Reports Server (NTRS)
Maung, Khin Maung; Trapathi, R. K.
1997-01-01
Transportation of ions and neutrons in matter is of direct interest in several technologically important and scientific areas, including space radiation, cosmic ray propagation studies in galactic medium, nuclear power plants and radiological effects that impact industrial and public health. For the proper assessment of radiation exposure, both reliable transport codes and accurate data are needed. Nuclear cross section data is one of the essential inputs into the transport codes. In order to obtain an accurate parametrization of cross section data, theoretical input is indispensable especially for processes where there is little or no experimental data available. In this grant period work has been done on the studies of the use of relativistic equations and their one-body limits. The results will be useful in choosing appropriate effective one-body equation for reaction calculations. Work has also been done to improve upon the data base needed for the transport codes used in the studies of radiation transport and shielding for space exploration and high speed flight transportation. A phenomenological model was developed for the total absorption cross sections valid for any system of charged and/or uncharged collision pairs for the entire energy range. The success of the model is gratifying. It is being used by other federal agencies, national labs and universities. A list of publications based on the work during the grant period is given below and copies are enclosed with this report.
COMPARISON OF MONTE CARLO METHODS FOR NONLINEAR RADIATION TRANSPORT
W. R. MARTIN; F. B. BROWN
2001-03-01
Five Monte Carlo methods for solving the nonlinear thermal radiation transport equations are compared. The methods include the well-known Implicit Monte Carlo method (IMC) developed by Fleck and Cummings, an alternative to IMC developed by Carter and Forest, an ''exact'' method recently developed by Ahrens and Larsen, and two methods recently proposed by Martin and Brown. The five Monte Carlo methods are developed and applied to the radiation transport equation in a medium assuming local thermodynamic equilibrium. Conservation of energy is derived and used to define appropriate material energy update equations for each of the methods. Details of the Monte Carlo implementation are presented, both for the random walk simulation and the material energy update. Simulation results for all five methods are obtained for two infinite medium test problems and a 1-D test problem, all of which have analytical solutions. Conclusions regarding the relative merits of the various schemes are presented.
Shift: A Massively Parallel Monte Carlo Radiation Transport Package
Pandya, Tara M; Johnson, Seth R; Davidson, Gregory G; Evans, Thomas M; Hamilton, Steven P
2015-01-01
This paper discusses the massively-parallel Monte Carlo radiation transport package, Shift, developed at Oak Ridge National Laboratory. It reviews the capabilities, implementation, and parallel performance of this code package. Scaling results demonstrate very good strong and weak scaling behavior of the implemented algorithms. Benchmark results from various reactor problems show that Shift results compare well to other contemporary Monte Carlo codes and experimental results.
Predicting radiative transport properties of plasma sprayed porous ceramics
NASA Astrophysics Data System (ADS)
Wang, B. X.; Zhao, C. Y.
2016-03-01
The typical yttria-stabilized zirconia material for making the thermal barrier coatings (TBCs) is intrinsically semitransparent to thermal radiation, and the unique disordered microstructures in TBCs make them surprisingly highly scattering. To quantitatively understand the influence of disordered microstructures, this paper presents a quantitative prediction on the radiative properties, especially the transport scattering coefficient of plasma sprayed TBC based on microstructure analysis and rigorous electromagnetic theory. The impact of the porosity, shape, size, and orientation of different types of voids on transport scattering coefficient is comprehensively investigated under the discrete dipole approximation. An inverse model integrating these factors together is then proposed to quantitatively connect transport scattering coefficient with microstructural information, which is also validated by available experimental data. Afterwards, an optimization procedure is carried out based on this model to obtain the optimal size and orientation distribution of the microscale voids to achieve the maximal radiation insulation performance at different operating temperatures, providing guidelines for practical coating design and fabrication. This work suggests that the current model is effective and also efficient for connecting scattering properties to microstructures and can be implemented as a quantitative tool for further studies like non-destructive infrared imaging as well as micro/nanoscale thermal design of TBCs.
NERO- a post-maximum supernova radiation transport code
NASA Astrophysics Data System (ADS)
Maurer, I.; Jerkstrand, A.; Mazzali, P. A.; Taubenberger, S.; Hachinger, S.; Kromer, M.; Sim, S.; Hillebrandt, W.
2011-12-01
The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a non-local thermodynamic equilibrium (NLTE) treatment of radiation transport has rarely been studied. In this paper, we present a new SN radiation transport code, NERO, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500 days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare NERO to the radiation transport code of Jerkstrand, Fransson & Kozma and to the nebular code of Mazzali et al. All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, NERO is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' Type Ia supernova (SN Ia) 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post-maximum epochs.
Efficient, Automated Monte Carlo Methods for Radiation Transport
Kong, Rong; Ambrose, Martin; Spanier, Jerome
2012-01-01
Monte Carlo simulations provide an indispensible model for solving radiative transport problems, but their slow convergence inhibits their use as an everyday computational tool. In this paper, we present two new ideas for accelerating the convergence of Monte Carlo algorithms based upon an efficient algorithm that couples simulations of forward and adjoint transport equations. Forward random walks are first processed in stages, each using a fixed sample size, and information from stage k is used to alter the sampling and weighting procedure in stage k + 1. This produces rapid geometric convergence and accounts for dramatic gains in the efficiency of the forward computation. In case still greater accuracy is required in the forward solution, information from an adjoint simulation can be added to extend the geometric learning of the forward solution. The resulting new approach should find widespread use when fast, accurate simulations of the transport equation are needed. PMID:23226872
Whicker, Jeffrey J.; Kirchner, Thomas B.; Breshears, David D.; Field, Jason P.
2012-03-27
The Technical Area 54 (TA-54) Area G disposal facility is used for the disposal of radioactive waste at Los Alamos National Laboratory (LANL). U.S. Department of Energy (DOE) Order 435.1 (DOE, 2001) requires that radioactive waste be managed in a manner that protects public health and safety and the environment. In compliance with that requirement, DOE field sites must prepare and maintain site-specific radiological performance assessments for facilities that receive waste after September 26, 1988. Sites are also required to conduct composite analyses for facilities that receive waste after this date; these analyses account for the cumulative impacts of all waste that has been (and will be) disposed of at the facilities and other sources of radioactive material that may interact with these facilities. LANL issued Revision 4 of the Area G performance assessment and composite analysis in 2008. In support of those analyses, vertical and horizontal sediment flux data were collected at two analog sites, each with different dominant vegetation characteristics, and used to estimate rates of vertical resuspension and wind erosion for Area G. The results of that investigation indicated that there was no net loss of soil at the disposal site due to wind erosion, and suggested minimal impacts of wind on the long-term performance of the facility. However, that study did not evaluate the potential for contaminant transport caused by the horizontal movement of soil particles over long time frames. Since that time, additional field data have been collected to estimate wind threshold velocities for initiating sediment transport due to saltation and rates of sediment transport once those thresholds are reached. Data such as these have been used in the development of the Vegetation Modified Transport (VMTran) model. This model is designed to estimate patterns and long-term rates of contaminant redistribution caused by winds at the site, taking into account the impacts of plant
Stormtime transport of ring current and radiation belt ions
Chen, M.W.; Schulz, M.; Lyons, L.R.; Gorney, D.J.
1993-04-01
This is an investigation of stormtime particle transport that leads to formation of the ring current. The method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants [mu]) in response to model substorm-associated impulses in the convection electric field. The simulation results are compared qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For [mu] approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of the model storm, and the authors find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For [mu] approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For [mu] approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion.
Atmospheric transport, clouds and the Arctic longwave radiation paradox
NASA Astrophysics Data System (ADS)
Sedlar, Joseph
2016-04-01
Clouds interact with radiation, causing variations in the amount of electromagnetic energy reaching the Earth's surface, or escaping the climate system to space. While globally clouds lead to an overall cooling radiative effect at the surface, over the Arctic, where annual cloud fractions are high, the surface cloud radiative effect generally results in a warming. The additional energy input from absorption and re-emission of longwave radiation by the clouds to the surface can have a profound effect on the sea ice state. Anomalous atmospheric transport of heat and moisture into the Arctic, promoting cloud formation and enhancing surface longwave radiation anomalies, has been identified as an important mechanism in preconditioning Arctic sea ice for melt. Longwave radiation is emitted equally in all directions, and changes in the atmospheric infrared emission temperature and emissivity associated with advection of heat and moisture over the Arctic should correspondingly lead to an anomalous signal in longwave radiation at the top of the atmosphere (TOA). To examine the role of atmospheric heat and moisture transport into the Arctic on TOA longwave radiation, infrared satellite sounder observations from AIRS during 2003-2014 are analyzed for summer (JJAS). Thermodynamic metrics are developed to identify months characterized by a high frequency of warm and moist advection into the Arctic, and segregate the 2003-14 time period into climatological and anomalously warm, moist summer months. We find that anomalously warm, moist months result in a significant TOA longwave radiative cooling, which is opposite the forcing signal that the surface experiences during these months. At the timescale of the advective events, 3-10 days, the TOA cooling can be as large as the net surface energy budget during summer. When averaged on the monthly time scale, and over the full Arctic basin (poleward of 75°N), summer months experiencing frequent warm, moist advection events are
Radiative transport in fluorescence-enhanced frequency domain photon migration.
Rasmussen, John C; Joshi, Amit; Pan, Tianshu; Wareing, Todd; McGhee, John; Sevick-Muraca, Eva M
2006-12-01
Small animal optical tomography has significant, but potential application for streamlining drug discovery and pre-clinical investigation of drug candidates. However, accurate modeling of photon propagation in small animal volumes is critical to quantitatively obtain accurate tomographic images. Herein we present solutions from a robust fluorescence-enhanced, frequency domain radiative transport equation (RTE) solver with unique attributes that facilitate its deployment within tomographic algorithms. Specifically, the coupled equations describing time-dependent excitation and emission light transport are solved using discrete ordinates (SN) angular differencing along with linear discontinuous finite-element spatial differencing on unstructured tetrahedral grids. Source iteration in conjunction with diffusion synthetic acceleration is used to iteratively solve the resulting system of equations. This RTE solver can accurately and efficiently predict ballistic as well as diffusion limited transport regimes which could simultaneously exist in small animals. Furthermore, the solver provides accurate solutions on unstructured, tetrahedral grids with relatively large element sizes as compared to commonly employed solvers that use step differencing. The predictions of the solver are validated by a series of frequency-domain, phantom measurements with optical properties ranging from diffusion limited to transport limited propagation. Our results demonstrate that the RTE solution consistently matches measurements made under both diffusion and transport-limited conditions. This work demonstrates the use of an appropriate RTE solver for deployment in small animal optical tomography. PMID:17278821
Applications of the Monte Carlo radiation transport toolkit at LLNL
NASA Astrophysics Data System (ADS)
Sale, Kenneth E.; Bergstrom, Paul M., Jr.; Buck, Richard M.; Cullen, Dermot; Fujino, D.; Hartmann-Siantar, Christine
1999-09-01
Modern Monte Carlo radiation transport codes can be applied to model most applications of radiation, from optical to TeV photons, from thermal neutrons to heavy ions. Simulations can include any desired level of detail in three-dimensional geometries using the right level of detail in the reaction physics. The technology areas to which we have applied these codes include medical applications, defense, safety and security programs, nuclear safeguards and industrial and research system design and control. The main reason such applications are interesting is that by using these tools substantial savings of time and effort (i.e. money) can be realized. In addition it is possible to separate out and investigate computationally effects which can not be isolated and studied in experiments. In model calculations, just as in real life, one must take care in order to get the correct answer to the right question. Advancing computing technology allows extensions of Monte Carlo applications in two directions. First, as computers become more powerful more problems can be accurately modeled. Second, as computing power becomes cheaper Monte Carlo methods become accessible more widely. An overview of the set of Monte Carlo radiation transport tools in use a LLNL will be presented along with a few examples of applications and future directions.
Aksenov, A.F.; Burnazyan, A.I.
1985-03-01
The purpose and application of the provisional standards for radiation safety of crew and passengers in civil aviation are given. The radiation effect of cosmic radiation in flight on civil aviation air transport is described. Standard levels of radiation and conditions of radiation safety are discussed.
NASA Technical Reports Server (NTRS)
Aksenov, A. F.; Burnazyan, A. I.
1985-01-01
The purpose and application of the provisional standards for radiation safety of crew and passengers in civil aviation are given. The radiation effect of cosmic radiation in flight on civil aviation air transport is described. Standard levels of radiation and conditions of radiation safety are discussed.
Los Alamos National Laboratory.
ERIC Educational Resources Information Center
Hammel, Edward F., Jr.
1982-01-01
Current and post World War II scientific research at the Los Alamos National Laboratory (New Mexico) is discussed. The operation of the laboratory, the Los Alamos consultant program, and continuation education, and continuing education activities at the laboratory are also discussed. (JN)
Compendium of Material Composition Data for Radiation Transport Modeling
Williams, Ralph G.; Gesh, Christopher J.; Pagh, Richard T.
2006-10-31
Computational modeling of radiation transport problems including homeland security, radiation shielding and protection, and criticality safety all depend upon material definitions. This document has been created to serve two purposes: 1) to provide a quick reference of material compositions for analysts and 2) a standardized reference to reduce the differences between results from two independent analysts. Analysts are always encountering a variety of materials for which elemental definitions are not readily available or densities are not defined. This document provides a location where unique or hard to define materials will be located to reduce duplication in research for modeling purposes. Additionally, having a common set of material definitions helps to standardize modeling across PNNL and provide two separate researchers the ability to compare different modeling results from a common materials basis.
Radiation safety aspects of commercial high-speed flight transportation
NASA Astrophysics Data System (ADS)
Wilson, John W.; Nealy, John E.; Cucinotta, Francis A.; Shinn, Judy L.; Hajnal, Ferenc; Reginatto, Marcel; Goldhagen, Paul
1995-05-01
High-speed commercial flight transportation is being studied for intercontinental operations in the 21st century, the projected operational characteristics for these aircraft are examined, the radiation environment as it is now known is presented, and the relevant health issues are discussed. Based on a critical examination of the data, a number of specific issues need to be addressed to ensure an adequate knowledge of the ionizing radiation health risks of these aircraft operations. Large uncertainties in our knowledge of the physical fields for high-energy neutrons and multiply-charged ion components need to be reduced. Improved methods for estimating risks in prenatal exposure need to be developed. A firm basis for solar flare monitoring and forecasting needs to be developed with means of exposure abatement.
Radiation safety aspects of commercial high-speed flight transportation
NASA Technical Reports Server (NTRS)
Wilson, John W.; Nealy, John E.; Cucinotta, Francis A.; Shinn, Judy L.; Hajnal, Ferenc; Reginatto, Marcel; Goldhagen, Paul
1995-01-01
High-speed commercial flight transportation is being studied for intercontinental operations in the 21st century, the projected operational characteristics for these aircraft are examined, the radiation environment as it is now known is presented, and the relevant health issues are discussed. Based on a critical examination of the data, a number of specific issues need to be addressed to ensure an adequate knowledge of the ionizing radiation health risks of these aircraft operations. Large uncertainties in our knowledge of the physical fields for high-energy neutrons and multiply-charged ion components need to be reduced. Improved methods for estimating risks in prenatal exposure need to be developed. A firm basis for solar flare monitoring and forecasting needs to be developed with means of exposure abatement.
A FORMALISM FOR COVARIANT POLARIZED RADIATIVE TRANSPORT BY RAY TRACING
Gammie, Charles F.; Leung, Po Kin
2012-06-20
We write down a covariant formalism for polarized radiative transfer appropriate for ray tracing through a turbulent plasma. The polarized radiation field is represented by the polarization tensor (coherency matrix) N{sup {alpha}{beta}} {identical_to} (a{sup {alpha}}{sub k} a*{sup {beta}}{sub k}), where a{sub k} is a Fourier coefficient for the vector potential. Using Maxwell's equations, the Liouville-Vlasov equation, and the WKB approximation, we show that the transport equation in vacuo is k{sup {mu}}{nabla}{sub {mu}} N{sup {alpha}{beta}} = 0. We show that this is equivalent to Broderick and Blandford's formalism based on invariant Stokes parameters and a rotation coefficient, and suggest a modification that may reduce truncation error in some situations. Finally, we write down several alternative approaches to integrating the transfer equation.
Radiation transport modeling using extended quadrature method of moments
NASA Astrophysics Data System (ADS)
Vikas, V.; Hauck, C. D.; Wang, Z. J.; Fox, R. O.
2013-08-01
The radiative transfer equation describes the propagation of radiation through a material medium. While it provides a highly accurate description of the radiation field, the large phase space on which the equation is defined makes it numerically challenging. As a consequence, significant effort has gone into the development of accurate approximation methods. Recently, an extended quadrature method of moments (EQMOM) has been developed to solve univariate population balance equations, which also have a large phase space and thus face similar computational challenges. The distinct advantage of the EQMOM approach over other moment methods is that it generates moment equations that are consistent with a positive phase space density and has a moment inversion algorithm that is fast and efficient. The goal of the current paper is to present the EQMOM method in the context of radiation transport, to discuss advantages and disadvantages, and to demonstrate its performance on a set of standard one-dimensional benchmark problems that encompass optically thin, thick, and transition regimes. Special attention is given in the implementation to the issue of realizability—that is, consistency with a positive phase space density. Numerical results in one dimension are promising and lay the foundation for extending the same framework to multiple dimensions.
High thermal-transport capacity heat pipes for space radiators
NASA Technical Reports Server (NTRS)
Carlson, Albert W.; Gustafson, Eric; Roukis, Susan L.
1987-01-01
This paper presents the results of performance tests of several dual-slot heat pipe test articles. The dual-slot configuration has a very high thermal transport capability and has been identified as a very promising candidate for the radiator system for the NASA Space Station solar dynamic power modules. Two six-foot long aluminum heat pipes were built and tested with ammonia and acetone. A 20-ft long heat pipe was also built and tested with ammonia. The test results have been compared with performance predictions. A thermal transport capacity of 2000 W at an adverse tilt of 1 in. and a 1000 W capacity at an adverse tilt of 2 in. were achieved on the 20-ft long heat pipe. These values are in close agreement with the predicted performance limits.
Radiation Transport in the Atmosphere - Sea Ice - Ocean System.
NASA Astrophysics Data System (ADS)
Jin, Zhonghai
1995-01-01
A comprehensive radiative transfer model for the coupled atmosphere-sea ice-ocean system has been developed. The theoretical work required for constructing such a coupled model is described first. This work extends the discrete ordinate method, which has been proven to be effective in studies of radiative transfer in the atmosphere, to solve the radiative transfer problem pertaining to a system consisting of two strata with different indices of refraction, such as the atmosphere-ocean system and the atmosphere -sea ice-ocean system. The relevant changes (as compared to the standard problem with constant index of refraction throughout the medium) in formulation and solution of the radiative transfer equation, including the proper application of interface and boundary conditions, are presented. This solution is then applied to the atmosphere -sea ice-ocean system to study the solar energy balance in this coupled system. The input parameters required by the model are observable physical properties (e.g., the profiles of temperature and gas concentrations in the atmosphere, and the profiles of temperature, density, and salinity in the ice). The atmosphere, sea ice and ocean are each divided into a sufficient number of layers in the vertical to adequately resolve changes in their optical properties. This model rigorously accounts for the multiple scattering and absorption by atmospheric molecules, clouds, snow and sea water, as well as inclusions in the sea ice, such as brine pockets and air bubbles. The effects of various factors on the solar energy distribution in the entire system have been studied quantitatively. These factors include the ice salinity and density variations, cloud microphysics as well as variations in melt ponds and snow cover on the ice surface. Finally, the coupled radiative transfer model is used to study the impacts of clouds, snow and ice algae on the light transport in sea ice and in the ocean, as well as to simulate spectral irradiance and
New Parallel computing framework for radiation transport codes
Kostin, M.A.; Mokhov, N.V.; Niita, K.; /JAERI, Tokai
2010-09-01
A new parallel computing framework has been developed to use with general-purpose radiation transport codes. The framework was implemented as a C++ module that uses MPI for message passing. The module is significantly independent of radiation transport codes it can be used with, and is connected to the codes by means of a number of interface functions. The framework was integrated with the MARS15 code, and an effort is under way to deploy it in PHITS. Besides the parallel computing functionality, the framework offers a checkpoint facility that allows restarting calculations with a saved checkpoint file. The checkpoint facility can be used in single process calculations as well as in the parallel regime. Several checkpoint files can be merged into one thus combining results of several calculations. The framework also corrects some of the known problems with the scheduling and load balancing found in the original implementations of the parallel computing functionality in MARS15 and PHITS. The framework can be used efficiently on homogeneous systems and networks of workstations, where the interference from the other users is possible.
Hierarchical Adaptive Solution of Radiation Transport Problems on Unstructured Grids
Dr. Cassiano R. E de Oliveira
2008-06-30
Computational radiation transport has steadily gained acceptance in the last decade as a viable modeling tool due to the rapid advancements in computer software and hardware technologies. It can be applied for the analysis of a wide range of problems which arise in nuclear reactor physics, medical physics, atmospheric physics, astrophysics and other areas of engineering physics. However, radiation transport is an extremely chanllenging computational problem since the governing equation is seven-deimensional (3 in space, 2 in direction, 1 in energy, and 1 in time) with a high degree of coupleing betwen these variables. If not careful, this relatively large number of independent variables when discretized can potentially lead to sets of linear equations of intractable size. Though parallel computing has allowed the solution of very large problems, avaliable computational resources will always be finite due to the fact that every more sophisticated multiphysics models are being demanded by industry. There is thus the pressing requirement to optimize the discretizations so as to minimize the effort and maximize the accuracy.
Compendium of Material Composition Data for Radiation Transport Modeling
McConn, Ronald J.; Gesh, Christopher J.; Pagh, Richard T.; Rucker, Robert A.; Williams III, Robert
2011-03-04
Introduction Meaningful simulations of radiation transport applications require realistic definitions of material composition and densities. When seeking that information for applications in fields such as homeland security, radiation shielding and protection, and criticality safety, researchers usually encounter a variety of materials for which elemental compositions are not readily available or densities are not defined. Publication of the Compendium of Material Composition Data for Radiation Transport Modeling, Revision 0, in 2006 was the first step toward mitigating this problem. Revision 0 of this document listed 121 materials, selected mostly from the combined personal libraries of staff at the Pacific Northwest National Laboratory (PNNL), and thus had a scope that was recognized at the time to be limited. Nevertheless, its creation did provide a well-referenced source of some unique or hard-to-define material data in a format that could be used directly in radiation transport calculations being performed at PNNL. Moreover, having a single common set of material definitions also helped to standardize at least one aspect of the various modeling efforts across the laboratory by providing separate researchers the ability to compare different model results using a common basis of materials. The authors of the 2006 compendium understood that, depending on its use and feedback, the compendium would need to be revised to correct errors or inconsistencies in the data for the original 121 materials, as well as to increase (per users suggestions) the number of materials listed. This 2010 revision of the compendium has accomplished both of those objectives. The most obvious change is the increased number of materials from 121 to 372. The not-so-obvious change is the mechanism used to produce the data listed here. The data listed in the 2006 document were compiled, evaluated, entered, and error-checked by a group of individuals essentially by hand, providing no library
Stockpile Stewardship: Los Alamos
McMillan, Charlie; Morgan, Nathanial; Goorley, Tom; Merrill, Frank; Funk, Dave; Korzekwa, Deniece; Laintz, Ken
2012-01-26
"Heritage of Science" is a short video that highlights the Stockpile Stewardship program at Los Alamos National Laboratory. Stockpile Stewardship was conceived in the early 1990s as a national science-based program that could assure the safety, security, and effectiveness of the U.S. nuclear deterrent without the need for full-scale underground nuclear testing. This video was produced by Los Alamos National Laboratory for screening at the Lab's Bradbury Science Museum in Los Alamos, NM and is narrated by science correspondent Miles O'Brien.
Stockpile Stewardship: Los Alamos
McMillan, Charlie; Morgan, Nathanial; Goorley, Tom; Merrill, Frank; Funk, Dave; Korzekwa, Deniece; Laintz, Ken
2014-08-12
"Heritage of Science" is a short video that highlights the Stockpile Stewardship program at Los Alamos National Laboratory. Stockpile Stewardship was conceived in the early 1990s as a national science-based program that could assure the safety, security, and effectiveness of the U.S. nuclear deterrent without the need for full-scale underground nuclear testing. This video was produced by Los Alamos National Laboratory for screening at the Lab's Bradbury Science Museum in Los Alamos, NM and is narrated by science correspondent Miles O'Brien.
Morgan C. White
2000-07-01
The fundamental motivation for the research presented in this dissertation was the need to development a more accurate prediction method for characterization of mixed radiation fields around medical electron accelerators (MEAs). Specifically, a model is developed for simulation of neutron and other particle production from photonuclear reactions and incorporated in the Monte Carlo N-Particle (MCNP) radiation transport code. This extension of the capability within the MCNP code provides for the more accurate assessment of the mixed radiation fields. The Nuclear Theory and Applications group of the Los Alamos National Laboratory has recently provided first-of-a-kind evaluated photonuclear data for a select group of isotopes. These data provide the reaction probabilities as functions of incident photon energy with angular and energy distribution information for all reaction products. The availability of these data is the cornerstone of the new methodology for state-of-the-art mutually coupled photon-neutron transport simulations. The dissertation includes details of the model development and implementation necessary to use the new photonuclear data within MCNP simulations. A new data format has been developed to include tabular photonuclear data. Data are processed from the Evaluated Nuclear Data Format (ENDF) to the new class ''u'' A Compact ENDF (ACE) format using a standalone processing code. MCNP modifications have been completed to enable Monte Carlo sampling of photonuclear reactions. Note that both neutron and gamma production are included in the present model. The new capability has been subjected to extensive verification and validation (V&V) testing. Verification testing has established the expected basic functionality. Two validation projects were undertaken. First, comparisons were made to benchmark data from literature. These calculations demonstrate the accuracy of the new data and transport routines to better than 25 percent. Second, the ability to
P1 Nonconforming Finite Element Method for the Solution of Radiation Transport Problems
NASA Technical Reports Server (NTRS)
Kang, Kab S.
2002-01-01
The simulation of radiation transport in the optically thick flux-limited diffusion regime has been identified as one of the most time-consuming tasks within large simulation codes. Due to multimaterial complex geometry, the radiation transport system must often be solved on unstructured grids. In this paper, we investigate the behavior and the benefits of the unstructured P(sub 1) nonconforming finite element method, which has proven to be flexible and effective on related transport problems, in solving unsteady implicit nonlinear radiation diffusion problems using Newton and Picard linearization methods. Key words. nonconforrning finite elements, radiation transport, inexact Newton linearization, multigrid preconditioning
Validation of a comprehensive space radiation transport code.
Shinn, J L; Cucinotta, F A; Simonsen, L C; Wilson, J W; Badavi, F F; Badhwar, G D; Miller, J; Zeitlin, C; Heilbronn, L; Tripathi, R K; Clowdsley, M S; Heinbockel, J H; Xapsos, M A
1998-12-01
The HZETRN code has been developed over the past decade to evaluate the local radiation fields within sensitive materials on spacecraft in the space environment. Most of the more important nuclear and atomic processes are now modeled and evaluation within a complex spacecraft geometry with differing material components, including transition effects across boundaries of dissimilar materials, are included. The atomic/nuclear database and transport procedures have received limited validation in laboratory testing with high energy ion beams. The codes have been applied in design of the SAGE-III instrument resulting in material changes to control injurious neutron production, in the study of the Space Shuttle single event upsets, and in validation with space measurements (particle telescopes, tissue equivalent proportional counters, CR-39) on Shuttle and Mir. The present paper reviews the code development and presents recent results in laboratory and space flight validation. PMID:11542474
Space radiation transport properties of polyethylene-based composites.
Kaul, R K; Barghouty, A F; Dahche, H M
2004-11-01
Composite materials that can serve as both effective shielding materials against cosmic-ray and energetic solar particles in deep space, as well as structural materials for habitat and spacecraft, remain a critical and mission enabling component in mission planning and exploration. Polyethylene is known to have excellent shielding properties due to its low density, coupled with high hydrogen content. Polyethylene-fiber reinforced composites promise to combine this shielding effectiveness with the required mechanical properties of structural materials. Samples of polyethylene-fiber reinforced epoxy matrix composite 1-5 cm thick were prepared at the NASA Marshall Space Flight Center and tested against a 500 MeV/nucleon Fe beam at the HIMAC facility of NIRS in Chiba, Japan. This paper presents measured and calculated results for the radiation transport properties of these samples. PMID:15644352
Radiation Transport Properties of Polyethylene-Fiber Composites
NASA Technical Reports Server (NTRS)
Kaul, Raj K.; Barghouty, A. F.; Dahche, H. M.
2003-01-01
Composite materials that can both serve as effective shielding materials against cosmic-ray and energetic solar particles in deep space as well as structural materials for habitat and spacecraft remain a critical and mission enabling piece in mission planning and exploration. Polyethylene is known to have excellent shielding properties due to its low density coupled with high hydrogen content. Polyethylene fiber reinforced composites promise to combine this shielding effectiveness with the required mechanical properties of structural materials. Samples of Polyethylene-fiber reinforced epoxy matrix composite 1-5 cm thick were prepared at NASA's Marshall Space Flight Center and tested against 500 MeV/nucleon Fe beam at the HIMAC facility of NIRS in Chiba, Japan. This paper presents measured and calculated results for the radiation transport properties of these samples.
Space radiation transport properties of polyethylene-based composites
NASA Technical Reports Server (NTRS)
Kaul, R. K.; Barghouty, A. F.; Dahche, H. M.
2004-01-01
Composite materials that can serve as both effective shielding materials against cosmic-ray and energetic solar particles in deep space, as well as structural materials for habitat and spacecraft, remain a critical and mission enabling component in mission planning and exploration. Polyethylene is known to have excellent shielding properties due to its low density, coupled with high hydrogen content. Polyethylene-fiber reinforced composites promise to combine this shielding effectiveness with the required mechanical properties of structural materials. Samples of polyethylene-fiber reinforced epoxy matrix composite 1-5 cm thick were prepared at the NASA Marshall Space Flight Center and tested against a 500 MeV/nucleon Fe beam at the HIMAC facility of NIRS in Chiba, Japan. This paper presents measured and calculated results for the radiation transport properties of these samples.
Charge transport properties of CdMnTe radiation detectors
Kim K.; Rafiel, R.; Boardman, M.; Reinhard, I.; Sarbutt, A.; Watt, G.; Watt, C.; Uxa, S.; Prokopovich, D.A.; Belas, E.; Bolotnikov, A.E.; James, R.B.
2012-04-11
Growth, fabrication and characterization of indium-doped cadmium manganese telluride (CdMnTe)radiation detectors have been described. Alpha-particle spectroscopy measurements and time resolved current transient measurements have yielded an average charge collection efficiency approaching 100 %. Spatially resolved charge collection efficiency maps have been produced for a range of detector bias voltages. Inhomogeneities in the charge transport of the CdMnTe crystals have been associated with chains of tellurium inclusions within the detector bulk. Further, it has been shown that the role of tellurium inclusions in degrading chargecollection is reduced with increasing values of bias voltage. The electron transit time was determined from time of flight measurements. From the dependence of drift velocity on applied electric field the electron mobility was found to be n = (718 55) cm2/Vs at room temperature.
Deriving Kubelka-Munk theory from radiative transport.
Sandoval, Christopher; Kim, Arnold D
2014-03-01
We derive Kubelka-Munk (KM) theory systematically from the radiative transport equation (RTE) by analyzing the system of equations resulting from applying the double spherical harmonics method of order one and transforming that system into one governing the positive- and negative-going fluxes. Through this derivation, we establish the theoretical basis of KM theory, identify all parameters, and determine its range of validity. Moreover, we are able to generalize KM theory to take into account general boundary sources and nonhomogeneous terms, for example. The generalized Kubelka-Munk (gKM) equations are also much more accurate at approximating the solution of the RTE. We validate this theory through comparison with numerical solutions of the RTE. PMID:24690662
Acceleration of a Monte Carlo radiation transport code
Hochstedler, R.D.; Smith, L.M.
1996-03-01
Execution time for the Integrated TIGER Series (ITS) Monte Carlo radiation transport code has been reduced by careful re-coding of computationally intensive subroutines. Three test cases for the TIGER (1-D slab geometry), CYLTRAN (2-D cylindrical geometry), and ACCEPT (3-D arbitrary geometry) codes were identified and used to benchmark and profile program execution. Based upon these results, sixteen top time-consuming subroutines were examined and nine of them modified to accelerate computations with equivalent numerical output to the original. The results obtained via this study indicate that speedup factors of 1.90 for the TIGER code, 1.67 for the CYLTRAN code, and 1.11 for the ACCEPT code are achievable. {copyright} {ital 1996 American Institute of Physics.}
The Subelement Sweeping method for radiation transport modeling on polygonal meshes
NASA Astrophysics Data System (ADS)
Gleicher, Frederick
A new memory efficient way of obtaining numerical solutions to the radiation transport equation on random polygon meshes is developed and analyzed. This method is called the Subelement Sweeping method, and is applied to the discrete ordinates form of the 1-D and 2-D mono-energetic transport equations. In the Subelement Sweeping method, the coarse mesh is first subdivided into triangular subelements, and the subelement mesh is then swept to obtain subelement angular flux solutions. As the subelements are swept the scalar flux for the scattering source is interpolated from the coarse mesh. Numerical solutions on the subelements are obtained with the linear discontinuous finite element method, and the resulting angular fluxes are projected back onto the coarse mesh and accumulated into new scalar flux values. The old subelement information is then thrown away allowing the method to be memory efficient. Formulas for the interpolation from the coarse mesh fluxes to the subelement mesh fluxes and for the projection from the subelement fluxes to the coarse mesh fluxes are derived by minimizing the squared error norm between coarse mesh and subelement scalar fluxes. Asymptotic analysis is carried out in 1-D, and the Subelement Sweeping method is shown to yield a valid diffusion discretization on the coarse mesh. Asymptotic analysis is also carried out in 2-D, and the Subelement Sweeping method is shown to have the diffusion limit for orthogonal quadrilateral meshes with some simplifying assumptions. The method was implemented for the 1-D slab geometry, and numerical experiments in 1-D show that the Subelement Sweeping method is at least third order accurate and has the diffusion limit. The method was also implemented in the Capsaicin framework within the Anaheim package for 2-D polygonal meshes at Los Alamos National Laboratory. Numerical experiments in 2-D show that the Subelement Sweeping method is at least second order accurate, and suggest strongly that the diffusion
Status of the MORSE multigroup Monte Carlo radiation transport code
Emmett, M.B.
1993-06-01
There are two versions of the MORSE multigroup Monte Carlo radiation transport computer code system at Oak Ridge National Laboratory. MORSE-CGA is the most well-known and has undergone extensive use for many years. MORSE-SGC was originally developed in about 1980 in order to restructure the cross-section handling and thereby save storage. However, with the advent of new computer systems having much larger storage capacity, that aspect of SGC has become unnecessary. Both versions use data from multigroup cross-section libraries, although in somewhat different formats. MORSE-SGC is the version of MORSE that is part of the SCALE system, but it can also be run stand-alone. Both CGA and SGC use the Multiple Array System (MARS) geometry package. In the last six months the main focus of the work on these two versions has been on making them operational on workstations, in particular, the IBM RISC 6000 family. A new version of SCALE for workstations is being released to the Radiation Shielding Information Center (RSIC). MORSE-CGA, Version 2.0, is also being released to RSIC. Both SGC and CGA have undergone other revisions recently. This paper reports on the current status of the MORSE code system.
An Improved Neutron Transport Algorithm for Space Radiation
NASA Technical Reports Server (NTRS)
Heinbockel, John H.; Clowdsley, Martha S.; Wilson, John W.
2000-01-01
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.
Stormtime transport of ring current and radiation belt ions
NASA Technical Reports Server (NTRS)
Chen, Margaret W.; Schulz, Michael; Lyons, L. R.; Gorney, David J.
1993-01-01
This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative
Guidelines for effective radiation transport for cable SGEMP modeling.
Drumm, Clifton Russell; Fan, Wesley C.; Turner, C. David
2014-07-01
This report describes experiences gained in performing radiation transport computations with the SCEPTRE radiation transport code for System Generated ElectroMagnetic Pulse (SGEMP) applications. SCEPTRE is a complex code requiring a fairly sophisticated user to run the code effectively, so this report provides guidance for analysts interested in performing these types of calculations. One challenge in modeling coupled photon/electron transport for SGEMP is to provide a spatial mesh that is sufficiently resolved to accurately model surface charge emission and charge deposition near material interfaces. The method that has been most commonly used to date to compute cable SGEMP typically requires a sub-micron mesh size near material interfaces, which may be difficult for meshing software to provide for complex geometries. We present here an alternative method for computing cable SGEMP that appears to substantially relax this requirement. The report also investigates the effect of refining the energy mesh and increasing the order of the angular approximation to provide some guidance on determining reasonable parameters for the energy/angular approximation needed for x-ray environments. Conclusions for -ray environments may be quite different and will be treated in a subsequent report. In the course of the energy-mesh refinement studies, a bug in the cross-section generation software was discovered that may cause under prediction of the result by as much as an order of magnitude for the test problem studied here, when the electron energy group widths are much smaller than those for the photons. Results will be presented and compared using cross sections generated before and after the fix. We also describe adjoint modeling, which provides sensitivity of the total charge drive to the source energy and angle of incidence, which is quite useful for comparing the effect of changing the source environment and for determining most stressing angle of incidence and
Modeling of external radiation from the transport of radionuclides across a canyon
Bowen, B.M.; Olsen, W.A.; Van Etten, D.; Chen, I-li
1986-01-01
The Los Alamos Meson Physics Facility (LAMPF) is an 800-million electron volt, l mA intensity linear proton accelerator used for studying subatomic particles at relativistic velocities. Routine operation of the accelerator results in the formation of short-lived air activation products, primarily in the beam stop section of LAMPF. This study presents the results of monitoring and modeling external radiation levels from LAMPF emissions at three locations during 1984. Measured radiation exposures are presented for all three locations during a 49-day period. Hourly radiation levels are calculated for all sites and compared with the prevalent wind patterns during the study period. Predicted daily levels are compared with measured values at all of the sites. Accuracy of the model is compared for day and night conditions. Annual model predictions are also compared with TLD measurements.
Stochastic methods for uncertainty quantification in radiation transport
Fichtl, Erin D; Prinja, Anil K; Warsa, James S
2009-01-01
The use of generalized polynomial chaos (gPC) expansions is investigated for uncertainty quantification in radiation transport. The gPC represents second-order random processes in terms of an expansion of orthogonal polynomials of random variables and is used to represent the uncertain input(s) and unknown(s). We assume a single uncertain input-the total macroscopic cross section-although this does not represent a limitation of the approaches considered here. Two solution methods are examined: The Stochastic Finite Element Method (SFEM) and the Stochastic Collocation Method (SCM). The SFEM entails taking Galerkin projections onto the orthogonal basis, which, for fixed source problems, yields a linear system of fully -coupled equations for the PC coefficients of the unknown. For k-eigenvalue calculations, the SFEM system is non-linear and a Newton-Krylov method is employed to solve it. The SCM utilizes a suitable quadrature rule to compute the moments or PC coefficients of the unknown(s), thus the SCM solution involves a series of independent deterministic transport solutions. The accuracy and efficiency of the two methods are compared and contrasted. The PC coefficients are used to compute the moments and probability density functions of the unknown(s), which are shown to be accurate by comparing with Monte Carlo results. Our work demonstrates that stochastic spectral expansions are a viable alternative to sampling-based uncertainty quantification techniques since both provide a complete characterization of the distribution of the flux and the k-eigenvalue. Furthermore, it is demonstrated that, unlike perturbation methods, SFEM and SCM can handle large parameter uncertainty.
Two dimensional (r-theta) transport model for synchrotron radiation of FRC plasma
NASA Astrophysics Data System (ADS)
Qerushi, Artan; Barnes, Dan; TAE Team
2013-10-01
A two dimensional (r-theta) transport model has been developed for describing the power loss in FRC reactor plasmas and the transport of energy due to synchrotron radiation as well as the transport of energy due to synchrotron radiation. The transport model uses 1d FRC equilibrium profiles and solves the equation of radiative transfer in two dimensions (r-theta) taking into account the absorption and emission of synchrotron radiation. Relativistic expressions are used for both the absorption and the emission coefficients of synchrotron radiation. The reflection of synchrotron radiation from metal walls is taken into account using the approach of Krajcik. The results of the two-dimensional calculations are compared with simpler 1d calculations, which use an approach developed by Dawson and Berk et al., and 0d calculations which use an approach developed by Trubnikov.
NASA Astrophysics Data System (ADS)
Los Alamos National Laboratory in New Mexico is calling for applications for postdoctoral appointments and research fellowships. The positions are available in geoscience as well as other scientific disciplines.The laboratory, which is operated by the University of California for the Department of Energy, awards J. Robert Oppenheimer Research Fellowships to scientists that either have or will soon complete doctoral degrees. The appointments are for two years, are renewable for a third year, and carry a stipend of $51,865 per year. Potential applicants should send a resume or employment application and a statement of research goals to Carol M. Rich, Div. 89, Human Resources Development Division, MS P290, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 by mid-November.
Dosimetry at the Los Alamos Critical Experiments Facility: Past, present, and future
Malenfant, R.E.
1993-10-01
Although the primary reason for the existence of the Los Alamos Critical Experiments Facility is to provide basic data on the physics of systems of fissile material, the physical arrangements and ability to provide sources of radiation have led to applications for all types of radiation dosimetry. In the broad definition of radiation phenomena, the facility has provided sources to evaluate biological effects, radiation shielding and transport, and measurements of basic parameters such as the evaluation of delayed neutron parameters. Within the last 15 years, many of the radiation measurements have been directed to calibration and intercomparison of dosimetry related to nuclear criticality safety. Future plans include (1) the new applications of Godiva IV, a bare-metal pulse assembly, for dosimetry (including an evaluation of neutron and gamma-ray room return); (2) a proposal to relocate the Health Physics Research Reactor from the Oak Ridge National Laboratory to Los Alamos, which will provide the opportunity to continue the application of a primary benchmark source to radiation dosimetry; and (3) a proposal to employ SHEBA, a low-enrichment solution assembly, for accident dosimetry and evaluation.
Coupled radiation transport/thermal analysis of the radiation shield for a space nuclear reactor
NASA Astrophysics Data System (ADS)
Barattino, W. J.
1985-07-01
A coupled radiation transport-heat transfer-stress analysis of the radiation shield for an SP-100 reactor was performed using a numerical code developed at the University of New Mexico and Sandia National Laboratory. For a fast reactor operating at 1.66 MW sub th, the energy deposited and resulting temperature distribution was determined for a shield consisting of tungsten and lithium hydride pressed into a stainless steel honeycomb matrix. While temperature feedback was shown to have a minor effect on energy deposition, the shielding configuration was found to have a major influence in meeting thermal requirements of the lithium hydride. It was shown that a shield optimized for radiation protection will fail because of melting. However, with minor modifications in the shield layering and material selection, the thermal integrity of the shield can be preserved. A shield design of graphite, depleted lithium hydride, tungsten, and natural lithium hydride was shown to satisfy neutron and gamma fluence requirements, maximum temperature limits, and minimize cracking in the LiH portion of the shield.
Deterministic methods in radiation transport. A compilation of papers presented February 4--5, 1992
Rice, A.F.; Roussin, R.W.
1992-06-01
The Seminar on Deterministic Methods in Radiation Transport was held February 4--5, 1992, in Oak Ridge, Tennessee. Eleven presentations were made and the full papers are published in this report, along with three that were submitted but not given orally. These papers represent a good overview of the state of the art in the deterministic solution of radiation transport problems for a variety of applications of current interest to the Radiation Shielding Information Center user community.
Implication of radiative forcing distribution for energy transport
NASA Astrophysics Data System (ADS)
Huang, Yi
2016-04-01
Radiative forcing of a homogeneous greenhouse gas can be very inhomogeneous because the forcing is dependent on other atmospheric and surface variables. In the case of doubling CO2, the mean instantaneous forcing at the top of the atmosphere is found to vary geographically and temporally from positive to negative values, with the range being more than three times the magnitude of the global mean value. The vertical temperature change across the atmospheric column (temperature lapse rate) is found to be the best single predictor for explaining forcing variation. In addition, the masking effects of clouds and water vapor also contribute to forcing inhomogeneity. A regression model that predicts forcing from geophysical variables is constructed. This model can explain more than 90% of the variance of the forcing. Applying this model to analyzing the forcing variation in the CMIP5 models, we find that inter-model discrepancy in CO2 forcing caused by model climatology leads to considerable discrepancy in their projected change in poleward energy transport.
NASA Astrophysics Data System (ADS)
Apruzese, J. P.; Giuliani, J. L.; Hansen, S. B.
2012-09-01
In modeling optically thick, high energy density plasmas (HEDP), radiation transport is of comparable importance to atomic physics, hydrodynamics, and other transport processes. Radiation transport theory is the framework for calculating radiation output based on atomic kinetics, and is thus critical to designing and diagnosing experiments in which radiation production is significant. Starting in the 1960s, the astrophysics community established benchmarks for computational radiation transport, based on a parameterized two-level-atom, and mostly applicable to media with very high optical depths and low collisional quenching of radiative transitions. The purpose of the work reported here is to establish a similar computational radiation transport benchmark that is more relevant to a laboratory HEDP environment. Our model consists of 8 levels of mostly K-shell Al ions, ranging from the ground state of Li-like Al to its bare nucleus. Rates connecting the levels are given by well-known, readily reproducible analytic prescriptions. The results presented consist of level populations as a function of position within the cylindrical medium, and emitted line profiles. For the plasma conditions considered, the magnitudes and spatial variations of the populations are a sensitive indicator of the accuracy of the radiation transport model, and are critical in calculating experimentally relevant quantities such as radiative powers and line ratios.
Roussin, R.W.
1993-03-01
From the very early days in its history Radiation Shielding Information Center (RSIC) has been involved with high energy radiation transport. The National Aeronautics and Space Administration was an early sponsor of RSIC until the completion of the Apollo Moon Exploration Program. In addition, the intranuclear cascade work of Bertini at Oak Ridge National Laboratory provided valuable resources which were made available through RSIC. Over the years, RSIC has had interactions with many of the developers of high energy radiation transport computing technology and data libraries and has been able to collect and disseminate this technology. The current status of this technology will be reviewed and prospects for new advancements will be examined.
Roussin, R.W.
1993-01-01
From the very early days in its history Radiation Shielding Information Center (RSIC) has been involved with high energy radiation transport. The National Aeronautics and Space Administration was an early sponsor of RSIC until the completion of the Apollo Moon Exploration Program. In addition, the intranuclear cascade work of Bertini at Oak Ridge National Laboratory provided valuable resources which were made available through RSIC. Over the years, RSIC has had interactions with many of the developers of high energy radiation transport computing technology and data libraries and has been able to collect and disseminate this technology. The current status of this technology will be reviewed and prospects for new advancements will be examined.
Radiative heat transport instability in a laser produced inhomogeneous plasma
Bychenkov, V. Yu.; Rozmus, W.
2015-08-15
A laser produced high-Z plasma in which an energy balance is achieved due to radiation emission and radiative heat transfer supports ion acoustic instability. A linear dispersion relation is derived, and instability is compared to the radiation cooling instability [R. G. Evans, Plasma Phys. Controlled Fusion 27, 751 (1985)]. Under conditions of indirect drive fusion experiments, the driving term for the instability is the radiative heat flux and, in particular, the density dependence of the radiative heat conductivity. A specific example of thermal Bremsstrahlung radiation source has been considered. This instability may lead to plasma jet formation and anisotropic x-ray generation, thus affecting inertial confinement fusion related experiments.
Serber, R.
1992-01-01
This book contains the 1943 lecture notes of Robert Serber. Serber was a protege of J. Robert Oppenheimer and member of the team that built the first atomic bomb - reveal what the Los Alamos scientists knew, and did not know, about the terrifying weapon they were building.
Investigation of stochastic radiation transport methods in random heterogeneous mixtures
NASA Astrophysics Data System (ADS)
Reinert, Dustin Ray
densely packed, optically thick kernels. A hybrid continuous energy Monte Carlo algorithm that combines homogeneous and explicit geometry models according to the energy dependent optical thickness was also developed. This resonance switch approach exhibited a remarkably high degree of accuracy in performing criticality calculations. The versatility of this hybrid modeling approach makes it an attractive acceleration strategy for a vast array of Monte Carlo radiation transport applications.
Radiation Transport Computation in Stochastic Media: Method and Application
NASA Astrophysics Data System (ADS)
Liang, Chao
Stochastic media, characterized by the stochastic distribution of inclusions in a background medium, are typical radiation transport media encountered in natural or engineering systems. In the community of radiation transport computation, there is always a demand of accurate and efficient methods that can account for the nature of the stochastic distribution. In this dissertation, we focus on methodology development for the radiation transport computation that is applied to neutronic analyses of nuclear reactor designs characterized by the stochastic distribution of particle fuel. Reactor concepts with the employment of a fuel design consisting of a random heterogeneous mixture of fissile material and non-fissile moderator are constantly proposed. Key physical quantities such as core criticality and power distribution, reactivity control design parameters, depletion and fuel burn-up need to be carefully evaluated. In order to meet these practical requirements, we first need to develop accurate and fast computational methods that can effectively account for the stochastic nature of double heterogeneity configuration. A Monte Carlo based method called Chord Length Sampling (CLS) method is considered to be a promising method for analyzing those TRISO-type fueled reactors. Although the CLS method has been proposed for more than two decades and much research has been conducted to enhance its applicability, further efforts are still needed to address some key research gaps that exist for the CLS method. (1) There is a general lack of thorough investigation of the factors that give rise to the inaccuracy of the CLS method found by many researchers. The accuracy of the CLS method depends on the optical and geometric properties of the system. In some specific scenarios, considerable inaccuracies have been reported. However, no research has been providing a clear interpretation of the reasons responsible for the inaccuracy in the reported scenarios. Furthermore, no any
A user's guide to radiation transport in ALEGRA-HEDP : version 4.6.
Mehlhorn, Thomas Alan; Brunner, Thomas A.
2005-01-01
This manual describes the input syntax to the ALEGRA radiation transport package. All input and output variables are defined, as well as all algorithmic controls. This manual describes the radiation input syntax for ALEGRA-HEDP. The ALEGRA manual[2] describes how to run the code and general input syntax. The ALEGRA-HEDP manual[13] describes the input for other physics used in high energy density physics simulations, as well as the opacity models used by this radiation package. An emission model, which is the lowest order radiation transport approximation, is also described in the ALEGRA-HEDP manual. This document is meant to be used with these other manuals.
Radiation inactivation studies of renal brush border water and urea transport
Verkman, A.S.; Dix, J.A.; Seifter, J.L.; Skorecki, K.L.; Jung, C.Y.; Ausiello, D.A.
1985-12-01
Radiation inactivation was used to determine the nature and molecular weight of water and urea transport pathways in brush border membrane vesicles (BBMV) isolated from rabbit renal cortex. BBMV were frozen to -50 degrees C, irradiated with 1.5 MeV electrons, thawed, and assayed for transport or enzyme activity. The freezing process had no effect on enzyme or transport kinetics. BBMV alkaline phosphatase activity gave linear ln(activity) vs. radiation dose plots with a target size of 68 +/- 3 kDa, similar to previously reported values. Water and solute transport were measured using the stopped-flow light-scattering technique. The rates of acetamide and osmotic water transport did not depend on radiation dose (0-7 Mrad), suggesting that transport of these substances does not require a protein carrier. In contrast, urea and thiourea transport gave linear ln(activity) vs. dose curves with a target size of 125-150 kDa; 400 mM urea inhibited thiourea flux by -50% at 0 and 4.7 Mrad, showing that radiation does not affect inhibitor binding to surviving transporters. These studies suggest that BBMV urea transport requires a membrane protein, whereas osmotic water transport does not.
The Los Alamos accelerator code group
Krawczyk, F.L.; Billen, J.H.; Ryne, R.D.; Takeda, Harunori; Young, L.M.
1995-05-01
The Los Alamos Accelerator Code Group (LAACG) is a national resource for members of the accelerator community who use and/or develop software for the design and analysis of particle accelerators, beam transport systems, light sources, storage rings, and components of these systems. Below the authors describe the LAACG`s activities in high performance computing, maintenance and enhancement of POISSON/SUPERFISH and related codes and the dissemination of information on the INTERNET.
Los Alamos low-level waste performance assessment status
Wenzel, W.J.; Purtymun, W.D.; Dewart, J.M.; Rodgers, J.E.
1986-06-01
This report reviews the documented Los Alamos studies done to assess the containment of buried hazardous wastes. Five sections logically present the environmental studies, operational source terms, transport pathways, environmental dosimetry, and computer model development and use. This review gives a general picture of the Los Alamos solid waste disposal and liquid effluent sites and is intended for technical readers with waste management and environmental science backgrounds but without a detailed familiarization with Los Alamos. The review begins with a wide perspective on environmental studies at Los Alamos. Hydrology, geology, and meteorology are described for the site and region. The ongoing Laboratory-wide environmental surveillance and waste management environmental studies are presented. The next section describes the waste disposal sites and summarizes the current source terms for these sites. Hazardous chemical wastes and liquid effluents are also addressed by describing the sites and canyons that are impacted. The review then focuses on the transport pathways addressed mainly in reports by Healy and Formerly Utilized Sites Remedial Action Program. Once the source terms and potential transport pathways are described, the dose assessment methods are addressed. Three major studies, the waste alternatives, Hansen and Rogers, and the Pantex Environmental Impact Statement, contributed to the current Los Alamos dose assessment methodology. Finally, the current Los Alamos groundwater, surface water, and environmental assessment models for these mesa top and canyon sites are described.
Roussin, R.W.
1994-10-01
From the very early days in its history RSIC has been involved with high energy radiation transport. The National Aeronautics and Space Administration was an early sponsor of RSIC until the completion of the Apollo Moon Exploration Program. In addition, the intranuclear cascade work of Bertini at Oak Ridge National Laboratory provided valuable resources which were made available through RSIC. Over the years, RSIC has had interactions with many of the developers of high energy radiation transport computing technology and data libraries and has been able to collect and disseminate this technology. The current status of this technology will be reviewed and prospects for new advancements will be examined.
Sandia Computational Engine for Particle Transport for Radiation Effects.
Energy Science and Technology Software Center (ESTSC)
2014-09-01
Version 00 The SCEPTRE code solves the linear Boltzmann transport equation for one-, two- and three-dimensional geometries. SCEPTRE is capable of handling any particle type for which multigroup-Legendre cross sections are available. However, the code is designed primarily to model the transport of photons, electrons, and positrons through matter. For efficiency and flexibility, SCEPTRE contains capability for both the first- and second-order forms of the Boltzmann transport equation.
Coupled Deterministic/Monte Carlo Simulation of Radiation Transport and Detector Response
Gesh, Christopher J.; Meriwether, George H.; Pagh, Richard T.; Smith, Leon E.
2005-09-01
The analysis of radiation sensor systems used to detect and identify nuclear and radiological weapons materials requires detailed radiation transport calculations. Two basic steps are required to solve radiation detection scenario analysis (RDSA) problems. First, the radiation field produced by the source must be calculated. Second, the response that the radiation field produces in a detector must be determined. RDSA problems are characterized by complex geometries, the presence of shielding materials, and large amounts of scattering (or absorption/re-emission). In this paper, we will discuss the use of the Attila code [2] for RDSA.
Radiation transport phenomena and modeling. Part A: Codes; Part B: Applications with examples
Lorence, L.J. Jr.; Beutler, D.E.
1997-09-01
This report contains the notes from the second session of the 1997 IEEE Nuclear and Space Radiation Effects Conference Short Course on Applying Computer Simulation Tools to Radiation Effects Problems. Part A discusses the physical phenomena modeled in radiation transport codes and various types of algorithmic implementations. Part B gives examples of how these codes can be used to design experiments whose results can be easily analyzed and describes how to calculate quantities of interest for electronic devices.
Ku, L.P.; Kolibal, J.G.; Liew, S.L.
1985-09-01
The computational models of the TFTR constructed for the radiation transport analysis for the Q approx. 1 demonstration are summarized and reviewed. These models can be characterized by the dimensionality required to describe the geometry, and by the numerical methods of solving the transport equation. Results obtained with these models in the test cell are compared and discussed.
Multigroup Three-Dimensional Direct Integration Method Radiation Transport Analysis Code System.
Energy Science and Technology Software Center (ESTSC)
1987-09-18
Version 00 TRISTAN solves the three-dimensional, fixed-source, Boltzmann transport equation for neutrons or gamma rays in rectangular geometry. The code can solve an adjoint problem as well as a usual transport problem. TRISTAN is a suitable tool to analyze radiation shielding problems such as streaming and deep penetration problems.
A Photon Free Method to Solve Radiation Transport Equations
Chang, B
2006-09-05
The multi-group discrete-ordinate equations of radiation transfer is solved for the first time by Newton's method. It is a photon free method because the photon variables are eliminated from the radiation equations to yield a N{sub group}XN{sub direction} smaller but equivalent system of equations. The smaller set of equations can be solved more efficiently than the original set of equations. Newton's method is more stable than the Semi-implicit Linear method currently used by conventional radiation codes.
Novel Parallel Numerical Methods for Radiation& Neutron Transport
Brown, P N
2001-03-06
In many of the multiphysics simulations performed at LLNL, transport calculations can take up 30 to 50% of the total run time. If Monte Carlo methods are used, the percentage can be as high as 80%. Thus, a significant core competence in the formulation, software implementation, and solution of the numerical problems arising in transport modeling is essential to Laboratory and DOE research. In this project, we worked on developing scalable solution methods for the equations that model the transport of photons and neutrons through materials. Our goal was to reduce the transport solve time in these simulations by means of more advanced numerical methods and their parallel implementations. These methods must be scalable, that is, the time to solution must remain constant as the problem size grows and additional computer resources are used. For iterative methods, scalability requires that (1) the number of iterations to reach convergence is independent of problem size, and (2) that the computational cost grows linearly with problem size. We focused on deterministic approaches to transport, building on our earlier work in which we performed a new, detailed analysis of some existing transport methods and developed new approaches. The Boltzmann equation (the underlying equation to be solved) and various solution methods have been developed over many years. Consequently, many laboratory codes are based on these methods, which are in some cases decades old. For the transport of x-rays through partially ionized plasmas in local thermodynamic equilibrium, the transport equation is coupled to nonlinear diffusion equations for the electron and ion temperatures via the highly nonlinear Planck function. We investigated the suitability of traditional-solution approaches to transport on terascale architectures and also designed new scalable algorithms; in some cases, we investigated hybrid approaches that combined both.
Brownell, J.H.; Bowers, R.L.
1997-04-01
The Los Alamos foil implosion program has the goal of producing an intense, high-energy density x-ray source by converting the energy of a magnetically imploded plasma into radiation and material energy. One of the methods for converting the plasma energy into thermal energy and radiation and utilizing it for experiments is called the flying radiation case (FRC). In this paper the authors shall model the FRC and provide a physical description of the processes involved. An analytic model of a planar FRC in the hydrodynamic approximation is used to describe the assembly and shock heating of a central cushion by a conducting liner driver. The results are also used to benchmark a hydrodynamics code for modeling an FRC. They then use a radiation-hydrodynamics computational model to explore the effects of radiation production and transport when a gold plasma assembles on a CH cushion. Results are presented for the structure and evolution of the radiation hohlraum.
Review of liquid metal heat pipe work at Los Alamos
Reid, R.S.; Merrigan, M.A.; Sena, J.T.
1990-01-01
A survey of space-power related liquid metal heat pipe work at Los Alamos National Laboratory is presented. Heat pipe development at Los Alamos has been on-going since 1963. Heat pipes were initially developed for thermionic nuclear-electrical power production in space. Since then Los Alamos has developed liquid metal heat pipes for numerous applications related to high temperature systems in both the space and terrestrial environments. Some of these applications include thermionic electrical generators, thermoelectric energy conversion (both in-core and direct radiation), thermal energy storage, hypersonic vehicle leading edge cooling, and heat pipe vapor laser cells. Some of the work performed at Los Alamos has been documented in internal reports that are often little-known. A representative description and summary of progress in space-related liquid metal heat pipe technology is provided followed by a reference section citing sources where these works may be found. 53 refs.
Recent Developments in Three Dimensional Radiation Transport Using the Green's Function Technique
NASA Technical Reports Server (NTRS)
Rockell, Candice; Tweed, John; Blattnig, Steve R.; Mertens, Christopher J.
2010-01-01
In the future, astronauts will be sent into space for longer durations of time compared to previous missions. The increased risk of exposure to dangerous radiation, such as Galactic Cosmic Rays and Solar Particle Events, is of great concern. Consequently, steps must be taken to ensure astronaut safety by providing adequate shielding. In order to better determine and verify shielding requirements, an accurate and efficient radiation transport code based on a fully three dimensional radiation transport model using the Green's function technique is being developed
Radiation transport and density effects in non-equilibrium plasmas
NASA Astrophysics Data System (ADS)
Fisher, Vladimir I.; Fisher, Dimitri V.; Maron, Yitzhak
2007-05-01
We describe a model for self-consistent computations of ionic level populations and the radiation field in transient non-equilibrium plasmas. In this model, the plasma density effects are described using the effective-statistical-weights (ESW) formalism based on the statistics of the microscopic environment of individual ions. In comparison to earlier work, the ESW formalism is expanded to a self-consistent treatment of the radiative transfer. For non-Maxwellian plasmas, the atomic-kinetics and radiative transfer computations may be performed for an arbitrary distribution of the free electrons. A plasma is presented by a finite number of cells, each with uniform thermodynamic parameters. The radiation field in each cell is computed by accounting for the radiation of entire plasma and of external sources. To demonstrate the predictions of the ESW approach and their difference from those of the traditional approach we apply the model to high-density plasmas. Based on hydrodynamic simulations of a laser-matter interaction, we use the model to analyze spectral line shapes, where the effects caused by the spatial dependence of the plasma flow velocity are demonstrated. In single-cell simulations, for acceleration of the computations, the model utilizes recently derived formula for the cell volume-average and direction-average specific intensity of radiation.
Transport analysis of high radiation and high density plasmas in the ASDEX Upgrade tokamak
NASA Astrophysics Data System (ADS)
Casali, L.; Bernert, M.; Dux, R.; Fischer, R.; Kallenbach, A.; Kurzan, B.; Lang, P.; Mlynek, A.; McDermott, R. M.; Ryter, F.; Sertoli, M.; Tardini, G.; Zohm, H.
2014-12-01
Future fusion reactors, foreseen in the "European road map" such as DEMO, will operate under more demanding conditions compared to present devices. They will require high divertor and core radiation by impurity seeding to reduce heat loads on divertor target plates. In addition, DEMO will have to work at high core densities to reach adequate fusion performance. The performance of fusion reactors depends on three essential parameters: temperature, density and energy confinement time. The latter characterizes the loss rate due to both radiation and transport processes. The DEMO foreseen scenarios described above were not investigated so far, but are now addressed at the ASDEX Upgrade tokamak. In this work we present the transport analysis of such scenarios.
On the Development of a Deterministic Three-Dimensional Radiation Transport Code
NASA Technical Reports Server (NTRS)
Rockell, Candice; Tweed, John
2011-01-01
Since astronauts on future deep space missions will be exposed to dangerous radiations, there is a need to accurately model the transport of radiation through shielding materials and to estimate the received radiation dose. In response to this need a three dimensional deterministic code for space radiation transport is now under development. The new code GRNTRN is based on a Green's function solution of the Boltzmann transport equation that is constructed in the form of a Neumann series. Analytical approximations will be obtained for the first three terms of the Neumann series and the remainder will be estimated by a non-perturbative technique . This work discusses progress made to date and exhibits some computations based on the first two Neumann series terms.
Ukhorskiy, A Y; Sitnov, M I; Millan, R M; Kress, B T; Smith, D C
2014-01-01
[1]Relativistic electron intensities in Earth's outer radiation belt can vary by multiple orders of magnitude on the time scales ranging from minutes to days. One fundamental process contributing to dynamic variability of radiation belt intensities is the radial transport of relativistic electrons across their drift shells. In this paper we analyze the properties of three-dimensional radial transport in a global magnetic field model driven by variations in the solar wind dynamic pressure. We use a test particle approach which captures anomalous effects such as drift orbit bifurcations. We show that the bifurcations lead to an order of magnitude increase in radial transport rates and enhance the energization at large equatorial pitch angles. Even at quiet time fluctuations in dynamic pressure, radial transport at large pitch angles exhibits strong deviations from the diffusion approximation. The radial transport rates are much lower at small pitch angle values which results in a better agreement with the diffusion approximation. PMID:26167431
NASA Astrophysics Data System (ADS)
Takabe, Hideaki
A brief review is given of the physics of radiation transport, a topic that is important in the study of astrophysics, laser-plasmas, divertor-plasmas, etc. In general, we must solve non-local thermodynamic equilibrium processes using an appropriate atomic model. The resultant data related to the spectral emissivity and opacity of partially ionized plasmas are then used to solve the radiation transfer equation. In this note, I briefly overview a variety of ways to carry out such a calculation. In addition, similarities and differences in the physical process between laser-plasmas and divertor-plasmas are briefly described.
The Roles of Transport and Wave-Particle Interactions on Radiation Belt Dynamics
NASA Technical Reports Server (NTRS)
Fok, Mei-Ching; Glocer, Alex; Zheng, Qiuhua
2011-01-01
Particle fluxes in the radiation belts can vary dramatically during geomagnetic active periods. Transport and wave-particle interactions are believed to be the two main types of mechanisms that control the radiation belt dynamics. Major transport processes include substorm dipolarization and injection, radial diffusion, convection, adiabatic acceleration and deceleration, and magnetopause shadowing. Energetic electrons and ions are also subjected to pitch-angle and energy diffusion when interact with plasma waves in the radiation belts. Important wave modes include whistler mode chorus waves, plasmaspheric hiss, electromagnetic ion cyclotron waves, and magnetosonic waves. We investigate the relative roles of transport and wave associated processes in radiation belt variations. Energetic electron fluxes during several storms are simulated using our Radiation Belt Environment (RBE) model. The model includes important transport and wave processes such as substorm dipolarization in global MHD fields, chorus waves, and plasmaspheric hiss. We discuss the effects of these competing processes at different phases of the storms and validate the results by comparison with satellite and ground-based observations. Keywords: Radiation Belts, Space Weather, Wave-Particle Interaction, Storm and Substorm
Energy distributions and radiation transport in uranium plasmas
NASA Technical Reports Server (NTRS)
Miley, G. H.; Bathke, C.; Maceda, E.; Choi, C.
1976-01-01
An approximate analytic model, based on continuous electron slowing, has been used for survey calculations. Where more accuracy is required, a Monte Carlo technique is used which combines an analytic representation of Coulombic collisions with a random walk treatment of inelastic collisions. The calculated electron distributions have been incorporated into another code that evaluates both the excited atomic state densities within the plasma and the radiative flux emitted from the plasma.
Radiation Effects on Transport and Bubble Formation in Silicate Glasses
Trifunac, A.D.; Shkrob, I.A.; Werst, D.W.
2001-12-31
Using advanced magnetic resonance spectroscopies and small-cluster modeling, atomic structure of radiation-induced point defects in alkali borate, silicate, and borosilicate glasses is fully characterized. It is shown that in boron-containing glasses, most of these point defects are electrons/holes trapped by cation/anion vacancies, such as O1 - - O3 + valence-alternation pairs. In microscopically phase-separated borosilicate glasses, radiation-induced defects are found to cluster at the interface between the borate and silicate phases. Reaction and diffusion dynamics of defect-annealing interstitial hydrogen atoms in boron and silica oxide glasses are studied. The yield of radiolytic O2 is estimated. This oxygen is shown to be the final product of triplet exciton decay. Plausible mechanisms for the oxygen bubble formation are put forward. Two practical conclusions relevant for the EMSP mission are made: First, the yield of radiolytic oxygen is shown to be too low to interfere with the storage of vitrified radioactive waste in the first 10 Kyr. Second, microscopic phase separation is demonstrated to increase both the chemical and radiation stability of borosilicate glass.
Analytical Radiation Transport Benchmarks for The Next Century
B.D. Ganapol
2005-01-19
Verification of large-scale computational algorithms used in nuclear engineering and radiological applications is an essential element of reliable code performance. For this reason, the development of a suite of multidimensional semi-analytical benchmarks has been undertaken to provide independent verification of proper operation of codes dealing with the transport of neutral particles. The benchmarks considered cover several one-dimensional, multidimensional, monoenergetic and multigroup, fixed source and critical transport scenarios. The first approach, called the Green's Function. In slab geometry, the Green's function is incorporated into a set of integral equations for the boundary fluxes. Through a numerical Fourier transform inversion and subsequent matrix inversion for the boundary fluxes, a semi-analytical benchmark emerges. Multidimensional solutions in a variety of infinite media are also based on the slab Green's function. In a second approach, a new converged SN method is developed. In this method, the SN solution is ''minded'' to bring out hidden high quality solutions. For this case multigroup fixed source and criticality transport problems are considered. Remarkably accurate solutions can be obtained with this new method called the Multigroup Converged SN (MGCSN) method as will be demonstrated.
MESTRN: A Deterministic Meson-Muon Transport Code for Space Radiation
NASA Technical Reports Server (NTRS)
Blattnig, Steve R.; Norbury, John W.; Norman, Ryan B.; Wilson, John W.; Singleterry, Robert C., Jr.; Tripathi, Ram K.
2004-01-01
A safe and efficient exploration of space requires an understanding of space radiations, so that human life and sensitive equipment can be protected. On the way to these sensitive sites, the radiation fields are modified in both quality and quantity. Many of these modifications are thought to be due to the production of pions and muons in the interactions between the radiation and intervening matter. A method used to predict the effects of the presence of these particles on the transport of radiation through materials is developed. This method was then used to develop software, which was used to calculate the fluxes of pions and muons after the transport of a cosmic ray spectrum through aluminum and water. Software descriptions are given in the appendices.
Radiation transport and energetics of laser-driven half-hohlraums at the National Ignition Facility
Moore, A. S. Graham, P.; Comley, A. J.; Foster, J.; Cooper, A. B. R.; Schneider, M. B.; MacLaren, S.; Lu, K.; Seugling, R.; Satcher, J.; Klingmann, J.; Marrs, R.; May, M.; Widmann, K.; Glendinning, G.; Castor, J.; Sain, J.; Baker, K.; Hsing, W. W.; Young, B.; and others
2014-06-15
Experiments that characterize and develop a high energy-density half-hohlraum platform for use in benchmarking radiation hydrodynamics models have been conducted at the National Ignition Facility (NIF). Results from the experiments are used to quantitatively compare with simulations of the radiation transported through an evolving plasma density structure, colloquially known as an N-wave. A half-hohlraum is heated by 80 NIF beams to a temperature of 240 eV. This creates a subsonic diffusive Marshak wave, which propagates into a high atomic number Ta{sub 2}O{sub 5} aerogel. The subsequent radiation transport through the aerogel and through slots cut into the aerogel layer is investigated. We describe a set of experiments that test the hohlraum performance and report on a range of x-ray measurements that absolutely quantify the energetics and radiation partition inside the target.
Los Alamos National Laboratory
Dogliani, Harold O
2011-01-19
The purpose of the briefing is to describe general laboratory technical capabilities to be used for various groups such as military cadets or university faculty/students and post docs to recruit into a variety of Los Alamos programs. Discussed are: (1) development and application of high leverage science to enable effeictive, predictable and reliability outcomes; (2) deter, detect, characterize, reverse and prevent the proliferation of weapons of mass destruction and their use by adversaries and terrorists; (3) modeling and simulation to define complex processes, predict outcomes, and develop effective prevention, response, and remediation strategies; (4) energetic materials and hydrodynamic testing to develop materials for precise delivery of focused energy; (5) materials cience focused on fundamental understanding of materials behaviors, their quantum-molecular properties, and their dynamic responses, and (6) bio-science to rapidly detect and characterize pathogens, to develop vaccines and prophylactic remedies, and to develop attribution forensics.
Kirk, B.L. )
1990-01-01
In nuclear applications, the conversion of mainframe software to the personal computer (PC) environment has seen an accelerated pace. Credit has to be extended to the software companies that have made the scientific language FORTRAN available on PCs. Not to be neglected are the scientists who dedicate their time in the conversion of codes and are challenged by the limited PC memory and disk space. The Radiation shielding Information Center (RSIC) at Oak Ridge National Laboratory has encouraged these developments, and the shielding community has cooperated by making these new tools available via RSIC. The PC codes in the shielding and radiation transport area are divided into five categories (these categories are not mutually exclusive): (1) gamma-ray scattering; (2) neutron and gamma-ray transport (also coupled); (3) environmental dose; (4) medical applications; and (5) reactor physics. Each category is discussed.
HETC radiation transport code development for cosmic ray shielding applications in space.
Townsend, L W; Miller, T M; Gabriel, Tony A
2005-01-01
In order to facilitate three-dimensional analyses of space radiation shielding scenarios for future space missions, the Monte Carlo radiation transport code HETC is being extended to include transport of energetic heavy ions, such as are found in the galactic cosmic ray spectrum in space. Recently, an event generator capable of providing nuclear interaction data for use in HETC was developed and incorporated into the code. The event generator predicts the interaction product yields and production angles and energies using nuclear models and Monte Carlo techniques. Testing and validation of the extended transport code has begun. In this work, the current status of code modifications, which enable energetic heavy ions and their nuclear reaction products to be transported through thick shielding, are described. Also, initial results of code testing against available laboratory beam data for energetic heavy ions interacting in thick targets are presented. PMID:16604614
Gleckler, P.J.; Randall, D.A.; Boer, G.
1994-03-01
This paper reports on energy fluxes across the surface of the ocean as simulated by fifteen atmospheric general circulation models in which ocean surface temperatures and sea-ice boundaries are prescribed. The oceanic meridional energy transport that would be required to balance these surface fluxes is computed, and is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean energy transport can be affected by the errors in simulated cloud-radiation interactions.
NASA Technical Reports Server (NTRS)
Hartung, Lin C.; Hassan, H. A.
1992-01-01
A moment method for computing 3D radiative transport in axisymmetric thermochemical nonequilibrium flows is developed. The method uses the P-1 approximation to reduce the governing system of integro-differential equations to a coupled set of partial differential equations. The numerical solution of these equations for realistic variations of the radiation properties is discussed. Representative results from the method are shown and compared to tangent slab calculations. The agreement between the transport methods is found to be about 10 percent in the stagnation region, with the difference increasing along the flank of the vehicle.
Basic knowledge on radiative and transport properties to begin in thermal plasmas modelling
Cressault, Y.
2015-05-15
This paper has for objectives to present the radiative and the transport properties for people beginning in thermal plasmas. The first section will briefly recall the equations defined in numerical models applied to thermal plasmas; the second section will particularly deal with the estimation of radiative losses; the third part will quickly present the thermodynamics properties; and the last part will concern the transport coefficients (thermal conductivity, viscosity and electrical conductivity of the gas or mixtures of gases). We shall conclude the paper with a discussion about the validity of these results the lack of data for some specific applications, and some perspectives concerning these properties for non-equilibrium thermal plasmas.
Towards a 3D Space Radiation Transport Code
NASA Technical Reports Server (NTRS)
Wilson, J. W.; Tripathl, R. K.; Cicomptta, F. A.; Heinbockel, J. H.; Tweed, J.
2002-01-01
High-speed computational procedures for space radiation shielding have relied on asymptotic expansions in terms of the off-axis scatter and replacement of the general geometry problem by a collection of flat plates. This type of solution was derived for application to human rated systems in which the radius of the shielded volume is large compared to the off-axis diffusion limiting leakage at lateral boundaries. Over the decades these computational codes are relatively complete and lateral diffusion effects are now being added. The analysis for developing a practical full 3D space shielding code is presented.
The implication of radiative forcing and feedback for meridional energy transport
NASA Astrophysics Data System (ADS)
Huang, Yi; Zhang, Minghong
2014-03-01
The distributions of radiative forcing and feedback in the Coupled Model Intercomparison Project phase 5 abrupt4xCO2 and Historical experiments are diagnosed, with a focus on their effects on the zonal mean structure of the top-of-the-atmosphere radiation anomalies and implications for the meridional energy transport. It is found that because the greenhouse gas longwave forcing peaks in the low latitudes, it reinforces the equator-to-pole net radiation gradient and accounts for the increase in the poleward energy transport in both hemispheres under global warming. The shortwave forcing by aerosol, ozone, etc. peaks in the Northern Hemisphere and instead implies an interhemispheric energy transport. Although the water vapor feedback also reinforces the equator-to-pole gradient of the net radiation, the temperature and albedo feedback act against it. The feedback tend to offset the zonal mean radiation anomaly caused by the forcing, although the overall feedback effect on the energy transport is rather uncertain, mainly due to the uncertainty in the cloud feedback.
A Radiation Transport Coupled Particle-In-Cell Model for Hg-Ar Discharges
NASA Astrophysics Data System (ADS)
Lee, Hae June; Verboncoeur, J. P.; Smith, H. B.; Parker, G. J.; Birdsall, C. K.
2000-10-01
We simulate a radial slice of the fluorescent lamp discharge in the positive column with a radiation transport coupled particle-in-cell (RT-PIC) code. In this model, the radiative and meta stable excited states of Hg-Ar mixture and their collisions as well as radiation transport are simulated by the fluid equations. The motions of electrons and ions and collisions with neutral or excited states are simulated by the conventional particle-in-cell method. We consider radiation transport of excited states using the Holstein equation[1] including the time varying nonuniform background gas density. The background gas density is calculated from the temperature profile by solving the heat transfer equation. The motion of charged particles are simulated by using the 1-D cylindrical particle-in-cell code, XPDC1[2]. Separate time scales are used for the charged particles, the excited states, and the neutral gas, respectively, and parallel processing can be used for the expensive calculation including radiation transport. This work was supported in part by General Electric Company contract GE-20000181. [1] T. Holstein, Phys. Rev. 72, 1213 (1947). [2] J. P. Verboncoeur, M. V. Alves, V. Vahedi, and C. K. Birdsall, Journal of Computational Physics 104(2), 321, (1993).
Experiences with radiation portal detectors for international rail transport
NASA Astrophysics Data System (ADS)
Stromswold, D. C.; McCormick, K.; Todd, L.; Ashbaker, E. D.; Evans, J. C.
2006-08-01
Radiation detectors monitored trains at two international borders to evaluate the performance of NaI(Tl) and plastic (polyvinyltoluene: PVT) gamma-ray detectors to characterize rail cargo. The detectors included a prototype NaI(Tl) radiation-portal-monitor panel having four large detectors (10-cm × 10-cm × 41-cm) and a PVT panel with a 41 cm × 173 cm × 3.8-cm detector. Spectral data from the NaI(Tl) and PVT detectors were recorded. Of particular emphasis was the identification of naturally occurring radioactive material (NORM) and the resultant frequency of nuisance alarms. For rail monitoring, the difficulty in stopping trains to perform secondary inspection on alarming cars creates a need for reliable identification of NORM during initial screening. Approximately 30 trains were monitored, and the commodities in individual railcars were ascertained from manifest information. At one test site, the trains carried inter-modal containers that had been unloaded from ships, and at the other site, the trains contained bulk cargo in tanker cars and hopper cars or individual items in boxcars or flatbeds. NORM encountered included potash, liquefied petroleum gas, fireworks, televisions, and clay-based products (e.g., pottery). Analysis of the spectral data included the use of the template-fitting portion of the program GADRAS developed at Sandia National Laboratories. For most of the NORM, the NaI(Tl) data produced a correct identification of the radionuclides present in the railcars. The same analysis was also used for PVT data in which the spectral information (no peaks but only gradual spectral changes including Compton edges) was limited. However, the PVT analysis provided correct identification of 40K and 226Ra in many cases.
Experiences with radiation portal detectors for international rail transport
Stromswold, David C.; McCormick, Kathleen R.; Todd, Lindsay C.; Ashbaker, Eric D.; Evans, J. C.
2006-08-30
Radiation detectors monitored trains at two international borders to evaluate the performance of NaI(Tl) and plastic (polyvinyltoluene: PVT) gamma-ray detectors to characterize rail cargo. The detectors included a prototype NaI(Tl) radiation-portal-monitor panel having four large detectors (10-cm × 10-cm × 41-cm) and a PVT panel with a 41 cm × 173 cm × 3.8-cm detector. Spectral data from the NaI(Tl) and PVT detectors were recorded. Of particular emphasis was the identification of naturally occurring radioactive material (NORM) and the resultant frequency of nuisance alarms. For rail monitoring, the difficulty in stopping trains to perform secondary inspection on alarming cars creates a need for reliable identification of NORM during initial screening. Approximately 30 trains were monitored, and the commodities in individual railcars were ascertained from manifest information. At one test site the trains carried inter-modal containers that had been unloaded from ships, and at the other site the trains contained bulk cargo or individual items in boxcars or flatbeds. NORM encountered included potash, liquefied petroleum gas, fireworks, televisions, and clay-based products (e.g., pottery). Analysis of the spectral data included the use of the template-fitting program GADRAS/FitToDB from Sandia National Laboratories. For much of the NORM the NaI(Tl) data produced a correct identification of the radionuclides present in the railcars. The same analysis was also used for PVT data in which the spectral information (no peaks but only gradual spectral changes including Compton edges) was limited. However, the PVT analysis provided correct identification of 40K and 226Ra in many cases.
A General Hybrid Radiation Transport Scheme for Star Formation Simulations on an Adaptive Grid
NASA Astrophysics Data System (ADS)
Klassen, Mikhail; Kuiper, Rolf; Pudritz, Ralph E.; Peters, Thomas; Banerjee, Robi; Buntemeyer, Lars
2014-12-01
Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodyanmics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.
Jacobs, D.B.; Berenski, C.J.; Spangler, R.A.; Jung, C.Y.
1987-06-15
The in situ assembly states of the glucose transport carrier protein in the plasma membrane and in the intracellular (microsomal) storage pool of rat adipocytes were assessed by studying radiation-induced inactivation of the D-glucose-sensitive cytochalasin B binding activities. High energy radiation inactivated the glucose-sensitive cytochalasin B binding of each of these membrane preparations by reducing the total number of the binding sites without affecting the dissociation constant. The reduction in total number of binding sites was analyzed as a function of radiation dose based on target theory, from which a radiation-sensitive mass (target size) was calculated. When the plasma membranes of insulin-treated adipocytes were used, a target size of approximately 58,000 daltons was obtained. For adipocyte microsomal membranes, we obtained target sizes of approximately 112,000 and 109,000 daltons prior to and after insulin treatment, respectively. In the case of microsomal membranes, however, inactivation data showed anomalously low radiation sensitivities at low radiation doses, which may be interpreted as indicating the presence of a radiation-sensitive inhibitor. These results suggest that the adipocyte glucose transporter occurs as a monomer in the plasma membrane while existing in the intracellular reserve pool either as a homodimer or as a stoichiometric complex with a protein of an approximately equal size.
A general hybrid radiation transport scheme for star formation simulations on an adaptive grid
Klassen, Mikhail; Pudritz, Ralph E.; Kuiper, Rolf; Peters, Thomas; Banerjee, Robi; Buntemeyer, Lars
2014-12-10
Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodyanmics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.
Dust Transport, Deposition and Radiative Effects Observed from MODIS
NASA Technical Reports Server (NTRS)
Kaufman, Y. J.; Koren, I.; Remer, L. A.; Tanre, D.; Ginoux, P.; Fan, S.
2003-01-01
Carlson (1977) used satellite (AVHRR) observation of dust episodes 3 estimate that 90 tg of dust are emitted from Africa (0-30 N) to the Atlantic Ocean between June and August. MODIS systematic measurements of aerosol optical thickness (AOT) and the fraction of the AOT (f) due to the fine mode (see Remer et al abstract), are used to derive the column concentration, flux and deposition of African dust over the Atlantic Ocean. The main data set is for 2001 but the results are consistent with MODIS measurements from 2002. The analysis first determines the properties of maritime baseline aerosol (AOT=0.06, f=0.5); followed by linear scaling of the dust AOT and the anthropogenic AOT, based on MODIS measured values of the fraction "f" being 0.9 for anthropogenic aerosol and 0.5 for dust. NCEP winds are used in the analysis and are evaluated against observed dust movements between the Terra and Aqua passes (see Koren et al. abstract). Monthly values of dust transport and deposition are calculated. Preliminary results show that 280 tg of dust are emitted annually from Africa to the Atlantic Ocean between 20s and 30N, with 40 tg returning to Africa and Europe between 30N and 50N. 85 tg reach the Americas, with 130-150 tg are deposited in the Atlantic Ocean. The results are compared with dust transport models that indicate 110-230 tg of dust being deposited in the Ocean. It is interesting to note that the early estimates of Carlson (1977) and Carlson & Prosper0 (1972) are very close to our estimate from MODIS of 100 tg for the same latitude range and monthly period.
Transport simulations of ITER with broad density profiles and high radiative fraction
NASA Astrophysics Data System (ADS)
Becker, G.
1995-08-01
Special versions of the 1.5-D BALDUR predictive transport code are used to explore the confinement in the ignited ITER EDA by self-consistent calculations. The code computes 2-D equilibria and solves 1-D transport equations in the bulk and scrape-off layer with empirical transport coefficients for the ohmic, L and ELMy H mode regimes. The emphasis is on scenarios with flat density profiles and high, fixed radiative power in the main chamber due to the seeded impurities argon and neon. It is shown that self-sustained steady state thermonuclear burn is achieved for 370 s and is compatible with the flat density profiles and strong radiative cooling. The necessary local energy and particle transport are presented. In the argon and neon scenarios, the required radiation corrected energy confinement times are 4.1 and 3.5 s, respectively, which are achievable according to the ITER ELMy H mode scaling. The advantage of neon originates from its smaller radiative loss within the separatrix of 37% of the total radiation in the main chamber, compared with 60% in the case of argon. A significant radiative loss from the confinement zone, mainly by bremsstrahlung, cannot be avoided. It raises the required energy confinement time and is the price to be paid for reduction of the divertor heat load by radiative cooling in the main chamber. In steady state, a helium fraction of 5% is computed. The fractions of helium, argon and neon and the resulting fuel dilution are considerably lower than with peaked density profiles
Advantages of Analytical Transformations in Monte Carlo Methods for Radiation Transport
McKinley, M S; Brooks III, E D; Daffin, F
2004-12-13
Monte Carlo methods for radiation transport typically attempt to solve an integral by directly sampling analog or weighted particles, which are treated as physical entities. Improvements to the methods involve better sampling, probability games or physical intuition about the problem. We show that significant improvements can be achieved by recasting the equations with an analytical transform to solve for new, non-physical entities or fields. This paper looks at one such transform, the difference formulation for thermal photon transport, showing a significant advantage for Monte Carlo solution of the equations for time dependent transport. Other related areas are discussed that may also realize significant benefits from similar analytical transformations.
NASA Astrophysics Data System (ADS)
Townsend, L. W.; Porter, J.; Spence, H. E.; Golightly, M. J.; Smith, S. S.; Schwadron, N.; Kasper, J. C.; Case, A. W.; Blake, J. B.; Mazur, J. E.; Looper, M. D.; Zeitlin, C. J.
2014-12-01
The Cosmic Ray Telescope for the Effects of Radiation (CRaTER) instrument on the Lunar Reconnaissance Orbiter (LRO) spacecraft measures the energy depositions by solar and galactic cosmic radiations in its silicon detectors. These energy depositions are converted to linear energy transfer (LET) spectra, which can contribute to benchmarking space radiation transport codes and also used to estimate doses for the Lunar environment. In this work the Monte Carlo transport code HETC-HEDS (High Energy Transport Code - Human Exploration and Development in Space) and the deterministic NASA space radiation transport code HZETRN2010 are used to estimate LET and dose contributions from the incident primary ions and their charged secondaries produced in nuclear collisions within the components of the CRaTER instrument. Comparisons of the calculated LET spectra with measurements of LET from the CRaTER instrument are made and clearly show the importance of including corrections to the calculated average energy deposition spectra in the silicon detectors using a Vavilov distribution function.
Frankel, J.I.
1997-09-01
This investigation used sysmbolic manipulation in developing analytical methods and general computational strategies for solving both linear and nonlinear, regular and singular integral and integro-differential equations which appear in radiative and mixed-mode energy transport. Contained in this report are seven papers which present the technical results as individual modules.
Time-dependent radiation transport in a one-dimensional medium
NASA Technical Reports Server (NTRS)
Nagel, W.; Meszaros, P.
1985-01-01
An analytic solution of the time-dependent radiation transport problem in a one-dimensional, stationary and homogeneous medium of finite thickness is presented. The solution is found by the method of images, and is compared with an eigenfunction expansion. Previous conjectures about the structure of such an expansion are clarified. The Green's function of this problem is also expanded in scattering orders.
Methods of treating complex space vehicle geometry for charged particle radiation transport
NASA Technical Reports Server (NTRS)
Hill, C. W.
1973-01-01
Current methods of treating complex geometry models for space radiation transport calculations are reviewed. The geometric techniques used in three computer codes are outlined. Evaluations of geometric capability and speed are provided for these codes. Although no code development work is included several suggestions for significantly improving complex geometry codes are offered.
Los Alamos National Laboratory Overview
Neu, Mary
2010-06-02
Mary Neu, Associate Director for Chemistry, Life and Earth Sciences at Los Alamos National Laboratory, delivers opening remarks at the "Sequencing, Finishing, Analysis in the Future" meeting in Santa Fe, NM
Warren, Kevin; Reed, Robert; Weller, Robert; Mendenhall, Marcus; Sierawski, Brian; Schrimpf, Ronald
2011-06-01
MRED (Monte Carlo Radiative Energy Deposition) is Vanderbilt University's Geant4 application for simulating radiation events in semiconductors. Geant4 is comprised of the best available computational physics models for the transport of radiation through matter. In addition to basic radiation transport physics contained in the Geant4 core, MRED has the capability to track energy loss in tetrahedral geometric objects, includes a cross section biasing and track weighting technique for variance reduction, and additional features relevant to semiconductor device applications. The crucial element of predicting Single Event Upset (SEU) parameters using radiation transport software is the creation of a dosimetry model that accurately approximates the net collected charge at transistor contacts as a function of deposited energy. The dosimetry technique described here is the multiple sensitive volume (MSV) model. It is shown to be a reasonable approximation of the charge collection process and its parameters can be calibrated to experimental measurements of SEU cross sections. The MSV model, within the framework of MRED, is examined for heavy ion and high-energy proton SEU measurements of a static random access memory.
Radiation transport and energetics of laser-driven half-hohlraums at the National Ignition Facility
Moore, A. S.; Cooper, A. B.R.; Schneider, M. B.; MacLaren, S.; Graham, P.; Lu, K.; Seugling, R.; Satcher, J.; Klingmann, J.; Comley, A. J.; Marrs, R.; May, M.; Widmann, K.; Glendinning, G.; Castor, J.; Sain, J.; Back, C. A.; Hund, J.; Baker, K.; Hsing, W. W.; Foster, J.; Young, B.; Young, P.
2014-06-01
Experiments that characterize and develop a high energy-density half-hohlraum platform for use in bench-marking radiation hydrodynamics models have been conducted at the National Ignition Facility (NIF). Results from the experiments are used to quantitatively compare with simulations of the radiation transported through an evolving plasma density structure, colloquially known as an N-wave. A half-hohlraum is heated by 80 NIF beams to a temperature of 240 eV. This creates a subsonic di usive Marshak wave which propagates into a high atomic number Ta_{2}O_{5} aerogel. The subsequent radiation transport through the aerogel and through slots cut into the aerogel layer is investigated. We describe a set of experiments that test the hohlraum performance and report on a range
NASA Technical Reports Server (NTRS)
Gleckler, P. J.; Randall, D. A.; Boer, G.; Colman, R.; Dix, M.; Galin, V.; Helfand, M.; Kiehl, J.; Kitoh, A.; Lau, W.
1995-01-01
This paper summarizes the ocean surface net energy flux simulated by fifteen atmospheric general circulation models constrained by realistically-varying sea surface temperatures and sea ice as part of the Atmospheric Model Intercomparison Project. In general, the simulated energy fluxes are within the very large observational uncertainties. However, the annual mean oceanic meridional heat transport that would be required to balance the simulated surface fluxes is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean heat transport can be affected by the errors in simulated cloud-radiation interactions. It is suggested that improved treatment of cloud radiative effects should help in the development of coupled atmosphere-ocean general circulation models.
Comparison of Stopping Power and Range Databases for Radiation Transport Study
NASA Technical Reports Server (NTRS)
Tai, H.; Bichsel, Hans; Wilson, John W.; Shinn, Judy L.; Cucinotta, Francis A.; Badavi, Francis F.
1997-01-01
The codes used to calculate stopping power and range for the space radiation shielding program at the Langley Research Center are based on the work of Ziegler but with modifications. As more experience is gained from experiments at heavy ion accelerators, prudence dictates a reevaluation of the current databases. Numerical values of stopping power and range calculated from four different codes currently in use are presented for selected ions and materials in the energy domain suitable for space radiation transport. This study of radiation transport has found that for most collision systems and for intermediate particle energies, agreement is less than 1 percent, in general, among all the codes. However, greater discrepancies are seen for heavy systems, especially at low particle energies.
Radiative or neutron transport modeling using a lattice Boltzmann equation framework
NASA Astrophysics Data System (ADS)
Bindra, H.; Patil, D. V.
2012-07-01
In this paper, the lattice Boltzmann equation (LBE)-based framework is used to obtain the solution for the linear radiative or neutron transport equation. The LBE framework is devised for the integrodifferential forms of these equations which arise due to the inclusion of the scattering terms. The interparticle collisions are neglected, hence omitting the nonlinear collision term. Furthermore, typical representative examples for one-dimensional or two-dimensional geometries and inclusion or exclusion of the scattering term (isotropic and anisotropic) in the Boltzmann transport equation are illustrated to prove the validity of the method. It has been shown that the solution from the LBE methodology is equivalent to the well-known Pn and Sn methods. This suggests that the LBE can potentially provide a more convenient and easy approach to solve the physical problems of neutron and radiation transport.
NASA Astrophysics Data System (ADS)
Smedley-Stevenson, Richard P.; McClarren, Ryan G.
2015-04-01
This paper attempts to unify the asymptotic diffusion limit analysis of thermal radiation transport schemes, for a linear-discontinuous representation of the material temperature reconstructed from cell centred temperature unknowns, in a process known as 'source tilting'. The asymptotic limits of both Monte Carlo (continuous in space) and deterministic approaches (based on linear-discontinuous finite elements) for solving the transport equation are investigated in slab geometry. The resulting discrete diffusion equations are found to have nonphysical terms that are proportional to any cell-edge discontinuity in the temperature representation. Based on this analysis it is possible to design accurate schemes for representing the material temperature, for coupling thermal radiation transport codes to a cell centred representation of internal energy favoured by ALE (arbitrary Lagrange-Eulerian) hydrodynamics schemes.
Smedley-Stevenson, Richard P.; McClarren, Ryan G.
2015-04-01
This paper attempts to unify the asymptotic diffusion limit analysis of thermal radiation transport schemes, for a linear-discontinuous representation of the material temperature reconstructed from cell centred temperature unknowns, in a process known as ‘source tilting’. The asymptotic limits of both Monte Carlo (continuous in space) and deterministic approaches (based on linear-discontinuous finite elements) for solving the transport equation are investigated in slab geometry. The resulting discrete diffusion equations are found to have nonphysical terms that are proportional to any cell-edge discontinuity in the temperature representation. Based on this analysis it is possible to design accurate schemes for representing the material temperature, for coupling thermal radiation transport codes to a cell centred representation of internal energy favoured by ALE (arbitrary Lagrange–Eulerian) hydrodynamics schemes.
Calculation of radiation therapy dose using all particle Monte Carlo transport
Chandler, W.P.; Hartmann-Siantar, C.L.; Rathkopf, J.A.
1999-02-09
The actual radiation dose absorbed in the body is calculated using three-dimensional Monte Carlo transport. Neutrons, protons, deuterons, tritons, helium-3, alpha particles, photons, electrons, and positrons are transported in a completely coupled manner, using this Monte Carlo All-Particle Method (MCAPM). The major elements of the invention include: computer hardware, user description of the patient, description of the radiation source, physical databases, Monte Carlo transport, and output of dose distributions. This facilitated the estimation of dose distributions on a Cartesian grid for neutrons, photons, electrons, positrons, and heavy charged-particles incident on any biological target, with resolutions ranging from microns to centimeters. Calculations can be extended to estimate dose distributions on general-geometry (non-Cartesian) grids for biological and/or non-biological media. 57 figs.
Calculation of radiation therapy dose using all particle Monte Carlo transport
Chandler, William P.; Hartmann-Siantar, Christine L.; Rathkopf, James A.
1999-01-01
The actual radiation dose absorbed in the body is calculated using three-dimensional Monte Carlo transport. Neutrons, protons, deuterons, tritons, helium-3, alpha particles, photons, electrons, and positrons are transported in a completely coupled manner, using this Monte Carlo All-Particle Method (MCAPM). The major elements of the invention include: computer hardware, user description of the patient, description of the radiation source, physical databases, Monte Carlo transport, and output of dose distributions. This facilitated the estimation of dose distributions on a Cartesian grid for neutrons, photons, electrons, positrons, and heavy charged-particles incident on any biological target, with resolutions ranging from microns to centimeters. Calculations can be extended to estimate dose distributions on general-geometry (non-Cartesian) grids for biological and/or non-biological media.
Atmospheric Ionizing Radiation and the High Speed Civil Transport. Chapter 1
NASA Technical Reports Server (NTRS)
Maiden, D. L.; Wilson, J. W.; Jones, I. W.; Goldhagen, P.
2003-01-01
Atmospheric ionizing radiation is produced by extraterrestrial radiations incident on the Earth's atmosphere. These extraterrestrial radiations are of two sources: ever present galactic cosmic rays with origin outside the solar system and transient solar particle events that are at times very intense events associated with solar activity lasting several hours to a few days. Although the galactic radiation penetrating through the atmosphere to the ground is low in intensity, the intensity is more than two orders of magnitude greater at commercial aircraft altitudes. The radiation levels at the higher altitudes of the High Speed Civil Transport (HSCT) are an additional factor of two higher. Ionizing radiation produces chemically active radicals in biological tissues that alter the cell function or result in cell death. Protection standards against low levels of ionizing radiation are based on limitation of excess cancer mortality or limitation of developmental injury resulting in permanent damage to the offspring during pregnancy. The crews of commercial air transport operations are considered as radiation workers by the EPA, the FAA, and the International Commission on Radiological Protection (ICRP). The annual exposures of aircrews depend on the latitudes and altitudes of operation and flight time. Flight hours have significantly increased since deregulation of the airline industry in the 1980's. The FAA estimates annual subsonic aircrew exposures to range from 0.2 to 9.1 mSv compared to 0.5 mSv exposure of the average nuclear power plant worker in the nuclear industry. The commercial aircrews of the HSCT may receive exposures above recently recommended allowable limits for even radiation workers if flying their allowable number of flight hours. An adequate protection philosophy for background exposures in HSCT commercial airtraffic cannot be developed at this time due to current uncertainty in environmental levels. In addition, if a large solar particle event
NASA Technical Reports Server (NTRS)
Miller, Jack; Heilbronn, Lawrence H.; Zeitlin, Cary J.; Wilson, John W.; Singleterry, Robert C., Jr.; Thibeault, Sheila Ann
2003-01-01
Mission crews in space outside the Earth s magnetic field will be exposed to high energy heavy charged particles in the galactic cosmic radiation (GCR). These highly ionizing particles will be a source of radiation risk to crews on extended missions to the Moon and Mars, and the biological effects of and countermeasures to the GCR have to be investigated as part of the planning of exploration-class missions. While it is impractical to shield spacecraft and planetary habitats against the entire GCR spectrum, biological and physical studies indicate that relatively modest amounts of shielding are effective at reducing the radiation dose. However, nuclear fragmentation in the shielding materials produces highly penetrating secondary particles, which complicates the problem: in some cases, some shielding is worse than none at all. Therefore the radiation transport properties of potential shielding materials need to be carefully investigated. One intriguing option for a Mars mission is the use of material from the Martian surface, in combination with chemicals carried from Earth and/or fabricated from elements found in the Martian atmosphere, to construct crew habitats. We have measured the transmission properties of epoxy-Martian regolith composites with respect to heavy charged particles characteristic of the GCR ions which bombard the Martian surface. The composites were prepared at NASA Langley Research Center using simulated Martian regolith, in the process also evaluating fabrication methods which could lead to technologies for in situ fabrication on Mars. Initial evaluation of the radiation shielding properties is made using radiation transport models developed at NASA-LaRC, and the results of these calculations are used to select the composites with the most favorable radiation transmission properties. These candidates are then evaluated at particle accelerators which produce beams of heavy charged particles representative in energy and charge of the radiation
NASA Technical Reports Server (NTRS)
Davis, Anthony B.; Marshak, Alexander
2010-01-01
The interplay of sunlight with clouds is a ubiquitous and often pleasant visual experience, but it conjures up major challenges for weather, climate, environmental science and beyond. Those engaged in the characterization of clouds (and the clear air nearby) by remote sensing methods are even more confronted. The problem comes, on the one hand, from the spatial complexity of real clouds and, on the other hand, from the dominance of multiple scattering in the radiation transport. The former ingredient contrasts sharply with the still popular representation of clouds as homogeneous plane-parallel slabs for the purposes of radiative transfer computations. In typical cloud scenes the opposite asymptotic transport regimes of diffusion and ballistic propagation coexist. We survey the three-dimensional (3D) atmospheric radiative transfer literature over the past 50 years and identify three concurrent and intertwining thrusts: first, how to assess the damage (bias) caused by 3D effects in the operational 1D radiative transfer models? Second, how to mitigate this damage? Finally, can we exploit 3D radiative transfer phenomena to innovate observation methods and technologies? We quickly realize that the smallest scale resolved computationally or observationally may be artificial but is nonetheless a key quantity that separates the 3D radiative transfer solutions into two broad and complementary classes: stochastic and deterministic. Both approaches draw on classic and contemporary statistical, mathematical and computational physics.
A hybrid transport-diffusion model for radiative transfer in absorbing and scattering media
Roger, M.; Caliot, C.; Crouseilles, N.; Coelho, P.J.
2014-10-15
A new multi-scale hybrid transport-diffusion model for radiative transfer is proposed in order to improve the efficiency of the calculations close to the diffusive regime, in absorbing and strongly scattering media. In this model, the radiative intensity is decomposed into a macroscopic component calculated by the diffusion equation, and a mesoscopic component. The transport equation for the mesoscopic component allows to correct the estimation of the diffusion equation, and then to obtain the solution of the linear radiative transfer equation. In this work, results are presented for stationary and transient radiative transfer cases, in examples which concern solar concentrated and optical tomography applications. The Monte Carlo and the discrete-ordinate methods are used to solve the mesoscopic equation. It is shown that the multi-scale model allows to improve the efficiency of the calculations when the medium is close to the diffusive regime. The proposed model is a good alternative for radiative transfer at the intermediate regime where the macroscopic diffusion equation is not accurate enough and the radiative transfer equation requires too much computational effort.
NASA Astrophysics Data System (ADS)
Davis, Anthony B.; Marshak, Alexander
2010-02-01
The interplay of sunlight with clouds is a ubiquitous and often pleasant visual experience, but it conjures up major challenges for weather, climate, environmental science and beyond. Those engaged in the characterization of clouds (and the clear air nearby) by remote sensing methods are even more confronted. The problem comes, on the one hand, from the spatial complexity of real clouds and, on the other hand, from the dominance of multiple scattering in the radiation transport. The former ingredient contrasts sharply with the still popular representation of clouds as homogeneous plane-parallel slabs for the purposes of radiative transfer computations. In typical cloud scenes the opposite asymptotic transport regimes of diffusion and ballistic propagation coexist. We survey the three-dimensional (3D) atmospheric radiative transfer literature over the past 50 years and identify three concurrent and intertwining thrusts: first, how to assess the damage (bias) caused by 3D effects in the operational 1D radiative transfer models? Second, how to mitigate this damage? Finally, can we exploit 3D radiative transfer phenomena to innovate observation methods and technologies? We quickly realize that the smallest scale resolved computationally or observationally may be artificial but is nonetheless a key quantity that separates the 3D radiative transfer solutions into two broad and complementary classes: stochastic and deterministic. Both approaches draw on classic and contemporary statistical, mathematical and computational physics.
A hybrid transport-diffusion model for radiative transfer in absorbing and scattering media
NASA Astrophysics Data System (ADS)
Roger, M.; Caliot, C.; Crouseilles, N.; Coelho, P. J.
2014-10-01
A new multi-scale hybrid transport-diffusion model for radiative transfer is proposed in order to improve the efficiency of the calculations close to the diffusive regime, in absorbing and strongly scattering media. In this model, the radiative intensity is decomposed into a macroscopic component calculated by the diffusion equation, and a mesoscopic component. The transport equation for the mesoscopic component allows to correct the estimation of the diffusion equation, and then to obtain the solution of the linear radiative transfer equation. In this work, results are presented for stationary and transient radiative transfer cases, in examples which concern solar concentrated and optical tomography applications. The Monte Carlo and the discrete-ordinate methods are used to solve the mesoscopic equation. It is shown that the multi-scale model allows to improve the efficiency of the calculations when the medium is close to the diffusive regime. The proposed model is a good alternative for radiative transfer at the intermediate regime where the macroscopic diffusion equation is not accurate enough and the radiative transfer equation requires too much computational effort.
Badhwar, G.D.; Konradi, A.; Cucinotta, F.A.; Braby, L.A.
1994-09-01
A new class of tissue-equivalent proportional counters has been flown on two space shuttle flights. These detectors and their associated electronics cover a lineal energy range from 0.4 to 1250 keV/{mu}m with a multichannel analyzer resolution of 0.1 keV/{mu}m from 0.4 to 20 keV/{mu} and 5 keV/{mu}m from 20 to 1250 keV/{mu}m. These detectors provide the most complete dynamic range and highest resolution of any technique currently in use. On one mission, one detector was mounted in the Shuttle payload bay and another older modele in the mid-deck, thus providing information on the depth dependence of the lineal energy spectrum. A detailed comparison of the observed lineal energy and calculated LET spectra for galacic cosmic radiation shows that, although the radiation transport models provide a rather accurate description of the dose ({+-}15%) and equivalent dose ({+-}15%), the calculations significantly underestimate the frequency of events below about 100 keV/{mu}m. This difference cannot be explained by the inclusion of the contribution of splash protons. The contribution of the secondary pions, kaons and electrons produce in the Shuttle shielding, if included in the radiation transport model, may explain these differences. There are also significant differences between the model predictions and observations above 1440 keV/{mu}m, particularly for 28.5{degrees} inclination orbit. 24 refs., 9 figs., 1 tab.
On the Transport and Radiative Properties of Plasmas with Small-Scale Electromagnetic Fluctuations
NASA Astrophysics Data System (ADS)
Keenan, Brett D.
Plasmas with sub-Larmor-scale ("small-scale") electromagnetic fluctuations are a feature of a wide variety of high-energy-density environments, and are essential to the description of many astrophysical/laboratory plasma phenomena. Radiation from particles, whether they be relativistic or non-relativistic, moving through small-scale electromagnetic turbulence has spectral characteristics distinct from both synchrotron and cyclotron radiation. The radiation, carrying information on the statistical properties of the turbulence, is also intimately related to the particle diffusive transport. We investigate, both theoretically and numerically, the transport of non-relativistic and transrelativistic particles in plasmas with high-amplitude isotropic sub-Larmor-scale magnetic turbulence---both with and without a mean field component---and its relation to the spectra of radiation simultaneously produced by these particles. Furthermore, the transport of particles through small-scale electromagnetic turbulence---under certain conditions---resembles the random transport of particles---via Coulomb collisions---in collisional plasmas. The pitch-angle diffusion coefficient, which acts as an effective "collision" frequency, may be substantial in these, otherwise, collisionless environments. We show that this effect, colloquially referred to as the plasma "quasi-collisionality", may radically alter the expected radiative transport properties of candidate plasmas. We argue that the modified magneto-optic effects in these plasmas provide an attractive, novel, diagnostic tool for the exploration and characterization of small-scale electromagnetic turbulence. Lastly, we speculate upon the manner in which quasi-collisions may affect inertial confinement fusion (ICF), and other laser-plasma experiments. Finally, we show that mildly relativistic jitter radiation, from laser-produced plasmas, may offer insight into the underlying electromagnetic turbulence. Here we investigate the
Grant, K.E.; Taylor, K.E.; Ellis, J.S.; Wuebbles, D.J.
1987-07-01
The authors have implemented a series of state of the art radiation transport submodels in previously developed one dimensional and two dimensional chemical transport models of the troposphere and stratosphere. These submodels provide the capability of calculating accurate solar and infrared heating rates. They are a firm basis for further radiation submodel development as well as for studying interactions between radiation and model dynamics under varying conditions of clear sky, clouds, and aerosols. 37 refs., 3 figs.
Environmental surveillance at Los Alamos during 1979
Not Available
1980-04-01
This report documents the environmental surveillance program conducted by the Los Alamos Scientific Laboratory (LASL) in 1979. Routine monitoring for radiation and radioactive or chemical substances was conducted on the Laboratory site and in the surrounding region to determine compliance with appropriate standards and permit early identification of possible undesirable trends. Results and interpretation of the data for 1979 on penetrating radiation, chemical and radiochemical quality of ambient air, surface and ground water, municipal water supply, soils and sediments, food, and airborne and liquid effluents are included. Comparisons with appropriate standards and regulations or with background levels from natural or other non-LASL sources provide a basis for concluding that environmental effects attributable to LASL operations are minor and cannot be considered likely to result in any hazard to the population of the area. Results of several special studies provide documentation of some unique environmental conditions in the LASL environs.
Environmental surveillance at Los Alamos during 1992
Kohen, K.; Stoker, A.; Stone, G.
1994-07-01
This report describes the environmental surveillance program at Los Alamos National Laboratory during 1992. The Laboratory routinely monitors for radiation and for radioactive and nonradioactive materials at (or on) Laboratory sites as well as in the surrounding region. LANL uses the monitoring results to determine compliance with appropriate standards and to identify potentially undesirable trends. Data were collected in 1992 to assess external penetrating radiation; quantities of airborne emissions and liquid effluents; concentrations of chemicals and radionuclides in ambient air, surface waters and groundwaters, municipal water supply, soils and sediments, and foodstuffs; and environmental compliance. Using comparisons with standards, regulations, and background levels, this report concludes that environmental effects from Laboratory operations are small and do not pose a demonstrable threat to the public, laboratory employees, or the environment.
Environmental surveillance at Los Alamos during 1989
Not Available
1990-12-01
This report describes the environmental surveillance program conducted by Los Alamos National Laboratory during 1989. Routine monitoring for radiation and radioactive or chemical materials is conducted on the Laboratory site as well as in the surrounding region. Monitoring results are used to determine compliance with appropriate standards and to permit early identification of potentially undesirable trends. Results and interpretation of data for 1989 cover external penetrating radiation; quantities of airborne emissions and effluents; concentrations of chemicals and radionuclides in ambient air, surface and ground waters, municipal water supply, soils and sediments, and foodstuffs; and environmental compliance. Comparisons with appropriate standards, regulations, and background levels provide the basis for concluding that environmental effects from Laboratory operations are small and do not pose a threat to the public, Laboratory employees, or the environment. 58 refs., 31 figs., 39 tabs.
Environmental surveillance at Los Alamos during 1995
1996-10-01
This report describes the environmental surveillance program at Los Alamos National Laboratory (LANL or the Laboratory) during 1995. The Laboratory routinely monitors for radiation and for radioactive and nonradioactive materials at (or on) Laboratory sites as well as in the surrounding region. LANL uses the monitoring result to determine compliance with appropriate standards and to identify potentially undesirable trends. Data were collected in 1995 to assess external penetrating radiation; quantities of airborne emissions and liquid effluents; concentrations of chemicals and radionuclides in ambient air, surface waters and groundwaters, municipal water supply, soils and sediments, and foodstuffs; and environmental compliance. Using comparisons with standards, regulations, and background levels, this report concludes that environmental effects from Laboratory operations are small and do not pose a demonstrable threat to the public, Laboratory employees, or the environment.
Environmental surveillance at Los Alamos during 1987
Not Available
1988-05-01
This report describes the environmental surveillance program conducted by Los Alamos National Laboratory during 1987. Routine monitoring for radiation and radioactive or chemical materials is conducted on the Laboratory site as well as in the surrounding region. Monitoring results are used to determine compliance with appropriate standards and to permit early identification of potentially undesirable trends. Results and interpretation of data for 1987 cover: external penetrating radiation; quantities of airborne emissions and liquid effluents; concentrations of chemicals and radionuclides in ambient air, surface and ground waters, municipal water supply, soils and sediments, and foodstuffs; and environmental compliance. Comparisons with appropriate standards, regulations, and background levels provide the basis for concluding that environmental effects from Laboratory operations are insignificant and do not pose a threat to the public, Laboratory employees, or the environment. 113 refs., 33 figs., 120 tabs.
Radiation Effects on Sorption and Mobilization of Radionuclides during Transport through Geosphere
Wang, L.M.; Ewing, R.C.; Hayes, K.F.
2002-03-14
Site restoration activities at DOE facilities and the permanent disposal of nuclear waste generated at DOE facilities involve working with and within various types and levels of radiation fields. Radionuclide decay and the associated radiation fields lead to physical and chemical changes that can degrade or enhance material properties. The principal sources of radiation at the DOE sites are the actinides and fission-products contained in high-level wastes currently in storage. Alpha-decay of the actinide elements and beta-decay of the fission products lead to atomic scale changes in the material (radiation damage and transmutation). During site restoration, materials will be exposed to radiation fields that exceed 104 rad/hr. The radiation exposure due to the release and sorption of long-lived actinides (e.g., 237Np) and fission products (e.g., 137Cs and 90Sr) may cause changes in important properties (e.g., cation exchange capacity) in geological materials (e.g., clays and zeolites) along transport pathways. Among these materials, clays and zeolites, which are expected to sorb and immobilize radionuclides, are known to be extremely susceptible to radiation-induced structure changes (e.g., bubble formation and solid state amorphization) through both collisional displacement and ionization processes. These changes will inevitably affect (either negatively or positively) the further sorption and the migration of radionuclides at waste sites (e.g., vadose zone at Hanford). Current models used for the longterm prediction of radionuclide transport have proven to be inadequate and unrealistic; however, these previous models did not take radiation effects into consideration.
The influence of resonance radiation transport on the contraction of a glow discharge in argon
NASA Astrophysics Data System (ADS)
Golubovskii, Yu B.; Maiorov, V. A.
2015-04-01
The role of resonance radiation transport in the contraction of a positive column in an argon glow discharge is studied numerically. The theory is based on the self-consistent solution of the ambipolar diffusion equation for electrons, the diffusion equation for metastable atoms and the Biberman-Holstein equation for resonance atoms. To calculate the ionization and excitation rates, the Boltzmann equation is solved in a local approximation taking into account elastic, inelastic and electron-electron collisions. A solution method for a boundary problem is developed which allows one to obtain a hysteresis of the parameters during a continuous transition from a diffuse mode to a contracted mode through an unstable branch. At small currents there is a diffuse discharge where the role of radiation transport is inessential because the radial distributions of electrons and excited atoms are close to the fundamental modes of the corresponding equations. Under these conditions, the traditional approximation of ‘effective lifetime’ is accurate enough. For a contracted discharge, this approximation is not applicable because the higher diffusion and radiation modes play a notable role and a more strict description of radiation transport is required. It is shown that, when radiation transport is taken into account, the width of a filament in a contracted discharge significantly exceeds that obtained in the traditional ‘effective lifetime’ approximation. The critical current, when the discharge abruptly turns into a contracted mode, is shifted towards higher current values. The results obtained in this paper can also relate to a discharge in other inert gases.
NASA Astrophysics Data System (ADS)
Bouchard, Hugo; Bielajew, Alex
2015-07-01
To establish a theoretical framework for generalizing Monte Carlo transport algorithms by adding external electromagnetic fields to the Boltzmann radiation transport equation in a rigorous and consistent fashion. Using first principles, the Boltzmann radiation transport equation is modified by adding a term describing the variation of the particle distribution due to the Lorentz force. The implications of this new equation are evaluated by investigating the validity of Fano’s theorem. Additionally, Lewis’ approach to multiple scattering theory in infinite homogeneous media is redefined to account for the presence of external electromagnetic fields. The equation is modified and yields a description consistent with the deterministic laws of motion as well as probabilistic methods of solution. The time-independent Boltzmann radiation transport equation is generalized to account for the electromagnetic forces in an additional operator similar to the interaction term. Fano’s and Lewis’ approaches are stated in this new equation. Fano’s theorem is found not to apply in the presence of electromagnetic fields. Lewis’ theory for electron multiple scattering and moments, accounting for the coupling between the Lorentz force and multiple elastic scattering, is found. However, further investigation is required to develop useful algorithms for Monte Carlo and deterministic transport methods. To test the accuracy of Monte Carlo transport algorithms in the presence of electromagnetic fields, the Fano cavity test, as currently defined, cannot be applied. Therefore, new tests must be designed for this specific application. A multiple scattering theory that accurately couples the Lorentz force with elastic scattering could improve Monte Carlo efficiency. The present study proposes a new theoretical framework to develop such algorithms.
Bouchard, Hugo; Bielajew, Alex
2015-07-01
To establish a theoretical framework for generalizing Monte Carlo transport algorithms by adding external electromagnetic fields to the Boltzmann radiation transport equation in a rigorous and consistent fashion. Using first principles, the Boltzmann radiation transport equation is modified by adding a term describing the variation of the particle distribution due to the Lorentz force. The implications of this new equation are evaluated by investigating the validity of Fano's theorem. Additionally, Lewis' approach to multiple scattering theory in infinite homogeneous media is redefined to account for the presence of external electromagnetic fields. The equation is modified and yields a description consistent with the deterministic laws of motion as well as probabilistic methods of solution. The time-independent Boltzmann radiation transport equation is generalized to account for the electromagnetic forces in an additional operator similar to the interaction term. Fano's and Lewis' approaches are stated in this new equation. Fano's theorem is found not to apply in the presence of electromagnetic fields. Lewis' theory for electron multiple scattering and moments, accounting for the coupling between the Lorentz force and multiple elastic scattering, is found. However, further investigation is required to develop useful algorithms for Monte Carlo and deterministic transport methods. To test the accuracy of Monte Carlo transport algorithms in the presence of electromagnetic fields, the Fano cavity test, as currently defined, cannot be applied. Therefore, new tests must be designed for this specific application. A multiple scattering theory that accurately couples the Lorentz force with elastic scattering could improve Monte Carlo efficiency. The present study proposes a new theoretical framework to develop such algorithms. PMID:26061045
Olson, Gordon L.
2015-09-24
One-dimensional models for the transport of radiation through binary stochastic media do not work in multi-dimensions. In addition, authors have attempted to modify or extend the 1D models to work in multidimensions without success. Analytic one-dimensional models are successful in 1D only when assuming greatly simplified physics. State of the art theories for stochastic media radiation transport do not address multi-dimensions and temperature-dependent physics coefficients. Here, the concept of effective opacities and effective heat capacities is found to well represent the ensemble averaged transport solutions in cases with gray or multigroup temperature-dependent opacities and constant or temperature-dependent heat capacities. Inmore » every case analyzed here, effective physics coefficients fit the transport solutions over a useful range of parameter space. The transport equation is solved with the spherical harmonics method with angle orders of n=1 and 5. Although the details depend on what order of solution is used, the general results are similar, independent of angular order.« less
Schneider, K.J.; Hostick, C.J.; Ross, W.A.; Peterson, R.W.; Smith, R.I.; Stiles, D.L.; Daling, P.M.; Weakley, S.A.; Grinde, R.B.; Young, J.R.
1987-11-01
This report contains a system study of estimated radiation doses to the public and workers resulting from the transport of spent fuel from commercial nuclear power reactors to a geologic repository. The report contains a detailed breakdown of activities and a description of time/distance/dose-rate estimates for each activity within the system. Collective doses are estimated for each of the major activities at the reactor site, in transit, and at the repository receiving facility. Annual individual doses to the maximally exposed individuals or groups of individuals are also estimated. A total of 17 alternatives and subalternatives to the postulated reference transportation system are identified, conceptualized, and their dose-reduction potentials and costs estimated. Resulting ratios of ..delta..cost/..delta..collective system dose for each alternative relative to the postulated reference transportation system are given. Most of the alternatives evaluated are estimated to provide both cost and dose reductions. Major reductions in transportation system dose and cost are estimated to result from using higher-capacity rail and truck casks, and particularly when replacing legalweight truck casks with ''advanced design'' overweight truck casks. The greatest annual dose reduction to the highest exposed individual workers (i.e., at the repository) is estimated to be achieved by using remote handling equipment for the cask handling operations at the repository. Additional shielding is also effective in reducing doses to both radiation workers at the reactor and repository and to transport workers. 69 refs., 36 figs., 156 tabs.
Comparing calvarial transport distraction with and without radiation and fat grafting.
Yuhasz, Mikell M; Koch, Felix P; Kwiatkowski, Anna; Young, Calvin; Clune, James; Travieso, Rob; Wong, Kenneth; Van Houten, Joshua; Steinbacher, Derek M
2014-10-01
The purpose of this study is to: a) assess transport distraction to reconstruct cranial defects in radiated and non-radiated fields b) examine adipose grafting's effect on the bony regenerate and overlying wound, and c) elucidate sources of bone formation during transport distraction osteogenesis. Twenty-three male New Zealand white rabbits (3 months; 3.5 kg) were used, 10 non-irradiated and 13 irradiated (17 treatment, 6 control) with a one-time fraction of 35 Gy. A 16 × 16 mm defect was abutted by a 10 × 16 mm transport disc 5 weeks after irradiation, and 11 animals were fat grafted at the distraction site. Latency (1 day), distraction (1.5 mm/day), and consolidation (4 weeks) followed. Fluorochromes were injected subcutaneously and microCT, fluorescence, and histology assessed. In distracted animals without fat grafting, bone density measured 701.87 mgHA/ccm and 2271.95 mgHA/ccm in irradiated and non-irradiated animals. In distracted animals with fat grafting, bone density measured 703.23 mgHA/ccm and 2254.27 mgHA/ccm in irradiated and non-irradiated animals. Fluorescence revealed ossification emanating from the dura, periosteum, and transport segment with decreased formation in irradiated animals. Transport distraction is possible for cranial reconstruction in irradiated fields but short-term osseous fill is significantly diminished. Adipose grafting enhances wound healing in previously irradiated fields but does not enhance ossification. PMID:24864072
Olson, Gordon L.
2015-09-24
One-dimensional models for the transport of radiation through binary stochastic media do not work in multi-dimensions. In addition, authors have attempted to modify or extend the 1D models to work in multidimensions without success. Analytic one-dimensional models are successful in 1D only when assuming greatly simplified physics. State of the art theories for stochastic media radiation transport do not address multi-dimensions and temperature-dependent physics coefficients. Here, the concept of effective opacities and effective heat capacities is found to well represent the ensemble averaged transport solutions in cases with gray or multigroup temperature-dependent opacities and constant or temperature-dependent heat capacities. In every case analyzed here, effective physics coefficients fit the transport solutions over a useful range of parameter space. The transport equation is solved with the spherical harmonics method with angle orders of n=1 and 5. Although the details depend on what order of solution is used, the general results are similar, independent of angular order.
Los Alamos High-Brightness Accelerator FEL (HIBAF) facility
Cornelius, W.D.; Bender, S.; Meier, K.; Thode, L.E.; Watson, J.M.
1989-01-01
The 10-/mu/m Los Alamos free-electron laser (FEL) facility is being upgraded. The conventional electron gun and bunchers have been replaced with a much more compact 6-MeV photoinjector accelerator. By adding existing parts from previous experiments, the primary beam energy will be doubled to 40 MeV. With the existing 1-m wiggler (/lambda//sub w/ = 2.7 cm) and resonator, the facility can produce photons with wavelengths from 3 to 100 /mu/m when lasing on the fundamental mode and produce photons in the visible spectrum with short-period wigglers or harmonic operation. After installation of a 150/degree/ bend, a second wiggler will be added as an amplifier. The installation of laser transport tubes between the accelerator vault and an upstairs laboratory will provide experimenters with a radiation-free environment for experiments. Although the initial experimental program of the upgraded facility will be to test the single accelerator-master oscillator/power amplifier configuration, some portion of the operational time of the facility can be dedicated to user experiments. 13 refs., 5 figs., 6 tabs.
A NEW MONTE CARLO METHOD FOR TIME-DEPENDENT NEUTRINO RADIATION TRANSPORT
Abdikamalov, Ernazar; Ott, Christian D.; O'Connor, Evan; Burrows, Adam; Dolence, Joshua C.; Loeffler, Frank; Schnetter, Erik
2012-08-20
Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck and Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.
Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.; Davidson, Gregory G.; Hamilton, Steven P.; Godfrey, Andrew T.
2015-12-21
This paper discusses the implementation, capabilities, and validation of Shift, a massively parallel Monte Carlo radiation transport package developed and maintained at Oak Ridge National Laboratory. It has been developed to scale well from laptop to small computing clusters to advanced supercomputers. Special features of Shift include hybrid capabilities for variance reduction such as CADIS and FW-CADIS, and advanced parallel decomposition and tally methods optimized for scalability on supercomputing architectures. Shift has been validated and verified against various reactor physics benchmarks and compares well to other state-of-the-art Monte Carlo radiation transport codes such as MCNP5, CE KENO-VI, and OpenMC. Somemore » specific benchmarks used for verification and validation include the CASL VERA criticality test suite and several Westinghouse AP1000® problems. These benchmark and scaling studies show promising results.« less
Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.; Davidson, Gregory G.; Hamilton, Steven P.; Godfrey, Andrew T.
2015-12-21
This paper discusses the implementation, capabilities, and validation of Shift, a massively parallel Monte Carlo radiation transport package developed and maintained at Oak Ridge National Laboratory. It has been developed to scale well from laptop to small computing clusters to advanced supercomputers. Special features of Shift include hybrid capabilities for variance reduction such as CADIS and FW-CADIS, and advanced parallel decomposition and tally methods optimized for scalability on supercomputing architectures. Shift has been validated and verified against various reactor physics benchmarks and compares well to other state-of-the-art Monte Carlo radiation transport codes such as MCNP5, CE KENO-VI, and OpenMC. Some specific benchmarks used for verification and validation include the CASL VERA criticality test suite and several Westinghouse AP1000^{®} problems. These benchmark and scaling studies show promising results.
Functional Size of Photosynthetic Electron Transport Chain Determined by Radiation Inactivation 1
Pan, Run Sun; Chien, Lee Feng; Wang, May Yun; Tsai, Mai Yu; Pan, Rong Long; Hsu, Ban Dar
1987-01-01
Radiation inactivation technique was employed to determine the functional size of photosynthetic electron transport chain of spinach chloroplasts. The functional size for photosystem I+II (H2O to methylviologen) was 623 ± 37 kilodaltons; for photosystem II (H2O to dimethylquinone/ferricyanide), 174 ± 11 kilodaltons; and for photosystem I (reduced diaminodurene to methylviologen), 190 ± 11 kilodaltons. The difference between 364 ± 22 (the sum of 174 ± 11 and 190 ± 11) kilodaltons and 623 ± 37 kilodaltons is partially explained to be due to the presence of two molecules of cytochrome b6/f complex of 280 kilodaltons. The molecular mass for other partial reactions of photosynthetic electron flow, also measured by radiation inactivation, is reported. The molecular mass obtained by this technique is compared with that determined by other conventional biochemical methods. A working hypothesis for the composition, stoichiometry, and organization of polypeptides for photosynthetic electron transport chain is proposed. PMID:16665649
Radice, David; Abdikamalov, Ernazar; Rezzolla, Luciano; Ott, Christian D.
2013-06-01
Recent work by McClarren and Hauck (2010) [31] suggests that the filtered spherical harmonics method represents an efficient, robust, and accurate method for radiation transport, at least in the two-dimensional (2D) case. We extend their work to the three-dimensional (3D) case and find that all of the advantages of the filtering approach identified in 2D are present also in the 3D case. We reformulate the filter operation in a way that is independent of the timestep and of the spatial discretization. We also explore different second- and fourth-order filters and find that the second-order ones yield significantly better results. Overall, our findings suggest that the filtered spherical harmonics approach represents a very promising method for 3D radiation transport calculations.
NASA Technical Reports Server (NTRS)
Hartung, Lin C.; Hassan, H. A.
1992-01-01
A moment method for computing 3-D radiative transport is applied to axisymmetric flows in thermochemical nonequilibrium. Such flows are representative of proposed aerobrake missions. The method uses the P-1 approximation to reduce the governing system of integro-di erential equations to a coupled set of partial di erential equations. A numerical solution method for these equations given actual variations of the radiation properties in thermochemical nonequilibrium blunt body flows is developed. Initial results from the method are shown and compared to tangent slab calculations. The agreement between the transport methods is found to be about 10 percent in the stagnation region, with the difference increasing along the flank of the vehicle.
Diagrammatic and asymptotic approaches to the origins of radiative transport theory: tutorial.
Cazé, A; Schotland, John C
2015-08-01
The radiative transport equation (RTE) is used widely to describe the propagation of multiply scattered light in disordered media. In this tutorial, we present two derivations of the RTE for scalar wave fields. The first derivation is based on diagrammatic perturbation theory, while the second stems from an asymptotic multiscale expansion. Although the two approaches are quite distinct mathematically, some common ground can be found and is discussed. PMID:26367292
A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes
NASA Technical Reports Server (NTRS)
Schnittman, Jeremy David; Krolik, Julian H.
2013-01-01
We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption mechanisms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code and also present results from a number of numerical tests to assess its accuracy and convergence.
A Monte Carlo Code for Relativistic Radiation Transport around Kerr Black Holes
NASA Astrophysics Data System (ADS)
Schnittman, Jeremy D.; Krolik, Julian H.
2013-11-01
We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption mechanisms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code and also present results from a number of numerical tests to assess its accuracy and convergence.
Visualization of a Deterministic Radiation Transport Model Using Standard Visualization Tools
James A. Galbraith; L. Eric Greenwade
2004-05-01
Output from a deterministic radiation transport code running on a CRAY SV1 is imported into a standard distributed, parallel, visualization tool for analysis. Standard output files, consisting of tetrahedral meshes, are imported to the visualization tool through the creation of a application specific plug-in module. Visualization samples are included, providing visualization of steady state results. Different plot types and operators are utilized to enhance the analysis and assist in reporting the results of the analysis.
A MONTE CARLO CODE FOR RELATIVISTIC RADIATION TRANSPORT AROUND KERR BLACK HOLES
Schnittman, Jeremy D.; Krolik, Julian H. E-mail: jhk@pha.jhu.edu
2013-11-01
We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption mechanisms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code and also present results from a number of numerical tests to assess its accuracy and convergence.
Extreme-ultraviolet radiation transport in small scale length laser-produced tin plasmas
NASA Astrophysics Data System (ADS)
Sequoia, Kevin Lamar Williams
The majority of the studies on laser-produced plasmas as an efficient extreme ultraviolet (EUV) light source have focused on relatively large plasmas produced at large laser facilities. However, to develop a commercially viable light source for EUV lithography, much smaller lasers and hence much smaller plasmas must be employed. Smaller plasmas behave quite differently than large plasmas in that the temperature and density are less uniform, and lateral expansion is more important. These differences affect the energy transport and, in particular, the radiation transport. This work studies the EUV radiation transport in small scale length tin plasmas, focusing on the effects of target geometry and laser pulse duration. Both planar and spherical tin targets were irradiated with an Nd:YAG laser operating at 1.064 microm. Conversion efficiency of laser light to 13.5 nm radiation (in-band), EUV emission spectrum, two-dimensional in-band emission profile, and the plasma electron density were measured experimentally. These measurements provide insight into where the laser is absorbed, where the in-band emission is produced, and how the radiation is transmitted. The plasma evolution in these experiments were simulated with a two-dimensional radiation hydrodynamic code, while the radiation transport and atomic kinetics where modeled with a collisional radiative code. Additional experiments were conducted using planar targets where the pulse duration was varied from 0.5 ns to 16 ns to understand the effects of laser pulse duration. It was found that the optimum plasma temperature for efficient generation and transmission of in-band emission is 20 eV. This is lower than the previously reported optimum temperature of 30 eV. The use of a 1.064 microm heating laser results in overheating of the plasma in a region that is much too dense to transmit the in-band emission. This overheating is necessary for the plasma to reach the optimum temperature in the region where the density is
Barattino, W.J.
1985-07-01
A coupled radiation transport-heat transfer-stress analysis of the radiation shield for an SP-100 reactor was performed using a numerical code developed at the University of New Mexico and Sandia National Laboratory. For a fast reactor operating at 1.66 MW/sub th/, the energy deposited and resulting temperature distribution was determined for a shield consisting of tungsten and lithium hydride pressed into a stainless steel honeycomb matrix. While temperature feedback was shown to have a minor effect on energy deposition, the shielding configuration was found to have a major influence in meeting thermal requirements of the lithium hydride. It was shown that a shield optimized for radiation protection will fail because of melting. However, with minor modifications in the shield layering and material selection, the thermal integrity of the shield can be preserved. A shield design of graphite, depleted lithium hydride, tungsten, and natural lithium hydride was shown to satisfy neutron and gamma fluence requirements, maximum temperature limits, and minimize cracking in the LiH portion of the shield.
A coupled radiation transport-thermal analysis of the radiation shield for an SP-100 type reactor
NASA Astrophysics Data System (ADS)
Barattino, William J.; El-Genk, Mohamed S.; McDaniel, Patrick J.
A coupled radiation transport-thermal analysis of the radiation shield for an SP-100 reactor was performed using finite element codes developed at the University of New Mexico and Sandia National Laboratories. For a fast reactor operating at 1.66 MWt, the energy deposited and resulting temperature distribution were determined for a shield consisting of tungsten and lithium hydride pressed into a stainless steel honeycomb matrix. While temperature feedback was shown to have a minor effect on energy deposition, the shielding configuration was found to have a major influence in meeting thermal requirements of the lithium hydride. It was shown that a shield optimized only for radiation protection will fail because of LiH melting. However, with minor modifications in the shield layering and material selection, the thermal integrity of the shield can be preserved. A shield design of graphite, depleted lithium hydride, tungsten, and natural lithium hydride was shown to satisfy neutron and gamma fluence requirements, and maximum temperature limits, and to minimize cracking in the LiH portion of the shield.
Discontinuous Galerkin Methods for the Two-Moment Model of Radiation Transport
NASA Astrophysics Data System (ADS)
Endeve, Eirik; Hauck, Cory
2016-03-01
We are developing computational methods for simulation of radiation transport in astrophysical systems (e.g., neutrino transport in core-collapse supernovae). Here we consider the two-moment model of radiation transport, where the energy density E and flux F - angular moments of a phase space distribution function - approximates the radiation field in a computationally tractable manner. We aim to develop multi-dimensional methods that are (i) high-order accurate for computational efficiency, and (ii) robust in the sense that the solution remains in the realizable set R = { (E , F) | E >= 0 and E - | F | >= 0 } (i.e., E and F are consistent with moments of an underlying distribution). Our approach is based on the Runge-Kutta discontinuous Galerkin method, which has many attractive properties, including high-order accuracy on a compact stencil. We present the physical model and numerical method, and show results from a multi-dimensional implementation. Tests show that the method is high-order accurate and strictly preserves realizability of the moments.
Explicit solutions of the radiative transport equation in the P{sub 3} approximation
Liemert, André Kienle, Alwin
2014-11-01
Purpose: Explicit solutions of the monoenergetic radiative transport equation in the P{sub 3} approximation have been derived which can be evaluated with nearly the same computational effort as needed for solving the standard diffusion equation (DE). In detail, the authors considered the important case of a semi-infinite medium which is illuminated by a collimated beam of light. Methods: A combination of the classic spherical harmonics method and the recently developed method of rotated reference frames is used for solving the P{sub 3} equations in closed form. Results: The derived solutions are illustrated and compared to exact solutions of the radiative transport equation obtained via the Monte Carlo (MC) method as well as with other approximated analytical solutions. It is shown that for the considered cases which are relevant for biomedical optics applications, the P{sub 3} approximation is close to the exact solution of the radiative transport equation. Conclusions: The authors derived exact analytical solutions of the P{sub 3} equations under consideration of boundary conditions for defining a semi-infinite medium. The good agreement to Monte Carlo simulations in the investigated domains, for example, in the steady-state and time domains, as well as the short evaluation time needed suggests that the derived equations can replace the often applied solutions of the diffusion equation for the homogeneous semi-infinite medium.