Search NASA⌕ Search

SEARCH · Search NASA

Results for “Convective mixing layer”

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.

352 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗
NASA NTRS1994

Direct numerical simulations and modeling of a spatially-evolving turbulent wake

Understanding of turbulent free shear flows (wakes, jets, and mixing layers) is important, not only for scientific interest, but also because of their appearance in numerous practical applications. Turbulent wakes, in particular, have recently received increased attention by researchers at NASA Langley. The turbulent wake generated by a two-dimensional airfoil has been selected as the test-case for detailed high-resolution particle image velocimetry (PIV) experiments. This same wake has also been chosen to enhance NASA's turbulence modeling efforts. Over the past year, the author has completed several wake computations, while visiting NASA through the 1993 and 1994 ASEE summer programs, and also while on sabbatical leave during the 1993-94 academic year. These calculations have included two-equation (K-omega and K-epsilon) models, algebraic stress models (ASM), full Reynolds stress closure models, and direct numerical simulations (DNS). Recently, there has been mutually beneficial collaboration of the experimental and computational efforts. In fact, these projects have been chosen for joint presentation at the NASA Turbulence Peer Review, scheduled for September 1994. DNS calculations are presently underway for a turbulent wake at Re(sub theta) = 1000 and at a Mach number of 0.20. (Theta is the momentum thickness, which remains constant in the wake of a two dimensional body.) These calculations utilize a compressible DNS code written by M. M. Rai of NASA Ames, and modified for the wake by J. Cimbala. The code employs fifth-order accurate upwind-biased finite differencing for the convective terms, fourth-order accurate central differencing for the viscous terms, and an iterative-implicit time-integration scheme. The computational domain for these calculations starts at x/theta = 10, and extends to x/theta = 610. Fully developed turbulent wake profiles, obtained from experimental data from several wake generators, are supplied at the computational inlet, along with appropriate noise. After some adjustment period, the flow downstream of the inlet develops into a fully three-dimensional turbulent wake. Of particular interest in the present study is the far wake spreading rate and the self-similar mean and turbulence profiles. At the time of this writing, grid resolution studies are underway, and a code is being written to calculate turbulence statistics from these wake calculations; the statistics will be compared to those from the ongoing PIV wake measurements, those of previous experiments, and those predicted by the various turbulence models. These calculations will lead to significant long-term benefits for the turbulence modeling effort. In particular, quantities such as the pressure-strain correlation and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. Present turbulence models do a very good job at prediction of the shape of the mean velocity and Reynolds stress profiles in a turbulent wake, but significantly underpredict the magnitude of the stresses and the spreading rate of the wake. Thus, the turbulent wake is an ideal flow for turbulence modeling research. By careful comparison and analysis of each term in the modeled Reynolds stress equations, the DNS data can show where deficiencies in the models exist; improvements to the models can then be attempted.

Cimbala, John M.↗
NASA NTRS1999

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗
NASA NTRS2003

Digital Particle Image Velocimetry (DPIV) Used for Space-Time Correlations in Nozzle Flow

An optical measurement technique known as Digital Particle Image Velocimetry (DPIV) was used previously to characterize the first- and second-order statistical properties of both cold and hot jet flows from externally mixed nozzles in NASA Glenn Research Center's Nozzle Acoustic Test Rig. In this technique, an electronic camera records particles entrained in a flow as a laser light sheet is pulsed at two instances in time. Correlation processing of the recorded particle image pairs yields the two-component velocity field across the imaged plane of the flow. The information acquired using DPIV is being used to improve our understanding of the decay of turbulence in jet flows-a critical element for understanding the acoustic properties of the flow. Recently, two independent DPIV systems were installed in Glenn's Small Hot Jet Acoustic Rig, enabling multiplane correlations in time and space. The data were collected over a range of different Mach numbers and temperature ratios. DPIV system 1 was fixed to a large traverse rig, and DPIV system 2 was mounted on a small traverse system mounted on the large traverse frame. The light sheets from the two DPIV systems were aligned to lie in the same axial plane, with DPIV system 2 being independently traversed downstream along the flow direction. For each measurement condition, the DPIV systems were started at a fully overlapping orientation. A polarization separation technique was used to avoid cross-talk between the two systems. Then, the DPIV systems fields were shifted axially apart, in successively increasing steps. The downstream DPIV system 2 was triggered at a short time delay after the upstream DPIV system 1, where the time delay was proportional to the convective flow velocity in the shear layer of the jet flow and the axial separation of the two DPIV systems. The acquired data were processed to obtain the instantaneous velocity vector maps over a range of time delays and spatial separations. The velocity fields from the different DPIV systems were then cross-correlated to determine the degree of correlation remaining in the flow as the downstream convection distance was increased. The new data provide Lagrangian measurements of the convective turbulent structures in the shear layer of an exhaust nozzle. These measurements, obtained in both cold and hot flows, will be used to validate and correct models for space-time velocity correlations-long a missing key to predicting jet noise.

