Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron source terms”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Asymptotic Giant Branch stars as a source of short-lived radioactive nuclei in the solar nebula

We carried out a theoretical evaluation of the contribution of Asymptotic Giant Branch (AGB) stars to some short-lived (10(exp 6) less than or equal to Tau-bar less than or equal to 2 x 10(exp 7) yr) isotopes in the Interstellar Medium (ISM) and in the early solar system using stellar model calculations for thermally pulsing evolutionary phases of low-mass stars. The yields of s-process nuclei in the convective He-shell for different neutron exposures tau(sub 0) were obtained, and AGB stars were shown to produce several radioactive nuclei (especially Pd-107, Pb-205, Fe-60, Zr-93, Tc-99, Cs-135, and Hf-182) in diferent amounts. Assuming either contamination of the solar nebula from a single AGB star or models for continuous injection and mixing from many stars into the ISM, we calculate the ratios of radioactive to stable nuclei at the epoch of the Sun's formation. The dilution factor between the AGB ejecta and the early solar system matter is obtained by matching the observed Pd-107/Pd-108 and depends on the value of tau(sub 0). It is found that small masses M(sub He) of He-shell material (10(exp -4)-10(exp -7) solar mass) enriched in s-process nuclei are sufficient to contaminate 1 solar mass of the ISM to produce the Pd-107 found in the early solar system. Predictions are made for all of the other radioactive isotopes. The optimal model to explain several observed radioactive species at different states of the proto-solar nebula involves a single AGB star with a low neutron exposure (tau(sub 0) = 0.03 mbarn(sup -1)) which contaminated the cloud with a dilution factor of M(sub He)/solar mass approximately 1.5 x 10(exp -4). This will also contribute newly synthesized stable s-process nuclei in the amount of approximately 10(exp -4) of their abundances already present in the proto-solar cloud. Variations in the degree of homogenization (approximately 30%) of the injected material may account for some of the small general isotopic anomalies found in meteorites. It is also found that Fe-60 is produced in small but significant quantities that may be sufficient to explain the observations if the time elapsed delta from the contamination of the ISM to the formation of protoplanetary bodies is not higher than delta = 5 x 10(exp 6) yr. If delta is longer, up to 10 x 10(exp 6) yr, this would require the single AGB star to experience enhanced neutron densities (n(sub n) approximately 3 x 10(exp 9)n/cu cm) in the s-processing zone in order to compensate for the branching at Fe-59. The alternative model of long-term continuous ejection of matter from many AGB stars does not appear to match the observations. We also estimate the Al-26 production from the H-shell and find that the Al-26 abundance in the early solar system may be readily explained in a self-consistent manner. Moreover, Al-26 from AGB stars may contribute substantially to the galactic Al-26 gamma-source, while no significant gamma-flux from Co-60 (deriving from Fe-60 decay) is to be expected.

Wasserburg, G. J.↗

A non-PRE double-peaked burst with oscillations: burning front propagation and stalling

Non-photospheric-radius-expansion (non-PRE) double-peaked bursts may be explained in terms of spreading (and temporary stalling) of thermonuclear flames from a rotational pole on the neutron star surface, as we argued in a previous study. Here we analyze Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) data of such a burst from the low mass X-ray binary (LMXB) system 4U 1636-536, and show that our model (with ignition at high latitudes) can qualitatively explain the observed burst profile, and spectral evolution. Moreover, the evolution of the source radius inferred from the data shows a strong signature of temporary stalling of the burning front, which is an essential ingredient of our model. This implies that an understanding of thermonuclear flame spreading on neutron stars can be achieved by a simultaneous study of the evolution of intensity and spectrum of these bursts. We also report the discovery of millisecond period brightness oscillations from this burst, which is the first such observation from a non-PRE double-peaked burst. Our model can explain the corresponding oscillation amplitude during the first (weaker) peak, and the absence of oscillations during the second peak. We discuss how observations of oscillations during non-PRE double-peaked bursts provide an additional t 001 for understanding thermonuclear flame spreading successfully.

Bhattacharyya, Sudip↗

Torque Reversal and Spin-Down of the Accretion-Powered Pulsar 4U 1626-67

