Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Modeling”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

A NICER view of PSR J0030+0451: Implications for the Dense Matter Equation of State

Both the mass and radius of the millisecond pulsar PSRJ0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA’s Neutron Star Interior Composition Explorer (NICER) mission. In this Letter we study the implications of the mass–radius inference reported for this source by Riley et al. for the dense matter equation of state (EoS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EoS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EoS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EoS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass–radius likelihood function derived by Riley et al. using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSRJ0030+0451.

G Raaijmakers↗
NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗
NASA NTRS2002

A model for interphase chromosomes and evaluation of radiation-induced aberrations

We have developed a theoretical model for evaluating radiation-induced chromosomal exchanges by explicitly taking into account interphase (G(0)/G(1)) chromosome structure, nuclear organization of chromosomes, the production of double-strand breaks (DSBs), and the subsequent rejoinings in a faithful or unfaithful manner. Each of the 46 chromosomes for human lymphocytes (40 chromosomes for mouse lymphocytes) is modeled as a random polymer inside a spherical volume. The chromosome spheres are packed randomly inside a spherical nucleus with an allowed overlap controlled by a parameter Omega. The rejoining of DSBs is determined by a Monte Carlo procedure using a Gaussian proximity function with an interaction range parameter sigma. Values of Omega and sigma have been found which yield calculated results of interchromosomal aberration frequencies that agree with a wide range of experimental data. Our preferred solution is one with an interaction range of 0.5 microm coupled with a relatively small overlap parameter of 0.675 microm, which more or less confirms previous estimates. We have used our model with these parameter values and with resolution or detectability limits to calculate yields of translocations and dicentrics for human lymphocytes exposed to low-LET radiation that agree with experiments in the dose range 0.09 to 4 Gy. Five different experimental data sets have been compared with the theoretical results. Essentially all of the experimental data fall between theoretical curves corresponding to resolution limits of 1 Mbp and 20 Mbp, which may reflect the fact that different investigators use different limits for sensitivity or detectability. Translocation yields for mouse lymphocytes have also been calculated and are in good agreement with experimental data from 1 cGy to 10 cGy. There is also good agreement with recent data on complex aberrations. Our model is expected to be applicable to both low- and high-LET radiation, and we include a sample prediction of the yield of interchromosomal rejoining in the dose range 0.22 Gy to 2 Gy of 1000 MeV/nucleon iron particles. This dose range corresponds to average particle traversals per nucleus ranging from 1.0 to 9.12.

Non-NASA Center↗
NASA NTRS1980

The chemical abundances of the Cassiopeia A fast-moving knots - Explosive nucleosynthesis on a minicomputer

A simplified nuclear reaction network for explosive nucleosynthesis calculations is described in which only the most abundant nuclear species and the most important reactions linking these species are considered. This scheme permits the exploration of many cases without excessive computational effort. Good agreement with previous calculations employing more complex reaction networks is obtained. This scheme is applied to the observed chemical abundances of the fast-moving knots in the supernova remnant Cassiopeia A and it is found that a wide range of initial conditions could yield the observed abundances. The abundances of four of the knots with significant and different amounts of elements heavier than oxygen are consistent with an origin in material of the same initial composition but processed at different peak temperatures and densities. Despite the observed high oxygen abundances and low abundances of light elements in the knots, they did not necessarily undergo incomplete oxygen burning; in fact, it is not even necessary that oxygen have been present in the initial composition. The agreement between the calculated and observed chemical abundances in Cas A and similar supernova remnants depends primarily upon the relevant nuclear physics and does not provide strong evidence in favor of any particular model of the supernova event.

Johnston, M. D.↗
NASA NTRS1985

Stopping relativistic Xe, Ho, Au and U nuclei in nuclear emulsions

Nuclei of Xe-54, Ho-67, Au-79 and U-92 accelerated at the Bevalac to energies between 1200 and 900 MeV/n were stopped in nuclear emulsions. The observed residual ranges were compared with those calculated from various models of energy loss and shown to be most consistent with a calculation that includes those higher order correction terms proposed previously to describe the energy loss of highly charged particles, for which the first Born approximation is not valid.

Waddington, C. J.↗
NASA NTRS1985

Stopping relativistic Xe, Ho, Au, and U nuclei in nuclear emulsions

Nuclei of Xe-54, Ho-67, Au-79 and U-92 accelerated at the Bevalac to energies between 1200 and 900 MeV/n were stopped in nuclear emulsions. The observed residual ranges were compared with those calculated from various models of energy loss and shown to be most consistent with a calculation that includes those higher order correction terms proposed previously to describe the energy loss of highly changed particles, for which the first Born approximation is not valid.

Waddington, C. J.↗
NASA NTRS1992

Measurements of LET spectra and comparison to models

Measurements of LET spectra for near-earth orbits with various inclinations and altitudes are presented. A comparison with calculated LET spectra shows that the contribution from direct ionizing galactic cosmic rays is well described by the models. An additional contribution to the spectra originates from stopping protons and from nuclear interactions of particles with material. In the case of an interaction, a large amount of energy is deposited in a small volume by target recoils or target fragments. These events will be called short range (SR) events. For a low-inclination orbit radiation belt protons are the main source of these events while galactic protons become more important when increasing the inclination to near-polar orbits. It is shown that the contribution of SR events for orbits with low altitude (324 km) and 57-deg inclination is comparable to that for an orbit with 28-deg inclination at a high altitude (510 km).

Wiegel, B.↗
NASA NTRS1995

HZETRN: Description of a free-space ion and nucleon transport and shielding computer program

The high-charge-and energy (HZE) transport computer program HZETRN is developed to address the problems of free-space radiation transport and shielding. The HZETRN program is intended specifically for the design engineer who is interested in obtaining fast and accurate dosimetric information for the design and construction of space modules and devices. The program is based on a one-dimensional space-marching formulation of the Boltzmann transport equation with a straight-ahead approximation. The effect of the long-range Coulomb force and electron interaction is treated as a continuous slowing-down process. Atomic (electronic) stopping power coefficients with energies above a few A MeV are calculated by using Bethe's theory including Bragg's rule, Ziegler's shell corrections, and effective charge. Nuclear absorption cross sections are obtained from fits to quantum calculations and total cross sections are obtained with a Ramsauer formalism. Nuclear fragmentation cross sections are calculated with a semiempirical abrasion-ablation fragmentation model. The relation of the final computer code to the Boltzmann equation is discussed in the context of simplifying assumptions. A detailed description of the flow of the computer code, input requirements, sample output, and compatibility requirements for non-VAX platforms are provided.

Wilson, John W.↗
NASA NTRS1994

Self-gravity and dissipation in polar rings

Studies of inclined rings inside galaxy potentials have mostly considered the influence of self-gravity and viscous dissipation separately. In this study, we construct models of highly inclined ('polar') rings in an external potential including both self-gravity and dissipation due to a drag force. We do not include pressure forces and thus ignore shock heating that dominates the evolution of gaseous rings inside strongly nonspherical potentials. We adopt an oblate spheroidal scale-free logarithmic potential with axis ratio q = 0.85 and an initial inclination of 80 deg for the self-gravitating rings. We find that stellar (dissipationless) rings suffer from mass loss during their evolution. Mass loss also drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than approximately 0.02 of the enclosed galaxy mass, rings remain bound and do not fall apart from differential precession. The rings precess at a constant rate for more than a precession period tau(sub p) finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self-gravity if the characteristic viscous inflow time scale tau(sub vi) is shorter than approximately 25(tau(sub p)). Rings with (tau(sub vi))/(tau(sub p)) less than or approximately equal to 25 collapse toward the nucleus of the potential within one precession period independent of the amount of self-gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.

Dubinski, John↗
NASA NTRS1963

Performance Capability of Single-Cavity Vortex Gaseous Nuclear Rockets

An analysis was made to determine the maximum powerplant thrust-to-weight ratio possible with a single-cavity vortex gaseous reactor in which all the hydrogen propellant must diffuse through a fuel-rich region. An assumed radial temperature profile was used to represent conduction, convection, and radiation heat-transfer effects. The effect of hydrogen property changes due to dissociation and ionization was taken into account in a hydrodynamic computer program. It is shown that, even for extremely optimistic assumptions of reactor criticality and operating conditions, such a system is limited to reactor thrust-to-weight ratios of about 1.2 x 10(exp -3) for laminar flow. For turbulent flow, the maximum thrust-to-weight ratio is less than 10(exp -3). These low thrusts result from the fact that the hydrogen flow rate is limited by the diffusion process. The performance of a gas-core system with a specific impulse of 3000 seconds and a powerplant thrust-to-weight ratio of 10(exp -2) is shown to be equivalent to that of a 1000-second advanced solid-core system. It is therefore concluded that a single-cavity vortex gaseous reactor in which all the hydrogen must diffuse through the nuclear fuel is a low-thrust device and offers no improvement over a solid-core nuclear-rocket engine. To achieve higher thrust, additional hydrogen flow must be introduced in such a manner that it will by-pass the nuclear fuel. Obviously, such flow must be heated by thermal radiation. An illustrative model of a single-cavity vortex system employing supplementary flow of hydrogen through the core region is briefly examined. Such a system appears capable of thrust-to-weight ratios of approximately 1 to 10. For a high-impulse engine, this capability would be a considerable improvement over solid-core performance. Limits imposed by thermal radiation heat transfer to cavity walls are acknowledged but not evaluated. Alternate vortex concepts that employ many parallel vortices to achieve higher hydrogen flow rates offer the possibility of sufficiently high thrust-to-weight ratios, if they are not limited by short thermal-radiation path lengths.

Ragsdale, Robert G.↗
NASA NTRS2000

Uncertainty Analysis of Decomposing Polyurethane Foam

Sensitivity/uncertainty analyses are necessary to determine where to allocate resources for improved predictions in support of our nation's nuclear safety mission. Yet, sensitivity/uncertainty analyses are not commonly performed on complex combustion models because the calculations are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, a variety of sensitivity/uncertainty analyses were used to determine the uncertainty associated with thermal decomposition of polyurethane foam exposed to high radiative flux boundary conditions. The polyurethane used in this study is a rigid closed-cell foam used as an encapsulant. Related polyurethane binders such as Estane are used in many energetic materials of interest to the JANNAF community. The complex, finite element foam decomposition model used in this study has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state decomposition front velocity calculated as the derivative of the decomposition front location versus time. An analytical mean value sensitivity/uncertainty (MV) analysis was used to determine the standard deviation by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation was essentially determined from a second derivative that was extremely sensitive to numerical noise. To minimize the numerical noise, 50-micrometer element dimensions and approximately 1-msec time steps were required to obtain stable uncertainty results. As an alternative method to determine the uncertainty and sensitivity in the decomposition front velocity, surrogate response surfaces were generated for use with a constrained Latin Hypercube Sampling (LHS) technique. Two surrogate response surfaces were investigated: 1) a linear surrogate response surface (LIN) and 2) a quadratic response surface (QUAD). The LHS techniques do not require derivatives of the response variable and are subsequently relatively insensitive to numerical noise. To compare the LIN and QUAD methods to the MV method, a direct LHS analysis (DLHS) was performed using the full grid and timestep resolved finite element model. The surrogate response models (LIN and QUAD) are shown to give acceptable values of the mean and standard deviation when compared to the fully converged DLHS model.

Hobbs, Michael L.↗
NASA NTRS1989

Fragmentation cross sections of 28Si at 14.5 GeV/nucleon

We used CR39 plastic nuclear track detectors (C12H18O7) in combination with automatic track measurement techniques to determine total charge changing and partial cross sections for the production of fragments of charge ZF = 6 to ZF = 13 in collisions of 28Si beam nuclei at 14.5 GeV/nucleon in targets H, C, CR39, (C18H38O), Al, Cu, Ag, and Pb. By application of factorization rules, measured partial cross sections are separated into pure nuclear and electromagnetic components. The cross sections for electromagnetic dissociation agree with theoretical models. The results are consistent with a Z2 dependence of virtual photon spectra.

NASA Discipline Number 04-10↗
NASA NTRSDate not supplied

Sensitivity of Hybrid NEP-Chemical Vehicle Mass to Assumptions for Crewed Opposition-Class Mars Missions

This research illustrates the impact of different technology choices and their associated performance on overall sizing for a hybrid Nuclear Electric Propulsion (NEP) and chemical crewed Mars transportation vehicle. A set of models for the trajectory and vehicle mass, power system specific mass(α), and power conversion performance enable the analysis. Vehicle mass required to close the mission is calculated as a function of α and specific impulse (Isp) for multiple choices of high-level assumptions including system alpha vs. power scaling and electric thruster type. System alpha is built up and compared for multiple cases of technology assumptions. A Brayton cycle model is used to illustrate the alpha sensitivities to assumptions in cycle design such as peak operating temperature and component performance. The analysis illustrates how technology performance sensitivities affect the alpha and thus vehicle sizing for a range of mission and technology assumptions.

Matthew Duchek↗
NASA NTRS2023

First Detection of X-Ray Polarization from the Accreting Neutron Star 4U 1820−303

This paper reports the first detection of polarization in the X-rays for atoll-source 4U 1820−303, obtained with the Imaging X-ray Polarimetry Explorer (IXPE) at 99.999% confidence level (CL). Simultaneous polarimetric measurements were also performed in the radio with the Australia Telescope Compact Array. The IXPE observations of 4U 1820−303 were coordinated with Swift X-ray Telescope, Neutron Star Interior Composition Explorer, and Nuclear Spectroscopic Telescope Array aiming to obtain an accurate X-ray spectral model covering a broad energy interval. The source shows a significant polarization above 4 keV, with a polarization degree of 2.0% ± 0.5% and a polarization angle of −55° ± 7° in the 4–7 keV energy range, and a polarization degree of 10% ± 2% and a polarization angle of −67° ± 7° in the 7–8 keV energy bin. This polarization also shows a clear energy trend with polarization degree increasing with energy and a hint for a position-angle change of ≃90° at 96% CL around 4 keV. The spectro-polarimetric fit indicates that the accretion disk is polarized orthogonally to the hard spectral component, which is presumably produced in the boundary/spreading layer. We do not detect linear polarization from the radio counterpart, with a 3σ upper limit of 50% at 7.25 GHz.

4U1820-303↗
NASA NTRS2011

PING Gamma Ray and Neutron Measurements of a Meter-Sized Carbonaceous Asteroid Analog

Determining the elemental composition of carbonaceous (spectral type C) asteroids is still one of the basic problems when studying these objects. The only main source of elemental composition information for asteroids is from their optical, NIR and IR properties, which include their spectral reflectance characteristics, albedo, polarization, and the comparison of optical spectroscopy with meteorite groups corresponding to asteroids of every spectral type. Unfortunately, these sources reflect observations from widely contrasting spatial scales that presently yield a void in the continuum of microscopic and macroscopic evidence, a lack of in situ measurement confirmation, and require deeper sensing techniques to discern the nature of these asteroids. The Probing In situ with Neutrons and Gamma rays (PING) instrument is ideally suited to address this problem because it can be used to determine the bulk elemental composition, H and C content, the average atomic weight and density of the surface and subsurface layers of C-type asteroids, and can provide measurements used to determine the difference between and distinguish between different types of asteroids. We are currently developing the PING instrument that combines gamma ray and neutron detectors with a 14 Me V pulsed neutron generator to determine the in-situ bulk elemental abundances and geochemistry of C-type asteroids with a spatial resolution of 1 m down to depths of tens of cm to 1 m. One aspect of the current work includes experimentally testing and optimizing PING on a known meter-sized Columbia River basalt C-type asteroid analog sample that has a similar composition and the same neutron response as that of a C-type asteroid. An important part of this effort focuses on utilizing timing measurements to isolate gamma rays produced by neutron inelastic scattering, neutron capture and delayed activation processes. Separating the gamma ray spectra by nuclear processes results in higher precision and sensitivity elemental composition measurements. Using gated data acquisition techniques allows for the unambiguous identification of gamma ray lines from different isotopes and nuclear processes, especially in situations when limited detector resolution results in overlapping gamma ray lines that cannot be individually resolved. In this paper, we will present the PING basalt layering experimental data, taken at the test facility at NASA Goddard Space Flight Center using the time tagged event-by-event data analysis technique, compared to our MCNPX computer simulation results for the C-type asteroid and basalt layering simulant models. Comparison of these data will show the advantages, validity, and measurement sensitivity of PING's nuclear interrogation methods to obtain more precise and sensitive in situ bulk elemental composition and density measurements of the subsurface of asteroids.

Bodnarik, J.↗
NASA NTRS1989

Biophysical aspects of heavy ion interactions in matter

The biological effects of high energy, high charge nuclei (HZE particles) occupy a central role in the management of space radiation hazards due to galactic cosmic rays. For the energy range of interest, the mean free path for nuclear interactions of these heavy ions is comparable to the thickness of the material traversed, and a significant fraction of stopping particles will undergo a nuclear reaction with the nuclei of the stopping material. Transport methods for HZE particles are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The basic physical properties of HZE particles will be reviewed and illustrated with the results of nuclear fragmentation experiments performed with 670A MeV neon ions incident on a water absorber and with measurements of multiple Coulomb scattering of uranium beams in copper. Finally, the extent to which physical measurements yield radiobiological predictions is illustrated for the example of neon.

Schimmerling, Walter↗
NASA NTRS1996

Model calculations of the radiation dose and LET spectra on LDEF and comparisons with flight data

Ionizing radiation environment models, a 3-D spacecraft mass model, and radiation transport codes have been used to predict the radiation dose and linear energy transfer (LET) spectra measured at various locations on the LDEF satellite. The predictions are compared with thermoluminescent dosimeter measurements of the trapped proton and electron doses and with LET spectra measured by plastic nuclear track detectors. The predicted vs observed comparisons indicate some of the uncertainties of present ionizing radiation environment models for low Earth-orbit missions.

LDEF Project↗