Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW - PLATES”

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 559 records · Page 31

Constraints on hydrothermal heat flux through the oceanic lithosphere from global heat flow

A significant discrepancy exists between the heat flow measured at the seafloor and the higher values predicted by thermal models of the cooling lithosphere. This discrepancy is generally interpreted as indicating that the upper oceanic crust is cooled significantly by hydrothermal circulation. The magnitude of this heat flow discrepancy is the primary datum used to estimate the volume of hydrothermal flow, and the variation in the discrepancy with lithospheric age is the primary constraint on how the hydrothermal flux is divided between near-ridge and off-ridge environments. The resulting estimates are important for investigation of both the thermal structure of the lithosphere and the chemistry of the oceans. We reevaluate the magnitude and age variation of the discrepancy using a global heat flow data set substantially larger than in earlier studies, and the GDHI (Global Depth and Heat Flow) model that better predicts the heat flow. We estimate that of the predicted global oceanic heat flux of 32 x 10(exp 12) W, 34% (11 x 10(exp 12) W) occurs by hydrothermal flow. Approximately 30% of the hydrothermal heat flux occurs in crust younger than 1 Ma, so the majority of this flux is off-ridge. These hydrothermal heat flux estimates are upper bounds, because heat flow measurements require sediment at the site and so are made preferentially at topographic lows, where heat flow may be depressed. Because the water temperature for the near-ridge flow exceeds that for the off-ridge flow, the near-ridge water flow will be even a smaller fraction of the total water flow. As a result, in estimating fluxes from geochemical data, use of the high water temperatures appropriate for the ridge axis may significantly overestimate the heat flux for an assumed water flux or underestimate the water flux for an assumed heat flux. Our data also permit improved estimates of the 'sealing' age, defined as the age where the observed heat flow approximately equals that predicted, suggesting that hydrothermal heat transfer has largely ceased. Although earlier studies suggested major differences in sealing ages for different ocean basins, we find that the sealing ages for the Atlantic, Pacific, and Indian oceans are similar and consistent with the sealing age for the entire data set, 65 +/- 10 Ma. The previous inference of a young (approximately 20 Ma) sealing age for the Pacific appears to have biased downward several previous estimates of the global hydrothermal flux. The heat flow data also provide indirect evidence for the mechanism by which the hydrothermal heat flux becomes small, which has often been ascribed to isolation of the igneous crust from seawater due to the hydraulic conductivity of the intervening sediment. We find, however, that even the least sedimented sites show the systematic increase of the ratio of observed to predicted heat flow with age, although the more sedimented sites have a younger sealing age. Moreover, the heat flow discrepancy persists at heavily sedimented sites until approximately 50 Ma. It thus appears that approximately 100-200 m of sediment is neither necessary nor sufficient to stop hydrothermal heat transfer. We therefore conclude that the age of the crust is the primary control on the fraction of heat transported by hydrothermal flow and that sediment thickness has a lesser effect. This inference is consistent with models in which hydrothermal flow decreases with age due to reduced crustal porosity and hence permeability.

Stein, Carol A.↗

Global Qualitative Flow-Path Modeling for Local State Determination in Simulation and Analysis

For qualitative modeling and analysis, a general qualitative abstraction of power transmission variables (flow and effort) for elements of flow paths includes information on resistance, net flow, permissible directions of flow, and qualitative potential is discussed. Each type of component model has flow-related variables and an associated internal flow map, connected into an overall flow network of the system. For storage devices, the implicit power transfer to the environment is represented by "virtual" circuits that include an environmental junction. A heterogeneous aggregation method simplifies the path structure. A method determines global flow-path changes during dynamic simulation and analysis, and identifies corresponding local flow state changes that are effects of global configuration changes. Flow-path determination is triggered by any change in a flow-related device variable in a simulation or analysis. Components (path elements) that may be affected are identified, and flow-related attributes favoring flow in the two possible directions are collected for each of them. Next, flow-related attributes are determined for each affected path element, based on possibly conflicting indications of flow direction. Spurious qualitative ambiguities are minimized by using relative magnitudes and permissible directions of flow, and by favoring flow sources over effort sources when comparing flow tendencies. The results are output to local flow states of affected components.

