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 343 records · Page 19

Acoustic pressures emanating from a turbomachine stage

A knowledge of the acoustic energy emission of each blade row of a turbomachine is useful for estimating the overall noise level of the machine and for determining its discrete frequency noise content. Because of the close spacing between the rotor and stator of a compressor stage, the strong aerodynamic interactions between them have to be included in obtaining the resultant flow field. This paper outlines a three-dimensional theory for determining the discrete frequency noise content of an axial compressor consisting of a rotor and a stator each with a finite number of blades. The lifting surface theory and the linearized equation of an ideal, nonsteady compressible fluid motion are used for thin blades of arbitrary cross section. The combined pressure field at a point of the fluid is constructed by linear addition of the rotor and stator solutions together with an interference factor obtained by matching them for net zero vorticity behind the stage.

Ramachandra, S. M.↗

A fully-coupled implicit method for thermo-chemical nonequilibrium air at sub-orbital flight speeds

A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.

Park, Chul↗

Educational Brief: Using Space for a Better Foundation on Earth Mechanics of Granular Materials

Soils are three-phase composite materials that consist of soil, solid particles, and voids filled with water and/or air. Based on the particle-size distribution, they are generally classified as fine-grained (clays and plastic silts) and coarse-grained soils (nonplastic silts, sand, and gravel). Soil's resistance to external loadings is mainly derived from friction between particles and cohesion. Friction resistance is due to particles' surface-to-surface friction, interlocking, crushing, rearrangement, and dilation (or expansion) during shearing. Cohesion can be due to chemical cementation between particles, electrostatic and electromagnetic forces, and soil-water reaction and equilibrium. The basic factor responsible for the strength of coarse-grained soils is friction. Cohesion can be ignored. This educational brief focuses on measuring shear strength of sands (typical example of coarse-grained soils) where, for the same material, packing density is a main factor to be considered when one asks about the shear strength value. As the external load is applied, the soil's resistance is attained through shearing resistance, which causes the soil volume to increase (expand) or decrease (compress) depending on the initial packing density.

Dooling, Dave↗

Working toward a three-dimensional fatigue closure model for surface cracks

The first reliable elastic fracture mechanics solutions for a surface crack in a plate were obtained by Newman and Raju. The authors, both from the Mechanics of Materials Branch at NASA-Langley, used a highly detailed finite element solution requiring substantial computational resources. Computers have since become more powerful and available; however, many important related problems remain computationally expensive. The problem of three-dimensional fatigue crack growth taking into account plasticity-induced crack closure is one such problem. It is the goal of this research to provide an efficient method to account for three-dimensional crack closure in fatigue. Newman developed a two-dimensional plasticity-induced crack closure model for center cracked specimens. This model requires iterations to determine both the contact solution at each growth step and the extent of the plastic zone at the crack tip. A three-dimensional version of this model would require obtaining these nonlinear variables all along the crack front. This model must be efficient so that repeated calculations can be performed during crack growth simulations. The highly versatile line spring model (LSM) with contact, fatigue, and plasticity will form the basis of the closure model. There are several required additions to past work to address the three-dimensional crack closure problem. Initially, these additions will include (1) an improved LSM to more accurately obtain the crack opening displacement, stress intensity factors, and elastic T-stress near the ends of the surface crack; (2) a method to determine the extent of the plastic zone all along the crack front; (3) a method to determine the contact zone given a perfectly plastic layer of material on the crack surfaces; (4) a method to determine the magnitude of the compressive contact stress; and (5) a way to implement the degree of constraint along the curved crack front. During the summer ASEE program an enhanced LSM was developed. A method similar to that of 'strip synthesis' first introduced by Fujimoto was used. Briefly, the crack opening displacements of 'slices' of the surface crack in a direction parallel to the plate surface are considered in addition to the standard LSM approach that makes use of springs obtained from slices perpendicular to the plate surface. This enhancement is necessary so that an accurate three-dimensional representation of quantities such as contact zone size, plastic zone size, stress intensity factors, T-stress, and crack opening displacement can be determined. By combining results of previous investigations with the LSM, the problem of three-dimensional crack closure will be addressed. In addition to a closure model, the enhanced LSM can be used for many other problems including interacting surface cracks and fatigue crack growth of a through crack with a curved crack front.

