Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

Computation of Anisotropic Bi-Material Interfacial Fracture Parameters and Delamination Creteria

This report documents the recent developments in methodologies for the evaluation of the integrity and durability of composite structures, including i) the establishment of a stress-intensity-factor based fracture criterion for bimaterial interfacial cracks in anisotropic materials (see Sec. 2); ii) the development of a virtual crack closure integral method for the evaluation of the mixed-mode stress intensity factors for a bimaterial interfacial crack (see Sec. 3). Analytical and numerical results show that the proposed fracture criterion is a better fracture criterion than the total energy release rate criterion in the characterization of the bimaterial interfacial cracks. The proposed virtual crack closure integral method is an efficient and accurate numerical method for the evaluation of mixed-mode stress intensity factors.

Chow, W-T.↗

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated-, the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered: one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of' the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows. pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow.

So, R. M. C.↗

Modeling the Dynamics of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distortion theory (RDT). The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined. Key words: transition/turbulence modeling, homogeneous shear, rapid distortion theory.

Thacker, W. D.↗

Airborne Sunphotometry of African Dust and Marine Boundary Layer Aerosols in PRIDE

The Puerto Rico Dust Experiment (PRIDE) was conducted during summer 2000 to study the radiative, microphysical and transport properties of Saharan dust in the Caribbean region. During PRIDE, NASA Ames Research Center's six-channel airborne autotracking sunphotometer (AATS-6) was operated aboard a Piper Navajo airplane based at Roosevelt Roads Naval Station on the northeast coast of Puerto Rico. AATS-6 measurements were taken during 21 science flights off the coast of Puerto Rico in the western Caribbean. Data were acquired within and above the Marine Boundary Layer (MBL) and the Saharan Aerosol Layer (SAL) up to 5.5 km altitude tinder a wide range of dust loadings. Aerosol optical depth (AOD) spectra and columnar water vapor (CWV) values have been calculated from the AATS-6 measurements by using sunphotometer calibration data obtained at Mauna Loa Observatory (3A kin ASL) before (May) and after (October) PRIDE. Mid-visible AOD values measured near the surface during PRIDE ranged from 0.07 on the cleanest day to 0.55 on the most turbid day. Values measured above the MBL were as high as 0.35; values above the SAL were as low as 0.01. The fraction of total column AOD due to Saharan dust cannot be determined precisely from AATS-6 AOD data alone due to the uncertainty in the extent of vertical mixing of the dust down through the MBL. However, analyses of ground-based and airborne in-situ aerosol sampling measurements and ground-based aerosol lidar backscatter data should yield accurate characterization of the vertical mixing that will enable calculation of the Saharan dust AOD component from the sunphotometer data. Examples will be presented showing measured AATS-6 AOD spectra, calculated aerosol extinction and water vapor density vertical profiles, and aerosol size distributions retrieved by inversion of the AOD spectra. Near sea-surface AOD spectra acquired by AATS-6 during horizontal flight legs at 30 m ASL are available for validation of AOD derived from coincident satellite sensor (TOMS, MODIS, MISR) measurements. AATS-6 AOD data acquired during numerous aircraft ascents and descents through the MBL and the SAL should permit atmospheric column closure analyses with respect to aerosol optical depth, extinction, and size distribution by comparison with coincident aircraft-based in-situ particle size distribution measurements and ground-based (Cabras Island, Puerto Rico) micropulse lidar aerosol backscatter measurements. The aerosol information derived from the column closure analyses can be used subsequently to calculate radiative flux changes, which can then be compared with coincident spectral flux measurements taken from aboard the aircraft with a solar spectral flux radiometer.

Livingston, John M.↗

Man-Made Closed Ecological Systems as Model of Natural Ecosystems and as Means to Provide High Quality of Human Life in Adverse Environment

