Search NASA⌕ Search

SEARCH · Search NASA

Results for “1-D 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 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 NTRS1992

Persistence and predictability in a perfect model

A realistic two-level GCM is used to examine the relationship between predictability and persistence. Predictability is measured by the average divergence of ensembles of solutions starting from perturbed initial conditions, and persistence is defined in terms of the autocorrelation function based on a single long-term model integration. The average skill of the dynamical forecasts is compared with the skill of simple persistence-based statistical forecasts. For initial errors comparable in magnitude to present-day analysis errors, the statistical forecast loses all skill after about one week, reflecting the lifetime of the lowest frequency fluctuations in the model. Large ensemble mean dynamical forecasts would be expected to remain skillful for about 3 wk. The disparity between the skill of the statistical and dynamical forecasts is greater for the higher frequency modes, which have little memory beyond 1 d, yet remain predictable for about 2 wk. The results are analyzed in terms of two characteristic time scales.

Schubert, Siegfried D.↗
NASA NTRS1989

Neural networks for data compression and invariant image recognition

An approach to invariant image recognition (I2R), based upon a model of biological vision in the mammalian visual system (MVS), is described. The complete I2R model incorporates several biologically inspired features: exponential mapping of retinal images, Gabor spatial filtering, and a neural network associative memory. In the I2R model, exponentially mapped retinal images are filtered by a hierarchical set of Gabor spatial filters (GSF) which provide compression of the information contained within a pixel-based image. A neural network associative memory (AM) is used to process the GSF coded images. We describe a 1-D shape function method for coding of scale and rotationally invariant shape information. This method reduces image shape information to a periodic waveform suitable for coding as an input vector to a neural network AM. The shape function method is suitable for near term applications on conventional computing architectures equipped with VLSI FFT chips to provide a rapid image search capability.

Gardner, Sheldon↗
NASA NTRS1988

A numerical simulation of the full two-dimensional electrothermal de-icer pad

The ability to predict the time-temperature history of electrothermal de-icer pads is important in the subsequent design of improved and more efficient versions. These de-icer pads are installed near the surface of aircraft components, for the specific purpose of removing accreted ice. The proposed numerical model can incorporate the full 2-D geometry through a section of a region (i.e., section of an airfoil), that current 1-D numerical codes are unable to do. Thus, the effects of irregular layers, curvature, etc., can now be accounted for in the thermal transients. Each layer in the actual geometry is mapped via a body-fitted coordinate transformation into uniform, rectangular computational grids. The relevant heat transfer equations are transformed and discretized. To model the phase change that might occur in any accreted ice, in an enthalpy formulation the phase change equations are likewise transformed and discretized. The code developed was tested against numerous classical numerical solutions, as well as against experimental de-icing data on a UH1H rotor blade obtained from the NASA Lewis Research Center. The excellent comparisons obtained show that this code can be a useful tool in predicting the performance of current de-icer models, as well as in the designing of future models.

Masiulaniec, Konstanty C.↗
NASA NTRS1991

Effect of wave-particle interaction on the outflow of ions at high latitudes

The objective of the research is to use the low-frequency wave spectrum measured by the Plasma Wave Instrument on the DE-1 spacecraft to include the wave-particle interaction (WPI) in the different polar wind models that are available at Utah State University. A Monte Carlo technique was used to simulate the ion diffusion in the velocity space due to scattering by the waves. This enabled us to study the effect of WPI on the magnitude as well as the composition of the outflow of the ionospheric ions. In particular, in the first year the 1-D semi-kinetic code developed by Barakat and Schunk (1983) was modified in order to include the effect of the WPI.

