ALBERT

Beschreibung: A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65 delta wing with different values of leading-edge radius. Although the geometry is quite simple, it poses a challenging problem for computing vortices originating from blunt leading edges. The second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the wind-tunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.

Beschreibung: Extensive pressure measurements and off-surface flow visualization were obtained on the forebody and strakes of the NASA F-18 High Alpha Research Vehicle (HARV) equipped with actuated forebody strakes. Forebody yawing moments were obtained by integrating the circumferential pressures on the forebody and strakes. Results show that large yawing moments can be generated with forebody strakes. At a 50 -angle-of-attack, deflecting one strake at a time resulted in a forebody yawing moment control reversal for small strake deflection angles. However, deflecting the strakes differentially about a 20 symmetric strake deployment eliminated the control reversal and produced a near linear variation of forebody yawing moment with differential strake deflection. At an angle of attack of 50 and for 0 and 20 symmetric strake deployments, a larger forebody yawing moment was generated by the forward fuselage (between the radome and the apex of the leading-edge extensions) than on the radome where the actuated forebody strakes were located. Cutouts on the flight vehicle strakes that were not on the wind tunnel models are believed to be responsible for deficits in the suction peaks on the flight radome pressure distributions and differences in the forebody yawing moments.

Beschreibung: The CFD modeling used has produced reasonably good global upper-surface pressure coefficient comparisons with measured flight data at both transonic and subsonic speeds at the angles of attack presented. Boundary layer comparisons showed the profiles to be reasonably well predicted inboard and under the primary vortex system. However, the secondary vortex profile was not well predicted either at the anticipated separation point or under the secondary vortex. Moreover, the flight data showed there to be a vortex/boundary-layer interaction that occurred in the vicinity of the secondary vortex. The spanwise distribution of local skin friction measured data was reasonably well predicted, especially away from the wing leading-edge. Lastly, predicted and measured flight-pressures, as well as flight-image data, for the F-16XL-1 airplane are now available via the World Wide Web.

Beschreibung: Surface and off-surface flow visualization techniques were used to visualize the three-dimensional vortex flows on the F-106 aircraft with vortex flaps installed. Results at angles of attack between 9 degrees to 18 degrees and Mach numbers from 0.3 to 0.9 are presented. A smoke flow vapor screen technique was used to document leading-edge vortex paths and sizes, while an oil flow technique was employed to provide detailed information on reattachment and separation line locations and other flow details. Results were obtained for two vortex flap deflection angles, 30 degrees and 40 degrees. Flow visualization revealed the existence of a multiple vortex system that had not previously been seen in subscale tests or predicted for this configuration. The vortex flap generated a leading-edge vortex system that reattached near the flap hinge over a wide angle of attack range. In addition to the primary vortex, flow visualization revealed the presence of several distinct vortices which traced a path from the vortex flap and then over the wing.

Beschreibung: A wind tunnel experiment was conducted in the NASA Langley Research Center (LaRC) 8-Foot Transonic Pressure Tunnel (TPT) to determine the effects of passive surface porosity on vortex flow interactions about a general research fighter configuration at subsonic and transonic speeds. Flow-through porosity was applied to a wing leading-edge extension (LEX) mounted to a 65 deg cropped delta wing model to promote large nose-down pitching moment increments at high angles of attack. Porosity decreased the vorticity shed from the LEX, which weakened the LEX vortex and altered the global interactions of the LEX and wing vortices at high angles of attack. Six-component forces and moments and wing upper surface static pressure distributions were obtained at free-stream Mach numbers of 0.50, 0.85, and 1.20, Reynolds number of 2.5(10(exp 6)) per foot, angles of attack up to 30 deg, and angles of sideslip to +/- 8 deg. The off-surface flow field was visualized in selected cross-planes using a laser vapor screen flow visualization technique. Test data were obtained with a centerline vertical tail and with alternate twin, wing-mounted vertical fins having 0 deg and 30 deg cant angles. In addition, the porosity of the LEX was compartmentalized to determine the sensitivity of the vortex-dominated aerodynamics to the location and level of porosity applied to the LEX.

Beschreibung: The influences of slow convective flow on droplet combustion, particularly in the low Reynolds number regime, have received very little attention in the past. Most studies in the literature are semi-empirical in nature and they were motivated by spray combustion applications in the moderate to high Reynolds number regime. None of the limited number of fundamental theoretical studies applicable to low Reynolds numbers have been verified by rigorous experimental data. Moreover, many unsteady phenomena associated with fluid-dynamic unsteadiness, such as impulsive starting or stopping of a burning droplet, or flow acceleration/deceleration effects, have not been investigated despite their importance in practical applications. In this study we investigate the effects of slow convection on droplet burning dynamics both experimentally and theoretically. The experimental portion of the study involves both ground-based experiments in the drop towers and future flight experiments on board the International Space Station. Heptane and methanol are used as test fuels, and this choice complements the quiescent-environment studies of the Droplet Combustion Experiment (DCE). An analytical model that employs the method of matched asymptotic expansions and uses the ratio of the convective velocity far from the droplet to the Stefan velocity at its surface as the small parameter for expansion has also been developed as a part of this investigation. Results from the ground-based experiments and comparison with the analytical model are presented in this report.

Beschreibung: Important issues: Mass gauging; Stability dynamics of disconnected capillary surfaces; Slow capillary driven flow (i.e. wicking); Long-term material property evolution in micro-g; Dumping problem with freezing of dump lines.

Beschreibung: The objective of this viewgraph presentation is to present a summary of computational methods developed at Ames Research Center for large scale fluid/structure interaction.

Beschreibung: This paper presents viewgraphs on a multiphase flow panel. The topics include: 1) Discussion of Priorities; 2) Critical Issues Reduced Gravity Instabilities; 3) Severely Limiting Phase Separation; 4) Severely-Limiting Phase Change; 5) Enhancements; 6) Awareness Instabilities; 7) Awareness; 8) Methods of Resolution; 9) 2008 Space Flight; 10) 2003-2008 Ground-Based Microgravity Facilities; 11) 2003-2008 Other; 12) 2009-2015 Space Flight; 13) 2009-2015 Ground-Based Microgravity Facilities; 14) 2009-2015 Other; and 15) 2016.

Beschreibung: The short term purpose of this research is to present a research plan and a roadmap developed for strategic research for the Office of Biological and Physical Research and the long term purpose is to conduct necessary ground-based and space-flight low gravity experiments, complemented by analyses, resulting in a documented framework for parameter prediction of needed by designers. This paper is presented in viewgraph form.

Beschreibung: The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Beschreibung: During the past two summers Professor Milanovic conducted Wind tunnel experiments on steady jets-in-cross-flow and synthetic jets. In her anticipated visit during the upcoming summer, she will continue and complete the research on synthetic jets involving 2-dimensional orifices of different aspect ratio as well as inclined slots. In addition, experiments will be conducted on pulsatile jets-in-cross-flow. The pulsation will be provided via an oscillating valve at controllable frequencies. The experiment will involve mainly hot-wire anemometer measurements in the low-speed wind tunnel. Overall goal will be to obtain database and investigate flow control strategies. The research will be of fundamental nature.

Beschreibung: We demonstrate a new variation of molecular-tagging velocimetry for hypersonic flows based on laser-induced fluorescence. A thin line of nitric-oxide molecules is excited with a laser beam and then, after a time delay, a fluorescence image of the displaced line is acquired. One component of velocity is determined from the time of flight. This method is applied to measure the velocity profile in a Mach 8.5 laminar, hypersonic boundary layer in the Australian National University s T2 free-piston shock tunnel. The single-shot velocity measurement uncertainty in the freestream was found to be 3.5%, based on 90% confidence. The method is also demonstrated in the separated flow region forward of a blunt fin attached to a flat plate in a Mach 7.4 flow produced by the Australian National University s T3 free-piston shock tunnel. The measurement uncertainty in the blunt fin experiment is approximately 30%, owing mainly to low fluorescence intensities, which could be improved significantly in future experiments. This velocimetry method is applicable to very high-speed flows that have low collisional quenching of the fluorescing species. It is particularly convenient in facilities where planar laser-induced fluorescence is already being performed.

Beschreibung: The spreading of a liquid on a solid surface is important for various practical processes, and contact-angle measurements provide an elegant method to characterize the interfacial properties of the liquid with the solid substrates. The complex physical processes occurring when a liquid contacts a solid play an important role in determining the performance of chemical processes and materials. Applications for these processes are in printing, coating, gluing, textile dyeing, and adhesives and in the pharmaceutical industry, biomedical research, adhesives, flat panel display manufacturing, surfactant chemistry, and thermal engineering.

Beschreibung: The interface between two fluids of different density can experience instability when gravity acts normal to the surface. The relatively well known Rayleigh-Taylor (RT) instability results when the gravity is constant with a heavy fluid over a light fluid. An impulsive acceleration applied to the fluids results in the Richtmyer-Meshkov (RM) instability. The RM instability occurs regardless of the relative orientation of the heavy and light fluids. In many systems, the passing of a shock wave through the interface provides the impulsive acceleration. Both the RT and RM instabilities result in mixing at the interface. These instabilities arise in a diverse array of circumstances, including supernovas, oceans, supersonic combustion, and inertial confinement fusion (ICF). The area with the greatest current interest in RT and RM instabilities is ICF, which is an attempt to produce fusion energy for nuclear reactors from BB-sized pellets of deuterium and tritium. In the ICF experiments conducted so far, RM and RT instabilities have prevented the generation of net-positive energy. The $4 billion National Ignition Facility at Lawrence Livermore National Laboratory is being constructed to study these instabilities and to attempt to achieve net-positive yield in an ICF experiment.

