Search NASA⌕ Search

SEARCH · Search NASA

Results for “freeze”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A model of the formation of high-salinity shelf water on polar continental shelves

This study presents a model of the flow and salinity fields forced by sea-surface salinity flux and wind stress curl and examines the processes that create High-Salinity Shelf Water (HSSW). To investigate the relative effects of wind stress and buoyancy forcing in HSSW production, the polynya freezing rate in the model is varied from 0.0 to 0.30 m/d, and the Ekman pumping derived from the wind stress curl is varied independently from 0.0 to 1.8 x 10 exp 6 m/s. The Ekman pumping was found to control the magnitude of the circulation, while the polynya freezing rate controlled the extent of salinization in the shelf water. The flux of HSSW increases linearly with increasing Ekman pumping above 0.3 x 10 exp -6 m/s. The modeled flux of HSSW and the flux of derived Botttom Water for the present estimates of the forcings (a peak freezing rate of 0.10 m/d and Ekman pumping of 0.2 x 10 exp -6 m/s) agree with the fluxes inferred from physical and chemical observations in the deep Weddel Sea by oceanographic field programs.

Grumbine, Robert W.↗
NASA NTRS1993

Sea ice radar signatures from ERS-1 SAR during late Summer and Fall in the Beaufort and Chukchi Seas

A study which examines ERS-1 C band SAR (Synthetic Aperture Radar) imagery of sea ice obtained in the Beaufort and Chukchi Seas from mid Summer through Fall freeze up and early Winter in 1991 is presented. Radar backscatter statistics of sea ice were obtained from the imagery, using common floes tracked through consecutive repeat images whenever possible. During the Summer months, strong fluctuations in ice signatures of several dB are observed over 2 to 3 day periods, which are found to be closely related to air temperature excursions above and below freezing that alters the phase of the ice surface. As air temperatures drop steadily below freezing in the Fall, the signatures of the pack ice increase in brightness and become more stable with time. Multiyear ice is distinguished from rough and smooth first year ice. There are also variations in the multiyear signatures with latitude. Large variations are seen in new ice and open water contained within leads which results in ambiguous classification.

Holt, Benjamin↗
NASA NTRS1993

Thermal and hydraulic considerations regarding the fate of water discharged by the outflow channels to the Martian northern plains

The identification of possible shorelines in the Martian northern plains suggests that the water discharged by the circum-Chryse outflow channels may have led to the formation of transient seas, or possibly even an ocean, covering as much as one-third of the planet. Speculations regarding the possible fate of this water have included local ponding and reinfiltration into the crust; freezing, sublimation, and eventual cold-trapping at higher latitudes; or the in situ survival of this now frozen water to the present day -- perhaps aided by burial beneath a protective cover of eolian sediment or lavas. Although neither cold-trapping at higher latitudes nor the subsequent freezing and burial of flood waters can be ruled out, thermal and hydraulic considerations effectively eliminate the possibility that any significant reassimilation of this water by local infiltration has occurred given climatic conditions resembling those of today. The arguments against the local infiltration of flood water into the northern plains are two-fold. First, given the climatic and geothermal conditions that are thought to have prevailed on Mars during the Late Hesperian (the period of peak outflow channel activity in the northern plains), the thickness of the cryosphere in Chryse Planitia is likely to have exceeded 1 km. A necessary precondition for the widespread occurrence of groundwater is that the thermodynamic sink represented by the cryosphere must already be saturated with ice. For this reason, the ice-saturated cryosphere acts as an impermeable barrier that effectively precludes the local resupply of subpermafrost groundwater by the infiltration of water discharged to the surface by catastraphic floods. Note that the problem of local infiltration is not significantly improved even if the cryosphere were initially dry, for as water attempts to infiltrate the cold, dry crust, it will quickly freeze, creating a seal that prevents any further infiltration from the ponded water above. The second argument against the local infiltration of flood water in the northern plains is based on hydraulic considerations. Repeated impacts have likely brecciated the Martian crust down to a depth of roughly 10 km. Given a value of permeability no greater than that inferred for the top 10 km of the Earth's crust (approximately 10(exp -2) darcies), a timescale as much as a billion years or more for the Martian groundwater system to achieve hydrostatic equilibrium, and the approximately 2-4 km elevation difference between the outflow channel source regions and the northern plains, the water confined beneath the frozen crust of the northern plains should have been under a significant hydraulic head. Thus, the existence of a hydraulic pathway between the ponded flood waters above the northern plains and the confined aquifer lying beneath it would not have led to the infiltration of flood water back into the crust, but rather the additional expulsion of groundwater onto the surface.