For its more than thirty year long history, the experimental creation of closed ecological systems has from its very sources been distinctly and strongly motivated by the development of human life-support systems for space. As the trend developed its fundamental significance and broad opportunities of terrestrial applications of the technologies under development were coming to the foreground. Nowadays, it can be argued that development of closed ecosystems is experimental foundation of a new branch of ecology biospherics, the goal of which is to comprehend the regularities of existence of the biosphere as a unique in the Universe (in that part of it that we know, at least) closed ecosystem. Closed technologies can be implemented in life-support systems under adverse conditions of life on the Earth - in Arctic and Antarctic latitudes, deserts, high mountains or deep in the ocean, as well as under the conditions of polluted water and air. In space where the environment is hostile for life all around the cell of life should be sealed and the life-support system as close to the ideally closed cyclic turnover of the matter as possible. Under terrestrial conditions designers should strive for maximum closure of the limiting factor: water - in deserts, oxygen - in high mountains, energy - in polar latitudes, etc. Essential closure of a life-support systems withstands also pollution of the environment by the wastes of human vital activity. This is of particular importance for the quarantine of visited planets, and on the Earth under the conditions of deficient heat in high latitudes and water in and areas. The report describes experimental ecosystem 'BIOS' and exohabitats being designed on its basis, which are adapted to various conditions, described capacities of the Center for Closed Ecosystems in Drasnoyarsk for international collaboration in research and education in this field.

Gitelson, I. I.↗

Cloud-Scale Numerical Modeling of the Arctic Boundary Layer

The research objective of this NASA grant-funded project was to determine in detail how large-scale processes. in combination with cloud-scale radiative, microphysical, and dynamical processes, govern the formation and multi-layered structure of Arctic stratus clouds. This information will be useful for developing and improving 1D (one dimensional) boundary layer models for the Arctic. Also, to quantitatively determine the effects of leads on the large-scale budgets of sensible heat, water vapor, and condensate in a variety of Arctic winter conditions. This information will be used to identify the most important lead-flux processes that require parameterization in climate models. Our approach was to use a high-resolution numerical model, the 2D (two dimensional) University of Utah Cloud Resolving Model (UU CRM), and its 1D version, the University of Utah Turbulence Closure Model (UU TCM), a boundary layer model based on third-moment turbulence closure, as well as a large-eddy simulation (LES) model originally developed by C.H. Moeng.

Kruegen, Steven K.↗

Shipboard Sunphotometer Measurements of Aerosol Optical Depth Spectra and Columnar Water Vapor During ACE-2, and Comparison with Selected Land, Ship, Aircraft, and Satellite Measurements

Analyses of aerosol optical depth (AOD) and columnar water vapor (CWV) measurements acquired with NASA Ames Research Center's six-channel Airborne Tracking Sunphotometer (AATS-6) operated aboard the R/V (research vehicle) Professor Vodyanitskiy during the second Aerosol Characterization Experiment (ACE-2) are discussed. Data are compared with various in situ and remote measurements for selected cases. The focus is on 10 July, when the Pelican airplane flew within 70 km of the ship near the time of a NOAA (National Oceanographic and Atmospheric Administration)-14/AVHRR (Advanced Very High Resolution Radiometer) satellite overpass and AOD measurements with the 14-channel Ames Airborne Tracking Sunphotometer (AATS-14) above the marine boundary layer (MBL) permitted calculation of AOD within the MBL from the AATS-6 measurements. A detailed column closure test is performed for MBL AOD on 10 July by comparing the AATS-6 MBL AODs with corresponding values calculated by combining shipboard particle size distribution measurements with models of hygroscopic growth and radiosonde humidity profiles (plus assumptions on the vertical profile of the dry particle size distribution and composition). Large differences (30-80% in the mid-visible) between measured and reconstructed AODs are obtained, in large part because of the high sensitivity of the closure methodology to hygroscopic growth models, which vary considerably and have not been validated over the necessary range of particle size/composition distributions. The wavelength dependence of AATS-6 AODs is compared with the corresponding dependence of aerosol extinction calculated from shipboard measurements of aerosol size distribution and of total scattering measured by a shipboard integrating nephelometer for several days. Results are highly variable, illustrating further the great difficulty of deriving column values from point measurements. AATS-6 CWV values are shown to agree well with corresponding values derived from radiosonde measurements during eight soundings on seven days and also with values calculated from measurements taken on 10 July with the AATS-14 and the University of Washington Passive Humidigraph aboard the Pelican.

Livingston, John M.↗

Gravitational Effects on Distortion in Sintering