Malin, Jane T.↗

Inhibition of the active lymph pump by flow in rat mesenteric lymphatics and thoracic duct

There are only a few reports of the influence of imposed flow on an active lymph pump under conditions of controlled intraluminal pressure. Thus, the mechanisms are not clearly defined. Rat mesenteric lymphatics and thoracic ducts were isolated, cannulated and pressurized. Input and output pressures were adjusted to impose various flows. Lymphatic systolic and diastolic diameters were measured and used to determine contraction frequency and pump flow indices. Imposed flow inhibited the active lymph pump in both mesenteric lymphatics and in the thoracic duct. The active pump of the thoracic duct appeared more sensitive to flow than did the active pump of the mesenteric lymphatics. Imposed flow reduced the frequency and amplitude of the contractions and accordingly the active pump flow. Flow-induced inhibition of the active lymph pump followed two temporal patterns. The first pattern was a rapidly developing inhibition of contraction frequency. Upon imposition of flow, the contraction frequency immediately fell and then partially recovered over time during continued flow. This effect was dependent on the magnitude of imposed flow, but did not depend on the direction of flow. The effect also depended upon the rate of change in the direction of flow. The second pattern was a slowly developing reduction of the amplitude of the lymphatic contractions, which increased over time during continued flow. The inhibition of contraction amplitude was dependent on the direction of the imposed flow, but independent of the magnitude of flow. Nitric oxide was partly but not completely responsible for the influence of flow on the mesenteric lymph pump. Exposure to NO mimicked the effects of flow, and inhibition of the NO synthase by N (G)-monomethyl-L-arginine attenuated but did not completely abolish the effects of flow.

NASA Discipline Musculoskeletal↗

Gas Flow Detection System

This system provides a portable means to detect gas flow through a thin-walled tube without breaking into the tubing system. The flow detection system was specifically designed to detect flow through two parallel branches of a manifold with only one inlet and outlet, and is a means for verifying a space shuttle program requirement that saves time and reduces the risk of flight hardware damage compared to the current means of requirement verification. The prototype Purge Vent and Drain Window Cavity Conditioning System (PVD WCCS) Flow Detection System consists of a heater and a temperature-sensing thermistor attached to a piece of Velcro to be attached to each branch of a WCCS manifold for the duration of the requirement verification test. The heaters and thermistors are connected to a shielded cable and then to an electronics enclosure, which contains the power supplies, relays, and circuit board to provide power, signal conditioning, and control. The electronics enclosure is then connected to a commercial data acquisition box to provide analog to digital conversion as well as digital control. This data acquisition box is then connected to a commercial laptop running a custom application created using National Instruments LabVIEW. The operation of the PVD WCCS Flow Detection System consists of first attaching a heater/thermistor assembly to each of the two branches of one manifold while there is no flow through the manifold. Next, the software application running on the laptop is used to turn on the heaters and to monitor the manifold branch temperatures. When the system has reached thermal equilibrium, the software application s graphical user interface (GUI) will indicate that the branch temperatures are stable. The operator can then physically open the flow control valve to initiate the test flow of gaseous nitrogen (GN2) through the manifold. Next, the software user interface will be monitored for stable temperature indications when the system is again at thermal equilibrium with the test flow of GN2. The temperature drop of each branch from its "no flow" stable temperature peak to its stable "with flow" temperature will allow the operator to determine whether a minimum level of flow exists. An alternative operation has the operator turning on the software only long enough to record the ambient temperature of the tubing before turning on the heaters and initiating GN2 flow. The stable temperature of the heated tubing with GN2 flow is then compared with the ambient tubing temperature to determine if flow is present in each branch. To help quantify the level of flow in the manifolds, each branch will be bench calibrated to establish its thermal properties using the flow detection system and different flow rates. These calibration values can then be incorporated into the software application to provide more detailed flow rate information.

Moss, Thomas↗

