Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase association”

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 451 records · Page 25

Stream structures in the outer heliosphere

Nineteen magnetic clouds are identified in the years from 1978 through 1982 and studied by the superimposed epoch method. The magnetic fluctuations, density, and temperature are enhanced ahead of the clouds preceded by shocks. Strong magnetic field intensities and low proton temperatures are observed in the clouds. A relatively large (2.5%) decrease in cosmic ray intensity is caused by the turbulent sheath behind an interplanetary shock ahead of a magnetic cloud. Only a small (0.5%) decrease in intensity is associated with the magnetic cloud itself. Magnetic clouds can produce geomagnetic activity with a decrease in the Dst index of the order of 100 gammas. The magnitude of the change in Dst index for the case when southward fields arrive first is comparable to that for the case of northward fields first, and the phase is such that geomagnetic activity is associated with the southward fields.

Whang, Y. C.↗

Magnetic clouds, geomagnetic disturbances, and cosmic ray decreases

Nineteen magnetic clouds are identified in the years from 1978 through 1982 and studied using superposed epoch analysis. The magnetic-field intensity, proton density and proton temperature are enhanced ahead of magnetic clouds that are preceded by shock, while strong magnetic-field intensity and low proton temperature are observed within the clouds. A relatively large (about 2.5 percent) decrease in cosmic-ray intensity is associated with magnetic clouds that are preceded by a shock, perhaps caused by the turbulent sheath behind an interplanetary shock ahead of the magnetic cloud, whereas only a small (0.5 percent) decrease in intensity is associated with the magnetic cloud itself. Magnetic clouds can produce geomagnetic activity with a decrease in Dst index of the order 100 gamma. The magnitude of the change in the Dst index for the case when southward fields arrive first is comparable to that for the case when northward fields arrive first, and the phase is such that geomagnetic activity is associated with southward fields.

Zhang, G.↗

Macroscopic aspects of interfacial reactions

The extent of interdiffusion and formation of new phases is determined by the constitution diagram of the alloy system, the interdiffusion coefficients of the phases present, and the thermal conditions (temperature and time) associated with the bonding process and/or subsequent use of the bonded structure. In many instance, the kinetics of interdiffusion and phase formation can be predicted from known parameters using numerical methods and computer techniques. Predictions are compared with experimentally determined parameters for a variety of metallurgical alloy systems.

Heckel, R. W.↗

Dispersive Elements for Enhanced Laser Gyroscopy and Cavity Stabilization

We analyze the effect of a highly dispersive element placed inside a modulated optical cavity on the frequency and amplitude of the modulation to determine the conditions for cavity self-stabilization and enhanced gyroscopic sensitivity. We find an enhancement in the sensitivity of a laser gyroscope to rotation for normal dispersion, while anomalous dispersion can be used to self-stabilize an optical cavity. Our results indicate that atomic media, even coherent superpositions in multilevel atoms, are of limited use for these applications, because the amplitude and phase filters work against one another, i.e., decreasing the modulation frequency increases its amplitude and vice-versa. On the other hand, for optical resonators the dispersion reversal associated with critical coupling enables the amplitude and phase filters to work together. We find that for over-coupled resonators, the absorption and normal dispersion on-resonance increase the contrast and frequency of the beat-note, respectively, resulting in a substantial enhancement of the gyroscopic response. Under-coupled resonators can be used to stabilize the frequency of a laser cavity, but result in a concomitant increase in amplitude fluctuations. As a more ideal solution we propose the use of a variety of coupled-resonator-induced transparency that is accompanied by anomalous dispersion.