During sintering a powder compact gains strength through low-temperature interparticle bonding, usually induced by solid-state surface diffusion, followed by further strength contributions from high-temperature densification. In cases where a liquid phase forms, sintering densification is accelerated and shape retention is sustained while open pores remain and contribute capillary forces. Unfortunately, sintering densification requires the compact become thermally softened to a point where creep strain rates reach levels near 10(exp -2)/s when the liquid forms. On the other hand, thermal softening of the powder compact substantially reduces the strength at high temperatures. Therefore, the in situ strength evolution during sintering is a primary focus to separate compact densification (as required for high performance) with minimized distortion (as required for net-shaping). With respect to gravitation effects on distortion during sintering there are two points of substantial weakness - prior to significant interparticle bonding and during final pore closure. This research is focused on understanding the competition among interparticle neck growth, densification, thermal softening, grain boundary wetting, capillary effects associated with liquid wetting and residual porosity, and gravity. Most surprising is the apparent role of gravity, where the deviatoric stress acting on the powder structure induces skeletal formation that reduces distortion. In contrast with theory, microgravity samples exhibit more distortion yet fail to fully densify. Results are presented on the experimental concepts supporting an emerging model of sintering strength evolution that enables understanding of both distortion and densification. The experiments have relied on tungsten heavy alloys, various combinations of dihedral angle, pore size, initial porosity, liquid:solid ratio, and heating rates. On Earth, the dominant factor with respect to distortion is the starting body heterogeneity. Current modeling efforts are seeking some means to uniformly predict the distortion based on a starting pore structure heterogeneity parameter. Densification is largely unaffected by the initial pore structure, but distortion is highly variable, suggesting that nonuniform pore closure might be a significant parameter during sintering. With respect to flight experiments, plans include removal of the solid body forces acting on the solid grains, allowing stabilization of the pore structure and examination of the buoyancy effects with regard to distortion. In microgravity there is the surprising result that compacts do not densify, yet distort a factor that is contrary to all current sintering models. Densification without distortion during liquid phase sintering was achieved by manipulating microstructure and its evolution during sintering. Microstructure parameters such as the solid volume fraction, dihedral angle, initial porosity, and pore size were varied to measure densification and distortion behavior during LPS using W-Ni-Cu alloys. Green compacts were formed using ethylene-bis-stearamide as a pore-forming agent with the amount of polymer controlling the initial porosity. Different initial pore sizes were generated by varying the polymer particle size. Dihedral angle was varied by changing the Ni:Cu ratio in the alloys. Finally, the solid volume fraction was adjusted via the tungsten content. Distortion was quantified using profiles determined with a coordinate measuring machine to calculate a distortion parameter. Sintering results showed that solid volume fraction and dihedral angle are the dominant factors on densification and distortion during liquid phase sintering. Distortion decreases with increasing solid volume fraction and dihedral angle, while initial porosity and pore size have no observable effect on distortion at nearly full densification. Various strategies emerge to improve distortion control in liquid phase sintering.

German, Randall M.↗

Performance of advanced trauma life support procedures in microgravity

BACKGROUND: Medical operations on the International Space Station will emphasize the stabilization and transport of critically injured personnel and so will need to be capable of advanced trauma life support (ATLS). METHODS: We evaluated the ATLS invasive procedures in the microgravity environment of parabolic flight using a porcine animal model. Included in the procedures evaluated were artificial ventilation, intravenous infusion, laceration closure, tracheostomy, Foley catheter drainage, chest tube insertion, peritoneal lavage, and the use of telemedicine methods for procedural direction. RESULTS: Artificial ventilation was performed and appeared to be unaltered from the 1-G environment. Intravenous infusion, laceration closure, percutaneous dilational tracheostomy, and Foley catheter drainage were achieved without difficulty. Chest tube insertion and drainage were performed with no more difficulty than in the 1-G environment due to the ability to restrain patient, operator and supplies. A Heimlich valve and Sorenson drainage system were both used to provide for chest tube drainage collection with minimal equipment, without the risk of atmospheric contamination, and with the capability to auto-transfuse blood drained from a hemothorax. The use of telemedicine in chest tube insertion was demonstrated to be useful and feasible. Peritoneal lavage using a percutaneous technique, although requiring less training to perform, was found to be dangerous in weightlessness due to the additional pressure of the bowel on the anterior abdominal wall creating a high risk of bowel perforation. CONCLUSIONS: The performance of ATLS procedures in microgravity appears to be feasible with the exception of diagnostic peritoneal lavage. Minor modifications to equipment and techniques are required in microgravity to effect surgical drainage in the presence of altered fluid dynamics, to prevent atmospheric contamination, and to provide for the restraint requirements. A parabolic simulation system was developed for equipment and procedure verification, physiological research, and possible crew medical officer training in the future.

Parabolic Flight↗

Prediction of Crack Growth under Variable-Amplitude Loading in Thin-Sheet 2024-T3 Aluminum Alloys