Hydrodynamics of Two-Phase Flows through Porous Media in Microgravity: Packed Bed Reactor Experiment onboard of the International Space Station

The objective of the Packed Bed Reactor Experiment was to develop a fundamental understanding of the hydrodynamics of two-phase flow with no phase transition through porous media in microgravity. This work describes two experimental campaigns flown on the International Space Station (ISS) that were designed to achieve the objective. This work presents the results on flow patterns observed, two-phase flow pressure drop in porous media, and the impact of capillary effects on pressure drop at various gas and liquid flow rates. This work is the first to present predictive correlations of the two-phase friction factor for the different regimes identified based on the superficial liquid and gas velocity in microgravity, the first to hypothesize the different regimes based on the change in slope of the pressure gradient versus flow rate plots, the first to address the effects of the capillarity on the pressure gradient and the first to measure column holdup and assess the effective porosity based on the packing and the trapped gas bubbles. Experimental data on pressure drop in gas (N2) – liquid (water) flow show the presence of four different flow regimes in microgravity. Dispersed bubble flow (DB) and pulse (P) flow regimes are detected at high liquid flow rates, whereas at low liquid flow rates “large bubble” or elongated bubble (LB) and “gas channeling” (GC) regimes are observed. For these different flow regimes, different two-phase flow friction factor and pressure gradient correlations are presented for the first time as a function of the gas and liquid modified Reynolds numbers and Suratman number. Within the viscous-capillary (V-C) regime, it is found that the capillary contribution is the dominant force that contributes to the pressure drop for the wetting case (glass). However, for the non-wetting packing (Teflon), the viscous contribution dominates. It was found that the gas hold-up and pressure drop are functions of bed history at low liquid and gas flow rates with the magnitude of the hysteresis decreasing with increasing flow rates. PBRE-2 results show that the capillary force is a strong function of the superficial liquid velocity but is a much weaker function of the superficial gas velocity and varies inversely with the particle diameter. Within the Viscous-Capillary (V-C) regime, over 90% of the pressure gradient is attributed to the capillary contribution in the gas continuous regime. However, in the large bubble regime, the viscous and capillary contributions were comparable. After the completion of PBRE-2 with glass spherical beads, another bed packed with alumina was installed. Pressure gradient data for the alumina packed bed (PBRE-Water Recovery) were transferred to the PI and were not analyzed by the authors of this work. The Packed Bed Reactor Experiment concluded its on-orbit operation after two successful campaigns of testing in 2017 and in 2021. The flight hardware was brought back from orbit and is being reconditioned for a series of future experiments referred to as PBRE-WRS in support of water recovery (WRS). The WRS series consists of testing different two-phase fluid system components and packed beds to assess their pressure gradient characteristics, which will be used for designing packed bed reactors for various applications relevant to life support. Although there is so much relevance of the experimental results obtained from the two PBRE campaign in microgravity to life support, these findings also apply to other applications that involve two phase flow in porous medium such as fuel cells, transport of nutrient to plants in space and other chemical and materials processing that involve two phase flows. These systems operate differently in microgravity because, due to the lack of buoyancy, the density difference between the phases becomes irrelevant and no longer leads to phase separation.

Packed Bed↗

Studies on Normal and Microgravity Annular Two Phase Flows

Two-phase gas-liquid flows occur in a wide variety of situations. In addition to normal gravity applications, such flows may occur in space operations such as active thermal control systems, power cycles, and storage and transfer of cryogenic fluids. Various flow patterns exhibiting characteristic spatial and temporal distribution of the two phases are observed in two-phase flows. The magnitude and orientation of gravity with respect to the flow has a strong impact on the flow patterns observed and on their boundaries. The identification of the flow pattern of a flow is somewhat subjective. The same two-phase flow (especially near a flow pattern transition boundary) may be categorized differently by different researchers. Two-phase flow patterns are somewhat simplified in microgravity, where only three flow patterns (bubble, slug and annular) have been observed. Annular flow is obtained for a wide range of gas and liquid flow rates, and it is expected to occur in many situations under microgravity conditions. Slug flow needs to be avoided, because vibrations caused by slugs result in unwanted accelerations. Therefore, it is important to be able to accurately predict the flow pattern which exists under given operating conditions. It is known that the wavy liquid film in annular flow has a profound influence on the transfer of momentum and heat between the phases. Thus, an understanding of the characteristics of the wavy film is essential for developing accurate correlations. In this work, we review our recent results on flow pattern transitions and wavy films in microgravity.