Beschreibung: Lightweight, strong, tough high-temperature materials are required to complement efficiency improvements for next-generation gas turbine engines that can operate with minimum cooling. Because of their low density, high-temperature strength, and high thermal conductivity, ceramics are being investigated as materials to replace the nickelbase superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass (ref. 1). To complement the effectiveness of the ceramics and their applicability for turbine engine applications, a parametric study using the finite element method is being carried out. The NASA Glenn Research Center remains very active in conducting and supporting a variety of research activities related to ceramic matrix composites through both experimental and analytical efforts (ref. 1). The objectives of this work are to develop manufacturing technology, develop a thermal and environmental barrier coating (TBC/EBC), develop an analytical modeling capability to predict thermomechanical stresses, and perform a minimal burner rig test on silicon nitride (Si3N4) and SiC/SiC turbine nozzle vanes under simulated engine conditions. Moreover, we intend to generate a detailed database of the material s property characteristics and their effects on structural response. We expect to offer a wide range of data since the modeling will account for other variables, such as cooling channel geometry and spacing. Comprehensive analyses have begun on a plate specimen with Si3N4 cooling holes.

Beschreibung: A high-performance spontaneous Raman scattering (SRS) system for measuring quantitative species concentration and temperature in high-pressure flames is now operational. The system is located in Glenn s Engine Research Building. Raman scattering is perhaps the only optical diagnostic technique that permits the simultaneous (single-shot) measurement of all major species (N2, O2, CO2, H2O, CO, H2, and CH4) as well as temperature in combustion systems. The preliminary data acquired with this new system in a 20-atm hydrogen-air (H2-air) flame show excellent spectral coverage, good resolution, and a signal-to-noise ratio high enough for the data to serve as a calibration standard. This new SRS diagnostic system is used in conjunction with the newly developed High- Pressure Gaseous Burner facility (ref. 1). The main purpose of this diagnostic system and the High-Pressure Gaseous Burner facility is to acquire and establish a comprehensive Raman-scattering spectral database calibration standard for the combustion diagnostic community. A secondary purpose of the system is to provide actual measurements in standardized flames to validate computational combustion models. The High-Pressure Gaseous Burner facility and its associated SRS system will provide researchers throughout the world with new insights into flame conditions that simulate the environment inside the ultra-high-pressure-ratio combustion chambers of tomorrow s advanced aircraft engines.

Beschreibung: The general theme of the research my NASA colleague and I have planned is "Optical and probe diagnostics applied to reacting flows". We plan to explore three major threads during the fellowship period. The first interrogates the flame synthesis of carbon nanotubes using aerosol catalysts. Having demonstrated the viability of the technique for nanotube synthesis, we seek to understand the details of this reacting system which are important to its practical application. Laser light scattering will reveal changes in particle size at various heights above the burner. Analysis of the flame gas by mass spectroscopy will reveal the chemical composition of the mixture. Finally, absorption measurements will map the nanotube concentration within the flow. The second thread explores soot oxidation kinetics. Despite the impact of soot on engine performance, fire safety and pollution, models for its oxidation are inhibited by uncertainty in the values of the oxidation rate. We plan to employ both optical and microscopic measurements to refine this rate. Cavity ring-down absorption measurements of the carbonaceous aerosol can provide a measure of the mass concentration with time and, hence, an oxidation rate. Spectroscopic and direct probe measurements will provide the temperature of the system needed for subsequent modeling. These data will be benchmarked against changes in soot nanostructures as revealed by transmission electron microscopic images from directly sampled material.

Beschreibung: Active flow control of boundary-layer separation using glow-discharge plasma actuators is studied experimentally. Separation is induced on a flat plate installed in a closed-circuit wind tunnel by a shaped insert on the opposite wall. The flow conditions represent flow over the suction surface of a modem low-pressure-turbine airfoil. The Reynolds number, based on wetted plate length and nominal exit velocity, is varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low (0.2%) and high (2.5%) free-stream turbulence intensities are set using passive grids. A spanwise-oriented phased-plasma-array actuator, fabricated on a printed circuit board, is surface-flush-mounted upstream of the separation point and can provide forcing in a wide frequency range. Static surface pressure measurements and hot-wire anemometry of the base and controlled flows are performed and indicate that the glow-discharge plasma actuator is an effective device for separation control.

Beschreibung: 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.

Beschreibung: Most combustion processes have, in some way or another, a recirculating flow field. This recirculation stabilizes the reaction zone, or flame, but an unnecessarily large recirculation zone can result in high nitrogen oxide (NOx) values for combustion systems. The size of this recirculation zone is crucial to the performance of state-of-the-art, low-emissions hardware. If this is a large-scale combustion process, the flow field will probably be turbulent and, therefore, three-dimensional. This research dealt primarily with flow fields resulting from lean direct injection (LDI) concepts, as described in Research & Technology 2001. LDI is a concept that depends heavily on the design of the swirler. The LDI concept has the potential to reduce NOx values from 50 to 70 percent of current values, with good flame stability characteristics. It is cost effective and (hopefully) beneficial to do most of the design work for an LDI swirler using computer-aided design (CAD) and computer-aided engineering (CAE) tools. Computational fluid dynamics (CFD) codes are CAE tools that can calculate three-dimensional flows in complex geometries. However, CFD codes are only beginning to correctly calculate the flow fields for complex devices, and the related combustion models usually remove a large portion of the flow physics.

Beschreibung: Knudsen effusion mass spectrometry (KEMS) allows the simultaneous determination of the identity and pressure of vapor species in equilibrium with a condensed phase as a function of temperature. This information can be used to determine the thermodynamic properties of materials. The partial pressure of species j in the cell is related to the measured intensity of the ion k formed from j, I(sup +) less than SABjk, and the absolute temperature T, where S(sub jk) is the sensitivity factor.

Beschreibung: Rotating instability is a phenomenon that occurs in the tip flow region of axial compressor stages during stable operation. It can be observed in highly staggered rotors with significant tip clearance and is strongest at high-load operating points where the characteristic levels off. In this condition, the single-stage fan under investigation radiates an audible, whistling tone, and wall pressure spectra in the vicinity of the rotor disk exhibit nonrotational components. The graph shows the spectrum of static pressure at a point on the endwall near the leading edge. A hump appears at roughly half of the blade passing frequency (BPF) and is characteristic of rotating instability. A computational model was developed at the NASA Glenn Research Center to investigate the mechanism behind this phenomenon. A three-dimensional steady Navier-Stokes code that has been successfully tested for a wide range of turbomachinery flows was modified to execute a time-accurate simulation of the full annulus of the compressor. At the inlet of the computational domain, the total pressure, total temperature, and two velocity components are specified. Since no unsteady measurements of static pressure or other flow variables were available downstream of the rotor, circumferentially averaged static pressure was specified on the shroud at the outlet of the computational domain. A three-dimensional view of the vortex from the numerical model is shown. Particle traces released near the leading edge tip have rolled up to illustrate the tip clearance vortex. Flow near the trailing edge is pushed forward by the axially reversed flow. It then interacts with the tip clearance flow and the incoming flow and results in the rotating instability vortex, the core of which is illustrated by total pressure shading on planes located successively downstream. The rotating instability vortex is formed periodically midway between the blades and moves toward the pressure side of the passage. The unsteady behavior of this vortex structure is the main mechanism of the rotating instability is shown. The numerical model can be used to detect any possible occurrence of rotating instability when the tip clearance increases during engine service.

Beschreibung: High-capacity cooling options remain limited for many small-scale applications such as microelectronic components, miniature sensors, and microsystems. A microelectromechanical system (MEMS) using a Stirling thermodynamic cycle to provide cooling or heating directly to a thermally loaded surface is being developed at the NASA Glenn Research Center to meet this need. The device can be used strictly in the cooling mode or can be switched between cooling and heating modes in milliseconds for precise temperature control. Fabrication and assembly employ techniques routinely used in the semiconductor processing industry. Benefits of the MEMS cooler include scalability to fractions of a millimeter, modularity for increased capacity and staging to low temperatures, simple interfaces, limited failure modes, and minimal induced vibration. The MEMS cooler has potential applications across a broad range of industries such as the biomedical, computer, automotive, and aerospace industries. The basic capabilities it provides can be categorized into four key areas: 1) Extended environmental temperature range in harsh environments; 2) Lower operating temperatures for electronics and other components; 3) Precision spatial and temporal thermal control for temperature-sensitive devices; and 4) The enabling of microsystem devices that require active cooling and/or temperature control. The rapidly expanding capabilities of semiconductor processing in general, and microsystems packaging in particular, present a new opportunity to extend Stirling-cycle cooling to the MEMS domain. The comparatively high capacity and efficiency possible with a MEMS Stirling cooler provides a level of active cooling that is impossible at the microscale with current state-of-the-art techniques. The MEMS cooler technology builds on decades of research at Glenn on Stirling-cycle machines, and capitalizes on Glenn s emerging microsystems capabilities.