Over 5 yr of hard X-ray (20-60 keV) monitoring of the 7.66 s accretion-powered pulsar 4U 1626-67 with the Compton Gamma Ray Observatory/BATSE large-area detectors has revealed that the neutron star is now steadily spinning down, in marked contrast to the steady spin-up and spin-down torques differ by only 15% with the neutron star spin changing on a timescale |v/v| approximately equals 5000 yr in both states. The current spin-down rate is itself decreasing on a timescale |v/v| approximately equals 26 yr. The long-term timing history shows small-amplitude variations on a 4000 day timescale, which are probably due to variations in the mass transfer rate. The pulsed 20-60 keV emission from 4U 1626-67 is well-fitted by a power-law spectrum with photon index gamma = 4.9 and a typical pulsed intensity of 1.5 x 10(exp -10) ergs cm (exp -2)s(exp -1). The low count rates with BATSE prohibited us from constraining the reported 42 minute binary orbit, but we can rule out long-period orbits in the range 2 days < or = P(orb) < or = 900 days. We compare the long-term torque behavior of 4U 1626-67 to other disk-fed accreting pulsars and discuss the implications of our results for the various theories of magnetic accretion torques. The abrupt change in the sign of the torque is difficult to reconcile with the extremely smooth spin-down now observed. The strength of the torque noise in 4U 1626-67, approximately 10(exp -22) Hz(exp 2)s(exp -2) Hz(exp -1), is the smallest ever measured for an accreting X-ray pulsar, and it is comparable to the timing noise seen in young radio pulsars. We close by pointing out that the core temperature and external torque (the two parameters potentially relevant to internal sources of timing noise) of an accreting neutron star are also comparable to those of young radio pulsars.

Chakrabarty, Deepto L.↗

LumaCam: a novel class of position-sensitive event mode particle detectors using scintillator screens

A new type of position-sensitive detectors is gaining attention in the neutron community. They are scintillator based detectors that detect the scintillation light on an individual photon basis via an image intensifier and a fast image sensor. Their readout operates in event mode i.e. it produces information about individual neutron interactions, reconstructed from the sensor data, thus enabling to achieve superior spatial and temporal resolutions compared to regular detectors. Although the development of current detectors is focused on neutrons, the concept is also applicable to the detection of other particles such as high-energy photons. This document provides a description on how these detectors are built, how they operate, and what their characteristics are. An example of a detector implementation based on a Timepix3 chip is described to illustrate the detector concept. This includes a detailed description of the algorithm that reconstructs the neutron interactions from the sensor data, one of the core components that sets it apart from established scintillator-based imaging detectors. Energy-resolved epithermal neutron radiography was performed at the ISIS EMMA beamline with this detector, illustrating some of the fundamental differences in the data that can be produced with the new type of detector compared to more established types of scintillator based neutron detectors. The term LumaCam is proposed to refer to this new class of position-sensitive event-mode detectors.

47 OTHER INSTRUMENTATION↗

Position-Dependent Neutron Time-of-Flight Deviation at VULCAN Diffractometer

In neutron diffraction, it is critical to precisely measure the lattice spacing, as it is an indicator of a material’s physical characteristics, such as lattice strain, thermal expansion, and phase structures. In time-of-flight (TOF) measurements, the lattice spacing is determined by the recorded TOF from a well-calibrated instrument. However, changes in neutron time-of-flight are sensitive to many factors, such as the alignment of instrument optics, temperature, sample positions, sample dimensions, internal strains, and chemical or physical heterogeneities at the grain level. At VULCAN (SNS, ORNL), which is a high-flux engineering neutron diffractometer, we used a 1-mm-diameter diamond powder sample to scan for changes in TOF by measuring d-spacing values at different sample positions under several configurations. The 2D map of the TOF deviation or d-spacing deviation, in terms of lattice shift/lattice strain, is reported. The change in TOF is dependent on the scanned location in the beam as well as on detector locations. The results are informative for experimental planning, data interpretation, and future instrument design.

36 MATERIALS SCIENCE↗

Extremum Seeking-Based Control System for Particle Accelerator Beam Loss Minimization

