SEARCH · Search NASA

Results for “Particle transport”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Cosmic Ray Propagation through the Magnetic Fields of the Galaxy with Extended Halo

In this project we perform theoretical studies of 3-dimensional cosmic ray propagation in magnetic field configurations of the Galaxy with an extended halo. We employ our newly developed Markov stochastic process methods to solve the diffusive cosmic ray transport equation. We seek to understand observations of cosmic ray spectra, composition under the constraints of the observations of diffuse gamma ray and radio emission from the Galaxy. The model parameters are directly are related to properties of our Galaxy, such as the size of the Galactic halo, particle transport in Galactic magnetic fields, distribution of interstellar gas, primary cosmic ray source distribution and their confinement in the Galaxy. The core of this investigation is the development of software for cosmic ray propagation models with the Markov stochastic process approach. Values of important model parameters for the halo diffusion model are examined in comparison with observations of cosmic ray spectra, composition and the diffuse gamma-ray background. This report summarizes our achievement in the grant period at the Florida Institute of Technology. Work at the co-investigator's institution, the University of New Hampshire, under a companion grant, will be covered in detail by a separate report.

Zhang, Ming↗

NASA NTRS1994

An outer gap model of high-energy emission from rotation-powered pulsars

We describe a refined calculation of high-energy emission from rotation-powered pulsars based on the outer gap model of Cheng, Ho, & Ruderman (1986 a, b). In this calculation, vacuum gaps form in regions near the speed-of-light cylinder of the pulsar magnetosphere along the boundary between the closed and open field line zones. We have improved upon previous efforts to model the spectra from these pulsars (e.g., Cheng et al. 1986b; Ho 1989) by following the variation in particle production and radiation properties with position in the outer gap. Curvature, synchotron, and inverse-Compton scattering fluxes vary significantly over the gap, and their interactions via photon-photon pair production build up the radiating charge populations at varying rates. We have also incorporated an approximate treatment of the transport particle and photon fluxes between gap emission zones. These effects, along with improved computations of the particle and photon distribution, provide very important modifications of the model gamma-ray flux. In paticular, we attempt to make specific predictions of pulse profile shapes and spectral variations as a function of pulse phase and suggest further extensions to the model which may provide accurate computations of the observed high-energy emissions.

Chiang, James↗

NASA NTRS1982

The acceleration of charged particles in interplanetary shock waves

Consideration of the theoretical and observational literature on energetic ion acceleration in interplanetary shock waves is the basis for the present discussion of the shock acceleration of the solar wind plasma and particle transport effects. It is suggested that ISEE data be used to construct data sets for shock events that extend continuously from solar wind to galactic cosmic ray energies, including data for electrons, protons, alphas and ions with Z values greater than 2.0, and that the temporal and spatial evolution of two- and three-dimensional particle distribution functions be studied by means of two or more spacecraft.

Pesses, M. E.↗

NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

A dynamical guiding-center simulation model is used to study the stormtime ion transport which leads to the formation of the ring current and diffusion in the radiation belts. Representative ions guiding-center motion in response to model storm-associated impulses in the convection electric field is traced for a range of ion mu values. The present numerical results are compared with previously formulated limiting idealization of particle transport in order to assess the limits of validity of these approximations. For ions having drift periods that exceed the duration of the main phase of the storm, their inward transport to form the stormtime ring current is appropriately described as direct convective access. For ions having drift periods comparable to the duration of the main phase of the storm, there is a transition between direct convective access and transport that resembles radial diffusion. Lower-energy ring-current ions at L of about 3 are freshly injected there from open adiabatic trajectories, whereas the higher-energy ring-current population consists of a mixture of freshly injected and previously trapped ions.

Chen, Margaret W.↗

NASA NTRS2022

Probing the High-energy Gamma-Ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light-curve Modeling

Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV γ-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to ∼100 GeV) and the Fermi Large Area Telescope to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current-sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second γ-ray light-curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the E field, as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