Beschreibung: Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.

Beschreibung: The generalization of a class of low-dissipative high order filter finite difference schemes for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been developed. The new scheme consists of a divergence free preserving high order spatial base scheme with a filter approach which can be divergence-free preserving depending on the type of filter operator being used, the method of applying the filter step, and the type of flow problem to be considered. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta * B) numerical error in the sense that no standard divergence cleaning is required. Performance evaluation of these variants, and the key role that the proper treatment of their corresponding numerical boundary conditions can play will be illustrated. Many levels of grid refinement and detailed comparison with several commonly used compressible MHD shock-capturing schemes will be sought. For certain MHD 2-D test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Beschreibung: An advanced thermal control system for the Burst Alert Telescope on the Swift satellite has been designed and an engineering test unit (ETU) has been built and tested in a thermal vacuum chamber. The ETU assembly consists of a propylene loop heat pipe, two constant conductance heat pipes, a variable conductance heat pipe (VCHP), which is used for rough temperature control of the system, and a radiator. The entire assembly was tested in a thermal vacuum chamber at NASA/GSFC in early 2002. Tests were performed with thermal mass to represent the instrument and with electrical resistance heaters providing the heat to be transferred. Start-up and heat transfer of over 300 W was demonstrated with both steady and variable condenser sink temperatures. Radiator sink temperatures ranged from a high of approximately 273 K, to a low of approximately 83 K, and the system was held at a constant operating temperature of 278 K throughout most of the testing. A novel LHP temperature control methodology using both temperature-controlled electrical resistance heaters and a small VCHP was demonstrated. This paper describes the system and the tests performed and includes a discussion of the test results.

Beschreibung: Modeling fuel slosh effects accurately and completely on spinning spacecraft has been a long standing concern within the aerospace community. Gyroscopic stiffness is obtained by spinning spacecraft as it is launched from one of the upper stages before it is placed in orbit. Unbalances in the spacecraft causes it to precess (wobble). The oscillatory motions are caused in the fuel due to this precession. This phenomenon is called 'fuel slosh' and the dynamic forces induced could adversely affect the stability and control of a spacecraft that spins about one of its minor moments of inertia axes. An equivalent mechanical model of fuel slosh is developed using springs and dampers that are connected to the rigid fuel tank. The stiffness and damping coefficients are the parameters that need to be identified in the fuel slosh mechanical model. This in turn is used in the spacecraft model to estimate the Nutation Time Constant (NTC) for the spacecraft. The experimental values needed for the identification of the model parameters are obtained from experiments conducted at Southwest Research Institute using the Spinning Slosh Test Rig (SSTR). The research focus is on developing models of different complexity and estimating the model parameters that will ultimately provide a more realistic Nutation Time Constant obtained through simulation.

Beschreibung: A computational and experimental study was conducted to investigate the effects of multiple injection ports in a two-dimensional, convergent-divergent nozzle, for fluidic thrust vectoring. The concept of multiple injection ports was conceived to enhance the thrust vectoring capability of a convergent-divergent nozzle over that of a single injection port without increasing the secondary mass flow rate requirements. The experimental study was conducted at static conditions in the Jet Exit Test Facility of the 16-Foot Transonic Tunnel Complex at NASA Langley Research Center. Internal nozzle performance was obtained at nozzle pressure ratios up to 10 with secondary nozzle pressure ratios up to 1 for five configurations. The computational study was conducted using the Reynolds Averaged Navier-Stokes computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. Internal nozzle performance was predicted for nozzle pressure ratios up to 10 with a secondary nozzle pressure ratio of 0.7 for two configurations. Results from the experimental study indicate a benefit to multiple injection ports in a convergent-divergent nozzle. In general, increasing the number of injection ports from one to two increased the pitch thrust vectoring capability without any thrust performance penalties at nozzle pressure ratios less than 4 with high secondary pressure ratios. Results from the computational study are in excellent agreement with experimental results and validates PAB3D as a tool for predicting internal nozzle performance of a two dimensional, convergent-divergent nozzle with multiple injection ports.

Beschreibung: It is well known that numerical warm season quantitative precipitation forecasts lack significant skill for numerous reasons. Some are related to the model--it may lack physical processes required to realistically simulate convection or the numerical algorithms and dynamics employed may not be adequate. Others are related to initialization-mesoscale features play an important role in convective initialization and atmospheric observation systems are incapable of properly depicting the three-dimensional stability structure at the mesoscale. The purpose of this study is to determine if a mesoscale model initialized with a diabatic initialization scheme can improve short-term (0 to 12h) warm season quantitative precipitation forecasts in the Southeastern United States. The Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory is used to diabatically initialize the Pennsylvania State University/National center for Atmospheric Research (PSUNCAR) Mesoscale Model version 5 (MM5). The SPORT Center runs LAPS operationally on an hourly cycle to produce analyses on a 15 km covering the eastern 2/3 of the United States. The 20 km National Centers for Environmental Prediction (NCEP) Rapid Update Cycle analyses are used for the background fields. Standard observational data are acquired from MADIS with GOES/CRAFT Nexrad data acquired from in-house feeds. The MM5 is configured on a 140 x 140 12 km grid centered on Huntsville Alabama. Preliminary results indicate that MM5 runs initialized with LAPS produce improved 6 and 12h QPF threat scores compared with those initialized with the NCEP RUC.

Beschreibung: Four Sunpower M87N Stirling-cycle cryocoolers will be used to extend the lifetime of the Alpha Magnetic Spectrometer-02 (AMS-02) experiment. The cryocoolers will be mounted to the AMS-02 vacuum case using a structure that will thermally and mechanically decouple the cryocooler from the vacuum case while providing compliance to allow force attenuation using a passive balancer system. The cryocooler drive is implemented using a 60Hz pulse duration modulated square wave. Details of the testing program, mounting assembly and drive scheme will be presented. AMS-02 is a state-of-the-art particle physics detector containing a large superfluid helium-cooled superconducting magnet. Highly sensitive detector plates inside the magnet measure a particle s speed, momentum, charge, and path. The AMS-02 experiment, which will be flown as an attached payload on the International Space Station, will study the properties and origin of cosmic particles and nuclei including antimatter and dark matter. Two engineering model cryocoolers have been under test at NASA Goddard since November 2001. Qualification testing of the engineering model cryocooler bracket assembly is near completion. Delivery of the flight cryocoolers to Goddard is scheduled for September 2003.

Beschreibung: It is well established that residual flows exist in contained liquid metal processes. In 1-g processing, buoyancy forces often drive these flows and their magnitudes can be substantial. It is also known that residual flows can exist during microgravity processing, and although greatly reduced in magnitude, they can influence the properties of the processed materials. Unfortunately, there are very few techniques to visualize flows in opaque, high temperature liquid metals, and those available are not easily adapted to flight investigation. In this study, a novel technique is developed that uses liquid tin as the model fluid and solid-state electrochemical cells constructed from Yttria-Stabilized Zirconia (YSZ) to establish and measure dissolved oxygen boundary conditions. The melt serves as a common electrode for each of the electrochemical cells in this design, while independent reference electrodes are maintained at the outside surfaces of the electrolyte. By constructing isolated electrochemical cells at various locations along the container walls, oxygen is introduced or extracted by imposing a known electrical potential or passing a given current between the melt and the reference electrode. This programmed titration then establishes a known oxygen concentration boundary condition at the selected electrolyte-melt interface. Using the other cells, the concentration of oxygen at the electrolyte-melt interface is also monitored by measuring the open-circuit potentials developed between the melt and reference electrodes. Thus the electrochemical cells serve to both establish boundary conditions for the passive tracer and sense its path. Rayleigh-Benard convection was used to validate the electrochemical approach to flow visualization. Thus, a numerical characterization of the second critical Rayleigh numbers in liquid tin was conducted for a variety of Cartesian aspect ratios. The extremely low Prandtl number of tin represents the lowest value studied numerically. Additionally, flow field oscillations are visualized and the effect of tilt on convecting systems is quantified. Experimental studies of the effect of convection in liquid tin are presented. Three geometries are studied: (1) double electrochemical cell with vertical concentration gradients; (2) double cell with horizontal concentration gradients; and (3) multiple cells with vertical temperature gradients. The first critical Rayleigh number transition is detected with geometry (1) and it is concluded that current measurements are not as affected by convection as EMF measurements. The system is compared with numerical simulations in geometry (2), and oscillating convection is detected with geometry (3).

Beschreibung: A magnetohydrodynamic model that examines the effect of rotating an electrically conducting cylinder with a uniform external magnetic field applied orthogonal to its axis is presented. Noting a simple geometry, it can be classified as a fundamental dynamo problem. For the case of an infinitely long cylinder, an analytical solution is obtained and analyzed in detail. A semi-analytical model was developed that considers a finite cylinder. Experimental data from a spinning brass wheel in the presence of Earth's magnetic field were compared to the proposed theory and found to fit well.