Particle accelerators throughout the world vary widely in terms of age and availability of advanced noninvasive diagnostics that provide varying levels of detail about the accelerated beams. Additionally, beam loss measurements and current monitors are ubiquitous in the accelerator community, they are noninvasive, and they are some of the most important metrics in terms of preventing damage or irradiation of beam pipes and equipment for high-energy machines. However, beam loss and current measurements are difficult to use for feedback tuning because of a lack of a known analytic relationship between scalar loss measurements throughout an accelerator and the hundreds of thousands of parameters that influence the beam and time variation of the beam source itself. In this work, we present a model-independent extremum seeking (ES) controller, which has been implemented for automated tuning and optimization of the Los Alamos Neutron Science Center (LANSCE) linear ion accelerator based only on scalar quantities, such as beam loss or flux measurements. We demonstrate the approach on various accelerator subsystems, including groups of radio frequency (RF) accelerating cavities, quadrupole magnets, bending magnets, and steering magnets. This tool is now available as a use graphical user interface (GUI) in the LANSCE control room, allowing beam physicists to choose arbitrary groups of parameters and beam loss or current monitors for adaptive tuning.

43 PARTICLE ACCELERATORS↗

Transformational challenge reactor design characteristics

The Transformational Challenge Reactor (TCR) program was conceived with the goal to reduce costs and time frames associated with advanced reactor deployment by leveraging developments in advanced manufacturing, advanced materials, data science, and rapid prototyping and testing. The final deliverable of the TCR program was to be an operational test of a novel reactor design. The TCR core design incorporates a dense tri-structural-isotropic/SiC fuel form and volumetrically efficient yttrium hydride moderator, both of which were manufactured and characterized under the TCR program. The TCR is a 3 MW{sub th} He-cooled experimental nuclear reactor designed to reach a total integrated burnup of less than 24 effective full-power hours to keep the radioactive source term to a very low level. TCR design process revealed a positive moderator coefficient; however, the negative doppler coefficients for the fuel and thermal expansion of fuel, moderator, and core support plate yield an overall negative reactivity coefficient. Calculated fuel element temperatures and stresses are well within safety margins. The maximum hypothetical accident (i.e., de-pressurized loss of forced cooling) yields only a modest increase in reactor temperatures that are all within safety margins. This paper summarizes the high-level TCR design characteristics, which were derived from neutronics, thermohydraulics, thermomechanics, and safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The 1979 X-ray outburst of Centaurus X-4

X-ray observations of the first major outburst of the classical transient X-ray source Centaurus X-4 since its discovery in 1969 are presented. The observations were obtained in May, 1979, with the all-sky monitor on board Ariel 5. The flare light curve is shown to exhibit many of the characteristics of other transients, including a double-peaked maximum, as well as significant, apparently random, variations and a lower peak flux and shorter duration than the 1969 event. Application of a standard epoch-folding technique to data corrected for linear decay trends indicates a possible source modulation at 0.3415 days (8.2 hours). Comparison of the results with previous other data on Cen X-4 and the characteristics of the soft X-ray transients allows a total X-ray output of approximately 3 x 10 to the 43rd ergs to be estimated, and reveals the duration and decay time of the 1979 Cen X-4 outburst to be the shortest yet observed from soft X-ray transients. The observations are explained in terms of episodic mass exchange from a late-type dwarf onto a neutron star companion in a relatively close binary system.

Kaluzienski, L. J.↗

Short-lived noble gas effluent trends from a research reactor

An understanding of anthropogenic sources of radioactive noble gases in the atmosphere is needed to enhance the discrimination ability of the International Monitoring System's sensors. These sources include commercial and research nuclear reactors and medical isotope production facilities. While abiding by local environmental ordinances these facilities all emit noble gas radioisotopes through normal operation. Here, this research presents measurements and analysis of noble gas isotopes ( 41 Ar, 135 Xe, 135m Xe, 137 Xe, 138 Xe, 87 Kr, 88 Kr, and 89 Kr) made directly at the stack of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. The Xe and Kr noble gases are concurrently observed with 41 Ar, a neutron activation product, when the reactor is operational. The magnitude of the Xe and Kr noble gases released is not constant over the HFIR cycle, but they temporally match the 41 Ar trend. An isotope activity ratio analysis of these shorter lived isotopes combined with the observation of the cycle's temporal trend helps understand the noble gas production mechanism at the HFIR. Isotopes with short half-lives are not useful for long-range environmental monitoring. However, these measurements could potentially be combined with atmospheric modeling to predict the background source term of the longer-lived Xe ratios at a monitoring station.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Outbursts and Orbit of the Accreting Pulsar GS 1843 - 02 = 2S 1845 - 024

