Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective conductivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 181 records · Page 10

Liquid metal micro heat pipes for space radiator applications

Micromachining is a chemical means of etching three-dimensional structures, typically in single-crystalline silicon. These techniques are leading toward what is coming to be referred to as MEMS (micro electro mechanical systems), where in addition to the ordinary two dimensional (planar) microelectronics, it is possible to build three-dimensional micromotors, electrically-actuated microvalves, hydraulic systems, and much more on the same microchip. These techniques become possible because of differential etching rates of various crystallographic planes and materials used for semiconductor microfabrication. The University of Cincinnati group in collaboration with NASA Lewis formed micro heat pipes in silicon by the above techniques. Work is ongoing at a modest level, but several essential bonding and packaging techniques have been recently developed. Currently, we have constructed and filled water/silicon micro heat pipes. Preliminary thermal tests of arrays of 125 micro heat pipes etched in a 1 inch x 1 inch x 250 micron silicon wafer have been completed. These pipes are instrumented with extremely small P-N junctions to measure their effective conductivity and their maximum operating power. A relatively simple one-dimensional model has been developed in order to predict micro heat pipes' operating characteristics. This information can be used to optimize micro heat pipe design with respect to length, hydraulic diameter, and number of pipes. Work is progressing on the fabrication of liquid-metal micro heat pipes. In order to be compatible with liquid metal (sodium or potassium), the inside of the micro heat pipes will be coated with a refractory metal (such as tungsten, molybdenum, or titanium).

Gerner, F. M.↗

The Center for Astrophysics in Antarctica

Nowhere on Earth are the infrared skies clearer, darker, or more stable than on the high Antarctic Plateau. At some wavelengths, Antarctic telescopes may be more than one to two orders of magnitude more efficient than at other sites. However, exploiting these advantages requires first addressing the formidable practical difficulties of working in the remote and frigid polar environment. This was the motivation for the Center for Astrophysical Research in Antarctica (CARA), one of twenty-five National Science Foundation Science and Technology Centers. At its inception, the Center organized its research into four projects. Three - AST/RO, COBRA, and SPIREX - address key problems in star formation, evolution of galaxies, and the distribution of matter in the early universe. They feature surveys which can be conducted effectively with moderate-size telescopes operated in a highly automated mode. They also explore the potential of the Antarctic Plateau for a broad range of astrophysical research over a spectral range extending from the near-infrared to millimeter wavelengths. A fourth, ATP, was created to obtain quantitative data on the qualities of the South Pole site and to plan for future scientific projects. During the next five years, AST/RO, COBRA, and SPIREX will become operational, and the Center will begin to build a second generation of telescopes which can address a broader range of problems and accommodate a larger community of users.

Pernic, Robert J.↗

Near Surface Vapor Bubble Layers in Buoyant Low Stretch Burning of Polymethylmethacrylate

Large-scale buoyant low stretch stagnation point diffusion flames over solid fuel (polymethylmethacrylate) were studied for a range of aerodynamic stretch rates of 2-12/ sec which are of the same order as spacecraft ventilation-induced stretch in a microgravity environment. An extensive layer of polymer material above the glass transition temperature is observed. Unique phenomena associated with this extensive glass layer included substantial swelling of the burning surface, in-depth bubble formation, and migration and/or elongation of the bubbles normal to the hot surface. The bubble layer acted to insulate the polymer surface by reducing the effective conductivity of the solid. The reduced in-depth conduction stabilized the flame for longer than expected from theory neglecting the bubble layer. While buoyancy acts to move the bubbles deeper into the molten polymer, thermocapillary forces and surface regression both act to bring the bubbles to the burning surface. Bubble layers may thus be very important in low gravity (low stretch) burning of materials. As bubbles reached the burning surface, monomer fuel vapors jetted from the surface, enhancing burning by entraining ambient air flow. Popping of these bubbles at the surface can expel burning droplets of the molten material, which may increase the fire propagation hazards at low stretch rates.

Olson, Sandra L.↗

Arc Jet Testing of the TIRS Cover Thermal Protection System for Mars Exploration Rover