The present paper is concerned with the application of a "plasticity-induced" crack closure model to study fatigue crack growth under various load histories. The model was based on the Dugdale model but modified to leave plastically deformed material in the wake of the advancing crack. The model was used to correlate crack growth rates under constant-amplitude loading and then used to predict crack growth under variable-amplitude and spectrum loading on thin-sheet 2024- T3 aluminum alloys. Predicted crack-opening stresses agreed well with test data from the literature. The crack-growth lives agreed within a factor of two for single and repeated spike overloads/underloads and within 20 percent for spectrum loading. Differences were attributed to fretting-product-debris-induced closure and three-dimensional affects not included in the model.

Newman, J. C., Jr.↗

Cardiogenic oscillation phase relationships during single-breath tests performed in microgravity

We studied the phase relationships of the cardiogenic oscillations in the phase III portion of single-breath washouts (SBW) in normal gravity (1 G) and in sustained microgravity (microG). The SBW consisted of a vital capacity inspiration of 5% He-1.25% sulfurhexafluoride-balance O2, preceded at residual volume by a 150-ml Ar bolus. Pairs of gas signals, all of which still showed cardiogenic oscillations, were cross-correlated, and their phase difference was expressed as an angle. Phase relationships between inspired gases (e.g., He) and resident gas (n2) showed no change from 1 G (211 +/- 9 degrees) to microG (163 +/- 7 degrees). Ar bolus and He were unaltered between 1 G (173 +/- 15 degrees) and microG (211 +/- 25 degrees), showing that airway closure in microG remains in regions of high specific ventilation and suggesting that airway closure results from lung regions reaching low regional volume near residual volume. In contrast, CO2 reversed phase with He between 1 G (332 +/- 6 degrees) and microG (263 +/- 27 degrees), strongly suggesting that, in microG, areas of high ventilation are associated with high ventilation-perfusion ratio (VA/Q). This widening of the range of VA/Q in microG may explain previous measurements (G.K. Prisk, A.R. Elliott, H.J.B. Guy, J.M. Kosonen, and J.B. West J. Appl. Physiol. 79: 1290-1298, 1995) of an overall unaltered range of VA/Q in microG, despite more homogeneous distributions of both ventilation and perfusion.

Non-NASA Center↗

Mapping and Monitoring Delmarva Fox Squirrel Habitat Using an Airborne LiDAR Profiler

Twenty five hundred thirty nine kilometers of airborne laser profiling and videography data were acquired over the state of Delaware during the summer of 2000. The laser ranging measurements and video from approximately one-half of that data set (1304 km) were analyzed to identify and locate forested sites that might potentially support populations of Delmarva fox squirrel (DFS, Sciurus niger cinereus). The DFS is an endangered species previously endemic to tall, dense, mature forests with open understories on the Eastern Shore of the Chesapeake Bay. The airborne LiDAR employed in this study can measure forest canopy height and canopy closure, but cannot measure or infer understory canopy conditions. Hence the LiDAR must be viewed as a tool to map potential, not actual, habitat. Fifty-three potentially suitable DFS sites were identified in the 1304 km of flight transect data. Each of the 53 sites met the following criteria according to the LiDAR and video record: (1 ) at least 120m of contiguous forest; (2) an average canopy height greater than 20m; (3) an average canopy closure of >80%; and (4) no roofs, impervious surface (e.g., asphalt, concrete), and/or open water anywhere along the 120m length of the laser segment. Thirty-two of the 53 sites were visited on the ground and measurements taken for a DFS habitat suitability model. Seventy eight percent of the sites (25 of 32) were judged by the model to be suited to supporting a DFS population. Twenty-eight of the 32 sites visited in the field were in forest cover types (hardwood, mixed wood, conifer, wetlands) according to a land cover GIS map. Of these, 23 (82%) were suited to support DFS. The remaining 4 sites were located in nonforest cover types - agricultural or residential areas. Two of the four, or 50% were suited to the DFS. All of the LiDAR flight data, 2539 km, were analyzed to

Nelson, Ross↗

Compressor Stall Recovery Through Tip Injection Assessed