Joseph, Paul F.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Accuracies of southwell and force/stiffness methods in the prediction of buckling strength of hypersonic aircraft wing tubular panels

Accuracies of the Southwell method and the force/stiffness (F/S) method are examined when the methods were used in the prediction of buckling loads of hypersonic aircraft wing tubular panels, based on nondestructive buckling test data. Various factors affecting the accuracies of the two methods were discussed. Effects of load cutoff point in the nondestructive buckling tests on the accuracies of the two methods were discussed in great detail. For the tubular panels under pure compression, the F/S method was found to give more accurate buckling load predictions than the Southwell method, which excessively overpredicts the buckling load. It was found that the Southwell method required a higher load cutoff point, as compared with the F/S method. In using the F/S method for predicting the buckling load of tubular panels under pure compression, the load cutoff point of approximately 50 percent of the critical load could give reasonably accurate predictions.

Ko, William L.↗

An implicit multigrid scheme for hypersonic strong-interaction flowfields

An implicit multigrid method for the compressible Navier-Stokes equations is applied to a series of two-dimensional, steady, laminar hypersonic flows over a compression ramp, including both attached and separated flow conditions. The algorithm uses upwind-biased differencing for the convective and pressure terms and central differencing for the shear-stress and heat-transfer terms. An implicit spatially factored scheme is used to advance the solution in terms on a given mesh, and a V-cycle multigrid strategy is used to accelerate convergence. The multigrid algorithm provided an order of magnitude decrease in the computational time required to obtain a converged position of the separation location and enabled a substantially grid-independent result to be obtained for each of the flows considered.

Thomas, James L.↗

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗

MESSENGER Observations of Extreme Magnetic Tail Loading and Unloading During its Third Flyby of Mercury: Substorms?

During MESSENGER's third flyby of Mercury on September 29, 2009, a variable interplanetary magnetic field produced a series of several minute enhancements of the tail magnetic field hy factors of approx. 2 to 3.5. The magnetic field flaring during these intervals indicates that they result from loading of the tail with magnetic flux transferred from the dayside magnetosphere. The unloading intervals were associated with plasmoids and traveling compression regions, signatures of tail reconnection. The peak tail magnetic flux during the smallest loading events equaled 30% of the magnetic flux emanating from Mercury, and may have reached 100% for the largest event. In this case the dayside magnetic shielding is reduced and solar wind flux impacting the surface may be greatly enhanced. Despite the intensity of these events and their similarity to terrestrial substorm magnetic flux dynamics, no energetic charged particles with energies greater than 36 keV were observed.

Slavin, James A.↗

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

Factors affecting the thermal shock behavior of yttria stabilized hafnia based graphite and tungsten composites.

Hafnia-based composites containing either graphite or tungsten were investigated as rocket nozzle throat inserts in solid propellant rocket engines. The thermal shock resistance of these materials is considered in terms of macroscopic thermal conductivity, thermal expansion, modulus of elasticity, and compressive fracture stress. The effect of degree of hafnia stabilization, density, and graphite or tungsten content upon these parameters is discussed. The variation of the ratio of elastic modulus to compressive fracture stress with density and its effect upon thermal shock resistance of these materials are discussed in detail.

Lineback, L. D.↗

Data Compression Algorithm Architecture for Large Depth-of-Field Particle Image Velocimeters

A large depth-of-field particle image velocimeter (PIV) is designed to characterize dynamic dust environments on planetary surfaces. This instrument detects lofted dust particles, and senses the number of particles per unit volume, measuring their sizes, velocities (both speed and direction), and shape factors when the particles are large. To measure these particle characteristics in-flight, the instrument gathers two-dimensional image data at a high frame rate, typically >4,000 Hz, generating large amounts of data for every second of operation, approximately 6 GB/s. To characterize a planetary dust environment that is dynamic, the instrument would have to operate for at least several minutes during an observation period, easily producing more than a terabyte of data per observation. Given current technology, this amount of data would be very difficult to store onboard a spacecraft, and downlink to Earth. Since 2007, innovators have been developing an autonomous image analysis algorithm architecture for the PIV instrument to greatly reduce the amount of data that it has to store and downlink. The algorithm analyzes PIV images and automatically reduces the image information down to only the particle measurement data that is of interest, reducing the amount of data that is handled by more than 10(exp 3). The state of development for this innovation is now fairly mature, with a functional algorithm architecture, along with several key pieces of algorithm logic, that has been proven through field test data acquired with a proof-of-concept PIV instrument.

