Search NASASearch

SEARCH · Search NASA

Results for “sound velocity anomalies”

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.

96 records · Page 5

Summary of recent NASA propeller research

Advanced high speed propellers offer large performance improvements for aircraft that cruise in the Mach 0.7 to 0.8 speed regime. At these speeds, studies indicate that there is a 15 to near 40 percent block fuel savings and associated operating cost benefits for advanced turboprops compared to equivalent technology turbofan powered aircraft. Recent wind tunnel results for five eight to ten blade advanced models are compared with analytical predictions. Test results show that blade sweep was important in achieving net efficiencies near 80 percent at Mach 0.8 and reducing nearfield cruise noise about 6 dB. Lifting line and lifting surface aerodynamic analysis codes are under development and some results are compared with propeller force and probe data. Also, analytical predictions are compared with some initial laser velocimeter measurements of the flow field velocities of an eight bladed 45 swept propeller. Experimental aeroelastic results indicate that cascade effects and blade sweep strongly affect propeller aeroelastic characteristics. Comparisons of propeller nearfield noise data with linear acoustic theory indicate that the theory adequately predicts nearfield noise for subsonic tip speeds, but overpredicts the noise for supersonic tip speeds.

Daniel C Mikkelson

Extensional Flow Convecting a Reactant Undergoing a First Order Homogeneous Reaction and Diffusional Mass Transfer From a Sphere at Low to Intermediate Peclet and Damkohler Numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 ≤ Pe ≤ 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 ≤ Da ≤ 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r,Θ) = Σ c n (r)P n (cosΘ) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions P n (cosΘ), the variable coefficient PDE can be reduced to a system of N+1 second order ODEs for the radial functions C n (r), n=0,1,2, ... ,N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N+ 1) first order differential equations as a two-point boundary value problem on 1 ≤ r ≤ r •• . Convergence of the expansion for a specific value of N can thus be established and provides "spectral" behavior as well as the full concentration field c(r,Θ).

N Y Shah

Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies

Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.

Mars2020

Behavior of Ullage Bubbles During Blowdown in Low-g Experiment (BUBBLE): Overview of a Cryogenic Tank Depressurization Test

Behavior of Ullage Bubbles during Blowdown in Low-g Experiment (BUBBLE) is a cryogenic tank depressurization drop tower experiment currently being designed at NASA Glenn Research Center. There is a need for liquid level rise data while venting a cryogenic propellant tank below the liquid saturation pressure in a reduced gravity environment. When a cryogenic tank is vented to a vapor pressure below the liquid saturation pressure, bubbles become entrained in the liquid causing the bulk liquid-vapor interface to rise. In reduced gravity, bubble rise velocity is diminished leading to a larger liquid level rise compared to a 1g environment. The purpose of the experiment is to gain further understanding of tank depressurization fluid physics and obtain high-fidelity data for model validation in 1g and reduced gravity environments. On-orbit cryogenic propellant tanks operating at high fill levels must efficiently manage venting operations to avoid the risk of liquid entrainment in the vent line, which could lead to asymmetric control thruster loads, freezing and clogging, and loss of liquid propellant. Validated models could be used to design settling and venting profiles to reduce risk and increase efficiency for cryogenic storage and transfer operations.

tank venting

A Comparison of Doppler Lidar Wind Sensors for Earth-Orbit Global Measurement Applications