Aerodynamic stability is a fundamental limit in the compressor design process. The development of robust techniques for increasing stability has several benefits: enabling higher loading and fewer blades, increasing safety throughout a mission, increasing tolerance to stage mismatch during part-speed operation and speed transients, and providing an opportunity to match stages at the compressor maximum efficiency point, thus reducing fuel burn. Mass injection upstream of the tip of a high-speed axial compressor rotor is a stability enhancement approach known to be effective in suppressing stall in tip-critical rotors if the injection is activated before stall occurs. This approach to stall suppression requires that a reliable stall warning system be available. Tests have recently been performed to assess whether steady injection can also be used to recover from fully developed stall. If mass injection is effective in recovering from stall quickly enough to avoid structural damage or loss of engine power, then a stall warning system may not be required. The stall recovery tests were performed on a transonic compressor rotor at its design tip speed of 1475 ft/sec using four injectors evenly spaced around the compressor case upstream of the rotor. The injectors were connected to an external air source. In an actual engine application, the injected air would be supplied with compressor bleed air. The injectors were isolated from the air source by a fast-acting butterfly valve. With the injectors turned off, the compressor was throttled into stall. Air injection was then activated with no change in throttle setting by opening the butterfly valve. The compressor recovered from stall at a fixed throttle setting with the aid of tip injection. The unsteady operating characteristic of the rotor was measured during these tests using high-response pressure sensors located upstream and downstream of the rotor. The figure shows the results, where the unsteady pressure and mass flow are superimposed on the steady operating characteristic. The total injected mass flow was equal to 1.3 percent of the compressor flow. The solid line with no solid squares on it denotes the operating point during the beginning of throttle closure and the initial drop into stall. The gray traces denote the operating point during an additional throttle closure that occurred over the next 1200 rotor revolutions (4 sec). The dashed line denotes the recovery from stall that occurred during 90 rotor revolutions (0.3 sec) after the injectors were activated with no change in throttle setting. Tip injection not only recovers the compressor from stall, but also restores the compressor to its pre-stall level of pressure rise. In contrast, standard stall recovery schemes such as compressor bleed, stator vane actuation, or engine throttle modulation result in a loss of pressure rise across the compressor, which results in a loss of engine power.

Suder, Ken L.↗

Adaptable Constrained Genetic Programming: Extensions and Applications

An evolutionary algorithm applies evolution-based principles to problem solving. To solve a problem, the user defines the space of potential solutions, the representation space. Sample solutions are encoded in a chromosome-like structure. The algorithm maintains a population of such samples, which undergo simulated evolution by means of mutation, crossover, and survival of the fittest principles. Genetic Programming (GP) uses tree-like chromosomes, providing very rich representation suitable for many problems of interest. GP has been successfully applied to a number of practical problems such as learning Boolean functions and designing hardware circuits. To apply GP to a problem, the user needs to define the actual representation space, by defining the atomic functions and terminals labeling the actual trees. The sufficiency principle requires that the label set be sufficient to build the desired solution trees. The closure principle allows the labels to mix in any arity-consistent manner. To satisfy both principles, the user is often forced to provide a large label set, with ad hoc interpretations or penalties to deal with undesired local contexts. This unfortunately enlarges the actual representation space, and thus usually slows down the search. In the past few years, three different methodologies have been proposed to allow the user to alleviate the closure principle by providing means to define, and to process, constraints on mixing the labels in the trees. Last summer we proposed a new methodology to further alleviate the problem by discovering local heuristics for building quality solution trees. A pilot system was implemented last summer and tested throughout the year. This summer we have implemented a new revision, and produced a User's Manual so that the pilot system can be made available to other practitioners and researchers. We have also designed, and partly implemented, a larger system capable of dealing with much more powerful heuristics.

Janikow, Cezary Z.↗

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters (or closwre models) for an initial thermai non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes currently being used at Glenn Research Center (GRC) for Stirling engine modeling are Fluent and CFD-ACE. The porous-media models available in each of these codes are equilibrium models, which assmne that the solid matrix and the fluid are in thermal equilibrium at each spatial location within the porous medium. This is believed to be a poor assumption for the oscillating-flow environment within Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, we non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location end time during the cycle. A NASA regenerator research grant has been providing experimental and computational results to support definition of various empirical coefficients needed in defining a noa-equilibrium, macroscopic, porous-media model (i.e., to define "closure" relations). The grant effort is being led by Cleveland State University, with subcontractor assistance from the University of Minnesota, Gedeon Associates, and Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the NASAlSunpower oscillating-flow test rig also provide experimentally based correlations that are useful in defining parameters for the porous-media model; these correlations are documented in Gedeon Associates' Sage Stirling-Code Manuals. These sources of experimentally based information were used to define the following terms and parameters needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity (including themal dispersion and estimate of tortuosity effects}, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity (including the effect of tortuosity) was also estimated. Determination of the porous-media model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber regenerator matrix. The non-equilibrium porous-media model presented is considered to be an initial, or "draft," model for possible incorporation in commercial CFD codes, with the expectation that the empirical parameters will likely need to be updated once resulting Stirling CFD model regenerator and engine results have been analyzed. The emphasis of the paper is on use of available data to define empirical parameters (and closure models) needed in a thermal non-equilibrium porous-media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates. However, it is anticipated that a thermal non-equilibrium model such as that presented here, when iacorporated in the CFD codes, will improve our ability to accurately model Stirling regenerators with CFD relative to current thermal-equilibrium porous-media models.

