ALBERT

Description: No abstract available

Description: An orifice element is commonly used in liquid rocket engine test facilities either as a flow metering device, a damper for acoustic resonance or to provide a large reduction in pressure over a very small distance in the piping system. While the orifice as a device is largely effective in stepping down pressure, it is also susceptible to a wake-vortex type instability that generates pressure fluctuations that propagate downstream and interact with other elements of the test facility resulting in structural vibrations. Furthermore in piping systems an unstable feedback loop can exist between the vortex shedding and acoustic perturbations from upstream components resulting in an amplification of the modes convecting downstream. Such was the case in several tests conducted at NASA as well as in the Ariane 5 strap-on P230 engine in a static firing test where pressure oscillations of 0.5% resulted in 5% thrust oscillations. Exacerbating the situation in cryogenic test facilities, is the possibility of the formation of vapor clouds when the pressure in the wake falls below the vapor pressure leading to a cavitation instability that has a lower frequency than the primary wake-vortex instability. The cavitation instability has the potential for high amplitude fluctuations that can cause catastrophic damage in the facility. In this paper high-fidelity multi-phase numerical simulations of an orifice element are used to characterize the different instabilities, understand the dominant instability mechanisms and identify the tonal content of the instabilities.

Description: A review of recently published coupled radiation and ablation capabilities involving the simulation of hypersonic flowfields relevant to Earth, Mars, or Venus entry is presented. The three fundamental mechanisms of radiation coupling are identified as radiative cooling, precursor photochemistry, and ablation-radiation interaction. The impact of these mechanisms are shown to be significant for a 3 m radius sphere entering Earth at hypothetical Mars return conditions (approximately 15 km/s). To estimate the influence precursor absorption on the radiative flux for a wide range of conditions, a simplified approach is developed that requires only the non-precursor solution. Details of a developed coupled ablation approach, which is capable of treating both massively ablating flowfields in the sublimation regime and weakly ablating diffusion Climited oxidation cases, are presented. A review of the two primary uncoupled ablation approximations, identified as the blowing correction and film coefficient approximations, is made and their impact for Earth and Mars entries is shown to be significant for recession and convective heating predictions. Fully coupled ablation and radiation simulations are presented for the Mars return sphere throughout its entire trajectory. Applying to the Mars return sphere the Pioneer- Venus heritage carbon phenolic heatshield, which has properties available in the open literature, the differences between steady state ablation and coupling to a material response code are shown to be significant.

Description: The SPHERES Slosh Experiment (SSE) is a free floating experimental platform developed for the acquisition of long duration liquid slosh data aboard the International Space Station (ISS). The data sets collected will be used to benchmark numerical models to aid in the design of rocket and spacecraft propulsion systems. Utilizing two SPHERES Satellites, the experiment will be moved through different maneuvers designed to induce liquid slosh in the experiment's internal tank. The SSE has a total of twenty-four thrusters to move the experiment. In order to design slosh generating maneuvers, a parametric study with three maneuvers types was conducted using the General Moving Object (GMO) model in Flow-30. The three types of maneuvers are a translation maneuver, a rotation maneuver and a combined rotation translation maneuver. The effectiveness of each maneuver to generate slosh is determined by the deviation of the experiment's trajectory as compared to a dry mass trajectory. To fully capture the effect of liquid re-distribution on experiment trajectory, each thruster is modeled as an independent force point in the Flow-3D simulation. This is accomplished by modifying the total number of independent forces in the GMO model from the standard five to twenty-four. Results demonstrate that the most effective slosh generating maneuvers for all motions occurs when SSE thrusters are producing the highest changes in SSE acceleration. The results also demonstrate that several centimeters of trajectory deviation between the dry and slosh cases occur during the maneuvers; while these deviations seem small, they are measureable by SSE instrumentation.

Description: This paper describes the development and testing of a scalable thermal control architecture for instruments, subsystems, or systems that must operate in severe space environments with wide variations in sink temperature. The architecture is comprised by linking one or more hot-side variable conductance heat pipes (VCHPs) in series with one or more cold-side loop heat pipes (LHPs). The VCHPs provide wide area heat acquisition, limited distance thermal transport, modest against gravity pumping, concentrated LHP startup heating, and high switching ratio variable conductance operation. The LHPs provide localized heat acquisition, long distance thermal transport, significant against gravity pumping, and high switching ratio variable conductance operation. Combining two variable conductance devices in series ensures very high switching ratio isolation from severe environments like the Earth's moon, where each lunar day spans 15 Earth days (270 K sink, with a surface-shielded/space viewing radiator) and each lunar night spans 15 Earth days (80-100 K radiative sink, depending on location). The single VCHP-single LHP system described herein was developed to maintain thermal control of International Lunar Network (ILN) anchor node lander electronics, but it is also applicable to other variable heat rejection space missions in severe environments. The LHPVCHP system utilizes a stainless steel wire mesh wick ammonia VCHP, a Teflon wick propylene LHP, a pair of one-third square meter high radiators (one capillary-pumped horizontal radiator and a second gravity-fed vertical radiator), a half-meter of transport distance, and a wick-bearing co-located flow regulator (CLFR) to allow operation with a hot (deactivated) radiator. The VCHP was designed with a small reservoir formed by extending the length of its stainless steel heat pipe tubing. The system was able to provide end-to-end switching ratios of 300-500 during thermal vacuum testing at ATK, including 3-5 W/K ON conductance and 0.01 W/K OFF conductance. The test results described herein also include an in-depth analysis of VCHP condenser performance to explain VCHP switching operation in detail. Future multi-VCHP/multi-LHP thermal management system concepts that provide scalability to higher powers/longer transport lengths are also discussed in the paper.

Description: It is known that the adverse effects of shock wave boundary layer interactions in high speed inlets include reduced total pressure recovery and highly distorted flow at the aerodynamic interface plane (AIP). This paper presents a design method for flow control which creates perturbations in geometry. These perturbations are tailored to change the flow structures in order to minimize shock wave boundary layer interactions (SWBLI) inside supersonic inlets. Optimizing the shape of two dimensional micro-size bumps is shown to be a very effective flow control method for two-dimensional SWBLI. In investigating the three dimensional SWBLI, a square duct is employed as a baseline. To investigate the mechanism whereby the geometric elements of the baseline, i.e. the bottom wall, the sidewall and the corner, exert influence on the flow's aerodynamic characteristics, each element is studied and optimized separately. It is found that arrays of micro-size bumps on the bottom wall of the duct have little effect in improving total pressure recovery though they are useful in suppressing the incipient separation in three-dimensional problems. Shaping sidewall geometry is effective in re-distributing flow on the side wall and results in a less distorted flow at the exit. Subsequently, a near 50% reduction in distortion is achieved. A simple change in corner geometry resulted in a 2.4% improvement in total pressure recovery.

Description: Nuclear power provides an enabling capability for NASA missions that might otherwise be constrained by power availability, mission duration, or operational robustness. NASA and the Department of Energy (DOE) are developing fission power technology to serve a wide range of future space uses. Advantages include lower mass, longer life, and greater mission flexibility than competing power system options. Kilowatt-class fission systems, designated "Kilopower," were conceived to address the need for systems to fill the gap above the current 100-Wclass radioisotope power systems being developed for science missions and below the typical 100-kWe-class reactor power systems being developed for human exploration missions. This paper reviews the current fission technology project and examines some Kilopower concepts that could be used to support future science missions or human precursors.

Description: Fixed packed bed reactors are compact, require minimum power and maintenance to operate, and are highly reliable. These features make this technology a highly desirable unit operation for long duration life support systems in space. NASA is developing an ISS experiment to address this technology with particular focus on water reclamation and air revitalization. Earlier research and development efforts funded by NASA have resulted in two hydrodynamic models which require validation with appropriate instrumentation in an extended microgravity environment. To validate these models, the instantaneous distribution of the gas and liquid phases must be measured.Electrical Capacitance Volume Tomography (ECVT) is a non-invasive imaging technology recently developed for multi-phase flow applications. It is based on distributing flexible capacitance plates on the peripheral of a flow column and collecting real-time measurements of inter-electrode capacitances. Capacitance measurements here are directly related to dielectric constant distribution, a physical property that is also related to material distribution in the imaging domain. Reconstruction algorithms are employed to map volume images of dielectric distribution in the imaging domain, which is in turn related to phase distribution. ECVT is suitable for imaging interacting materials of different dielectric constants, typical in multi-phase flow systems. ECVT is being used extensively for measuring flow variables in various gas-liquid and gas-solid flow systems. Recent application of ECVT include flows in risers and exit regions of circulating fluidized beds, gas-liquid and gas-solid bubble columns, trickle beds, and slurry bubble columns. ECVT is also used to validate flow models and CFD simulations. The technology is uniquely qualified for imaging phase concentrations in packed bed reactors for the ISS flight experiments as it exhibits favorable features of compact size, low profile sensors, high imaging speed, and flexibility to fit around columns of various shapes and sizes. ECVT is also safer than other commonly used imaging modalities as it operates in the range of low frequencies (1 MHz) and does not radiate radioactive energy. In this effort, ECVT is being used to image flow parameters in a packed bed reactor for an ISS flight experiment.

Description: Air plasma radiation in Local Thermodynamic Equilibrium (LTE) within cylindrical geometries is studied with an application towards modeling the radiative transfer inside arc-constrictors, a central component of constricted-arc arc jets. A detailed database of spectral absorption coefficients for LTE air is formulated using the NEQAIR code developed at NASA Ames Research Center. The database stores calculated absorption coefficients for 1,051,755 wavelengths between 0.04 m and 200 m over a wide temperature (500K to 15 000K) and pressure (0.1 atm to 10.0 atm) range. The multi-group method for spectral reduction is studied by generating a range of reductions including pure binning and banding reductions from the detailed absorption coefficient database. The accuracy of each reduction is compared to line-by-line calculations for cylindrical temperature profiles resembling typical profiles found in arc-constrictors. It is found that a reduction of only 1000 groups is sufficient to accurately model the LTE air radiation over a large temperature and pressure range. In addition to the reduction comparison, the cylindrical-slab formulation is compared with the finite-volume method for the numerical integration of the radiative flux inside cylinders with varying length. It is determined that cylindrical-slabs can be used to accurately model most arc-constrictors due to their high length to radius ratios.

