Results for “Droplet effective radius” — Page 5

NASA NTRS1994

Simulations of the effects of water vapor, cloud liquid water, and ice on AMSU moisture channel brightness temperatures

Radiative transfer simulations are performed to determine how water vapor and nonprecipitating cloud liquid water and ice particles within typical midlatitude atmospheres affect brightness temperatures T(sub B)'s of moisture sounding channels used in the Advanced Microwave Sounding Unit (AMSU) and AMSU-like instruments. The purpose is to promote a general understanding of passive top-of-atmosphere T(sub B)'s for window frequencies at 23.8, 89.0, and 157.0 GHz, and water vapor frequencies at 176.31, 180.31, and 182.31 GHz by documenting specific examples. This is accomplished through detailed analyses of T(sub B)'s for idealized atmospheres, mostly representing temperate conditions over land. Cloud effects are considered in terms of five basic properties: droplet size distribution, phase, liquid or ice water content, altitude, and thickness. Effects on T(sub B) of changing surface emissivity also are addressed. The brightness temperature contribution functions are presented as an aid to physically interpreting AMSU T(sub B)'s. Both liquid and ice clouds impact the T(sub B)'s in a variety of ways. The T(sub B)'s at 23.8 and 89 GHz are more strongly affected by altostratus liquid clouds than by cirrus clouds for equivalent water paths. In contrast, channels near 157 and 183 GHz are more strongly affected by ice clouds. Higher clouds have a greater impact on 157- and 183-GHz T(sub B)'s than do lower clouds. Clouds depress T(sub B)'s of the higher-frequency channels by suppressing, but not necessarily obscuring, radiance contributions from below. Thus, T(sub B)'s are less closely associated with cloud-top temperatures than are IR radiometric temperatures. Water vapor alone accounts for up to 89% of the total attenuation by a midtropospheric liquid cloud for channels near 183 GHz. The Rayleigh approximation is found to be adequate for typical droplet size distributions; however, Mie scattering effects from liquid droplets become important for droplet size distribution functions with modal radii greater than 20 micrometers near 157 and 183 GHz, and greater than 30-40 micrometers at 89 GHz. This is due mainly to the relatively small concentrations of droplets much larger than the mode radius. Orographic clouds and tropical cumuli have been observed to contain droplet size distributions with mode radii in the 30-40 micrometers range. Thus, as new instruments bridge the gap between microwave and infrared to frequencies even higher than 183 GHz, radiative transfer modelers are cautioned to explicitly address scattering characteristics of such clouds.

Muller, Bradley M.↗

NASA NTRS1995

Pinatubo Aerosol Global-to-Micro-Scale Evolution: A Unified Picture From Space, Air, and Ground Measurements

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations In showing a dominant particle composition of H2SO4-H2O mature with H2SO4 weight traction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all attitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained In high latitude winters and at the tropical tropopause. Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra for comparison to direct measurements by impactors and optical counters, All three techniques paint a generally consistent picture of the evolution of R(sub eff) the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the pre-eruption value of 0.1 to 0.2 microns because both small (r less than 0.2 microns) and large (r greater than 0.6 microns) particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns, reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 does optical depth spectra begin to show a clear return to the preeruption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- and 1.0-kilometer optical depth derived from the AVHRR and SAGE satellite measurements. It also sets a context for evaluating remaining uncertainties in each of these satellite data products. We also show how the effects of wavelength-dependent refractive index on backscatter spectra can influence particle sizes retrieved from multiwavelength lidar measurements.

Russell, Philip B.↗

NASA NTRS1995

A Unified Picture of Pinatubo Aerosol Global-to Micro-Scale Evolution, From Space, Air, and Ground Measurements

