Search NASA⌕ Search

SEARCH · Search NASA

Results for “1-D fluid flow”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

A comparison of numerical solutions to the inviscid equations of fluid motion

The second-order MacCormack method and the third-order Rusanov and Kutler-Warming-Lomax methods are applied to the inviscid 1-d Burgers' equation, wedge flow, and the problem of shock reflection from a rigid boundary. The numerical solution in each case is compared to the exact solution and the quantitative estimates of accuracy are obtained. Results of this study show that the third-order Kutler-Warming-Lomax or the Rusanov methods tuned for minimum dissipation or minimum dispersion provide the most accurate solution in each of the examples considered.

Anderson, D. A.↗
NASA NTRS1990

Prediction of the ullage gas thermal stratification in a NASP vehicle propellant tank experimental simulation using FLOW-3D

As part of the National Aero-Space Plane (NASP) project, the multi-dimensional effects of gravitational force, initial tank pressure, initial ullage temperature, and heat transfer rate on the 2-D temperature profiles were studied. FLOW-3D, a commercial finite difference fluid flow model, was used for the evaluation. These effects were examined on the basis of previous liquid hydrogen experimental data with gaseous hydrogen pressurant. FLOW-3D results were compared against an existing 1-D model. In addition, the effects of mesh size and convergence criteria on the analytical results were investigated. Suggestions for future modifications and uses of FLOW-3D for modeling of a NASP tank are also presented.

Hardy, Terry L.↗
NASA NTRS1998

1-D Numerical Analysis of RBCC Engine Performance

An RBCC engine combines air breathing and rocket engines into a single engine to increase the specific impulse over an entire flight trajectory. Considerable research pertaining to RBCC propulsion was performed during the 1960's and these engines were revisited recently as a candidate propulsion system for either a single-stage-to-orbit (SSTO) or two-stage-to-orbit (TSTO) launch vehicle. There are a variety of RBCC configurations that had been evaluated and new designs are currently under development. However, the basic configuration of all RBCC systems is built around the ejector scramjet engine originally developed for the hypersonic airplane. In this configuration, a rocket engine plays as an ejector in the air-augmented initial acceleration mode, as a fuel injector in scramjet mode and the rocket in all rocket mode for orbital insertion. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in RBCC propulsion systems. The objective of the present research was to develop a transient 1-D numerical model that could be used to predict flow behavior throughout a generic RBCC engine following a flight path.

Han, Samuel S.↗
NASA NTRS1999

1-D Numerical Analysis of ABCC Engine Performance

ABCC engine combines air breathing and rocket engine into a single engine to increase the specific impulse over an entire flight trajectory. Except for the heat source, the basic operation of the ABCC is similar to the basic operation of the RBCC engine. The ABCC is intended to have a higher specific impulse than the RBCC for single stage Earth to orbit vehicle. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in ABCC propulsion system. The objective of the present research was to develop a transient 1-D numerical model using conservation of mass, linear momentum, and energy equations that could be used to predict flow behavior throughout a generic ABCC engine following a flight path. At specific points during the development of the 1-D numerical model a myriad of tests were performed to prove the program produced consistent, realistic numbers that follow compressible flow theory for various inlet conditions.

Holden, Richard↗
NASA NTRS1998

Physics of Boundaries and Their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of January. During this time, we primarily worked on further analysis of the results presented at the AGU as well as writing them up for publication. Using large scale simulations, we showed that the magnetopause during the southward IMF case is quite irregular with varying thickness, and has a complex flow pattern owing to the nonlinear effects of the convective flow superimposed on the flows generated in the reconnection layer. We used inflow-outflow boundary conditions to examine the kinetic nature of the discontinuities that are formed in the reconnection layer and concluded that nonlocal effects play a major role in the formation of such discontinuities and can alter their properties from the usual structures expected from 1-D simulations or from fluid theories. Finally, we used our 3-D simulations to examine the nonlinear interaction of the tearing mode with the Kelvin-Helmholtz instability. We showed that this interaction leads to the generation of a large core field which is observed both in the magnetotail as well as the magnetopause.

Omidi, Nojan↗
NASA NTRS1993

SAFSIM overview

