Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 541 records · Page 30

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, C. M.↗

dc power control for a liquid-fed resistojet

A simple breadboard power controller was designed and demonstrated for a new liquid-fed water resistojet. The 1-piece laboratory model thruster has an integrated vaporizer/superheater using a single heating element. Heater temperature was maintained at or near a preset reference value with the closed loop controller providing pulse width modulated (PWM) dc power into the thruster heater. A combined thruster, temperature readout, PWM transfer function was experimentally determined. This transfer function was used to design a proportional plus integral controller that demonstrated zero steady state error, conservative stability margins and adequate transient response to step changes in propellant flow rate, input voltage and temperature reference. Initial turn-on temperature overshoot from room temperature to a 650 C setpoint was 80 C. In addition, EMI was alleviated by reducing heater dI/dt and dV/dt using a simple diode-inductor-capacitor network. Based on limited initial tests, thruster preheat with no propellant flow was necessary to achieve stable system operation during startup. Breadboard power efficiency was 99 percent at 1 kW, and component mass was 0.4 kg excluding the power loss and mass of an input filter required for spacecraft integration.

Gruber, Robert P.↗

Integration of Biological, Physical/Chemical and Energy Efficient Systems in the CELSS Antarctic Analog: Performance of Prototype Systems and Issues for Life Support

The Controlled Ecological Life Support System (CELSS) Antarctic Analog Project (CAAP) is a joint endeavor between the National Science Foundation, Office of Polar Programs (NSF-OPP), and the National Aeronautics and Space Administration (NASA). The fundamental objective is to develop, deploy, and operate a testbed of advanced life support technologies at the Amundsen-Scott South Pole Station that enable the objectives of both the NSF and NASA. The functions of food production, water purification, and waste treatment, recycle, and reduction provided by CAAP will improve the quality of life for the South Pole inhabitants, reduce logistics dependence, enhance safety, and minimize environmental impacts associated with human presence on the polar plateau. Because of the analogous technical, scientific, and mission features with Planetary missions, such as a mission to Mars, CAAP provides NASA with a method for validating technologies and overall approaches to supporting humans. Prototype systems for waste treatment, water recycle, resource recovery and crop production are being evaluated in a testbed at Ames Research Center. The combined performance of these biological and physical/chemical systems as an integrated function in support of the human habitat will be discussed. Overall system performance will be emphasized. The effectiveness and efficiency of component technologies will be discussed in the context of energy and mass flow within the system and contribution to achieving a mass and energy conservative system. Critical to the discussion are interfaces with habitat functions outside of the closed-loop life support: the ability of the system to satisfy the life support requirements of the habitat and the ability to define input requirements. The significance of analog functions in relation to future Mars habitats will be discussed.

Bubenheim, David L.↗

Development of A Maximum Anti-Slosh Baffle Pressure Load Model

The sloshing of propellants can affect the stability of a spacecraft and the integrity of the tank structure. Undesirable sloshing can be controlled by the addition of anti-slosh baffles, and the spacing and configuration of baffles are driven by damping requirements. The structural design of the baffle is determined after consideration of many factors, such as the strength and rigidity needed to support the baffle for its lifetime. Therefore, knowledge of distributed pressure loading is important for detailed structural design. In addition, the resultant force and moment produced by the distributed pressure are of direct importance to the design of a vehicle’s control system. Previous experimental investigations have been conducted to determine the liquid pressure loads and slosh damping associated with a rigid ring baffle. The results suggested that when the nondimensional velocity parameter is larger than 1.0, the theories agree with the test. However, when the velocity parameter is less than 1.0, all theories are nonconservative and under-predict the pressure loads. The present study has derived a maximum pressure load on the slosh baffle based on the energy conservation principle. It is verified from CFD that pressure in the slosh flow field can be decomposed into static and transient components. The CFD results confirm that there is a phase shift in pressure across the baffle, which depends on the fluid damping. Higher damping leads to a higher phase shift. The CFD investigation further verifies the proposed theory: the maximum pressure load occurs when the phase shift is 90 degrees. A comparison of the present computational results to the previous comprehensive experimental data validates the maximum pressure theory. When the baffle is submerged, the maximum pressure theory envelopes all the experimental data points.

CFD↗

Investigation of Extreme Leading-Edge Roughness on Thick Low-Drag Airfoils to Indicate Those Critical to Separation

Several airfoils, Including a conventional NACA 23021 and some low-drag airfoils for which the thickness had been increased to the point that they were considered doubtfully conservative with respect to separation, were investigated as smooth airfoils and after the application of a standard roughness. The results show some of the airfoils to be critical to separation resulting from such flow disturbances. It is concluded, pending the further investigation of separation difficulties, that airfoil sections falling definitely within the conservative range should be used.

AIRFOILS - NACA 65,3-418 - CHAR & MOMENTS - LOW-DR↗

