Search NASA⌕ Search

SEARCH · Search NASA

Results for “compressibility factor”

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 577 records · Page 32

Implicit time-marching solution of the Navier-Stokes equations for thrust reversing and thrust vectoring nozzle flows

An implicit finite volume method is investigated for the solution of the compressible Navier-Stokes equations for flows within thrust reversing and thrust vectoring nozzles. Thrust reversing nozzles typically have sharp corners, and the rapid expansion and large turning angles near these corners are shown to cause unacceptable time step restrictions when conventional approximate factorization methods are used. In this investigation these limitations are overcome by using second-order upwind differencing and line Gauss-Siedel relaxation. This method is implemented with a zonal mesh so that flows through complex nozzle geometries may be efficiently calculated. Results are presented for five nozzle configurations including two with time varying geometries. Three cases are compared with available experimental data and the results are generally acceptable.

Imlay, S. T.↗

Deep space communications, weather effects, and error control

Deep space telemetry is and will remain signal-to-noise limited and vulnerable to interference. A need exists to increase received signal power and decrease noise. This includes going to Ka-band in the mid-1990's to increase directivity. The effects of a wet atmosphere can increase the noise temperature by a factor of 5 or more, even at X-band, but the order of magnitude increase in average data rate obtainable at Ka-band relative to X-band makes the increased uncertainty a good trade. Lowbit error probabilities required by data compression are available both theoretically and practically with coding, at an infinitesimal power penalty rather than the 10 to 15 dB more power required to reduce error probabilities without coding. Advances are coming rapidly in coding, as with the new constraint-length 15 rate 1/4 convolutional code concatenated with the already existing Reed-Solomon code to be demonstrated on Galileo. In addition, high density spacecraft data storage will allow selective retransmissions, even from the edge of the Solar System, to overcome weather effects. In general, deep space communication was able to operate, and will continue to operate, closer to theoretical limits than any other form of communication. These include limits in antenna area and directivity, system noise temperature, coding efficiency, and everything else. The deep space communication links of the mid-90's and beyond will be compatible with new instruments and compression algorithms and represent a sensible investment in an overall end-to-end information system design.

Posner, Edward C.↗

Shock interference heating in scramjet engines

Experimental and analytical research sponsored by the NASA Langley Research center and the NASP Structures Technology Maturation Program to define critical aerothermal loads for the NASP engine is summarized. Presented is a review of (1) shock-shock interaction on the engine cowl leading edge that results in a supersonic jet impinging on the leading edge surface and causes the heat transfer rate to be amplified by a factor of 30 or more over the undisturbed (no shock interaction) flow stagnation point heat transfer rate, (2) the effectiveness of supersonic film cooling with and without the effects of an impinging oblique shock wave, and (3) oblique shock impingement in an axial compression corner.

Wieting, Allan R.↗

Effects of bleed-hole geometry and plenum pressure on three-dimensional shock-wave/boundary-layer/bleed interactions

A numerical study was performed to investigate 3D shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. This study was focused on how bleed-hole geometry and pressure ratio across bleed holes affect the bleed rate and the physics of the flow in the vicinity of the holes. The aspects of the bleed-hole geometry investigated include angle of bleed hole and the number of bleed holes. The plenum/freestream pressure ratios investigated range from 0.3 to 1.7. This study is based on the ensemble-averaged, 'full compressible' Navier-Stokes (N-S) equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the ensemble-averaged N-S equations were obtained by an implicit finite-volume method using the partially-split, two-factored algorithm of Steger on an overlapping Chimera grid.

Chyu, Wei J.↗

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗

Multiple grid problems on concurrent-processing computers

Three computer codes were studied which make use of concurrent processing computer architectures in computational fluid dynamics (CFD). The three parallel codes were tested on a two processor multiple-instruction/multiple-data (MIMD) facility at NASA Ames Research Center, and are suggested for efficient parallel computations. The first code is a well-known program which makes use of the Beam and Warming, implicit, approximate factored algorithm. This study demonstrates the parallelism found in a well-known scheme and it achieved speedups exceeding 1.9 on the two processor MIMD test facility. The second code studied made use of an embedded grid scheme which is used to solve problems having complex geometries. The particular application for this study considered an airfoil/flap geometry in an incompressible flow. The scheme eliminates some of the inherent difficulties found in adapting approximate factorization techniques onto MIMD machines and allows the use of chaotic relaxation and asynchronous iteration techniques. The third code studied is an application of overset grids to a supersonic blunt body problem. The code addresses the difficulties encountered when using embedded grids on a compressible, and therefore nonlinear, problem. The complex numerical boundary system associated with overset grids is discussed and several boundary schemes are suggested. A boundary scheme based on the method of characteristics achieved the best results.

