Nuclear thermal rocket testbed restart feasibility study, NLV-019-21, Year 1 of 1
Presentation to be presented during the Site-Directed Research and Development (SDRD) program FY 2021 review meeting (Webex), September 22–23, 2021.
SEARCH · Engineering Papers
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.
Presentation to be presented during the Site-Directed Research and Development (SDRD) program FY 2021 review meeting (Webex), September 22–23, 2021.
Nuclear thermal rockets are currently NASA’s preferred option for use in a manned mission to Mars in the 2040s. The communication delay between an Earth ground station and a spacecraft heading toward Mars can be up to 20 min. Therefore, controlling the nuclear rocket engine would require either a full-time reactor operator on the mission or an autonomous control system for the reactor. The latter idea of making space nuclear reactors fully autonomous has drawn more interest from stakeholders, but such an autonomous control system must be rigorously tested and validated before it is certified for human use. The cost of a full ground test for a space nuclear reactor is tremendous, so a nonnuclear mock reactor test bed was created to test and validate control elements and control algorithms for space nuclear reactors. The test bed consists of control element hardware that inputs physical measurement data into a reactor emulator to produce the reactor’s performance under steady-state, transient, and fault conditions. The control element hardware consists of six full-sized control drums equipped with servo drives and motors and is instrumented with optical encoders, resolvers, and torque sensors for drum movement characterization. In addition to the drums, a two-phase flow loop was designed and built to mimic the valves and turbomachinery associated with the propellant flow through a nuclear thermal rocket engine; components such as pressure sensors, flow meters, thermocouples, and tachometers are instrumented throughout the loop to characterize the fluid flow, valve, and turbomachinery behavior of the system. The data from the physical hardware (e.g., drum position, propellant flow rates) are input to a nuclear reactor simulator to determine the actual nuclear reactor parameters, and the data are sent back to a control algorithm to complete the control loop. The ability to conduct numerous tests of the control systems and autonomous algorithms can help validate the instrumentation and control aspects for a space nuclear reactor for every possible fault situation.
The fundamental approach to nuclear physics was prepared to introduce basic reactor principles to various groups of non-nuclear technical personnel associated with NERVA Test Operations. NERVA Test Operations functions as the field test group for the Nuclear Rocket Engine Program. Nuclear Engine for Rocket Vehicle Application (NERVA) program is the combined efforts of Aerojet-General Corporation as prime contractor, and Westinghouse Astronuclear Laboratory as the major subcontractor, for the assembly and testing of nuclear rocket engines. Development of the NERVA Program is under the direction of the Space Nuclear Propulsion Office, a joint agency of the U. S. Atomic Energy Commission and the National Aeronautics and Space Administration. This report is being reprinted for use in the U. S. Atomic Energy Commission and National Aeronautics and Space Administration educational and technology utilization programs.
Abstract The development and experimentation of a solid propulsion concept consisting of separate hydrogen‐rich solid propellant and solid oxidizer grains is presented. This system consists of an extremely fuel‐rich solid propellant, which, once ignited, undergoes a self‐sustaining decomposition to release fuel‐rich product gases, which are particularly rich in hydrogen gas. These relatively high‐temperature gases subsequently react with a solid oxidizer grain stored downstream. Progress in high‐nitrogen materials has allowed for the synthesis of fuel‐rich compounds that can store significant amounts of hydrogen, with little or no bound oxygen. With these new developments in propellant chemistry, a unique segregated propulsion system has been developed, creating an all‐solid propulsion system that combusts analogously to a reverse hybrid rocket system. Because of the physical separation of fuel and oxidizer, a higher level of safety is achieved which concurrently allows for the utilization of higher energy ingredients without the penalty of higher sensitivity. Theoretical performance calculations and experimental data have shown the potential for this system to compete with existing composite rocket propellants while having significant benefits in terms of safety, toxicity, and mission space. This paper presents a detailed overview of the conception, development, and testing of this propulsion system.
The mixing dynamics of injected propellants is a key factor in determining the ignition performance and combustion-instability response of rocket engines, internal combustion engines, and gas turbines. A key source of uncertainty in the prediction of phase-exchange dynamics is the immiscibility of the injected propellants at supercritical pressure conditions, which induces liquid/vapor phase separation and surface-tension dynamics. While experimental observations indicate the presence of liquid/vapor interfacial structures, this interfacial dynamics is typically neglected in numerical analyses. To address this issue, the objective of the present study is to systematically evaluate the importance of species immiscibility on the phase-exchange dynamics of cryogenic LOX/GH2 mixing layers at typical rocket engine injection conditions. This is accomplished by comparing simulations of (i) a recently developed interface-capturing Regularized-Interface Method (RIM) formulation, and (ii) the commonly employed Diffuse-Interface Method (DIM) formulation. Analysis shows that the interfacial dynamics significantly impact the atomization and mixing of the propellants in the near-injector region, which is not captured by the DIM formulation. The findings of this study extend to other immiscible injection systems, such as LOX/kerosene and hydrocarbon-fuel/air, thereby highlighting the importance of resolving species immiscibility in simulating high-pressure combustion engines.
Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.
The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.
Digital image correlation (DIC) is a popular, noncontacting technique to measure full-field deformation by using cameras to track the motion of an applied surface pattern. Because it is noncontacting, DIC can be performed for extreme temperature applications (e.g., hot-fire rocket testing of carbon composite rocket nozzles) under harsh conditions during which bonded gauges are damaged. Speckle pattern inversion is a phenomenon that sometimes occurs while performing high-temperature DIC. During speckle pattern inversion, portions of the surface pattern that were initially darker at room temperature (e.g., graphite) may emit more light due to blackbody radiation than the portions that were initially paler, thereby producing images in which the pattern appears inverted at high temperature relative to the initial pattern at room temperature. This phenomenon can prevent the correlation algorithm from being able to resolve the displacements between images. This work compares three methods to mitigate speckle pattern inversion: (A) the subtraction method, a recently-published technique in which two high-temperature images are subtracted to remove unwanted light; (B) the filtering method, a popular technique in which optical bandpass filters screen out unwanted light; and (C) the histogram rescaling method, a proposed new method that pairs a color camera with a blue light source and uses information from the green sensor of the camera to correct against inversion in the blue sensor through postprocessing. The histogram rescaling method is shown to successfully eliminate speckle pattern inversion and has the added advantages that it does not require quasi-static loading to be able to compensate for speckle pattern inversion, nor does it impose thick-glass distortions caused by the optical filter.
Here we investigate the application of ground-coupled airwaves observed by seismoacoustic stations at local to near-regional scales to detect signals of interest and determine back-azimuth information. Ground-coupled airwaves are created from incident pressure waves traveling through the atmosphere that couple to the earth and transmit as a seismic wave with retrograde elliptical motion. Previous studies at sub-local scales (<10 km from a source of interest) found the back-azimuth to the source could be accurately determined from seismoacoustic signals recorded by acoustic and 3-component seismic sensors spatially separated on the order of 10 to 150 m. The potential back-azimuth directions are estimated from the coherent signals between the acoustic and vertical seismic data, via a propagation-induced phase shift of the seismoacoustic signal. A unique solution is then informed by the particle motion of the 3-component seismic station, which was previously found to be less accurate than the seismoacoustic-sensor method. We investigate the applicability of this technique to greater source-receiver distances, from 50-100 km and up to 400 km, which contains pressure waves with tropospheric and stratospheric ray paths, respectively. Specifically, we analyze seismoacoustic sources with ground truth from rocket motor fuel elimination events at the Utah Test and Training Range (UTTR) as well as a 2020 rocket launch in Southern California. From these sources we observe evidence that while coherent signals can be seen from both sources on multiple seismoacoustic station pairs, the determined ground-coupled airwave back-azimuths are more complicated than results at more local scales. Our findings suggest more complex factors including incidence angle, coupling location, subsurface material, and atmospheric propagation effects need to be fully investigated before the ground-coupled airwave back-azimuth determination method can be applied or assessed at these further distances.
The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.
This work explores deriving transmissibility functions for a missile from a measured location at the base of the fairing to a desired location within the payload. A pressure on the outside of the fairing and the rocket motor’s excitation creates an acceleration at a measured location and a desired location. Typically, the desired location is not measured. In fact, it is typical that the payload may change, but measured acceleration at the base of the fairing is generally similar to previous test flights. Given this knowledge, it is desired to use a finite-element model to create a transmissibility function which relates acceleration from the previous test flight’s measured location at the base of the fairing to acceleration at a location in the new payload. Four methods are explored for deriving this transmissibility, with the goal of finding an appropriate transmissibility when both the pressure and rocket motor excitation are equally present. In conclusion, these methods are assessed using transient results from a simple example problem, and it is found that one of the methods gives good agreement with the transient results for the full range of loads considered.
Abstract: Boron-loaded propellants have the potential to dramatically increase the performance of solid fuel ramjets, ducted rockets, and hybrid rocket engines. However, difficult ignition of boron decreases the combustion efficiency of these propellants. One approach to solving this problem involves the use of magnesium diboride, MgB2, which ignites easier than boron. Magnesium tetraboride, MgB4, potentially offers greater energetic performance as B has a higher energy density than Mg. However, the effect of the higher boron/metal ratio on the ignition and combustion is unclear. Nanoscale MgB2 particles and quasi 2D structures are promising propellant ingredients, but the oxidation and combustion properties of nanoscale MgB4 remain unknown. Nanoscale magnesium borides are also of interest as precursors for the synthesis of magnesium borohydride, Mg(BH4)2, a promising hydrogen storage material, but hydrogenation of MgB4 has not been studied yet. The objectives of the present work included synthesis, purification, and high-energy ball milling of MgB2 and MgB4 powders as well as investigation of their hydrogen uptake, thermal decomposition, oxidation, and combustion. The powders were fabricated by combustion synthesis and by heating in a tube furnace. The latter method was superior in the synthesis of MgB4. Oxide impurities in the synthesized powders were removed by acid leaching. Nanoscale powders were obtained by ball-mill exfoliation. The hydrogen intake of the obtained magnesium borides was examined at 700 bar and 300 ℃ and did not reveal any advantage of MgB4 over MgB2. Their thermal decomposition and oxidation were investigated with thermogravimetric analysis (TGA), while their combustion was studied using laser ignition and high-speed video recording. TGA has confirmed prior observations of multistep decomposition of magnesium borides, where each step involves formation of a boride with a higher B/Mg ratio and evaporation of formed magnesium. The oxidation rates of the borides are much higher than that of boron at temperatures over 1200 °C for MgB2 and over 900 °C for MgB4. The burning rates of non-milled MgB₂ and MgB₄ powders were much higher than for the used submicron boron. Milling the MgB₂ and MgB₄ powders further increased their burning rates. The milled MgB4 burned 7.5 times faster than submicron boron.
SAND2022-12904 O Initial Atmospheric Transport (IAT) predicts the initial puff/plume rise resulting from either a solid rocket propellant or liquid rocket fuel fire. The code generates initial conditions for subsequent atmospheric transport calculations. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Stoichiometric mixing length L s of reacting coaxial jet flames is a critical scaling parameter for liquid rocket engine combustors. Previous studies have shown that L s for shear coaxial flames can be scaled like their non-reacting counterparts using a non-dimensional momentum flux ratio J. In addition, stoichiometric mixing lengths of reacting and nonreacting coaxial jets collapse upon a single line by altering J using an effective outer flow gas density. This effective density is calculated from a modified version of the equivalence principle, originally developed by Tacina and Dahm [1, 2] and accounts for the effects of heat release on mixing. However, previous studies also required a second nonphysical scaling constant S c for the reacting jets, which is not predicted by the equivalence principle [3]. It was originally hypothesized that S c is attributed to the limitation of hydroxyl (OH) planar laser-induced fluorescence, which only infers L s . Direct quantitative measurement of conserved scalar fields using conventional optical diagnostics is difficult due to the lack of a tracer that easily fluoresces, survives high temperature oxygen flames, and is not dominated by quenching effects. To measure a conserved scalar field, this work implements x-ray fluorescence of Kr and Ar tracers to obtain quantitative mixture fraction fields. From these mixture fraction fields, stoichiometric mixing lengths for two CH 4 /O 2 flames are calculated and scaled against nonreacting coaxial mixing lengths using the equivalence principle. By directly measuring the stoichiometric mixing length, it is established that the additional constant is a byproduct of the OH measurement technique and the equivalence principle fully captures the scaling. Finally, comparison with high-fidelity simulation of the flame further supports this conclusion. In addition to further strengthening this scaling method, this work represents the first use of x-ray fluorescence to make quantitative conserved scalar measurements in turbulent flames.
A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.
The purpose of this letter report is to document the mitigation of adverse effects of a proposed undertaking that would demolish Building 25-3124 in Area 25 of the Nevada National Security Site (NNSS) in compliance with the terms of the 2024 Programmatic Agreement among the U.S. Department of Energy and the Nevada State Historic Preservation Officer and the Advisory Council on Historic Preservation Concerning the Protection of Historic Properties on the Nevada National Security Site, Nye County, Nevada (NNSS PA). Building 25-3124 (State Historic Preservation Office [SHPO] Resource # B19009) is a contributing element to the Nuclear Rocket Development Station (NRDS) Historic District (SHPO Resource # D424), which has been determined eligible for listing in the National Register of Historic Places (NRHP) under Criteria A, B, C, and D (Reno et al. 2023). The proposed demolition constitutes an adverse effect. The U.S. Department of Energy (DOE) National Nuclear Security Administration Nevada Field Office (NNSA/NFO) consulted with the SHPO on an adverse effect finding for the undertaking and notified the SHPO of its intent to use the standard mitigation in the NNSS PA on December 11, 2024 (Julian to Reed). The SHPO concurred on December 30, 2024 (Reed to Julian). Therefore, NNSA/NFO has prepared mitigation pursuant to the standard mitigation in the NNSS PA Appendix D.I.A for contributing elements to recorded, NRHP-eligible districts.
The purpose of this letter report is to document the mitigation of adverse effects of a proposed undertaking that would demolish Building 25-3153 in Area 25 of the Nevada National Security Site (NNSS) in compliance with the terms of the 2024 Programmatic Agreement among the U.S. Department of Energy and the Nevada State Historic Preservation Officer and the Advisory Council on Historic Preservation Concerning the Protection of Historic Properties on the Nevada National Security Site, Nye County, Nevada (NNSS PA). Building 25-3153 (State Historic Preservation Office [SHPO] Resource # B19004) is a contributing element to the Nuclear Rocket Development Station (NRDS) Historic District (SHPO Resource # D424), which has been determined eligible for listing in the National Register of Historic Places (NRHP) under Criteria A, B, C, and D (Reno et al. 2023). The proposed demolition constitutes an adverse effect. The U.S. Department of Energy (DOE) National Nuclear Security Administration Nevada Field Office (NNSA/NFO) consulted with the SHPO on an adverse effect finding for the undertaking and notified the SHPO of its intent to use the standard mitigation in the NNSS PA on December 11, 2024 (Julian to Reed). The SHPO concurred on December 30, 2024 (Reed to Julian). Therefore, NNSA/NFO has prepared mitigation pursuant to the standard mitigation in the NNSS PA Appendix D.I.A for contributing elements to recorded, NRHP-eligible districts.
The purpose of this letter report is to document the mitigation of adverse effects of a proposed undertaking that would demolish Building 25-4314 in Area 25 of the Nevada National Security Site (NNSS) in compliance with the terms of the 2024 Programmatic Agreement among the U.S. Department of Energy and the Nevada State Historic Preservation Officer and the Advisory Council on Historic Preservation Concerning the Protection of Historic Properties on the Nevada National Security Site, Nye County, Nevada (NNSS PA). Building 25-4314 (State Historic Preservation Office [SHPO] Resource # B19016) is a contributing element to the Nuclear Rocket Development Station (NRDS) Historic District (SHPO Resource # D424), which has been determined eligible for listing in the National Register of Historic Places (NRHP) under Criteria A, B, C, and D (Reno et al. 2023). The proposed demolition constitutes an adverse effect. The U.S. Department of Energy (DOE) National Nuclear Security Administration Nevada Field Office (NNSA/NFO) consulted with the SHPO on an adverse effect finding for the undertaking and notified the SHPO of its intent to use the standard mitigation in the NNSS PA on December 11, 2024 (Julian to Reed). The SHPO concurred on December 30, 2024 (Reed to Julian). Therefore, NNSA/NFO has prepared mitigation pursuant to the standard mitigation in the NNSS PA Appendix D.I.A for contributing elements to recorded, NRHP-eligible districts.