Wernet, Mark P.↗
NASA NTRS2008

Using Wind and Temperature Fields to Study Dehydration Mechanisms in the Tropical Tropopause Layer

The tropics are the main region for troposphere-to-stratosphere transport (TST) of air. One of the dominant mechanisms that control tropical TST of water vapor is freeze-drying by the cold tropical tropopause. This mechanism is supported by evidence from satellite observations of the "tape recorder", where seasonal changes in stratospheric water vapor are in phase with seasonal changes in tropopause temperatures in the tropics. Over the last few years, however, the concept of the tropical tropopause has evolved from a single material surface to a layer called the Tropical Tropopause Layer (TTL). A recent hypothesis on dehydration mechanisms suggests that dehydration and entry point into the stratosphere are not always co-located (Holton and Gettelman, 2001). Instead, dehydration can occur during horizontal advection through Lagrangian 'cold pools', or coldest regions along a parcel's trajectory, as air ascends within the TTL while the entry point into the stratosphere occurs at a different geographical location. In this study, we investigate the impact that these Lagrangian cold pools have on TTL moisture. For this purpose, we use in situ measurements of TTL water vapor obtained aboard NASA's WB-57 aircraft over the Eastern Tropical Pacific, and we compare these measurements to minimum saturation water vapor mixing ratios obtained from three-dimensional backward trajectory calculations. Aircraft measurements show frequent unsaturated conditions, which suggest that the entry value of stratospheric water vapor in this region was not set by local saturation conditions. Trajectory calculations, driven by both ECMWF operational analysis and reanalysis winds and temperature fields, are used to explore the impact (e.g., geographical location, timing, dehydration magnitude) of the Lagrangian cold pools intercepted by the parcels sampled by the aircraft. We find noteworthy differences in the location of the Lagrangian cold pools using the two ECMWF data sets, namely influence of the Western Tropical Pacific region when using operational analysis fields versus influence of the Eastern Tropical Pacific and South America regions when using reanalysis fields. These results have a significant impact on our scientific conclusions on dehydration mechanisms affecting the air sampled by the aircraft, given that these regions have different thermodynamic and convective properties.

Pittman, Jasna↗
NASA NTRS2000

Stereo Imaging Velocimetry of Mixing Driven by Buoyancy Induced Flow Fields

Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.

Duval, W. M. B.↗
NASA NTRS2000

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗
NASA NTRS2004

Florida Thunderstorms: A Faucet of Reactive Nitrogen to the Upper Troposphere