Clifford, S. M.↗
NASA NTRS1994

Influence of convection on microstructure

The primary motivation for this research was to determine the cause for space processing altering the microstructure of some eutectics, especially the MnBi-Bi eutectic. Four primary hypotheses were to be tested under this current grant: (1) A fibrous microstructure is much more sensitive to convection than a lamellar microstructure, which was assumed in our prior theoretical treatment. (2) An interface with one phase projecting out into the melt is much more sensitive to convection than a planar interface, which was assumed in our prior theoretical treatment. (3) The Soret effect is much more important in the absence of convection and has a sufficiently large influence on microstructure that its action can explain the flight results. (4) The microstructure is much more sensitive to convection when the composition of the bulk melt is off eutectic. These hypotheses were tested. It was concluded that none of these can explain the Grumman flight results. Experiments also were performed on the influence of current pulses on MnBi-Bi microstructure. A thorough review was made of all experimental results on the influence of convection on the fiber spacing in rod eutectics, including results from solidification in space or at high gravity, and use of mechanical stirring or a magnetic field. Contradictory results were noted. The predictions of models for convective influences were compared with the experimental results. Vigorous mechanical stirring appears to coarsen the microstructure by altering the concentration field in front of the freezing interface. Gentle convection is believed to alter the microstructure of a fibrous eutectic only when it causes a fluctuating freezing rate with a system for which the kinetics of fiber branching differs from that for fiber termination. These fluctuations may cause the microstructure to coarsen or to become finer, depending on the relative kinetics of these processes. The microstructure of lamellar eutectics is less sensitive to freezing rate fluctuations and to gentle convection.

Wilcox, William R.↗
NASA NTRS1994

The tropical water and energy cycles in a cumulus ensemble model. Part 1: Equilibrium climate

A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Sui, C. H.↗
NASA NTRS1994

Analysis of ERS 1 synthetic aperture radar data of frozen lakes in northern Montana and implications for climate studies

Lakes that freeze each winter are good indicators of regional climate change if key parameters, such as freeze-up and breakup date and maximum ice thickness, are measured over a decade-scale time frame. Synthetic aperture radar (SAR) satellite data have proven to be especially useful for measurement of climatologically significant parameters characteristic of frozen lakes. In this paper, five lakes in Glacier National Park, Montana, have been studied both in the field and using Earth Remote-Sensing Satellite (ERS) 1 SAR data during the 1992-1993 winter. The lakes are characterized by clear ice, sometimes with tubular or rounded bubbles, and often with a layer of snow ice on top of the clear ice. They are also often snow covered. Freeze-up is detected quite easily using ERS 1 SAR data as soon as a thin layer of ice forms. The effect of snow ice on the backscatter is thought to be significant but is, as yet, undetermined. On the five lakes studied, relative backscatter was found to increase with ice thickness until a maximum was reached in February. Breakup, an often ill-defined occurrence, is difficult to detect because surface water causes the SAR signal to be absorbed, thus masking the ice below. Comparison of the bubble structure of thaw lakes in northern Alaska with lakes in northern Montana has shown that the ice structure is quite different, and this difference may contribute to differential SAR signature evolution in the lakes of the two areas.

Hall, Dorothy K.↗
NASA NTRS1995

Getting granite dikes out of the source region