Now, there are four Doppler lidar configurations which are being promoted for the measurement of tropospheric winds: (1) the coherent CO2 Lidar, operating in the 9 micrometer region using a pulsed, atmospheric pressure CO2 gas discharge laser transmitter, and heterodyne detection; (2) the coherent Neodymium doped YAG or Glass Lidar, operating at 1.06 micrometers, using flashlamp or diode laser optical pumping of the solid state laser medium, and heterodyne detection; (3) the Neodymium doped YAG/Glass Lidar, operating at the doubled frequency (at 530 nm wavelength), again using flashlamp or diode laser pumping of the laser transmitter, and using a high resolution tandem Fabry-Perot filter and direct detection; and (4) the Raman shifted Xenon Chloride Lidar, operating at 350 nm wavelength, using a pulsed, atmospheric pressure XeCl gas discharge laser transmitter at 308 nm, Raman shifted in a high pressure hydrogen cell to 350 nm in order to avoid strong stratospheric ozone absorption, also using a high resolution tandem Fabry-Perot filter and direct detection. Comparisons of these four systems can include many factors and tradeoffs. The major portion of this comparison is devoted to efficiency. Efficiency comparisons are made by estimating the number of transmitted photons required for a single pulse wind velocity estimate of + or - 1 m/s accuracy in the middle troposphere, from an altitude of 800 km, which is assured to be reasonable for a polar orbiting platform.

Robert T Menzies

Timing the Flames: Geostationary Satellite Detection of Diurnally Shifting Stubble Burning in Northwestern India

Post-monsoon open-field stubble burning in northwestern (NW) India—a key agricultural region known as the “breadbasket”—is a longstanding practice used to clear fields. Satellite observations spanning over two decades have revealed significant upward trends in crop production, vegetative greenness, and the frequency of post-harvest fires, with this last contributing to hazardous air quality during the peak burning season (mid-October to mid-November). Since 2022, thermal anomaly data from Aqua-MODIS and SNPP-VIIRS sensors have shown a sharp decline in reported fire events—an observation that contrasts starkly with the concurrent rise in regional aerosol loading detected from space. This apparent discrepancy became particularly pronounced in 2024–2025, prompting a closer examination using high-temporal-resolution imagery from the Advanced Meteorological Imager (AMI) on the geostationary satellite GEO-KOMPSAT-2A. These observations revealed a clear spike in fire-related signals occurring around and after 4:00 p.m. local time, i.e., outside the typical noon to 2:00 p.m. detection window of the MODIS and VIIRS. A fire detection algorithm exploiting the fire-sensitive shortwave-infrared 3.8 μm signal and its contrast to 11.2 μm infrared observations is designed to adopt AMI observations and applied to its multi-year observations (2019–2025). The resulting fire dataset unambiguously shows a gradual shift in stubble burning activity toward the late afternoon hours beginning in 2022 which is underreported by polar-orbiting satellites. The orbital drift of NASA’s MODIS sensor on the Aqua platform allows detection of some of the gradually shifting fires during afternoon hours, but the MODIS still misses a large number of fires occurring around and after 4 p.m. The AMI’s relatively coarse spatial resolution (~4 km), a consequence of its slant viewing geometry over NW India, imposes inherent limitations on quantifying the full extent of fire occurrences. The operational air quality forecasting models currently assimilate satellite fire detections predominantly captured during early afternoon overpasses of the MODIS and VIIRS. The temporal shift in fire activity complicates such forecast, leading to a substantial underestimation of emissions. Intense stubble burning and the resulting air pollution highlight the need for effective crop residue management practices for mitigating the frequency of open biomass burning and thereby reducing episodic degradation of air quality and its associated public health and economic impacts.

post-monsoon stubble burning; northwestern India;

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.

Trailing Edge

Engineering the Interface: Advanced Surface Technologies for Lunar Dust Management and Equipment Longevity