Balakotaiah, V.↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗

Thin-Film Air-Mass-Flow Sensor of Improved Design Developed

Researchers at the NASA Glenn Research Center have developed a new air-mass-flow sensor to solve the problems of existing mass flow sensor designs. NASA's design consists of thin-film resistors in a Wheatstone bridge arrangement. The resistors are fabricated on a thin, constant-thickness airfoil to minimize disturbance to the airflow being measured. The following photograph shows one of NASA s prototype sensors. In comparison to other air-mass-flow sensor designs, NASA s thin-film sensor is much more robust than hot wires, causes less airflow disturbance than pitot tubes, is more accurate than vane anemometers, and is much simpler to operate than thermocouple rakes. NASA s thin-film air-mass-flow sensor works by converting the temperature difference seen at each leg of the thin-film Wheatstone bridge into a mass-flow rate. The following figure shows a schematic of this sensor with air flowing around it. The sensor operates as follows: current is applied to the bridge, which increases its temperature. If there is no flow, all the arms are heated equally, the bridge remains in balance, and there is no signal. If there is flow, the air passing over the upstream legs of the bridge reduces the temperature of the upstream legs and that leads to reduced electrical resistance for those legs. After the air has picked up heat from the upstream legs, it continues and passes over the downstream legs of the bridge. The heated air raises the temperature of these legs, increasing their electrical resistance. The resistance difference between the upstream and downstream legs unbalances the bridge, causing a voltage difference that can be amplified and calibrated to the airflow rate. Separate sensors mounted on the airfoil measure the temperature of the airflow, which is used to complete the calculation for the mass of air passing by the sensor. A current application for air-mass-flow sensors is as part of the intake system for an internal combustion engine. A mass-flow sensor is used to provide accurate information about the amount of air entering the engine so that the amount of fuel can be adjusted to give the most efficient combustion. The ideal mass-flow sensor would be a rugged design that minimizes the disturbance to the flow stream and provides an accurate reading of both smooth and turbulent flows; NASA's design satisfies these requirements better than any existing design. Most of the mass-flow sensors used today are the hot wire variety. Hot wires can be fragile and cannot accurately measure a turbulent or reversing flow, which is often encountered in an intake manifold. Other types of mass-flow sensors include pitot tubes, vane anemometers, and thermocouple rakes-all of which suffer from some type of performance problem. Because it solves these performance problems while maintaining a simple design that lends itself to low-cost manufacturing techniques, NASA s thin-film resistance temperature detector air-mass-flow sensor should lead to more widespread use of mass-flow sensors.

Fralick, Gustave C.↗

Terrestrial Analogue Studies from ISRO’s Venus Mission Perspective: Polarimetric Radar Properties of Hawaiian Lava Flows