Description: Updated information on the status of the UH-60A Airloads wind tunnel test data, including recent work performed by NASA.

Description: Lateral nozzle forces are known to cause severe structural damage to any new rocket engine in development. Currently there is no fully coupled computational tool to analyze this fluid/structure interaction process. The objective of this study was to develop a fully coupled aeroelastic modeling capability to describe the fluid/structure interaction process during the transient nozzle operations. The aeroelastic model composes of three components: the computational fluid dynamics component based on an unstructured-grid, pressure-based computational fluid dynamics formulation, the computational structural dynamics component developed in the framework of modal analysis, and the fluid-structural interface component. The developed aeroelastic model was applied to the transient nozzle startup process of the Space Shuttle Main Engine at sea level. The computed nozzle side loads and the axial nozzle wall pressure profiles from the aeroelastic nozzle are compared with those of the published rigid nozzle results, and the impact of the fluid/structure interaction on nozzle side loads is interrogated and presented.

Description: The overall goal of this study was achieved: Replicated the numerical assessment performed by Chen et. al. (2005). Displayed the ability of Thermal Desktop to be coupled with thermo-electrochemical analysis techniques. such that the local heat generated on the cells is a function of the model itself using logic blocks and arrays. Differences in the TD temperature vs. depth of discharge profiles and Chen's was most likely due to differences in two primary areas: Contact regions and conductance values. Differences in density and specific heat values. center dot The model results are highly dependent on the accuracy of the material properties with respect to the multiple layers of an individual cell.

Description: The turbulent flows through a plane channel and a channel with a constriction (2-D hill) are numerically simulated using DNS and RANS calculations. The Navier-Stokes equations in the DNS are solved using a higher order kinetic energy preserving central schemes and a fifth order accurate upwind biased WENO scheme for the space discretization. RANS calculations are performed using the NASA code CFL3D with the komega SST two-equation model and a full Reynolds stress model. Using DNS, the magnitudes of different terms that appear in the enstrophy equation are evaluated. The results show that the dissipation and the diffusion terms reach large values at the wall. All the vortex stretching terms have similar magnitudes within the buffer region. Beyond that the triple correlation among the vorticity and strain rate fluctuations becomes the important kinematic term in the enstrophy equation. This term is balanced by the viscous dissipation. In the separated flow, the triple correlation term and the viscous dissipation term peak locally and balance each other near the separated shear layer region. These findings concur with the analysis of Tennekes and Lumley, confirming that the energy transfer terms associated with the small-scale dissipation and the fluctuations of the vortex stretching essentially cancel each other, leaving an equation for the dissipation that is governed by the large-scale motion.

Description: The functional equivalence of the unstructured grid code FUN3D to the the structured grid code LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. Examples from an existing suite of regression tests are used to demonstrate the functional equivalence, encompassing various thermochemical models and vehicle configurations. Algorithm modifications required for the node-based unstructured grid code (FUN3D) to reproduce functionality of the cell-centered structured code (LAURA) are also documented. Challenges associated with computation on tetrahedral grids versus computation on structured-grid derived hexahedral systems are discussed.

Description: No abstract available

Description: The fully implicit ablation and thermal response code simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid.This work describes new modeling capabilities that are added to a special version of code. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. In the first group, an Orion flight environment for a proposed lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas model.

Description: With the era of full utilization of the ISS now upon us, this presentation will discuss some of the highest-priority areas for two-phase flow systems with thermal control applications. These priorities are guided by recommendations of a 2011 NRC Decadal Survey report, Recapturing a Future for Space Exploration, Life and Physical Sciences for a New Era as well as an internal NASA exercise in response to the NRC report conducted in early 2012. Many of these proposals are already in various stages of development, while others are still conceptual.

Description: Thin Film Heat Flux Sensors: Design and Methodology: (1) Heat flux is one of a number of parameters, together with pressure, temperature, flow, etc. of interest to engine designers and fluid dynamists, (2) The measurement of heat flux is of interest in directly determining the cooling requirements of hot section blades and vanes, and (3)In addition, if the surface and gas temperatures are known, the measurement of heat flux provides a value for the convective heat transfer coefficient that can be compared with the value provided by CFD codes.

Description: Upward flame spread experiments were conducted on a thin fabric cloth consisting of 75% cotton and 25% fiberglass. The sample is sandwiched symmetrically with stainless steel plates with the exposed width varying between 2 to 8.8 cm from test to test and 〉1.5m tall. The bottom edge was ignited resulting in a symmetric two sided flame. For the narrower samples (. 5cm), two sided flame growth would proceed until reaching some limiting value (15-30 cm depending on sample width). Fluctuation or instability of the flame base on one side would initially become visible and then the flame base would retreat downstream and cause extinguishment on one side. Detailed examination of the still images shows that the fuel continues to vaporize from the extinguished side due to the thermally thin nature of the fuel. But, due to the remaining inert fiberglass mesh, which acts as a flashback arrestor, the extinguished side was not able to be reignited by the remaining flame. The remaining flame would then shrink in length due to the reduced heat transfer to the solid to a shorter length. The one-sided flame will spread stably with a constant speed and a constant flame length to the end of the sample. A constant length flame implies that the pyrolysis front and the burnt out fronts move at the same speed. For the wider samples (. 7cm), no one-sided extinction is observed. Two-sided flames spread all the way to the top of the sample. For these wider widths, the flames are still growing and have not reached their limiting length if it exists. Care was taken to minimize the amount of non-symmetries in the experimental configuration. Repeated tests show that blow-off can occur on either side of the sample. The flame growth is observed to be very symmetric during the growth phase and grew to significant length (〉10cm) before extinction of the flame on one side. Our proposed explanation of this unusual phenomenon (i.e. stronger two ]sided flame cannot exist but weaker one-sided flame can) is as follows: The observed one-sided extinction is a blow- off induced by buoyant entrainment. It is known that the flammable diffusion flame regime is bounded by quenching and blow ]off limits when varying incoming air velocity. The narrowest samples tested (between 2 and 5 cm) begin within the flammable range, but as the flame grows, the buoyancy driven air velocity increases at the neighborhood of the flame base. The initially stable flame crosses the extinguishment boundary resulting in a flame blow-off. When one-side of the flame extinguishes, the remaining side shrinks due to the reduced heat transfer to the solid. This reduces the induced velocity and the flame becomes stable. It is proposed that this may have implications to upward flame growth beyond this experiment.