Beschreibung: A theoretical and experimental study is presented on the stability of solutal convection of a magnetized fluid in the presence of a magnetic field. The total force on the fluid is derived and equilibrium positions are computed establishing the field necessary to counter fluid buoyancy. The requirements for stability are developed and compared with experiments with a paramagnetic fluid. The experiments are in good agreement not only with the theoretical predictions for equilibrium but also verify the stability theory which predicts both horizontal and vertical stability. Analogous to results for levitation, the theory indicates that solutal convection in paramagnetic fluids cannot be completely stabilized while that in diamagnetic liquid are possible.

Beschreibung: Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Beschreibung: Progress can be reported in two areas related to characterizing the properties of cusp diamagnetic cavities. Laboratory terrella experiments have been conducted for the purpose of using neutral gas excitation as a tracer of trapped electron populations in the presence of two dipoles that are used to develop a magnetic cusp topology. Figure 1 and 2 show top and side views of two configurations. Dipole trapped electron populations appear as the two luminous annular rings. Other populations are the most intense regions are shown. Interspersed between these regions are narrow regions that represent the topological cusps in these configurations. That they contain luminous gas is evidence for cusp trapping similar to what we believe exists in the terrestrial magnetosphere. The asymmetry of these cusp regions as seen in Figure 1 is the result of a relative tilt between the two dipoles suggestive of what would be expected in space. It is in these regions that particle observations were sought, so as to validate the realization of proposed and laboratory achieved trapping in a diamagnetic cusp. Figure 3 shows particle trajectories in a modeled cusp magnetic topology for three particle energies. Blue, green, and red traces correspond to increasing energies. Due to factors discussed outside of this final report, a thorough exploration of relevant satellite observations have not been achieved.

Beschreibung: The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Beschreibung: Development of hydrogen in sealed silica glass ampoules during annealing at elevated temperatures was investigated. The dependence of hydrogen pressure in the ampoules as a function of time, for different temperatures and ampoule parameters was measured. The process was modeled assuming chemical solution of hydrogen according to the reaction: silica + H2 = H- Si= + H-O-Si=. The equilibrium constant of the reaction was determined by fitting the theoretical curves to the experimental data. The Gibbs function for this reaction was estimated at deltaG = -25.8 + 54T.

Beschreibung: Volumetric forces resulting from local density variations and gravitational acceleration cause buoyancy induced convective motion in melts and solutions. Solutal buoyancy is a result of concentration differences in an otherwise isothermal fluid. If the fluid also exhibits variations in magnetic susceptibility with concentration then convection control by external magnetic fields can be hypothesized. Magnetic control of thermal buoyancy induced convection in ferrofluids (dispersions of ferromagnetic particles in a carrier fluid) and paramagnetic fluids have been demonstrated. Here we show the nature of magnetic control of solutal buoyancy driven convection of a paramagnetic fluid, an aqueous solution of Manganese Chloride hydrate. We predict the critical magnetic field required for balancing gravitational solutal buoyancy driven convection and validate it through a simple experiment. We demonstrate that gravity driven flow can be completely reversed by a magnetic field but the exact cancellation of the flow is not possible. This is because the phenomenon is unstable. The technique can be applied to crystal growth processes in order to reduce convection and to heat exchanger devices for enhancing convection. The method can also be applied to impose a desired g-level in reduced gravity applications.

Beschreibung: Direct and large-eddy simulations of the interaction between the wake of a circular cylinder and a flat-plate boundary layer are conducted. Two Reynolds numbers are examined. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re(sub theta) achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established.

Beschreibung: We analyze the velocity and temperature fields at steady state due to thermocapillary convection around a gas bubble that is stationary in a liquid. A linear temperature field is imposed in the undisturbed liquid. Our interest is in investigating the effect of convective transport of momentum and energy on the velocity and temperature fields. We assume the pertinent physical properties to be constant, and that buoyant convection is negligible. Suitably defined Reynolds and Marangoni numbers are assumed to be small compared with unity. When both the Reynolds and Marangoni numbers are set equal to zero, a solution can be found. In this solution, far from the bubble, the velocity field decays as the inverse of the distance from the bubble, and the disturbance temperature field decays as the inverse of the square of this distance. We now attempt to obtain a solution when the Reynolds number is zero, but the Marangoni number is small, but non-zero, by a perturbation expansion in the Marangoni number. When the temperature field is expanded in a regular perturbation series in the Marangoni number, we show that the problem for the first correction field is ill-posed. The governing equation for this perturbation field contains an inhomogeneity, and the corresponding particular solution neither decays far from the bubble, nor can be canceled by a homogeneous solution. Additional information is included in the original extended abstract.

Beschreibung: A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.

Beschreibung: This report discusses the: 1. Development coupled electro/fluid/structural lumped-element model (LEM) of a prototypical flow-control actuator. 2. Validation the coupled electro/fluid/structural dynamics lumped-element models. 3. Development simple, yet effective, design tools for actuators. 4. Development structural dynamic models that accurately characterize the dynamic response of piezoelectric flap actuators using the Finite Element Method (FEW as well as analytical methods. 5. Perform a parametric study of a piezo-composite flap actuator. 6.Develop an optimization scheme for maximizing the actuator performance.

Beschreibung: The objective of this research is to improve and implement the filtered mass density function (FDF) methodology for large eddy simulation (LES) of high-speed reacting turbulent flows. We have just completed Year 1 of this research. This is the Final Report on our activities during the period: January 1, 2003 to December 31, 2003. 2002. In the efforts during the past year, LES is conducted of the Sandia Flame D, which is a turbulent piloted nonpremixed methane jet flame. The subgrid scale (SGS) closure is based on the scalar filtered mass density function (SFMDF) methodology. The SFMDF is basically the mass weighted probability density function (PDF) of the SGS scalar quantities. For this flame (which exhibits little local extinction), a simple flamelet model is used to relate the instantaneous composition to the mixture fraction. The modelled SFMDF transport equation is solved by a hybrid finite-difference/Monte Carlo scheme.

Beschreibung: We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.

Beschreibung: This paper presents the results obtained from a recent ethane heat pipe program. Three identical ethane heat pipes were tested individually, and then two selected heat pipes were tested collectively in their system configuration. Heat transport, thermal conductance, and non-condensable gas tests were performed on each heat pipe. To gain insight into the reflux operation as seen at spacecraft level ground testing, the test fixture was oriented in a vertical configuration. The system level test included a computer-controlled heater designed to emulate the heat load generated at the thermoelectric cooler interface. The system performance was successfully characterized for a wide range of environmental conditions while staying within the operating limits.

Beschreibung: A silicon-based microhotplate tin oxide (SnO2) gas sensor integrated into a polymer-based microfluidic system for monitoring of contaminants in water systems is presented. This device is designed to sample a water source, control the sample vapor pressure within a microchannel using integrated resistive heaters, and direct the vapor past the integrated gas sensor for analysis. The sensor platform takes advantage of novel technology allowing direct integration of discrete silicon chips into a larger polymer microfluidic substrate, including seamless fluidic and electrical interconnects between the substrate and silicon chip.

Beschreibung: Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.

Beschreibung: The space-time conservation-element and solution-element method is employed to numerically study the near-field screech-tone noise of a typical underexpanded circular jet issuing from a sonic nozzle. Both axisymmetric and fully three-dimensional computations are carried out. The self-sustained feedback loop is properly simulated. The computed shock-cell structure, acoustic wave length, screech-tone frequency, and sound-pressure levels are in good agreement with existing experimental results.

Beschreibung: The flow fields of synthetic jets in a cross-flow from orifices of different geometry are investigated. The geometries include a straight, a tapered, a pitched and a cluster of nine orifices, all having the same cross-sectional area through which the perturbation is discharged into the cross-flow. The strength of the jet from the tapered orifice in comparison to that from the straight one is found to be only slightly enhanced. The flow field from the cluster of orifices, when viewed a few equivalent diameters downstream, is similar to that from the single orifice. However, the penetration is somewhat lower in the former case due to the increased mixing of the distributed jets with the cross-flow. The penetration for the pitched configuration is the lowest, as expected. The jet trajectories for the straight and pitched orifices are well represented by correlation equations available for steady jets-in-cross-flow. Distributions of streamwise velocity, vorticity as well as turbulence intensity are documented for various cases. In addition, distributions of phase-averaged velocity and vorticity for the cylindrical and the clustered orifices are presented providing an insight into the flow dynamics.

Beschreibung: Contents include the following: 1. High frequency, low amplitude temperature oscillations: LHP operation - governing equations; interactions among LHP components; factors affecting low amplitude temperature oscillations. 2. Test results. 3. Conclusions.

Beschreibung: This paper presents viewgraphs on the low frequency high amplitude temperature oscillations observed in loop heat pipe operations. The topics include: 1) Proposed Theory; 2) Test Loop and Test Results; and 3) Effects of Various Parameters. The author also presents a short summary on the conditiions that must be met in order to sustain a low frequency high amplitude temperature oscillation.