We present observations of a series of 10 outbursts of pulsed hard X-ray flux from the transient 10.6 mHz accreting pulsar GS 1843-02, using the Burst and Transient Source Experiment on the Compton Gamma Ray Observatory. These outbursts occurred regularly every 242 days, coincident with the ephemeris of the periodic transient GRO J1849-03, which has recently been identified with the SAS 3 source 2S 1845-024. Our pulsed detection provides the first clear identification of GS 1843-02 with 2S 1845-024. We present a pulse timing analysis that shows that the 2S 1845-024 outbursts occur near the periastron passage of the neutron star's highly eccentric (e = 0.88 +/- 0.01) 242.18 +/- 0.01 day period binary orbit about a high-mass (M(sub c) greater than 7 solar mass) companion. The orbit and transient outburst pattern strongly suggest that the pulsar is in a binary system with a Be star. Our observations show a long-term spin-up trend, with most of the spin-up occurring during the outbursts. From the measured spin-up rates and inferred luminosities we conclude that an accretion disk is present during the outbursts.

Finger, Mark H.↗

Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection (RR): ▪ Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in RR. ▪ The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research. ▪ RR plays a key role in global models of high-energy sources, in terms of both global-scale layers and reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools.

acceleration of particles↗

Feasibility of Isotopically Tailored Low Activation Materials

Materials employed in future fusion reactors and advanced fission reactors are susceptible to activation and transmutation due to fast neutron flux, which drives the need for Low Activation Materials for survivability, maintenance, and long-term disposal considerations. This paper investigates the feasibility of using ORNL’s planned enrichment technologies to produce isotopically enriched or depleted materials for use in future reactor construction. We focus on identifying isotopes of interest and performing estimates of the required scale of enrichment to impact future design and operation.

07 ISOTOPE AND RADIATION SOURCES↗

Design and Characterization of an Optically Segmented Single Volume Scatter Camera Module

Here, the optically segmented single volume scatter camera (OS-SVSC) aims to image neutron sources for nuclear nonproliferation applications using the kinematic reconstruction of elastic double-scatter events. We report on the design, construction, and calibration of one module of a new prototype. The module includes 16 EJ-204 organic plastic scintillating bars individually wrapped in Teflon tape, each measuring 0.5cx0.5cmx20cm . The scintillator array is coupled to two custom silicon photomultiplier (SiPM) boards consisting of a 2x8 array of SensL J-Series-60035 SiPMs, which are read out by a custom 16 channel DRS4 based digitizer board. The electrical crosstalk between SiPMs within the electronics chain is measured as 0.76%±0.11% among all 16 channels. We report the detector response of one module including interaction position, time, and energy, using two different optical coupling materials: EJ-560 silicone rubber optical coupling pads and EJ-550 optical coupling grease. We present results in terms of the overall mean and standard deviation of the z -position reconstruction and interaction time resolutions for all 16 bars in the module. We observed the 1σz -position resolution for gamma interactions in the 0.3–0.4 MeVee range to be 2.24 cm ± 1.10 cm and 1.45 cm ± 0.19 cm for silicone optical coupling pad and optical grease, respectively. The observed 1σ interaction time resolution is 265 ps ± 29 ps and 235 ps ± 10 ps for silicone optical coupling pad and optical grease, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High Order Implicit Residual-Based Spatial Discretization Error Estimation for S N Neutron Transport