Eberhardt, D. S.↗

Modeling of transient heat pipe operations

An analysis of the steady, compressible, one-dimensional, laminar flow of sodium vapor is presented for a case of a flat plate-type heat pipe with asymmetrical boundary conditions. In addition, shear stress at the liquid-vapor interface, variations of vapor quality, and momentum and energy factors are considered. A similarity solution for a semiporous channel is used to provide the velocity profile at cross sections.

Colwell, Gene T.↗

Alterations in crack-tip deformation during variable-amplitude fatigue crack growth

Fatigue crack growth rate changes caused by single overloads and an overload followed by an underload have been studied in the high-strength aluminum alloy 7091. Crack-tip plasticity parameters were measured at each step in the loading sequence using the stereoimaging technique. Effective stress-intensity factor was measured with high resolution immediately at the crack tip both before an overload and during the subsequent growth rate retardation period. Crack-tip opening displacement and strain were altered in the same way as growth rate following the overload/underload events, and similitude with respect to crack-tip plasticity was preserved during the growth retardation period. Effective stress-intensity factor was found to correlate well with crack growth rate. Strains within the plastic zone were converted to stress at each step in load change. These stresses were then summed to determine the zone of residual stress caused by the overload. The region of compressive residual stress ahead of the crack tip was found to increase with overload severity but was limited to yield stress in magnitude. An underload shrank the spatial extent of the residual stresses caused by the overload.

Davidson, David L.↗

Boundary-layer stability analysis of Langley Research Center 8-foot LFC experimental data

An analytical study of linear-amplifying instabilities of a laminar boundary layer as found in the experimental data of the LaRC/8-foot laminar-flow control (LFC) experiment was completed and the results are presented. The LFC airfoil used for this experiment was a swept, supercritical design which removed suction air through spanwise slots. The amplification of small disturbances by linear processes on a swept surface such as this can be due to either Tollmien-Schlichting (TS) and/or crossflow (CF) mechanisms. This study consists of the examination of these two instabilities by both the commonly used incompressible (SALLY and MARIA) analysis and the more involved compressible (COSAL) analysis. A wide range of experimental test conditions with variations in Mach number, Reynolds number, and suction distributions were available for this study. Experimentally determined transition locations were found from thin-film techniques and were used to correlate the n-factors at transition for the range of test cases.

Berry, Scott↗

Prediction of steady and unsteady asymmetric vortical flows around cones

Steady and unsteady, supersonic asymmetric vortical flows and their passive control around circular and noncircular cones are considered in this paper. These problems are formulated by using the unsteady, compressible, single and double, thin-layer. Navier-Stokes equations. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme, either in a pseudotime stepping or in an accurate-time stepping. An implicit, approximately-factored, central-difference finite-volume scheme has also been used to validate some applications of the upwind scheme. Steady asymmetric vortical flows have been predicted by using random and controlled disturbances for circular and noncircular cones. Unsteady asymmetric vortex-shedding flows have also been predicted, for the first time, using time-accurate solutions, for circular and noncircular cones. Control of flow asymmetry have been demonstrated computationally, for the first time, by inserting a vertical fin the leeward plane of geometric symmetry.

Kandil, Osama A.↗

Creep Properties of NiAl-1Hf Single Crystals Re-Investigated