An overview of the systems analysis flow simulator (SAFSIM) computer program is provided. SAFSIM is being developed at Sandia National Laboratories and is currently funded by the Air Force Space Nuclear Thermal Propulsion (SNTP) Program. SAFSIM is a general purpose, Fortran computer program to simulate the integrated performance of complex systems involving fluid mechanics, heat transfer, and reactor dynamics. SAFSIM provides sufficient versatility to allow the engineering simulation of almost any system. SAFSIM is based on a 1-D finite element model and provides the analyst with approximate solutions to complex problems.

Dobranich, Dean↗
NASA NTRS1999

Inlet-Compressor Analysis using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. Isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (1-D) Euler inlet simulation showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition. To aid in this understanding and to help evaluate possible boundary conditions, an inlet-engine experiment was conducted at the University of Cincinnati. The interaction of acoustic pulses, generated in the inlet, with the engine were investigated. Because of the availability of experimental data for validation, it was decided to simulate the experiment using CFD. The philosophy here is that the inlet-engine system is best simulated by coupling (existing) specialized CFD component-codes. The objectives of this work were to aid in a better understanding of inlet-compressor interaction physics and the formulation of a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets. Previous simulations have used 1-D Euler engine simulations in conjunction with 1-D Euler and axisymmetric Euler inlet simulations. This effort is a first step toward CFD simulation of an entire engine by coupling multidimensional component codes.

Cole, Gary↗
NASA NTRS2001

Thermal Fault Tolerance Analysis of Carbon Fiber Rope Barrier Systems for Use in the Reusable Solid Rocket Motor ( RSRM) Nozzle Joints

Carbon Fiber Rope (CFR) thermal barrier systems are being considered for use in several RSRM (Reusable Solid Rocket Motor) nozzle joints as a replacement for the current assembly gap close-out process/design. This study provides for development and test verification of analysis methods used for flow-thermal modeling of a CFR thermal barrier subject to fault conditions such as rope combustion gas blow-by and CFR splice failure. Global model development is based on a 1-D (one dimensional) transient volume filling approach where the flow conditions are calculated as a function of internal 'pipe' and porous media 'Darcy' flow correlations. Combustion gas flow rates are calculated for the CFR on a per-linear inch basis and solved simultaneously with a detailed thermal-gas dynamic model of a local region of gas blow by (or splice fault). Effects of gas compressibility, friction and heat transfer are accounted for the model. Computational Fluid Dynamic (CFD) solutions of the fault regions are used to characterize the local flow field, quantify the amount of free jet spreading and assist in the determination of impingement film coefficients on the nozzle housings. Gas to wall heat transfer is simulated by a large thermal finite element grid of the local structure. The employed numerical technique loosely couples the FE (Finite Element) solution with the gas dynamics solution of the faulted region. All free constants that appear in the governing equations are calibrated by hot fire sub-scale test. The calibrated model is used to make flight predictions using motor aft end environments and timelines. Model results indicate that CFR barrier systems provide a near 'vented joint' style of pressurization. Hypothetical fault conditions considered in this study (blow by, splice defect) are relatively benign in terms of overall heating to nozzle metal housing structures.

Clayton, J. Louie↗
NASA NTRS1992

An extended Lagrangian method for subsonic flows

It is well known that fluid motion can be specified by either the Eulerian of Lagrangian description. Most of Computational Fluid Dynamics (CFD) developments over the last three decades have been based on the Eulerian description and considerable progress has been made. In particular, the upwind methods, inspired and guided by the work of Gudonov, have met with many successes in dealing with complex flows, especially where discontinuities exist. However, this shock capturing property has proven to be accurate only when the discontinuity is aligned with one of the grid lines since most upwind methods are strictly formulated in 1-D framework and only formally extended to multi-dimensions. Consequently, the attractive property of crisp resolution of these discontinuities is lost and research on genuine multi-dimensional approach has just been undertaken by several leading researchers. Nevertheless they are still based on the Eulerian description.

Liou, Meng-Sing↗
NASA NTRS2000

Internal Flow Analysis of Large L/D Solid Rocket Motors

