Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Numerical modeling”

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 55 records · Page 3

Numerical Model of Flame Spread Over Solids in Microgravity: A Supplementary Tool for Designing a Space Experiment

The recently developed numerical model of concurrent-flow flame spread over thin solids has been used as a simulation tool to help the designs of a space experiment. The two-dimensional and three-dimensional, steady form of the compressible Navier-Stokes equations with chemical reactions are solved. With the coupled multi-dimensional solver of the radiative heat transfer, the model is capable of answering a number of questions regarding the experiment concept and the hardware designs. In this paper, the capabilities of the numerical model are demonstrated by providing the guidance for several experimental designing issues. The test matrix and operating conditions of the experiment are estimated through the modeling results. The three-dimensional calculations are made to simulate the flame-spreading experiment with realistic hardware configuration. The computed detailed flame structures provide the insight to the data collection. In addition, the heating load and the requirements of the product exhaust cleanup for the flow tunnel are estimated with the model. We anticipate that using this simulation tool will enable a more efficient and successful space experiment to be conducted.

Shih, Hsin-Yi

Numerical Modeling and Testing of an Inductively-Driven and High-Energy Pulsed Plasma Thrusters

Pulsed Plasma Thrusters (PPTs) are advanced electric space propulsion devices that are characterized by simplicity and robustness. They suffer, however, from low thrust efficiencies. This summer, two approaches to improve the thrust efficiency of PPTs will be investigated through both numerical modeling and experimental testing. The first approach, an inductively-driven PPT, uses a double-ignition circuit to fire two PPTs in succession. This effectively changes the PPTs configuration from an LRC circuit to an LR circuit. The LR circuit is expected to provide better impedance matching and improving the efficiency of the energy transfer to the plasma. An added benefit of the LR circuit is an exponential decay of the current, whereas a traditional PPT s under damped LRC circuit experiences the characteristic "ringing" of its current. The exponential decay may provide improved lifetime and sustained electromagnetic acceleration. The second approach, a high-energy PPT, is a traditional PPT with a variable size capacitor bank. This PPT will be simulated and tested at energy levels between 100 and 450 joules in order to investigate the relationship between efficiency and energy level. Arbitrary Coordinate Hydromagnetic (MACH2) code is used. The MACH2 code, designed by the Center for Plasma Theory and Computation at the Air Force Research Laboratory, has been used to gain insight into a variety of plasma problems, including electric plasma thrusters. The goals for this summer include numerical predictions of performance for both the inductively-driven PPT and high-energy PFT, experimental validation of the numerical models, and numerical optimization of the designs. These goals will be met through numerical and experimental investigation of the PPTs current waveforms, mass loss (or ablation), and impulse bit characteristics.

Parma, Brian

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen

Comparison of Numerical Model Results with Diffusion Flames in Microgravity

The effects of gravity on methane diffusion flames were studied using a finite difference numerical model. Variables under consideration included velocity, pressure, temperature, enthalpy, reactive species, and inert carrier gas. Results obtained showed that flames in zero gravity were wider and taller than those in normal gravity. The use of a six reaction scheme produced similar predictions of flame width and height compared to those of the single reaction case. The numerical models predicted a considerable difference in the flame base region between normal gravity and low gravity. This region is considered to be beyond the boundaries of current analytical models.

Hyer, Paul V.

Numerical modeling of two-dimensional confined flows

A numerical model of two-dimensional confined flows is presented. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying conservation principles directly to the individual streams. A listing of a computer code based on this approach in FORTRAN 4 language is presented. The code computes two dimensional compressible turbulent flows in ducts when the duct area along the flow is specified and the pressure gradient is unknown.

Greywall, M. S.

A numerical model for low-latitude ionospheric TEC