Smith, David D.↗

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from the three level-2 TRMM standard profile algorithms for the 1998 annual cycle. Results show that diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are producing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas, rainfall peaks during daytime; however, important exceptions are noted. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening-early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall, there is a secondary MLA maximum coexisting with the LE-EM maximum at latitudes dominated by large-scale convergence and deep convection. Analogously, there is a preponderancy for an LE EM maximum over land coexisting with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large-scale circulation. Neither of the secondary maxima exhibit phase behavior that can be considered semidiurnal in nature. Diurnal rainfall variability over the ocean associated with large-scale convection is clearly an integral component of the general circulation. Phase analysis reveals differences in regional and seasonal features of the diurnal cycle, indicating that underlying forcing mechanisms differ from place to place. This is underscored by the appearance of secondary ocean maxima in the presence of large-scale convection, along with other important features. Among these, there are clear-cut differences between the diurnal variability of seasonal rainfall over the mid-Pacific and Indian Ocean Basins. The mid-Pacific exhibits double maxima in spring and winter but only LE-EM maxima in summer and autumn, while the Indian Ocean exhibits double maxima in spring and summer and only an LE-EM maximum in autumn and winter. There are also evident daytime maxima within the major large-scale marine stratocumulus regions off the west coasts of continents. The study concludes with a discussion concerning how the observational evidence either supports or repudiates possible forcing mechanisms that have been suggested to explain diurnal rainfall variability.

Yang, Song↗

Influence of the continuous and dispersed phases on the symmetry of a gas turbine air-blast atomizer

Current trends in liquid-fueled practical combustion systems are leaving less tolerance for fuel injection deficiencies such as poor spray field symmetry. The present paper evaluates the symmetry of the flowfield produced by a practical airblast atomizer. Specifically, the influence of both the continuous phase and dispersed phase on the spray field symmetry is assessed. In the present case, asymmetry in volume flux is associated principally with disparities in the injection of the dispersed phase, which is manifested by a maldistribution of larger drops. Asymmetries observed in the continuous phase without the dispersed phase are reduced in magnitude by the presence of the dispersed phase, but still contribute to asymmetry in radial spread of the dispersed phase.

Mcdonell, V. G.↗

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

UV spectroscopy of Z Chamaeleontis. I - Time dependent dips in superoutburst

Extensive IUE observations of the dwarf nova Z Cha during the 1987 April superoutburst and IUE-Exosat observations during the 1985 July superoutburst are presented. The UV light curve shows two dips when folded on the orbital period. One dip, at orbital phase 0.8 becomes shallower as the superoutburst progresses and can be associated with decreasing mass transfer rate from the secondary star. The other dip, at orbital phase 0.15, appears after the development of the superhump and some days after the occurrence of the largest dip at phase 0.8. The continuum flux distribution during superoutbursts is fainter and redder than in low-inclination dwarf novae during superoutbursts. This is interpreted in terms of the extended vertical disk structure which occults the hot inner parts of the disk with the development of a 'cool' bulge on the edge of the disk at orbital phase 0.75. Details of the behaviour of the UV emission lines during eclipse and away from eclipse are discussed.

Harlaftis, E. T.↗

Satellite-based Estimates of Ambient Air Pollution and Global Variations in Childhood Asthma Prevalence

Background: The effect of ambient air pollution on global variations and trends in asthma prevalence is unclear. Objectives: Our goal was to investigate community-level associations between asthma prevalence data from the International Study of Asthma and Allergies in Childhood (ISAAC) and satellite-based estimates of particulate matter with aerodynamic diameter < 2.5 microm (PM2.5) and nitrogen dioxide (NO2), and modelled estimates of ozone. Methods: We assigned satellite-based estimates of PM2.5 and NO2 at a spatial resolution of 0.1deg × 0.1deg and modeled estimates of ozone at a resolution of 1deg × 1deg to 183 ISAAC centers. We used center-level prevalence of severe asthma as the outcome and multilevel models to adjust for gross national income (GNI) and center- and country-level sex, climate, and population density. We examined associations (adjusting for GNI) between air pollution and asthma prevalence over time in centers with data from ISAAC Phase One (mid-1900s) and Phase Three (2001-2003). Results: For the 13- to 14-year age group (128 centers in 28 countries), the estimated average within-country change in center-level asthma prevalence per 100 children per 10% increase in center-level PM2.5 and NO2 was -0.043 [95% confidence interval (CI): -0.139, 0.053] and 0.017 (95% CI: -0.030, 0.064) respectively. For ozone the estimated change in prevalence per parts per billion by volume was -0.116 (95% CI: -0.234, 0.001). Equivalent results for the 6- to 7-year age group (83 centers in 20 countries), though slightly different, were not significantly positive. For the 13- to 14-year age group, change in center-level asthma prevalence over time per 100 children per 10% increase in PM2.5 from Phase One to Phase Three was -0.139 (95% CI: -0.347, 0.068). The corresponding association with ozone (per ppbV) was -0.171 (95% CI: -0.275, -0.067). Conclusion: In contrast to reports from within-community studies of individuals exposed to traffic pollution, we did not find evidence of a positive association between ambient air pollution and asthma prevalence as measured at the community level.