Traditionally, Solid Rocket Motor (SRM) internal ballistic performance has been analyzed and predicted with either zero-dimensional (volume filling) codes or one-dimensional ballistics codes. One dimensional simulation of SRM performance is only necessary for ignition modeling, or for motors that have large length to port diameter ratios which exhibit an axial "pressure drop" during the early burn times. This type of prediction works quite well for many types of motors, however, when motor aspect ratios get large, and port to throat ratios get closer to one, two dimensional effects can become significant. The initial propellant grain configuration for the Space Shuttle Reusable Solid Rocket Motor (RSRM) was analyzed with 2-D, steady, axi-symmetric computational fluid dynamics (CFD). The results of the CFD analysis show that the steady-state performance prediction at the initial burn geometry, in general, agrees well with 1-D transient prediction results at an early time, however, significant features of the 2-D flow are captured with the CFD results that would otherwise go unnoticed. Capturing these subtle differences gives a greater confidence to modeling accuracy, and additional insight with which to model secondary internal flow effects like erosive burning. Detailed analysis of the 2-D flowfield has led to the discovery of its hidden 1-D isentropic behavior, and provided the means for a thorough and simplified understanding of internal solid rocket motor flow. Performance parameters such as nozzle stagnation pressure, static pressure drop, characteristic velocity, thrust and specific impulse are discussed in detail and compared for different modeling and prediction methods. The predicted performance using both the 1-D codes and the CFD results are compared with measured data obtained from static tests of the RSRM. The differences and limitations of predictions using ID and 2-D flow fields are discussed and some suggestions for the design of large L/D motors and more critically, motors with port to throat ratios near one, are covered.

Laubacher, Brian A.↗
NASA NTRS1990

Carbon deposition model for oxygen-hydrocarbon combustion. Task 6: Data analysis and formulation of an empirical model

The formation and deposition of carbon (soot) was studied in the Carbon Deposition Model for Oxygen-Hydrocarbon Combustion Program. An empirical, 1-D model for predicting soot formation and deposition in LO2/hydrocarbon gas generators/preburners was derived. The experimental data required to anchor the model were identified and a test program to obtain the data was defined. In support of the model development, cold flow mixing experiments using a high injection density injector were performed. The purpose of this investigation was to advance the state-of-the-art in LO2/hydrocarbon gas generator design by developing a reliable engineering model of gas generator operation. The model was formulated to account for the influences of fluid dynamics, chemical kinetics, and gas generator hardware design on soot formation and deposition.

Makel, Darby B.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗
NASA NTRS1989

Thermodynamic and fluid mechanic analysis of rapid pressurization in a dead-end tube

Three models have been applied to very rapid compression of oxygen in a dead-ended tube. Pressures as high as 41 MPa (6000 psi) leading to peak temperatures of 1400 K are predicted. These temperatures are well in excess of the autoignition temperature (750 K) of teflon, a frequently used material for lining hoses employed in oxygen service. These findings are in accord with experiments that have resulted in ignition and combustion of the teflon, leading to the combustion of the stainless steel braiding and catastrophic failure. The system analyzed was representative of a capped off-high-pressure oxygen line, which could be part of a larger system. Pressurization of the larger system would lead to compression in the dead-end line, and possible ignition of the teflon liner. The model consists of a large plenum containing oxygen at the desired pressure (500 to 6000 psi). The plenum is connected via a fast acting valve to a stainless steel tube 2 cm inside diameter. Opening times are on the order of 15 ms. Downstream of the valve is an orifice sized to increase filling times to around 100 ms. The total length from the valve to the dead-end is 150 cm. The distance from the valve to the orifice is 95 cm. The models describe the fluid mechanics and thermodynamics of the flow, and do not include any combustion phenomena. A purely thermodynamic model assumes filling to be complete upstream of the orifice before any gas passes through the orifice. This simplification is reasonable based on experiment and computer modeling. Results show that peak temperatures as high as 4800 K can result from recompression of the gas after expanding through the orifice. An approximate transient model without an orifice was developed assuming an isentropic compression process. An analytical solution was obtained. Results indicated that fill times can be considerably shorter than valve opening times. The third model was a finite difference, 1-D transient compressible flow model. Results from the code show the recompression effect but predict much lower peak temperatures than the thermodynamic model.

