Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flux adjustment”

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 253 records · Page 14

A model of stationary gravity wave breakdown with convective adjustment

A steady WKB model of gravity wave propagation including convective adjustment is used to investigate approximations used in various gravity-wave parameterization schemes. First, it is shown that estimates of the wave breaking height assuming a single horizontal wavenumber gravity wave can lead to errors if the topography is not sinusoidal. Second, the model results show that the assumption that wave growth ceases with the onset of convection or shear instability is an oversimplification. Since convection appears first over a very limited spatial region of the wave field, the wave is initially unaffected by turbulent mixing. However, when the convection zone spreads over a large portion of the wave field the amplitude is constrained. Estimates of the heat flux by breaking gravity waves are used to develop a simple parameterization of the vertical diffusion in terms of the Reynolds stress.

Schoeberl, Mark R.↗

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗

Contactless Growth of ZnSe Single Crystals by Physical Vapor Transport

ZnSe crystals were grown by self-seeded physical vapor transport (PVT) technique in the horizontal configuration. The source materials were heat treated by H2 reduction to remove the oxide followed by baking under dynamic vacuum to adjust the source composition toward that of congruent sublimation. Contactless growth of ZnSe single crystals have been performed consistently using three different source materials. The crystals grew away from the wall during the later stage of the growth with large (110) facets tend to align parallel to the gravity direction. The Scanning Electron Micrography (SEM) micrographs and the Atomic Force Microscope (AFM) images showed that large (110) terraces and steps dominate the as-grown facets. The measured residual gas pressures in the processed ampoules agree well among various source materials and the major components were CO and H2. No preferred growth direction was found. The one-dimensional diffusion model on the mass flux of a multi-species PVT system was employed to analyze the conditions for contactless growth. The calculated thermal profile for supersaturation is very close to the thermal profile measured inside the empty furnace bore in the region of contactless growth. The effects of convective flows in the vapor phase inside the ampoule on the growth processes are discussed.

Su, Ching-Hua↗

Growth Rate of Deep Convective System Cloud Shields: Satellite Observations and Km-Scale Radiative Convective Equilibrium Simulations

Deep convection gives rise to large upper level clouds that strongly interact with radiation and are important to the climate energy budget. From an object-oriented perspective, these individual deep cloud systems are characterized by a well depicted cloud shield life cycle, starting with small cloud extents that grow at varying rates before decaying and vanishing. A simple formulation of the growth rate of the cloud shield has been proposed that links together the growth rate on the convective part of the cloud, the mass flux of both the convective and stratiform parts of the cluster and a simple removal sink term (Elsaesser et al., 2022). In this presentation we first show using a suite of satellite observations (infrared from geostationary satellites, GPM radar, etc.) that the functional form of the proposed equation is well suited to quantify the shield growth rate. We then focus on RCE simulations, with deep cloud system objects post processed, to explore the relative role of each term of the growth rate budget. Three different models are used in the same RCEMIP-like configurations. The results show that the budget equation works very well for each model, although the time constants require model-dependent adjustments. We will further show in Vienna the commonalities and the specificities of each model.

Deep convection↗

A one- and two-layer model for estimating evapotranspiration with remotely sensed surface temperature and ground-based meteorological data over partial canopy cover

Surface-air temperature differences are commonly used in a bulk resistance equation for estimating sensible heat flux (H), which is inserted in the one-dimensional energy balance equation to solve for the latent heat flux (LE) as a residual. Serious discrepancies between estimated and measured LE have been observed for partial-canopy-cover conditions, which are mainly attributed to inappropriate estimates of H. To improve the estimates of H over sparse canopies, one- and two-layer resistance models that account for some of the factors causing poor agreement are developed. The utility of the two models is tested with remotely sensed and micrometeorological data for a furrowed cotton field with 20 percent cover and a dry soil surface. It is found that the one-layer model performs better than the two-layer model when a theoretical bluff-body correction for heat transfer is used instead of an empirical adjustment; otherwise, the two-layer model is better.

Kustas, William P.↗

Facility Designed and Built to Investigate the Combined Effects of Contaminant and Atomic Oxygen on a Light-Transmitting Surface

