SEARCH · Search NASA
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.
Computer program VARI-QUIR 3 provides Gauss-Seidel type of solution with inner and outer iterations for steady-state, multigroup, two-dimensional neutron diffusion equations. The program has no restrictions on any of the input parameters such as the number of groups, regions, or materials.
The extent to which suitable solutions may be obtained for one physics problem and two engineering type problems is traced. NASTRAN appears to be a practical tool to solve one-group steady-state neutron diffusion equations. Transient diffusion analysis may be performed after new levels that allow time-dependent temperature calculations are developed. NASTRAN piecewise linear anlaysis may be applied to solve those plasticity problems for which a smooth stress-strain curve can be used to describe the nonlinear material behavior. The accuracy decreases when sharp transitions in the stress-strain relations are involved. Improved NASTRAN usefulness will be obtained when nonlinear material capabilities are extended to axisymmetric elements and to include provisions for time-dependent material properties and creep analysis. Rigid formats 3 and 5 proved to be very convenient for the buckling and normal-mode analysis of a nuclear fuel element.
A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.
A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.
Three mechanisms that promote the loss of magnetic flux from an isolated neutron star - Ohmic decay, ambipolar diffusion, and Hall drift - are investigated. Equations of motions are solved for charged particles in the presence of a magnetic field and a fixed background of neutrons, while allowing for the creation and destruction of particles by weak interactions. Although these equations apply to normal neutrons and protons, the present interpretations of their solutions are extended to cover cases of neutron superfluidity and proton superconductivity. The equations are manipulated to prove that, in the presence of a magnetic force, the charged particles cannot be simultaneously in magnetostatic equilibrium and chemical equilibrium with the neutrons. The application of the results to real neutron stars is discussed.
The high energy proton radiation belts of Saturn are energetically dominated by the source from cosmic ray albedo neutron decay (CRAND), trapping of protons from beta decay of neutrons emitted from galactic cosmic ray nuclear interactions with the main rings. These belts were originally discovered in wide gaps between the A-ring, Janus/Epimetheus, Mimas, and Enceladus. The narrow F and G rings significant affected the CRAND protons but did not produce total depletion. Voyager 2 measurements subsequently revealed an outermost CRAND proton belt beyond Enceladus. Although the source rate is small, the trapping times limited by radial magnetospheric diffusion are very long, about ten years at peak measured flux inwards of the G ring, so large fluxes can accumulate unless otherwise limited in the trapping region by neutral gas, dust, and ring body interactions. One proposed final extension of the Cassini Orbiter mission would place perikrone in a 3000-km gap between the inner D ring and the upper atmosphere of Saturn. Experience with CRAND in the Earth's inner Van Allen proton belt suggests that a similar innermost belt might be found in this comparably wide region at Saturn. Radial dependence of magnetospheric diffusion, proximity to the ring neutron source, and northward magnetic offset of Saturn's magnetic equator from the ring plane could potentially produce peak fluxes several orders of magnitude higher than previously measured outside the main rings. Even brief passes through such an intense environment of highly penetrating protons would be a significant concern for spacecraft operations and science observations. Actual fluxes are limited by losses in Saturn's exospheric gas and in a dust environment likely comparable to that of the known CRAND proton belts. The first numerical model of this unexplored radiation belt is presented to determine limits on peak magnitude and radial profile of the proton flux distribution.
It has been shown previously (Anath et al., 1973 and Kane, 1974) that 20 to 25% of days, the diffusion component of the cosmic-ray neutron diurnal anisotropy is directed more than 30 degrees away from the ecliptic projection of the interplanetary magnetic field averaged over the same 24 hours. A number of explanations for this deviation are discussed and it is concluded that transverse gradient drifts due to gradients perpendicular to the ecliptic are likely, that diurnal variations in the diffusion component of the neutron anisotropy may affect results from single stations and that the 24 hour mean interplanetary magnetic field may not be the field appropriate to the streaming equation at neutron monitor energies.
Results of the numerical solution of the anistoropic diffusion equation are presented. The modulation depth of galactic cosmic rays is defined by the degree of curvature of the neutral current sheet in the heliosphere. The effect of the regular interplanetary magnetic field (IMF) on cosmic ray anisotropy in the period of solar activity minimum (in 1976) is analyzed by the data of the neutron super-monitors of the world network, and the heliolatitudinal gradient and cosmic ray diffusion coefficient are defined.
This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.
A method is presented for determining the critical size of a cylindrical reactor by a one-dimensional group-diffusion method extended to a two-dimensional solution by prescribing values of axial leakage based on assumed flux levels. The neutron energy spectrum is divided into nine groups, and the reactor into four concentric cylinders. A computing machine program for an IBM 650 computer with attachments and a method for using the program with nine, or less, groups and four, or less, regions are given.
The Diviner lunar radiometer instrument aboard the Lunar Reconnaissance Orbiter mission has revealed large areas of lunar polar terrain with surface temperatures well below 100K. At these temperatures, the sublimation rate of water ice is well below 1 mm per billion years. In contrast, the loss rate at 120K is more than 1 meter of ice in that time consequently volatiles delivered to the coldest locations can be trapped for over 1 Ga, but will be quickly lost from warmer locales. Here we investigate the loss or retention of a layer of ice-bearing regolith at the lunar south poe, assuming contemporary surface temperature conditions and no other loss processes. We use an analytic solution for the one-dimensional diffusion equation of water ice, assuming an isothermal regolith with pore space comparable to mean grain size, 75 micrometers. Only the top meter of soil is assumed to be ice-bearing. We can then calculate the history of ice content with time based on local temperature, and predict what the epithermal neutron output would be in the presence of such a concentration of hydrogen. We compare the present, observed distribution of hydrogen with what one would expect from the temperature-dependent loss or retention of ice for various times since emplacement.
The relative abundances of seven constituent nuclei, He4, C12, O16, Ne20, Mg24, Si28, and Fe56, are calculated as a function of time for neutron star atmospheres within which exist magnetic fields of the order of 10 to the 13th G. The opacity, equation of state of the electrons, and cooling rate of the magnetic star are discussed, and it is shown to be a reasonable approximation to assume an atmosphere to be isothermal. The effects of particle diffusion are included in the nuclear reaction network. Computations are performed both for a constant mass atmosphere and for an atmosphere in which mass is being ejected. It is found that the final abundances are model-independent, as long as the initial model contains predominantly He4. The relative abundances are compared to the cosmic ray spectrum. For both the constant-mass and mass-loss atmospheres, nucleosynthesis proceeds virtually completely to Fe56. However the outermost layers of the envelope, in which no mass is being ejected, are composed almost entirely of He4 with trace amounts of Fe56. After the loss of about 10 to the 21st g, only Fe56 is ejected from atmospheres expelling mass.
A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material show that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux but increases sharply in the presence of a sustainer electric field.
A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material shows that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux.
The experiments and theories to explain the high-energy protons trapped in the earth's radiation belt are reviewed. The theory of cosmic ray albedo neutron decay injection of protons into the radiation belt is discussed. Radial diffusion and change in the earth's dipole moment are considered along with losses of protons by ionization and nuclear collision. It is found that the measured albedo neutron escape current is sufficient to supply trapped protons above 30 MeV. The theoretical calculations of the trapped protons are in agreement with the measurements for L less than or equal to 1.7 both on and off the equator. For L greater than or equal to 1.7, additional trapped proton differential energy measurements should be made before the theory can be adequately tested. It appears that an additional loss mechanism such as pitch angle scattering may be required.
Most plausible scenarios for the formation of planets around pulsars end with a disk of gas around the pulsar. The supplicant author then points to the solar system to bolster faith in the miraculous transfiguration of gas into planets. We here investigate this process of transfiguration. We derive analytic sequences of quasi-static disks which give good approximations to exact solutions of the disk diffusion equation with realistic opacity tables. These allow quick and efficient surveys of parameter space. We discuss the outward transfer of mass in accretion disks and the resulting timescale constraints, the effects of illumination by the central source on the disk and dust within it, and the effects of the widely different elemental compositions of the disks in the various scenarios, and their extensions to globular clusters. We point out where significant uncertainties exist in the appropriate grain opacities, and in the effect of illumination and winds from the neutron star.
On the basis of the model and experimental investigations, the spatial distribution of cosmic ray anisotropy for the different epochs of solar activity is studied. A solution is offered to the anisotropic diffusion equation with regard to the electromagnetic conditions for the periods of minimum and maximum solar activity and cosmic ray particle drift in a regular interplanetary magnetic field. It is shown that the long period changes amplitude and phase of the diurnal variations of cosmic rays is limited not only by convection and diffusion of particles but also by the drift effect before and after the sing change of the general magnetic field of the Sun. The calculated model is compared with the results obtained on the basis of an analysis of the experimental data from the neutron super monitor station, and it is shown that for the periods when the lines of magnetic force of the Sun come from the Northern Hemisphere the phase of the first harmonic diurnal variation is shifted forwards to an earlier time.
Steady state, spherically symmetric solutions of the cosmic-ray transport equation describing the acceleration of energetic particles in galactic accretion flows onto neutron stars, black holes, white dwarfs, and protostars are studied. The results indicate that astrophysical accretion flows can be partitioned into distinct classes depending upon whether the accretion rate lies above or below a critical value of a few times 10 to the -7th stellar masses/yr. When the particle transport is convection-dominated, both classes of accretion flows exhibit a spectral index appropriate for first-order Fermi acceleration at a plane shock in the absence of losses. As the particle transport becomes diffusion-dominated, both classes show a break and precipitous falloff in the particle spectrum due to the escape of these particles from the accretion flow. The precise nature of the spectrum depends on the relationship between the particle momentum and the spatial diffusion coefficient.