NiAl-1Hf single crystals have been shown to be quite strong at 1027 C, with strength levels approaching those of advanced Ni-based superalloys. Initial testing, however, indicated that the properties might not be reproducible. Study of the 1027 C creep behavior of four different NiAl-1Hf single-crystal ingots subjected to several different heat treatments indicated that strength lies in a narrow band. Thus, we concluded that the mechanical properties are reproducible. Recent investigations of the intermetallic NiAl have confirmed that minor alloying additions combined with single-crystal growth technology can produce elevated temperature strength levels approaching those of Ni-based superalloys. For example, General Electric alloy AFN 12 {Ni-48.5(at.%) Al-0.5Hf-1Ti-0.05Ga} has a creep rupture strength equivalent to Rene 80 combined with a approximately 30-percent lower density, a fourfold improvement in thermal conductivity, and the ability to form a self-protective alumina scale in aggressive environments. Although the compositions of strong NiAl single crystals are relatively simple, the microstructures are complex and vary with the heat treatment and with small ingot-toingot variations in the alloy chemistry. In addition, initial testing suggested a strong dependence between microstructure and creep strength. If these observations were true, the ability to utilize NiAl single-crystal rotating components in turbine machinery could be severely limited. To investigate the possible limitations in the creep response of high-strength NiAl single crystals, the NASA Glenn Research Center at Lewis Field initiated an in depth investigation of the effect of heat treatment on the microstructure and subsequent 1027 C creep behavior of [001]-oriented NiAl-1Hf with a nominal chemistry of Ni-47.5Al-1Hf-0.5Si. This alloy was selected since four ingots, grown over a number of years and possessing slightly different compositions, were available for study. Specimens taken from the ingots were subjected to several heat treatment schedules, examined by transmission electron microscopy, and tested in both compression and tension. An example of the microstructure found in a [001]-oriented NiAl-1Hf specimen after a solution treatment at 1317 C for 50 hr followed by air cooling is illustrated in the image on the left, where the NiAl matrix contains a uniform distribution of nanometer-scale Gphase (Ni16Hf6Si7) precipitates. Other heat treating schedules produced microstructures with nanometer-sized G-phase cubes and plates or, in an extreme case, produced a microstructure with all the G-phase converted to Heusler (Ni2AlHf) particles. The results of 1027 C creep strength and strain rate testing are illustrated which summarizes data from tensile and compressive testing of samples cut from all four NiAl-1Hf ingots and subjected to a variety of heat treatment schedules. With one exception, all the strength values lie in a narrow band that spans six orders of magnitude in strain rate. The only factor that produced results outside of this band was the heat treatment schedule that dissolved all the G-phase and replaced it with Heusler precipitates. The results portrayed in this figure lead to the important practical conclusion that the elevated-temperature creep properties of NiAl-1Hf single crystals are reproducible and are not affected by small variations in alloy chemistry from ingot to ingot or by different initial distributions of G-phase in the heat-treated alloy. The only variable in this study that produced a significant and delerious effect on mechanical strength was a post-solution heat treatment that lead to the complete disappearance of the G-phase in favor of Heusler precipitates.

Whittenberger, J. Daniel↗

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗

Geomagnetically Induced Currents Caused by Interplanetary Shocks with Different Impact Angles and Speeds

The occurrence of geomagnetically induced currents (GICs) poses serious threats to modern technological infrastructure. Large GICs result from sharp variations of the geomagnetic field (dB∕dt) caused by changes of large-scale magnetospheric and ionospheric currents. Intense dB∕dt perturbations are known to occur often in high-latitude regions as a result of storm time substorms. Magnetospheric compressions usually caused by interplanetary shocks increase the magnetopause current leading to dB∕dt perturbations more evident in midlatitude to low-latitude regions, while they increase the equatorial electrojet current leading to dB∕dt perturbations in dayside equatorial regions. We investigate the effects of shock impact angles and speeds on the subsequent dB∕dt perturbations with a database of 547 shocks observed at the L1 point. By adopting the threshold of dB∕dt = 100 nT/min, identified as a risk factor to power systems, we find that dB∕dt generally surpasses this threshold when following impacts of high-speed and nearly frontal shocks in dayside high-latitude locations. The same trend occurs at lower latitudes and for all nightside events but with fewer high-risk events. Particularly, we found nine events in equatorial locations with dB∕dt > 100 nT/min. All events were caused by high-speed and nearly frontal shock impacts and were observed by stations located around noon local time. These high-risk perturbations were caused by sudden strong and symmetric magnetospheric compressions, more effectively intensifying the equatorial electrojet current, leading to sharp dB∕dt perturbations. We suggest that these results may provide insights for GIC forecasting aiming at preventing degradation of power systems due to GICs.

Oliveira, D. M.↗

COSAL: A black-box compressible stability analysis code for transition prediction in three-dimensional boundary layers