Monica Barnard↗

NASA NTRS1993

Ion acceleration in impulsive solar flares

Nonrelativistic spectra of protons and ions accelerated in impulsive solar flares are derived using more realistic turbulence power spectra. The calculation is based on a particle transport equation extracted from a second step acceleration model containing stochastic acceleration. The turbulence model is generalized to waves with a small angle to the magnetic field vector and to turbulence power spectra with spectral indices s smaller than 2. Due to the occurrence of impulsive flares at low coronal heights, Coulomb losses at the dense coronal plasma and diffusive particle escape are taken into account. The ion spectra show deviations from long-duration spectra near the Coulomb barrier, where the losses become maximal. The Z-squared/A-dependence of the Coulomb losses leads to spectral variations for different ions. We present a method to estimate the turbulence parameters and injection conditions of the flare particles using ion ratios like Fe/O of impulsive flares.

Steinacker, Jurgen↗

NASA NTRS2017

Transverse eV Ion Heating by Random Electric Field Fluctuations in the Plasmasphere

Charged particle acceleration in the Earth inner magnetosphere is believed to be mainly due to the local resonant wave-particle interaction or particle transport processes. However, the Van Allen Probes have recently provided interesting evidence of a relatively slow transverse heating of eV ions at distances about 2-3 Earth radii during quiet times. Waves that are able to resonantly interact with such very cold ions are generally rare in this region of space, called the plasmasphere. Thus, non-resonant wave-particle interactions are expected to play an important role in the observed ion heating. We demonstrate that stochastic heating by random transverse electric field fluctuations of whistler (and possibly electromagnetic ion cyclotron) waves could explain this weak and slow transverse heating of H+ and O+ ions in the inner magnetosphere. The essential element of the proposed model of ion heating is the presence of trains of random whistler (hiss) wave packets, with significant amplitude modulations produced by strong wave damping, rapid wave growth, or a superposition of wave packets of different frequencies, phases, and amplitudes. Such characteristics correspond to measured characteristics of hiss waves in this region. Using test particle simulations with typical wave and plasma parameters, we demonstrate that the corresponding stochastic transverse ion heating reaches 0.07-0.2 eV/h for protons and 0.007-0.015 eV/h for O+ ions. This global temperature increase of the Maxwellian ion population from an initial Ti approx. 0.3 eV could potentially explain the observations.

Artemyev, A. V.↗

NASA NTRS2023

On the Onset Delays of Solar Energetic Electrons and Protons: Evidence for a Common Accelerator

The processes responsible for the acceleration of solar energetic particles (SEPs) are still not well understood, including whether SEP electrons and protons are accelerated by common or separate processes. Using a numerical particle transport model that includes both pitch-angle and perpendicular spatial diffusion, we simulate, among other quantities, the onset delay for MeV electrons and protons and compare the results to observations of SEPs from widely separated spacecraft. Such observations have previously been interpreted, in a simple scenario assuming no perpendicular diffusion, as evidence for different electron and proton sources. We show that, by assuming a common particle source together with perpendicular diffusion, we are able to simultaneously reproduce the onset delays for both electrons and protons. We argue that this points toward a common accelerator for these particles. Moreover, a relatively broad particle source is required in the model to correctly describe the observations. This is suggestive of diffusive shock acceleration occurring at large shock structures playing a significant role in the acceleration of these SEPs.

R D Strauss↗

NASA NTRS2024

Double Power-law Formation by Sequential Particle Acceleration

