Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reservoir flow model”

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

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

A basic code for the prediction of transient three-dimensional turbulent flowfields

A primitive pressure-velocity finite difference code has been developed to predict transient three-dimensional turbulent flow. The code is a simplified yet effective prediction procedure for use by persons with little experience in computational fluid dynamics, and into which user-oriented complexities can easily be added. The method is based on the transient two-dimensional Los Alamos SOLA prediction technique for laminar flows. Turbulence is simulated by means of a two-equation turbulence model; species diffusion and buoyancy are also included. Two applications of the code are presented to local destratification near the release structure of a reservoir and to the deflection of a jet entering normally into a uniform cross-flow. Predicted results exhibit good agreement with experimental data, showing that a useful characterization of fully three-dimensional flows is now available.

Busnaina, A. A.↗

Nonequilibrium Flow Expansion Experiment Around a Blunted Cone

Eleven tests were carried out in DLR's high enthalpy tunnel in Goettingen (HEG), Germany, for reservoir conditions ranging from 10 to 23 MJ/kg and a free stream Mach number of approximately 10. A blunted cone model with a cylindrical afterbody (sting) was investigated. To obtain information on the influence of defined parameters on the body back flow, especially of the reacting gas, the heat transfer rate along the body contour was measured with fast response surface thermocouples on the forebody and sensitive thin film heat transfer gauges on the base and sting of the model. Flow visualization with a holographic interferometry system was provided. Tests are described, and a preliminary interpretation of the observed flow effects are given.

Kastell, Dirk↗

Propulsion System Development for the Iodine Satellite (iSAT) Demonstration Mission

The development and testing of a 200-W iodine-fed Hall thruster propulsion system that will be flown on a 12-U CubeSat is described. The switch in propellant from more traditional xenon gas to solid iodine yields the advantage of high density, low pressure propellant storage but introduces new requirements that must be addressed in the design and operation of the propulsion system. The thruster materials have been modified from a previously-flown xenon Hall thruster to make it compatible with iodine vapor. The cathode incorporated into this design additionally requires little or no heating to initiate the discharge, reducing the power needed to start the thruster. The feed system produces iodine vapor in the propellant reservoir through sublimation and then controls the flow to the anode and cathode of the thruster using a pair of proportional flow control valves. The propellant feeding process is controlled by the power processing unit, with feedback control on the anode flow rate provided through a measure of the thruster discharge current. Thermal modeling indicates that it may be difficult to sufficiently heat the iodine if it loses contact with the propellant reservoir walls, serving to motivate future testing of that scenario to verify the modeling result and develop potential mitigation strategies. Preliminary, short-duration materials testing has thus-far indicated that several materials may be acceptable for prolonged contact with iodine vapor, motivating longer-duration testing. A propellant loading procedure is presented that aims to minimize the contaminants in the feed system and propellant reservoir. Finally, an 80-hour duration test being performed to gain experience operating the thruster over long durations and multiple restarts is discussed.

Polzin, Kurt A.↗

A comparison of economic evaluation models as applied to geothermal energy technology

Several cost estimation and financial cash flow models have been applied to a series of geothermal case studies. In order to draw conclusions about relative performance and applicability of these models to geothermal projects, the consistency of results was assessed. The model outputs of principal interest in this study were net present value, internal rate of return, or levelized breakeven price. The models used were VENVAL, a venture analysis model; the Geothermal Probabilistic Cost Model (GPC Model); the Alternative Power Systems Economic Analysis Model (APSEAM); the Geothermal Loan Guarantee Cash Flow Model (GCFM); and the GEOCOST and GEOCITY geothermal models. The case studies to which the models were applied include a geothermal reservoir at Heber, CA; a geothermal eletric power plant to be located at the Heber site; an alcohol fuels production facility to be built at Raft River, ID; and a direct-use, district heating system in Susanville, CA.

Ziman, G. M.↗

Experimental and Computational Studies of the Flow Over a Sting Mounted Planetary Probe Configuration

