Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gas compression”

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 667 records · Page 37

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗

Modelling of Spherical Gas Bubble Oscillations and Sonoluminescence

The discovery of single-bubble sonoluminescence has led to a renewed interest in the forced radial oscillations of gas bubbles. Many of the more recent studies devoted to this topic have used several simplifications in the modelling, and in particular in accounting for liquid compressibility and thermal processes in the bubble. In this paper the significance of these simplifications is explored by contrasting the results of Lohse and co-workers with those of a more detailed model. It is found that, even though there may be little apparent difference between the radius-versus time behaviour of the bubble as predicted by the two models, quantities such as the spherical stability boundary and the threshold for rectified diffusion are affected in a quantitatively significant way. These effects are a manifestation of the subtle dependence upon dissipative processes of the phase of radial motion with respect to the driving sound field. The parameter space region, where according to the theory of Lohse and co-workers, sonoluminescence should be observable, is recalculated with the new model and is found to be enlarged with respect to the earlier estimate. The dependence of this parameter region on sound frequency is also illustrated.

Prosperetti, A.↗

Stirring Up an Elastic Fluid: Critical Viscosity of Xenon-2 (CVX-2)

Whipped cream stays in place even when turned upside down. Yet it readily flows through the nozzle of a spray can to reach the dessert plate. This demonstrates the phenomenon of shear thinning that is important to many industrial and physical processes. Paints, film emulsions, and other complex solutions that are highly viscous under normal conditions but become thin and flow easily under shear forces. A simple fluid, such as water, does not exhibit shear thinning under normal conditions. Very close to the liquid-vapor critical point, where the distinction between liquid and vapor disappears, the fluid becomes more complex and is predicted to display shear thinning. At the critical point, xenon atoms interact over long distances in a classical model of cooperative phenomena. Physicists rely on this system to learn how long-range order arises. The Critical Viscosity of Xenon Experiment (CVX-2) will measure the viscous behavior of xenon, a heavy inert gas used in flash lamps and ion rocket engines, at its critical point. Although it does not easily combine with other chemicals, its viscosity at the critical point can be used as a model for a range of fluids. Viscosity originates from the interactions of individual molecules. It is so complicated that, except for the simplest gas, it cannot be calculated accurately from theory. Tests with critical fluids can provide key data, but are limited on Earth because critical fluids are highly compressed by gravity. CVX-2 employs a tiny metal screen vibrating between two electrodes in a bath of critical xenon. The vibrations and how they dampen are used to measure viscosity. CVX flew on STS-85 (1997), where it revealed that, close to the critical point, the xenon is partly elastic: it can 'stretch' as well as flow. For STS-107, the hardware has been enhanced to determine if critical xenon is a shear-thinning fluid.

Berg, Robert F.↗

Critical Technologies for Hypersonic Vehicle Development

The objective of the present RTO/VKI Lecture Series is to focus attention on the technologies that play a critical role in hypersonic vehicle development. Owing to the demanding flight environment, successful development of hypersonic vehicles will require progress in understanding and surmounting specific technical problems in ground test, numerical simulation and flight test. Following the introductory remarks stating the objectives of the Lecture Series, a historical perspective on hypersonic vehicle development was given to establish the context and to highlight the emerging role of numerical simulation in vehicle design. A thorough review of ongoing vehicle development and flight-test programs completed the introductory portion. To provide sufficient background information, lectures were presented on compressible flows of perfect and imperfect gases, with special attention given to shock-wave boundary-layer interactions, laminar to turbulent transition, and high temperature gas properties. Advancements in computational modeling capabilities for hypersonic vehicle design were addressed, and the ongoing need for validation data, as well as a consistent validation methodology, where highlighted. Specific critical technology areas, including trajectory aerothermal environment definition, system trade studies, propulsion systems, aerodynamic control laws, flight instrumentation, thermal protection systems and flight experiments were all addressed in separate lectures given by international experts. Finally, the potential of plasma devices to ameliorate or avoid critical aerothermal phenomena was presented, providing a good example of the type of creative effort required for achieving progress. In all of the presentations the needs of future hypersonic flight vehicle development programs were addressed.