We combine a variety of measurements to develop a composite picture of the post-Pinatubo aerosol and assess the consistency and uncertainties of the measurement and retrieval techniques. Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements are in accord with theoretical calculations in showing a dominant particle composition of H2SO4-H2O mixture, with H2SO4 weight fraction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are: (1) the presence of volcanic ash at all altitudes initially and in a layer just above the tropopause until at least March 1992, and (2) much smaller H2SO4 weight fractions at the low temperatures attained in high latitude winters and at the tropical tropopause, Laboratory spectroscopy and theoretical calculations yield wavelength- and temperature-dependent refractive indices for the dominant H2SO4-H2O droplets. These in turn permit derivation of particle size spectra from measured optical depth spectra, for comparison to direct measurements by impactors and optical counters. All three techniques paint a generally consistent picture of the evolution of R(sub eff), the effective, or area-weighted, particle radius. In the first month after the eruption, although particle numbers increased by orders of magnitude, R(sub eff) was similar to the preemption value of 0.1 to 0.2 microns, because both small (r less than 0.2 microns) and large (r greater than 0.6 micron particles increased in number. Over the next 3-6 months, R(sub eff) increased to about 0.5 microns reflecting particle growth through condensation and coagulation. In general, R(sub eff) continued to increase for about a year after the eruption. Extinction spectra computed from in situ size distribution measurements are consistent with optical depth measurements, which show spectra with maxima initially at wavelengths of 0.42 microns or less, and thereafter progressively increasing to between 0.78 and 1 micron. Not until 1993 do optical depth spectra begin to show a clear return to the preemption signature of maximizing at the shortest visible wavelengths or in the near UV. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between the global maps of 0.5- 1.0- micron optical depth derived from the AVHRR and SAGE satellite measurements.

Russell, Philip B.↗

NASA NTRS1983

An analytical evaluation of the icing properties of several low and medium speed airfoils

A droplet trajectory computer code is used to analyze the water droplet impingement characteristics of several low and medium speed airfoils. The results are reported in terms of the airfoil maximum impingement efficiency, and limits of impingement. The airfoils are analyzed for angles of attack of -4 to 12 degrees, and modified inertia parameters of 0.01 to 1.0. These data are then used to evaluate the effect of airfoil geometry: leading edge radius, maximum thickness, maximum camber, and angle of attack on airfoil icing characteristics. Airfoils from the NACA four digit series to the new NASA NLF airfoils are considered. The strongest correlations were seen between airfoil leading edge radius and maximum local collection efficiency, thickness and total collection efficiency, and camber and the limits of impingement.

Bragg, M. B.↗

NASA NTRS2004

Backscattering enhancement with a finite beam width for millimeter - wavelength weather radars

Backscattering enhancement from random hydrometeors should increase as wavelengths of radars reach millimeter regions. For 95 GHz radars, the reflectivity of backscattering is expected to increase by 2 dB, due to multiple scattering including backscattering enhancement, for water droplets of diameter of 1 mm with a density of 5×103 m−3. Previous theoretical studies of backscattering enhancement considered infinitely extending plane waves. In this paper, we expand the theory to spherical waves with a Gaussian antenna pattern, including depolarizing effects. While the differences from the plane wave results are not great when the optical thickness is small, as the latter increases the differences become significant, and essentially depend on the ratio of radar footprint radius to the mean free path of hydrometeors. In this regime, for a radar footprint that is smaller than the mean free path, the backscattering-enhancement reflectivity corresponding to spherical waves is significantly less pronounced than in the case of the plane wave theory. Hence this reduction factor must be taken into account when analyzing radar reflectivity factors for use in remote sensing applications.

Im, Eastwood↗

NASA NTRSDate not supplied

Surface Radiation Budget (SRB) Release 2.5 Longwave cloud props 3 hourly Data in Native Format (SRB_REL2.5_LW_CLDPROPS_3HRLY)

The Surface Radiation Budget (SRB) data set contain 47 surface and cloud properties including cloud fraction, cloud visible optical depth, cloud top temperature, cloud top pressure, cloud base pressure, and cloud water and ice content for each of 5 cloud types and in total. In addition to these cloud properties, the parameters also include surface pressure, total column ozone, day/night flag, and precipitable water. All perameters were measured at three hourly intervals for each day for the entire globe between 07/01/1983 and 06/30/2005 and serve as inputs to the Global Energy and Water-Cycle Experiment (GEWEX) longwave algorithm (GLW). [Location=GLOBAL] [Temporal_Coverage: Start_Date=1998-01-01; Stop_Date=2005-06-30] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1 degree; Longitude_Resolution=Ranges from 1 degree (tropics and subtropics) to 120 degrees (the poles).; Temporal_Resolution=3 hourly; Temporal_Resolution_Range=3 hourly].

CLOUD TOP TEMPERATURE↗

NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

NASA NTRS1998

Undercooling, Liquid Separation and Solidification of Cu-Co Alloys