Air pollution↗

Aeroacoustic Computations of a Transonic Truss-Braced Wing Aircraft: Part 2 – Acoustic Signature and Noise Source Identification

High-fidelity, time-dependent simulations of a Boeing-designed, transonic, truss-braced-wing aircraft in cruise (clean) and landing configurations are leveraged to generate synthetic microphone-phased-array data for airframe noise prediction and assessment. These data sets are used to compute source localization (beamform) maps to determine the location and strength of primary and secondary airframe noise sources associated with this unique configuration. The synthetic phased-array implementation mimics the setup of a flight test. As this study is ongoing, preliminary integrated far-field spectra for the cruise configuration obtained at multiple spatial resolutions revealed significant tonal content that lacked convergence with increased resolution. The origin of several of these tones and their unusual convergence behavior was traced to the larger-than-normal trailing-edge thickness of the “as-tested” cruise model being simulated. Reducing the trailing-edge thickness to more realistic values eliminated most of the tones at low to moderate frequencies and improved spectrum convergence significantly. Applying lessons learned from the cruise simulations, several modifications to the geometry of the landing configuration were made and are described in this work. Results from permeable and solid Ffowcs-Williams and Hawkings surfaces at two different spatial resolutions (coarse and medium) are used to illustrate the major noise sources and determine convergence of the CLEAN integrated noise levels for the entire aircraft as well as major subcomponents. We demonstrate that the low-frequency content of the far-field spectrum is dominated by noise generated from the main landing gear, while the medium- and high-frequency content is dominated by the wing-leading-edge Krueger flaps. Since analysis of the acoustic maps for the landing configuration revealed several clusters of multiple sources along the wing leading edge, “high resolution” processing of the array data was used to distinguish more accurately the locations of sources.

Transonic Truss-Braced Wing↗

Performing spectroscopic and specific heat studies of improper ferroelectrics

The results of infrared measurements on Ni-Br, Cu-Cl, and Fe-I boracite improper ferroelectrics and far infrared measurements of Ni-Br boracite are presented. The boracites have the general formula X3B7O3Y, where X = divalent metal and Y = halogen. They undergo a first order phase transition from a high temperature paraelectric phase with cubic symmetry to a ferroelectric phase with orthorhombic symmetry. The boracites are "improper ferroelectrics" since the spontaneous polarization is not the primary order parameter in the cubic-orthorhombic phase transition. Current understanding of these materials is that the primary order parameter is associated with a doubly degenerate zone-boundary phonon in the cubic phase. The degenerate critical modes become homogeneous and split into the A sub 1 and A sub 2 modes in the orthorhombic phase, doubling the volume of the primitive cell. An harmonic coupling between the softing A sub 1 and a low frequency A sub 1 optic mode induces a spontaneous polarization as a secondary effect in the ferroelectric phase. This secondary non-critical nature of the ferroelectric mode earns these materials the "improper" title and is responsible for their unique properties and high figure of merit in detector use.