Spectral double power laws are common in solar high-energy phenomena such as flares and interplanetary energetic-electron events. However, the physical mechanism that produces the changes in power-law index within a single spectrum is unclear. We developed a fully analytical method of forming single power-law spectra from sequential acceleration of particles orbiting inside and hopping between simulated large-scale magnetic islands formed by flare reconnection. Here, we extend the analytical method to the formation of double power-law spectra by assuming sequential acceleration in two successive regions with different acceleration and particle-transport rates. The resulting spectral distribution is continuous and smooth, with a flattening at low energies, two power-law regions at mid-energies, and a steep rollover at high energies. The model provides analytical expressions for the spectral indices, all energy breaks, and normalization constants as functions of just three physical parameters of each acceleration region: (1) the energy gain in each accelerator, (2) the percentage of particles transferred between accelerators, and (3) the number of accelerators visited. One of the most salient predictions of our work is that the spectral index at high (low) energies is determined by the parameters of the first "seed" (second) acceleration region. By constructing the spectral distribution through an iterative analytical process, the evolution toward a double power law is easily characterized and explained. Our analytical model provides tools to interpret space- and ground-based observations from RHESSI, FOXSI, NuSTAR, Solar Orbiter/STIX, EOVSA, and future high-energy missions.

Solar flares↗

NASA NTRS2005

Modeling Solar Flare Hard X-ray Images and Spectra Observed with RHESSI

Observations obtained with the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of a flare on February 20, 2002 indicate a hard X-ray (HXR) coronal source at or near the top of a flare loop (called a HXR looptop source). The existence of the HXR looptop source suggests that magnetic reconnection, which is believed to power flares, occurs above the loop. In order to explain this HXR looptop source, I created a steady-state particle transport model, in which high-energy electrons are continuously injected at the top of a semicircular flare loop. Based on the simulation results, I find that the model predictions are consistent with the RHESSI observations in many respects, but the spectrum of the looptop source obtained from the model is steeper than that from the RHESSI data. This suggests that, instead of being accelerated above the loop as generally believed, the particles might be accelerated in the looptop itself. RHESSI-observations of three other homologous flares that occurred between April 14 and 16, 2002, provide strong evidence for the presence of a large scale current sheet above a flare loop, which is the basis of standard flare models. The most convincing finding is the presence of the temperature distribution of a separate coronal source above the flare loops: the hotter part of the coronal source was located lower in altitude the cooler part. Together with the fact that the hotter flare loops are higher than the cooler loops, the observations support the existence of a large-scale current sheet between the top of the flare loops and the coronal source above. Blob-Like sources along a line above the loop in the decay phase of the April 15, 2002, flare, which are suggestive of magnetic islands initiated by the tearing-mode instability, and the observation of a cusp structure in microwaves, further support the presence of the current sheet. The observations of the three homologous flares reveal two other features which are beyond the predictions of the standard flare models: the downward motion of flare loops in the early impulsive phase of each flare, and an initially stationary coronal source above the loops. These features me believed to be related to the formation and development of a current sheet. In particular, the downward loop motion seem to be a common phenomenon in flares, suggesting the necessity for modifications to the existing standard flare. models. Finally, thanks to the broad energy coverage of the RHESSI spectra, a low-energy cutoff of 28 (plus or minus 2) keV in the nonthermal electron distribution was determined for the April 15, 2002, flare. As a result, the energy carried by the nonthermal electrons is found to be comparable to the thermal energy of the flare, but one order of magnitude larger than the kinetic energy of the associated coronal mass ejection. The method used to deduce the electron low-energy cutoff will be useful in the analyses of similar events.

Sui, Linhui↗

NASA NTRS2009

Coupled Neutron Transport for HZETRN

Exposure estimates inside space vehicles, surface habitats, and high altitude aircrafts exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETC-HEDS, FLUKA, and MCNPX, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light particle transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Slaba, Tony C.↗

NASA NTRS1999

Laboratory Simulations of Martian and Venusian Aeolian Processes