A fast computer code COSAL for transition prediction in three dimensional boundary layers using compressible stability analysis is described. The compressible stability eigenvalue problem is solved using a finite difference method, and the code is a black box in the sense that no guess of the eigenvalue is required from the user. Several optimization procedures were incorporated into COSAL to calculate integrated growth rates (N factor) for transition correlation for swept and tapered laminar flow control wings using the well known e to the Nth power method. A user's guide to the program is provided.

Malik, M. R.↗

Calculation of aerodynamic forces on a propeller in pitch or yaw

An analysis was made to determine the applicability of existing propeller theory and the theory of oscillating airfoils to the problem of determining the magnitude of the forces on propellers in pitch or yaw. Strip calculations using compressible airfoil characteristics were first made as though steady-state conditions existed successively at several blade positions of the propeller blades during one revolution. A theory of oscillating airfoils in pulsating incompressible potential flow was then considered from which it was possible to determine factors which would modify the steady-state forces.

PROPELLERS - PITCH AND YAW↗

Experimental Investigation on Metallic Droplet Behavior in Molten BOF Slag

In order to better understand the metallic droplet behavior during a slag treatment process, a physical modeling based on the similarity principle was performed in a transparent scaled-down vessel at room temperature. Paraffin oil, 20 wt pct copper sulfate solution, and compressed air were used to simulate the molten slag, metallic droplet, and carrier gas, respectively. The droplets injected into paraffin oil during the experiment were captured by a high speed camera and were analyzed by Image Pro Plus software to obtain the droplet size distribution. The critical droplet size in the physical modeling and slag treatment process is quantitatively correlated. The results show that droplet breakage phenomenon is dominant over its coalescence in the current industrial practice, and droplet breakage is enhanced with increasing gas flow rate and/or lance depth. No significant effect of the nozzle configuration was found on the droplet breakage and coalescence. The droplet size distribution varies with the lance position. Gas flow rate and lance depth are the most important factors for droplet breakage, the extent of which can be reduced through a proper selection of the operational conditions. A linear relationship between the droplet size and the input energy flux is obtained.

36 MATERIALS SCIENCE↗

Additively manufactured β-Ti5553 with laser powder bed fusion: Microstructures and mechanical properties of bulk and lattice parts

Ti5553 (Ti-5Al-5Mo-5V-3Cr wt%) is a titanium alloy widely used for its high strength-to-weight ratio and good formability at elevated temperatures. Unlike Ti-6Al-4V, Ti5553 does not undergo martensitic transformation, preventing cracking of brittle martensite upon rapid cooling. This makes it a strong candidate for additive manufacturing (AM), particularly laser powder bed fusion (L-PBF). L-PBF offers the unique opportunity to make fine lattice structures to reduce component weight. Despite the growing field of AM, there have been limited studies on L-PBF Ti5553 lattices and how their properties differ from the bulk. Here the present work addresses this knowledge gap by investigating microstructures and properties of L-PBF bulk and lattice parts and the effect of post L-PBF heat treatments. Electron microscopy and mechanical testing show that the high dislocation density formed during L-PBF increases bulk part’s yield strength by approximately 100 MPa compared to the conventional alloy. Digital image correlation during compression testing of octet truss lattices reveals a layer-by-layer failure mode. Compared to the bulk, the lattice contains copious ω nanoprecipitation, weaker <001> texture, smaller average grain sizes, and larger content of high-angle grain boundaries. These features elicit differences in Taylor factor distributions for the lattice depending on load direction, underlining challenges in predicting lattice mechanical response based on bulk properties. By examining the processing-structure-property relationships in the bulk and lattice, the present results delineate their microstructural and mechanical differences and establish a benchmark for the future design applications of L-PBF Ti5553.

36 MATERIALS SCIENCE↗

Three dimensional unsteady aerodynamics and aeroelastic response of advanced turboprops

A method for the prediction of steady and unsteady aerodynamic loads and aeroelastic response of advanced turboprops is presented. The aerodynamic analysis uses three dimensional unsteady linearized compressible flow theory to compute the blade pressure distribution. The aeroelastic analysis is based on a normal mode representation of the structure. The method is applicable to both conventional and advanced turbo-prop configurations, provided that blade stall and transonic shock waves are not important factors. Aerodynamic results are presented which validate the model in various limits by comparisons to alternative theories and experimental data. Finally, results of a stability analysis of an advanced turboprop are given, with comparisons to measurements made at NASA Lewis Research Center.

Williams, M. H.↗