Search

An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time�-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.

The isolator pseudo-shock provides necessary compression within a dual-mode scramjet engine and buffers the engine system against unstart. Quasi-1D flux-conserved models are the state-of-the-art reduced-order model for optimization and online control of dual-mode scramjet engines. The stability and efficacy of this modeling approach is evaluated against data from a combustor-driven direct-connect experiment. The experiment exhibited strong combustor-driven unsteadiness that produced upstream propagating weak shocks into the isolator, interacting with the pseudo-shock. While this configuration resulted in unsteadiness that is atypical of standard operation, the experiment provided an opportunity to evaluate the modeling techniques in highly transient states. Such transients could occur during maneuvering or result from unexpected combustor events. A flexible quasi-1D formulation, the Fievet flux-conserved model is fit using Bayesian inference in the laboratory and shock-stationary reference frames. Model performance is analyzed using the Bayesian posteriors and model evaluations over the measured shock train speed range. It is concluded that to produce consistent isolator pressure profile estimates in this unsteady environment, the model must be implemented in a shock-stationary reference frame. Implementing this conclusion in model-based engine controllers may reduce needed unstart safety margins and increase maximum performance.

A digital model of an OxEon Energy methanation reactor (with a cooling jacket) planned for use in a NASA designed In-Situ Resource Utilization (ISRU) propellant production system was developed by the UTEP Aerospace Center using STAR-CCM+ software, a Siemens Computational Fluid Dynamics (CFD) package. The reactor digital model created includes two primary components: (i) a cooling jacket that maintains the required reactor chamber surface temperature and (ii) a reactor core chamber where the exothermic methanation reactions occur. The StarCCM+ digital model has two objectives: (i) maintain the surface temperature of the reactor below a set value so that the extra heat accumulation does not melt or crack the chamber, and (ii) produce the desired amount of methane (CH4) according to test cases completed by OxEon Energy for NASA JSC. The computational digital model was validated with the test data enabling the cooling jacket scale to be predicted based on the reactor's operating temperature. The StarCCM+ simulation of the digital model shows a good agreement with the cooling jacket performance and the CH4 production. One test case for the cooling jacket maintained the reactor wall temperature around 310 °C was validated and the reactor produced 185 g/hr CH4.

We present the results from a direct numerical simulation of a spanwise-periodic turbulent flow past a Gaussian bump. The problem setup is designed to investigate the interaction of an incoming turbulent boundary layer with the strong favorable and adverse pressure gradients generated by the Gaussian bump as the flow passes over it at a Reynolds number of 340000 based on the bump height, or 4 million based on the bump length. The statistical results from the present simulation are compared against our earlier results at a Reynolds number of 2 million. An internal layer, which forms beneath the strongly accelerated boundary layer over the windward side of the bump, is found to generate its near-wall turbulence stress peaks in closer proximity of the wall in the higher Reynolds-number case. Furthermore, the logarithmic layer of the higher Reynolds-number boundary layer appears more resistant to changes induced by strong acceleration and surface curvature effects over the same region. Despite a nearly identical flow separation point in the two flows, the detached shear layer grows at a faster rate and subsequently reattaches at an earlier point in the higher Reynolds-number flow. The surface pressure and skin-friction distributions over the attached flow region compare well against the corresponding experimental data for both flows. However, some differences appear in the separated flow region, which are attributed to the three-dimensionality of the experimental model setup that is not included in the simulation owing to the spanwise periodic assumption. Comparisons with the stereoscopic particle image velocimetry measurements on the central plane of the experimental model over the windward side of the bump show reasonable overall agreement in the mean velocity components, but the turbulence stress components do not agree well at some streamwise locations. Comparisons over the leeward side of the bump show that the mean separated shear layer in the simulation is tilted significantly more toward the wall than the experimental shear layer on the central plane. This mismatch in the mean shear layer orientation is due to the experimental model three-dimensionality and tunnel end-wall effects, which are not modeled in the present spanwise-periodic simulation.

This work develops a Computational Fluid Dynamics (CFD) model of a Commercial-Off-The-Shelf (COTS) condenser used for ground-based testing of an In-Situ Resource Utilization (ISRU) propellant production system. STAR-CCM+ was used to simulate water condensing through the condenser downstream of a methanation reactor. The computational data of the digital model is validated with test cases using nitrogen and water in ambient conditions. After validation with the test data, the digital model can be used to predict scale at higher production targets. The StarCCM+ simulation of the digital model shows good agreement with the test data provided. From the simulation, the condensation rate calculated was 18.3 g/hr and the process gas mass flow rate was 709.2 g/hr. The test data measured an average 16.9 g/hr condensation rate and 693.6 g/hr process gas mass flow rate. The simulation shows a 7.3% error in the water condensation and 2.2% error in the process gas at the outlet.

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Estimation of boundary layer edge or thickness is an essential aspect of utilizing high fidelity aerothermodynamic flow solutions. For viscous flow over a flat plate, the boundary layer thickness is usually defined as the distance from the wall at which the flow velocity (zero at the wall) has increased through some profile to 99% or 99.5% of the free stream velocity. In real gas flows over aerospace vehicles, a total enthalpy profile is employed for edge calculations. The importance of knowing boundary layer thickness and its relation to possible transition from laminar to turbulent flow is illustrated by Space Shuttle mission STS-114, when the appearance of tile gap filler protrusions on the nose was deemed serious enough that a crew member was called upon to exit the vehicle and remove the potential threat prior to its descent from orbit. The BLAYER utility used as part of post-processing Orion/MPCV flow solutions is a generalization of the earlier BLAYER_RESULTS utility developed following the STS-107 Columbia accident. BLAYER handles a variable number of gas species (not just 5), allows for more than one temperature and for optional extra flow quantities, and omits optional handling of Orbiter tile datasets.