Description: The present effort aims to strengthen modeling work conducted at the NASA Ames Research Center by measuring the critical plasma electron characteristics within and slightly outside of an arc jet plasma column. These characteristics are intended to give physical insights while assisting in the formulation of boundary conditions to validate full scale simulations. Single and triple Langmuir probes have been used to achieve estimates of the electron temperature (T(sub e)), electron number density (n(sub e)) and plasma potential (outside of the plasma column) as probing location is varied radially from the flow centerline. Both the electron temperature and electron number density measurements show a large dependence on radial distance from the plasma column centerline with T(sub e) approx. = (3 - 12 eV and n(sub e) approx. = 10(exp 12) - 10(exp 14)/cu cm.

Description: This work uses a discontinuous-Galerkin spectral-element method (DGSEM) to solve the compressible Navier-Stokes equations [1{3]. The inviscid ux is computed using the approximate Riemann solver of Roe [4]. The viscous fluxes are computed using the second form of Bassi and Rebay (BR2) [5] in a manner consistent with the spectral-element approximation. The method of lines with the classical 4th-order explicit Runge-Kutta scheme is used for time integration. Results for polynomial orders up to p = 15 (16th order) are presented. The code is parallelized using the Message Passing Interface (MPI). The computations presented in this work are performed using the Sandy Bridge nodes of the NASA Pleiades supercomputer at NASA Ames Research Center. Each Sandy Bridge node consists of 2 eight-core Intel Xeon E5-2670 processors with a clock speed of 2.6Ghz and 2GB per core memory. On a Sandy Bridge node the Tau Benchmark [6] runs in a time of 7.6s.

Description: Thermal protection materials and systems (TPS) protect vehicles from the heat generated when entering a planetary atmosphere. NASA has developed many TPS systems over the years for vehicle ranging from planetary probes to crewed vehicles. The goal for all TPS is efficient and reliable performance. Efficient means using the right material for the environment and minimizing the mass of the heat shield without compromising safety. Efficiency is critical if the payload such as science experiments is to be maximized on a particular vehicle. Reliable means that we understand and can predict performance of the material. Although much characterization and testing of materials is performed to qualify and certify them for flight, it is not possible to completely recreate the reentry conditions in test facilities, and flight-testing

Description: Analysis tools are needed to investigate aerodynamic performance of Variable-Speed Power Turbines (VSPT) for rotorcraft applications. The VSPT operates at low Reynolds numbers (transitional flow) and over a wide range of incidence. Previously, the capability of a published three-equation turbulence model to predict accurately the transition location for three-dimensional heat transfer problems was assessed. In this paper, the results of a post-diction exercise using a three-dimensional flow in a transonic linear cascade comprising VSPT blading are presented. The measured blade pressure distributions and exit total pressure and flow angles for two incidence angles corresponding to cruise (i = 5.8deg) and takeoff (i = -36.7deg) were used for this study. For the higher loading condition of cruise and the negative incidence condition of takeoff, overall agreement with data may be considered satisfactory but areas of needed improvement are also indicated.

Description: The goal of this study was to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC) installed in Node 3. Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow enable identifying the paths of water transport. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 5-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain.

Description: Five different central difference schemes, based on a conservative differencing form of the Kennedy and Gruber skew-symmetric scheme, were compared with six different upwind schemes based on primitive variable reconstruction and the Roe flux. These eleven schemes were tested on a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem and a turbulent channel flow problem. The central schemes were generally very accurate and stable, provided the grid stretching rate was kept below 10%. As near-DNS grid resolutions, the results were comparable to reference DNS calculations. At coarser grid resolutions, the need for an LES SGS model became apparent. There was a noticeable improvement moving from CD-2 to CD-4, and higher-order schemes appear to yield clear benefits on coarser grids. The UB-7 and CU-5 upwind schemes also performed very well at near-DNS grid resolutions. The UB-5 upwind scheme does not do as well, but does appear to be suitable for well-resolved DNS. The UF-2 and UB-3 upwind schemes, which have significant dissipation over a wide spectral range, appear to be poorly suited for DNS or LES.

Description: The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel ow problem.

Description: This paper compares computational and experimental data from the Aeroelastic Prediction Workshop (AePW) held in April 2012. This workshop was designed as a series of technical interchange meetings to assess the state of the art of computational methods for predicting unsteady flowfields and static and dynamic aeroelastic response. The goals are to provide an impartial forum to evaluate the effectiveness of existing computer codes and modeling techniques to simulate aeroelastic problems and to identify computational and experimental areas needing additional research and development. Three subject configurations were chosen from existing wind-tunnel data sets where there is pertinent experimental data available for comparison. Participant researchers analyzed one or more of the subject configurations, and results from all of these computations were compared at the workshop.

Description: A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

Description: The Aeroelastic Prediction Workshop brought together an international community of computational fluid dynamicists as a step in defining the state of the art in computational aeroelasticity. The Rectangular Supercritical Wing (RSW) was chosen as the first configuration to study due to its geometric simplicity, perceived simple flow field at transonic conditions and availability of an experimental data set containing forced oscillation response data. Six teams performed analyses of the RSW; they used Reynolds-Averaged Navier-Stokes flow solvers exercised assuming that the wing had a rigid structure. Both steady-state and forced oscillation computations were performed by each team. The results of these calculations were compared with each other and with the experimental data. The steady-state results from the computations capture many of the flow features of a classical supercritical airfoil pressure distribution. The most dominant feature of the oscillatory results is the upper surface shock dynamics. Substantial variations were observed among the computational solutions as well as differences relative to the experimental data. Contributing issues to these differences include substantial wind tunnel wall effects and diverse choices in the analysis parameters.

Description: Measurements have been conducted at the University of Virginia Supersonic Combustion Facility in configuration C of the dual-mode scramjet. This is a continuation of previously published works on configuration A. The scramjet is hydrogen fueled and operated at two equivalence ratios, one representative of the scram mode and the other of the ram mode. Dual-pump CARS was used to acquire the mole fractions of the major species as well as the rotational and vibrational temperatures of N2. Developments in methods and uncertainties in fitting CARS spectra for vibrational temperature are discussed. Mean quantities and the standard deviation of the turbulent fluctuations at multiple planes in the flow path are presented. In the scram case the combustion of fuel is completed before the end of the measurement domain, while for the ram case the measurement domain extends into the region where the flow is accelerating and combustion is almost completed. Higher vibrational than rotational temperature is observed in those parts of the hot combustion plume where there is substantial H2 (and hence chemical reaction) present.

Description: A technique for predicting supersonic pitch damping of blunt axisymmetric bodies from static CFD data is presented. The contributions to static pitching moment due to forebody and aftbody pressure distributions are broken out and considered separately. The one-dimension moment equation is cast to model the separate contributions from forebody and aftbody pressures with no traditional damping term included. The aftbody contribution to pitching moment is lagged by a phase angle of the natural oscillation period. This lag represents the time for aftbody wake structures to equilibrate while the body is oscillation. The characteristic equation of this formulation indicates that the lagged backshell moment adds a damping moment equivalent in form to a constant pitch damping term. CFD calculations of the backshell's contribution to the static pitching moment for a range of angles-of-attack is used to predict pitch damping coefficients. These predictions are compared with ballistic range data taken of the Mars Exploration Rover (MER) capsule and forced oscillation data of the Mars Viking capsule. The lag model appears to capture dynamic stability variation due to backshell geometry as well as Mach number.

Description: Non-intrusive hydroxyl radical (OH) planar laser-induced fluorescence (PLIF) measurements were obtained in configuration C of the University of Virginia supersonic combustion experiment. The combustion of hydrogen fuel injected through an unswept compression ramp into a supersonic cross-flow was imaged over a range of streamwise positions. Images were corrected for optical distortion, variations in the laser sheet profile, and different camera views. Results indicate an effect of fuel equivalence ratio on combustion zone shape and local turbulence length scale. The streamwise location of the reaction zone relative to the fuel injector was also found to be sensitive to the fuel equivalence ratio. The flow boundary conditions in the combustor section, which are sensitive to the fuel flow rate, are believed to have caused this effect. A combination of laser absorption and radiative trapping effects are proposed to have caused asymmetry observed in the images. The results complement previously published OH PLIF data obtained for configuration A along with other non-intrusive measurements to form a database for computational fluid dynamics (CFD) model validation.

Description: The effects of geometric changes on the amplification of first mode instability waves in an external supersonic boundary layer were investigated using numerical techniques. Boundary layer stability was analyzed at Mach 6 conditions similar to freestream conditions obtained in quiet ground test facilities so that results obtained in this study may be applied to future test article design to measure first mode instability waves. The DAKOTA optimization software package was used to optimize an axisymmetric geometry to maximize the amplification of the waves at first mode frequencies as computed by the 2D STABL hypersonic boundary layer stability analysis tool. First, geometric parameters such as nose radius, cone half angle, vehicle length, and surface curvature were examined separately to determine the individual effects on the first mode amplification. Finally, all geometric parameters were allowed to vary to produce a shape optimized to maximize the amplification of first mode instability waves while minimizing the amplification of second mode instability waves. Since first mode waves are known to be most unstable in the form of oblique wave, the geometries were optimized using a broad range of wave frequencies as well as a wide range of oblique wave angles to determine the geometry that most amplifies the first mode waves. Since first mode waves are seen most often in flows with low Mach numbers at the edge of the boundary layer, the edge Mach number for each geometry was recorded to determine any relationship between edge Mach number and the stability of first mode waves. Results indicate that an axisymmetric cone with a sharp nose and a slight flare at the aft end under the Mach 6 freestream conditions used here will lower the Mach number at the edge of the boundary layer to less than 4, and the corresponding stability analysis showed maximum first mode N factors of 3.

Description: Heat flux sensors are routinely used in arc jet facilities to determine heat transfer rates from plasma plume. The goal of this study is to assess the impact of surface composition changes on these heat flux sensors. Surface compositions can change due to oxidation and material deposition from the arc jet. Systematic surface analyses of the sensors were conducted before and after exposure to plasma. Currently copper is commonly used as surface material. Other surface materials were studied including nickel, constantan gold, platinum and silicon dioxide. The surfaces were exposed to plasma between 0.3 seconds and 3 seconds. Surface changes due to oxidation as well as copper deposition from the arc jets were observed. Results from changes in measured heat flux as a function of surface catalycity is given, along with a first assessment of enthalpy for these measurements. The use of cupric oxide is recommended for future heat flux measurements, due to its consistent surface composition arc jets.

Description: The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Description: Future space exploration missions require advanced thermal control systems (TCS) to dissipate heat from spacecraft, rovers, or habitats operating in environments that can vary from extremely hot to extremely cold. A lightweight, reliable TCS is being developed to effectively control cabin and equipment temperatures under widely varying heat loads and ambient temperatures. The system uses freeze-tolerant radiators, which eliminate the need for a secondary circulation loop or heat pipe systems. Each radiator has a self-regulating variable thermal conductance to its ambient environment. The TCS uses a nontoxic, water-based working fluid that is compatible with existing lightweight aluminum heat exchangers. The TCS is lightweight, compact, and requires very little pumping power. The critical characteristics of the core enabling technologies were demonstrated. Functional testing with condenser tubes demonstrated the key operating characteristics required for a reliable, freeze-tolerant TCS, namely (1) self-regulating thermal conductance with short transient responses to varying thermal loads, (2) repeatable performance through freeze-thaw cycles, and (3) fast start-up from a fully frozen state. Preliminary coolant tests demonstrated that the corrosion inhibitor in the water-based coolant can reduce the corrosion rate on aluminum by an order of magnitude. Performance comparison with state-of-the-art designs shows significant mass and power saving benefits of this technology.

Description: Reactor fuel rod surface area that is perpendicular to coolant flow direction (+S) i.e. perpendicular to the P creates areas of coolant stagnation leading to increased coolant temperatures resulting in localized changes in fluid properties. Changes in coolant fluid properties caused by minor increases in temperature lead to localized reductions in coolant mass flow rates leading to localized thermal instabilities. Reductions in coolant mass flow rates result in further increases in local temperatures exacerbating changes to coolant fluid properties leading to localized thermal runaway. Unchecked localized thermal runaway leads to localized fuel melting. Reactor designs with randomized flow paths are vulnerable to localized thermal instabilities, localized thermal runaway, and localized fuel melting.

Description: A review of literature associated with Pebble Bed and Particle Bed reactor core research has revealed a systemic problem inherent to reactor core concepts which utilize randomized rather than structured coolant channel flow paths. For both the Pebble Bed and Particle Bed Reactor designs; case studies reveal that for indeterminate reasons, regions within the core would suffer from excessive heating leading to thermal runaway and localized fuel melting. A thermal Computational Fluid Dynamics model was utilized to verify that In both the Pebble Bed and Particle Bed Reactor concepts randomized coolant channel pathways combined with localized high temperature regions would work together to resist the flow of coolant diverting it away from where it is needed the most to cooler less resistive pathways where it is needed the least. In other words given the choice via randomized coolant pathways the reactor coolant will take the path of least resistance, and hot zones offer the highest resistance. Having identified the relationship between randomized coolant channel pathways and localized fuel melting it is now safe to assume that other reactor concepts that utilize randomized coolant pathways such as the foam core reactor are also susceptible to this phenomenon.

Description: NASA is developing new portable life support system (PLSS) technologies, which it is demonstrating in an unmanned ground based prototype unit called PLSS 2.0. One set of technologies within the PLSS provides suitable ventilation to an astronaut while on an EVA. A new component within the ventilation gas loop is a liquid-to-gas heat exchanger to transfer excess heat from the gas to the thermal control system s liquid coolant loop. A unique bench top prototype heat exchanger was built and tested for use in PLSS 2.0. The heat exchanger was designed as a counter-flow, compact plate fin type using stainless steel. Its design was based on previous compact heat exchangers manufactured by United Technologies Aerospace Systems (UTAS), but was half the size of any previous heat exchanger model and one third the size of previous liquid-to-gas heat exchangers. The prototype heat exchanger was less than 40 cubic inches and weighed 2.57 lb. Performance of the heat exchanger met the requirements and the model predictions. The water side and gas side pressure drops were less 0.8 psid and 0.5 inches of water, respectively, and an effectiveness of 94% was measured at the nominal air side pressure of 4.1 psia.

Description: Future exploration missions require the development of a new liquid cooling garment (LCG) that offers greater system reliability, is more comfortable, and maximizes thermal performance. To inform the development of a future LCG a thermal performance test was conducted to evaluate three factors: (1) the effect of the thermal comfort undergarment (TCU) on tactile and thermal comfort, (2) the comparable thermal performance of an CSAFE developed engineering evaluation unit (EEU) LCG, which uses a commercial-off-the-shelf (COTS) wicking garment as the base, and (3) the performance of a torso or upper body only LCG configuration to evaluate a proposed auxiliary loop configuration. To evaluate the thermal performance of each configuration a metabolic suit test was conducted, utilizing suited subjects to generate metabolic heat by walking on a treadmill at various speeds. Three (3) test subjects of similar height and weight produced a metabolic load for five tests by either resting (300-600 BTU/hr), walking at a slow pace (1200 BTU/hr), and walking at a brisk pace (2200 BTU/hr). During the test, data was collected that would allow us to track the heat transfer to the LCG and ventilation system to determine the thermal performance of the LCG configurations. Four different test configurations were tested, with one configuration tested twice. The test results show that the CSAFE EEU LCG and EMU LCG had comparable performance. The testing also showed that an auxiliary loop LCG, sized similarly to the shirt-only configuration, should provide adequate cooling for contingency scenarios. Finally, the testing showed the previous analysis that assumed a UA deterioration from the TCU was too conservative and the TCU may prove to be acceptable for future development with additional analysis and testing.

Description: No abstract available

Description: A validation study of Computational Fluid Dynamics (CFD) for Supersonic Retropropulsion (SRP) was conducted using three Navier-Stokes flow solvers (DPLR, FUN3D, and OVERFLOW). The study compared results from the CFD codes to each other and also to wind tunnel test data obtained in the NASA Ames Research Center 90 70 Unitary PlanWind Tunnel. Comparisons include surface pressure coefficient as well as unsteady plume effects, and cover a range of Mach numbers, levels of thrust, and angles of orientation. The comparisons show promising capability of CFD to simulate SRP, and best agreement with the tunnel data exists for the steadier cases of the 1-nozzle and high thrust 3-nozzle configurations.

Description: Planar laser-induced fluorescence (PLIF) of naturally occurring nitric oxide (NO) has been used to obtain instantaneous flow visualization images, and to make both radial and axial velocity measurements in the HYMETS (Hypersonic Materials Environmental Test System) 400 kW arc-heated wind tunnel at NASA Langley Research Center. This represents the first application of NO PLIF flow visualization in HYMETS. Results are presented at selected facility run conditions, including some in a simulated Earth atmosphere (75% nitrogen, 20% oxygen, 5% argon) and others in a simulated Martian atmosphere (71% carbon dioxide, 24% nitrogen, 5% argon), for specific bulk enthalpies ranging from 6.5 MJ/kg to 18.4 MJ/kg. Flow visualization images reveal the presence of large scale unsteady flow structures, and indicate nitric oxide fluorescence signal over more than 70% of the core flow for specific bulk enthalpies below about 11 MJ/kg, but over less than 10% of the core flow for specific bulk enthalpies above about 16 MJ/kg. Axial velocimetry was performed using molecular tagging velocimetry (MTV). Axial velocities of about 3 km/s were measured along the centerline. Radial velocimetry was performed by scanning the wavelength of the narrowband laser and analyzing the resulting Doppler shift. Radial velocities of +/- 0.5 km/s were measured.

Description: Under the Advanced Exploration Systems (AES) Spacecraft Fire Safety Demonstration Project (SFSDP), as a risk mitigation activity in support of the development of a large-scale fire demonstration experiment in microgravity, flame-spread tests were conducted in normal gravity on thin, cellulose-based fuels in a sealed chamber. The primary objective of the tests was to measure pressure rise in a chamber as sample material, burning direction (upward/downward), total heat release, heat release rate, and heat loss mechanisms were varied between tests. A Design of Experiments (DOE) method was imposed to produce an array of tests from a fixed set of constraints and a coupled response model was developed. Supplementary tests were run without experimental design to additionally vary select parameters such as initial chamber pressure. The starting chamber pressure for each test was set below atmospheric to prevent chamber overpressure. Bottom ignition, or upward propagating burns, produced rapid acceleratory turbulent flame spread. Pressure rise in the chamber increases as the amount of fuel burned increases mainly because of the larger amount of heat generation and, to a much smaller extent, due to the increase in gaseous number of moles. Top ignition, or downward propagating burns, produced a steady flame spread with a very small flat flame across the burning edge. Steady-state pressure is achieved during downward flame spread as the pressure rises and plateaus. This indicates that the heat generation by the flame matches the heat loss to surroundings during the longer, slower downward burns. One heat loss mechanism included mounting a heat exchanger directly above the burning sample in the path of the plume to act as a heat sink and more efficiently dissipate the heat due to the combustion event. This proved an effective means for chamber overpressure mitigation for those tests producing the most total heat release and thusly was determined to be a feasible mitigation strategy to incorporate into the microgravity experiment.

Description: This paper considers the small amplitude unsteady motion of an inviscid non-heat conducting compressible fluid on a transversely sheared mean flow. It extends a previous result given in Goldstein (1978(b) and 1979(a)) which shows that the hydrodynamic component of the motion is determined by two arbitrary convected quantities in the absence of solid surfaces or other external sources. The result is important because it can be used to specify appropriate boundary conditions for unsteady surface interaction problems on transversely sheared mean flows in the same way that the vortical component of the Kovasznay (1953) decomposition is used to specify these conditions for surface interaction problems on uniform mean flows. But unlike the Kovasznay (1953) case the arbitrary convected quantities no longer bear a simple relation to the physical variables. One purpose of this paper is to derive a formula that relates these quantities to the (physically measurable) vorticity and pressure fluctuations in the flow.

Description: This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintain ability by eliminating the requirement for problem dependent recompilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multi-physics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the Fun3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Description: We report a high-speed laser diagnostics technique based on point-wise linear Raman spectroscopy for measuring the frequency content of a CH4-air premixed flame stabilized behind a circular bluff body. The technique, which primarily employs a Nd:YLF pulsed laser and a fast image-intensified CCD camera, successfully measures the time evolution of scalar parameters (N2, O2, CH4, and H2O) in the vortex-induced flame instability at a data rate of 1 kHz. Oscillation of the V-shaped flame front is quantified through frequency analysis of the combustion species data and their correlations. This technique promises to be a useful diagnostics tool for combustion instability studies.

Description: Supplemental information includes input and output data files for Generalized Fluid System Simulation Program Version 6.0.

Description: NASA Glenn Research Center (GRC) has developed the LERCHP code. The PC-based LERCHP code can be used to predict the steady-state performance of heat pipes, including the determination of operating temperature and operating limits which might be encountered under specified conditions. The code contains a vapor flow algorithm which incorporates vapor compressibility and axially varying heat input. For the liquid flow in the wick, Darcy s formula is employed. Thermal boundary conditions and geometric structures can be defined through an interactive input interface. A variety of fluid and material options as well as user defined options can be chosen for the working fluid, wick, and pipe materials. This report documents the current effort at GRC to update the LERCHP code for operating in a Microsoft Windows (Microsoft Corporation) environment. A detailed analysis of the model is presented. The programming architecture for the numerical calculations is explained and flowcharts of the key subroutines are given

Description: No abstract available

Description: Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This proposal describes an approach to validate the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft. The research described here is absolutely cutting edge. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional"validation by test only'' mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computationaf Fluid Dynamics can be used to veritY these requirements; however, the model must be validated by test data. The proposed research project includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT and OPEN FOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid . . . Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. To date, the author is the only person to look at the uncertainty in the entire computational domain. For the flow regime being analyzed (turbulent, threedimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Description: The effects of fin-leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9deg wedge generated a planar oblique shock at 16.7deg to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0deg (normal to the free-stream) to 15deg and 25deg swept forward. Global temperature data was obtained from the surface of the fused silica fins through phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite one-dimensional method.

Description: Advances in the understanding of multiphase flow characteristics under variable gravity conditions will ultimately lead to improved and as of yet unknown process designs for advanced space missions. Such novel processes will be of paramount importance to the success of future manned space exploration as we venture into our solar system and beyond. In addition, because of the ubiquitous nature and vital importance of biological and environmental processes involving airwater mixtures, knowledge gained about fundamental interactions and the governing properties of these mixtures will clearly benefit the quality of life here on our home planet. The techniques addressed in the current research involving multiphase transport in porous media and gas-liquid phase separation using capillary pressure gradients are also a logical candidate for a future International Space Station (ISS) flight experiment. Importantly, the novel and potentially very accurate Lattice-Boltzmann (LB) modeling of multiphase transport in porous media developed in this work offers significantly improved predictions of real world fluid physics phenomena, thereby promoting advanced process designs for both space and terrestrial applications.This 3-year research effort has culminated in the design and testing of a zero-g demonstration prototype. Both the hydrophilic (glass) and hydrophobic (Teflon) media Capillary Pressure Gradient (CPG) cartridges prepared during the second years work were evaluated. Results obtained from ground testing at 1-g were compared to those obtained at reduced gravities spanning Martian (13-g), Lunar (16-g) and zero-g. These comparisons clearly demonstrate the relative strength of the CPG phenomena and the efficacy of its application to meet NASAs unique gas-liquid separation (GLS) requirements in non-terrestrial environments.LB modeling software, developed concurrently with the zero-g test effort, was shown to accurately reproduce observed CPG driven gas-liquid separation phenomena. The design and fabrication of a micropost plate-lamina Hele-Shaw (HS) cell was performed which served as a computationally attainable geometric structure facilitating direct comparison between physical phenomena observed in our laboratory and the LB software predictions.

Description: The optimum radiator configuration in hot lunar thermal environments is one in which the radiator is parallel to the ground and has no view to the hot lunar surface. However, typical spacecraft configurations have limited real estate available for top-mounted radiators, resulting in a desire to use the spacecraft's vertically oriented sides. Vertically oriented, flat panel radiators will have a large view factor to the lunar surface, and thus will be subjected to significant incident lunar infrared heat. Consequently, radiator fluid temperatures will need to exceed approximately 325 K (assuming standard spacecraft radiator optical properties) in order to provide positive heat rejection at lunar noon. Such temperatures are too high for crewed spacecraft applications in which a heat pump is to be avoided. A recent study of vertically oriented radiator configurations subjected to lunar noon thermal environments led to the discovery of a novel radiator concept that yielded positive heat rejection at lower fluid temperatures. This radiator configuration, called the Intense Thermal Infrared Reflector (ITIR), has exhibited superior performance to all previously analyzed concepts in terms of heat rejection in the lunar noon thermal environment. A key benefit of ITIR is the absence of louvers or other moving parts and its simple geometry (no parabolic shapes). ITIR consists of a specularly reflective shielding surface and a diffuse radiating surface joined to form a horizontally oriented V-shape (shielding surface on top). The point of intersection of these surfaces is defined by two angles, those which define the tilt of each surface with respect to the local horizontal. The optimum set of these angles is determined on a case-by-case basis. The idea assumes minimal conductive heat transfer between shielding and radiating surfaces, and a practical design would likely stack sets of these surfaces on top of one another to reduce radiator thickness.

Description: A process is provided for producing a salt pill for use in very low temperature adiabatic demagnetization refrigerators (ADRs). The method can include providing a thermal bus in a housing. The thermal bus can include an array of thermally conductive metal conductors. A hydrated salt can be grown on the array of thermally conductive metal conductors. Thermal conductance can be provided to the hydrated salt.

Description: The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependant flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 25 demonstrated example problems.

Description: A short term effort (3-months) was undertaken to demonstrate the feasibility of a novel method to locally control the heat transfer rate and demonstrate the potential to achieve a turndown ratio of approximately 10:1. The technology had to be demonstrated to be at a TRL of 2-3, with a plan to advance it to a TRL 5-6. Here, we show that the concept recently developed in our laboratory, namely the pinned-contact, double droplet switch made by overfilling a hole drilled in a suitable substrate can be implemented with a low-melting temperature metal. When toggled near a second substrate, a liquid bridge can be reversibly connected or disconnected, on demand. We have shown experimentally that liquid-metal bridge switches can be made from gallium with a suitable choice of substrate materials, activation strategies, and control techniques. Individual as well as arrays of gallium bridge switches were shown to be feasible and can be robustly controlled. The very short response time of the bridge connection and disconnection (on the order of 1 millisecond) provides for utility in a wide range of applications. The liquid bridge switches may be controlled actively or passively. We have shown through computations and analysis that liquid bridge switches provide locally large turndown ratios (on the order of 103:1), so a relatively sparse packing of them would be needed to obtain the desired turndown ratio of 10:1. For the laboratory demonstrations, pressure activation was utilized. Simple designs for a passive control strategy are presented which are highly attractive for several reasons, including i) large turndown ratio, ii) no solid-moving parts, and iii) stable operation. Finally, we note that passive systems do not require any electronics for their control. This along with the relatively small molecular weight of candidate materials for the system, makes for a robust design outside of Earth?s magnetic field, where spacecraft are subject to significant radiation bombardment.