A need exists to investigate changes in the transparency of a light-transmitting surface during simultaneous exposure to a contaminant and an atomic oxygen (AO) flux. This mechanism may be responsible for the degradation of the light-transmitting surfaces of both photovoltaic cells and photodiodes currently in use on many low-Earth-orbit spacecraft. To address this need, researchers from the Electro-Physics Branch of the NASA Glenn Research Center at Lewis Field built such a test system within their AO beam facility. This facility produces an effective AO flux of 1.4 10(exp 16) atoms/sq cm/sec and contains a three-axis positioning system that provides the motion capability necessary for test operations. During testing, a target surface is held directly within the AO beam and close to two contaminant effusion cells. The effusion cells are shielded from the AO beam, and the outgassing contaminant is constrained to move across the target surface when heat is applied to either of the reservoirs. A light source is periodically moved over the target surface, and the transmitted light intensity is checked with a photodiode located below the target. This light source is also periodically checked with a separate photodiode, which is protected from contamination and AO exposure, to allow adjustments necessary to maintain a consistent light intensity.

Sechkar, Edward A.↗

Lunar Neutral Exposphere Properties from Pickup Ion Analysis

Composition and structure of neutral constituents in the lunar exosphere can be determined through measurements of phase space distributions of pickup ions borne from the exosphere [1]. An essential point made in an early study [ 1 ] and inferred by recent pickup ion measurements [2, 3] is that much lower neutral exosphere densities can be derived from ion mass spectrometer measurements of pickup ions than can be determined by conventional neutral mass spectrometers or remote sensing instruments. One approach for deriving properties of neutral exospheric source gasses is to first compare observed ion spectra with pickup ion model phase space distributions. Neutral exosphere properties are then inferred by adjusting exosphere model parameters to obtain the best fit between the resulting model pickup ion distributions and the observed ion spectra. Adopting this path, we obtain ion distributions from a new general pickup ion model, an extension of a simpler analytic description obtained from the Vlasov equation with an ion source [4]. In turn, the ion source is formed from a three-dimensional exospheric density distribution, which can range from the classical Chamberlain type distribution to one with variable exobase temperatures and nonthermal constituents as well as those empirically derived. The initial stage of this approach uses the Moon's known neutral He and Na exospheres to deriv e He+ and Na+ pickup ion exospheres, including their phase space distributions, densities and fluxes. The neutral exospheres used are those based on existing models and remote sensing studies. As mentioned, future ion measurements can be used to constrain the pickup ion model and subsequently improve the neutral exosphere descriptions. The pickup ion model is also used to estimate the exosphere sources of recently observed pickup ions on KAGUYA [3]. Future missions carrying ion spectrometers (e.g., ARTEMIS) will be able to study the lunar neutral exosphere with great sensitivity, yielding the necessary ion velocity spectra needed to further analysis of parent neutral exosphere properties.

Hartle, R. E.↗

Can Terrestrial Microbes Grow on Mars?

The theme for AbSciCon 2012 is "Exploring Life: Past and Present, Near and Far." The conference will address our current understanding of life - from processes at the molecular level to those which operate at planetary scales. Studying these aspects of life on Earth provides an essential platform from which to examine the potential for life on other worlds, both within our solar system and beyond. Mars exhibits a variety of extreme environments characterized by high UV and ionizing radiation flux, low pressure anoxic atmosphere, scarce or absent liquid water, extreme low temperatures, etc. The ability of terrestrial microorganisms to survive and adapt to the Mars environment has profound implications for astrobiology, planetary protection, and Mars life detection missions. At the NASA Ames Synthetic Biology Initiative, we believe that synthetic biology has the potential to revolutionize human space exploration. As such, the initiative is dedicated to applying the tools and techniques of synthetic biology to space exploration and astrobiology. Biological solutions will be invaluable for space exploration because they are not resource intensive, and they are versatile and self-renewing. An understanding of how to work with DNA in an unfavorable environment is paramount to utilizing biological tools on space missions. Furthermore, the ability to adjust life to the parameters of Mars is vital both to discovering what life on Mars might look like, and to using biological tools under such conditions. As a first step, we need an energy-efficient, low cost means of transporting, storing, and protecting genomic DNA, DNA parts, and whole microbial strains. Our goal is to develop and demonstrate viable and superior alternatives to standard DNA storage methods, which can be optimized to the conditions of space exploration, using synthetic biology as a tool. This includes protocols and kit designs for easy and repeatable DNA and strain recovery from protective storage conditions. We are constructing newly engineered genetic parts for different valuable host organisms, designed to increased long-term survival and functional retention. These methods should be applied for DNA and strain storage and transportation. In parallel, we seek inspiration from natural organisms that have developed means for survival in extreme environmental conditions. We are utilizing novel techniques for analysis of lipid biomarkers in the Antarctic Dry Valleys in order to identify resident microbes in the Antarctic soil and permafrost, as well as biomarker fossils of organisms that survived in the valleys in ages past. Through the identification of these life forms, we hope to understand and draw on new biological tools and strategies for synthetic biological applications on Mars.