Calculation of flow about two-dimensional bodies by means of the velocity-vorticity formulation on a staggered grid

A method for calculating the incompressible viscous flow about two-dimensional bodies, utilizing the velocity-vorticity form of the Navier-Stokes equations using a staggered-grid formulation is presented. The solution is obtained by employing an alternative-direction implicit method for the solution of the block tridiagonal matrix resulting from the finite-difference representation of the governing equations. The boundary vorticity and the conservation of mass are calculated implicitly as a part of the solution. The mass conservation is calculated to machine zero for the duration of the computation. Calculations for the flow about a circular cylinder, a 2-pct thick flat plate at 90-deg incidence, an elliptic cylinder at 45-deg incidence, and a NACA 0012, with and without a deflected flap, at - 90-deg incidence are performed and compared with the results of other numerical investigations.

Stremel, Paul M.↗

Nonequilibrium thermo-chemical calculations using a diagonal implicit scheme

A recently developed computer program for hypersonic vehicle flow analysis is described. The program uses a diagonal implicit algorithm to solve the equations of viscous flow for a gas in thermochemical nonequilibrium. The diagonal scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The program uses multiple zones of grids patched together and includes radiation wall and rarefied gas boundary conditions. Solutions are presented for hypersonic flows of air and hydrogen air mixtures.

Imlay, Scott T.↗

Three dimensional hypersonic flow simulations with the CSCM implicit upwind Navier-Stokes method

The Bardina and Lombard (1985) three-dimensional CSCM Navier-Stokes method is presently extended to the simulation of complex hypersonic reentry vehicle external flows at angle of attack. The robust stability of the method derives from the combination of conservative implicit upwind flux difference splitting with a three-dimensional diagonally-dominant approximate factorization and relaxation scheme and characteristic-based implicit boundary approximations. The method's efficiency derives from an implicit symmetric Gauss-Seidel 'method of planes' relaxation scheme with alternating directional space marching sweeps along the flow coordinate direction.

Bardina, Jorge↗

A New Test Method for Material Flammability Assessment in Microgravity and Extraterrestrial Environments

The standard oxygen consumption (cone) calorimeter (described in ASTM E 1354 and NASA STD 6001 Test 2) is modified to provide a bench-scale test environment that simulates the low velocity buoyant or ventilation flow generated by or around a burning surface in a spacecraft or extraterrestrial gravity level. The Equivalent Low Stretch Apparatus (ELSA) uses an inverted cone geometry with the sample burning in a ceiling fire (stagnation flow) configuration. For a fixed radiant flux, ignition delay times for characterization material PMMA are shown to decrease by a factor of three at low stretch, demonstrating that ignition delay times determined from normal cone tests significantly underestimate the risk in microgravity. The critical heat flux for ignition is found to be lowered at low stretch as the convective cooling is reduced. At the limit of no stretch, any heat flux that exceeds the surface radiative loss at the surface ignition temperature is sufficient for ignition. Regression rates for PMMA increase with heat flux and stretch rate, but regression rates are much more sensitive to heat flux at the low stretch rates, where a modest increase in heat flux of 25 kW/m2 increases the burning rates by an order of magnitude. The global equivalence ratio of these flames is very fuel rich, and the quantity of CO produced in this configuration is significantly higher than standard cone tests. These results [2] demonstrate the ELSA apparatus allows us to conduct normal gravity experiments that accurately and quantifiably evaluate a material s flammability characteristics in the real-use environment of spacecraft or extra-terrestrial gravitational acceleration. These results also demonstrate that current NASA STD 6001 Test 2 (standard cone) is not conservative since it evaluates materials flammability with a much higher inherent buoyant convective flow.

Olson, S. L.↗

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented.

Johnson, G. M.↗

Review of research on wave mean flow interactions using FGGE data

Comparison is made between the conventional Eulerian diagnostics of wave mean flow interactions in the atmosphere and the diagnostics introduced by Andrews and McIntyre (1976). In the latter, the meridional circulation is represented by a residual circulation, and the eddy torque in the zonal momentum equation by the Eliassen-Palm (1960) flux divergence (or equivalently, by the poleward eddy flux of geostrophic potential vorticity). It is found that, while the new diagnostics are useful in studying stratospheric warmings and transports of conservative quantities, the conventional diagnostics are more useful in studying global energetics and wave mean flow interactions in the troposphere and lower stratosphere. Case studies diagnosed with FGGE IIIb data are presented to illustrate this point. Kung and Tanaka's (1983) comparisons of the atmospheric energy cycle for SOP-1 and SOP-2, estimated with the use of the ECMWF and GFDL FGGE IIIb data bases, are used to assess the uncertainty in calculations of energetics based on the use of FGGE IIIb data. These results demonstrate the high degree of model dependence in the FGGE IIb data sets.

Pfeffer, R. L.↗

A numerical study of flow phenomena within a vortex sink rate sensor