A wind tunnel test program was conducted to obtain experimental convective heating data on models of the Low Earth Orbital Flight Test of an Inflatable Decelerator (LOFTID) vehicle in support of the flight test program. Wind tunnel testing was conducted at Mach 6 and Mach 10 at unit Reynolds numbers from 0.5×106/ft to 4.0×106/ft on 6.0inch diameter aeroshell models. Global heating data were obtained through phosphor thermography measurements on the front and back of the aeroshell and on the payload. For all test conditions, attached, laminar flow was definitively produced on the front of the aeroshell, while in the wake the flow was separated, unsteady and possibly transitional or turbulent. A complementary computational study was performed to assess the accuracy of heating predictions for LOFTID wake flow environments. A new grid topology was developed for the study to allow the use of structured grid solvers for complex, open-backed geometries such as that of LOFTID. Close agreement between measurements and predictions was obtained on the front of the aeroshell. Comparisons were more challenging in the wake because of both the extremely low heating rates and the accompanying high uncertainties and the unsteady nature of the separated wake flow simulations. While the differences between heating predictions and measurements for the wake flow were higher than for the forebody flow, these differences were used to assist in the definition of conservative uncertainty margins for predictions of the flight environments for the successful flight test of the LOFTID vehicle in November 2022.

One of the scientific goals of NASA’s Artemis missions is to return regolith samples from the lunar surface. To maintain scientific integrity of these samples, including any volatiles that may exist, cryogenic stowage is desired. FROSTE (Frozen Return Of Samples To Earth) is the element of the Moon to Mars architecture that is tasked with developing and advancing concepts related to this goal. Both active and passive solutions are being considered, but active solutions (cryo-coolers) have high mass and cost penalties. Phase change materials (PCMs) are commonly used as passive thermal management solutions, particularly for sensitive scientific samples. However, there are currently no PCMs capable of providing cooling at the cryogenic temperatures needed to enable lunar sample return with volatiles intact (<120 K). FROSTE has identified two possible candidate materials that could fill this gap. Testing is required to verify the thermal performance of these materials, as well as gain experience with their handling characteristics. We are conducting this testing at the Marshall Space Flight Center’s Space Environmental Effects laboratory, in vacuum, with liquid nitrogen providing cooling ability below the freezing point of the materials. In this talk, we will present an overview of the FROSTE program element and share results from the initial development testing of two cryogenic phase change materials.

Emerging concepts for fission surface power and nuclear electric propulsion necessitate lightweight, mechanically robust, and thermally efficient heat rejection radiators. State-of-the-art intermediate-temperature (~400 K) composite radiator assemblies have been developed based on titanium-water heat pipes bonded to metal, graphite, and carbon-fiber-based panels. NASA has identified a need for new radiator concepts that can operate at even higher temperatures (500 – 600 K), minimize thermal resistances and thermal stress failures at bond interfaces, and approach areal densities of 2 – 3 kg m-2. To meet these needs, our team is developing additively manufactured (AM) radiator panels with integral branching wicking heat pipe networks. Water is selected as the working fluid for this temperature range. Based on simulations and thermal vacuum experiments, these branching embedded heat pipe networks can efficiently distribute heat over panels for finned surface efficiencies of ηf >70% at TH = 500 K input heat. This paper first presents laser powder-bed fusion AM strategies to produce embedded porous structures for wicking heat pipes in Inconel 718 and titanium alloys (commercially pure and Ti-6Al-4V alloys). Post-build chemical and thermal treatments are described that yield hydrophilic wicking surfaces for operation with water. Transient rate-of-rise experiments with water and acetone are reported that yield estimates for AM wick porosity (ϵ), permeability (K), and effective pore radius (rpore). Based on the wick characterization results, small prototype radiator panels (75 × 125 mm) with integrated heat pipe networks were manufactured. Heat rejection performance data are presented from cold thermal vacuum testing, with heat input temperatures up to ~510 K. Future efforts will focus on improving heat pipe performance, optimizing radiator mass, and evaluating larger panels to assess scalability.

Spacecraft nuclear electric propulsion, as well as other emerging technologies, can benefit from heat pipes with optimal performance in the intermediate temperature, 450–700 K, regime. In this work, a 1D model was developed to investigate the performance of halide working fluids coupled with 3D printed-ceramic shells to operate at 500–600 K and interface with a nuclear electric propulsion system. An initial conceptual design and performance limits of a 3D printed-AlN heat pipe radiator with AlBr3 as the working fluid is presented. An areal density of 5.2 kg/m2 and mass per unit heat performance of 1.00kg/kW is estimated.

