SEARCH · Search NASA

Results for “quasi-one-dimensional”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1991

The coupling of radiative transfer to quasi 1-D flows with thermochemical nonequilibrium

Quasi-one-dimensional nonequilibrium nozzle flows with coupled radiative transfer are considered. The strongly coupled formulation of radiation and flowfield leads to a governing set of integro-differential equations. A fully implicit numerical method using the full matrix inversion or block iteration methods is presented to solve these equations. The nonequilibrium gas model consists of two chemical species, molecular and atomic nitrogen. The thermodynamic state of the gas is described by two temperatures, translational-rotational and vibrational, and the thermal radiation is assumed to be governed by the vibrational temperature. In radiative transfer, gases are assumed to be absorbing and emitting, and a detailed spectral dependency of the absorption coefficient is prescribed for a non-gray gas. The numerical solutions of strongly radiating nonequilibrium flows are presented for both gray and non-gray gases.

Gokcen, Tahir↗

NASA NTRS1991

Nonlinear analysis of longitudinal-mode liquid propellant rocket combustion instability

A computational technique for the nonlinear analysis of longitudinal-mode liquid propellant rocket combustion instability is examined based on the unsteady, quasi-one-dimensional Euler equations with appropriate source terms introduced to account for interphase transport coupling with the spray. The method is first assessed for unsteady, nonreacting, isentropic duct flow with specified admittances at the outflow and inflow boundaries. For small amplitude disturbances, numerical results based on a sufficient number of grid points compare favorably with predictions from a small-disturbance linearized analysis.

Jeng, San-Mou↗

NASA NTRS1990

Real-time simulation of supersonic inlets

A previously published real-time simulation algorithm, the matrix stability region placement (MSRP) method, is used to simulate a small perturbation model of the NASA Lewis Mach 2.5 40-60 mixed compression inlet. The model is representative of high-speed internal flow propulsion systems which can be approximated as quasi-one-dimensional flows. The resulting system of equations, which is stiff, is also simulated by the second-order Adam-Bashforth method. It is shown that the MSRP method can be used to simulate small perturbation models of high-speed internal flow propulsion systems in real time.

Mossayebi, F.↗

NASA NTRS1991

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant.

Hartley, Tom T.↗

NASA NTRS1991

Physical lumping methods for developing linear reduced models for high speed propulsion systems

In gasdynamic systems, information travels in one direction for supersonic flow and in both directions for subsonic flow. A shock occurs at the transition from supersonic to subsonic flow. Thus, to simulate these systems, any simulation method implemented for the quasi-one-dimensional Euler equations must have the ability to capture the shock. In this paper, a technique combining both backward and central differencing is presented. The equations are subsequently linearized about an operating point and formulated into a linear state space model. After proper implementation of the boundary conditions, the model order is reduced from 123 to less than 10 using the Schur method of balancing. Simulations comparing frequency and step response of the reduced order model and the original system models are presented.

Immel, S. M.↗

NASA NTRS1991

Modeling of linear isentropic flow systems

A modeling approach for linear isentropic flow systems based on the quasi-one-dimensional Euler equations of non-viscous, compressible flow are presented. Such systems are representative of certain high speed propulsion systems. Accurate models useful in control system studies are developed. A supersonic inlet is considered, and the resulting set of partial differential equations with boundary conditions is solved for a linear transfer matrix using Laplace transforms.

Sarantopoulos, Athan D.↗

NASA NTRS1991

Transient multidimensional analysis of nonconventional heat pipes with uniform and nonuniform heat distributions

A numerical analysis of transient heat pipe performance including nonconventional heat pipes with nonuniform heat distributions is presented. A body-fitted grid system was applied to a three-dimensional wall and wick model, which was coupled with a transient compressible quasi-one-dimensional vapor flow model. The numerical results were first compared with experimental data from cylindrical heat pipes with good agreement. Numerical calculations were then made for a leading edge heat pipe with localized high heat fluxes. Performance characteristics different from conventional heat pipes are illustrated and some operating limits concerning heat pipe design are discussed.

Cao, Y.↗

NASA NTRS1991

Mach 4 and Mach 8 axisymmetric nozzles for a high-enthalpy shock tunnel