Tew, Roy↗

Hybrid LES/RANS Simulation of Transverse Sonic Injection into a Mach 2 Flow

A computational study of transverse sonic injection of air and helium into a Mach 1.98 cross-flow is presented. A hybrid large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) turbulence model is used, with the two-equation Menter baseline (Menter-BSL) closure for the RANS part of the flow and a Smagorinsky-type model for the LES part of the flow. A time-dependent blending function, dependent on modeled turbulence variables, is used to shift the closure from RANS to LES. Turbulent structures are initiated and sustained through the use of a recycling / rescaling technique. Two higher-order discretizations, the Piecewise Parabolic Method (PPM) of Colella and Woodward, and the SONIC-A ENO scheme of Suresh and Huyhn are used in the study. The results using the hybrid model show reasonably good agreement with time-averaged Mie scattering data and with experimental surface pressure distributions, even though the penetration of the jet into the cross-flow is slightly over-predicted. The LES/RANS results are used to examine the validity of commonly-used assumptions of constant Schmidt and Prandtl numbers in the intense mixing zone downstream of the injection location.

Boles, John A.↗

Deflection Analysis of the Space Shuttle External Tank Door Drive Mechanism

Upon observing an abnormal closure of the Space Shuttle s External Tank Doors (ETD), a dynamic model was created in MSC/ADAMS to conduct deflection analyses for assessing whether the Door Drive Mechanism (DDM) was subjected to excessive additional stress, and more importantly, to evaluate the magnitude of the induced step or gap with respect to shuttle s body tiles. To model the flexibility of the DDM, a lumped parameter approximation was used to capture the compliance of individual parts within the drive linkage. These stiffness approximations were then validated using FEA and iteratively updated in the model to converge on the actual distributed parameter equivalent stiffnesses. The goal of the analyses is to determine the deflections in the mechanism and whether or not the deflections are in the region of elastic or plastic deformation. Plastic deformation may affect proper closure of the ETD and would impact aero-heating during re-entry.

Tosto, Michael A.↗

Large Eddy Simulation of Transonic Flow Field in NASA Rotor 37

The current paper reports on numerical investigations on the flow characteristics in a transonic axial compressor, NASA Rotor 37. The flow field was used previously as a CFD blind test case conducted by American Society of Mechanical Engineers in 1994. Since the CFD blind-test exercise, many numerical studies on the flow field in the NASA Rotor 37 have been reported. Although steady improvements have been reported in both numerical procedure and turbulence closure, it is believed that all the important aspects of the flow field have not been fully explained with numerical studies based on the Reynolds Averaged Navier-Stokes (RANS) solution. Experimental data show large dip in total pressure distribution near the hub at downstream of the rotor at 100% rotor speed. Most original numerical solutions from the blind test exercise did not predict this total pressure deficit correctly. This total pressure deficit at the rotor exit was attributed to a hub corner flow separation by the author. Several subsequent numerical studies with different turbulence closure model also calculated this dip in total pressure rise. Also, several studies attributed this total pressure deficit to a small leakage flow coming from the hub in the test article. As the experimental study cannot be repeated, either explanation cannot be validated. The primary purpose of the current investigation is to investigate the transonic flow field with both RANS and a Large Eddy Simulation (LES). The RANS approach gives similar results presented at the original blind test exercise. Although the RANS calculates higher overall total pressure rise, the total pressure deficit near the hub is calculated correctly. The numerical solution shows that the total pressure deficit is due to a hub corner flow separation. The calculated pressure rise from the LES agrees better with the measured total pressure rise especially near the casing area where the passage shock interacts with the tip clearance vortex and flow becomes unsteady due to this interaction. The LES simulation also calculates the total pressure rise deficit near the hub and it agrees well with the measured data.

Hah, Chunill↗