Heat transfer in high-temperature, high-porosity, flexible refractory ceramic fibrous insulation felts is investigated. Heat transfer in these insulation materials consists of combined gas conduction, solid conduction, and radiation modes, with the precise theoretical modeling of the latter two modes being formidable. A semi-empirical model that requires inverse methods and steady-state thermal test data to infer some of the required model parameters is further developed in this study, and applied to five insulation materials for temperatures between 300 K and 1900 K. The steady-state thermal test setup at NASA Langley Research Center with recent modifications to increase its testing capability to 1900 K is discussed. Design considerations to ensure one-dimensional heat transfer in the test setup are described. Test data and corresponding thermal models for alumina and zirconia-based fibrous insulation felts are presented. Furthermore, test data and thermal models on two fibrous insulation samples containing additives to further suppress either radiation or gas conduction modes of heat transfer are presented. Previously published alumina-based insulation data are also re-processed using the updated modeling methodology. The significance of various heat transfer modes in typical insulation samples is discussed and used to provide general guidance on optimum insulation layups.

The NASA Juncture Flow Experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This report presents and discusses the results of a recent experiment with the Juncture Flow Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of the test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model angle of incidence. Laser Doppler velocimetry (LDV) measurements were made at three model angles of incidence (0° : fully attached, 1° : incipient separation, and 5° : separated flow) and for each angle, mean-flow and Reynolds-stress data were obtained on the fuselage and at several locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and include model and tunnel wall static pressures, tunnel wall boundary-layer profiles, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. A companion set of Reynolds-averaged Navier-Stokes (RANS) CFD simulations were also performed and those results were used to support the discussion of the experimental test results.

GlennICE (Glenn Icing Computational Environment) is a computational tool designed to calculate ice growth on complex three-dimensional geometries using the input from a user-supplied computational fluid dynamics (CFD) solution for the geometry of interest. The NASA John H. Glenn Research Center at Lewis Field is developing this tool to aid those evaluating, designing and certifying aircraft, engines, and aircraft components for flight in icing conditions. This domestically available software is being developed to enable the introduction of new icing physics into a computational environment in a manner that is open for evaluation and eventual use by industry, academia, and other government organizations.

A model for solid rocket motor residual thrust is built from the concept that hot nozzle phenolic or carbon materials provide a radiating source of heat flux to the insulation, leading to continued off-gassing by pyrolysis. This model is 1. Fitted to flight datasets and scaled to the Mars Ascent Vehicle second stage motor (MAV SRM2), and 2. Derived with MAV SRM2 design and material estimates. The two methods show a potential range of modeling uncertainty for use by mission analysts. Further loss methods that could reduce the maximum residual impulse are identified but left to further work to evaluate.

NASA is working to define the architecture needed for a Mars exploration campaign. Initial analysis assumptions allow for pre-deployment of essential cargo and equipment to support a crew landing, including the pre-positioning of a Mars Ascent Vehicle (MAV). This MAV is likely to represent the largest single payload that must be landed on the Mars surface. Its size would be influenced by the amount of mass that state-of-the-art Entry, Descent, and Landing (EDL) systems would be capable of placing on Mars. One possible method of increasing the usable size of the MAV without exceeding available EDL capabilities is to land the MAV without ascent propellant on board. Following such a method may necessitate a strategy to acquire sufficient ascent propellant to allow a crew to safely depart the Martian surface. This paper describes a conceptual return propellant strategy that uses a liquid transportation skid, or pallet, to be used in conjunction with a rover mobility system to transfer water across the Martian surface from a source point to an in-situ resource utilization (ISRU) plant that would use this water as a feedstock to generate oxygen and methane to enable launch of the MAV. Design considerations, concept of operations, and rover energetics will be discussed in this paper.

NASA’s Space Launch System (SLS) is powered by four Aerojet-Rocketdyne RS-25 engines, previously referred to as Space Shuttle Main Engines (SSMEs), with two RSRMV solid rocket boosters (SRBs) derived from the Shuttle booster as seen in Fig. 1. Ignition and the start transient of the SRBs, each of which generate 3.6 million pounds of thrust, create a series of large magnitude pressure waves separate from launch acoustics that have the potential to damage the vehicle if unmitigated. This portion of the SLS launch environment consists of igniter shock (IS), ignition overpressure (IOP), and duct overpressure (DOP).

In this paper, we report progress in the development of economically high-order flux-solution-reconstruction (FSR) schemes, which are second-order accurate on general unstructured grids but achieve high-order accuracy when a grid is regular, i.e., has the same stencil (with the same spacing) throughout the domain. Two variants of the FSR schemes are discussed: chain-rule-flux-solution reconstruction (CFSR) and quadratic-form-flux-solution reconstruction (QFSR), where the former is based on the chain rule and the latter on the flux reconstruction expressed as a function of solution variables. These schemes are tested for flows with shock waves with a limiter incorporated in the flux and solution reconstructions. Improved results, compared with second-order methods are demonstrated for inviscid and viscous flows with smooth grids.

Spanwise-periodic computation of a turbulent flow past a two-dimensional smooth ramp geometry is performed in the form of a direct numerical simulation. The Reynolds number based on the ramp height is about 147,000. A straight section that precedes the smooth ramp allows the incoming turbulent boundary layer to grow under a weak favorable pressure gradient. The boundary layer introduced at the domain inlet has a momentum-thickness based Reynolds number of 2000. The turbulent boundary layer nearing the ramp first interacts with a relatively stronger favorable pressure gradient, followed by a strong adverse pressure gradient. Consequently, the boundary layer experiences a modest acceleration before decelerating and separating. Analysis of the data over this region hints at the formation of an internal layer beneath the accelerated boundary layer. The analysis also reveals that this internal layer forms the origin of the free shear layer that emerges in the deceleration region and separates. The streamwise extent of the separated region is comparable to the ramp length, while the viscous layer thickness near reattachment is about the same as the ramp height; hence, the boundary layer undergoing separation and subsequent reattachment in the present configuration experiences its thickness being amplified by about tenfold. The reattached flow continues to develop further under a diminishing pressure gradient in the recovery region in a similar fashion to a zero pressure gradient turbulent boundary layer.