Beschreibung: This work examined the wake aerodynamics of a single helicopter rotor blade with several tip shapes operating on a hover test stand. Velocity field measurements were conducted using three-component laser Doppler velocimetry (LDV). The objective of these measurements was to document the vortex velocity profiles and then extract the core properties, such as the core radius, peak swirl velocity, and axial velocity. The measured test cases covered a wide range of wake-ages and several tip shapes, including rectangular, tapered, swept, and a subwing tip. One of the primary differences shown by the change in tip shape was the wake geometry. The effect of blade taper reduced the initial peak swirl velocity by a significant fraction. It appears that this is accomplished by decreasing the vortex strength for a given blade loading. The subwing measurements showed that the interaction and merging of the subwing and primary vortices created a less coherent vortical structure. A source of vortex core instability is shown to be the ratio of the peak swirl velocity to the axial velocity deficit. The results show that if there is a turbulence producing region of the vortex structure, it will be outside of the core boundary. The LDV measurements were supported by laser light-sheet flow visualization. The results provide several benchmark test cases for future validation of theoretical vortex models, numerical free-wake models, and computational fluid dynamics results.

Beschreibung: A finite volume based network analysis procedure has been applied to model unsteady flow without and with heat transfer. Liquid has been modeled as compressible fluid where the compressibility factor is computed from the equation of state for a real fluid. The modeling approach recognizes that the pressure oscillation is linked with the variation of the compressibility factor; therefore, the speed of sound does not explicitly appear in the governing equations. The numerical results of chilldown process also suggest that the flow and heat transfer are strongly coupled. This is evident by observing that the mass flow rate during 90-second chilldown process increases by factor of ten.

Beschreibung: Numerical simulations have been completed for a variety of designs for a 90 deg elbow duct. The objective is to identify a design that minimizes the dynamic load entering a LOX turbopump located at the elbow exit. Designs simulated to date indicate that simpler duct geometries result in lower losses. Benchmark simulations have verified that the compressible flow codes used in this study are applicable to these incompressible flow simulations.

Beschreibung: Reynolds-averaged Navier-Stokes calculations are used to investigate porous side-edge treatment as a passive means for flap noise reduction. Steady-state simulations are used to infer effects of the treatment on acoustically relevant features of the mean flow near the flap side edge. Application of the porous treatment over a miniscule fraction of the wetted flap area (scaling with the flap thickness) results in significantly weaker side-edge vortex structures via modification of the vortex initiation and roll-up processes. At high flap deflections, the region of axial flow reversal associated with the breakdown of the side-edge vortex is also eliminated, indicating an absence of vortex bursting in the presence of the treatment. Potential ramifications of the mean-flow modifications for flap-noise reduction are examined in the light of lessons learned from recent studies on flap noise. Computations confirm that any noise reduction benefit via the porous treatment would be achieved without compromising the aerodynamic effectiveness of the flap. Results of the parameter study contribute additional insight into the measured data from the 7x10 wind tunnel at NASA Ames and provide preliminary guidance for specifying optimal treatment characteristics in terms of treatment location, spatial extent, and flow resistance of the porous skin.

Beschreibung: As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Beschreibung: This paper presents the Turbine Air-Flow Test (TAFT) rig computational fluid dynamics (CFD) results for test matrix. The topics include: 1) TAFT Background; 2) Design Point CFD; 3) TAFT Test Plan and Test Matrix; and 4) CFD of Test Points. This paper is in viewgraph form.

Beschreibung: This report consists of two published papers, 'Computation of Magnetohydrodynamic Flows Using an Iterative PNS Algorithm' and 'Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm'.

Beschreibung: A computational and experimental study was conducted to investigate the effects of multiple injection ports in a two-dimensional, convergent-divergent nozzle, for fluidic thrust vectoring. The concept of multiple injection ports was conceived to enhance the thrust vectoring capability of a convergent-divergent nozzle over that of a single injection port without increasing the secondary mass flow rate requirements. The experimental study was conducted at static conditions in the Jet Exit Test Facility of the 16-Foot Transonic Tunnel Complex at NASA Langley Research Center. Internal nozzle performance was obtained at nozzle pressure ratios up to 10 with secondary nozzle pressure ratios up to 1 for five configurations. The computational study was conducted using the Reynolds Averaged Navier-Stokes computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. Internal nozzle performance was predicted for nozzle pressure ratios up to 10 with a secondary nozzle pressure ratio of 0.7 for two configurations. Results from the experimental study indicate a benefit to multiple injection ports in a convergent-divergent nozzle. In general, increasing the number of injection ports from one to two increased the pitch thrust vectoring capability without any thrust performance penalties at nozzle pressure ratios less than 4 with high secondary pressure ratios. Results from the computational study are in excellent agreement with experimental results and validates PAB3D as a tool for predicting internal nozzle performance of a two dimensional, convergent-divergent nozzle with multiple injection ports.

Beschreibung: During the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of CFD discipline. Even though incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient, CFD tools become indispensable in fluid engineering for incompressible and low speed flow. This paper is intended to review some of the successes made possible by advances in computational technologies during the same period, and discuss some of the current challenges.

Beschreibung: A generalized multi-phase formulation for cavitation in fluids operating at temperatures elevated relative to their critical temperatures is presented. The thermal effects and the accompanying property variations due to phase change are modeled rigorously. Thermal equilibrium is assumed and fluid thermodynamic properties are specified along the saturation line using the NIST-12 databank. Fundamental changes in the physical characteristics of the cavity when thermal effects become pronounced are identified; the cavity becomes more porous, the interface less distinct, and has increased entrainment when temperature variations are present. Quantitative estimates of temperature and pressure depressions in both liquid nitrogen and liquid hydrogen were computed and compared with experimental data of Hord for hydrofoils. Excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable. Liquid nitrogen cavities were consistently found to be in thermal equilibrium while liquid hydrogen cavities exhibited small, but distinct, non-equilibrium effects.

Beschreibung: A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment using 13 computers located at 4 different geographical sites. Process automation and web-based access to the database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The database consists of forces, moments, and solution files obtained by varying the Mach number, angle of attack, and sideslip angle. The forces and moments compare well with experimental data. Stability derivatives are also computed using a monotone cubic spline procedure. Flow visualization and three-dimensional surface plots are used to interpret and characterize the nature of computed flow fields.

Beschreibung: In this study 3-D Computational Fluid Dynamic (CFD) runs have been made for the Space Shuttle RSRM using 2 different grids and 4 different turbulent models, which were the Standard KE, the RNG KE, the Realizable KE, and the Reynolds stress model. The RSRM forward segment consists of 11 fins. By taking advantage of the forward fin symmetry only half of one fin along the axis had to be used in making the grid. This meant that the 3-D model consisted of a pie slice that encompassed 1/22nd of the motor circumference and went along the axis of the entire motor. The 3-D flow patterns in the forward fin region are of particular interest. Close inspection of these flow patterns indicate that 2 counter-rotating axial vortices emerge from each submerged solid propellant fin. Thus, the 3-D CFD analysis allows insight into complicated internal motor flow patterns that are not available from the simpler 2-D axi-symmetric studies. In addition, a comparison is made between the 3-D bore pressure drop and the 2-D axi-symmetric pressure drop.

Beschreibung: A four-parameter, three-level, central composite experiment design has been used to optimize the configuration of an annular scramjet injector geometry using computational fluid dynamics. The computational fluid dynamic solutions played the role of computer experiments, and response surface methodology was used to capture the simulation results for mixing efficiency and total pressure recovery within the scramjet flowpath. An optimization procedure, based upon the response surface results of mixing efficiency, was used to compare the optimal design configuration against the target efficiency value of 92.5%. The results of three different optimization procedures are presented and all point to the need to look outside the current design space for different injector geometries that can meet or exceed the stated mixing efficiency target.

Beschreibung: A concerted effort is underway at NASA Langley Research Center to create a benchmark for Computational Fluid Dynamic (CFD) codes. both unstructured and structured, against a data set for the hump model with actuation. The hump model was tested in the NASA Langley 0.3-m Transonic Cryogenic Tunnel. The CFD codes used for the analyses are the FUN2D (Full Unstructured Navier-Stokes 2-Dimensional) code, the structured TLNS3D (Thin-Layer Navier-Stokes 3-Dimensional) code, and the structured CFL3D code, all developed at NASA Langley. The current investigation uses the time-accurate Reynolds-Averaged Navier-Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. Two-dimensional computational results verified that steady blowing and suction and oscillatory suction/blowing can be used to significantly reduce the separated flow region on the model. Discrepancies do exist between the CFD results and experimental data in the region downstream of the slot with the largest differences in the oscillatory cases. Overall, the structured CFD codes exhibited similar behavior with each other for a wide range of control conditions, with the unstructured FUN2D code showing moderately different results in the separated flow region for the suction and oscillatory cases.

Beschreibung: The focus of this research is on supercritical C7H16/N2 and O2/H2 mixing layers undergoing transitions to turbulence. The C7H16/N2 system serves as a simplified model of hydrocarbon/air systems in gas-turbine and diesel engines; the O2/H2 system is representative of liquid rocket engines. One goal of this research is to identify ways of controlling area production to increase disintegration of fluids and enhance combustion in such engines. As used in this research, "area production" signifies the fractional rate of change of surface area oriented perpendicular to the mass-fraction gradient of a mixing layer. In the study, a database of transitional states obtained from direct numerical simulations of the aforementioned mixing layers was analyzed to investigate global layer characteristics, phenomena in regions of high density-gradient magnitude (HDGM), irreversible entropy production and its relationship to the HDGM regions, and mechanisms leading to area production.