A numerical model of total electron content (TEC) at low latitudes has been constructed from the Faraday rotation data recorded at a number of locations in India during the period Oct. 1975-July 1976 when the geostationary satellite ATS-6 was located at 35 deg E. In all, 36 coefficients for each season are required to represent the model. The model can be used for satellite tracking systems in estimating quickly the range, range rate and angular refraction errors. The present model is an improvement over the previously constructed model based on orbiting satellite data.

Sethia, G.

An improved numerical model for wave rotor design and analysis

A numerical model has been developed which can predict both the unsteady flows within a wave rotor and the steady averaged flows in the ports. The model is based on the assumptions of one-dimensional, unsteady, and perfect gas flow. Besides the dominant wave behavior, it is also capable of predicting the effects of finite tube opening time, leakage from the tube ends, and viscosity. The relative simplicity of the model makes it useful for design, optimization, and analysis of wave rotor cycles for any application. This paper discusses some details of the model and presents comparisons between the model and two laboratory wave rotor experiments.

Paxson, Daniel E.

An improved numerical model for wave rotor design and analysis

A numerical model has been developed which can predict both the unsteady flows within a wave rotor and the steady averaged flows in the ports. The model is based on the assumptions of one-dimensional, unsteady, and perfect gas flow. Besides the dominant wave behavior, it is also capable of predicting the effects of finite tube opening time, leakage from the tube ends, and viscosity. The relative simplicity of the model makes it useful for design, optimization, and analysis of wave rotor cycles for any application. This paper discusses some details of the model and presents comparisons between the model and two laboratory wave rotor experiments.

Paxson, Daniel E.

Implementation of a Sage-Based Stirling Model Into a System-Level Numerical Model of the Fission Power System Technology Demonstration Unit

The Fission Power System (FPS) project is developing a Technology Demonstration Unit (TDU) to verify the performance and functionality of a subscale version of the FPS reference concept in a relevant environment, and to verify component and system models. As hardware is developed for the TDU, component and system models must be refined to include the details of specific component designs. This paper describes the development of a Sage-based pseudo-steady-state Stirling convertor model and its implementation into a system-level model of the TDU.

Briggs, Maxwell H.

Testing and numerical modeling of hypervelocity impact damaged Space Station multilayer insulation

Results are presented of experiments measuring the degradation of the insulating capabilities of the multilayer insulation (MLI) of the Space Station Freedom, when subjected to hypervelocity impact damage. A simple numerical model was developed for use in an engineering design environment for quick assessment of thermal effect of the impact. The model was validated using results from thermal vacuum tests on MLI with simulated damage. The numerical model results agreed with experimental data.

Rule, William K.

Numerical models for cometary dust environments

The development of a numerical model for the dust environment near a comet is described. The major components of the model are a simple treatment of the acceleration of the dust by the gas outflow from an isothermal icy nucleus, a detailed treatment of the dust particle trajectories controlled by solar gravity and radiation pressure, including the formation of the dust envelope and tail, and a simple treatment of the radiation from the dust in the visible and thermal infrared. The results of recent comet observations are used to estimate parameter values for the model, including dust-to-gas ratios by light and by mass, the distribution of dust and gas emission on the nuclear surface, the distributions of dust particle size, density, mass, and radiative properties, and other quantities. Gas and dust production rates as functions of time are derived from observed light curves using the simple photometric theory. Sample results are presented for the dust concentration, flux and fluence along the trajectory of a spacecraft which flies by Comet P/Halley in 1986.

Divine, N.

Numerical models of extragalactic radio sources

When supercomputer-implemented numerical simulations analyzing the nonlinear physics inherent in the hydrodynamic and MHD equations are applied to extragalactic radio sources, many of the complex structures observed on telescopic images are reproduced. Attention is presently given to recently obtained results from 2D and 3D numerical simulations of the formation and evolution of extended radio morphologies; these numerical models allow the exploration of such physical phenomena as the role of magnetic fields in the dynamics and emissivity of extended radio galaxies, intermittent outflow from the cores of active galaxies, fluid-jet instabilities, and the bending of collimated outflows by motion through the intergalactic medium.