With rapidly increased worldwide interest in landing on the moon, the issue of Plume Surface Interactions (PSI) is gaining attention. Hazards posed by lander plume induced dust and debris, as well as landing site deformation can be mitigated when better understood through predictive simulations. As simulation enabling computational power continues to increase, hybrid Computational Fluid Dynamics (CFD)/Engineering models provide the immediate ability to conduct parametric/trade studies driving design decisions for landers. Reduced order erosion engineering models apply correlations of the surface erosion rate to the plume induced surface forces with the correlations anchored to flight observations from Apollo LM landings. The MSFC propulsion fluid dynamics branch has developed such a hybrid model for predicting the plume induced viscous erosive regression of the Lunar surface beneath a landing vehicle. The (Descent Interpolated Gas Granular Erosion Model) DIGGEM was originally implemented as a post processing tool to calculate induced erosion rates through vehicle descent using CFD solutions of the vehicle plume at several fixed altitudes. This capability has since been advanced in the Loci/Chem-DIGGEM model to allow transient, moving vehicle, fully coupled viscous erosion modeling of vehicle descent PSI. This model has recently been used to make preflight predictions of the erosive regression of the ground beneath a Commercial Lunar Payload Services (CLPS) vehicle in support of measurements to be made by the Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) instrument.

With rapidly increased worldwide interest in landing on the moon, the issue of Plume Surface Interactions (PSI) is gaining attention. Hazards posed by lander plume induced dust and debris, as well as landing site deformation can be mitigated when better understood through predictive simulations. As simulation enabling computational power continues to increase, hybrid Computational Fluid Dynamics (CFD)/Engineering models provide the immediate ability to conduct parametric/trade studies driving design decisions for landers. Reduced order erosion engineering models apply correlations of the surface erosion rate to the plume induced surface forces with the correlations anchored to flight observations from Apollo LM landings. The MSFC propulsion fluid dynamics branch has developed such a hybrid model for predicting the plume induced viscous erosive regression of the Lunar surface beneath a landing vehicle. The (Descent Interpolated Gas Granular Erosion Model) DIGGEM was originally implemented as a post processing tool to calculate induced erosion rates through vehicle descent using CFD solutions of the vehicle plume at several fixed altitudes. This capability has since been advanced in the Loci/Chem-DIGGEM model to allow transient, moving vehicle, fully coupled viscous erosion modeling of vehicle descent PSI. This model has recently been used to make preflight predictions of the erosive regression of the ground beneath a Commercial Lunar Payload Services (CLPS) vehicle in support of measurements to be made by the Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) instrument.

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Measuring fluid mass in micro gravity, where fluid behavior is dominated by fluid properties, is a challenging problem. To address this problem engineers at NASA are developing a capacitance-based, mass-fraction gauge for vessels containing two-phase fluids. The vessel volume is enclosed with an array of electrodes, and a unique set of capacitance measurements of the enclosed volume are made between the electrodes. The capacitance measurements are scaled with appropriate weighting factors derived from Laplace's Equation to compensate for the highly non-uniform electric fields inside the measurement volume and achieve a greater level of mass fraction accuracy.

This manual describes the installation and execution of FUN3D version 14.0.2, including optional dependent packages. FUN3D is a suite of computational fluid dynamics simulation and design tools that uses mixed-element unstructured grids in a large number of formats, including structured multiblock and overset grid systems. A discretely-exact adjoint solver may be used for formal design optimization, error estimation, and mesh adaptation. FUN3D also offers a reacting, real-gas capability and provides GPU acceleration of many common simulation options.

As the advancement of In-Situ Resource Utilization concepts and systems continue to develop, applicable technology development and maturation continues in parallel. While there are many different ways to use the resources found on other bodies, one of the most prevalent suggested applications is the manufacturing of propellants. One of the key technologies for the eventual use of these propellant production based ISRU systems is the liquefaction and storage of the produced propellants. The most mentioned propellant combinations include oxygen-hydrogen and oxygen-methane. The liquefaction of oxygen in these systems will be different than oxygen liquefaction systems on Earth, which mainly revolves around the distillation of air. These systems have been developed conceptually and many of the components have been previously tested or are in development. However, the need to demonstrate the system level operations still exists. The demonstration of a prototypical oxygen liquefaction system using tube-on-tank broad area cooling was completed to better understanding system level operations during liquefaction activities. Demonstration testing included system performance determination, constant liquefaction demonstrations, and transient liquefaction demonstrations. The demonstrations showed the operational capabilities of the tube-on-tank system with an integrated cryocooler. Additional testing explored subsurface vs ullage introduction of the gaseous oxygen flow stream as well as demonstrating a novel fiber optic sensor that measured the temperature gradients along the fluid center line within the tank.