Rothschild, Lynn↗

Black Carbon and Precipitation: an Energetics Perspective

Black carbon aerosols (BC) influence precipitation through a range of processes. The climate response to the presence of BC is however highly dependent on its vertical distribution. Here, we analyze the changes in the energy budget and precipitation impacts of adding a layer of BC at a range of altitudes in two independent global climate models. The models are run with atmosphere‐only and slab ocean model setup to analyze both fast and slow responses, respectively. Globally, precipitation changes are tightly coupled to the energy budget. We decompose the precipitation change into contributions from absorption of solar radiation, atmospheric longwave radiative cooling, and sensible heat flux at the surface. We find that for atmosphere‐only simulations, BC rapidly suppress precipitation, independent of altitude, mainly because of strong atmospheric absorption. This reduction is offset by increased atmospheric radiative longwave cooling and reduced sensible heat flux at the surface, but not of sufficient magnitude to prevent reduced precipitation. On longer timescales, when the surface temperature is allowed to respond, we find that the precipitation increase associated with surface warming can compensate for the initial reduction, particularly for BC in the lower atmosphere. Even though the underlying processes are strikingly similar in the two models, the resulting change in precipitation and temperature by BC differ quite substantially.

black carbon aerosols↗

Regional Deformation Studies with GRACE and GPS

GRACE data indicate large seasonal variations in gravity that have been shown to be to be related to climate-driven fluxes of surface water. Seasonal redistribution of surface mass deforms the Earth, and our previous study using GRACE data demonstrate that annual radial deformations of +/-13 mm in the region of Amazon River Basin were observed by both GRACE and ten GPS sites in the region. For the GRACE determinations, we estimate in a least-squares solution for each Stokes coefficient parameters that represent the amplitudes of the annual variation. We then filter these parameters based on a statistical test that uses the scatter of the postfit residuals. We demonstrate by comparison to the GPS amplitudes that this method is more accurate, for this region, than Gaussian smoothing. Our model for the temporal behavior of the gravity coefficients includes a rate term, and although the time series are noisy, the glacial isostatic adjustment signal over Hudson s Bay can be observed. .

Davis, J. L.↗

Characteristics of Whipple Shield Performance in the Shatter Regime

