Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPRESSIBLE FLUID”

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 379 records · Page 21

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗

Euler Technology Assessment for Preliminary Aircraft Design: Compressibility Predictions by Employing the Cartesian Unstructured Grid SPLITFLOW Code

The objective of the second phase of the Euler Technology Assessment program was to evaluate the ability of Euler computational fluid dynamics codes to predict compressible flow effects over a generic fighter wind tunnel model. This portion of the study was conducted by Lockheed Martin Tactical Aircraft Systems, using an in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaption of the volume grid during the solution to resolve high-gradient regions. The SPLITFLOW code predictions of configuration forces and moments are shown to be adequate for preliminary design, including predictions of sideslip effects and the effects of geometry variations at low and high angles-of-attack. The transonic pressure prediction capabilities of SPLITFLOW are shown to be improved over subsonic comparisons. The time required to generate the results from initial surface data is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.

Finley, Dennis B.↗

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗

Sphere based fluid systems

Systems are described for using multiple closely-packed spheres. In one system for passing fluid, a multiplicity of spheres lie within a container, with all of the spheres having the same outside diameter and with the spheres being closely nested in one another to create multiple interstitial passages of a known size and configuration and smooth walls. The container has an inlet and outlet for passing fluid through the interstitial passages formed between the nested spheres. The small interstitial passages can be used to filter out material, especially biological material such as cells in a fluid, where the cells can be easily destroyed if passed across sharp edges. The outer surface of the spheres can contain a material that absorbs a constitutent in the flowing fluid, such as a particular contamination gas, or can contain a catalyst to chemically react the fluid passing therethrough, the use of multiple small spheres assuring a large area of contact of these surfaces of the spheres with the fluid. In a system for storing and releasing a fluid such as hydrogen as a fuel, the spheres can include a hollow shell containing the fluid to be stored, and located within a compressable container that can be compressed to break the shells and release the stored fluid.

Elleman, Daniel D.↗

Method and Apparatus for Predicting Unsteady Pressure and Flow Rate Distribution in a Fluid Network

A method and apparatus for analyzing steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics, external body forces such as gravity and centrifugal force and conjugate heat transfer. In some embodiments, a graphical user interface provides for the interactive development of a fluid network simulation having nodes and branches. In some embodiments, mass, energy, and specific conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. In some embodiments, contained herein are data objects for computing thermodynamic and thermophysical properties for fluids. In some embodiments, the systems of equations describing the fluid network are solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods.

Majumdar, Alok K.↗

Measurement of the Specific Heat Using a Gravity Cancellation Approach

The specific heat at constant volume C(sob V) of a simple fluid diverges near its liquid-vapor critical point. However, gravity-induced density stratification due to the divergence of isothermal susceptibility hinders the direct comparison of the experimental data with the predictions of renormalization group theory. In the past, a microgravity environment has been considered essential to eliminate the density stratification. We propose to perform specific heat measurements of He-3 on the ground using a method to cancel the density stratification. A He-3 fluid layer will be heated from below, using the thermal expansion of the fluid to cancel the hydrostatic compression. A 6% density stratification at a reduced temperature of 10(exp -5) can be cancelled to better than 0.1% with a steady 1.7 micro K temperature difference across a 0.05 cm thick fluid layer. A conventional AC calorimetry technique will be used to determine the heat capacity. The minimized bulk density stratification with a relaxation time 6500 sec at a reduced temperature of 10(exp -5) will stay unchanged during 1 Hz AC heating. The smear of the specific heat divergence due to the temperature difference across the cell is about 0.1% at a reduced temperature of 10(exp -6). The combination of using High Resolution Thermometry with a 0.5 n K temperature resolution in the AC technique and the cancellation of the density stratification will enable C(sub V) to be measured down to a reduced temperature of 10(exp -6) with less than a 1% systematic error.

Zhong, Fang↗

Computation of Pressurized Gas Bearings Using CE/SE Method