While oxygen will probably be the first fluid liquefied from in-situ developed propellants, hydrogen will also be a sought-after commodity on both the Lunar and Martian surfaces. General hydrogen liquefaction processes are not significantly affected by the atmosphere on Mars and similar processes can be applied on both surfaces. However, hydrogen will require different processes than oxygen (which can be directly liquefied into its final use tank) due to the different spin states the hydrogen nucleus (i.e. parahydrogen and orthohydrogen) can have as well as the significant sensible energy that has to be removed prior to the actual phase change. Due to the low efficiency of 20 K class cryocoolers and high sensible heat removal requirements at higher temperatures, significant benefits can be obtained by precooling the hydrogen gas in a variety of manners. The benefits of three concepts are explored to take advantage of different precooling technics as well as accounting for sensible, latent, and para to ortho conversion heating requirements. Process concepts are shown that can reduce the base input power requirements for heat removal by more than 67% and the radiator area by similar amounts. A majority of the technology development for these processes has already begun, both within NASA and externally, with the few gaps being identified.

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics, the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal., 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer., 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy. This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

This final report documents work performed by ATA Engineering, Inc., (ATA) to develop computational models and analyze the coupled fluid-structure response of two types of noise treatments applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The first treatment is a slat-gap filler (SGF), which closes the gap between the suction surfaces of a deployed slat and an aircraft main wing, and the second treatment is a slat-cove filler (SCF), which replaces the recirculating flow on the slat cove with a surface that promotes flow attachment. The representative airframe chosen for this work was NASA’s High-Lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs were assessed for their static aeroelastic deflection as well as their dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL-CSD) tools. Specifically, fluid-structure interaction (FSI) problems were solved computationally using the CFD solver Loci/CHEM and the NL-CSD solver Abaqus. As a part of the overall project, a similar capability was implemented using the CFD solver FUN3D coupled to Abaqus, although that work is documented in a separate report and that FSI framework was not used to analyze any of NASA’s SGF and SCF configurations. The technical approach consisted of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design was analyzed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability, whereas the SCF was analyzed in 2D only. SGF static aeroelastic response predictions focused on characterizing the deformed shape, with maximum displacements predicted to be on the order of magnitude of the technology element panel thickness. SGF dynamic aeroelastic response predictions used Partial Floquet analysis of the temporal evolution of selected nodal displacements to quantify the sign and magnitude of aeroelastic damping. Results suggest that the CRM-HL operating conditions would result in a dynamically stable response. The simulated dynamic pressure was also increased up to a factor of about four, and resulting responses suggest that predicted dynamic stability would be achieved with some margin.

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m2s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m2s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

This report provides turbulent flow statistics made using particle image velocimetry (PIV) for several internally mixed exhaust systems with external plugs (Plug20 designs). The tests were conducted at the NASA Glenn Research Center’s Aero-Acoustic Propulsion Laboratory (AAPL) in May and June 2022. The corresponding acoustic data for the configurations and flows presented in this report can be found in an earlier Technical Memorandum, NASA/TM-20210010291, “Plug20 Test Report.” The Plug20 nozzle family was designed to represent commercial supersonic vehicles in the 2025 to 2035 timeframe with a bypass ratio around 4, having an internal mixer, with and without an external plug to improve nacelle boattail drag and shock. The flow conditions measured represent engine cycles for commercial supersonic aircraft at takeoff, having pressure ratios of 2.0 and heated fan and core streams, measured both with and without a flight stream surrounding the plume. A few unheated flow conditions were also tested to tie results to other baseline test results. The PIV measurements in both streamwise and cross-stream planes are described, along with brief analysis to visually catalog the contents and confirm their validity. In this analysis, the streamwise PIV is found to match benchmark datasets within its experimental error. However, the cross-stream turbulent quantities measured by two systems do not match each other or the streamwise measurements in amplitude—caution is required when using the cross-stream data. Data acquired and included in this report are (1) velocity statistics of the plumes in both a streamwise and several cross-stream planes, and (2) the hot nozzle geometry. Tabular forms of the data and computer-aided design (CAD) files of the hot geometry are presented as supplemental files described in Appendix A (TM-20230006779SUPPL1.zip) and Appendix B (TM-20230006779SUPPL2.zip).

An experimental rig was designed and constructed at the NASA Glenn Research Center to measure the alternating current (AC) losses of round, multifilament, superconducting (SC), magnesium diboride (MgB2) wire in the temperature range from 20 to 30 K. The rig resembles a permanent magnet (PM) motor with the stator windings omitted. Experiment packs containing the superconductor samples will be placed in the annular space between the rotor and the stator back iron. The original plan was to measure the AC losses by calorimetry using nonboiling liquid hydrogen as the cooling fluid. The present plan is to use helium gas (GHe) under a few bar pressure instead. This report summarizes the calculations that indicate that reasonable results can be obtained, even though the mass density of the GHe will be only a fraction of that of liquid hydrogen. It is necessary to show that the heat transfer coefficient from SC wire to GHe can be high enough, that the temperature rise of the He between the experiment inlet and outlet can be small enough, and that sufficient He flow can be driven by a commercially available He fan. Both straight wire samples to validate loss models and coils to simulate motor stator coils can be tested in the rig. The lowest superconductor temperature that is feasible to reach in the rig is estimated. The temperature range available for samples will be wider with the gaseous He system than it would have been with liquid hydrogen.