Through the Artemis program, NASA intends to develop a sustainable human foothold on the Moon, ultimately paving the way for crewed exploration of Mars. The Moon's hostile environment poses numerous obstacles, including exposure to radiation, temperature extremes, micrometeoroid threats, and particularly the persistent problem of lunar dust. Lunar dust impacts nearly every aspect of surface operations through adhesion and abrasion mechanisms, with contamination from anthropogenic activities (landing, rovers) far outweighing natural phenomena. Multiple adhesion pathways contribute to surface contamination in the lunar environment, including van der Waals forces, electrostatic forces, chemical reaction, and magnetic forces from elemental iron deposits. Sharp asperities from micrometeoroid bombardment and atmospheric absence increase interaction potential and enable mechanical interlocking. Low cohesion between dust particles exacerbates these challenges, as minimal interaction potential between dust and nearby surfaces overcomes particle cohesion, causing contamination. Lunar dust adhesion mitigation technologies can be categorized as either active, requiring external energy, or passive, relying on intrinsic material properties. Ultrasonic and electrodynamic technologies have been developed to the highest technology readiness level for active approaches. Passive strategies primarily focus on surface chemistry and topography modifications. At NASA Langley Research Center, approaches include surface migration agents to reduce surface energy, topographical modification using laser ablation patterning, and tailored surface conductivity to reduce intrinsic adhesion force. Performance has been evaluated using custom-built ultrasonic and centrifuge instruments. Plume-surface interactions from lunar landers can propel micrometer-sized particles at velocities up to 1000 m s-1.8 These particles pose risks to landers, habitats and infrastructure, leading to erosion, degradation, and reduced component lifespan. A panel recovered from Surveyor III was determined to have been severely abraded because of lunar dust displaced from the Apollo 12 lunar module that landed 160 m away. The performance of metallic surfaces has been evaluated via high velocity single particle impact using the laser-induced project impact test (LIPIT) facility at the University of Utah. Peridynamics modeling, a form of continuum mechanics that uses a nonlocal approach enabling greater simulation capabilities of crack initiation and fracture, has also been utilized to gain greater insight into material response during impact events. Lunar dust contamination challenges extend to power generation systems and moving equipment. Cables, rotation stages, and other mechanisms may experience limited range of motion and reduced lifetime due to dust infiltration. NASA Langley Research Center has evaluated traditional aerospace alloys, softgoods, wear resistant ceramics, and several polymer and polymer composite materials. Test methods have included traditional techniques like Taber abrasion testing, as well as designed test configurations developed in the DUSTE (dust, ultraviolet radiation, and space thermal environmental) chamber that reproduce mechanism functions in operational environment. Beyond laboratory experiments, several flight experiments have been conducted. Materials were exposed to the low Earth orbit environment on the Materials International Space Station Experiment (MISSE) and to the lunar surface environment through the Aegis Aerospace Regolith Adherence Characterization (RAC) payload and the Honeybee Robotics PlanetVac payload. Determining lunar dust's impact on surface exploration and habitation requires comprehensive experimental and computational capabilities combined with lessons learned from initial lunar activities. Identifying the greatest environmental challenges and developing mitigation technologies provides the clearest path toward successfully, expeditiously, and efficaciously completing NASA's mission. This presentation will discuss ongoing efforts at NASA Langley Research Center and collaborator contributions to these critical objectives.

Surface Engineering

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations

Lubricating Properties of Ceramic-Bonded Calcium Fluoride Coatings on Nickel-Base Alloys From 75° to 1900° F

The endurance life and the friction coefficient of ceramic-bonded calcium fluoride (CaF2) coatings on nickel-base alloys were determined at temperatures from 75° to 1900° F. The specimen configuration consisted of a hemispherical rider (3/16-in. rad.) sliding against the flat surface of a rotating disk. Increasing the ambient temperature (up to 1500° F) or the sliding velocity generally reduced the friction coefficient and improved coating life. Base-metal selection was critical above 1500° F. For instance, cast Inconel sliding against coated Inconel X was lubricated effectively to 1500° F, but at 1600° F severe blistering of the coatings occurred. However, good lubrication and adherence were obtained for René 41 sliding against coated René 41 at temperatures up to 1900° F; no blisters developed, coating wear life was fairly good, and the rider wear rate was significantly lower than for the unlubricated metals . Friction coefficients were 0.12 at 1500° F, 0.15 at 1700° F, and 0.17 at 1800° and 1900° F. Because of its ready availability, Inconel X appears to be the preferred substrate alloy for applications in which the temperature does not exceed 1500° F. René 41 would have to be used in applications involving higher temperatures. Improved coating life was derived by either preoxidizing the substrate metals prior to the coating application or by applying a very thin (<0.0002 in.) burnished and sintered overlay to the surface of the coating. Preoxidation did not affect the friction coefficient. The overlay generally resulted i n a higher friction coefficient than that obtained without the overlay. The combination of both modifications resulted in longer coating life and in friction coefficients intermediate between those obtained with either modification alone.