Description: Computational Fluid Dynamics (CFD) is a powerful and versatile tool simulating fluid and thermal environments of launch and re-entry vehicles alike. Where it excels in power and accuracy, however, it lacks in speed. An alternative tool for this purpose is known as MINIVER, an aeroheating code widely used by NASA and within the aerospace industry. Capable of providing swift, reasonably accurate approximations of the fluid and thermal environment of launch vehicles, MINIVER is used where time is of the essence and accuracy need not be exact. However, MINIVER is an old, aging tool: running on a user-unfriendly, legacy command-line interface, it is difficult for it to keep pace with more modem software tools. Florida Institute of Technology was tasked with the construction of a new Graphical User Interface (GUI) that implemented the legacy version's capabilities and enhanced them with new tools and utilities. This thesis provides background to the legacy version of the program, the progression and final version of a modem user interface, and benchmarks to demonstrate its usefulness.

Description: The principal purpose of this report is to describe a validation process that can be performed in part on the ground prior to launch, and in space for the IVGEN system. The general approach taken is derived from standard pharmaceutical industry validation schemes modified to fit the special requirements of in-space usage.

Description: Planar, semiconductor heat arrays have been previously proposed and developed; however, this design makes use of a novel, microscale black silicon wick structure that provides increased capillary pumping pressure of the internal working fluid, resulting in increased effective thermal conductivity of the device, and also enables operation of the device in any orientation with respect to the gravity vector. In a heat pipe, the efficiency of thermal transfer from the case to the working fluid is directly proportional to the surface area of the wick in contact with the fluid. Also, the primary failure mechanism for heat pipes operating within the temperature range of interest is inadequate capillary pressure for the return of fluid from the condenser to the wick. This is also what makes the operation of heat pipes orientation-sensitive. Thus, the two primary requirements for a good wick design are a large surface area and high capillary pressure. Surface area can be maximized through nanomachined surface roughening. Capillary pressure is largely driven by the working fluid and wick structure. The proposed nanostructure wick has characteristic dimensions on the order of tens of microns, which promotes menisci of very small radii. This results in the possibility of enormous pumping potential due to the inverse proportionality with radius. Wetting, which also enhances capillary pumping, can be maximized through growth of an oxide layer or material deposition (e.g. TiO2) to create a superhydrophilic surface.