Detection of present and past volcanism on Venus is one of the major goals of the proposed ISRO’s Venus orbiter mission. The S-band, high-resolution (40 m/pixel) fully Polarimetric Synthetic Aperture Radar (PolSAR) instrument on this mission [1] may have the capabilities to detect volcanism within the mission lifetime using repeated PolSAR imagery, and possibly SAR Interferometry (InSAR, experimental mode). The global mapping efforts of Venus by Magellan mission has enabled comprehensive mapping of lava flows and indicated that Venus has an extensive history of volcanism, the ages of which are largely unknown (e.g. [2]). While Magellan radar emissivity data (e.g. [3, 4]) and Near Infrared emissivity data from VIRTIS [5] suggested the occurrence of recent episodes of volcanic activity (e.g. Maat Mons and Ganis Chasma), as well as the presence of stratigraphically young lava flows (e.g. Idunn Mons), measuring changes in radar backscatter amplitude alone to identify lava flows has many challenges. Although very large changes in the shape of the terrain can be observed in radar backscatter amplitude changes (e.g. [6]), smaller, or relatively flat lava flows are difficult to detect. Previous terrestrial studies suggest that PolSAR and InSAR techniques are very effective for mapping lava flows (e.g. [7-9]), and can be used when changes cannot be distinguished in radar backscatter images. We use the unvegetated lava flows on Hawaiʻi island as a terrestrial analogue to study Venus lava flows for the following reasons: (a) It is extensively studied at several wavelengths commonly used in remote sensing studies (including PolSAR and InSAR methods); and (b) it is a volcanically active area with new lava flows frequently covering older emplaced flows. To investigate the surface roughness, texture, and fine-grained mantling associated with Mauna Loa and Kilauea lava flows, we utilize C- and L-band PolSAR datasets obtained from RISAT-1A (EOS-4) and ALOS PALSAR missions respectively. In particular, we will use the quad-polarized backscatter and polarimetric parameters to characterize the texture of the terrestrial lava flows to understand whether the Venus crust is continuously disrupted during flow emplacement. While some previous studies (e.g. [10, 11]) suggested that surface roughness of most of the Venus flows is comparable to that of terrestrial pāhoehoe flows, other studies indicated that fractal dimensions of some large lava flows on Venus imply high eruption rates which favour the formation of a’a flows (e.g. [12]). We will also analyse the terrestrial flows for the presence/absence of pyroclastic mantling as radar-bright diffuse deposits near the summit regions of some coronae on Venus have been proposed to be young pyroclastics, and possible evidence of a renewed epoch of mantle volcanism that taps into deeper volatiles [13]. A recent study using EOS-4 RISAT-1 data of a part of fresh Mauna Loa lava flows (2022 eruption) emphasizes the ability of fully polarimetric SAR data to understand the diversity of physical properties (e.g. texture and morphology) associated with them (Sreejith et al. 2024); and we will apply similar methods to the PolSAR data obtained from ISRO’s Venus mission for our proposed objectives.

Sriram S Bhiravarasu↗

Extensive lava flow fields on Venus: Preliminary investigation of source elevation and regional slope variations

Large-volume lava flow fields have been identified on Venus, the most areally extensive of which are known as fluctus and have been subdivided into six morphologic types. Sheetlike flow fields (Type 1) lack the numerous, closely spaced, discrete lava flow lobes that characterize digitate flow fields. Transitional flow fields (Type 2) are similar to sheetlike flow fields but contain one or more broad flow lobes. Digitate flow fields are divided further into divergent (Types 3-5) and subparallel (Type 6) classes on the basis of variations in the amount of downstream flow divergence. As a result of our previous analysis of the detailed morphology, stratigraphy, and tectonic associations of Mylitta Fluctus, we have formulated a number of questions to apply to all large flow fields on Venus. In particular, we would like to address the following: (1) eruption conditions and style of flow emplacement (effusion rate, eruption duration), (2) the nature of magma storage zones (presence of neutral buoyancy zones, deep or shallow crustal magma chambers), (3) the origin of melt and possible link to mantle plumes, and (4) the importance of large flow fields in plains evolution. To answer these questions we have begun to examine variations in flow field dimension and morphology; the distribution of large flow fields in terms of elevation above the mean planetary radius; links to regional tectonic or volcanic structures (e.g., associations with large shield edifices, coronae, or rift zones); statigraphic relationships between large flow fields, volcanic plains, shields, and coronae; and various models of flow emplacement in order to estimate eruption parameters. In this particular study, we have examined the proximal elevations and topographic slopes of 16 of the most distinctive flow fields that represent each of the 6 morphologic types.

Magee-Roberts, K.↗

Gravitational Effects on Flow Instability and Transition in Low Density Jets

Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Agrawal A. K.↗

The Summer 1997 Eruption at Pillan Patera on Io: Implications for Ultrabasic Lava Flow Emplacement

