Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Equation of State”

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 73 records · Page 4

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star’s equatorial circumferential radius Re and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSRJ0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are Re = 13.02+1.24-1.06 km andM = 1.44+0.15-0.14 Mmoon (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of Re and M that are consistent with those found in this work. We show that our measurements of Re and M for PSRJ0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

M C Miller↗

Dynamic Nuclear Polarization in Samarium Doped Lanthanum Magnesium Nitrate

The dynamic nuclear polarization of hydrogen nuclei by the solid effect in single crystals of samarium doped lanthanum magnesium nitrate (Sm:LMN) was studied theoretically and experimentally. The equations of evolution governing the dynamic nuclear polarization by the solid effect were derived in detail using the spin temperature theory and the complete expression for the steady state enhancement of the nuclear polarization was calculated. Experimental enhancements of the proton polarization were obtained for eight crystals at 9.2 GHz and liquid helium temperatures. The samarium concentration ranged from 0.1 percent to 1.1 percent as determined by X-ray fluorescence. A peak enhancement of 181 was measured for a 1.1 percent Sm:LMN crystal at 3.0 K. The maximum enhancements extrapolated with the theory using the experimental data for peak enhancement versus microwave power and correcting for leakage, agree with the ideal enhancement (240 in this experiment) within experimental error for three of the crystals.

Byvik, C. E.↗

Simplified model of a volumetric direct nuclear pumped He-3-Ar laser

A physical model for the volumetric nuclear pumped He-3-Ar laser is described and the source terms arising from the He-3(n, p)H-3 reaction as well as the rate equations are given. Steady state solutions are derived for the densities of the species vs the fraction of argon. The dominant pumping mechanism is found to be collisional radiative recombination of the argon atomic ions and subsequent cascading into the upper laser level. The reactions following recombination are described and an expression is derived for the difference in densities of the argon 3d and 4p levels. Comparison with experiments shows good agreement.

Harries, W. L.↗

Elliptic integral solutions to a class of space flight optimization problems

This paper is initially concerned with the minimum-time, exoatmospheric flight of a rocket with constant thrust acceleration magnitude, as in the cases of nuclear and solar electric propulsion. Gravitational acceleration is assumed to be a constant scalar multiple of the radius vector, plus a correction term which is a given function of time. The solution to the state equations is obtained in terms of elliptic integrals. A method is presented for the solution of the two-point boundary-condition problem associated with orbital transfer. At most, the latter method requires iteration upon final time, angle of injection, and two other parameters which are bounded. An example problem is provided which involves a rocket with very low thrust and a spiraling trajectory of many revolutions, but an altitude change of only several hundred miles above the earth. Finally, the original elliptic integral solution is extended to a larger class of low and intermediate thrust problems with constant thrust magnitude, mass decreasing with time, and an inverse square gravitational force.

Andrus, J. F.↗

Gamma-ray lines from neutron stars as probes of fundamental physics

The detection of gamma-ray lines produced at the surface of neutron stars will serve to test both the strong and gravitational interactions under conditions unavailable in terrestrial laboratories. Observation of a single redshifted gamma-ray line, combined with an estimate of the mass of the star will serve as a strong constraint on allowable equations of state of matter at supernuclear densities. Detection of two redshifted lines arising from different physical processes at the neutron star surface can provide a test of the strong principle of equivalence. Expected fluxes of nuclear gamma-ray lines from accreting neutron stars were calculated, including threshold, radiative transfer and redshift effects. The most promising probes of neutron star structure are the deuterium formation line and the positron annihilation line. Detection of sharp redshifted gamma-ray lines from X-ray sources such as Cyg X-1 would argue strongly in favor of a neutron star rather than black hole identification for the object.

Brecher, K.↗

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗

Surface composition of magnetic neutron stars.

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.

Rosen, L. C.↗

Study of the extra-ionic electron distributions in semi-metallic structures by nuclear quadrupole resonance techniques

A straightforward self-consistent method was developed to estimate solid state electrostatic potentials, fields and field gradients in ionic solids. The method is a direct practical application of basic electrostatics to solid state and also helps in the understanding of the principles of crystal structure. The necessary mathematical equations, derived from first principles, were presented and the systematic computational procedure developed to arrive at the solid state electrostatic field gradients values was given.

Murty, A. N.↗

Implementation of variable time step stochastic dynamics for electronically inelastic gas-surface collisions

A variable time step algorithm has been implemented for solving the stochastic equations of motion for gas-surface collisions. It has been tested for a simple model of electronically inelastic collisions with an insulator surface in which the phonon manifold acts as a heat bath and electronic states are localized. In addition to reproducing the accurate nuclear dynamics of the surface atoms, numerical calculations have shown the algorithm to yield accurate ensemble averages of physical observables such as electronic transition probabilities and total energy loss of the gas atom to the surface. This new algorithm offers a gain in efficieny of up to an order of magnitude compared to fixed time step integration.

Garrett, Bruce C.↗

Rigid rotators

The two-particle, steady-state Schroedinger equation is transformed to center of mass and internuclear distance vector coordinates, leading to the free particle wave equation for the kinetic energy motion of the molecule and a decoupled wave equation for a single particle of reduced mass moving in a spherical potential field. The latter describes the vibrational and rotational energy modes of the diatomic molecule. For fixed internuclear distance, this becomes the equation of rigid rotator motion. The classical partition function for the rotator is derived and compared with the quantum expression. Molecular symmetry effects are developed from the generalized Pauli principle that the steady-state wave function of any system of fundamental particles must be antisymmetric. Nuclear spin and spin quantum functions are introduced and ortho- and para-states of rotators, along with their degeneracies, are defined. Effects of nuclear spin on entropy are deduced. Next, rigid polyatomic rotators are considered and the partition function for this case is derived. The patterns of rotational energy levels for nonlinear molecules are discussed for the spherical symmetric top, for the prolate symmetric top, for the oblate symmetric top, and for the asymmetric top. Finally, the equilibrium energy and specific heat of rigid rotators are derived.

Source record↗

An Analog Computer Study of the Low-Frequency Flow Dynamics of Two Nuclear-Rocket Cold-Flow Engine Systems

An analog computer study is made of the low-frequency flow dynamics observed with liquid hydrogen flowing in two different nuclear-rocket cold-flow (unfueled) engine systems. The prediction of the low-frequency response depends largely upon the generation of the proper mass-flow storage in the system. This is achieved by relating, in the model lumping, the average values of fluid enthalpy and a certain thermodynamic fluid property. A two-phase state equation is derived for thermodynamic equilibrium between the gas and liquid phases. Parametric studies are made which demonstrate the effect of the pump characteristics on the flow damping.

Fox, Harry W., Jr.↗

Surface structure of neutron stars with high magnetic fields

The equation of state of cold dense matter in strong magnetic fields is calculated in the Thomas-Fermi and Thomas-Fermi-Dirac approximations. For use in the latter calculation, a new expression is derived for the exchange energy of the uniform electron gas in a strong magnetic field. Detailed calculations of the density profile in the surface region of a neutron star are described for a variety of equations of state, and these show that the surface density profile is strongly affected by the magnetic field, irrespective of whether or not matter in a magnetic field has a condensed state bound with respect to isolated atoms. It is also shown that, as a consequence of the field dependence of the screening potential, magnetic fields can significantly increase nuclear reaction rates.

Fushiki, I.↗

Vibrational Frequencies and Spectroscopic Constants for 1(sup 3)A' HNC and 1(sup 3)A' HOC+ from High-Accuracy Quartic Force Fields

The spectroscopic constants and vibrational frequencies for the 1(sup 3)A' states of HNC, DNC, HOC+, and DOC+ are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schroedinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1(sup 3)A' HCN and HCO+ isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC+ as well as astronomical observations of chemically active astrophysical environments.

Astrochemistry↗

Correlations in alpha-alpha scattering and semi-classical optical models

We show the equivalence of semi-classical solutions to optical model coupled-channel equations derived from Watson's form of the nucleus-nucleus multiple-scattering series to the Glauber multiple-scattering series. A second-order solution to the semi-classical coupled-channel elastic amplitude is shown to be nearly equivalent to a second-order optical-phase-shift approximation to the Glauber amplitude if the densities of all nuclear excited states are approximated by the ground-state density. Using the Jastrow method to model the two-body density we find an average excited-state density to be of negligible importance in the double-scattering region of alpha-alpha scattering.

NASA Discipline Number 04-10↗

Cold equation of state in a strong magnetic field - Effects of inverse beta-decay

The influence of a high magnetic field (B is greater than 10 exp 12 G) on the degenerate matter equation of state appropriate to a neutron star is studied. The regime dominated by relativistic electrons up to the neutron drip density is highlighted. The equilibrium matter composition and equation of state, allowing for inverse beta-decay. Two different equilibrium models are determined: an ideal neutron-proton-electron (npe) gas and the more realistic model of Baym, Pethick, and Sutherland (1971) consisting of a Coulomb lattice of heavy nuclei embedded in an electron gas. For a sufficiently high field strength, the magnetic field has an appreciable effect, changing the adiabatic index of the matter and the nuclear transition densities. The influence of a strong field on some simple nonequilibrium processes, including beta-decay and inverse beta-decay (electron capture) is also considered. The effects produced by the magnetic field are mainly due to the changes in the transverse electron quantum orbits and the allowed electron phase space induced by the field.

Lai, Dong↗

Investigation A: Numerical solutions of driven vortices of binary fluid in a cylindrical geometry

The steady state laminar motion of a viscous, incompressible and binary fluid is studied for a rotating flow in a cylindrical geometry. The mathematical model employed is a cold flow simulation of the fluid mechanics of the light-bulb concept of the gaseous core nuclear engine. A numerical treatment is developed for the rotating flow which includes a description of the nuclear fuel addition. The problem is formulated with the complete Navier-Stokes equations in order to show the interaction between the fuel addition, the main flow, and the boundary layer flow in an accurate manner. The results presented show holdup of the nuclear fuel for the case of steady fuel addition.

Shin, Y. W.↗

Modeling of transient heat pipe operation

The use of heat pipes is being considered as a means of reducing the peak temperature and large thermal gradients at the leading edges of reentry vehicles and hypersonic aircraft and in nuclear reactors. In the basic cooling concept, the heat pipe covers the leading edge, a portion of the lower wing surface, and a portion of the upper wing surface. Aerodynamic heat is mainly absorbed at the leading edge and transported through the heat pipe to the upper and lower wing surface, where it is rejected by thermal radiation and convection. Basic governing equations are written to determine the startup, transient, and steady state performance of a haet pipe which has initially frozen alkali-metal as the working fluid.

Colwell, Gene T.↗

Path Toward a Unifid Geometry for Radiation Transport

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 widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