On May 5, 1961, Alan B. Shepard became the first American in space (K. Mars); since then, NASA has made leaps and bounds in space exploration and continues to develop new technologies. Through their X-Hab program, NASA has teamed up with the National Space Grant Foundation to fund our student-centered project with one main goal: Move That Goo! Liquid amine absorption is a method by which CO2 is removed from an environment and is intended to be the primary method for CO2 removal for NASA’s sustained human presence on extraterrestrial bodies such as the Moon and Mars. Liquid amines can be highly viscous when absorbing CO2, thus the team needed to produce a design that could move a high viscosity fluid while staying under the power requirement of 1 kW and allowing a few square meters in a closed system. Prompted by NASA and with the support of our faculty mentors, we have been tasked with the research, development, production, and testing of a system able to move a high viscosity liquid amine which requires adequate air exposure to absorb carbon dioxide from a habitat atmosphere in low gravity. The system is split into two main groups, fluids, and materials, with two separate teams of students assigned to their half of the design. The fluids team is responsible for everything in the design that directly touches the fluid. The materials team is responsible for the test stand to which the fluids team’s design is attached to as well as additive manufacturing material selection. In the first semester, the team focused on brainstorming ideas, design selection, budget creation, and CAD modeling. Multiple concepts were generated but the team decided on a design to get the high viscosity fluid to move around a rectangular channel using a screw conveyor mechanism. The liquid must be moved by the turning screw from the thermal chamber going around a rectangular trough back to the thermal chamber. The screw and the trough geometry will be important to allow for the high viscosity fluid to flow, utilizing surface tension to maintain contact to the trough as seen in Figure 1 below.

I worked as a NASA Intern during the Summer 2023 term in the DS-00 division under the supervision of my mentor, Dr. Jonathan Pitt. Our goal was to build on our previous work from 2022 to develop a plume strike estimation tool using a prescribed physics methodology and model plume impingement effects while considering the dynamics of a rendezvous, operations, proximity, and docking (RPOD) maneuver. This tool supports previously configured CFD-DSMC calculations by allowing for rapid analysis of initial designs using a low-fidelity source flow model. Engineers can then use the high-fidelity CFD-DSMC tool to consolidate results as they work towards finalizing a design. This year’s project was focused on developing a software application that other engineers would be using in their analysis. Thus, the user’s experience was considered in the development of this application. Proper documentation, testability, and modularity of the codebase was our priority. For example, the project included auto documentation procedures to start building towards a User Manual, while also including dedicated demonstration cases for more explicit communication of functionality. Also, this project included a framework for testing the source code for future developments. Additionally, care was taken to develop the code using an Object-Oriented Programming approach. Thus, allowing for a modular extensibility of functionality in anticipation of future developments. The core work of this project was developing an algorithm that would transform the visiting vehicle and associated thruster data according to the kinematics described in the jet firing history. It would then calculate the estimated plume strikes on two of the target vehicles and write data accordingly into a VTK file. Summer work is to conclude by developing and presenting a PowerPoint slide deck at the intern exit briefing on August 11th, 2023. Once the model for simple plume strike calculations is developed and tested there are several avenues to explore to continue development of this tool. These are also discussed in this report.

This paper describes the development of a thermal management system (TMS) concept design and analysis software package called Heat Transport System Simulation (HeaTSSPy). Built within Python using the OpenMDAO framework, HeaTSSPy can be used to size and optimize an active (using liquid/air heat exchangers) or passive (using finned heat sinks) TMS. The package makes use of modular TMS elements that allow for the creation of different system architectures and includes components such as heat sinks, heat exchangers, liquid pumps, fans, ducts, air inlets, air nozzles, and liquid pipes. Modeling methods for these components include a combination of physics-based analytical and empirical equations that relate component sizing criteria to system performance. The HeatSSPy heat sink methods are fully detailed within this paper, while the heat exchanger methods are described and referenced from previous work. This paper also uses high-fidelity simulation to validate two different methods of calculating thermal resistance using CFD results. Once the methods are fully described, the code is exercised to compare an active TMS with that of a passive TMS. Design criteria for this study include rejected heat, system altitude, Mach number, and ambient temperature. These criteria are used to develop a TMS system with estimated performance metrics such as weight, drag, and operational power. Results of this paper show the crossover point when a passive system begins to weigh more than an active system for a given heat rejection.

Thrusters are used to maneuver spacecraft in space. Spacecraft maneuvering and attitude control is achieved by accelerating the high temperature and pressure exhaust gases, which are the result of chemical combustion in the case of chemical engines, through the rocket nozzle. Following the end of the nozzle a plume is created containing combustion gases and high velocity droplets of unburnt propellant that with time can erode the surfaces it strikes. A simulation that can accurately model a thruster’s plume can help us reduce the erosion of components struck by the plume. The creation of this simulation can be broken into several parts, and this paper focuses on finding out what distribution of particles along a plane in the nozzle, where combustion can be assumed to be complete, creates the observed particle distribution in the plume. Findings show that unburnt droplets are mostly concentrated in the center of the nozzle or uniformly dispersed within the nozzle.

NASA’s return to the moon brings about many challenges, including issues with survival on the Lunar surface. Marshall Space Flight Center (MSFC) has been developing a hybrid thermal control system that can be utilized for various surface assets that must survive in extreme lunar environments. This scalable system is targeted for human-rated systems and utilizes a combination of a- pumped fluid loop (PFL) and a loop heat pipe (LHP) with thermal control valve (TCV). Pumped fluid loops have a long history of use in human rated systems. They can collect large amounts of waste heat and transport it over long distances in rovers, habitats, and other systems. By utilizing a non-toxic working fluid in the habitable volume, the PFL is easily serviceable by the crew and does not pose a risk during any unexpected leaks or failures. The addition of a LHP for the exterior heat transport and rejection adds several benefits to the system for extreme environment survival. The quantity of LHP scan be tailored to optimize heat rejection for different systems and allows for larger adiative surfaces. By utilizing a TCV in combination with the LHP, there can be high heat transfer during the daytime with minimal heat transfer during the night. The TCV passively controls the amount of heat transfer through the LHP based on the environmental temperature. A prototype system has been tested in bench top and thermal vacuum conditions. Results are presented showing the feasibility of the system concept.