Burns, Jack O.

Numerical modeling of preburner flowfield

This work is intended to numerically predict the flowfields inside the preburner of the Space Shuttle Main Engine. The computer code (FDNS) based on pressure correction method is modified and adapted with an elliptic grid generator. The original configuration of the preburner in conjunction with downstream gas turbines has been simplified geometrically and numerically modeled at its full power in this work. The computational results are presented and qualitatively discussed with test data collected in NASA/MSFC.

Chow, A. S.

Noninvasive assessment of mitral inertness: clinical results with numerical model validation

Inertial forces (Mdv/dt) are a significant component of transmitral flow, but cannot be measured with Doppler echo. We validated a method of estimating Mdv/dt. Ten patients had a dual sensor transmitral (TM) catheter placed during cardiac surgery. Doppler and 2D echo was performed while acquiring LA and LV pressures. Mdv/dt was determined from the Bernoulli equation using Doppler velocities and TM gradients. Results were compared with numerical modeling. TM gradients (range: 1.04-14.24 mmHg) consisted of 74.0 +/- 11.0% inertial forcers (range: 0.6-12.9 mmHg). Multivariate analysis predicted Mdv/dt = -4.171(S/D (RATIO)) + 0.063(LAvolume-max) + 5. Using this equation, a strong relationship was obtained for the clinical dataset (y=0.98x - 0.045, r=0.90) and the results of numerical modeling (y=0.96x - 0.16, r=0.84). TM gradients are mainly inertial and, as validated by modeling, can be estimated with echocardiography.

NASA Discipline Cardiopulmonary

Five years' central Pacific sea level from in situ array, satellite altimeter and numerical model

Temporal and spatial fetures of central equatorial Pacific Ocean sea-level variation appear similar, in measurements from two very different systems (one in the ocean and one carried on a satellite), and in results from a numerical model of the region. In particular, there is an interannual cycle: during El Nino, Kelvin waves appear at the equator, and the sea-surface ridge associated with the equatorial current system shifts southward; in non-El Nino years, instability waves appear at 6 deg N (strongest around the end of each calendar year), and the ridge shifts to the north. This three-way comparison gives support to both measurement systems and to the numerical model.

Donohue, Kathleen A.

Numerical Modeling of Boiling in a Heated Tube

This paper presents numerical models of boiling in a heated tube using the Generalized Fluid System Simulation Program (GFSSP), a finite-volume-based general-purpose flow network code developed at NASA/Marshall Space Flight Center. The heated tube is discretized into a one-dimensional array of nodes and branches to represent the flow of liquid and vapor in a tube with a prescribed pressure differential. The solid wall is also discretized into solid nodes and conductors to allow for heat transfer between the wall and the fluid. The conservation equations of mass, momentum, and energy of the fluid are solved simultaneously with the energy conservation equation for the solid wall. Two experimental configurations of fluid flowing in a vertical tube have been simulated, one with water and the other with liquid hydrogen. This paper compares experimental data with numerical predictions based on four different published correlations for boiling heat transfer coefficients. Three of these correlations are applicable to the saturated vertical flow conditions of the experiments. One of them is applicable to film boiling and has been used for the liquid hydrogen experiment, which was in film boiling regime. For the case of boiling water, the predictions of wall temperatures using the boiling heat transfer correlations agreed well with the experimental results. However, in the case of boiling hydrogen larger discrepancies were observed between the experimental data and numerical predictions.

Majumdar, Alok

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.

Numerical modeling of D-mappings with applications to chemical kinetics

Numerical modeling of D-mappings was studied and applied to solving nonlinear stiff systems. These mappings were locally linearized for convergence analysis, and some applications were made to chemical kinetics. The technique avoids using multistep implicit codes that require inversion of Jacobian matrices, but depends on the Jacobians for its convergence analysis.

Dey, S. K.