Coleman, L. B.↗

Variation in Sleep Duration and Circadian Phase by Duty Start Time Among Short-Haul Commercial Airline Pilots

Prior studies have confirmed that commercial airline pilots experience circadian phase shifts and short sleep duration following travel with layovers in different time zones. Few studies have examined the impact of early and late starts on the sleep and circadian phase of airline pilots who return to their domicile after each duty period. We recruited 44 pilots (4 female) from a short-haul commercial airline to participate in a study examining sleep and circadian phase over four duty schedules (baseline, early starts, mid-day starts, late starts). Each duty schedule was five days long, separated by three rest days. Participants completed the rosters in the same order. Sleep outcomes were estimated using wrist-borne actigraphy (Actiware Software, Respironics, Bend, OR) and daily sleep diaries. Thirteen participants volunteered to collect urine samples for the assessment of 6-sulfatoxymelatonin (aMT6s). Urine samples were collected in four hourly bins during the day and eight-hourly bins during sleep episodes, for 24 hours immediately following each experimental duty schedule. The aMT6s results were fit to a cosine in order to obtain the acrophase to estimate circadian phase. Univariate statistics were calculated for acrophase changes, schedule start times and sleep times. All statistical analyses were computed using SAS software (Cary, IN).The mean duty start times varied as expected (baseline 10:17 [ 3:50], early starts 5:24 [ 0:30], mid-day starts 13:52 [ 1:20], late starts 16:33 [ 1:33]). Actigraphy-derived sleep duration varied between schedule types, with the shortest average sleep durations occurring during the early starts and night duty. Mean circadian phase was similar during each schedule type (baseline 26:14 [ 3:22], early starts 25:29 [2:13], mid-day starts 26:20 [ 3:16], late starts 25:49 [ 4:28]), but there were wide inter-individual differences in phase shifting from the beginning to the end of the study, with phase shifts ranging from a 5.98 hour phase advance to an 11.34 hour phase delay. Our preliminary findings suggest that early and late starts are associated with reduced sleep duration. The dispersion in inter-individual differences in circadian phase across shifting duty schedules should be further evaluated to determine how schedule-induced phase shifts may affect operational performance.

Flynn-Evans, Erin↗

Interaction of Light and Ethylene on Stem Gravitropism

The major objective of this study was to evaluate light-regulated ethylene production during gravitropic bending in etiolated pea stems. Previous investigations indicated that ethylene production increases after gravistimulation and is associated with the later (counter-reactive) phase of bending. Additionally, changes in the counter-reaction and locus of curvature during gravitropism are greatly influenced by red light and ethylene production. Ethylene production may be regulated by the levels of available precursor (1-aminocyclopropane-l-carboxylic acid, ACC) via its synthesis, conjugation to malonyl-ACC or glutamyl-ACC, or oxidation to ethylene. The regulation of ethylene production by quantifying ACC and conjugated ACC levels in gravistimulated pea stemswas examined. Also measured was the changes in protein and enzyme activity associated with gravitropic curvature by electrophoretic and spectrophotometric techniques. An image analysis system was used to visualize and quantify enzymatic activity and transcriptional products in gravistimulated and red-light treated etiolated pea stem tissues.

Harrison, Marcia A.↗

Three-Phased Wake Vortex Decay

A detailed parametric study is conducted that examines vortex decay within turbulent and stratified atmospheres. The study uses a large eddy simulation model to simulate the out-of-ground effect behavior of wake vortices due to their interaction with atmospheric turbulence and thermal stratification. This paper presents results from a parametric investigation and suggests improvements for existing fast-time wake prediction models. This paper also describes a three-phased decay for wake vortices. The third phase is characterized by a relatively slow rate of circulation decay, and is associated with the ringvortex stage that occurs following vortex linking. The three-phased decay is most prevalent for wakes imbedded within environments having low-turbulence and near-neutral stratification.

Proctor, Fred H.↗

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↗

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↗