The storing and transfer of cryogenic propellants is an enabling technology for NASA as it moves to expand future missions into low Earth orbit (LEO) and beyond. The transfer of cryogenic propellants is of particular interest because it has never been successfully demonstrated in microgravity. Several ground tests have been conducted in past decades examining various tank chilldown and fill methods. Creating accurate models using these historical datasets is a vital step toward developing appropriate modeling tools to form pre-test predictions for eventual flight and ground tests of various novel injector concepts. The goal of this work was to create a working model in Thermal Desktop using data from the CRYOTE-2 experiments whose purpose was to demonstrate the chill and fill process in a spherical receiver tank parametrically examining the effect of injection method, receiver tank initial wall temperature and initial state, supply pressure, and mass flow rate on the NVF process. These tests were conducting using liquid nitrogen with different injectors with different injection patterns. The focus of the current work is on one particular spray nozzle. The pressure rise inside the tank, the fill level, and the change in the outside wall temperature were all tracked and compared with test data. The Thermal Desktop model was able to accurately predict each of these metrics for this nozzle within 11.0%, 22.0%, and 7.0%, respectively.

The storing of cryogenic propellants is an enabling technology for NASA as it moves to expand future missions into low Earth orbit (LEO) and beyond. The Robotic Refueling Mission–3 (RRM3) experiment launched in 2018 and was designed to demonstrate tank self-pressurization during storage and active cooling of liquid methane in microgravity. Creating accurate models using existing datasets from flight experiments such as RRM3 is a vital step toward developing appropriate modeling tools to form pre-test predictions for future flight systems. The goal of this work was to create a working model in Thermal Desktop of the RRM3 source Dewar during the microgravity self-pressurization experiment. The pressure rise, temperature distribution of the fluid inside the tank, and the wall temperatures were all tracked and compared with test data. Initial results show good agreement between each of these metrics and the experimental results. This work is ongoing and will yield further insights into the thermal behavior of cryogens in microgravity.

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) activity has been investigating concepts for the liquefaction of cryogenic fluids produced on the surface of the Moon and Mars. Liquefaction testing of oxygen in a scaled lander tank using an integrated, industrial cryocooler system was completed. The testing covered the determination of the nominal performance operation of the system, constant liquefaction performance, and transient liquefaction performance. The goal of the testing was to demonstrate cryogenic liquefaction operational capabilities on the Lunar and Martian surfaces for landers and In Situ Resource Utilization. The testing met all key performance parameter threshold values including liquefaction rate (demonstrated 1.6 kg/hr compared to threshold of 1.1 kg/hr), number of variables as transient (goal of 3: liquefaction flow, environmental temperature, and cryocooler input power – all demonstrated), and number of fill levels (goal of 3: < 5%, 50%, and 90% - all demonstrated). Key testing results from oxygen liquefaction testing will be discussed including sensitivities provided by the analysis team.

The Extended Stroke, Reverse Operation, Differential Thermal Expansion (DTE) Thermal Switch (TS) for Extreme Variable Environments uses materials with both positive and negative coefficients of thermal expansion (CTE) arranged in concentric cylinders so that the TS operates passively, opening a gap between interfaces to create a low thermal contact when cold, and closing the gap between interfaces to create a high thermal contact when hot. This paper details the testing of the TS at Marshall Space Flight Center (MSFC) including the vibration (vibe) test, pre-vibe thermal vacuum (t-vac) cycle test, post-vibe t-vac cycle test, and long duration t-vac cycle test. The data collected from the pre- and post-vibe tests are compared to show the impact of the vibe test on performance. The data collected from the long duration test is used to calculate the total thermal contact conductance through the TS interfaces for both ‘on’ (TS closed, during hot plateau) and ‘off’ (TS open, during cold plateau) configurations, and the turndown ratio is also calculated. Challenges experienced throughout testing are also discussed.

This paper examines the use of actively controlled louvers to attenuate UV and dust, as well as mitigate freezing concerns for human spacecraft radiators during Artemis missions. Artemis missions to the lunar orbit or surface will expose the radiators to high energy UV radiation and dust, which will degrade the radiator’s coating emissivity and consequently reduce heat rejection performance. In addition, subfreezing environmental temperatures during transit to lunar orbit and nighttime on lunar south pole can rupture coolant tubes, reduce heat rejection performance, and worst-case scenario result in a Loss of Mission (LOM). Louver technology would be a promising solution to maintaining radiator performance and integrity for Artemis missions, but heritage louvers are passively controlled. This technology needs maturing to active control, or motor actuation, to achieve faster thermal response times. Actively controlled louver design considerations are discussed in this paper. The analysis that follows shows actively controlled louvers can attenuate high energy UV radiation and dust, as well as protect the coolant from freezing.