Galileo data and numerical modeling were used to investigate the summer 1997 eruption at Pillan Patera on Io. This event, now defined as 'Pillanian' eruption style, included a high-temperature (greater than 1600 C), possibly ultrabasic, 140-km-high plume eruption that deposited dark, orthopyroxene-rich pyroclastic material over greater than 125,000 sq km, followed by emplacement of dark flow-like material over greater than 3100 sq km to the north of the caldera. We estimate that the high-temperature, energetic episode of this eruption had a duration of 52- 167 days between May and September 1997, with peak eruption temperatures around June 28, 1997. Galileo 20 m/pixel images of part of the Pillan flow field show a widespread, rough, pitted surface that is unlike any flow surface we have seen before. We suggest that th.s surface may have resulted from (1) a fractured lava crust formed during rapid, low-viscosity lava surging, perhaps including turbulent flow emplacement; (2) disruption of the lava flow by explosive interaction with a volatile-rich substrate: or (3) a combination of 1 and 2 with or without accumulation of pyroclastic materials on the surface. Well-developed flow lobes are observed, suggesting that this is a relatively distal part of the flow field. Shadow measurements at flow margins indicate a thickness of approx. 8-10 m. We have modeled the emplacement of putative ultrabasic flows from the summer 1997 Pillan eruption using constraints from new Galileo data. Results suggest that either laminar sheet flows or turbulent channelized flows could have traveled 50-150 km on a flat. unobstructed surface, which is consistent with the estimated length of the Pillan flow field (approx. 60 km). Our modeling suggests low thermal erosion rates (less than 0.1 m/d), and that the formation of deep (greater than 20 m) erosion channels was unlikely, especially distal to the source. We calculate a volumetric flow rate of approx. 2-7 x l0(exp 3) cu m/s, which is greater than those for typical Mauna Loa/Kilauea flows but comparable to those for the (1783) Laki eruption and the inferred flow rates of the Roza flows in the Columbia River flood basalts. The differences in ultrabasic eruption styles on Earth and Io appear to be controlled by the different eruption environments: Plumes at sites of ultrabasic eruptions on Io suggest strong magma-volatile: interactions on a low-gravity body lacking an atmosphere, whereas the geology at sites of komatiite eruptions on Earth suggest mostly submarine emplacement of thick flows with a pronounced lack of subaerial explosive activity.

Williams, David A.↗

The Summer 1997 Eruption at Pillan Patera on Io: Implications for Ultrabasic Lava Flow Emplacement

Galileo data and numerical modeling were used to investigate the summer 1977 eruption at Pillan Patera on Io. This event, now defined as "Pillanian" eruption style, included a high-temperature (greater than 1600 C), possible ultrabasic , 140-km-high plume eruption that deposited dark, orthopyroxene-rich pyroclastic material over greater than 125,000 sq km, followed by emplacement of dark flow-like material over greater than 3100 sq km to the north of the caldera. We estimate that the high-temperature, energetic episode of this eruption had a duration of 52 - 167 days between May and September 1997, with peak eruption temperatures around June 28, 1997. Galileo 20 m/pixel images of part of the Pillan flow field show a wide-spread, rough, pitted surface that is unlike any flow surface we have seen before. We suggest that this surface may have resulted from: 1. A fractured lava crust formed during rapid, low-viscosity lava surging, perhaps including turbulent flow emplacement. 2. Disruption of the lava flow by explosive interaction with a volatile-rich substrate. or 3. A combination of 1 and 2 with or without accumulation of pyroclastic material on the surface. Well-developed flow lobes are observed, suggesting that this is a relatively distant part of the flow field.Shadow measurements at flow margins indicate a thickness of-8 - 10 m. We have modeled the emplacement of putative ultrabasic flow from the summer 1997 Pillan eruption using constraints from new Galileo data. Results suggest that either laminar sheet flows or turbulent channelized flows could have traveled 50 - 150 km on a flat, unobstructed surface, which is consistent with the estimated length of the Pillan flow field (approx. 60 km). Our modeling suggests low thermal erosion rates (less than 4.1 m/d), and that the formation of deep (greater than 20 m) erosion channels was unlikely, especially distal to the source. We calculate a volumetric flow rate of approx. 2 - 7 x 10(exp 3)cu m/s, which is greater than those for typical Mauna Loa/Kilaueaq flows but comparable to those for the (1783) Laki eruption and the inferred flow rates of the Roza flows in the Columbia River flood basalts. The differences in ultrabasic eruption styles on Earth and Io appear to be controlled by the different eruption environments; Plumes at sites of ultrabasic eruptions on Io suggest strong magma-volatile interactions on a low-gravity body lacking an atmosphere, whereas the geology at sites of komatiite eruptions on Earth suggest mostly submarine emplacement of thick flows with a pronounced lack of subaerial explosive activity.