Description: The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Description: A free piston Stirling machine including a thermal buffer tube extending from the machine's expansion space and surrounded by its heat rejector and its regenerator, a displacer cylinder extending from the thermal buffer tube to the compression space and surrounded by the heat rejecting heat exchanger, and a displacer that reciprocates within an excursion limit that extends into the regenerator by no more than 20% of the length of the regenerator during normal operation and preferably within excursion limits that are substantially the length of the heat rejector.

Description: A prediction of whether a point on a computer-generated surface is adjacent to laminar or turbulent flow is made using a transition prediction technique. A plurality of boundary-layer properties at the point are obtained from a steady-state solution of a fluid flow in a region adjacent to the point. A plurality of instability modes are obtained, each defined by one or more mode parameters. A vector of regressor weights is obtained for the known instability growth rates in a training dataset. For each instability mode in the plurality of instability modes, a covariance vector is determined, which is the covariance of a predicted local growth rate with the known instability growth rates. Each covariance vector is used with the vector of regressor weights to determine a predicted local growth rate at the point. Based on the predicted local growth rates, an n-factor envelope at the point is determined.

Description: A method and apparatus for the measurement of thermal conductivity combines the following capabilities: 1) measurements of very small specimens; 2) measurements of specimens with thermal conductivity on the same order of that as air; and, 3) the ability to use air as a reference material. Care is taken to ensure that the heat flow through the test specimen is essentially one-dimensional. No attempt is made to use heated guards to minimize the flow of heat from the hot plate to the surroundings. Results indicate that since large correction factors must be applied to account for guard imperfections when specimen dimensions are small, simply measuring and correcting for heat from the heater disc that does not flow into the specimen is preferable.

Description: A parallel-plate system collects data used to determine viscosity. A first plate is coupled to a translator so that the first plate can be moved along a first direction. A second plate has a pendulum device coupled thereto such that the second plate is suspended above and parallel to the first plate. The pendulum device constrains movement of the second plate to a second direction that is aligned with the first direction and is substantially parallel thereto. A force measuring device is coupled to the second plate for measuring force along the second direction caused by movement of the second plate.