This paper summarizes the arc jet test results of the Mars Exploration Rover (MER) Silicone Impregnated Reusable Ceramic Ablator (SIRCA) Transverse Impulse Rocket System (TIRS) Cover test series in the Panel Test Facility (PTF) at NASA Ames Research Center (ARC). NASA ARC performed aerothermal environment analyses, TPS sizing and thermal response analyses, and arc jet testing to evaluate the MER SIRCA TIRS Cover design and interface to the aeroshell structure. The primary objective of this arc jet test series was to evaluate specific design details of the SIRCA TIRS Cover interface to the MER aeroshell under simulated atmospheric entry heating conditions. Four test articles were tested in an arc jet environment with various sea] configurations. The test condition was designed to match the predicted peak flight heat load at the gap region between the SIRCA and the backshell TPS material, SLA-561S, and resulted in an over-test (with respect to heat flux and heat load) for the apex region of the SIRCA TIRS Cover. The resulting pressure differential was as much as twenty times that predicted for the flight case, depending on the location, and there was no post-test visual evidence of over-heating or damage to the seal, bracket, or backshell structure. The exposed titanium bolts were in good condition at post-test and showed only a small amount of oxidation at the leading edge locations. Repeatable thermocouple data were obtained and SIRCA thermal response analyses were compared to applicable thermocouple data. For the apex region of the SIRCA TIRS Cover, a one-dimensional thermal response prediction proved overly conservative, as there were strong multi-dimensional conduction effects evident from the thermocouple data. The one-dimensional thermal response prediction compared well with the thermocouple data for the leading edge "lip" region at the bolt location. In general, the test results yield confidence in the baseline seal design to prevent hot gas ingestion at the bracket and composite aeroshell structure interface.

Szalai, Christine E.↗

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

Analysis of ISS Plasma Interaction

To date, the International Space Station (ISS) has been one of the largest objects flown in lower earth orbit (LEO). The ISS utilizes high voltage solar arrays (160V) that are negatively grounded leading to pressurized elements that can float negatively with respect to the plasma. Because laboratory measurements indicate a dielectric breakdown potential difference of 80V, arcing could occur on the ISS structure. To overcome the possibility of arcing and clamp the potential of the structure, two Plasma Contactor Units (PCUs) were designed, built, and flown. Also a limited amount of measurements of the floating potential for the present ISS configuration were made by a Floating Potential Probe (FPP), indicating a minimum potential of 24 Volts at the measurement location. A predictive tool, the ISS Plasma Interaction Model (PIM) has been developed accounting for the solar array electron collection, solar array mast wire and effective conductive area on the structure. The model has been used for predictions of the present ISS configuration. The conductive area has been inferred based on available floating potential measurements. Analysis of FPP and PCU data indicated distribution of the conductive area along the Russian segment of the ISS structure. A significant input to PIM is the plasma environment. The International Reference Ionosphere (IRI 2001) was initially used to obtain plasma temperature and density values. However, IRI provides mean parameters, leading to difficulties in interpretation of on-orbit data, especially at eclipse exit where maximum charging can occur. This limits our predicative capability. Satellite and Incoherent Scatter Radar (ISR) data of plasma parameters have also been collected. Approximately 130,000 electron temperature (Te) and density (Ne) pairs for typical ISS eclipse exit conditions have been extracted from the reduced Langmuir probe data flown aboard the NASA DE-2 satellite. Additionally, another 18,000 Te and Ne pairs of ISR data from several radar locations around the globe were used to assure consistency of the satellite data. PIM predictions for ISS charging made with this data correlated very well with FPP data, indicating that the general physics of spacecraft charging with high voltage solar arrays have been captured. The predictions also provided the probabilities of occurrences for ISS charging. These probabilities give a numerical measure of the number of times when the ISS will approach or exceed the vehicle plasma hazard conditions for each configuration. In this paper we shall present the interaction mechanisms between the ISS and the surrounding plasma and give an overview of the PIM components. PIM predictions are compared with available data followed by a discussion of the variability of plasma parameters and the conductive area on the ISS. The ISS PIM will be further tested and verified as data from the Floating Potential Measurement Unit become available, and construction of the ISS continues.