A numerical investigation of a viscous incompressible fluid with constant properties within a vortex sink rate sensor is undertaken. A computer program was developed, for obtaining the numerical solution of the Navier-Stokes equations, assuming laminar flow, having general prescribed inlet conditions and axisymmetric boundary conditions. Computational results for various viscous flows and assorted boundary conditions have been obtained. In the vortex chamber, the radial momentum is conserved, whereas the axial momentum is negligible. In the sink region, the radial and tangential velocities, as well as the tangential vorticities, have significantly high values. In the sink tube, the radial momentum is negligible. The radial momentum within the vortex chamber is converted to the axial momentum in the sink tube. A change in the flow rate had an appreciable effect within the sensor and particularly near the sink tube entrance. As the flow rate was increased, both the tangential velocity and tangential vorticity increased rapidly. At the higher flow rates, vortices were produced at the corner of the entrance section of the sink tube.

Patel, D. K.↗

Development of prediction techniques for multi-jet thermal ground flow fields and fountain formation

Avoiding detrimental ground interaction is important for practical V/STOL aircraft. This paper reports recent developments in a numerical method for estimating thermal ground footprints. Upwash and fountain formation for arbitrarily oriented jet arrangements is predicted. Flow asymmetry due to roll, pitch, differential thrust or ground inclination is included. The prediction methodology uses simple inviscid relations for energy and momentum conservation along with an empirical entrainment law, applied in independent sectors of the wall jet and upwash. Asymmetrical stagnation line prediction is compared with experiment. Detailed flow measurements for a three-jet interaction are also presented.

Siclari, M. J.↗

Modified ICED-ALE method for astrogeophysical plasma flows

The Implicit-Continuous-Eulerian-Difference Mesh-Arbitrary-Lagrangian-Eulerian (ICED-ALE) algorithm of Brackbill and Pracht (1973) is modified for the study of astrophysical plasma flows in which dynamical effects are important. In the present study the general-energy-conservation law is directly applied to the iteration process, with the total (kinetic-, specific-internal-, and magnetic-) energy density being obtained implicitly at the end of the process. An example is computed in which the convergence speed of the latter method is substantially faster than that of the former. The initializing calculation, or explicit phase, in which the electric current density, magnetic diffusion of field, energy augmentation, and zero-order approximation of the flow velocity are given, is described. Consideration is given to the iteration process, or implicit phase, from which the exact Lagrangian solution for energy density, velocity, and a magnetic field is obtained.

Wu, S. T.↗

Application of viscous-inviscid interaction methods to transonic turbulent flows

Two different viscous-inviscid interaction schemes were developed for the analysis of steady, turbulent, transonic, separated flows over axisymmetric bodies. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity approach. In the semi-inverse interaction scheme, the viscous and inviscid equations are solved in an explicitly separate manner and the displacement thickness distribution is iteratively updated by a simple coupling algorithm. In the simultaneous interaction method, local solutions of viscous and inviscid equations are treated simultaneously, and the displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure. The inviscid flow region is described by a direct finite-difference solution of a velocity potential equation in conservative form. The potential equation is solved on a numerically generated mesh by an approximate factorization (AF2) scheme in the semi-inverse interaction method and by a successive line overrelaxation (SLOR) scheme in the simultaneous interaction method. The boundary-layer equations are used for the viscous flow region. The continuity and momentum equations are solved inversely in a coupled manner using a fully implicit finite-difference scheme.

Lee, D.↗

Vector processor algorithms for transonic flow calculations

This paper discusses a number of algorithms for solving the transonic full-potential equation in conservative form on a vector computer, such as the CDC STAR-100 or the CRAY-1. Recent research with the 'artificial density' method for transonics has led to development of some new iteration schemes which take advantage of vector-computer architecture without suffering significant loss of convergence rate. Several of these more promising schemes are described and 2-D and 3-D results are shown comparing the computational rates on the STAR and CRAY vector computers, and the CYBER-175 serial computer. Schemes included are: (1) Checkerboard SOR, (2) Checkerboard Leapfrog, (3) odd-even vertical line SOR, and (4) odd-even horizontal line SOR.

South, J. C., Jr.↗

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

A time-accurate algorithm for chemical non-equilibrium viscous flows at all speeds

A time-accurate, coupled solution procedure is described for the chemical nonequilibrium Navier-Stokes equations over a wide range of Mach numbers. This method employs the strong conservation form of the governing equations, but uses primitive variables as unknowns. Real gas properties and equilibrium chemistry are considered. Numerical tests include steady convergent-divergent nozzle flows with air dissociation/recombination chemistry, dump combustor flows with n-pentane-air chemistry, nonreacting flow in a model double annular combustor, and nonreacting unsteady driven cavity flows. Numerical results for both the steady and unsteady flows demonstrate the efficiency and robustness of the present algorithm for Mach numbers ranging from the incompressible limit to supersonic speeds.

Shuen, J.-S.↗