The performance of two axisymmetric nozzles which were designed to produce uniform, parallel flow with nominal Mach numbers of 4 and 8 is examined. A free-piston-driven shock tube was used to supply the nozzle with high-temperature, high-pressure test gas. The inviscid design procedure treated the nozzle expansion in two stages. Close to the nozzle throat, the nozzle wall was specified as conical and the gas flow was treated as a quasi-one-dimensional chemically-reacting flow. At the end of the conical expansion, the gas was assumed to be calorically perfect, and a contoured wall was designed (using method of characteristics) to convert the source flow into a uniform and parallel flow at the end of the nozzle. Performance was assessed by measuring Pitot pressures across the exit plane of the nozzles and, over the range of operating conditions examined, the nozzles produced satisfactory test flows. However, there were flow disturbances in the Mach 8 nozzle flow that persisted for significant times after flow initiation.

Jacobs, P. A.↗

NASA NTRS1992

Vibrational relaxation of anharmonic oscillators in expanding flows

Although the Landau-Teller vibrational model accurately predicts vibrational excitation process in post-shock and compressing flows, it under-predicts the rate of de-excitation in cooling and expanding flows. In the present paper, detailed calculations of the vibrational relaxation process of N2 and CO in cooling flows are conducted. A coupled set of vibrational transition rate equations and quasi-one-dimensional fluid dynamic equations is solved. Multiple quantum level transition rates are computed using SSH theory. The SSH transition rate results are compared with available experimental data and other theoretical models. Vibration-vibration exchange collisions are responsible for some vibrational relaxation acceleration in situations of high vibrational temperature and low translational temperature. The present results support the relaxation mechanisms proposed by Bray and by Treanor Rich and Rehm.

Ruffin, Stephen M.↗

NASA NTRS1993

Adjoint methods for aerodynamic wing design

A model inverse design problem is used to investigate the effect of flow discontinuities on the optimization process. The optimization involves finding the cross-sectional area distribution of a duct that produces velocities that closely match a targeted velocity distribution. Quasi-one-dimensional flow theory is used, and the target is chosen to have a shock wave in its distribution. The objective function which quantifies the difference between the targeted and calculated velocity distributions may become non-smooth due to the interaction between the shock and the discretization of the flowfield. This paper offers two techniques to resolve the resulting problems for the optimization algorithms. The first, shock-fitting, involves careful integration of the objective function through the shock wave. The second, coordinate straining with shock penalty, uses a coordinate transformation to align the calculated shock with the target and then adds a penalty proportional to the square of the distance between the shocks. The techniques are tested using several popular sensitivity and optimization methods, including finite-differences, and direct and adjoint discrete sensitivity methods. Two optimization strategies, Gauss-Newton and sequential quadratic programming (SQP), are used to drive the objective function to a minimum.

Grossman, Bernard↗

NASA NTRS1993

Source term evaluation for combustion modeling

A modification is developed for application to the source terms used in combustion modeling. The modification accounts for the error of the finite difference scheme in regions where chain-branching chemical reactions produce exponential growth of species densities. The modification is first applied to a one-dimensional scalar model problem. It is then generalized to multiple chemical species, and used in quasi-one-dimensional computations of shock-induced combustion in a channel. Grid refinement studies demonstrate the improved accuracy of the method using this modification. The algorithm is applied in two spatial dimensions and used in simulations of steady and unsteady shock-induced combustion. Comparisons with ballistic range experiments give confidence in the numerical technique and the 9-species hydrogen-air chemistry model.

Sussman, Myles A.↗

NASA NTRS1993

Time-dependent simulation of reflected-shock/boundary layer interaction

An initial numerical/experimental investigation has been done to better understand multi-dimensional flow phenomena inside pulse facilities. Time-dependent quasi-one-dimensional and axisymmetric numerical simulations of complete shock tube flow are compared with experimental pressure traces recorded at the NASA Ames electric-arc driven shock tube facility (from cold driver shots). Of particular interest is the interaction between the reflected shock wave and the boundary layer. Evidence of the shock bifurcation caused by this interaction is clearly seen in the present experimental data. The axisymmetric simulations reproduce this phenomenon and demonstrate how this interaction can provide a mechanism for driver gas to contaminate the stagnation region. The simulations incorporate finite-rate chemistry, a moving mesh and laminar viscosity.

Wilson, Gregory J.↗

NASA NTRS1993

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗

NASA NTRS1991

Modeling of linear isentropic flow systems

Sarantopoulos, Athan D.↗

NASA NTRS1993

Analysis and design of an ultrahigh temperature hydrogen-fueled MHD generator