Beschreibung: This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Beschreibung: Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Beschreibung: Immiscible fluid layers are found in a host of applications ranging from materials processing, for example the use of encapsulants in float zone crystal growth technique and a buffer layer in industrial Czochralski growth of crystals to prevent Marangoni convection, to heat transfer phenomena in day-to-day processes like the presence of air pockets in heat exchangers. In the microgravity and space processing realm, the exploration of other planets requires the development of enabling technologies in several fronts. The reduction in the gravity level poses unique challenges for fluid handling and heat transfer applications. The present work investigates the efficacy of controlling thermal convective flow using magnetic fluids and magnetic fields. The setup is a two-layer immiscible liquids system with one of the fluids being a diluted ferrofluid (super paramagnetic nano particles dispersed in carrier fluid). Using an external magnetic field one can essentially dial in a volumetric force - gravity level, on the magnetic fluid and thereby affect the system thermo-fluid behavior. The paper will describe the experimental and numerical modeling approach to the problem and discuss results obtained to date.

Beschreibung: Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Beschreibung: Experiments to suppress thermocapillary oscillations using high-frequency vibrations were performed on float-zones. Such a float-zone is formed by melting one end of a vertically held sodium nitrate-barium nitrate crystal rod in contact with a hot surface at the top. In the experiments, when thermocapillary oscillation occurred, the bottom end of the rod was vibrated at a high frequency to generate fine ripples on the melt surface, driving a streaming flow in the opposite direction to that of the thermocapillary convection. It was observed that by generating a sufficiently strong streaming flow the thermocapillary flow can be offset enough such that the associated thermocapillarity oscillations can be quenched.

Beschreibung: Imagine a sequence of equidistant steps parallel to one another on a crystal face growing from solution. Steps may be one or part of a lattice spacing high. With such ideal morphology, the most perfect material is expected to be produced. Sometimes this happens. Much more often, however, steps form bunches, acquire wavy shapes and interlace. With time, these spontaneous perturbations may be increased in amplitude, i.e. the stepped surface becomes morphologically unstable. Interferometric studies of Y.G. Kuznetsov, L.N. Rashkovich, P.G. Vekilov and N.A. Booth and the author demonstrated that the step bunching depends at least on direction and rate of solution flow over the interface, supersaturation and presence of impurities. In particular, if solution and step flow are parallel to one another, the steps convene into bunches. If these flows are antiparallel, the bunches disappear. As the step train propagates, the overall bunch height often increases infinitely. However, we have found recently that if the flow rate is large enough, (approximately 1 meter per second), the bunch height may be limited. All these phenomena come from interaction between steps. The interaction mechanisms are still not fully understood. I plan to overview major approaches to the still unresolved problem on how these dissipative structures on growing crystal face appear and evolve.

Beschreibung: Turbulent fluxes of momentum, moisture, and heat at the air-sea interface are essential for climate studies. Version 2 Goddard Satellite-based Surface Turbulent Fluxes (GSSTF2) has been derived from the Special Sensor Microwave/Imager (SSM/I) radiance measurements. This dataset, covering the period July 1987-December 2000 over global oceans, has a spatial resolution of 1 deg x 1 deg lat-long and a temporal resolution of 1 day. Turbulent fluxes are derived from the SSM/I surface winds and surface air humidity, as well as the 2-m air and sea surface temperatures (SST) of the NCEP/NCAR reanalysis, using a bulk aerodynamic algorithm based on the surface layer similarity theory. The GSSTF2 bulk flux model, and retrieved daily wind stress, latent heat flux, wind speed, and surface air humidity validate well with ship observations of ten field experiments over the tropical and midlatitude oceans during 1991-99. The global distributions of 1988-2000 annual- and seasonal-mean turbulent fluxes show reasonable patterns related to the atmospheric general circulation and seasonal variations. Latent heat fluxes and related input parameters over global oceans during 1992-93 have been compared among GSSTF1 (version 1), GSSTF2, HOAPS (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data), NCEP/NCAR reanalysis, and one based on COADS (Comprehensive Ocean-Atmosphere Data Set). Our analyses suggest that the GSSTF2 latent heat flux, surface air humidity, surface wind, and SST are quite realistic compared to the other four flux datasets examined. However, significant differences are found among these five flux datasets. The GSSTF2, available at http://daac.gsfc.nasa.gov/CAMPAIGN_DOCS/hydrology/hd_gsstf2.O.html, is useful for climate studies.

Beschreibung: Containerless processing is an important tool for materials research. The freedom from a crucible allows processing of liquid materials in a metastable undercooled state, as well as allowing processing of high temperature and highly reactive melts. Electrostatic levitation (ESL) is a containerless method which provides a number of unique advantages, including the ability to process non-conducting materials, the ability to operate in ultra-high vacuum or at moderate gas pressure (approx. = 5 atm), and the decoupling of positioning force from sample heating. ESL also has the potential to reduce internal flow velocities below those possible with electromagnetic, acoustic, or aero-acoustic techniques. In electrostatic levitation, the acceleration of gravity (or residual acceleration in reduced gravity) is opposed by the action of an applied electric field on a charged sample. Microgravity allows electrostatic levitation to work even more effectively. The ESL facility at NASA s Marshall Space Flight Center is in use for materials research and thermophysical property measurement by a number of different internal and external investigators. Results from the recent CDDF studies on the high energy X-ray beamline at the Advanced Photon Source of Argonne National Laboratory will be presented. The Microgravity Research Program supports the facility.

Beschreibung: We investigate the convective dynamics of evaporating pools of volatile liquids using an ultra-sensitive thermal imaging camera. During evaporation, there are significant convective flows inside the liquid due to Marangoni forces. We find that Marangoni convection during evaporation can dramatically affect the evaporation rates of volatile liquids. A simple heat balance model connects the convective velocities and temperature gradients to the evaporation rates.

Beschreibung: Low-dissipative high order filter finite difference methods for shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (del (raided dot) B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Beschreibung: Low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta . B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Beschreibung: The capabilities of the Dense Media Radiative Transfer model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP) to reproduce measured radiometric data were tested using the University of Tokyo Ground Based Microwave Radiometer (GBMR-7) during the third Intensive Observation Period (IOP3) of the NASA Cold-land Processes Field Experiment (CLPX). The data were collected at the Local-Scale Observation Site (LSOS), a 0.8-ha study site consisting of two open meadows separated by trees. Intensive measurements were also made of snow depth and temperature, density, and grain size profiles. A DMRT model is needed to describe radiative transfer in a medium such as snow because the assumption of independent scattering used in classical radiative transfer theory (CRT) is not valid. Validation of the DMRT approach requires a relationship between measured snow grain size and the DMRT approximation of snow grain radius as spherical particles with a mean radius of the log-normal particle-size distribution. This relationship is very important for a better understanding of snow modelling and for practical applications. DMRT simulations were compared with observations of microwave brightnesses at 18.7, 36.5 and 89-GHz (V and H polarizations) collected on February-1 9-25, 2003. Observation angles ranged from 30\deg to 70\deg. Model inputs included measured snow parameters except mean grain size. The average snow temperature, fractional volume and depth were held constant, together with the ice and soil permittivities. The minimum and maximum measured mean grain sizes were used to test the capabilities of the DMRT to reproduce the brightnesses as upper and lower limits. The sensitivity to the largest and smallest measured grain size in the three classes of minimum, medium and maximum observed grain sizes was also investigated. DMRT particle sizes yielding a best-fit to the experimental data for each date were computed. Results show that the measured brightnesses fall within the range of simulated brightnesses using the smallest and largest measured grain size values. The DMRT best-fit radii are comparable to the average radii for the medium observed grain sizes.

Beschreibung: A new, interferometer-based system for measuring thermal expansion to an absolute accuracy of 20 ppb or better at cryogenic temperatures has been developed. Data from NIST Copper SRM 736 measured from room temperature to 15 K will be presented along with data from many other materials including beryllium, ULE, Zerodur, and composite materials. Particular attention will be given to a study by the Space Optics Manufacturing Technology Center (SOMTC) investigating the variability of ULE and beryllium materials used in the AMSD program Approximately 20 samples of each material, tested from room temperature to below 30 K are compared as a function of billet location.

Beschreibung: Experiments on solutal convection in a paramagnetic fluid were conducted in a strong magnetic field gradient using a dilute solution of Manganese Chloride. The observed flows indicate that the magnetic field can completely counter the settling effects of gravity locally and are consistent with the theoretical predictions presented.