With the flyby of the Neptune system by Voyager, the preliminary exploration of the Solar System was accomplished. Data have been returned for all major planets and satellites except the Pluto system. Results show that the surfaces of terrestrial planets and satellites have been subjected to a wide variety of geological processes. On solid- surface planetary objects having an atmosphere, aeolian processes are important in modifying their surfaces through the redistribution of fine-grained material by the wind. Bedrock may be eroded to produce particles and the particles transported by wind for deposition in other areas. This process operates on Earth today and is evident throughout the geological record. Aeolian processes also occur on Mars, Venus, and possibly Titan and Triton, both of which are outer planet satellites that have atmospheres. Mariner 9 and Viking results show abundant wind-related landforms on Mars, including dune fields and yardangs (wind-eroded hills). On Venus, measurements made by the Soviet Venera and Vega spacecraft and extrapolations from the Pioneer Venus atmospheric probes show that surface winds are capable of transporting particulate materials and suggest that aeolian processes may operate on that planet as well. Magellan radar images of Venus show abundant wind streaks in some areas, as well as dune fields and a zone of possible yardangs. The study of planetary aeolian processes must take into account diverse environments, from the cold, low-density atmosphere of Mars to the extremely hot, high- density Venusian atmosphere. Factors such as threshold wind speeds (minimum wind velocity needed to move particles), rates of erosion and deposition, trajectories of windblown particles, and aeolian flow fields over various landforms are all important aspects of the problem. In addition, study of aeolian terrains on Earth using data analogous to planetary data-collection systems is critical to the interpretation of spacecraft information and places constraints on results from numerical models and laboratory simulations.

Greeley, Ronald↗

NASA NTRS2023

Robotic and Crewed Mars Missions Increasing the Demand for Planetary Protection Technology Needs

Planetary protection (PP) policy seeks to avoid harmful contamination by limiting biological and relevant organic contamination from spacecraft as well as preventing adverse changes to Earth’s biosphere when extraterrestial samples are brought back to Earth. The PP policy at NASA was updated in 2021 (NPR 8715.24) and 2022 (NASA-STD-8719.27) to enable missions by expanding the decades old prescriptive requirements to allow for an option of adopting performance-based requirements that are objectivesdriven, risk-informed and case-assured. In parallel, the final PP knowledge gap workshop was completed representing the international consensus on the key areas to be considered in developing crew PP policy. These knowledge gaps focused on key technology development areas in 1) microbial and human health monitoring, 2) technical and operations needed for contamination control and 3) natural transport of contamination on Mars. As robotic missions start to implement performance-based approaches and research and technology efforts commence to inform crew policy the demand for data quality driven verification and validation in relevant space environments. Examples of the types of testing that is envisioned includes test as you fly validation and verification of decontamination systems in a relevant on-orbit and Mars environment, developing lethality curves of terrestrial organisms to further our understanding of the biocidal impacts of Mars and the space environment, and particle transport model validation and verification. Thus, the PP discipline has identified the need for groundbased space environments to perform preliminary testing as validation and verification of flight systems and to advance the technology readiness level prior to further testing on-orbit or lunar environments to prepare for Mars.

J. Nick Benardini↗

NASA NTRSDate not supplied

Robotic and Crewed Mars Missions Increasing the Demand for Planetary Protection Technology Needs

Planetary protection (PP) policy seeks to avoid harmful contamination by limiting biological and relevant organic contamination from spacecraft as well as preventing adverse changes to Earth’s biosphere when extraterrestrial samples are brought back to Earth. The PP policy at NASA was updated in 2021 (NPR 8715.24) and 2022 (NASA-STD-8719.27) to enable missions by expanding the decades old prescriptive requirements to allow for an option of adopting performance-based requirements that are objectives-driven, risk-informed and case-assured. In parallel, the final PP knowledge gap workshop was completed representing the international consensus on the key areas to be considered in developing crew PP policy. These knowledge gaps focused on key technology development areas in 1) microbial and human health monitoring, 2) technical and operations needed for contamination control and 3) natural transport of contamination on Mars. As robotic missions start to implement performance-based approaches and research and technology efforts commence to inform crew policy the demand for data quality driven verification and validation in relevant space environments. Examples of the types of testing that is envisioned includes test as you fly validation and verification of decontamination systems in a relevant on-orbit and Mars environment, developing lethality curves of terrestrial organisms to further our understanding of the biocidal impacts of Mars and the space environment, and particle transport model validation and verification. Thus, the PP discipline has identified the need for ground-based space environments to perform preliminary testing as validation and verification of flight systems and to advance the technology readiness level prior to further testing on-orbit or lunar environments to prepare for Mars.