Research and Technology Organisation (RTO)↗

General Flow-Solver Code for Turbomachinery Applications

Phantom is a computer code intended primarily for real-fluid turbomachinery problems. It is based on Corsair, an ideal-gas turbomachinery code, developed by the same authors, which evolved from the ROTOR codes from NASA Ames. Phantom is applicable to real and ideal fluids, both compressible and incompressible, flowing at subsonic, transonic, and supersonic speeds. It utilizes structured, overset, O- and H-type zonal grids to discretize flow fields and represent relative motions of components. Values on grid boundaries are updated at each time step by bilinear interpolation from adjacent grids. Inviscid fluxes are calculated to third-order spatial accuracy using Roe s scheme. Viscous fluxes are calculated using second-order-accurate central differences. The code is second-order accurate in time. Turbulence is represented by a modified Baldwin-Lomax algebraic model. The code offers two options for determining properties of fluids: One is based on equations of state, thermodynamic departure functions, and corresponding state principles. The other, which is more efficient, is based on splines generated from tables of properties of real fluids. Phantom currently contains fluid-property routines for water, hydrogen, oxygen, nitrogen, kerosene, methane, and carbon monoxide as well as ideal gases.

Dorney, Daniel↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschens law that takes into account the flow of ambient gas past electrode surfaces. Paschens law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschens law when the mean velocity is zero. A new theoretical formulation of Paschens law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

discharge↗

Least-squares finite element methods for compressible Euler equations

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L-sq-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L-sq method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H1-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.

Jiang, Bo-Nan↗

Universality in the compressive behavior of solids

It was discovered that the isothermal equation of state for solids in compression is a simple, universal form. This single form accurately describes the pressure and bulk modulus as a function of volume for tonic, metallic, covalent, and rare gas solids.

Vinet, P.↗

Comment on 'The semiannual variation of great geomagnetic storms and the postshock Russell-McPherron effect preceding coronal mass ejecta' by N. U. Crooker, E. W. Cliver and B. T. Tsurutani

It is proposed by Crooker et al. (1992) that for a subgroup of great geomagnetic storms, for which the associated strong southward IMF (B(S)) fields reside in the postshock plasma, preceding the driver gas of coronal mass ejections, such strong B(S) fields result from a 'major increase in the Russell-McPherron polarity effect, through a systematic pattern of compression and draping' of the Archimedean field in the x-y plane. The critics test the scenario proposed by Crooker et al., namely, that the Russell-McPherron polarity effect is a major contribution to the semiannual variable of intense geomagnetic storms. It is found by the critics that in the cases studied there is little difference between the B(S) values as measured in geocentric solar ecliptic and geocentric solar magnetospheric coordinates, and it is concluded that the Russell-McPherron mechanism cannot explain by itself the seasonal dependence of intense storms, for which the variation is the largest. Crooker et al. present arguments to show that the combined preshock and postshock Russell-McPherron effect remains the sole cause of the semiannual variation of great storm occurrence.

Gonzalez, W. D.↗

Cryogenic Fuel Tank Draining Analysis Model

One of the technological challenges in designing advanced hypersonic aircraft and the next generation of spacecraft is developing reusable flight-weight cryogenic fuel tanks. As an aid in the design and analysis of these cryogenic tanks, a computational fluid dynamics (CFD) model has been developed specifically for the analysis of flow in a cryogenic fuel tank. This model employs the full set of Navier-Stokes equations, except that viscous dissipation is neglected in the energy equation. An explicit finite difference technique in two-dimensional generalized coordinates, approximated to second-order accuracy in both space and time is used. The stiffness resulting from the low Mach number is resolved by using artificial compressibility. The model simulates the transient, two-dimensional draining of a fuel tank cross section. To calculate the slosh wave dynamics the interface between the ullage gas and liquid fuel is modeled as a free surface. Then, experimental data for free convection inside a horizontal cylinder are compared with model results. Finally, cryogenic tank draining calculations are performed with three different wall heat fluxes to demonstrate the effect of wall heat flux on the internal tank flow field.