Description: A method for making a microporous membrane comprises the steps of: providing a plurality of carbon nanotubes having a hollow interior diameter of 20 Angstroms or less; sonicating the plurality of carbon nanotubes utilizing a solution comprising deionized, distilled water and a surfactant that coats at least one of the plurality of carbon nanotubes; collecting the coated carbon nanotubes; forming a matrix that supports the plurality of carbon nanotubes; embedding the coated carbon nanotubes into the matrix; rinsing the coated nanotubes to remove at least a portion of the surfactant; curing the nanotube-matrix assembly; and cutting the nanotube-matrix assembly to a particular thickness so as to open the ends of the embedded nanotubes. The hollow interiors of the plurality of embedded carbon nanotubes comprise the pores of the microporous membrane.

Description: The Mars Science Laboratory heatshield was designed to withstand a fully turbulent heat pulse based on test results and computational analysis on a pre-flight design trajectory. Instrumentation on the flight heatshield measured in-depth temperatures in the thermal protection system. The data indicate that boundary layer transition occurred at 5 of 7 thermocouple locations prior to peak heating. Data oscillations at 3 pressure measurement locations may also indicate transition. This paper presents the heatshield temperature and pressure data, possible explanations for the timing of boundary layer transition, and a qualitative comparison of reconstructed and computational heating on the as-flown trajectory. Boundary layer Reynolds numbers that are typically used to predict transition are compared to observed transition at various heatshield locations. A uniform smooth-wall transition Reynolds number does not explain the timing of boundary layer transition observed during flight. A roughness-based Reynolds number supports the possibility of transition due to discrete or distributed roughness elements on the heatshield. However, the distributed roughness height would have needed to be larger than the pre-flight assumption. The instrumentation confirmed the predicted location of maximum turbulent heat flux near the leeside shoulder. The reconstructed heat flux at that location is bounded by smooth-wall turbulent calculations on the reconstructed trajectory, indicating that augmentation due to surface roughness probably did not occur. Turbulent heating on the downstream side of the heatshield nose exceeded smooth-wall computations, indicating that roughness may have augmented heating. The stagnation region also experienced heating that exceeded computational levels, but shock layer radiation does not fully explain the differences.

Description: Current CFD modeling techniques are known to do a poor job of predicting the mixing rate and persistence of slot film flow in co-annular flowing ducts with relatively small velocity differences but large thermal gradients. A co-annular test was devised to empirically determine the mixing rate of slot film flow in a constant area circular duct (D approx. 1ft, L approx. 10ft). The axial rate of wall heat-up is a sensitive measure of the mixing rate of the two flows. The inflow conditions were varied to simulate a variety of conditions characteristic of moderate by-pass ratio engines. A series of air temperature measurements near the duct wall provided a straightforward means to measure the axial temperature distribution and thus infer the mixing rate. This data provides a characterization of the slot film mixing rates encountered in typical jet engine environments. The experimental geometry and entrance conditions, along with the sensitivity of the results as the entrance conditions vary, make this a good test for turbulence models in a regime important to modern air-breathing propulsion research and development.

Description: As part of the Nuclear Systems Office Fission Surface Power Technology Demonstration Unit (TDU) project, a reactor simulator test loop (RxSim) was design & built to perform integrated testing of the TDU components. In particular, the objectives of RxSim testing was to verify the operation of the core simulator, the instrumentation and control system, and the ground support gas and vacuum test equipment. In addition, it was decided to include a thermal test of a cold trap purification design and a pump performance test at pump voltages up to 150 V since the targeted mass flow rate of 1.75 kg/s was not obtained in the RxSim at the originally constrained voltage of 120 V. This paper summarizes RxSim testing. The gas and vacuum ground support test equipment performed effectively in NaK fill, loop pressurization, and NaK drain operations. The instrumentation and control system effectively controlled loop temperature and flow rates or pump voltage to targeted settings. The cold trap design was able to obtain the targeted cold temperature of 480 K. An outlet temperature of 636 K was obtained which was lower than the predicted 750 K but 156 K higher than the cold temperature indicating the design provided some heat regeneration. The annular linear induction pump (ALIP) tested was able to produce a maximum flow rate of 1.53 kg/s at 800 K when operated at 150 V and 53 Hz. Keywords: fission, space power, nuclear, liquid metal, NaK.

Description: A database compilation of hypersonic shock-wave/turbulent boundary layer experiments is provided. The experiments selected for the database are either 2D or axisymmetric, and include both compression corner and impinging type SWTBL interactions. The strength of the interactions range from attached to incipient separation to fully separated flows. The experiments were chosen based on criterion to ensure quality of the datasets, to be relevant to NASA's missions and to be useful for validation and uncertainty assessment of CFD Navier-Stokes predictive methods, both now and in the future. An emphasis on datasets selected was on surface pressures and surface heating throughout the interaction, but include some wall shear stress distributions and flowfield profiles. Included, for selected cases, are example CFD grids and setup information, along with surface pressure and wall heating results from simulations using current NASA real-gas Navier-Stokes codes by which future CFD investigators can compare and evaluate physics modeling improvements and validation and uncertainty assessments of future CFD code developments. The experimental database is presented tabulated in the Appendices describing each experiment. The database is also provided in computer-readable ASCII files located on a companion DVD.

Description: Three major projects were undertaken during the Summer 2013 USRP Internship: (A) assisting the cTAPS group with component and pressure vessel system analyses and documentation, (B) designing a hoisting fixture for a solid rocket motor, (C) finding an alternative to removing the DOT rated gaseous nitrogen tank from the roof for hydrostatic testing. Hypergolic Material Assessments (HMAs) and pressure calculations were performed on components of pressure systems. Additionally, component information was logged in the Standard Parts Database to provide a location where system designers can find information regarding components, including their specifications and compatibility with fluids. A hoisting fixture was designed to hoist a solid rocket motor and meets the specifications related to stress and size. However, there are issues with the fixtures bolt head allotment, the bolt spacing, and the complexity of the part. Finally, calculations were performed on an expiring DOT rated gaseous nitrogen tank in an attempt to re-rate it per ASME standards. This was unsuccessful so other options are being explored for the tank. While much progress was made on all three projects, there is still work to be performed on each project to achieve the desired results.

Description: Fluid-flow simulation over a computer-generated aircraft surface is generated using inviscid and viscous simulations. A fluid-flow mesh of fluid cells is obtained. At least one inviscid fluid property for the fluid cells is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. A set of intersecting fluid cells that intersects the aircraft surface are identified. One surface mesh polygon of the surface mesh is identified for each intersecting fluid cell. A boundary-layer prediction point for each identified surface mesh polygon is determined. At least one boundary-layer fluid property for each boundary-layer prediction point is determined using the at least one inviscid fluid property of the corresponding intersecting fluid cell and a boundary-layer simulation that simulates fluid viscous effects. At least one updated fluid property for at least one fluid cell is determined using the at least one boundary-layer fluid property and the inviscid fluid simulation.

Description: A fluid-flow simulation over a computer-generated aircraft surface is generated using a diffusion technique. The surface is comprised of a surface mesh of polygons. A boundary-layer fluid property is obtained for a subset of the polygons of the surface mesh. A pressure-gradient vector is determined for a selected polygon, the selected polygon belonging to the surface mesh but not one of the subset of polygons. A maximum and minimum diffusion rate is determined along directions determined using a pressure gradient vector corresponding to the selected polygon. A diffusion-path vector is defined between a point in the selected polygon and a neighboring point in a neighboring polygon. An updated fluid property is determined for the selected polygon using a variable diffusion rate, the variable diffusion rate based on the minimum diffusion rate, maximum diffusion rate, and angular difference between the diffusion-path vector and the pressure-gradient vector.

Description: Future manned missions to Mars are expected to greatly increase the space vehicle's size, weight, and heat dissipation requirements. An effective means to reducing both size and weight is to replace single-phase thermal management systems with two-phase counterparts that capitalize upon both latent and sensible heat of the coolant rather than sensible heat alone. This shift is expected to yield orders of magnitude enhancements in flow boiling and condensation heat transfer coefficients. A major challenge to this shift is a lack of reliable tools for accurate prediction of two-phase pressure drop and heat transfer coefficient in reduced gravity. Developing such tools will require a sophisticated experimental facility to enable investigators to perform both flow boiling and condensation experiments in microgravity in pursuit of reliable databases. This study will discuss the development of the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS), which was initiated in 2012 in collaboration between Purdue University and NASA Glenn Research Center. This facility was recently tested in parabolic flight to acquire condensation data for FC-72 in microgravity, aided by high-speed video analysis of interfacial structure of the condensation film. The condensation is achieved by rejecting heat to a counter flow of water, and experiments were performed at different mass velocities of FC-72 and water and different FC-72 inlet qualities. It is shown that the film flow varies from smooth-laminar to wavy-laminar and ultimately turbulent with increasing FC-72 mass velocity. The heat transfer coefficient is highest near the inlet of the condensation tube, where the film is thinnest, and decreases monotonically along the tube, except for high FC-72 mass velocities, where the heat transfer coefficient is enhanced downstream. This enhancement is attributed to both turbulence and increased interfacial waviness. One-ge correlations are shown to predict the average condensation heat transfer coefficient with varying degrees of success, and a recent correlation is identified for its superior predictive capability, evidenced by a mean absolute error of 21.7%.

Description: An aeroheating measurement test campaign of an Apollo capsule model with laminar and turbulent boundary layer was performed in the free-piston shock tunnel HIEST at JAXA Kakuda Space Center. A 250mm-diameter 6.4%-scaled Apollo CM capsule model made of SUS-304 stainless steel was applied in this study. To measure heat flux distribution, the model was equipped with 88 miniature co-axial Chromel-Constantan thermocouples on the heat shield surface of the model. In order to promote boundary layer transition, a boundary layer trip insert with 13 "pizza-box" isolated roughness elements, which have 1.27mm square, were placed at 17mm below of the model geometric center. Three boundary layer trip inserts with roughness height of k=0.3mm, 0.6mm and 0.8mm were used to identify the appropriate height to induce transition. Heat flux records with or without roughness elements were obtained for model angles of attack 28 under stagnation enthalpy between H(sub 0)=3.5MJ/kg to 21MJ/kg and stagnation pressure between P(sub 0)=14MPa to 60MPa. Under the condition above, Reynolds number based on the model diameter was varied from 0.2 to 1.3 million. With roughness elements, boundary layer became fully turbulent less than H(sub 0)=9MJ/kg condition. However, boundary layer was still laminar over H(sub 0)=13MJ/kg condition even with the highest roughness elements. An additional experiment was also performed to correct unexpected heat flux augmentation observed over H(sub 0)=9MJ/kg condition.