Reddell, Brandon↗

Assessment of CFD Hypersonic Turbulent Heating Rates for Space Shuttle Orbiter

Turbulent CFD codes are assessed for the prediction of convective heat transfer rates at turbulent, hypersonic conditions. Algebraic turbulence models are used within the DPLR and LAURA CFD codes. The benchmark heat transfer rates are derived from thermocouple measurements of the Space Shuttle orbiter Discovery windward tiles during the STS-119 and STS-128 entries. The thermocouples were located underneath the reaction-cured glass coating on the thermal protection tiles. Boundary layer transition flight experiments conducted during both of those entries promoted turbulent flow at unusually high Mach numbers, with the present analysis considering Mach 10{15. Similar prior comparisons of CFD predictions directly to the flight temperature measurements were unsatisfactory, showing diverging trends between prediction and measurement for Mach numbers greater than 11. In the prior work, surface temperatures and convective heat transfer rates had been assumed to be in radiative equilibrium. The present work employs a one-dimensional time-accurate conduction analysis to relate measured temperatures to surface heat transfer rates, removing heat soak lag from the flight data, in order to better assess the predictive accuracy of the numerical models. The turbulent CFD shows good agreement for turbulent fuselage flow up to Mach 13. But on the wing in the wake of the boundary layer trip, the inclusion of tile conduction effects does not explain the prior observed discrepancy in trends between simulation and experiment; the flight heat transfer measurements are roughly constant over Mach 11-15, versus an increasing trend with Mach number from the CFD.

Wood, William A.↗

Real-Time Deposition Monitor for Ultrathin Conductive Films

A device has been developed that can be used for the real-time monitoring of ultrathin (2 or more) conductive films. The device responds in less than two microseconds, and can be used to monitor film depositions up to about 60 thick. Actual thickness monitoring capability will vary based on properties of the film being deposited. This is a single-use device, which, due to the very low device cost, can be disposable. Conventional quartz/crystal microbalance devices have proven inadequate to monitor the thickness of Pd films during deposition of ultrathin films for hydrogen sensor devices. When the deposited film is less than 100 , the QCM measurements are inadequate to allow monitoring of the ultrathin films being developed. Thus, an improved, high-sensitivity, real-time deposition monitor was needed to continue Pd film deposition development. The new deposition monitor utilizes a surface acoustic wave (SAW) device in a differential delay-line configuration to produce both a reference response and a response for the portion of the device on which the film is being deposited. Both responses are monitored simultaneously during deposition. The reference response remains unchanged, while the attenuation of the sensing path (where the film is being deposited) varies as the film thickness increases. This device utilizes the fact that on high-coupling piezoelectric substrates, the attenuation of an SAW undergoes a transition from low to very high, and back to low as the conductivity of a film on the device surface goes from nonconductive to highly conductive. Thus, the sensing path response starts with a low insertion loss, and as a conductive film is deposited, the film conductivity increases, causing the device insertion loss to increase dramatically (by up to 80 dB or more), and then with continued film thickness increases (and the corresponding conductivity increases), the device insertion loss goes back down to the low level at which it started. This provides a continuous, real-time monitoring of film deposition. For use with different films, the device would need to be calibrated to provide an understanding of how film thickness is related to film conductivity, as the device is responding primarily to conductivity effects (and not to mass loading effects) in this ultrathin film regime.

Hines, Jacqueline↗

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗

A Systems Engineering Approach to Architecture Development

Architecture development is conducted prior to system concept design when there is a need to determine the best-value mix of systems that works collectively in specific scenarios and time frames to accomplish a set of mission area objectives. While multiple architecture frameworks exist, they often require use of unique taxonomies and data structures. In contrast, this presentation characterizes architecture development using terminology widely understood within the systems engineering community. Using a notional civil space architecture example, it employs a multi-tier framework to describe the enterprise level architecture and illustrates how results of lower tier, mission area architectures integrate into the enterprise architecture. It also presents practices for conducting effective mission area architecture studies, including establishing the trade space, developing functions and metrics, evaluating the ability of potential design solutions to meet the required functions, and expediting study execution through the use of iterative design cycles.

Architecture Development↗

NASAs Evolvable Mars Campaign: Mars Moons Robotic Precursor

Human exploration missions to the moons of Mars are being considered within NASA's Evolvable Mars Campaign (EMC) as an intermediate step for eventual human exploration and pioneering of the surface of Mars. A range of mission architectures is being evaluated in which human crews would explore one or both moons for as little as 14 days or for as long as 500 days with a variety of orbital and surface habitation and mobility options being considered. Relatively little is known about the orbital, surface, or subsurface characteristics of either moon. This makes them interesting but challenging destinations for human exploration missions during which crewmembers must be able to effectively conduct scientific exploration without being exposed to undue risks due to radiation, dust, micrometeoroids, or other hazards. A robotic precursor mission to one or both moons will be required to provide data necessary for the design and operation of subsequent human systems and for the identification and prioritization of scientific exploration objectives. This paper identifies and discusses considerations for the design of such a precursor mission based on current human mission architectures. Objectives of a Mars' moon precursor in support of human missions are expected to include: 1) identifying hazards on the surface and the orbital environment at up to 50-km distant retrograde orbits; 2) collecting data on physical characteristics for planning of detailed human proximity and surface operations; 3) performing remote sensing and in situ science investigations to refine and focus future human scientific activities; and 4) prospecting for in situ resource utilization. These precursor objectives can be met through a combination or remote sensing (orbital) and in-situ (surface) measurements. Analysis of spacecraft downlink signals using radio science techniques would measure the moon's mass, mass distribution, and gravity field, which will be necessary to enable trajectory planning. Laser altimetry would precisely measure the moon's shape and improve the accuracy of radio science measurements. A telescopic imaging camera would map the moon at submeter resolution and photograph selected areas of interest at subcentimeter resolution and a visible and near-infrared (0.4-3.0 mm) imaging spectrograph would produce a global map of mineral composition variations at a resolution of tens of meters and maps of selected areas of interest at meter resolution. Additional remote sensing capabilities could include a thermal infrared imager (heat flow, thermal inertia, and grain size distributions), a gamma-ray and neutron detector (atomic composition), a ground-penetrating radar (internal structure), and a magnetometer and Langmuir probe (magnetic properties and plasma field). Once on the surface of Phobos or Deimos, necessary instrumentation would include a penetrometer (regolith compressive strength), a motion-imagery camera (to observe the penetrometer tests before, during, and after contact), a dust-adhesion witness plate and camera (dust levitation), a microimager (dust particle sizes and shapes), and an alpha-proton-X-ray, X-ray fluorescence, Mossbauer, or Raman spectrometer (atomic and mineral composition of surface materials) and an optional temperature probe (regolith thermal properties). A variety of robotic mission design options to enable both orbital and surface measurements are being considered that include fully integrated and modular approaches. In-situ measurements from at least one surface location would be required, with additional measurement locations possible through use of multiple landers, through propulsive relocation of a single lander, or through electromechanical surface translation by a walking or hopping lander vehicle, which could also serve to evaluate such mobility capabilities for subsequent human missions. Preliminary orbital analysis suggests that remote sensing would likely be performed while in a distant retrograde orbit around the target moon. Mission design options to enable characterization of both Mars’ moons in a single mission are also being studied.

