SEARCH · Search NASA
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.
Advanced Stirling-cycle cryogenic apparatus suitable for cooling sensitive infrared detectors to very low temperatures. Working fluid in refrigerator helium. Working fluid compressed and circulated by three piezoelectrically actuated diaphragm pumps offering advantages of greater reliability, relative simplicity, and lower cost.
Topologically Independent Grid, Euler Refinement (TIGER) computer program solves Euler equations of three-dimensional, unsteady flow of inviscid, compressible fluid by numerical integration on unstructured hexahedral coordinate grid refined where necessary to resolve shocks and other details. Hexahedral cells subdivided, each into eight smaller cells, as needed to refine computational grid in regions of high flow gradients. Grid Interactive Refinement and Flow-Field Examination (GIRAFFE) computer program written in conjunction with TIGER program to display computed flow-field data and to assist researcher in verifying specified boundary conditions and refining grid.
Cryogenic damper seals operating close to the liquid-vapor region (near the critical point or slightly su-cooled) are likely to present two-phase flow conditions. Under single phase flow conditions the mechanical energy conveyed to the fluid increases its temperature and causes a phase change when the fluid temperature reaches the saturation value. A bulk-flow analysis for the prediction of the dynamic force response of damper seals operating under two-phase conditions is presented as: all-liquid, liquid-vapor, and all-vapor, i.e. a 'continuous vaporization' model. The two phase region is considered as a homogeneous saturated mixture in thermodynamic equilibrium. Th flow in each region is described by continuity, momentum and energy transport equations. The interdependency of fluid temperatures and pressure in the two-phase region (saturated mixture) does not allow the use of an energy equation in terms of fluid temperature. Instead, the energy transport is expressed in terms of fluid enthalpy. Temperature in the single phase regions, or mixture composition in the two phase region are determined based on the fluid enthalpy. The flow is also regarded as adiabatic since the large axial velocities typical of the seal application determine small levels of heat conduction to the walls as compared to the heat carried by fluid advection. Static and dynamic force characteristics for the seal are obtained from a perturbation analysis of the governing equations. The solution expressed in terms of zeroth and first order fields provide the static (leakage, torque, velocity, pressure, temperature, and mixture composition fields) and dynamic (rotordynamic force coefficients) seal parameters. Theoretical predictions show good agreement with experimental leakage pressure profiles, available from a Nitrogen at cryogenic temperatures. Force coefficient predictions for two phase flow conditions show significant fluid compressibility effects, particularly for mixtures with low mass content of vapor. Under these conditions, an increase on direct stiffness and reduction of whirl frequency ratio are shown to occur. Prediction of such important effects will motivate experimental studies as well as a more judicious selection of the operating conditions for seals used in cryogenic turbomachinery.
Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.
Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.
The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.
To assess engine performance during the testing of Space Shuttle Main Engines (SSMEs), the design of an optimal altitude diffuser is studied for future Space Transportation Systems (STS). For other Space Transportation Systems, rocket propellant using kerosene is also studied. Methane and dodecane have similar reaction schemes as kerosene, and are used to simulate kerosene combustion processes at various temperatures. The equations for the methane combustion mechanism at high temperature are given, and engine combustion is simulated on the General Aerodynamic Simulation Program (GASP). The successful design of an altitude diffuser depends on the study of a sub-scaled diffuser model tested through two-dimensional (2-D) flow-techniques. Subroutines given calculate the static temperature and pressure at each Mach number within the diffuser flow. Implementing these subroutines into program code for the properties of 2-D compressible fluid flow determines all fluid characteristics, and will be used in the development of an optimal diffuser design.
The pulsing inertial oscillation (PIO) model is a nonlinear, time-dependent, translating vortex solution of the inviscid, compressible fluid dynamic equations in the middle troposphere. The translation of this vortex during a pulse is strikingly similar to that of a supercell storm -- a rotating thunderstorm that can generate tornadoes and hail. Two studies were performed to test the hypothesis that some supercell storms are manifestations of a PIO pulse. The first study applied the model to an intense interior draft whose buoyancy was bounded by a temperature excess of +/- 12 K. The peak updraft speed achieved was 41.5 m/ s and the peak Rossby number was 92.9. The study also pointed to an advanced concept for attaining higher values. The second study applied the PIO model to a supercell storm as a whole and succeeded in replicating its bulk properties, such as mesocyclonic circulation, net mass and moisture influxes, and time track. This study also identified a critical feature of the PIO model that could be tested against storm data: The average vertical draft is downward before the turn in the storm track and upward afterwards. In the conventional theory, the average vertical draft is upward from storm inception until dissipation. These differing draft predictions were compared with the best available data, which are surface mesonetwork data. These data were found to support the PIO model. However, surface data alone are not conclusive, and further measurements are warranted.
In the present report the theory of free turbulence propagation and the boundary layer theory are developed for a plane-parallel free stream of a compressible fluid. In constructing the theory use was made of the turbulence hypothesis by Taylor (transport of vorticity) which gives best agreement with test results for problems involving heat transfer in free jets.
In the present paper, the motion of a gas in a plane-parallel Laval nozzle in the neighborhood of the transition from subsonic to supersonic velocities is studied. In a recently published paper, F. I. Frankl, applying the holograph method of Chaplygin, undertook a detailed investigation of the character of the flow near the line of transition from subsonic to supersonic velocities. From the results of Tricomi's investigation on the theory of differential equations of the mixed elliptic-hyperbolic type, Frankl introduced as one of the independent variables in place of the modulus of the velocity, a certain specially chosen function of this modulus. He thereby succeeded in explaining the character of the flow at the point of intersection of the transition line and the axis of symmetry (center of the nozzle) and in studying the behavior of the stream function in the neighborhood of this point by separating out the principal term having, together with its derivatives, the maximum value as compared with the corresponding corrections. This principal term is represented in Frankl's paper in the form of a linear combination of two hypergeometric functions. In order to find this linear combination, it is necessary to solve a number of boundary problems, which results in a complex analysis. In the investigation of the flow with which this paper is concerned, a second method is applied. This method is based on the transformation of the equations of motion to a form that may be called canonical for the system of differential equations of the mixed elliptic-hyperbolic type to which the system of equations of the motion of an ideal compressible fluid refers. By studying the behavior of the integrals of this system in the neighborhood of the parabolic line, the principal term of the solution is easily separated out in the form of a polynomial of the third degree. As a result, the computation of the transitional part of the nozzle is considerably simplified.
The flow about a conical body of an ideal compressible fluid is considered. Assume that the velocity of the oncoming flow at infinity W is directed along the z-axis. The system of Cartesian coordinates x, y, z with origin at the vertex of the cone O is shown. From the considerations,of the dimensional theory, it may be found that along any ray issuing from O the components of the velocity u, v, W+w along the coordinate axes will maintain a constant value. It is further assumed that the conical body has such shape and disposition relative to the flow that u, v, and w are small in comparison with W.
An airbag system for elevating the fuselage of an aircraft off a landing surface a sufficient degree to allow for emergency egress of passengers and crew through ventral emergency exit doors. An airbag assembly made up of a plurality of independent airbags is disposed within the aircraft. When activated, the airbag system deploys the airbags external of the aircraft that elevate the fuselage of the aircraft a sufficient degree to allow for utilizing the ventral emergency exit doors on the fuselage to enable evacuating the passengers and crew. An activation mechanism is connected to the inflation.devices associated with each of the airbags. The activation mechanism generates an electrical signal which activates the inflation devices, which in turn fill the airbags with a compressed fluid, thus expanding the airbags and lifting the fuselage. A crew member initiates the activation of the airbag system through one or more switches.
This paper describes a method to estimate key aerodynamic parameters of single and multistage axial and centrifugal compressors. This mean-line compressor code COMDES provides the capability of sizing single and multistage compressors quickly during the conceptual design process. Based on the compressible fluid flow equations and the Euler equation, the code can estimate rotor inlet and exit blade angles when run in the design mode. The design point rotor efficiency and stator losses are inputs to the code, and are modeled at off design. When run in the off-design analysis mode, it can be used to generate performance maps based on simple models for losses due to rotor incidence and inlet guide vane reset angle. The code can provide an improved understanding of basic aerodynamic parameters such as diffusion factor, loading levels and incidence, when matching multistage compressor blade rows at design and at part-speed operation. Rotor loading levels and relative velocity ratio are correlated to the onset of compressor surge. NASA Stage 37 and the three-stage NASA 74-A axial compressors were analyzed and the results compared to test data. The code has been used to generate the performance map for the NASA 76-B three-stage axial compressor featuring variable geometry. The compressor stages were aerodynamically matched at off-design speeds by adjusting the variable inlet guide vane and variable stator geometry angles to control the rotor diffusion factor and incidence angles.
A high speed centrifugal compressor for compressing fluids includes a permanent magnet synchronous motor (PMSM) having a hollow shaft, the being supported on its ends by ball bearing supports. A permanent magnet core is embedded inside the shaft. A stator with a winding is located radially outward of the shaft. The PMSM includes a rotor including at least one impeller secured to the shaft or integrated with the shaft as a single piece. The rotor is a high rigidity rotor providing a bending mode speed of at least 100,000 RPM which advantageously permits implementation of relatively low-cost ball bearing supports.
A paper presents a theoretical investigation of subsonic and supersonic effects in a Bose-Einstein condensate (BEC). The BEC is represented by a time-dependent, nonlinear Schroedinger equation that includes terms for an external confining potential term and a weak interatomic repulsive potential proportional to the number density of atoms. From this model are derived Madelung equations, which relate the quantum phase with the number density, and which are used to represent excitations propagating through the BEC. These equations are shown to be analogous to the classical equations of flow of an inviscid, compressible fluid characterized by a speed of sound (g/Po)1/2, where g is the coefficient of the repulsive potential and Po is the unperturbed mass density of the BEC. The equations are used to study the effects of a region of perturbation moving through the BEC. The excitations created by a perturbation moving at subsonic speed are found to be described by a Laplace equation and to propagate at infinite speed. For a supersonically moving perturbation, the excitations are found to be described by a wave equation and to propagate at finite speed inside a Mach cone.
A computer program solves the Navier- Stokes equations governing the flow of a viscous, compressible fluid in an axisymmetric or two-dimensional (2D) setting. To obtain solutions more accurate than those generated by prior such programs that utilize regular and/or fixed computational meshes, this program utilizes unstructured (that is, irregular triangular) computational meshes that are automatically adapted to solutions. The adaptation can refine to regions of high change in gradient or can be driven by a novel residual minimization technique. Starting from an initial mesh and a corresponding data structure, the adaptation of the mesh is controlled by use of minimization functional. Other improvements over prior such programs include the following: (1) Boundary conditions are imposed weakly; that is, following initial specification of solution values at boundary nodes, these values are relaxed in time by means of the same formulations as those used for interior nodes. (2) Eigenvalues are limited in order to suppress expansion shocks. (3) An upwind fluctuation-splitting distribution scheme applied to inviscid flux requires fewer operations and produces less artificial dissipation than does a finite-volume scheme, leading to greater accuracy of solutions.
This paper considers the small amplitude unsteady motion of an inviscid non-heat conducting compressible fluid on a transversely sheared mean flow. It extends a previous result given in Goldstein (1978(b) and 1979(a)) which shows that the hydrodynamic component of the motion is determined by two arbitrary convected quantities in the absence of solid surfaces or other external sources. The result is important because it can be used to specify appropriate boundary conditions for unsteady surface interaction problems on transversely sheared mean flows in the same way that the vortical component of the Kovasznay (1953) decomposition is used to specify these conditions for surface interaction problems on uniform mean flows. But unlike the Kovasznay (1953) case the arbitrary convected quantities no longer bear a simple relation to the physical variables. One purpose of this paper is to derive a formula that relates these quantities to the (physically measurable) vorticity and pressure fluctuations in the flow.
This work explores the impact of adding volume dynamics to an electrified turbofan engine model. The zero-dimensional single-input single-output volume elements capture the physics of fluid compressibility of flow in the engine core gas path turbomachinery during transient operation with greater fidelity. Some background regarding the motivation for this effort is provided. A model of the volume element is developed from first principles. The baseline engine model is described and an explanation is given of the modifications made to incorporate the volume elements into the new model. Details of the simulated inputs, designed to elicit exaggerated dynamic responses and to facilitate comparison between the baseline and modified engine models, are explicated. The results of these simulations are presented and discussed, demonstrating that the transient performance and operability of the modified engine model, measured by parameters such as shaft speeds, thrust, and high-pressure compressor (HPC) stall margin, are impacted by the additional dynamics. Commentary is provided regarding potential explanations for the differences in behavior exhibited by the two models. Finally, suggestions are offered for future work that might improve the volume models developed in this effort and leverage the improved fidelity enabled by the modeling of turbomachinery volume dynamics.