Leslie, Ian H.↗
NASA NTRS1995

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.↗
NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗
NASA NTRS1989

A finite difference treatment of Stokes-type flows: Preliminary report

The equations Laplacian operator omega = 0, (1.1a) and omega = Laplacian operator Chi, (1.1b) describe, in suitable units, 2-D Stokes flow of an incompressible fluid occupying a domain D in which omega is the vorticity and Chi is the stream function. The flow is uniquely determined by specifying the velocity on the boundary B of D, a condition which leads to specifying the stream function Chi and its normal derivative Chi sub n on B. A mathematically similar problem arises in describing the equilibrium of a flat plate in structural mechanics where a related 1-D problem by finite difference or finite element methods is to introduce effective methods for imposing the boundary conditions through which (1.1a) is coupled to (1.1b). These models thus provide a simple starting point for examining the general treatment of boundary conditions for more general time dependent Navier-Stokes incompressible flows. For the purpose of discussion it is assumed that D is a square domain. A standard finite difference method to solve (1.1) is to introduce a uniform grid and then use standard five point finite difference operators to express each equation in (1.1). At any point on the boundary B a value of Chi is specified by the boundary conditions but a value of omega at the same boundary mesh point will also be required to complete the computation. Methods are discussed which overcome the difficulty in solving these problems.

Rose, M. E.↗
NASA NTRS1990

One- and two-dimensional Stirling machine simulation using experimentally generated flow turbulence models

Investigations of one- and two-dimensional (1- or 2-D) simulations of Stirling machines centered around experimental data generated by the U. of Minnesota Mechanical Engineering Test Rig (METR) are covered. This rig was used to investigate oscillating flows about a zero mean with emphasis on laminar/turbulent flow transitions in tubes. The Space Power Demonstrator Engine (SPDE) and in particular, its heater, were the subjects of the simulations. The heater was treated as a 1- or 2-D entity in an otherwise 1-D system. The 2-D flow effects impacted the transient flow predictions in the heater itself but did not have a major impact on overall system performance. Information propagation effects may be a significant issue in the simulation (if not the performance) of high-frequency, high-pressure Stirling machines. This was investigated further by comparing a simulation against an experimentally validated analytic solution for the fluid dynamics of a transmission line. The applicability of the pressure-linking algorithm for compressible flows may be limited by characteristic number (defined as flow path information traverses per cycle); this warrants further study. Lastly the METR was simulated in 1- and 2-D. A two-parameter k-w foldback function turbulence model was developed and tested against a limited set of METR experimental data.

Goldberg, Louis F.↗
NASA NTRS1993

Low Reynolds number kappa-epsilon and empirical transition models for oscillatory pipe flow and heat transfer

Stirling engine heat exchangers are shell-and-tube type with oscillatory flow (zero-mean velocity) for the inner fluid. This heat transfer process involves laminar-transition turbulent flow motions under oscillatory flow conditions. A low Reynolds number kappa-epsilon model, (Lam-Bremhorst form), was utilized in the present study to simulate fluid flow and heat transfer in a circular tube. An empirical transition model was used to activate the low Reynolds number k-e model at the appropriate time within the cycle for a given axial location within the tube. The computational results were compared with experimental flow and heat transfer data for: (1) velocity profiles, (2) kinetic energy of turbulence, (3) skin friction factor, (4) temperature profiles, and (5) wall heat flux. The experimental data were obtained for flow in a tube (38 mm diameter and 60 diameter long), with the maximum Reynolds number based on velocity being Re(sub max) = 11840, a dimensionless frequency (Valensi number) of Va = 80.2, at three axial locations X/D = 16, 30 and 44. The agreement between the computations and the experiment is excellent in the laminar portion of the cycle and good in the turbulent portion. Moreover, the location of transition was predicted accurately. The Low Reynolds Number kappa-epsilon model, together with an empirical transition model, is proposed herein to generate the wall heat flux values at different operating parameters than the experimental conditions. Those computational data can be used for testing the much simpler and less accurate one dimensional models utilized in 1-D Stirling Engine design codes.

Bauer, Christopher↗