This paper summarizes the results of a series of experimental studies in the LENS shock tunnel and computations with DSMC and Navier Stokes codes which have been made to examine the aerothermal and flowfield characteristics of the flow over a sting-supported planetary probe configuration in hypervelocity air and nitrogen flows. The experimental program was conducted in the LENS hypervelocity shock tunnel at total enthalpies of 5and 10 MJkg for a range of reservoir pressure conditions from 70 to 500 bars. Heat transfer and pressure measurements were made on the front and rear face of the probe and along the supporting sting. High-speed and single shot schlieren photography were also employed to examine the flow over the model and the time to establish the flow in the base recirculation region. Predictions of the flowfield characteristics and the distributions of heat transfer and pressure were made with DSMC codes for rarefied flow conditions and with the Navier-Stokes solvers for the higher pressure conditions where the flows were assumed to be laminar. Analysis of the time history records from the heat transfer and pressure instrumentation on the face of the probe and in the base region indicated that the base flow was fully established in under 4 milliseconds from flow initiation or between 35 and 50 flow lengths based on base height. The measurements made in three different tunnel entries with two models of identical geometries but with different instrumentation packages, one prepared by NASA Langley and the second prepared by CUBRC, demonstrated good agreement between heat transfer measurements made with two different types of thin film and coaxial gage instrumentation. The measurements of heat transfer and pressure to the front face of the probe were in good agreement with theoretical predictions from both the DSMC and Navier Stokes codes. For the measurements made in low density flows, computations with the DSMC code were found to compare well with the pressure and heat transfer measurements on the sting, although the computed heat transfer rates in the recirculation region did not exhibit the same characteristics as the measurements. For the 10MJkg and 500 bar reservoir match point condition, the measurements and heat transfer along the sting from the first group of studies were in agreement with the Navier Stokes solutions for laminar conditions. A similar set of measurements made in later tests where the model was moved to a slightly different position in the test section indicated that the boundary layer in the reattachment compression region was close to transition or transitional where small changes in the test environment can result in larger than laminar heating rates. The maximum heating coefficients on the sting observed in the present studies was a small fraction of similar measurements obtained at nominally the same conditions in the HEG shock tunnel, where it is possible for transition to occur in the base flow, and in the low enthalpy studies conducted in the NASA Langley high Reynolds number Mach 10 tunnel where the base flow was shown to be turbulent. While the hybrid Navier- StokedDMSC calculations by Gochberg et al. (Reference 1) suggested that employing the Navier- Stokes calculations for the entire flowfield could be seriously in error in the base region for the 10 MJkg, 500 bar test case, similar calculations performed by Cornell, presented here, do not.

Holden, Michael S.↗

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗

Fundamentals of chemistry modeling applicable to a vectorized particle simulation

This paper describes the fundamentals of extending the vectorized particle simulation method derived by Baganoff and McDonald (1990), McDonald and Baganoff (1988), and McDonald (1989) for modeling chemically reacting flows. Details of reaction mechanics per reaction are presented, with particular attention given to the quantum nature of the vibrational mode. The models of reactive flows developed here were verified through a simulation of a superheated diatomic gas relaxing thermochemically to equilibrium in a reservoir.

Haas, Brian L.↗

Transient Analysis of a Magnetic Heat Pump

An experimental heat pump that uses a rare earth element as the refrigerant is modeled using NASTRAN. The refrigerant is a ferromagnetic metal whose temperature rises when a magnetic field is applied and falls when the magnetic field is removed. The heat pump is used as a refrigerator to remove heat from a reservoir and discharge it through a heat exchanger. In the NASTRAN model the components modeled are represented by one-dimensional ROD elements. Heat flow in the solids and fluid are analyzed. The problem is mildly nonlinear since the heat capacity of the refrigerant is temperature-dependent. One simulation run consists of a series of transient analyses, each representing one stroke of the heat pump. An auxiliary program was written that uses the results of one NASTRAN analysis to generate data for the next NASTRAN analysis.

Schroeder, E. A.↗

Simulation Software