The space-time conservation element and solution element (CE/SE) method is extended to compute compressible viscous flows in pressurized thin fluid films. This numerical scheme has previously been used successfully to solve a wide variety of compressible flow problems, including flows with large and small discontinuities. In this paper, the method is applied to calculate the pressure distribution in a hybrid gas journal bearing. The formulation of the problem is presented, including the modeling of the feeding system. the numerical results obtained are compared with experimental data. Good agreement between the computed results and the test data were obtained, and thus validate the CE/SE method to solve such problems.

Cioc, Sorin↗

Transport Coefficients in Rotating Weakly Compressible Turbulence

Analytical studies of compressible turbulence have found that compressible velocity fluctuations create both effective fluid transport properties and an effective equation of state. This paper investigates the effects of rotation on compressible turbulence. It is shown that rotation modifies the transport properties of compressible turbulence by replacing the turbulence time scale by a rotational time scale, much as rotation modifies the transport properties of incompressible turbulence. But thermal equilibrium properties are modified in a more complex manner. Two regimes are possible: one dominated by incompressible fluctuations, in which the sound speed is modified as it is in non-rotating turbulence, and a rotation dominated regime in which the sound speed enhancement is rotation dependent. The dimensionless parameter which discriminates between regimes is identified. In general, rotation is found to suppress the effects of compressibility. A novel feature of the present analysis is the use of a non-Kolmogorov steady state as the reference state of turbulence. introduction of such steady states expands the power and utility of analytical turbulence closures to a wider range of problems.

Rubinstein, Robert↗

Motion of particles with inertia in a compressible free shear layer