During the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) enhanced mixing ratios of nitric oxide were measured in the anvils of thunderstorms and in clear air downwind of storm systems on flights of a Wl3-57F high-altitude aircraft. Mixing ratios greater than l0 - 20 times background were readily observed over distances of 25-120 km due to lightning activity. In many of the Florida storms deposition of NO occurred up to near the tropopause but major deposition usually occurred 1 - 2 km below the tropopause, or mostly within the visible anvil volume formed prior to storm decay. Observations from two storms of very different anvil size and electrical activity allowed estimates of the total mass of NO, vented to the middle and upper troposphere. Using the cloud-to ground (CG) flash accumulations from the National Lightning Detection Network, climatological intra-cloud (IC) to CG ratios, and assuming that CG and IC flashes were of equivalent efficiency for NO production, the ranges of production per flash for a moderate-sized and a large storm were (0.51 - 1.0) x l0(exp 26) and (2.3 - 3.1) x 10(exp 26) molecules NO/flash, respectively. Using the recently determined average global flash rate of 44 8, a gross extrapolation of these two storms to represent possible global annual production rates yield 1.6 - 3.2 and 7.3 - 9.9 Tg(N)/yr, respectively. If the more usual assumption is made that IC efficiency is l/l0th that of CG activity, the ranges of production for the moderate-sized and large storm were (1.3 - 2.7) x l0(exp 26) and (6.0 - 8.1) x l0(exp 26) molecules NO/CG flash, respectively. The estimates from the large storm may be high because there is indirect evidence that the IC/CG ratio was larger than would be derived from climatology. These two storms and others studied did not have flash rates that scaled as approx. H(sup 5) where H is the cloud top altitude. The observed CG flash accumulations and NO(x) mass production estimate for the month of July over the Florida area were compared with a representative 3D global Chemistry-Transport Model (CTMJ that uses the Price et al. lightning parameterization. For two land grid points representing the Florida peninsula the model compared well with the observations: CG flash rates were low by only a factor of approx. 2. When the model grid points included the coastal regions of Florida the flash accumulations were lower than observed by a factor of 3.4 - 4.6. It is recommended that models using the Price et al. parameterization allow any global coastal grid point to maintain the land rather than the marine flash rate parameterization. The convection in this CTM underestimated the actual cloud top heights over Florida by 1 - 2 km and thus the total lightning flash rates and the altitude range of reactive nitrogen deposition. Broad scale (20 - 120 km) median mixing ratios of NO within anvils over Florida were significantly larger than in storms previously investigated over Colorado and New Mexico.

Ridley, B.↗
NASA NTRS1999

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗
NASA NTRS2010

Reevaluating the Role of Saharan Air Layer in Atlantic Tropical Cyclogenesis and Evolution

The existence of the Saharan air layer (SAL), a layer of warm, dry, dusty air that frequently moves westward off of the Saharan desert of Africa and over the tropical Atlantic Ocean, has long been appreciated. As air moves over the desert, it is strongly heated from below, producing a very hot air mass at low levels. Because there is no moisture source over the Sahara, the rise in temperature causes a sharp drop in relative humidity, thus drying the air. In addition, the warm air produces a very strong jet of easterly flow in the middle troposphere called the African easterly jet that is thought to play a critical role in hurricane formation. In recent years, there has been an increased focus on the impact that the SAL has on the formation and evolution of hurricanes in the Atlantic. However, the nature of its impact remains unclear, with some researchers arguing that the SAL amplifies hurricane development and with others arguing that it inhibits it. The argument for positively influencing hurricane development is based upon the fact that the African easterly jet provides an energy source for the waves that eventually form hurricanes and that it leads to rising motion south of the jet that favors the development of deep thunderstorm clouds. The potential negative impacts of the SAL include 1) low-level vertical wind shear associated with the African easterly jet; 2) warm SAL air aloft, which increases thermodynamic stability and suppresses cloud development; and 3) dry air, which produces cold downdrafts in precipitating regions, thereby removing energy needed for storm development. As part of this recent focus on the SAL and hurricanes (which motivated a 2006 NASA field experiment), there has been little emphasis on the SAL s potential positive influences and almost complete emphasis on its possible negative influences, almost to the point of claims that the SAL is the major suppressing influence on hurricanes in the Atlantic. In this study, multiple NASA satellite data sets (TRMM, MODIS, CALIPSO, and AIRS/AMSU) and National Centers for Environmental Prediction global analyses are used to see if the proposed negative influences deserve all of the attention they have recently received. The results show that storms generally form on the southern side of the African jet, where favorable background rotation is high. The jet often helps to form the northern side of the storms and is typically stronger in storms that intensify than those that weaken, suggesting that jet-induced vertical wind shear is not a negative influence on developing storms. Warm SAL air is confined to regions north of the jet and generally does not impact the tropical cyclone precipitation south of the jet. A comparison of the environments of strongly strengthening storms and of weakening storms shows no differences in SAL structure, indicating that the SAL has little influence on whether storms weaken or intensify. The large-scale flow at upper levels above the SAL was found to be most important, with the environment of strengthening storms having very little vertical wind shear and also favoring more expansive outflow from the storm. The SAL is shown to occur in a large-scale environment that is already characteristically dry as a result of large-scale subsidence (sinking air motions). Strong surface heating and deep dry convective mixing enhance dryness at low levels, but moisten the air at midlevels. Therefore, mid-to-upper-level dryness is not a defining characteristic of the SAL, but is instead a signature of subsidence. As a result, we conclude that the SAL is not the major negative influence on hurricanes that recent studies have emphasized. It is just one of many possible influences and can be both positive and negative.

Braun, Scott A.↗