Beschreibung: One of the power systems under consideration for nuclear electric propulsion or as a planetary surface power source is a heatpipe-cooled reactor coupled to a Brayton cycle. In this system, power is transferred from the heatpipes to the Brayton gas via a heat exchanger attached to the heatpipes. This paper discusses the fluid, thermal and structural analyses that were performed in support of the design of the heat exchanger to be tested in the SAFE-100 experimental program at Marshall Space Flight Center. A companion paper, "Mechanical Design and Fabrication of a SAFE-100 Heat Exchanger for use in NASA s Advanced Propulsion Thermal-hydraulic Simulator", presents the fabrication issues and prototyping studies that, together with these analyses, led to the development of this heat exchanger. An important consideration throughout the design development of the heat exchanger was its capability to be utilized for higher power and temperature applications. This paper also discusses this aspect of the design and presents designs for specific applications that are under consideration.

Beschreibung: Re-circulating flow in molten metal columns can be conveniently induced by the axisymmetric traveling magnetic field. A number of applications can benefit from this technique, such as mixing under microgravity environment, or.crysta1 growth from metallic melts. For small magnetic field excitations, the flow is laminar and stationary. As the imposed field increases, a more complex flow will set up in the cylindrical column. Conditions for stable laminar flow are of importance for practical applications. In this work, a linear stability analysis is performed in order to determine the onset of the bifurcation in the system. Here the analysis is restricted to the axisymmetric modes and the low-frequency regime.

Beschreibung: Velocity fluctuations and particle concentrations are studied in a liquid fluidized bed to investigate the origin of steady state sedimentation. Both the velocity fluctuations and the particle concentrations are found to strongly depend on height. A flux balance model shows why the bed is stable: velocity fluctuations drive a downward particle flux that is compensated by an upward particle flux stemming from an excess flow velocity due to the stratification in concentration. Our results show that in steady state the magnitudes of the fluctuations are related to the degree of stratification.

Beschreibung: The study examined the effect of natural convection in applying the modified laser flash method to measure thermal properties of semiconductor melts. Common laser flash method uses a laser pulse to heat one side of a thin circular sample and measures the temperature response of the other side. Thermal diffusivity can be calculations based on a heat conduction analysis. For semiconductor melt, the sample is contained in a specially designed quartz cell with optical windows on both sides. When laser heats the vertical melt surface, the resulting natural convection can introduce errors in calculation based on heat conduction model alone. The effect of natural convection was studied by CFD simulations with experimental verification by temperature measurement. The CFD results indicated that natural convection would decrease the time needed for the rear side to reach its peak temperature, and also decrease the peak temperature slightly in our experimental configuration. Using the experimental data, the calculation using only heat conduction model resulted in a thermal diffusivity value is about 7.7% lower than that from the model with natural convection. Specific heat capacity was about the same, and the difference is within 1.6%, regardless of heat transfer models.

Beschreibung: The Lorentz force induced by a traveling magnetic field (TMF) in a cylindrical container has been calculated. The force can be used to control flow in dectrically conducting melts and the direction of the magnetic field and resulting flow can be reversed. A TMF can be used to partially cancel flow driven by buoyancy. The penetration of the field into the cylinder decreases as the frequency increases, and there exists an optimal value of frequency for which the resulting force is a maximum. Expressions for the Lorentz force in the limiting cases of low frequency and infinite cylinder are also given and compared to the numerical calculations.

Beschreibung: The ancillary components for Unitized Regenerative Fuel Cell (URFC) Energy Storage System are being developed at the NASA Glenn Research Center. This URFC system is unique in that it uses the surface area of the hydrogen and oxygen storage tanks as radiating heat surfaces for overall thermal control of the system. The waste heat generated by the URFC stack during charging and discharging is transferred from the cell stack to the surface of each tank by loop heat pipes. The heat pipes are coiled around each tank and covered with a thin layer of thermally conductive layer of carbon composite. The thin layer of carbon composite acts as a fin structure that spreads the heat away from the heat pipe and across the entire tank surface. Two different sized commercial grade composite tanks were constructed with integral heat pipes and tested in a thermal vacuum chamber to examine the feasibility of using the storage tanks as system radiators. The storage radiators were subjected to different steady-state heat loads and varying heat load profiles. The surface emissivity and specific heat capacity of each tank were calculated. The results were incorporated into a model that simulates the performance of similar radiators using lightweight, space rated carbon composite tanks.

Beschreibung: The focus of this project was to study the physical processes that govern tachocline dynamics and structure. Specific features explored included stratification, shear, waves, and toroidal and poloidal background fields. In order to address recent theoretical work on anisotropic mixing and dynamics in the tachocline, we were particularly interested in such anisotropic mixing for the specific tachocline processes studied. Transition to turbulence often shapes the largest-scale features that appear spontaneously in a flow during the development of turbulence. The resulting large-scale straining field can control the subsequent dynamics; therefore, anticipation of the large-scale straining field that results for individual realizations of the transition to turbulence can be important for subsequent dynamics, flow morphology, and transport characteristics. As a result, we paid particular attention to the development of turbulence in the stratified and sheared environment of the tachocline. This is complicated by the fact that the linearly stability of sheared MHD flows is non-self-adjoint, implying that normal asymptotic linear stability theory may not be relevant.

Beschreibung: Pulsed combustion appears to have the potential to provide for rapid fuel/air mixing, compact and economical combustors, and reduced exhaust emissions. The objective of this experiment (PuFF, for Pulsed-Fully Flames) is to increase the fundamental understanding of the fuel/air mixing and combustion behavior of pulsed, turbulent diffusion flames by conducting experiments in microgravity. In this research the fuel jet is fully-modulated (i.e., completely shut off between pulses) by an externally controlled valve system. This gives rise to drastic modification of the combustion and flow characteristics of flames, leading to enhanced fuel/air mixing compared to acoustically excited or partially-modulated jets. Normal-gravity experiments suggest that the fully-modulated technique also has the potential for producing turbulent jet flames significantly more compact than steady flames with no increase in exhaust emissions. The technique also simplifies the combustion process by avoiding the acoustic forcing generally present in pulsed combustors. Fundamental issues addressed in this experiment include the impact of buoyancy on the structure and flame length, temperatures, radiation, and emissions of fully-modulated flames.

Beschreibung: Schiller, Ross, and Sirignano (1996) studied ignition and flame spread above liquid fuels initially below the flashpoint temperature by using a two-dimensional computational fluid dynamics code that solves the coupled equations of both the gas and the liquid phases. Pulsating flame spread was attributed to the establishment of a gas-phase recirculation cell that forms just ahead of the flame leading edge because of the opposing effect of buoyancy-driven flow in the gas phase and the thermocapillary-driven flow in the liquid phase. Schiller and Sirignano (1996) extended the same study to include flame spread with forced opposed flow in the gas phase. A transitional flow velocity was found above which an originally uniform spreading flame pulsates. The same type of gas-phase recirculation cell caused by the combination of forced opposed flow, buoyancy-driven flow, and thermocapillary-driven concurrent flow was responsible for the pulsating flame spread. Ross and Miller (1998) and Miller and Ross (1998) performed experimental work that corroborates the computational findings of Schiller, Ross, and Sirignano (1996) and Schiller and Sirignano (1996). Cai, Liu, and Sirignano (2002) developed a more comprehensive three-dimensional model and computer code for the flame spread problem. Many improvements in modeling and numerical algorithms were incorporated in the three-dimensional model. Pools of finite width and length were studied in air channels of prescribed height and width. Significant three-dimensional effects around and along the pool edge were observed. The same three-dimensional code is used to study the detailed effects of pool depth, pool width, opposed air flow velocity, and different levels of air oxygen concentration (Cai, Liu, and Sirignano, 2003). Significant three-dimensional effects showing an unsteady wavy flame front for cases of wide pool width are found for the first time in computation, after being noted previously by experimental observers (Ross and Miller, 1999). Regions of uniform and pulsating flame spread are mapped for the flow conditions of pool depth, opposed flow velocity, initial pool temperature, and air oxygen concentration under both normal and microgravity conditions. Details can be found in Cai et al. (2002, 2003). Experimental results recently performed at NASA Glenn of flame spread across a wide, shallow pool as a function of liquid temperature are also presented here.

Beschreibung: It is well known that buoyancy has a major influence on the flow structure of turbulent nonpremixed jet flames. For example, previous studies have shown that transitional and turbulent jet flames exhibit flame lengths that are as much as a factor of two longer in microgravity than in normal gravity. The objective of this study is to extend these previous studies by investigating both mean and fluctuating characteristics of turbulent nonpremixed jet flames under three different gravity levels (1 g, 20 mg and 100 micrograms). This work is described in more detail elsewhere. In addition, we have recently initiated a new study into the effects of buoyancy on turbulent nonpremixed jet flames in cross-flow (JFICF). Buoyancy has been observed to play a key role in determining the centerline trajectories of such flames.6 The objective of this study is to use the low gravity environment to study the effects of buoyancy on the turbulent characteristics of JFICF.