Description: We report the progress in the development of the chilldown model for rapid cryogenic loading system developed at KSC. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The two-phase flow model of the chilldown is approximated as one-dimensional homogeneous fluid flow with no slip condition for the interphase velocity. The model is built using commercial SINDAFLUINT software. The results of numerical predictions are in good agreement with the experimental time traces. The obtained results pave the way to the application of the SINDAFLUINT model as a verification tool for the design and algorithm development required for autonomous loading operation.

Description: Flame spread tests have been conducted over polymethylmethacrylate (PMMA) samples in San Diego State University's Narrow Channel Apparatus (SDSU NCA). The Narrow Channel Apparatus (NCA) has the ability to suppress buoyant flow in horizontally spreading flames, and is currently being investigated as a possible replacement or complement to NASA's current material flammability test standard for non-metallic solids, NASA-STD-(I)-6001B Test 1. The buoyant suppression achieved with a NCA allows for tests to be conducted in a simulated microgravity atmosphere-a characteristic that Test 1 lacks since flames present in Test 1 are buoyantly driven. The SDSU NCA allows for flame spread tests to be conducted with varying opposed flow oxidizer velocities, oxygen percent by volume, and total pressure. Also, since the test sample is placed symmetrically between two confining plates so that there is a gap above and below the sample, this gap can be adjusted. This gap height adjustment allows for a compromise between heat loss from the flame to the confining boundaries and buoyancy suppression achieved by those boundaries. This article explores the effect gap height has on the flame spread rate for 75 m thick PMMA at 1 atm pressure and 21% oxygen concentration by volume in the SDSU NCA. Flame spread results from the SDSU NCA for thin cellulose fuels have previously been compared to results from tests in actual microgravity at various test conditions with the same sample materials and were found to be in good agreement. This article also presents results from the SDSU NCA for PMMA at 1 atm pressure, opposed oxidizer velocity ranging from 3 to 35 cm/s, oxygen concentration by volume at 21%, 30 %, and 50% and fuel thicknesses of 50 and 75 m. These results are compared to results obtained in actual microgravity for PMMA obtained at the 4.5s drop tower of MGLAB in Gifu, Japan, and the 5.2s drop tower at NASA's Zero-Gravity Research Facility in Cleveland, OH. This comparison confirms that at 1 atm pressure, the SDSU NCA successfully simulates microgravity for not only thin cellulose fuels, but also for thin PMMA sheets as well. This further supports the idea that the NCA is a viable option to complement or replace NASA's Test 1 for material flammability testing. Tests with thick fuels will be conducted in the future to further characterize the SDSU NCA.

Description: This paper describes the development and testing of a scalable thermal management architecture for instruments, subsystems, or systems that must operate in severe space environments with wide variations in sink temperature. The architecture involves a serial linkage of one or more hot-side variable conductance heat pipes (VCHPs) to one or more cold-side loop heat pipes (LHPs). The VCHPs provide wide area heat acquisition, limited distance thermal transport, modest against gravity pumping, concentrated LHP startup heating, and high switching ratio variable conductance operation. The LHPs provide localized heat acquisition, long distance thermal transport, significant against gravity pumping, and high switching ratio variable conductance operation. The single-VCHP, single-LHP system described herein was developed to maintain thermal control of a small robotic lunar lander throughout the lunar day-night thermal cycle. It is also applicable to other variable heat rejection space missions in severe environments. Operationally, despite a 60-70% gas blocked VCHP condenser during ON testing, the system was still able to provide 2-4 W/K ON conductance, 0.01 W/K OFF conductance, and an end-to-end switching ratio of 200-400. The paper provides a detailed analysis of VCHP condenser performance, which quantified the gas blockage situation. Future multi-VCHP/multi-LHP thermal management system concepts that provide power/transport length scalability are also discussed.

Description: An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6- DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Description: This study presents a detailed analysis of a sweeping jet (fluidic oscillator) actuator. The sweeping jet actuator promises to be a viable flow control actuator candidate due to its simple, no moving part structure and its high momentum, spatially oscillating flow output. Hot-wire anemometer and particle image velocimetry measurements were carried out with an emphasis on understanding the actuator flow field in a quiescent environment. The time averaged, fluctuating, and instantaneous velocity measurements are provided. A modified actuator concept that incorporates high-speed solenoid valves to control the frequency of oscillation enabled phase averaged measurements of the oscillating jet. These measurements reveal that in a given oscillation cycle, the oscillating jet spends more time on each of the Coanda surfaces. In addition, the modified actuator generates four different types of flow fields, namely: a non oscillating downward jet, a non oscillating upward jet, a non oscillating straight jet, and an oscillating jet. The switching from an upward jet to a downward jet is accomplished by providing a single pulse from the solenoid valve. Once the flow is switched, the flow stays there until another pulse is received. The oscillating jet is compared with a non oscillating straight jet, which is a typical planar turbulent jet. The results indicate that the oscillating jet has a higher (5 times) spreading rate, more flow entrainment, and higher velocity fluctuations (equal to the mean velocity).

Description: Direct numerical simulation is used to investigate the transition induced by threedimensional isolated roughness elements in a supersonic boundary layer at a free stream Mach number of 3.5. Simulations are performed for two different configurations: one is a square planform roughness and the other is a diamond planform roughness. The mean-flow calculations show that the roughness induces counter rotating streamwise vortices downstream of the roughness. These vortices persist for a long distance downstream and lift the low momentum fluid from the near wall region and place it near the outer part of the boundary layer. This forms highly inflectional boundary layer profiles. These observations agree with recent experimental observations. The receptivity calculations showed that the amplitudes of the mass-flux fluctuations near the neutral point for the diamond shape roughness are the same as the amplitude of the acoustic disturbances. They are three times smaller for the square shape roughness.

Description: For thermal protection system (heat shield) design for space vehicle entry into earth and other planetary atmospheres, it is essential to know the augmentation of the heat flux due to vehicle surface roughness. At the NASA Ames Hypervelocity Free Flight Aerodynamic Facility (HFFAF) ballistic range, a campaign of heat flux studies on rough models, using infrared camera techniques, has been initiated. Several phenomena can interfere with obtaining good heat flux data when using this measuring technique. These include leakage of the hot drive gas in the gun barrel through joints in the sabot (model carrier) to create spurious thermal imprints on the model forebody, deposition of sabot material on the model forebody, thereby changing the thermal properties of the model surface and unknown in-barrel heating of the model. This report presents developments in launch techniques to greatly reduce or eliminate these problems. The techniques include the use of obturator cups behind the launch package, enclosed versus open front sabot designs and the use of hydrogen gas in the launch tube. Attention also had to be paid to the problem of the obturator drafting behind the model and impacting the model. Of the techniques presented, the obturator cups and hydrogen in the launch tube were successful when properly implemented

Description: There have been few discussions on using Computational Fluid Dynamics (CFD) without experimental validation. Pairing experimental data, uncertainty analysis, and analytical predictions provides a comprehensive approach to verification and is the current state of the art. With pressed budgets, collecting experimental data is rare or non-existent. This paper investigates and proposes a method to perform CFD uncertainty analysis only from computational data. The method uses current CFD uncertainty techniques coupled with the Student-T distribution to predict the heat transfer coefficient over a at plate. The inputs to the CFD model are varied from a specified tolerance or bias error and the difference in the results are used to estimate the uncertainty. The variation in each input is ranked from least to greatest to determine the order of importance. The results are compared to heat transfer correlations and conclusions drawn about the feasibility of using CFD without experimental data. The results provide a tactic to analytically estimate the uncertainty in a CFD model when experimental data is unavailable

Description: Several analytical mechanical slosh models for a cylindrical tank with flat bottom are reviewed. Even though spacecrafts use cylinder shaped tanks, most of those tanks usually have elliptical domes. To extend the application of the analytical models for a cylindrical tank with elliptical domes, the modified slosh parameter models are proposed in this report by mapping an elliptical dome cylindrical tank to a flat top/bottom cylindrical tank while maintaining the equivalent liquid volume. For the low Bond number case, the low-g slosh models were also studied. Those low-g models can be used for Bond number 〉 10. The current low-g slosh models were also modified to extend their applications for the case that liquid height is smaller than the tank radius. All modified slosh models are implemented in MATLAB m-functions and are collected in the developed MST (Mechanical Slosh Toolbox).

Description: The application of a circulation control system for high Reynolds numbers was experimentally validated with the Fundamental Aerodynamic Subsonic Transonic Modular Active Control semi-span model in the NASA Langley National Transonic Facility. This model utilized four independent flow paths to modify the lift and thrust performance of a representative advanced transport type of wing. The design of the internal flow paths highlights the challenges associated with high Reynolds number testing in a cryogenic pressurized wind tunnel. Weight flow boundaries for the air delivery system were identified at mildly cryogenic conditions ranging from 0.1 to 10 lbm/sec. Results from the test verified system performance and identified solutions associated with the weight-flow metering system that are linked to internal perforated plates used to achieve flow uniformity at the jet exit.