Nickel alloy

Resistance of Silicon Nitride Turbine Components to Erosion and Hot Corrosion/Oxidation Attack

Silicon nitride turbine components are under intensive development by AlliedSignal to enable a new generation of higher power density auxiliary power systems. In order to be viable in the intended applications, silicon nitride turbine airfoils must be designed for survival in aggressive oxidizing combustion gas environments. Erosive and corrosive damage to ceramic airfoils from ingested sand and sea salt must be avoided. Recent engine test experience demonstrated that NT154 silicon nitride turbine vanes have exceptional resistance to sand erosion, relative to superalloys used in production engines. Similarly, NT154 silicon nitride has excellent resistance to oxidation in the temperature range of interest - up to 1400 C. Hot corrosion attack of superalloy gas turbine components is well documented. While hot corrosion from ingested sea salt will attack silicon nitride substantially less than the superalloys being replaced in initial engine applications, this degradation has the potential to limit component lives in advanced engine applications. Hot corrosion adversely affects the strength of silicon nitride in the 850 to 1300 C range. Since unacceptable reductions in strength must be rapidly identified and avoided, AlliedSignal and the NASA Lewis Research Center have pioneered the development of an environmental life prediction model for silicon nitride turbine components. Strength retention in flexure specimens following 1 to 3300 hour exposures to high temperature oxidation and hot corrosion has been measured and used to calibrate the life prediction model. Predicted component life is dependent upon engine design (stress, temperature, pressure, fuel/air ratio, gas velocity, and inlet air filtration), mission usage (fuel sulfur content, location (salt in air), and times at duty cycle power points), and material parameters. Preliminary analyses indicate that the hot corrosion resistance of NT154 silicon nitride is adequate for AlliedSignal's initial engine applications. Protective coatings and/or inlet air filtration may be required to achieve required ceramic component lives in more aggressive environments.

Nitride

Petermann Glacier on the Brink: Progress, Challenges and Insights

Petermann Glacier, the largest marine-terminating glacier in northern Greenland based on catchment area and ice discharge, plays a key role in regulating ice discharge from the Greenland Ice Sheet into the Arctic Ocean. With an upstream catchment connected to the ice sheet interior via a deep subglacial canyon, its future stability has major implications for sea level rise. In this review, we synthesize recent advances in understanding Petermann’s dynamics across three critical interfaces: the ice–ocean, ice–atmosphere, and ice–bed boundaries. At the surface, observations show that reanalysis products underestimate air temperatures and melt, while regional climate models diverge significantly in their estimates of surface mass balance, underscoring the need for improved in situ data and models. At the ocean boundary, enhanced basal melting driven by both subglacial runoff and Atlantic water intrusions is identified as the dominant driver of recent mass loss of Petermann Glacier, with continued warming posing a serious threat to the stability of the floating tongue. At the bed, new geophysical synthesis reveals complex geology, likely spatial variability in geothermal heat flux, and the influence of the megacanyon on seasonal hydrology and velocity fluctuations. Petermann’s mass balance has been negative in recent decades without corresponding flow acceleration. However, the glacier has undergone significant calving events, and the current rifting that began in September 2025 highlights its vulnerability. This upcoming calving event underscores the timeliness of this review, as it will put Petermann Glacier’s terminus at its most retreated position since records began in 1923. The anticipated retreat of the ice tongue also reduces buttressing and brings the terminus closer to the grounding zone, and modeling studies suggest that calving within 12 km of the grounding zone could potentially trigger dynamic retreat, accelerating ice discharge, and a doubling of flow speeds. We conclude that improved observations, sustained monitoring of oceanographic and atmospheric properties, and high-resolution modeling are critical to constraining projections of Petermann Glacier’s future and its role in the stability of the Greenland Ice Sheet.