J Nick Benardini↗

NASA NTRS1996

Investigation of Natural and Man-Made Radiation Effects on Crews on Long Duration Space Missions

Over the past several years, NASA has studied a variety of mission scenarios designed to establish a permanent human presence on the surface of Mars. Nuclear electric propulsion (NEP) is one of the possible elements in this program. During the initial stages of vehicle design work, careful consideration must be given to not only the shielding requirements of natural space radiation, but to the shielding and configuration requirements of the on-board reactors. In this work, the radiation transport code MCNP has been used to make initial estimates of crew exposures to reactor radiation fields for a specific manned NEP vehicle design. In this design, three 25 MW(sub th), scaled SP-100-class reactors are shielded by three identical shields. Each shield has layers of beryllium, tungsten, and lithium hydride between the reactor and the crew compartment. Separate calculations are made of both the exiting neutron and gamma fluxes from the reactors during beginning-of-life, full-power operation. This data is then used as the source terms for particle transport in MCNP. The total gamma and neutron fluxes exiting the reactor shields are recorded and separate transport calculations are then performed for a 10 g/sq cm crew compartment aluminum thickness. Estimates of crew exposures have been assessed for various thicknesses of the shield tungsten and lithium hydride layers. A minimal tungsten thickness of 20 cm is required to shield the reactor photons below the 0.05 Sv/y man-made radiation limit. In addition to a 20-cm thick tungsten layer, a 40-cm thick lithium hydride layer is required to shield the reactor neutrons below the annual limit. If the tungsten layer is 30-cm thick, the lithium hydride layer should be at least 30-cm thick. These estimates do not take into account the photons generated by neutron interactions inside the shield because the MCNP neutron cross sections did not allow reliable estimates of photon production in these materials. These results, along with natural space radiation shielding estimates calculated by NASA Langley Research Center, have been used to provide preliminary input data into a new Macintosh-based software tool. A skeletal version of this tool being developed will allow rapid radiation exposure and risk analyses to be performed on a variety of Lunar and Mars missions utilizing nuclear-powered vehicles.

Bolch, Wesley E.↗

NASA NTRS2009

Sources of Sodium in the Lunar Exosphere: Modeling Using Ground-Based Observations of Sodium Emission and Spacecraft Data of the Plasma

Observations of the equatorial lunar sodium emission are examined to quantify the effect of precipitating ions on source rates for the Moon's exospheric volatile species. Using a model of exospheric sodium transport under lunar gravity forces, the measured emission intensity is normalized to a constant lunar phase angle to minimize the effect of different viewing geometries. Daily averages of the solar Lyman alpha flux and ion flux are used as the input variables for photon-stimulated desorption (PSD) and ion sputtering, respectively, while impact vaporization due to the micrometeoritic influx is assumed constant. Additionally, a proxy term proportional to both the Lyman alpha and to the ion flux is introduced to assess the importance of ion-enhanced diffusion and/or chemical sputtering. The combination of particle transport and constrained regression models demonstrates that, assuming sputtering yields that are typical of protons incident on lunar soils, the primary effect of ion impact on the surface of the Moon is not direct sputtering but rather an enhancement of the PSD efficiency. It is inferred that the ion-induced effects must double the PSD efficiency for flux typical of the solar wind at 1 AU. The enhancement in relative efficiency of PSD due to the bombardment of the lunar surface by the plasma sheet ions during passages through the Earth's magnetotail is shown to be approximately two times higher than when it is due to solar wind ions. This leads to the conclusion that the priming of the surface is more efficiently carried out by the energetic plasma sheet ions.

Sarantos, Menelaos↗

NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

NASA NTRS2008

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