The effects of the inertia of a particle on its flow-tracking accuracy and particle dispersion are studied using direct numerical simulations of 2D compressible free shear layers in convective Mach number (Mc) range of 0.2 to 0.6. The results show that particle response is well characterized by tau, the ratio of particle response time to the flow time scales (Stokes' number). The slip between particle and fluid imposes a fundamental limit on the accuracy of optical measurements such as LDV and PIV. The error is found to grow like tau up to tau = 1 and taper off at higher tau. For tau = 0.2 the error is about 2 percent. In the flow visualizations based on Mie scattering, particles with tau more than 0.05 are found to grossly misrepresent the flow features. These errors are quantified by calculating the dispersion of particles relative to the fluid. Overall, the effect of compressibility does not seem to be significant on the motion of particles in the range of Mc considered here.

Samimy, M.↗

Heat transfer head for a Stirling cycle machine

A common heat acceptor is provided between opposed displacers in a Stirling cycle machine. It includes two sets of open channels in separate fluid communications with the expansion spaces of the receptive cyclinders. The channels confine movement of working fluid in separate paths that extend between the expansion space of one cylinder and the compression space of the other. The method for operating the machine involves alternatively directing working fluid from the expansion space of each cylinder in a fluid path leading to the compression space of the other cylinder and from the compression space of each cylinder in a fluid path leading to the expansion space of the other cylinder.

Emigh, Stuart G.↗

Simulation of tail buffet using delta wing-vertical tail configuration

Computational simulation of the vertical tail buffet problem is accomplished using a delta wing-vertical tail configuration. Flow conditions are selected such that the wing primary-vortex cores experience vortex breakdown and the resulting flow interacts with the vertical tail. This multidisciplinary problem is solved successively using three sets of equations for the fluid flow, aeroelastic deflections and grid displacements. For the fluid dynamics part, the unsteady, compressible, full Navier-Stokes equations are solved accurately in time using an implicit, upwind, flux-difference splitting, finite-volume scheme. For the aeroelastic part, the aeroelastic equation for bending vibrations is solved accurately in time using the Galerkin method and the four-stage Runge-Kutta scheme. The grid for the fluid dynamics computations is updated every few time steps using a third set of interpolation equations. The computational application includes a delta wing of aspect ratio 1 and a rectangular vertical tail of aspect ratio 2, which is placed at 0.5 root-chord length downstream of the wing trailing edge. The wing angle of attack is 35 deg and the flow Mach number and Reynolds number are 0.4 and 10,000, respectively.

Kandil, Osama A.↗

Investigation of Compressibility Effect for Aeropropulsive Shear Flows

Rocket Based Combined Cycle (RBCC) engines operate within a wide range of Mach numbers and altitudes. Fundamental fluid dynamic mechanisms involve complex choking, mass entrainment, stream mixing and wall interactions. The Propulsion Research Center at the University of Alabama in Huntsville is involved in an on- going experimental and numerical modeling study of non-axisymmetric ejector-based combined cycle propulsion systems. This paper attempts to address the modeling issues related to mixing, shear layer/wall interaction in a supersonic Strutjet/ejector flow field. Reynolds Averaged Navier-Stokes (RANS) solutions incorporating turbulence models are sought and compared to experimental measurements to characterize detailed flow dynamics. The effect of compressibility on fluids mixing and wall interactions were investigated using an existing CFD methodology. The compressibility correction to conventional incompressible two- equation models is found to be necessary for the supersonic mixing aspect of the ejector flows based on 2-D simulation results. 3-D strut-base flows involving flow separations were also investigated.

Balasubramanyam, M. S.↗

Numerical comparison of convective heat transfer augmentation devices used in cooling channels of hypersonic vehicles

Hypersonic vehicles are exposed to extreme thermal conditions compared to subsonic aircraft; therefore, some level of thermal management is required to protect the materials used. Normally, hypersonic vehicles experience the highest temperatures in the nozzle throat, and aircraft and propulsion system leading edges. Convective heat transfer augmentation techniques can be used in the thermal management system to increase heat transfer of the cooling channels in those areas. The techniques studied in this report are pin-fin, offset-fin, ribbed and straight roughened channel. A smooth straight channel is used as the baseline for comparing the techniques. SINDA '85, a lumped parameter finite difference thermal analyzer, is used to model the channels. Subroutines are added to model the fluid flow assuming steady one dimensional compressible flow with heat addition and friction. Correlations for convective heat transfer and friction are used in conjunction with the fluid flow analysis mentioned. As expected, the pin-fin arrangement has the highest heat transfer coefficient and the largest pressure drop. All the other devices fall in between the pin-fin and smooth straight channel. The selection of the best heat augmentation method depends on the design requirements. A good approach may be a channel using a combination of the techniques. For instance, several rows of pin-fins may be located at the region of highest heat flux, surrounded by some of the other techniques. Thus, the heat transfer coefficient is maximized at the region of highest heat flux while the pressure drop is not excessive.

Maldonado, Jaime J.↗

New Type of the Interface Evolution in the Richtmyer-Meshkov Instability

We performed systematic theoretical and numerical studies of the nonlinear large-scale coherent dynamics in the Richtmyer-Meshkov instability for fluids with contrast densities. Our simulations modeled the interface dynamics for compressible and viscous uids. For a two-fluid system we observed that in the nonlinear regime of the instability the bubble velocity decays and its surface attens, and the attening is accompanied by slight oscillations. We found the theoretical solution for the system of conservation laws, describing the principal influence of the density ratio on the motion of the nonlinear bubble. The solution has no adjustable parameters, and shows that the attening of the bubble front is a distinct property universal for all values of the density ratio. This property follows from the fact that the RM bubbles decelerate. The theoretical and numerical results validate each other, describe the new type of the bubble front evolution in RMI, and identify the bubble curvature as important and sensitive diagnostic parameter.

Abarzhi, S. I.↗

Actively suspended counter-rotating machine

A counter-rotating machine, such as a positive displacement pump having a pair of meshed, non-contacting helical screws (10,12), subjects its rotating members to axial and radial thrust forces when used for such purposes as compression of liquid or gaseous phase fluids while transporting them through a pump cavity (11,13). Each helical screw (10,12) has a shaft (17,17') which is actively suspended at opposite ends (11a,11b) of the pump cavity by a servo-controlled magnetic bearing assembly (19) and a servo-controlled rotary drive motor (20). Both bearing assemblies and drive motors are mounted on the outside of the pump cavity (11,13). Opto-electric angular position sensors (250) provide synchronization between radial orientation of the drive motors. The bearing assemblies and drive motors conjugately provide axial stabilization and radial centering of the helical screws during volumetric compression of aspirated liquid or gaseous phase fluids.

Studer, Philip A.↗

Effects of leading-edge flap oscillation on unsteady delta wing flow and rock control

The isolated and interdisciplinary problems of unsteady fluid dynamics and rigid-body dynamics and control of delta wings with and without leading-edge flap oscillation are considered. For the fluid dynamics problem, the unsteady, compressible, thin-layer Navier-Stokes (NS) equations, which are written relative to a moving frame of reference, are solved along with the unsteady, linearized, Navier-displacement (ND) equations. The NS equations are solved for the flowfield using an implicit finite-volume scheme. The ND equations are solved for the grid deformation, if the leading-edge flaps oscillate, using an ADI scheme. For the dynamics and control problem, the Euler equation of rigid-body rolling motion for a wing and its flaps are solved interactively with the fluid dynamics equations for the wing-rock motion and subsequently for its control. A four-stage Runge-Kutta scheme is used to explicitly integrate the dynamics equation.

Kandil, Osama A.↗

Compressibility effects on rotor forces in the leakage path between a shrouded pump impeller and its housing

A modified approach to Childs' previous work on fluid-structure interaction forces in the leakage path between an impeller shroud and its housing is presented in this paper. Three governing equations consisting of continuity, path-momentum, and circumferential-momentum equations were developed to describe the leakage path inside a pump impeller. Radial displacement perturbations were used to solve for radial and circumferential force coefficients. In addition, impeller-discharge pressure disturbances were used to obtain pressure oscillation responses due to precessing impeller pressure wave pattern. Childs' model was modified from an incompressible model to a compressible barotropic-fluid model (the density of the working fluid is a function of the pressure and a constant temperature only). Results obtained from this model yielded interaction forces for radial and circumferential force coefficients. Radial and circumferential forces define reaction forces within the impeller leakage path. An acoustic model for the same leakage path was also developed. The convective, Coriolis, and centrifugal acceleration terms are removed from the compressible model to obtain the acoustics model. A solution due to impeller discharge pressure disturbances model was also developed for the compressible and acoustics models. The results from these modifications are used to determine what effects additional perturbation terms in the compressible model have on the acoustic model. The results show that the additional fluid mechanics terms in the compressible model cause resonances (peaks) in the force coefficient response curves. However, these peaks only occurred at high values of inlet circumferential velocity ratios greater than 0.7. The peak pressure oscillation was shown to occur at the wearing ring seal. Introduction of impeller discharge disturbances with n = 11 diametral nodes showed that maximum peak pressure oscillations occurred at nondimensional precession frequencies of f = 6.4 and f = 7.8 for this particular pump. Bolleter's results suggest that for peak pressure oscillations to occur at the wearing ring seal, the nondimensional excitation frequency should be on the order of f = 2.182 for n = 11. The resonances found in this research do not match the excitation frequencies predicted by Bolleter. At the predicted peak excitation frequencies given by Bolleter, the compressible model shows an attenuation of the pressure oscillations at the seal exit. The compressibility of the fluid does not have a significant influence on the model at low values of nondimensional excitation frequency. At high values of nondimensional frequency, the effects of compressibility become more significant. For the acoustic analysis, the convective, Coriolis, and centrifugal acceleration terms do affect the results to a limited extent for precession excitation and to a large extent for a pressure excitation when the fluid operates at relatively high Mach numbers.

Cao, Nhai The↗