Greer, Donald↗

Development of Efficient Real-Fluid Model in Simulating Liquid Rocket Injector Flows

The characteristics of propellant mixing near the injector have a profound effect on the liquid rocket engine performance. However, the flow features near the injector of liquid rocket engines are extremely complicated, for example supercritical-pressure spray, turbulent mixing, and chemical reactions are present. Previously, a homogeneous spray approach with a real-fluid property model was developed to account for the compressibility and evaporation effects such that thermodynamics properties of a mixture at a wide range of pressures and temperatures can be properly calculated, including liquid-phase, gas- phase, two-phase, and dense fluid regions. The developed homogeneous spray model demonstrated a good success in simulating uni- element shear coaxial injector spray combustion flows. However, the real-fluid model suffered a computational deficiency when applied to a pressure-based computational fluid dynamics (CFD) code. The deficiency is caused by the pressure and enthalpy being the independent variables in the solution procedure of a pressure-based code, whereas the real-fluid model utilizes density and temperature as independent variables. The objective of the present research work is to improve the computational efficiency of the real-fluid property model in computing thermal properties. The proposed approach is called an efficient real-fluid model, and the improvement of computational efficiency is achieved by using a combination of a liquid species and a gaseous species to represent a real-fluid species.

Cheng, Gary↗

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗

Supersonic flow calculation using a Reynolds-stress and an eddy thermal diffusivity turbulence model

A second-order model for the velocity field and a two-equation model for the temperature field are used to calculate supersonic boundary layers assuming negligible real gas effects. The modeled equations are formulated on the basis of an incompressible assumption and then extended to supersonic flows by invoking Morkovin's hypothesis, which proposes that compressibility effects are completely accounted for by mean density variations alone. In order to calculate the near-wall flow accurately, correction functions are proposed to render the modeled equations asymptotically consistent with the behavior of the exact equations near a wall and, at the same time, display the proper dependence on the molecular Prandtl number. Thus formulated, the near-wall second order turbulence model for heat transfer is applicable to supersonic flows with different Prandtl numbers. The model is validated against flows with different Prandtl numbers and supersonic flows with free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.3. Among the flow cases considered, the momentum thickness Reynolds number varies from approximately 4,000 to approximately 21,000. Good correlation with measurements of mean velocity, temperature, and its variance is obtained. Discernible improvements in the law-of-the-wall are observed, especially in the range where the big-law applies.

Sommer, T. P.↗

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschen's law that takes into account the flow of ambient gas past electrode surfaces. Paschen's law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschen's law when the mean velocity is zero. A new theoretical formulation of Paschen's law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

electrostatic↗

Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)

The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.

Adam Johnson↗

Calculation of the Pressure Distribution on Bodies of Revolution in the Subsonic Flow of a Gas: Axially Symmetrical Flow - Part 1

The present report concerns a method of computing the velocity and pressure distributions on bodies of revolution in axially symmetrical flow in the subsonic range. The differential equation for the velocity potential Phi of a compressible fluid motion is linearized tn the conventional manner, and then put in the form Delta(Phi) = 0 by affine transformation. The quantity Phi represents the velocity potential of a fictitious incompressible flow, for which a constant superposition of sources by sections is secured by a method patterned after von Karman which must comply with the boundary condition delta(phi)/delta(n) = 0 at the originally specified contour. This requirement yields for the "pseudo-stream function" psi a differential equation which must be fulfilled for as many points on the contour as source lengths are assumed. In this manner, the problem of defining the still unknown source intensities is reduced to the solution of an inhomogeneous equation system. The pressure distribution is then determined with the aid of Bernoulli's equation and adiabatic equation of state. Lastly, the pressure distributions in compressible and incompressible medium are compared on a model problem.

Bilharz, Herbert↗