Large undercooling can induce not only various solidification pathways, but also a precursor reaction, or liquid separation. This paper deals with the latter effect of undercooling using examples of the Cu-Co system which has a flattened liquidus. Bulk Cu-Co alloys (about 7mm diameter) at compositions ranging from 10 to 90 wt pct Co were highly undercooled using a fluxing technique. Except for Cu-90 wt pct Co, liquid separation was directly observed as undercooling exceeded a critical value depending on the composition. It was also confirmed by a microstructural transition from dendrites to droplets above the critical undercooling. Finally, theoretical calculations regarding the metastable miscibility boundary and maximum droplet radius were made to analyze the experimental results.

Robinson, M. B.↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

Modeling and Simulation of Tank Pressure Control using Zero-Boiloff 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↗

NASA NTRS2013

A Simple Model of Global Aerosol Indirect Effects

Most estimates of the global mean indirect effect of anthropogenic aerosol on the Earth's energy balance are from simulations by global models of the aerosol lifecycle coupled with global models of clouds and the hydrologic cycle. Extremely simple models have been developed for integrated assessment models, but lack the flexibility to distinguish between primary and secondary sources of aerosol. Here a simple but more physically based model expresses the aerosol indirect effect (AIE) using analytic representations of cloud and aerosol distributions and processes. Although the simple model is able to produce estimates of AIEs that are comparable to those from some global aerosol models using the same global mean aerosol properties, the estimates by the simple model are sensitive to preindustrial cloud condensation nuclei concentration, preindustrial accumulation mode radius, width of the accumulation mode, size of primary particles, cloud thickness, primary and secondary anthropogenic emissions, the fraction of the secondary anthropogenic emissions that accumulates on the coarse mode, the fraction of the secondary mass that forms new particles, and the sensitivity of liquid water path to droplet number concentration. Estimates of present-day AIEs as low as 5 W/sq m and as high as 0.3 W/sq m are obtained for plausible sets of parameter values. Estimates are surprisingly linear in emissions. The estimates depend on parameter values in ways that are consistent with results from detailed global aerosol-climate simulation models, which adds to understanding of the dependence on AIE uncertainty on uncertainty in parameter values.

Clouds↗

NASA NTRS1981

The formation and growth of sulfate aerosols in the stratosphere

The critical compositions, critical radii, and nucleation rates of sulfuric acid droplets formed through homogeneous heteromolecular nucleation process in the stratosphere under a variety of ambient conditions are calculated. Results are plotted as a function of H2SO4 concentration and H2O concentration so that characteristics of sulfate aerosols can be quickly estimated whenever environmental conditions are known. The effect of the uncertainty in the saturated vapor pressure of H2SO4 on the estimated properties of sulfate aerosols is elucidated. Plottings are presented that illustrate the relative importance of new particle formation by homogeneous nucleation vs growth by heteromolecular condensation on pre-existing nuclei in the stratosphere. The effects of ambient conditions, such as the concentrations of H2SO4 and H2O, the average stratospheric particle radius, and the number concentration of pre-existing nuclei on the relative importance of these two competing processes are discussed.

Yue, G. K.↗

NASA NTRSDate not supplied

Look-up Table for Lidar and Polarimeter Liquid Water Cloud Studies

- Combined lidar and polarimeter retrievals of aerosol, cloud and ocean microphysical properties involve single scattering cloud calculations that are time consuming. - We create a look-up table to speed up calculations for water droplets in the atmosphere. The gained increase in computational efficiency is up to 10^(4). - In our new Lorenz-Mie look-up table we tabulate the light scattering by an ensemble of homogeneous isotropic spheres at wavelengths starting from 0.35 µm. - The look-up table covers liquid water cloud particles with radii in the range of 0.001—500 µm. - The covered complex refractive indices range from 1.25 to 1.36 for the real part and from 0 to 0.001 for the imaginary part. - Using look-up table, we can precisely compute inherent optical properties for the particle size distributions ranging up to 100 µm for the effective radius and up to 0.6 for the effective variance. - We test wavelengths from 0.35 to 2.3 µm and find that the P12 element of the normalized scattering matrix is precise to within 1% for 96.7% of cases; the absorption and backscatter coefficients, and lidar ratio are precise to within 1% for 99% of cases; the P11 , P33 and P34 elements of the normalized scattering matrix, the asymmetry parameter, the extinction and scattering coefficients are precise to within 1% for all cases. - We provide an example of using the look-up table with in situ measurements to determine agreement with remote sensing. - The look-up table along with examples of its use soon will be freely available at https://science.larc.nasa.gov/polarimetry

Eduard Chemyakin↗

NASA NTRS1983

Retrieval of stratospheric aerosol size distribution from atmospheric extinction of solar radiation at two wavelengths