Gernhardt, Michael L.↗

A Systems Engineering Approach to Architecture Development

Architecture development is often conducted prior to system concept design when there is a need to determine the best-value mix of systems that works collectively in specific scenarios and time frames to accomplish a set of mission area objectives. While multiple architecture frameworks exist, they often require use of unique taxonomies and data structures. In contrast, this paper characterizes architecture development using terminology widely understood within the systems engineering community. Using a notional civil space architecture example, it employs a multi-tier framework to describe the enterprise level architecture and illustrates how results of lower tier, mission area architectures integrate into the enterprise architecture. It also presents practices for conducting effective mission area architecture studies, including establishing the trade space, developing functions and metrics, evaluating the ability of potential design solutions to meet the required functions, and expediting study execution through the use of iterative design cycles.

Multi-tier Framework↗

A Systems Engineering Approach to Architecture Development

Architecture development is often conducted prior to system concept design when there is a need to determine the best-value mix of systems that works collectively in specific scenarios and time frames to accomplish a set of mission area objectives. While multiple architecture frameworks exist, they often require use of unique taxonomies and data structures. In contrast, this paper characterizes architecture development using terminology widely understood within the systems engineering community. Using a notional civil space architecture example, it employs a multi-tier framework to describe the enterprise level architecture and illustrates how results of lower tier, mission area architectures integrate into the enterprise architecture. It also presents practices for conducting effective mission area architecture studies, including establishing the trade space, developing functions and metrics, evaluating the ability of potential design solutions to meet the required functions, and expediting study execution through the use of iterative design cycles