A coupled gas dynamics/radiative heat transfer analysis of partially ionized hydrogen, in local thermodynamic equilibrium, flowing through an ultrahigh temperature (10,000-20,000 K) magnetohydrodynamic (MHD) generator is performed. Gas dynamics are modeled by a set of quasi-one-dimensional, nonlinear differential equations which account for friction, convective and radiative heat transfer, and the interaction between the ionized gas and applied magnetic field. Radiative heat transfer is modeled using nongray, absorbing-emitting 2D and 3D P-1 approximations which permit an arbitrary variation of the spectral absorption coefficient with frequency. Gas dynamics and radiative heat transfer are coupled through the energy equation and through the temperature- and density-dependent absorption coefficient. The resulting nonlinear elliptic problem is solved by iterative methods. Design of such MHD generators as onboard, open-cycle, electric power supplies for a particular advanced airbreathing propulsion concept produced an efficient and compact 128-MWe generator characterized by an extraction ratio of 35.5 percent, a power density of 10,500 MWe/cu m, and a specific (extracted) energy of 324 MJe/kg of hydrogen. The maximum wall heat flux and total wall heat load were 453 MW/sq m and 62 MW, respectively.

Moder, Jeffrey P.↗

NASA NTRS1993

Shape optimization governed by the Euler equations using an adjoint method

A numerical approach for the treatment of optimal shape problems governed by the Euler equations is discussed. Focus is on flows with embedded shocks. A very simple problem is considered: the design of a quasi-one-dimensional Laval nozzle. A cost function and a set of Lagrange multipliers are introduced to achieve the minimum. The nature of the resulting costate equations is discussed. A theoretical difficulty that arises for cases with embedded shocks is pointed out and solved. Finally, some results are given to illustrate the effectiveness of the method.

Iollo, Angelo↗

NASA NTRS1994

An analysis code for the Rapid Engineering Estimation of Momentum and Energy Losses (REMEL)

Nonideal behavior has traditionally been modeled by defining efficiency (a comparison between actual and isentropic processes), and subsequent specification by empirical or heuristic methods. With the increasing complexity of aeropropulsion system designs, the reliability of these more traditional methods is uncertain. Computational fluid dynamics (CFD) and experimental methods can provide this information but are expensive in terms of human resources, cost, and time. This report discusses an alternative to empirical and CFD methods by applying classical analytical techniques and a simplified flow model to provide rapid engineering estimates of these losses based on steady, quasi-one-dimensional governing equations including viscous and heat transfer terms (estimated by Reynold's analogy). A preliminary verification of REMEL has been compared with full Navier-Stokes (FNS) and CFD boundary layer computations for several high-speed inlet and forebody designs. Current methods compare quite well with more complex method results and solutions compare very well with simple degenerate and asymptotic results such as Fanno flow, isentropic variable area flow, and a newly developed, combined variable area duct with friction flow solution. These solution comparisons may offer an alternative to transitional and CFD-intense methods for the rapid estimation of viscous and heat transfer losses in aeropropulsion systems.

Dechant, Lawrence J.↗

NASA NTRS1994

A numerical analysis of shipboard and coastal zone color scanner time series of new production within Gulf Stream cyclonic eddies in the South Atlantic Bight

Eddy-induced upwelling occurs along the western edge of the Gulf Stream between Cape Canaveral, Florida, and Cape Hatteras, North Carolina, in the South Atlantic Bight (SAB). Coastal zone color scanner images of 1-km resolution spanning the period April 13-21, 1979, were processed to examine these eddy features in relation to concurrent shipboard and current/temperature measurements at moored arrays. A quasi-one-dimensional (z), time dependent biological model, using only nitrate as a nutrient source, has been combined with a three-dimensional physical model in an attempt to replicate the observed phytoplankton field at the northward edge of an eddy. The model is applicable only to the SAB south of the Charleston Bump, at approximately 31.5 deg N, since no feature analogous to the bump exists in the model bathymetry. The modeled chlorophyll, nitrate, and primary production fields of the euphotic zone are very similar to those obtained from the satellite and shipboard data at the leading edges of the observed eddies south of the Charleston Bump. The horizontal and vertical simulated fluxes of nitrate and chlorophyll show that only approximately 10% of the upwelled nitrate is utilized by the phytoplankton of the modeled grid box on the northern edge of the cyclone, while approximately 75% is lost horizontally, with the remainder still in the euphotic zone after the 10-day period of the model. Loss of chlorophyll due to sinking is very small in this strong upwelling region of the cyclone. The model is relatively insensitive to variations in the sinking parameterization and the external nitrate and chlorophyll fields but is very sensitive to a reduction of the maximum potential growth rate to half that measured. Given the success of this model in simulating the new production of the selcted upwelling region, other upwelling regions for which measurements or successful models of physical and biological quantities and rates exist could be modeled similarly.

Pribble, J. Raymond↗