A parametric study is presented of the effects of the size distribution and composition of stratospheric aerosols on the value R, the ratio of the extinction of solar radiation at wavelength 0.45 micron to the extinction at wavelength 1.0 micron. The aerosol size distributions under study comprise nine analytical expressions, including most of the stratospheric aerosol models used by investigators. The aerosol compositions considered are the supercooled sulfuric-acid droplets, with different weight percentages of H2SO4 in the aerosol. It is found that whereas R is not very sensitive to either the composition or the radii limits of the stratospheric aerosols under consideration, it is quite sensitive to the value of the variable parameter governing the mode radius of the size distribution of aerosol particles. The results of the parametric study underlie the method proposed here of retrieving the size distribution of aerosol particles in the stratosphere from the experimental results of the extinction of solar radiation at two wavelengths.

Yue, G. K.↗

NASA NTRS1953

Maximum Evaporation Rates of Water Droplets Approaching Obstacles in the Atmosphere Under Icing Conditions

When a closed body or a duct envelope moves through the atmosphere, air pressure and temperature rises occur ahead of the body or, under ram conditions, within the duct. If cloud water droplets are encountered, droplet evaporation will result because of the air-temperature rise and the relative velocity between the droplet and stagnating air. It is shown that the solution of the steady-state psychrometric equation provides evaporation rates which are the maximum possible when droplets are entrained in air moving along stagnation lines under such conditions. Calculations are made for a wide variety of water droplet diameters, ambient conditions, and flight Mach numbers. Droplet diameter, body size, and Mach number effects are found to predominate, whereas wide variation in ambient conditions are of relatively small significance in the determination of evaporation rates. The results are essentially exact for the case of movement of droplets having diameters smaller than about 30 microns along relatively long ducts (length at least several feet) or toward large obstacles (wings), since disequilibrium effects are then of little significance. Mass losses in the case of movement within ducts will often be significant fractions (one-fifth to one-half) of original droplet masses, while very small droplets within ducts will often disappear even though the entraining air is not fully stagnated. Wing-approach evaporation losses will usually be of the order of several percent of original droplet masses. Two numerical examples are given of the determination of local evaporation rates and total mass losses in cases involving cloud droplets approaching circular cylinders along stagnation lines. The cylinders chosen were of 3.95-inch (10.0+ cm) diameter and 39.5-inch 100+ cm) diameter. The smaller is representative of icing-rate measurement cylinders, while with the larger will be associated an air-flow field similar to that ahead of an airfoil having a leading-edge radius comparable with that of the cylinder. It is found that the losses are less than 5 percent. It is concluded that such losses are, in general, very small (less than 1 percent) in the case of smaller obstacles (of icing-rate measurement- cylinder size); the motional dynamics are such, however, that exceptions will occur by reason of failure of very small droplets (moving along stagnation lines) to impinge upon obstacle surfaces. In such cases, the droplets will evaporate completely.

Lowell, H. H.↗

NASA NTRS2004

Parameterization of the Vertical Variability of Tropical Cirrus Cloud Microphysical and Optical Properties

Cloud Integrating Nephelometers (CIN) were flown on the U. North Dakota Citation aircraft and the NASA WB-57 aircraft for the purpose of measuring in-situ the optical extinction coefficient and the asymmetry parameter (g) at a wavelength of 635 nm of primarily ice particles encountered during the NASA CRYSTAL-FACE study of large cumulus clouds (Cu) and their anvils found in the southern Florida region. The probes performance was largely successful and produced archived data for vertical profiles of extinction, asymmetry parameter, and effective radius (Re), the latter being obtained by combining CIN and CVI (total water; Oregon State U.) measurements. Composites of the CIN and CVI data describing the average microphysical and optical behavior of the Cu and their anvils showed the following: The extinction increases with height as a result of the size of the particles also decreasing with height as shown by the Re measurements; near the top of anvils the size of the primary ice particles is about 10-um radius; and the value of g does not vary significantly with height and has a mean value of about 0.73 consistent with the idea that ambient ice crystals are primarily of complex shape and reflect solar radiation more efficiently than particles of pristine crystal shape. Other observations include: The g measurements were found to be an indicator of the phase of the cloud permitting identification of the clouds with water droplets, rain, and ice; visual ranges as small as several tens of meters were occasionally found in "extinction cores" that coincided with strong updraft cores; and comparison of the cloud probes on the Citation showed significant disagreement.

Gerber, Hermann E.↗