Whether a dike can propagate far from a magma reservoir depends upon the competition between the rate at which propagation widens the dike and the rate at which freezing constricts the aperture available for magma flow. Various formulations are developed for a viscous fluid at temperature T(sub m) intruding a growing crack in an elastic solid. The initial solid temperature equals T(sub m) at the source and decreases linearly with distance from the source. If T(sub m) is the unique freezing temperature of the fluid, dike growth is initially self-similar and an essentially exact solution is obtained; if T(sub m) is above the solidus temperature, the solution is approximate but is designed to overestimate the distance the dike may propagate. The ability of a dike to survive thermally depends primarily upon a single parameter that is a measure of the ratio of the dike frozen margin thickness to elastic thickness. Perhaps more intuitively, one may define a minimum distance from the essentially solid reservoir wall to the point at which the host rock temperature drops below the solidus, necessary for dikes to propagate far into subsolidus rock. It is concluded that for reasonable material properties and source conditions, most basalt dikes will have little difficulty leaving the source region, but most rhyolite dikes will be halted by freezing soon after the magma encounters rock at temperatures below the magma solidus. While these results can explain why granitic dikes are common near granitic plutons but rare elsewhere, the potentially large variation in magmatic systems makes it premature to rule out the possibility that most granites are transported through the crust in dikes. Nonetheless, these results highlight difficulties with such proposals and suggest that it may also be premature to rule out the possibility that most granite plutons ascend as more equidimensional bodies.

Rubin, Allan M.↗
NASA NTRS1996

Ice Detector and Deicing Fluid Effectiveness Monitoring System

An ice detector and deicing fluid effectiveness monitoring system for an aircraft is disclosed. The ice detection portion is particularly suited for use in flight to notify the flight crew of an accumulation of ice on an aircraft lifting and control surfaces, or helicopter rotors, whereas the deicing fluid effectiveness monitoring portion is particularly suited for use on the ground to notify the flight crew of the possible loss of the effectiveness of the deicing fluid. The ice detection portion comprises a temperature sensor and a parallel arrangement of electrodes whose coefficient of coupling is indicative of the formation of the ice, as well as the thickness of the formed ice. The fluid effectiveness monitoring portion comprises a temperature sensor and an ionic-conduction cell array that measures the conductivity of the deicing fluid which is indicative of its concentration and, thus, its freezing point. By measuring the temperature and having knowledge of the freezing point of the deicing fluid, the fluid effectiveness monitoring portion predicts when the deicing fluid may lose its effectiveness because its freezing point may correspond to the temperature of the ambient.

Seegmiller, H. Lee B.↗
NASA NTRS1996

Dynamics of the Molten Contact Line

In contrast to the ordinary contact line problem, virtually no information is available on the similar problem associated with a molten material spreading on a solid which is below the melt's fusion point. The latter is a more complex problem which heat transfer and solidification take place simultaneously with spreading, and requires answers not only for the hot melt's advance speed over the cold solid as a function of contact angle, but also for how one is to predict the point of the molten contact line's arrest by freezing. This issues are of importance in evolving methods of materials processing. The purpose of our work is to develop, based on both experiments and theory, an understanding of the dynamic processes that occur when a molten droplet touches a subcooled solid, spreads partly over it by capillary action, and freezes. We seek answers to the following basic questions. First, what is the relationship between the melt's contact line speed and the apparent (dynamic) contact angle? Secondly, at what point will the contact line modon be arrested by freezing? The talk will describe three components of our work: (1) deposition experiments with small molten droplets; (2) investigation of the dynamics of the molten contact line by means of a novel forced spreading method; and (3) an attempt to provide a theoretical framework for answering the basic questions posed above.

Sonin, Ain A.↗
NASA NTRS1995

Thermal Characterization of the Universal Multizone Crystallizator

The Universal Multizone Crystallizator (UMC) is a special apparatus for crystal growth under terrestrial and microgravity conditions. The use of twenty-five zones allows the, UMC to be used for several normal freezing growth techniques. The thermal profile is electronically translated along the stationary sample by systematically reducing the power to the control zones. Elimination of mechanical translation devices increases the systems reliability while simultaneously reducing the size and weight. This paper addresses the UMC furnace design, sample cartridge, typical thermal profiles and corresponding power requirements necessary for two normal freezing techniques: dynamic gradient freeze and zone melting crystal growth.

Watring, D. A.↗
NASA NTRS1997

DRA/NASA/ONERA Collaboration on Icing Research: Prediction of Airfoil Ice Accretion - Part 2