Bos, Brent↗

Thermal Effects on the Bearing Behavior of Composite Joints

Thermal effects on the pin-bearing behavior of an IM7/PET15 composite laminate are studied comprehensively. A hypothesis presents factors influencing a change in pin-bearing strength with a change in temperature for a given joint design. The factors include the change in the state of residual cure stress, the material properties, and the fastener fit with a change in temperature. Experiments are conducted to determine necessary lamina and laminate material property data for the IM7/PET15 being utilized in this study. Lamina material properties are determined between the assumed stress free temperature of 460 F down to -200 F. Laminate strength properties are determined for several lay-ups at the operating temperatures of 350 F, 70 F, and -200 F. A three-dimensional finite element analysis model of a composite laminate subject to compressive loading is developed. Both the resin rich layer located between lamina and the thermal residual stresses present in the laminate due to curing are determined to influence the state of stress significantly. Pin-bearing tests of several lay-ups were conducted to develop an understanding on the effect of temperature changes on the pin-bearing behavior of the material. A computational study investigating the factors influencing pin-bearing strength was performed. A finite element model was developed and used to determine the residual thermal cure stresses in the laminate containing a hole. Very high interlaminar stress concentrations were observed two elements away from the hole boundary at all three operating temperatures. The pin-bearing problem was modeled assuming a rigid frictionless pin and restraining only radial displacements at the hole boundary. A uniform negative pressure load was then applied to the straight end of the model. A solution, where thermal residual stresses were combined with the state of stress due to pin-bearing loads was evaluated. The presence of thermal residual stresses intensified the interlaminar stresses predicted at the hole boundary in the pin-bearing problem. This dissertation shows that changes in material properties drives pin-bearing strength degradation with increasing temperature.

Walker, Sandra Polesky↗

Sound scattering from atmospheric turbulence

Sound propagation through a turbulent flow field with nonvanishing mean-flow velocity has been studied. The effects of compressibility on sound scattering are discussed. Using an analysis describing moving sources in a fixed frame, particular results are obtained for the far-field scattered sound intensity due to uniformly moving scattering turbulence. A formula is obtained showing the factor by which the scattered sound intensity is amplified due to the motion of the turbulence. In addition, this analysis was extended to include wave propagating through mean flow with constant gradient.

Huang, M. N.↗

Effect of wall suction on the stability of compressible subsonic flows over smooth two-dimensional backward-facing steps

The effect of suction on the stability of compressible flows over backward-facing steps is investigated. Mach numbers up to 0.8 are considered. The results show that continuous suction stabilizes the flow outside the separation bubble, but it destabilizes the flow inside it. Nevertheless, the overall N factor decreases as the suction level increases due to the considerable reduction of the separation bubble. For the same suction flow rate, properly distributed suction strips stabilize the flow more than continuous suction. The size of the separation bubble, and hence its effect on the instability can be considerably reduced by placing strips with high suction velocities in the separation region.

Al-Maaitah, Ayman A.↗

Elevated temperature compressive properties of reaction milled NiAl-AlN and Zr-doped NiAl-AlN composites

Previous studies of a single lot of NiAl powder which had been ground under high intensity conditions in liquid nitrogen (cryomilling) indicated that this processing leads to a high strength, elevated temperature NiAl-AlN composite. Because this was the first known example of the use of the reaction milling process to produce a high temperature composite, the reproducibility of this technique was unknown. Two additional lots of NIAl powder and a lot of a Zr-doped NiAl powder have been cryomilled, and analyses indicate that AlN was formed within a NiAl matrix in all three cases. Compression testing between 1200 K and 1400 K has shown that the deformation resistance of these heats is similar to that of the first lot of NiAl-AlN; thus cryomilling can improve the creep resistance of NiAl by a factor of six. Based on this work, it is concluded that cryomilling of NiAl powder to form high temperature, high strength NiAl-AlN composites is a reproducible process.

Whittenberger, J. D.↗