For models with a limitation on the overall node count, the typical approach to Multi-Layer Insulation (MLI) modeling may generate nodes that are necessary for the analysis, but do not represent components of particular interest. This leaves fewer nodes that can be utilized to model components of greater importance than the MLI. A common approach to modeling MLI is to include a separate MLI node representing the outer layer of the insulation and a radiative coupling based on the area multiplied by an effective emissivity. Therefore, wherever insulation is included, one node is needed for the underlying surface and another node for the insulation. Since many spacecraft and instruments include MLI covering a sizable portion of their designs, this may result in a considerable number of nodes being used for MLI. An alternate method to MLI modeling was developed that eliminates the MLI node, while still preserving the effect of the insulation for the underlying surface, thereby increasing the available nodes that could be used elsewhere in the model. This approach relies on preserving the baseline reflectivity, while reducing the absorptivity (based on the blanket effective emittance) and including a transparency. An inactive second surface is placed just behind the base surface that fully absorbs any energy that is transmitted without including its effect in the model. In essence, this approach applies only the energy that makes it through the blanket to the underlying surface. This method was tested out on the Roman Space Telescope model in local areas in preparation for its use in the generation of a launch model, which is constrained in the allowable node count. This paper documents the performance of the method and presents a comparison between the One-Node MLI method and the traditional two node MLI approach.

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

The Generalized Fluid System Simulation Program (GFSSP) computer code is enhanced with the capability to model condensation in the presence of noncondensable gases. Condensation in the presence of noncondensables is modeled using the Couette flow film (stagnant film) model. Experimental data on the condensation of water vapor in downward flow of air-water vapor and helium-water vapor mixtures in vertical tubes are compared with the predictions of GFSSP. The comparisons show that GFSSP can predict the experimental data well.

Liquid nitrogen is often used as a substitute for oxygen testing due to safety concerns, and the general similarity between their fluid properties. During the Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) testing, an opportunity arose to compare liquid nitrogen and liquid oxygen behavior using the same test hardware in similar test conditions. Comparative testing would verify whether the system response to nitrogen and oxygen would behave in a similar manner. Tests that were investigated include system boil-off heat load determination, autogenous pressurization, cryocooler loop initiation, and cryocooler loop operations. Results from testing are shown and compared to verify that liquid nitrogen and liquid oxygen tests yield similar system responses, and that nitrogen can be used as a substitute for oxygen in developmental tests at cryogenic temperatures.

The use of cryogenic propellants has and will continue to play an integral role in manned-space exploration due the high specific impulses offered and its ubiquity through in-situ resource utilization. But guaranteeing vapor-free transfer of such low-surface tension liquids is difficult for traditional capillary-action propellant management devices (PMDs). Screen-channel liquid acquisition devices and compliant origami bladders are potential solutions, but to quantify the benefits these technologies offer, this paper present a cast study analyzing the performance metrics of an orthodox vane PMD as a benchmark for comparison. A 0.15 m3 liquid hydrogen tank at 20.3 K and 103 kPa was selected for study. Then assuming no body forces and no heat transfer (for simplicity), a steady-state 1-D differential equation was numerically solved in tandem with four possible wetted area configurations to yield an expulsion efficiency for a theoretically maximum inputted expulsion flow rate. Maximization of the flow rate was constrained by the onset of choked flow. Additional inputs, including vane height and vane number, were parametrically varied between 0.1 – 10 cm and 4 – 28 vanes, respectively.

Understanding two-phase cryogenic propellant behavior is key to enabling technologies for future spaceflight missions. Developing accurate models of two-phase flow phenomena, particularly flow boiling in the heating configuration is relevant to the propellant transfer process both in microgravity and on other planetary surfaces. Currently there is a need for more accurate, direct cryogenic data anchored models for various boiling phenomena. Recently, universal correlations for cryogens flowing in heated tubes have been developed for a wide variety of fluids, thermodynamic conditions, and various regimes across the boiling curve, and have been patched to provide a smooth, continuous predictive curve. This paper demonstrates how these correlations have been ported into Thermal Desktop to improve predictive performance. To test the validity of the new correlations, Thermal Desktop is validated against a historical dataset of flow boiling experiments in the heating configuration using liquid helium. Based on results, the new correlations show a substantial improvement over the original built-in flow boiling correlations in Thermal Desktop in predicting the wall temperature as a function of preponderant parameters for this quantum fluid.

Available experimental data dealing with critical heat flux (CHF) of liquid hydrogen (LH2), liquid methane (LCH4), and liquid oxygen (LO2) in pool and flow boiling are compiled. The compiled data are compared with widely used correlations. Experimental pool boiling CHF data for the aforementioned cryogens are scarce. Based on only 25 data points found in five independent sources, the correlation of Sun and Lienhard (1970) is recommended for predicting the pool CHF of LH2. Only two experiments with useful CHF data for the pool boiling of LCH4 could be found. Four different correlations including the correlation of Lurie and Noyes (1964) can predict the pool boiling CHF of LCH4 within a factor of two for more than 70% of the data. Furthermore, based on the 19 data points taken from only two available sources, the correlation of Sun and Lienhard (1970) is recommended for the prediction of pool CHF of LO2. Flow boiling CHF data for LH2 could be found in seven experimental studies, five of them from the same source. Based on the 91 data points, it is suggested that the correlation of Katto and Ohno (1984) be used to predict the flow CHF of LH2. No useful data could be found for flow boiling CHF of LCH4 or LO2. The available databases for flow boiling of LCH4 and LO2 are generally deficient in all boiling regimes. This deficiency is particularly serious with respect to flow boiling.

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.