This report presents results from a joint study by DRA, NASA, and ONERA for the purpose of comparing, improving, and validating the aircraft icing computer codes developed by each agency. These codes are of three kinds: (1) water droplet trajectory prediction, (2) ice accretion modeling, and (3) transient electrothermal deicer analysis. In this joint study, the agencies compared their code predictions with each other and with experimental results. These comparison exercises were published in three technical reports, each with joint authorship. DRA published and had first authorship of Part 1 - Droplet Trajectory Calculations, NASA of Part 2 - Ice Accretion Prediction, and ONERA of Part 3 - Electrothermal Deicer Analysis. The results cover work done during the period from August 1986 to late 1991. As a result, all of the information in this report is dated. Where necessary, current information is provided to show the direction of current research. In this present report on ice accretion, each agency predicted ice shapes on two dimensional airfoils under icing conditions for which experimental ice shapes were available. In general, all three codes did a reasonable job of predicting the measured ice shapes. For any given experimental condition, one of the three codes predicted the general ice features (i.e., shape, impingement limits, mass of ice) somewhat better than did the other two. However, no single code consistently did better than the other two over the full range of conditions examined, which included rime, mixed, and glaze ice conditions. In several of the cases, DRA showed that the user's knowledge of icing can significantly improve the accuracy of the code prediction. Rime ice predictions were reasonably accurate and consistent among the codes, because droplets freeze on impact and the freezing model is simple. Glaze ice predictions were less accurate and less consistent among the codes, because the freezing model is more complex and is critically dependent upon unsubstantiated heat transfer and surface roughness models. Thus, heat transfer prediction methods used in the codes became the subject for a separate study in this report to compare predicted heat transfer coefficients with a limited experimental database of heat transfer coefficients for cylinders with simulated glaze and rime ice shapes. The codes did a good job of predicting heat transfer coefficients near the stagnation region of the ice shapes. But in the region of the ice horns, all three codes predicted heat transfer coefficients considerably higher than the measured values. An important conclusion of this study is that further research is needed to understand the finer detail of of the glaze ice accretion process and to develop improved glaze ice accretion models.

Wright, William B.↗
NASA NTRS1997

Determination of Stress Intensity Factor Distributions for "Interface" Cracks in Incompressible, Dissimilar Materials

The present study was undertaken in order to develop test methods and procedures for measuring the variation of the stress intensity factor through the thickness in bimaterial specimens containing cracks within and parallel to the bond line using the frozen stress photoelastic method. Since stress freezing materials are incompressible above critical temperature, and since thick plates are to be employed which tend to produce a state of plane strain near the crack tip, the interface near tip fracture equations reduce to the classic form for homogeneous materials. Moreover, zero thickness interfaces do not exist when materials are bonded together. It was decided early on that it would be important to insure a uniform straight and accurate crack tip region through the thickness of the body to reduce scatter in the SIF distribution through the thickness. It was also observed that rubberlike materials which were desired to be modeled exhibited significant tip blunting prior to crack extension and that some blunting of the tip would provide a more realistic model. It should be noted that, in normal stress freezing photoelastic work, it is considered good practice to avoid utilizing data near bond lines in photoelastic models due to the bond line stresses which inevitably develop when two parts are bonded together. Thus, the present study involves certain exploratory aspects in deviating from standard practice in stress freezing work. With the above ideas in mind, several different test methods were investigated and are described in the following sections and appendices. The geometry selected for the program was a thick, edge cracked specimen containing a bond line.

Smith, C. W.↗
NASA NTRS2000

Experimental Results from the Thermal Energy Storage-2 (TES-2) Flight Experiment