Various NASA Small Business Innovation Research grants from Marshall Space Flight Center, Langley Research Center and Ames Research Center were used to develop the 'kernel' of COMCO's modeling and simulation software, the PHLEX finite element code. NASA needed it to model designs of flight vehicles; one of many customized commercial applications is UNISIM, a PHLEX-based code for analyzing underground flows in oil reservoirs for Texaco, Inc. COMCO's products simulate a computational mechanics problem, estimate the solution's error and produce the optimal hp-adapted mesh for the accuracy the user chooses. The system is also used as a research or training tool in universities and in mechanical design in industrial corporations.

Source record↗

Thermal and dynamical regimes of single- and two-phase magmatic flow in dikes

The coupling between thermal and dynamical regimes of single- and two-phase magmatic flow in dikes, due to temperature-dependent viscosity and dissipation, was investigated using finite element calculations of magma flow in dikelike channels with length-to-width ratios of 1000:1 or more. Solutions of the steady state equations governing magma flow are obtained for a variety of conditions ranging from idealized plane-parallel models to cases involving nonparallel geometry and two-phase flows. The implications of the numerical simulations for the dynamics of flow in a dike-reservoir system and the consequences of dike entrance conditions on magmatic storage are discussed. Consideration is also given to an unmixing/self-lubrication mechanism which may be important for the lubrication of silicic magmas rising to the earth's surface in mixed magma ascent scenarios, which naturally segregates magma mixtures of two components with differing viscosities to minimize the driving pressure gradient.

Carrigan, Charles R.↗

Engineering Model Propellant Feed System Development for an Iodine Hall Thruster Demonstration Mission

CUBESATS are relatively new spacecraft platforms that are typically deployed from a launch vehicle as a secondary payload, providing low-cost access to space for a wide range of end-users. These satellites are comprised of building blocks having dimensions of 10x10x10 cu cm and a mass of 1.33 kg (a 1-U size). While providing low-cost access to space, a major operational limitation is the lack of a propulsion system that can fit within a CubeSat and is capable of executing high (Delta)v maneuvers. This makes it difficult to use CubeSats on missions requiring certain types of maneuvers (i.e. formation flying, spacecraft rendezvous). Recently, work has been performed investigating the use of iodine as a propellant for Hall-effect thrusters (HETs) 2 that could subsequently be used to provide a high specific impulse path to CubeSat propulsion. 3, 4 Iodine stores as a dense solid at very low pressures, making it acceptable as a propellant on a secondary payload. It has exceptionally high ρIsp (density times specific impulse), making it an enabling technology for small satellite near-term applications and providing the potential for systems-level advantages over mid-term high power electric propulsion options. Iodine flow can also be thermally regulated, subliming at relatively low temperature (< 100 C) to yield I2 vapor at or below 50 torr. At low power, the measured performance of an iodine-fed HET is very similar to that of a state-of-the-art xenon-fed thruster. Just as importantly, the current-voltage discharge characteristics of low power iodine-fed and xenon-fed thrusters are remarkably similar, potentially reducing development and qualifications costs by making it possible to use an already-qualified xenon-HET PPU in an iodine-fed system. Finally, a cold surface can be installed in a vacuum test chamber on which expended iodine propellant can deposit. In addition, the temperature doesn't have to be extremely cold to maintain a low vapor pressure in the vacuum chamber (it is under 10(exp −6) torr at −75 C), making it possible to 'cryopump' the propellant with lower-cost recirculating refrigerant-based systems as opposed to using liquid nitrogen or low temperature gaseous helium cryopanels. In the present paper, we describe the design and testing of the engineering model propellant feed system for iSAT (see Fig. 1). The feed system is based around an iodine propellant reservoir and two proportional control valves (PFCVs) that meter the iodine flow to the cathode and anode. The flow is split upstream of the PFCVs to both components can be fed from a common reservoir. Testing of the reservoir is reported to demonstrate that the design is capable of delivering the required propellant flow rates to operate the thruster. The tubing and reservoir are fabricated from hastelloy to resist corrosion by the heated gaseous iodine propellant. The reservoir, tubing, and PFCVs are heated to ensure the sublimed propellant will not re-deposit within the feed system. Heating is accomplished through a number of individual zones to control the overall power expended on heating the system and insulation is employed to minimize the amount of power used to heat the system prior to thruster operation.