Between the onset of projectile fragmentation and the assumption of rear wall failure due to an impulsive load, multi-wall ballistic limit equations are linearly interpolated to provide reasonable yet conservative predictions of perforation thresholds with conveniently simple mathematics. Although low velocity and hypervelocity regime predictions are based on analytical expressions, there is no such scientific foundation for predictions in the intermediate (or shatter) regime. As the debris flux in low earth orbit (LEO) becomes increasingly dominated by manmade pollution, the profile of micrometeoroid and orbital debris (MMOD) risk shifts continually towards lower velocities. For the International Space Station (ISS), encounter velocities below 7 km/s now constitute approximately 50% of the penetration risk. Considering that the transition velocity from shatter to hypervelocity impact regimes described by common ballistic limit equations (e.g. new non-optimum Whipple shield equation [1]) occurs at 7 km/s, 50% of station risk is now calculated based on failure limit equations with little analytical foundation. To investigate projectile and shield behavior for impact conditions leading to projectile fragmentation and melt, a series of hypervelocity impact tests have been performed on aluminum Whipple shields. In the experiments projectile diameter, bumper thickness, and shield spacing were kept constant, while rear wall thickness was adjusted to determine spallation and perforation limits at various impact velocities and angles. The results, shown in Figure 1 for normal and 45 impacts, demonstrated behavior that was not sufficiently described by the simplified linear interpolation of the NNO equation (also shown in Figure 1). Hopkins et al. [2] investigated the performance of a nominally-identical aluminum Whipple shield, identifying the effects of phase change in the shatter regime. The results (conceptually represented in Figure 2) were found to agree well with those obtained in this study at normal incidence, suggesting that shielding performance in the shatter regime could be well described by considering more complex phase conditions than currently implemented in most BLEs. Furthermore, evidence of these phase effects were found in the oblique test results, providing the basis for an empirical description of these effects that can be applied in MMOD risk assessment software. In this paper, results of the impact experiments are presented, and characteristics of target damage are evaluated. A comparison of intermediate velocity impact failure mechanisms in current BLEs are discussed and compared to the findings of the experimental study. Risk assessment calculations have been made on a simplified structure using currently implemented penetration equations and predicted limits from the experimental program, and the variation in perceived mission risk is discussed. It was found that ballistic limit curves that explicitly incorporated phase change effects within the intermediate regime lead to a decrease in predicted MMOD risk for ISS-representative orbits. When considered for all Whipple-based shielding configurations onboard the ISS, intermediate phase change effects could lead to significant variations in predicted mission risk.

Ryan, Shannon↗

The Martian atmospheric planetary boundary layer stability, fluxes, spectra, and similarity