This work demonstrates our novel residual source spatial discretization error estimator (LeR/TEAD) for a DGFEM-1 discretization and assesses it along with two contemporary estimators, Ragusa and Wang's h -refinement estimator (RW) and Duo, Azmy, and Zikatanov's explicit residual-based estimator (DAZ), on a suite of Method of Manufactured Solutions (MMS) 2D problems and three realistic problem geometries. LeR/TE-AD is attractive because it directly estimates the local error in the angular flux, as opposed to a mere indicator of the error's behavior, on the same mesh and method order as the original numerical solution, thus typically being less computationally intensive than a refinement-based method. On the MMS suite, LeR/TE-AD consistently displayed a reduced performance versus its DGFEM-0 results in terms of accuracy and precision metrics, though it was not typically grossly inaccurate. This is attributed to the irregularities in the true solution across singular characteristics limiting the local accuracy of the numerical flux solution, leading to poor derivative approximations used in the residual approximations. The error transport problem then spreads the error in the residual to nearby cells, causing a greater degree of imprecision that did not afflict DAZ or RW. In testing the estimators on realistic problem geometries, however, LeR/TE-AD fared better. In practice, the true error is much larger in non-idealized geometries like in MMS, and a superlinear true solution means that RW and DAZ are not beneficially biased for DGFEM-1 error estimation. LeR/TE-AD was typically first or second in accuracy, primarily competing with RW, but the latter usually consumed 2-4 times the computational time as LeR/TE-AD, and requires a solution with four times as many unknowns. Furthermore, RW and LeR/TE-AD can be used to compute direct estimates of the error in any quantity of interest that is based on the angular ux solution, such as the fission rate density in a fuel pin, whereas DAZ requires a heuristic extension due to its norm-based nature.

97 MATHEMATICS AND COMPUTING↗

Calibration strategy of the PROSPECT-II detector with external and intrinsic sources

This paper presents an energy calibration scheme for an upgraded reactor antineutrino detector for the Precision Reactor Oscillation and Spectrum Experiment (PROSPECT). The PROSPECT collaboration is preparing an upgraded detector, PROSPECT-II (P-II), to advance capabilities for the investigation of fundamental neutrino physics, fission processes and associated reactor neutrino flux, and nuclear security applications. P-II will expand the statistical power of the original PROSPECT (P-I) dataset by at least an order of magnitude. The new design builds upon previous P-I design and focuses on improving the detector robustness and long-term stability to enable multi-year operation at one or more sites. The new design optimizes the fiducial volume by elimination of dead space previously occupied by internal calibration channels, which in turn necessitates the external deployment. In this paper, we describe a calibration strategy for P-II. Here, the expected performance of externally deployed calibration sources is evaluated using P-I data and a well-benchmarked simulation package by varying detector segmentation configurations in the analysis. The proposed external calibration scheme delivers a compatible energy scale model and achieves comparable performance with the inclusion of an additional AmBe neutron source, in comparison to the previous internal arrangement. Most importantly, the estimated uncertainty contribution from the external energy scale calibration model meets the precision requirements of the P-II experiment.

47 OTHER INSTRUMENTATION↗

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Calculations of cosmogenic nuclide production rates in the Earth's atmosphere and their inventories

The production rates of cosmogenic isotopes in the Earth's atmosphere and their resulting terrestrial abundances have been calculated, taking into account both geomagnetic and solar-modulatory effects. The local interstellar flux was assumed to be that of Garcia-Munoz, et al. Solar modulation was accounted for using the heliocentric potential model and expressed in terms of the Deep River neutron monitor count rates. The geomagnetic field was presented by vertical cutoffs calculated by Shea and Smart and the non-vertical cutoffs calculated using ANGRI. The local interstellar particle flux was first modulated using the heliocentric potential field. The modulated cosmic-ray fluxes reaching the earth's orbit then interacted with the geomagnetic field as though it were a high-pass filter. The interaction of the cosmic radiation with the Earth's atmosphere was calculated utilizing the Bolztmann transport equation. Spallation cross sections for isotope production were calculated using the formalism of Silberberg and Tsao and other cross sections were taken from standard sources. Inventories were calculated by accounting from the variation in solar modulation and geomagnetic field strength with time. Results for many isotope, including C-14, Be-7 and Be-10 are in generally good agreement with existing data. The C-14 inventory, for instance, amounts to 1.75/sq cm(e)/s, in excellent agreement with direct estimates.

Obrien, K.↗