Beschreibung: The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Beschreibung: Lightweight, flexible probes have been developed for measuring temperatures at multiple locations in tanks that contain possibly pressurized cryogenic fluids. If the fluid in a given tank is subcritical (that is, if it consists of a liquid and its vapor), then in one of two modes of operation, the temperature measurements made by a probe of this type can be used to deduce the approximate level of the liquid. The temperature sensors are silicon diodes located at intervals along a probe. If the probe is to be used to measure a temperature gradient along a given axis in the tank, then the probe must be mounted along that axis. In the non-liquid-level-sensing temperature-measurement mode, a constant small electric current is applied to each diode and the voltage across the diode . a known function of the current and temperature . is measured as an indication of its temperature. For the purpose of this measurement, "small electric current" signifies a current that is not large enough to cause a significant increase in the measured temperature. More specifically, the probe design calls for a current of 10 A, which, in the cryogenic temperature range of interest, generates heat at a rate of only about 0.01 mW per diode. In the liquid-level-sensing mode, one applies a larger current (30 mA) to each diode so as to heat each diode appreciably (with a power of about 36 mW in the temperature range of interest). Because the liquid cools the diode faster than does the vapor, the temperature of the diode is less when diode is immersed in the liquid than when it is above the surface of the liquid. Thus, the temperature (voltage) reading from each diode can be used to determine whether the liquid level is above or below the diode, and one can deduce that the liquid level lies between two adjacent diodes, the lower one of which reads a significantly lower temperature. The aforementioned techniques for measuring temperature and deducing liquid level are not new. What is new here are the designs of the probes and of associated external electronic circuitry. In each probe, the diodes and the lead wires are embedded in a strong, lightweight, flexible polyimide strip. Each probe is constructed as an integral unit that includes a multipin input/output plug or socket for solderless connection of the lead wires to the external circuitry. The polyimide strip includes mounting tabs with holes that can accommodate rivets, screws, or other fasteners. Alternatively, a probe can be mounted by use of an epoxy. A probe can be manufactured to almost any length or width, and the diodes can be embedded at almost any desired locations along and across the polyimide strip. In designing a probe for a specific application, one seeks a compromise between (1) minimizing the number of diodes in order to minimize the complexity of input/output connections and external electronic circuitry while (2) using enough diodes to obtain the required precision. Optionally, to minimize spurious heating of the cryogenic fluid, the external circuitry can be designed to apply power to the probe only during brief measurement intervals. Assuming that the external circuitry is maintained at a steady temperature, a power-on interval of only a few seconds is sufficient to obtain accurate data on temperatures and/or the height of the liquid/vapor interface.

Beschreibung: This paper proposes a task graph (TG) model to represent a single discrete step of multi-block overset grid computational fluid dynamics (CFD) applications. The TG model is then used to not only balance the computational workload across the overset grids but also to reduce inter-grid communication costs. We have developed a set of task assignment heuristics based on the constraints inherent in this class of CFD problems. Two basic assignments, the smallest task first (STF) and the largest task first (LTF), are first presented. They are then systematically costs. To predict the performance of the proposed task assignment heuristics, extensive performance evaluations are conducted on a synthetic TG with tasks defined in terms of the number of grid points in predetermined overlapping grids. A TG derived from a realistic problem with eight million grid points is also used as a test case.

Beschreibung: Currently, a wide variety of flow phenomena are addressed with numerical simulations. Many flow solvers are optimized to simulate a limited spectrum of flow effects effectively, such as single parts of a flow system, but are either inadequate or too expensive to be applied to a very complex problem. As an example, the flow through a gas turbine can be considered. In the compressor and the turbine section, the flow solver has to be able to handle the moving blades, model the wall turbulence, and predict the pressure and density distribution properly. This can be done by a flow solver based on the Reynolds-Averaged Navier-Stokes (RANS) approach. On the other hand, the flow in the combustion chamber is governed by large scale turbulence, chemical reactions, and the presence of fuel spray. Experience shows that these phenomena require an unsteady approach. Hence, for the combustor, the use of a Large Eddy Simulation (LES) flow solver is desirable. While many design problems of a single flow passage can be addressed by separate computations, only the simultaneous computation of all parts can guarantee the proper prediction of multi-component phenomena, such as compressor/combustor instability and combustor/turbine hot-streak migration. Therefore, a promising strategy to perform full aero-thermal simulations of gas-turbine engines is the use of a RANS flow solver for the compressor sections, an LES flow solver for the combustor, and again a RANS flow solver for the turbine section.

Beschreibung: This paper is divided into four parts. First, the level set/vortex sheet method for three-dimensional two-phase interface dynamics is presented. Second, the LSS model for the primary breakup of turbulent liquid jets and sheets is outlined and all terms requiring subgrid modeling are identified. Then, preliminary three-dimensional results of the level set/vortex sheet method are presented and discussed. Finally, conclusions are drawn and an outlook to future work is given.

Beschreibung: Fragmentation plays an important role in a variety of physical, chemical, and geological processes. Examples include atomization in sprays, crushing of rocks, explosion and impact of solids, polymer degradation, etc. Although each individual action of fragmentation is a complex process, the number of these elementary actions is large. It is natural to abstract a simple 'effective' scenario of fragmentation and to represent its essential features. One of the models is the fragmentation under the scaling symmetry: each breakup action reduces the typical length of fragments, r (right arrow) alpha r, by an independent random multiplier alpha (0 〈 alpha 〈 1), which is governed by the fragmentation intensity spectrum q(alpha), integral(sup 1)(sub 0) q(alpha)d alpha = 1. This scenario has been proposed by Kolmogorov (1941), when he considered the breakup of solid carbon particle. Describing the breakup as a random discrete process, Kolmogorov stated that at latest times, such a process leads to the log-normal distribution. In Gorokhovski & Saveliev, the fragmentation under the scaling symmetry has been reviewed as a continuous evolution process with new features established. The objective of this paper is twofold. First, the paper synthesizes and completes theoretical part of Gorokhovski & Saveliev. Second, the paper shows a new application of the fragmentation theory under the scale invariance. This application concerns the turbulent cascade with intermittency. We formulate here a model describing the evolution of the velocity increment distribution along the progressively decreasing length scale. The model shows that when the turbulent length scale gets smaller, the velocity increment distribution has central growing peak and develops stretched tails. The intermittency in turbulence is manifested in the same way: large fluctuations of velocity provoke highest strain in narrow (dissipative) regions of flow.

Beschreibung: With the renewed interest in Cartesian gridding methodologies for the ease and speed of gridding complex geometries in addition to the simplicity of the control volumes used in the computations, it has become important to investigate ways of extending the existing Cartesian grid solver functionalities. This includes developing methods of modeling the viscous effects in order to utilize Cartesian grids solvers for accurate drag predictions and addressing the issues related to the distributed memory parallelization of Cartesian solvers. This research presents advances in two areas of interest in Cartesian grid solvers, viscous effects modeling and MPI parallelization. The development of viscous effects modeling using solely Cartesian grids has been hampered by the widely varying control volume sizes associated with the mesh refinement and the cut cells associated with the solid surface. This problem is being addressed by using physically based modeling techniques to update the state vectors of the cut cells and removing them from the finite volume integration scheme. This work is performed on a new Cartesian grid solver, NASCART-GT, with modifications to its cut cell functionality. The development of MPI parallelization addresses issues associated with utilizing Cartesian solvers on distributed memory parallel environments. This work is performed on an existing Cartesian grid solver, CART3D, with modifications to its parallelization methodology.

Beschreibung: Here we report multiple harmonic theoretical solutions for a complete system of conservation laws, which describe the large-scale coherent dynamics in RTI and RMI for fluids with a finite density ratio in the general three-dimensional case. The analysis yields new properties of the bubble front dynamics. In either RTI or RMI, the obtained dependencies of the bubble velocity and curvature on the density ratio differ qualitatively and quantitatively from those suggested by the models of Sharp (1984), Oron et al. (2001), and Goncharov (2002). We show explicitly that these models violate the conservation laws. For the first time, our theory reveals an important qualitative distinction between the dynamics of the RT and RM bubbles.

Beschreibung: The original aim of this work was to shed some light on the physics of turbulence over rough walls using large-eddy simulations and the suboptimal-control wall boundary conditions introduced by Nicoud et al. It was hoped that, if that algorithm was used to fit the mean velocity profile of the simulations to that of a rough-walled channel, instead of to a smooth one, the wall stresses introduced by the control algorithm would give some indication of what aspects of rough walls are most responsible for the modification of the flow in real turbulence. It was similarly expected that the structure of the resulting velocity fluctuations would share some of the characteristics of rough-walled flows, thus again suggesting what is intrinsic and what is accidental in the effect of geometric wall roughness. A secondary goal was to study the effect of 'unphysical' boundary conditions on the outside flow by observing how a relatively major change of the target velocity profile, and therefore presumably of the applied wall stresses, modifies properties such as the dominant length scales of the velocity fluctuations away from the wall. As will be seen below, this secondary goal grew more important during the course of the study, which was carried out during a short summer visit of the first two authors to the CTR. It became clear that there are open questions about the way in which the control algorithm models the boundary conditions, even for smooth walls, and that these questions make the physical interpretation of the results difficult. Considerable more work in that area seems to be needed before even relatively advanced large-eddy simulations, such as these, can be used to draw conclusions about the physics of wall-bounded turbulent flows. The numerical method is the same as in Nicoud et al. The modifications introduced in the original code are briefly described in section 2, but the original paper should be consulted for a full description of the algorithm. The results are presented in section 3 and summarized in section 4. The elementary properties of turbulence over rough walls which are used in the text have been taken from recent reviews such as Raupach et al. or Jimenez.