Thermal Energy Storage-2 (TES-2) is a flight experiment that flew on the Space Shuttle Endeavour (STS-72), in January 1996. TES-2 originally flew with TES-1 as part of the OAST-2 Hitchhiker payload on the Space Shuttle Columbia (STS-62) in early 1994. The two experiments, TES-1 and TES-2 were identical except for the fluoride salts to be characterized. TES-1 provided data on lithium fluoride (LiF), TES-2 provided data on a fluoride eutectic (LiF/CaF2). Each experiment was a complex autonomous payload in a Get-Away-Special payload canister. TES-1 operated flawlessly for 22 hr. Results were reported in a paper entitled, Effect of Microgravity on Materials Undergoing Melting and Freezing-The TES Experiment, by David Namkoong et al. A software failure in TES-2 caused its shutdown after 4 sec of operation. TES-1 and 2 were the first experiments in a four experiment suite designed to provide data for understanding the long duration microgravity behavior of thermal energy storage salts that undergo repeated melting and freezing. Such data have never been obtained before and have direct application for the development of space-based solar dynamic (SD) power systems. These power systems will store energy in a thermal energy salt such as lithium fluoride or a eutectic of lithium fluoride/calcium difluoride. The stored energy is extracted during the shade portion of the orbit. This enables the solar dynamic power system to provide constant electrical power over the entire orbit. Analytical computer codes were developed for predicting performance of a space-based solar dynamic power system. Experimental verification of the analytical predictions were needed prior to using the analytical results for future space power design applications. The four TES flight experiments were to be used to obtain the needed experimental data. This paper will address the flight results from the first and second experiments, TES-1 and 2, in comparison to the predicted results from the Thermal Energy Storage Simulation (TESSIM) analytical computer code. An analysis of the TES-2 data was conducted by Cleveland State University Professor, Mounir Ibrahim. TESSIM validation was based on two types of results; temperature history of various points on the containment vessel and TES material distribution within the vessel upon return from flight. The TESSIM prediction showed close comparison with the flight data. Distribution of the TES material within the vessel was obtained by a tomography imaging process. The frozen TES material was concentrated toward the colder end of the canister. The TESSIM prediction indicated a similar pattern. With agreement between TESSIM and the flight data, a computerized representation was produced to show the movement and behavior of the void during the entire melting and freezing cycles.

Tolbert, Carol↗
NASA NTRS2000

The Cirrus Parcel Model Comparison Project

The cirrus Parcel Model Comparison Project involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 centimeters per second, respectively. Five models are participating in the project. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found amongst the models for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Systematic bias exists between results of a model based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies. This information is necessary, but not sufficient, to construct consistent formulae for the two approaches. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period.

Lin, Ruei-Fong↗
NASA NTRS2000

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗
NASA NTRS2001

Ice Nucleation in Deep Convection

The processes controlling production of ice crystals in deep, rapidly ascending convective columns are poorly understood due to the difficulties involved with either modeling or in situ sampling of these violent clouds. A large number of ice crystals are no doubt generated when droplets freeze at about -40 C. However, at higher levels, these crystals are likely depleted due to precipitation and detrainment. As the ice surface area decreases, the relative humidity can increase well above ice saturation, resulting in bursts of ice nucleation. We will present simulations of these processes using a large-eddy simulation model with detailed microphysics. Size bins are included for aerosols, liquid droplets, ice crystals, and mixed-phase (ice/liquid) hydrometers. Microphysical processes simulated include droplet activation, freezing, melting, homogeneous freezing of sulfate aerosols, and heterogeneous ice nucleation. We are focusing on the importance of ice nucleation events in the upper part of the cloud at temperatures below -40 C. We will show that the ultimate evolution of the cloud in this region (and the anvil produced by the convection) is sensitive to these ice nucleation events, and hence to the composition of upper tropospheric aerosols that get entrained into the convective column.

Jensen, Eric↗
NASA NTRS2001

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗
NASA NTRS2002

Simulation of TRMM Microwave Imager Brightness Temperature using Precipitation Radar Reflectivity for Convective and Stratiform Rain Areas over Land

Rain is highly variable in space and time. In order to measure rainfall over global land with satellites, we need observations with very high spatial resolution and frequency in time. On board the Tropical Rainfall Measuring Mission (TRMM) satellite, the Precipitation Radar (PR) and Microwave Imager (TMI) are flown together for the purpose of estimating rain rate. The basic method to estimate rain from PR has been developed over the past several decades. On the other hand, the TMI method of rain estimation is still in the state development, particularly over land. The objective of this technical memorandum is to develop a theoretical framework that helps relate the observations made by these two instruments. The principle result of this study is that in order to match the PR observations with the TMI observations in convective rain areas, a mixed layer of graupel and supercooled water drops above the freezing level is needed. On the other hand, to match these observations in the stratiform region, a layer of snowflakes with appropriate densities above the freezing level, and a melting layer below the freezing level, are needed. This understanding can lead to a robust rainfall estimation technique from the microwave radiometer observations.

Prabhakara, C.↗