Barakat, Abdallah R.↗
NASA NTRS1993

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem (`Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗
NASA NTRS1992

The 3D rocket combustor acoustics model

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Analyses including radial and hub baffles and absorbers can be performed in one, two, and three dimensions. Pressure and velocity oscillations calculated using this procedure are presented for the SSME to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Comparisons are made between the results obtained utilizing 1-D, 2-D, and 3-D assumptions with regards to capturing the physical phenomena of interest and computational requirements.

Priem, Richard J.↗
NASA NTRS1989

A first scramjet study

A variety of related scramjet engine topics are examined. The flow is assumed to be 1-D, the gas is thermally and calorically perfect, and focus is on low hypersonic Mach numbers. The thrust and lift of an exposed half nozzle, which is used on the aerospace plane, is evaluated as well as a fully confined nozzle. A rough estimate of the drag of an aerospace plane is provided. Thermal effects and shock waves are next discussed. A parametric scramjet model is then presented based on the influence coefficient method, which evaluates the dominant scramjet processes. The independent parameters are the ratio of specific heats, a nondimensional heat addition parameter, and four Mach numbers. The total thrust generated by the combustor and nozzle is shown to be independent of the heat release distribution and the combustor exit Mach number, providing thermal choking is avoided. An operating condition for the combustor is found that maximizes the thrust. An alternative condition is explored when this optimum is no longer realistic. This condition provides a favorable pressure gradient and a reasonable area ratio for the combustor. Parametric results based on the model is provided.

Emanuel, George↗
NASA NTRS1991

The Navy/NASA Engine Program (NNEP89): A user's manual

An engine simulation computer code called NNEP89 was written to perform 1-D steady state thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a set of standard components are connected at execution time to simulate almost any turbine engine configuration that the user could imagine. The code was used to simulate a wide range of engine cycles from turboshafts and turboprops to air turborockets and supersonic cruise variable cycle engines. Off design performance is calculated through the use of component performance maps. A chemical equilibrium model is incorporated to adequately predict chemical dissociation as well as model virtually any fuel. NNEP89 is written in standard FORTRAN77 with clear structured programming and extensive internal documentation. The standard FORTRAN77 programming allows it to be installed onto most mainframe computers and workstations without modification. The NNEP89 code was derived from the Navy/NASA Engine program (NNEP). NNEP89 provides many improvements and enhancements to the original NNEP code and incorporates features which make it easier to use for the novice user. This is a comprehensive user's guide for the NNEP89 code.

Plencner, Robert M.↗
NASA NTRS1992

High temperature furnace modeling and performance verifications

Analytical, numerical, and experimental studies were performed on two classes of high temperature materials processing sources for their potential use as directional solidification furnaces. The research concentrated on a commercially available high temperature furnace using a zirconia ceramic tube as the heating element and an Arc Furnace based on a tube welder. The first objective was to assemble the zirconia furnace and construct parts needed to successfully perform experiments. The 2nd objective was to evaluate the zirconia furnace performance as a directional solidification furnace element. The 3rd objective was to establish a data base on materials used in the furnace construction, with particular emphasis on emissivities, transmissivities, and absorptivities as functions of wavelength and temperature. A 1-D and 2-D spectral radiation heat transfer model was developed for comparison with standard modeling techniques, and were used to predict wall and crucible temperatures. The 4th objective addressed the development of a SINDA model for the Arc Furnace and was used to design sample holders and to estimate cooling media temperatures for the steady state operation of the furnace. And, the 5th objective addressed the initial performance evaluation of the Arc Furnace and associated equipment for directional solidification. Results of these objectives are presented.

Smith, James E., Jr.↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1993

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗
NASA NTRS1992

Two-photon annihilation of thermal pairs in strong magnetic fields

The annihilation spectrum of pairs with 1-D thermal distributions in the presence of a strong magnetic field is calculated. Numerical analysis of the spectrum are performed for mildly relativistic temperatures and for different angles of emission with respect to field lines. Teragauss magnetic fields are assumed so that conditions are typical of gamma ray burst and pulsar environments. The spectra at each viewing angle reveal asymmetric line profiles that are signatures of the magnetic broadening and red shifting of the line: these asymmetries are more prominent for small viewing angles. Thermal Doppler broadening tends to dominate in the right wing of the line and obscures the magnetic broadening more at high temperatures and smaller viewing angles. This angular dependence of the line asymmetry may prove a valuable diagnostic tool. For low temperatures and magnetic field strengths, useful analytic expressions are presented for the line width, and also for the annihilation spectrum at zero viewing angle. The results presented find application in gamma ray burst and pulsar models, and may prove very helpful in deducing field strengths and temperatures of the emission regions of these objects from line observations made by Compton GRO and future missions.

Baring, Matthew G.↗
NASA NTRS1989

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗
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 NTRS1992

Singularity classification as a design tool for multiblock grids

A major stumbling block in interactive design of 3-D multiblock grids is the difficulty of visualizing the design as a whole. One way to make this visualization task easier is to focus, at least in early design stages, on an aspect of the grid which is inherently easy to present graphically, and to conceptualize mentally, namely the nature and location of singularities in the grid. The topological behavior of a multiblock grid design is determined by what happens at its edges and vertices. Only a few of these are in any way exceptional. The exceptional behaviors lie along a singularity graph, which is a 1-D construct embedded in 3-D space. The varieties of singular behavior are limited enough to make useful symbology on a graphics device possible. Furthermore, some forms of block design manipulation that appear appropriate to the early conceptual-modeling phase can be accomplished on this level of abstraction. An overview of a proposed singularity classification scheme and selected examples of corresponding manipulation techniques is presented.

Jones, Alan K.↗
NASA NTRS1992

Density measurements in an expanding flow using holographic interferometry

A nonequilibrium expansion of nitrogen through a 2-dimensional nozzle of a reflected shock tunnel is investigated with laser holography, time resolved pressure measurements and emission spectroscopy to characterize the flow field for vibrational relaxation measurements. A 2D computation is used to simulate the nozzle flow field. Synthetic holograms are generated using the computed density profiles and are compared with the experimental holograms. Also, using a quasi 1-D approximation, populations for each vibrational level are computed by (1) solving the master equations and (2) using the Landau-Teller formulation. Computational accuracy of both models is assessed against the experimental data.

Sharma, Surendra P.↗