Architecture Development↗

Hall Thruster Thermal Modeling and Test Data Correlation

HERMeS - Hall Effect Rocket with Magnetic Shielding. Developed through a joint effort by NASA/GRC and the Jet Propulsion Laboratory (JPL). Design goals: High power (12.5 kW) high Isp (3000 sec), high efficiency (> 60%), high throughput (10,000 kg), reduced plasma erosion and increased life (5 yrs) to support Asteroid Redirect Robotic Mission (ARRM). Further details see "Performance, Facility Pressure Effects and Stability Characterization Tests of NASAs HERMeS Thruster" by H. Kamhawi and team. Hall Thrusters (HT) inherently operate at elevated temperatures approx. 600 C (or more). Due to electric magnetic (E x B) fields used to ionize and accelerate propellant gas particles (i.e., plasma). Cooling is largely limited to radiation in vacuum environment.Thus the hardware components must withstand large start-up delta-T's. HT's are constructed of multiple materials; assorted metals, non-metals and ceramics for their required electrical and magnetic properties. To mitigate thermal stresses HT design must accommodate the differential thermal growth from a wide range of material Coef. of Thermal Expansion (CTEs). Prohibiting the use of some bolted/torqued interfaces.Commonly use spring loaded interfaces, particularly at the metal-to-ceramic interfaces to allow for slippage.However most component interfaces must also effectively conduct heat to the external surfaces for dissipation by radiation.Thus contact pressure and area are important.

Model↗

Coastal California Water Resources II: Utilizing NASA Earth Observations to Detect and Assess the Impacts of Estuarine Breach Events for Improved Coastal Wetland Monitoring and Management

Estuaries are dynamic environments that provide a host of vital ecosystem services. California’s Marine Life Protection Act protects such ecosystems by creating Marine Protected Areas. California has approximately 440,000 acres of estuarine habitats as well as 23 Estuarine Marine Protected Areas (EMPAs); thus, in situ data collection is often difficult due to time and resource constraints. This project used remote sensing to gather data that examined the health and dynamics of California EMPAs in order to supplement ground-based field measurements. Through the use of Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS), Sentinel-2 Multispectral Instrument (MSI) and Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), this project assessed mouth state, inundation extent, turbidity, temperature, and tidal measurements for observable estuaries. The Normalized Water Difference Index from Sentinel-2 MSI captured estuary mouth state and inundation extent. Landsat 8 OLI and Sentinel-2 MSI detected differences in water quality metrics that correlated to changes in estuary mouth state (i.e., open or closed). The team’s California Estuary Assessment (CEA) tool in Google Earth Engine was successful in analyzing estuary mouth state, inundation, and water quality. It was most effective when breach events were larger than 10 meters in resolution, water surface was smooth, and imagery was unimpeded by algae or sun glint. The CEA tool will allow the partners, the Ocean Protection Council, Central Coast Wetlands Group, Southern California Coastal Water Research Project, and University of California Los Angeles (UCLA) and Davis (UCD), to better understand estuary dynamics and more effectively conduct in situ estuary monitoring.

Sarah Payne↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

Cameron J. Hines↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