This is the first analysis of the high frequency data from the Viking lander and spectra of wind, in the Martian atmospheric surface layer, along with the diurnal variation of the height of the mixed surface layer, are calculated for the first time for Mars. Heat and momentum fluxes, stability, and z(sub O) are estimated for early spring, from a surface temperature model and from Viking Lander 2 temperatures and winds at 44 deg N, using Monin-Obukhov similarity theory. The afternoon maximum height of the mixed layer for these seasons and conditions is estimated to lie between 3.6 and 9.2 km. Estimations of this height is of primary importance to all models of the boundary layer and Martian General Circulation Models (GCM's). Model spectra for two measuring heights and three surface roughnesses are calculated using the depth of the mixed layer, and the surface layer parameters and flow distortion by the lander is also taken into account. These experiments indicate that z(sub O), probably lies between 1.0 and 3.0 cm, and most likely is closer to 1.0 cm. The spectra are adjusted to simulate aliasing and high frequency rolloff, the latter caused both by the sensor response and the large Kolmogorov length on Mars. Since the spectral models depend on the surface parameters, including the estimated surface temperature, their agreement with the calculated spectra indicates that the surface layer estimates are self consistent. This agreement is especially noteworthy in that the inertial subrange is virtually absent in the Martian atmosphere at this height, due to the large Kolmogorov length scale. These analyses extend the range of applicability of terrestrial results and demonstrate that it is possible to estimate the effects of severe aliasing of wind measurements, to produce a models which agree well with the measured spectra. The results show that similarity theory developed for Earth applies to Mars, and that the spectral models are universal.

Tillman, James E.↗

Improved Spectral and Angular Characterization of Libya-4 and Dome-C for consistent MODIS and VIIRS Radiometric Scaling

The NASA CERES observed SW and LW fluxes are utilized by the climate community for monitoring the Earth’s energy imbalance and for climate model validation. To facilitate seamless flux and cloud properties across MODIS, VIIRS, and geostationary (GEO) imagers the CERES project intercalibrates the imagers directly comparing the coincident ray-matched analogous channel radiances. For CERES edition 4 products, the MODIS, VIIRS and GEO imagers were radiometrically scaled to the Aqua-MODIS C5 calibration reference. There are no direct comparison opportunities between NPP and NOAA20 VIIRS, since they are positioned a half an orbit apart. The previous usage of Aqua-MODIS as a transfer radiometer is coming to an end, since the Aqua orbit is slowly drifting towards the terminator and will be deorbited in 2026. The CERES project will utilize the Libya-4 and Dome-C invariant targets to radiometrically scale between MODIS and VIIRS reflective solar bands. The Libya-4 and Dome-C Earth invariant targets will be characterized by the repeat cycle angular configuration. The target spectral band adjustment factors (SBAF) will be derived using EMIT and DESIS hyper-spectral observations onboard the ISS and will be compared with the existing SCIAMACHY, GOME-2, and Hyperion SBAFs. The DESIS and EMIT SBAF approach can incorporate the future CLARREO CPF SI traceable hyper-spectral measurements allowing the targets to be referenced to an absolute calibration reference. Utilizing the BRDF and atmospheric corrected characterization of one sensor with another will provide the scaling factors. Consistent scaling factors between the two invariant target validates the method.

David R. Doelling↗

Constructing the Coronal Magnetic Field: by Correlating Parameterized Magnetic Field Lines with Observed Coronal Plasma Structures

The reconstruction of the coronal magnetic field is carried out using a perturbation procedure. A set of magnetic field lines generated from magnetogram data is parameterized and then deformed by varying the parameterized values. The coronal fluxtubes associated with this field are adjusted until the correlation between the field lines and the observed coronal loops is maximized. A mathematical formulation is described which ensures (1) that the normal component of the photospheric field remains unchanged, (2) that the field is given in the entire corona, (3) that the field remains divergence free, and (4) that electrical currents are introduced into the field. It is demonstrated that a simple radial parameterization of a potential field, comprising a radial stretching of the field, can provide a match for a simple bipolar active region, AR 7999, which crossed the central meridian on 1996 Nov 26. At a coronal height of 30 km, the resulting magnetic field is a non-force free magnetic field with the maximum Lorentz force being on the order of 2.6 x 10(exp -9) dyn resulting from an electric current density of $0.13 mu A/ sq m. This scheme is an important tool in generating a magnetic field solution consistent with the coronal flux tube observations and the observed photospheric magnetic field.

Gary, G. A.↗

Evaporation/SST Sensitivity Over the Tropical Oceans During ENSO Events as Estimated from the da Silva, Young, Levitus Surface Marine Data Set

The da Silva, Young and Levitus Surface Marine Atlas, based on observations from the Comprehensive Ocean Atmosphere Data Set (COADS) Release 1, has been used to investigate the relationship between evaporation and sea-surface temperature (SST) over the global oceans. For the period 1950 to 1987 SST, surface latent heat flux, and other related variables have been filtered to minimize data uncertainties and to focus upon interannual variations associated with warm (El Nino) and cold (La Nina) ENSO events. Compositing procedures have enabled identification of systematic variations in latent heat fluxes accompanying these events and the relationship to spatial anomalies in ocean surface wind speed and humidity. The evaporation response associated with ENSO sea surface temperature (SST) variability is systematic in nature and composed of offsetting contributions from the surface wind and humidity variations. During warm events exceeding 1.0 S.D. delta SST, increases in the surface humidity deficit, delta(qs-qa), between the surface and 2m height dominate regions of positive SST anomalies and lead to increases in evaporation of almost 2 Wm (exp -2) at deltaSST = 0.23 K. Despite the increases in specific humidity, relative humidity decreases slightly in regions of elevated SSTs. For the most part, variations in wind speed are consistent with previous investigations. Weakening of the equatorial easterlies (and generation of westerlies) between 160 degrees E and 140 degrees W dominates during the early phases of warm events. Elevated wind speeds in adjacent subtropical regions and in the eastern equatorial Pacific subsequently develop too. The net contribution of these winds, which reflect adjustments in Hadley and Walker circulation components is toward reduced evaporation. Results for cold periods are approximately similar, but opposite in sign to warm events, though evidence of different temporal evolution is noted.

Robertson, F. R.↗

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↗