Dominik Fahrner

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Manned Lunar Landing Via Rendezvous

In any mission description, the vehicles, the flight profiles, and the astrionics hardware to im- plement the mission are all tightly interwoven things. A final result evolves only after many iterations to the solution are made. This paper will describe one of these iterations in the Saturn C-5 Earth Orbit Rendezvous approach to the Manned Lunar Landing Program. Since the iter- ation to be described in an nthone, there exists some basis for the hope that the perturbation from the final solution is small. This paper is not concerned with the landing itself, but only with those operations leading to injection of the space craft into the lunar trans- fer trajectory. However, as is to be expected, it is the target conditions which set the pace for the overall operation. The entire operation must be sized to culminate at a time and place which places the lunar target in an attainable position. The procedure would call for a burst of activity lasting over a relatively short time as compared to the long and extensive preparations leading up to it. The activity must be aimed at the opening of the lunar "launch window". Figure 1 illustrates the variation in the velocity increment required to launch a vehicle into a lunar transfer trajectory from a 485 kilometer earth orbit. The minima are at irregularly spaced intervals and are a func- tion of the inclination of the lunar and the earth satellite planes and of the position of the moon in its orbit around the earth (i. e. , the day of the month). In an operations analysis these spacings will influence the number of vehicles on the launch pad (primary and back-up) , their state of readi- ness, the firing rate, and also the flight profile to be chosen. Whether it is decided to go by "con- necting" or by "tanking" mode, the objective must be to get the spacecraft in the launch ready state at the opening of one of these launch windows. It may be desired that the first vehicle be capable of remaining in a functionally capable state even after bridging one or more of the gaps between the windows. This consideration will influence the design of the vehicles as well as the operational modes to be designed into the flight control hard- ware. For example, a sleep switch may be de- sirable from the standpoint of savings in battery we ight.

LUNAR LANDING

Investigation of Bipropellant Plume-Induced Contamination Effects on Coverglass Materials

Contamination and degradation of external spacecraft materials by unburned and partially combusted species from bipropellant thruster plumes has long been observed as a key component of the induced space environment. Space shuttle flight experiments and returned flight hardware from the International Space Station (ISS) have both experienced microscopic impact features induced by high-velocity thruster plume droplets. Analytical results have shown that droplet impingement angle relative to a receiving surface plays a key role in the surface damage. Although impacts with normal impingement angles contribute more severely to surface degradation than highly oblique angles, surface effects at higher impingement angles should not be dismissed. Thruster plume-induced materials degradation is a complex phenomenon that depends on a variety of parameters, including but not limited to material type, system temperature and pressure, plume composition, and thruster firing specifications such as number of pulses, pulse duration, and sample distance from the thruster. For space applications, attaining the vacuum pressure and temperature conditions necessary for flight-like plume expansion and exposure conditions is not a trivial task. The German Aerospace Center (Deutsches Zentrum für Luft- und Raumfahrt, DLR) is a facility uniquely capable of simulating such conditions. Test coupons were exposed to bipropellant thruster firings under high vacuum at the DLR facility. Percent area coverage (PAC) and droplet size distributions were evaluated for the uncoated and coated solar array coverglass materials over a range of impingement angles (0̊ to 75̊). A post-test imaging workflow was developed that aimed to quantify changes in sample surface morphology obtained from scanning electron microscopy (SEM) images using the Image Processing and Analysis in Java (ImageJ) tool; an opensource image processing software. The goal was to create a framework through which to evaluate the effect of bipropellant-induced PAC and droplet size distribution on solar array coverglass optical transmission losses. Understanding this relationship is important because optical transmission losses are known to lead to current reduction in solar power generation systems. In addition to the development of surface characterization workflows, valuable lessons learned as they pertain to future investigations and experiments will be discussed. The authors hope that sharing these lessons will facilitate more utilization of DLR’s unique capabilities as well as open the conversation for how best to address experimental characterization of flight-like plume expansion and its impacts on materials surface degradation effects.

Gateway

Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer

With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.

Microgravity