Williams, David A.↗

Analytical and numerical studies of the thermocapillary flow in a uniformly floating zone

The microgravity environment of an orbiting vehicle permits crystal growth experiments in the presence of greatly reduced buoyant convection in the liquid melt. Crystals grown in ground-based laboratories do not achieve their potential properties because of dopant variations caused by flow in the melt. The floating zone crystal growing system is widely used to produce crystals of silicon and other materials. However, in this system the temperature gradient on the free sidewall surface of the melt is the source of a thermocapillary flow which does not disappear in the low-gravity environment. The idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow to the melt sidewall leaving the interior of the melt passive is examined. A cylinder of fluid with an axial temperature gradient imposed on the cylindrical sidewall is considered. A half zone and the linearized, axisymmetric flow in the absence of crystal growth is examined. Rotation is found to confine the linear thermocapillary flow. A simplified model is extended to a full zone and both linear and nonlinear thermocapillary flows are studied theoretically. Analytical and numerical methods are used for the linear flows and numerical methods for the nonlinear flows. It was found that the linear flows in the full zone have more complicated and thicker boundary layer structures than in the half zone, and that these flows are also confined by the rotation. However, for the simplified model considered and for realistic values for silicon, the thermocapillary flow is not linear. The fully nonlinear flow is strong and unsteady (a weak oscillation is present) and it penetrates the interior. Some non-rotating flow results are also presented. Since silicon as a large value of thermal conductivity, one would expect the temperature fields to be determined by conduction alone. This is true for the linear and weakly nonlinear flows, but for the stronger nonlinear flow the results show that temperature advection is also important. Uniform rotation may still be a means of confining the flow and the results obtained define the procedure to be used to examine this hypothesis.

Fowlis, W. W.↗

