SEARCH · Engineering Papers
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.
A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.
Circular-to-rectangular transition ducts are used as exhaust system components of aircraft with rectangular exhaust nozzles. Often, the incoming flow of these transition ducts includes a swirling velocity component remaining from the gas turbine engine. Previous transition duct studies have either not included inlet swirl or when inlet swirl was considered, only overall performance parameters were evaluated. Circular-to-rectangular transition duct flows with and without inlet swirl were explored in order to understand the effect of inlet swirl on the transition duct flow field and to provide detailed duct flow data for comparison with numerical code predictions. A method was devised to create a swirling, solid body rotational flow with minimal associated disturbances. Coefficients based on velocities and total and static pressures measured incross stream planes at four axial locations within the transition duct, along with surface static pressure measurements and surface oil film visualization, are presented for both nonswirling and swirling incoming flow. In both cases the inlet centerline Mach number was 0.35. The Reynolds number based on the inlet centerline velocity and duct inlet diameter was 1,547,000 for nonswirling and 1,366,000 for swirling flow. The maximum swirl angle was 15.6 deg. Two pair of counter-rotating side wall vortices appeared in the duct flow without inlet swirl. These vortices were absent in the swirling incoming flow cases.
Circular-to-rectangular transition duct flows with and without inlet swirl were investigated experimentally in order to determine the effect of inlet swirl on the transition duct flow field and to provide detailed duct flow data for comparison with numerical code predictions. Coefficients based on detailed measurements of velocity, total pressure and static pressure, acquired in four cross stream planes within a circular-to-rectangular transition duct, with and without inlet swirl, are presented, as are surface static pressure and surface oil film visualization results.
The PARC3D code was selected by the authors to analyze a variety of complex and high-speed flow configurations. Geometries considered for code validation include ramps and corner flows, which are characteristic of inlets and nozzles. Flows with Mach numbers of 3-14 were studied. Both two- and three-dimensional experimental data for shock-boundary-layer interaction were considered to validate the code. A detailed comparison of various flow parameters with available experimental data is presented; agreement between the solutions and the experimental data in terms of pitot pressure profiles, yaw-angle distributions, static pressures, and skin friction is found to be very good. In addition, two- and three-dimensional flow calculations were performed for a hypersonic nozzle. Comparison of the wall pressure results with the published solutions is made for the two-dimensional case.
One possible low speed propulsion system for the National Aerospace Plane is a liquid air cycle engine (LACE). The LACE system uses the heat sink in the liquid hydrogen propellant to liquefy air in a heat exchanger which is then pumped up to high pressure and used as the oxidizer in a hydrogen liquid air rocket. The inlet airstream must be dehumidified or moisture could freeze on the cryogenic heat exchangers and block them. The main objective of this research has been to develop a computer simulation of the cold tube/antifreeze-spray water alleviation system and to verify the model with experimental data. An experimental facility has been built and humid air tests were conducted on a generic heat exchanger to obtain condensing data for code development. The paper describes the experimental setup, outlines the method of calculation used in the code, and presents comparisons of the calculations and measurements. Cause of discrepancies between the model and data are explained.
This paper presents an overview of high angle-of-attack tests of a full-scale F/A-18 in the 80- by 120-Foot Wind Tunnel of the National Full-Scale Aerodynamic Complex at NASA Ames Research Center at Moffett Field, California. A production aircraft was tested over an angle-of-attack range of 18 to 50 deg and at wind speeds of up to 100 knots. These tests had three primary test objectives. Pneumatic and mechanical forebody flow control devices were tested at full-scale and shown to produce significant yawing moments for lateral control of the aircraft at high angles of attack. Mass flow requirements for the pneumatic system were found to scale with freestream density and speed rather than freestream dynamic pressure. Detailed measurements of the pressures buffeting the vertical tail were made and spatial variations in the buffeting frequency were found. The LEX fence was found to have a significant effect on the frequency distribution on the outboard surface of the vertical fin. In addition to the above measurements, an extensive set of data was acquired for the validation of computational fluid dynamics codes and for comparison with flight test and small-scale wind tunnel test results.
Computational Fluid Dynamics (CFD), because of its unique ability to predict complex three-dimensional flows, is being applied with increasing frequency in the aerospace industry. Currently, no consistent code validation procedure is applied within the industry. Such a procedure is needed to increase confidence in CFD and reduce risk in the use of these codes as a design and analysis tool. This final contract report defines classifications for three levels of code validation, directly relating the use of CFD codes to the engineering design cycle. Evaluation criteria by which codes are measured and classified are recommended and discussed. Criteria for selecting experimental data against which CFD results can be compared are outlined. A four phase CFD code validation procedure is described in detail. Finally, the code validation procedure is demonstrated through application of the REACT CFD code to a series of cases culminating in a code to data comparison on the Space Shuttle Main Engine High Pressure Fuel Turbopump Impeller.
As the size and complexity of three dimensional volume grids increases, there is a growing need for fast and efficient 3D volumetric elliptic grid solvers. Present day solvers are limited by computational speed and do not have all the capabilities such as interior volume grid clustering control, viscous grid clustering at the wall of a configuration, truncation error limiters, and convergence optimization residing in one code. A new volume grid generator, 3DMAGGS (Three-Dimensional Multi-Block Advanced Grid Generation System), which is based on the 3DGRAPE code, has evolved to meet these needs. This is a manual for the usage of 3DMAGGS and contains five sections, including the motivations and usage, a GRIDGEN interface, a grid quality analysis tool, a sample case for verifying correct operation of the code, and a comparison to both 3DGRAPE and GRIDGEN3D. Since it was derived from 3DGRAPE, this technical memorandum should be used in conjunction with the 3DGRAPE manual (NASA TM-102224).
The physics of the circulation controlled airfoil is complex and poorly understood, particularly with regards to jet stall, which is the eventual breakdown of lift augmentation by the jet at some sufficiently high blowing rate. The present paper describes the numerical simulation of stalled and unstalled flows over a two-dimensional circulation controlled airfoil using a fully implicit Navier-Stokes code, and the comparison with experimental results. Mach numbers of 0.3 and 0.5 and jet total to freestream pressure ratios of 1.4 and 1.8 are investigated. The Baldwin-Lomax and k-epsilon turbulence models are used, each modified to include the effect of strong streamline curvature. The numerical solutions of the post-stall circulation controlled airfoil show a highly regular unsteady periodic flowfield. This is the result of an alternation between adverse pressure gradient and shock induced separation of the boundary layer on the airfoil trailing edge.
Interpretation of the vast data flow from the recent and ongoing Mars missions requires extensive modeling activity to support the investigations and conclusions drawn from the observations. We have initiated a campaign to compare atmospheric radiation schemes for the lower Martian atmosphere with leastcompromise reference line-by-line calculations (Meadows and Crisp 1996, Crisp 1986). By using fixed reference conditions we are introducing an intercomparison framework to enable atmospheric modeling teams to compare their modeling results with the reference line-by-line calculations, and to improve their models. The atmosphere of Mars is thin and therefore its response to local radiative forcing is strong when compared to that of the Earth s troposphere. Airborne dust and also the water vapor give significant contributions to the radiative fluxes. Thus models of the Martian atmosphere should have fairly accurate radiative transfer (RT) algorithms. This is specifically important for the net radiative flux at the surface, which is the main driver for the surface temperature evolution, since the surface turbulent fluxes are small due to the thin atmosphere (e.g. Haberle et al. 1993, Savij arvi 1991a, Savij arvi 1991b, Savij arvi 1999). However, the accuracy of the RT schemes has not been rigorously tested in many cases as reference calculations have been rare. We introduce a set of reference cases for Martian atmospheric conditions with profiles for temperature, pressure, dust optical depth (Conrath 1975, Clancy et al. 1995, OckertBell et al. 1997, Smith et al. 2000), water, CO2,O3,CO and O2 (Table 1). The reference atmosphere uses a diurnally and globally averaged thermal structure derived from Mariner 9 IRIS observations acquired during late southern summer. Thermal structures for nighttime and daytime conditions will be included shortly. Computationally expensive reference simulations were performed by applying a spectrum resolving (line-by-line multiple scattering) model, SRM, to give the least-compromise base for comparisons. We call for other groups to join in future comparisons. This resembles the International Comparison of Radiation Codes for Climate Models (ICRCCM; e.g. Ellingson and Fouquart 1991, Ellingson et al. 1991). The ICRCCM led to many improvements in the Earth GCM radiation codes. The intercomparison framework presented in this paper is a step toward that direction for Mars. We have already five modeling teams participating in this intercomparison exercice. Additional teams investigating the Martian atmosphere are cordially invited to join this effort of comparing the radiation schemes in fixed conditions. The results achieved will be reported and some additional features to be included in the intercomparison will be discussed.
Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.
NASA and Industry are performing vehicle studies of configurations with low sonic boom pressure signatures. The computational analyses of modern configuration designs have matured to the point where there is confidence in the prediction of the pressure signature from the front of the vehicle, but uncertainty in the aft signatures with often greater boundary layer effects and nozzle jet pressures. Wind tunnel testing at significantly lower Reynolds numbers than in flight and without inlet and nozzle jet pressures make it difficult to accurately assess the computational solutions of flight vehicles. A wind tunnel test in the NASA Ames 9- by 7-Foot Supersonic Wind Tunnel from Mach 1.6 to 2.0 will be used to assess the effects of shocks from components passing through nozzle jet plumes on the sonic boom pressure signature and provide datasets for comparison with CFD codes. A large number of high-fidelity numerical simulations of wind tunnel test models with a variety of shock generators that simulate horizontal tails and aft decks have been studied to provide suitable models for sonic boom pressure measurements using a minimally intrusive pressure rail in the wind tunnel. The computational results are presented and the evolution of candidate wind tunnel models is summarized and discussed in this paper.
This paper investigates the effects of axial liner displacement on modal attenuation and scattering. Tests are conducted with a uniform liner mounted in the NASA Langley Curved Duct Test Rig, with a multimode source at flow conditions of Mach 0.0 and 0.3. The impedance of the liner is educed from data acquired in the NASA Langley Grazing Flow Impedance Tube. This impedance is then used as input to the CDUCT-LaRC propagation code to enable comparisons of predicted and measured modal content upstream and downstream of the liner as the axial displacement is varied. Results show that the liner position has a clear effect on the sound field as the number of modes is increased from two to six. The comparison of measured and predicted modal content is very good when only a few modes are present. For frequencies in which more modes are present, this comparison is less favorable very near resonance, but improves for frequencies away from resonance.
The overall objective of the current effort at NASA GRC is to evaluate, develop, and apply methodologies suitable for modeling intra-engine trace chemical changes over post combustor flow path relevant to the pollutant emissions from aircraft engines. At the present time, the focus is the high pressure turbine environment. At first, the trace chemistry model of CNEWT were implemented into GLENN-HT as well as NCC. Then, CNEWT, CGLENN-HT, and NCC were applied to the trace species evolution in a cascade of Cambridge University's No. 2 rotor and in a turbine vane passage. In general, the results from these different codes provide similar features. However, the details of some of the quantities of interest can be sensitive to the differences of these codes. This report summaries the implementation effort and presents the comparison of the No. 2 rotor results obtained from these different codes. The comparison of the turbine vane passage results is reported elsewhere. In addition to the implementation of trace chemistry model into existing CFD codes, several pre/post-processing tools that can handle the manipulations of the geometry, the unstructured and structured grids as well as the CFD solutions also have been enhanced and seamlessly tied with NCC, CGLENN-HT, and CNEWT. Thus, a complete CFD package consisting of pre/post-processing tools and flow solvers suitable for post-combustor intra-engine trace chemistry study is assembled.
Hypersonic flows over straight and bent biconics are calculated for a range of freestream conditions in which the gas behind the shock is treated as either perfect or real. The Parabolized-Navier-Stokes (PNS) equations form the basis of the approximation scheme. Good comparisons with experimental data for pressures, forces and moments, heat transfer, and oil-flow patterns serve to validate the perfect-gas version of the code. Circumferential velocity vector plots further aid in the interpretation of leeside oil-flow patterns. A variable-effective-gamma (VEG) option is implemented for the real-gas calculations. Gamma, now defined as the ratio of enthalpy to internal energy, is determined from a locally valid linear relation in enthalpy and pressure at every mesh point which in turn is calculated from a benchmark equilibrium code. The VEG option is easily incorporated into a host PNS code because it uses the underlying perfect-gas structure of the code. Comparisons to experimental data for heat transfer and shock shape in high enthalpy air have been obtained using the VEG option.
Computational fluid dynamics (CFD) software that solves the Reynolds-averaged Navier-Stokes (RANS) equations has been in routine use for more than a quarter of a century. It is currently employed not only for basic research in fluid dynamics, but also for the analysis and design processes in many industries worldwide, including aerospace, automotive, power generation, chemical manufacturing, polymer processing, and petroleum exploration. A key feature of RANS CFD is the turbulence model. Because the RANS equations are unclosed, a model is necessary to describe the effects of the turbulence on the mean flow, through the Reynolds stress terms. The turbulence model is one of the largest sources of uncertainty in RANS CFD, and most models are known to be flawed in one way or another. Alternative methods such as direct numerical simulations (DNS) and large eddy simulations (LES) rely less on modeling and hence include more physics than RANS. In DNS all turbulent scales are resolved, and in LES the large scales are resolved and the effects of the smallest turbulence scales are modeled. However, both DNS and LES are too expensive for most routine industrial usage on today's computers. Hybrid RANS-LES, which blends RANS near walls with LES away from walls, helps to moderate the cost while still retaining some of the scale-resolving capability of LES, but for some applications it can still be too expensive. Even considering its associated uncertainties, RANS turbulence modeling has proved to be very useful for a wide variety of applications. For example, in the aerospace field, many RANS models are considered to be reliable for computing attached flows. However, existing turbulence models are known to be inaccurate for many flows involving separation. Research has been ongoing for decades in an attempt to improve turbulence models for separated and other nonequilibrium flows. When developing or improving turbulence models, both verification and validation are important steps in the process. Verification insures that the CFD code is solving the equations as intended (no errors in the implementation). This is typically done either through the method of manufactured solutions (MMS) or through careful step-by-step comparisons with other verified codes. After the verification step is concluded, validation is performed to document the ability of the turbulence model to represent different types of flow physics. Validation can involve a large number of test case comparisons with experiments, theory, or DNS. Organized workshops have proved to be valuable resources for the turbulence modeling community in its pursuit of turbulence modeling verification and validation. Workshop contributors using different CFD codes run the same cases, often according to strict guidelines, and compare results. Through these comparisons, it is often possible to (1) identify codes that have likely implementation errors, and (2) gain insight into the capabilities and shortcomings of different turbulence models to predict the flow physics associated with particular types of flows. These are valuable lessons because they help bring consistency to CFD codes by encouraging the correction of faulty programming and facilitating the adoption of better models. They also sometimes point to specific areas needed for improvement in the models. In this paper, several recent workshops are summarized primarily from the point of view of turbulence modeling verification and validation. Furthermore, the NASA Langley Turbulence Modeling Resource website is described. The purpose of this site is to provide a central location where RANS turbulence models are documented, and test cases, grids, and data are provided. The goal of this paper is to provide an abbreviated survey of turbulence modeling verification and validation efforts, summarize some of the outcomes, and give some ideas for future endeavors in this area.
11 tests were selected from two microcomputer-based performance test batteries because previously these tests exhibited rapid stability (less than 10 min, of practice) and high retest reliability efficiencies (r greater than 0.707 for each 3 min. of testing). The battery was administered three times to each of 108 college students (48 men and 60 women) and a factor analysis was performed. Two of the three identified factors appear to be related to information processing ("encoding" and "throughput/decoding"), and the third named an "output/speed" factor. The spatial, memory, and verbal tests loaded on the "encoding" factor and included Grammatical Reasoning, Pattern Comparison, Continuous Recall, and Matrix Rotation. The "throughput/decoding" tests included perceptual/numerical tests like Math Processing, Code Substitution, and Pattern Comparison. The output speed factor was identified by Tapping and Reaction Time tests. The Wonderlic Personnel Test was group administered before the first and after the last administration of the performance tests. The multiple Rs in the total sample between combined Wonderlic as a criterion and less than 5 min. of microcomputer testing on Grammatical Reasoning and Math Processing as predictors ranged between 0.41 and 0.52 on the three test administrations. Based on these results, the authors recommend a core battery which, if time permits, would consist of two tests from each factor. Such a battery is now known to permit stable, reliable, and efficient assessment.
The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.