Description: The NASA Glenn Research Center (GRC) has been testing 100 We class, free-piston Stirling convertors for potential use in Stirling Radioisotope Power Systems (RPS) for space science and exploration missions. Free-piston Stirling convertors are capable of achieving a 38% conversion efficiency, making Stirling attractive for meeting future power system needs in light of the shrinking U.S. plutonium fuel supply. Convertors currently on test include four Stirling Technology Demonstration Convertors (TDCs), manufactured by the Stirling Technology Company (STC), and six Advanced Stirling Convertors (ASCs), manufactured by Sunpower, Inc. Total hours of operation is greater than 514,000 hours (59 years). Several tests have been initiated to demonstrate the functionality of Stirling convertors for space applications, including: in-air extended operation, thermal vacuum extended operation. Other tests have also been conducted to characterize Stirling performance in anticipated mission scenarios. Data collected during testing has been used to support life and reliability estimates, drive design changes and improve quality, and plan for expected mission scenarios. This paper will provide a summary of convertors tested at NASA GRC and discuss lessons learned through extended testing.

Description: The verification and validation of the Launch Ascent and Vehicle Aerodynamics (LAVA) computational fluid dynamics (CFD) solver is presented. A modern strategy for verification and validation is described incorporating verification tests, validation benchmarks, continuous integration and version control methods for automated testing in a collaborative development environment. The purpose of the approach is to integrate the verification and validation process into the development of the solver and improve productivity. This paper uses the Method of Manufactured Solutions (MMS) for the verification of 2D Euler equations, 3D Navier-Stokes equations as well as turbulence models. A method for systematic refinement of unstructured grids is also presented. Verification using inviscid vortex propagation and flow over a flat plate is highlighted. Simulation results using laminar and turbulent flow past a NACA 0012 airfoil and ONERA M6 wing are validated against experimental and numerical data.

Description: The effects of fin leading-edge radius and sweep angle on peak heating rates due to shock-shock interactions were investigated in the NASA Langley Research Center 20-inch Mach 6 Air Tunnel. The fin model leading edges, which represent cylindrical leading edges or struts on hypersonic vehicles, were varied from 0.25 inches to 0.75 inches in radius. A 9deg wedge generated a planar oblique shock at 16.7deg to the flow that intersected the fin bow shock, producing a shock-shock interaction that impinged on the fin leading edge. The fin angle of attack was varied from 0deg (normal to the free-stream) to 15deg and 25deg swept forward. Global temperature data was obtained from the surface of the fused silica fins using phosphor thermography. Metal oil flow models with the same geometries as the fused silica models were used to visualize the streamline patterns for each angle of attack. High-speed zoom-schlieren videos were recorded to show the features and temporal unsteadiness of the shock-shock interactions. The temperature data were analyzed using one-dimensional semi-infinite as well as one- and two-dimensional finite-volume methods to determine the proper heat transfer analysis approach to minimize errors from lateral heat conduction due to the presence of strong surface temperature gradients induced by the shock interactions. The general trends in the leading-edge heat transfer behavior were similar for the three shock-shock interactions, respectively, between the test articles with varying leading-edge radius. The dimensional peak heat transfer coefficient augmentation increased with decreasing leading-edge radius. The dimensional peak heat transfer output from the two-dimensional code was about 20% higher than the value from a standard, semi-infinite onedimensional method.

Description: The thermal efficiency of a SparkJet is evaluated by measuring the impulse response of a pendulum subject to a single spark discharge. The SparkJet is attached to the end of a pendulum. A laser displacement sensor is used to measure the displacement of the pendulum upon discharge. The pendulum motion is a function of the fraction of the discharge energy that is channeled into the heating of the gas (i.e., increasing the translational-rotational temperature). A theoretical perfect gas model is used to estimate the portion of the energy from the heated gas that results in equivalent pendulum displacement as in the experiment. The earlier results from multiple runs for different capacitances of C = 3, 5, 10, 20, and 40(micro)F demonstrate that the thermal efficiency decreases with higher capacitive discharges.1 In the current paper, results from additional run cases have been included and confirm the previous results

Description: A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Description: Preinjection of fuel on the forebody of an airbreathing vehicle is a proposed method to gain access to hypervelocity flight Mach numbers. However, this creates the possibility of autoignition either near the wall or in the core of the flow, thereby consuming fuel prematurely as well as increasing the amount of pressure drag on the vehicle. The computational fluid dynamics code VULCAN was used to conduct three dimensional simulations of the reacting flow in the vicinity of hydrogen injectors on a flat plate at conditions relevant to a Mach 12 notional flight vehicle forebody to determine the location where autoignition occurs. Active wall cooling strategies were formulated and simulated in response to regions of autoignition. It was found that tangential film cooling using hydrogen or helium were both able to nearly or completely eliminate wall autoignition in the flow domain of interest.

Description: Propellants flow through specialized mechanical hardware that is designed for effective and safe ignition and sustained combustion of the propellants. By integrating a micro-fluidic porous media element between a propellant feed source and the combustion chamber, an effective and reliable propellant injector head may be implemented that is capable of withstanding transient combustion and detonation waves that commonly occur during an ignition event. The micro-fluidic porous media element is of specified porosity or porosity gradient selected to be appropriate for a given propellant. Additionally the propellant injector head design integrates a spark ignition mechanism that withstands extremely hot running conditions without noticeable spark mechanism degradation.

Description: A variety of increasingly intricate container geometries are under consideration for the passive manipulation of liquids aboard spacecraft where the impact of gravity may be neglected. In this study we examine the mean curvature of a liquid volume confined within a radial array of disconnected vanes of infinite extent. This particular geometry possesses a number of desirable characteristics relevant to waste water treatment aboard spacecraft for life support. It is observed that under certain conditions the slender shape of the free surface approaches an asymptote, which can be predicted analytically using new hybrid boundary conditions proposed herein. This contribution represents possibly the final extension of what has been referred to as the method of de Lazzer et al. (1996). The method enables the integration of the Young-Laplace equation over a domain with its boundaries, including the wetted portion of the solid boundaries, symmetry planes, and circular arcs representing free surfaces at the center plane of the liquid body. Asymptotic solutions at several limits are obtained and the analysis is confirmed with numerical computations.

Description: Since 1960, the Arc Jet Complex at NASA Ames Research Center has been a source of long-duration, high-enthalpy flow for materials testing with application to the thermal protection of aerospace vehicle components. From their inception the facilities have played an integral role supporting many of NASA's space flight programs and numerous DoD projects. In recent years advancements in computational fluid dynamics (CFO) have made the resultant models a valuable tool for assessing and predicting performance, however, the inherent limitation of models to compensate for a dissociated, transitionally high temperature nonequilibrated flowfield have made further measurements necessary. The use of electrostatic probe diagnostics within similarly harsh plasma environments in previous studies have been met with much success. In this study, the use of a single Langmuir probe was implemented to characterize the plasma parameters of interest as they vary radially within a large volume of the plume. Classical Langmuir probe theory was applied to achieve first order estimates of the heavy particle temperature, the ratio T(sub i)/T(sub e), and the ionization fraction. As expected, both the electron temperature and electron density measurements show a dependence on radial distance from the plume centerline, with electron density profiles showing the largest dependence. This paper aims to validate and strengthen modeling work conducted at the NASA Ames Research Center by measuring the critical plasma parameters (electron temperature, electron density, and plasma potential) within the arc plume of a subscale arc jet. These parameters are intended to give physical insight into the flow characteristics while providing the necessary boundary conditions to validate full scale simulations.

Description: The development of computational fluid dynamics over the last few decades has yielded enormous successes and capabilities being routinely employed today; however there remain some open problems to be properly resolved-some are fundamental in nature and some resolvable by operational changes. These two categories are distinguished and broadly explored previously. One, that belongs to the former, is the so-called overheating problem, especially in rarefying flow. This problem up to date still dogs every method known to the author; a solution to it remains elusive. The study in this paper concludes that: (1) the entropy increase is quantitatively linked to the increase in the temperature increase, (2) it is argued that the overheating is inevitable in the current shock capturing or traditional finite difference framework, and (3) a simple hybrid method is proposed that removes the overheating problem in the rarefying problems, but also retains the property of accurate shock capturing. This remedy (enhancement of current numerical methods) can be included easily in the present Eulerian codes.

Description: Ultrasonic probes were designed, fabricated and tested for high temperature health monitoring system. The goal of this work was to develop the health monitoring system that can determine the height level of the condensed water through the pipe wall at high temperature up to 250 deg while accounting for the effects of surface perturbation. Among different ultrasonic probe designs, 2.25 MHz probes with air backed configuration provide satisfactory results in terms of sensitivity, receiving reflections from the target through the pipe wall. A series of tests were performed using the air-backed probes under irregular conditions, such as surface perturbation and surface disturbance at elevated temperature, to qualify the developed ultrasonic system. The results demonstrate that the fabricated air-backed probes combined with advanced signal processing techniques offer the capability of health monitoring of steam pipe under various operating conditions.

Description: No abstract available

Description: NASA's Environmentally Responsible Aviation (ERA) Project is working with industry to develop the fuel flexible combustor technologies for a new generation of low-emissions engine targeted for the 2020 timeframe. These new combustors will reduce nitrogen oxide (NOx) emissions to half of current state-of-the-art (SOA) combustors, while simultaneously reducing noise and fuel burn. The purpose of the low NOx fuel-flexible combustor research is to advance the Technology Readiness Level (TRL) and Integration Readiness Level (IRL) of a low NOx, fuel flexible combustor to the point where it can be integrated in the next generation of aircraft. To reduce project risk and optimize research benefit NASA chose to found two Phase 1 contracts. The first Phase 1 contracts went to engine manufactures and were awarded to: General Electric Company, and Pratt & Whitney Company. The second Phase 1 contracts went to fuel injector manufactures Goodrich Corporation, Parker Hannifin Corporation, and Woodward Fuel System Technology. In 2012, two sector combustors were tested at NASA's ASCR. The results indicated 75% NOx emission reduction below the 2004 CAEP/6 regulation level.