Polzin, Kurt A.↗

Mantle devolatilization and rheology in the framework of planetary evolution

We investigate the thermal history of an Earth-like planet with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow. To calculate the spreading rate, we assume that the heat flow under the mid-ocean ridges is double the average mantle heat flow. The rate of regassing also depends on the seafloor spreading rate as well as on other factors like the efficiency of volatile recycling through island arc volcanism. Both mechanisms (de- and regassing) are coupled self-consistently with the help of the parameterized convection model under implementation of a temperature and volatile-content-dependent mantle viscosity.

Franck, S.↗

Nature of convection-stabilized dc arcs in dual-flow nozzle geometry. I - The cold flow field and dc arc characteristics. II - Optical diagnostics and theory

Steady-state low-current air arcs in a dual-flow nozzle system are studied experimentally. The cold flow field with no arc is investigated using a 12.7-mm diameter dual-flow nozzle in a steady-flow facility. Mach number and mass flux distributions are determined for various nozzle-pressure ratios and nozzle-gap spacing. It is found that the shock waves in the converging-diverging nozzles result in a decrease in overal resistance by about 15 percent. Also, Schlieren and differential interferometry techniques are used to visualize the density gradients within the arc plasma and thermal mantle. Both optical techniques reveal a laminar arc structure for a reservoir pressure of 1 atm at various current levels. Experimentally determined axial static pressure and cold-flow mass flux rate distributions and a channel-flow model with constant arc temperatre are used to solve the energy integral for the arc radius as a function of axial distance. The arc electric field strength, voltage, resistance, and power are determined with Ohm's law and the total heat transfer is related to arc power.

Serbetci, Ilter↗

Temperature Oscillation in a Loop Heat Pipe with Gravity Assist

ATLAS Laser Thermal Control System (LTCS) thermal vacuum testing where the condenser-radiator was placed in a vertical position, it was found that the loop heat pipe (LHP) reservoir required much more control heater power than the analytical model had predicted. The required control heater power was also higher than the liquid subcooling entering the reservoir using the measured temperatures and the calculated mass flow rate based on steady state LHP operation. This presentation describes the investigation of the LHP behaviors under a gravity assist mode with a very cold radiator sink temperature and a large thermal mass attached to the evaporator. It is concluded that gravity caused the cold liquid to drop from the condenser-radiator to the reservoir, resulting in a rapid decrease of the reservoir temperature. When the reservoir temperature was increasing, a reverse flow occurred in the liquid line, carrying warm liquid to the condenser-radiator. Both events consumed the reservoir control heater power. The fall and rise of the reservoir temperature also caused the net heat input to the evaporator to vary due to the release and storage of the sensible heat of the thermal mass. The combination of these effects led to a persistent reservoir temperature oscillation and a repeated influx of cold liquid from the condenser. This was the root cause of the extraordinary high control heater power requirement in the LTCS TV test. Without gravity assist, such a persistent temperature oscillation will not be present.

000↗

The NATA code: Theory and analysis, volume 1

A computer program for calculating quasi-one-dimensional gas flow in axisymmetric and two-dimensional nozzles and rectangular channels is presented. Flow is assumed to start from a state of thermochemical equilibrium at a high temperature in an upstream reservoir. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. Electronic nonequilibrium effects can be included using a two-temperature model. An approximate laminar boundary layer calculation is given for the shear and heat flux on the nozzle wall. Boundary layer displacement effects on the inviscid flow are considered also. Chemical equilibrium and transport property calculations are provided by subroutines. The code contains precoded thermochemical, chemical kinetic, and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It provides calculations of the stagnation conditions on axisymmetric or two-dimensional models, and of the conditions on the flat surface of a blunt wedge. The primary purpose of the code is to describe the flow conditions and test conditions in electric arc heated wind tunnels.

Bade, W. L.↗

Calibration of Axisymmetric and Quasi-1D Solvers for High Enthalpy Nozzles

The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.

Papadopoulos, P. E.↗