[Low-Frequency Flow Oscillation

The results of the research conducted under this grant are presented in detail in three Master theses, by Heinrich, Balow, and Broeren. Additional analysis of the experimental data can be found in two AIAA Journal articles and two conference papers. Citations for all of the studies' publications can be found in the bibliography which is attached. The objective of Heinrich's study was to document the low-frequency flow oscillation on the LRN-1007 airfoil, which had been previously observed at low Reynolds number, to determine its origin, and explore the phenomenon at higher Reynolds number. Heinrich performed detailed flow visualization on the airfoil using surface fluorescent oil and laser-sheet off-body visualization. A large leading-edge separation bubble and trailing-edge separation was identified on the airfoil just prior to the onset of the unsteady stall flow oscillation. From the laser-sheet data, the unsteady flow appeared as a massive boundary-layer separation followed by flow reattachment. Hot-wire data were taken in the wake to identify the presence of the flow oscillation and the dominant frequency. The oscillation was found in the flow from a Reynolds number of 0.3 to 1.3 x 10 exp 6. The Strouhal number based on airfoil projected height was nominally 0.02 and increased slightly with increasing Reynolds number and significantly with increasing airfoil angle of attack. Balow focused his research on the leading-edge separation bubble which was hypothesized to be the origin of the low-frequency oscillation. Initially, experimental measurements in the bubble at the onset of the low-frequency oscillation were attempted to study the characteristics of the bubble and explain possible relationships to the shear-layer-flapping phenomena. Unfortunately, the bubble proved to be extremely sensitive to the probe interference and it drastically reduced the size of the bubble. These detailed measurements were then abandoned by Balow. However, this led to a series of tests where the leading-edge bubble and trailing-edge separation were altered and the affect on the flow-oscillation studied. Balow found that by tripping the airfoil boundary-layer with "zigzag" tape ahead of bubble separation, the bubble was effectively eliminated mid the oscillation suppressed. Wake survey drag measurements showed a drastic reduction in airfoil drag when the bubble and oscillation were eliminated. Using the "zigzag" tape, the trailing-edge separation was moved downstream approximately 5 percent chord. This was found to reduce the amplitude of the oscillation, particularly in the onset stage at low angle of attack (around 14 degrees). Through detailed analysis of the wake behind the airfoil during the unsteady flow oscillation, Balow provided a better understanding of the wake flowfield. Broeren studied the oscillating flowfield in detail at Reynolds number equal 3 x 10 exp 5 and an angle of attack of 15 degrees using laser Doppler velocimetry (LDV). Two-dimensional LDV data were acquired at 687 grid points above the model upper surface while hot-wire data were taken simultaneously in the wake. Using the hot-wire signal, the LDV data were phase averaged into 24 bins to represent a single ensemble average of one oscillation cycle. The velocity data showed a flowfield oscillation that could be divided into three flow regimes. In the first regime, the flow over the airfoil was completely separated initially, the flowfield reattached from the leading edge and the reattachment point moved downstream with increasing time or phase. Broeren referred to this as the reattachment regime. The bubble development regime followed, where a leading-edge separation bubble formed at the leading edge and grew with increasing time. During the initial part of this regime the trailing-edge separation continued to move downstream. However, during the last 30 degrees of phase the trailing-edge separation moved rapidly forward and appeared to merge with the leading-edge bubble. During the third regime, the separation regime, the flow was segmented from the airfoil leading edge and did not reattach to the airfoil surface. The reverse flow was seen to grow in vertical extent up from the model surface as the phase increased. Next reattachment began again at the leading edge signaling the start of the reattachment regime, and so the cycle continued. From Broeren's work, the details of the unsteady flowfield over the airfoil were seen for the first time. From this research a great deal has been learned about the low-frequency flow oscillation which naturally occurs on the LRN-1007 airfoil near stall. The oscillation was seen to persist at higher Reynolds number, the dependence of the Strouhal number on angle of attack and Reynolds number were discovered, the critical role played by the laminar bubble was shown and the entire upper surface flowfield during a flow oscillation cycle was measured and analyzed. What still eludes understanding is the scaling of the flow oscillation and why certain airfoils, such as the LRN, have a very strong low-frequency mode and other airfoils exhibit no organized low-frequency oscillation at all.

Bragg, Michael B.↗

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan Brian G.↗

Criterion for Identifying Vortices in High-Pressure Flows

A study of four previously published computational criteria for identifying vortices in high-pressure flows has led to the selection of one of them as the best. This development can be expected to contribute to understanding of high-pressure flows, which occur in diverse settings, including diesel, gas turbine, and rocket engines and the atmospheres of Jupiter and other large gaseous planets. Information on the atmospheres of gaseous planets consists mainly of visual and thermal images of the flows over the planets. Also, validation of recently proposed computational models of high-pressure flows entails comparison with measurements, which are mainly of visual nature. Heretofore, the interpretation of images of high-pressure flows to identify vortices has been based on experience with low-pressure flows. However, high-pressure flows have features distinct from those of low-pressure flows, particularly in regions of high pressure gradient magnitude caused by dynamic turbulent effects and by thermodynamic mixing of chemical species. Therefore, interpretations based on low-pressure behavior may lead to misidentification of vortices and other flow structures in high-pressure flows. The study reported here was performed in recognition of the need for one or more quantitative criteria for identifying coherent flow structures - especially vortices - from previously generated flow-field data, to complement or supersede the determination of flow structures by visual inspection of instantaneous fields or flow animations. The focus in the study was on correlating visible images of flow features with various quantities computed from flow-field data.

Bellan, Josette↗

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