Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation codes”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Three-Dimensional Simulation of Traveling-Wave Tube Cold-Test Characteristics Using MAFIA

The three-dimensional simulation code MAFIA was used to compute the cold-test parameters - frequency-phase dispersion, beam on-axis interaction impedance, and attenuation - for two types of traveling-wave tube (TWT) slow-wave circuits. The potential for this electromagnetic computer modeling code to reduce the time and cost of TWT development is demonstrated by the high degree of accuracy achieved in calculating these parameters. Generalized input files were developed for ferruled coupled-cavity and TunneLadder slow-wave circuits. These files make it easy to model circuits of arbitrary dimensions. The utility of these files was tested by applying each to a specific TWT slow-wave circuit and comparing the results with experimental data. Excellent agreement was obtained.

Kory, Carol L.↗
NASA NTRS1995

Importance of real gas effects in entry technology

Critical issues concerning real gas effects in hypervelocity flows where the thermochemical nonequilibrium effects are pronounced, are described. Emphasis is on the development and validation of benchmark analysis tools. Approaches to develop and/or enhance phenomenological models, including test techniques to acquire databases, and incorporate them into computational flow field simulation codes are described. The performance characteristics of shock tubes, ballistic ranges and arc heated wind tunnels, state of the art diagnostics are included. The following aerothermodynamic phenomena are discussed: aerodynamic parameters, viscous interactions, turbulent transition, forebody heating/heat transfer, radiative heating, Lee-Base flows and low density real gas.

Deiwert, George S.↗
NASA NTRS1994

Implosive accretion and outbursts of active galactic nuclei

A model and simulation code have been developed for time-dependent axisymmetric disk accretion onto a compact object including for the first time the influence of an ordered magnetic field. The accretion rate and radiative luminosity of the disk are naturally coupled to the rate of outflow of energy and angular momentum in magnetically driven (+/- z) winds. The magnetic field of the wind is treated in a phenomenological way suggested by self-consistent wind solutions. The radial accretion speed u(r, t) of the disk matter is shown to be the sum of the usual viscous contribution and a magnetic contribution proportional to r(exp 3/2)(B(sub p exp 2))/sigma, where B(sub p)(r,t) is the poloidal field threading the disk and sigma(r,t) is the disk's surface mass density. An enhancement or variation in B(sub p) at a large radial distance leads to the formation of a soliton-like structure in the disk density, temperature, and B-field which propagates implosively inward. The implosion gives a burst in the power output in winds or jets and a simultaneous burst in the disk radiation. The model is pertinent to the formation of discrete fast-moving components in jets observed by very long baseline interferometry. These components appear to originate at times of optical outbursts of the active galactic nucleus.

Lovelace, R. V. E.↗
NASA NTRS1994

Quantitative measurements of vaporization, burst ionization, and emission characteristics of shaped charge barium releases

Intensity-calibrated color video recordings of three barium-shaped charge injections in the ionopshere were used to determine the initial ionization, the column density corresponding to unity optical depth, and the yield of vaporized barium in the fast jet. It was found that the initial ionization at the burst was less than 1% and that 0% burst ionization was consistent with the observations. Owing to the Doppler shift, the column density for optical thickness in the neutral barium varies somewhat according to the velocity distribution. For the cases examined here, the column density was 2-5 x 10(exp 10) atoms/sq cm. This value, which occurred 12 to 15 s after release, should be approximately valid for most shaped charge experiments. The yield was near 30% (15% in the fast jet) for two of the releases and was somewhat lower in the third, which also had a lower peak velocity. This study also demonstrated the applicability of the computer simulation code developed for chemical releases by Stenbaek-Nielsen and provided experimental verification of the Doppler-corrected emission rates calculated b Stenbaek-Nielsen (1989).

Hoch, Edward L.↗
NASA NTRS1994

Cable Discharge System for fundamental detonator studies

Sandia National Laboratories has recently completed the modification and installation of a cable discharge system (CDS) which will be used to study the physics of exploding bridgewire (EBW) detonators and exploding foil initiators (EFI or slapper). Of primary interest are the burst characteristics of these devices when subjected to the constant current pulse delivered by this system. The burst process involves the heating of the bridge material to a conductive plasma and is essential in describing the electrical properties of the bridgewire foil for use in diagnostics or computer models. The CDS described herein is capable of delivering up to an 8000 A pulse of 3 micron duration. Experiments conducted with the CDS to characterize the EBW and EFI burst behavior are also described. In addition, the CDS simultaneous VISAR capability permits updating the EFI electrical Gurney analysis parameters used in our computer simulation codes. Examples of CDS generated data for a typical EFI and EBW detonator are provided.

Peevy, Gregg R.↗
NASA NTRS1996

Active Sonic Boom Control

A theory and simulation code are developed to study nonsteady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of flight path. An aircraft in nonsteady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender nonsteady sources but requires no far field approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two nonsteady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and convert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on the average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical

Crow, Steven C.↗
NASA NTRS1997

Validation of an Accurate Three-Dimensional Helical Slow-Wave Circuit Model

The helical slow-wave circuit embodies a helical coil of rectangular tape supported in a metal barrel by dielectric support rods. Although the helix slow-wave circuit remains the mainstay of the traveling-wave tube (TWT) industry because of its exceptionally wide bandwidth, a full helical circuit, without significant dimensional approximations, has not been successfully modeled until now. Numerous attempts have been made to analyze the helical slow-wave circuit so that the performance could be accurately predicted without actually building it, but because of its complex geometry, many geometrical approximations became necessary rendering the previous models inaccurate. In the course of this research it has been demonstrated that using the simulation code, MAFIA, the helical structure can be modeled with actual tape width and thickness, dielectric support rod geometry and materials. To demonstrate the accuracy of the MAFIA model, the cold-test parameters including dispersion, on-axis interaction impedance and attenuation have been calculated for several helical TWT slow-wave circuits with a variety of support rod geometries including rectangular and T-shaped rods, as well as various support rod materials including isotropic, anisotropic and partially metal coated dielectrics. Compared with experimentally measured results, the agreement is excellent. With the accuracy of the MAFIA helical model validated, the code was used to investigate several conventional geometric approximations in an attempt to obtain the most computationally efficient model. Several simplifications were made to a standard model including replacing the helical tape with filaments, and replacing rectangular support rods with shapes conforming to the cylindrical coordinate system with effective permittivity. The approximate models are compared with the standard model in terms of cold-test characteristics and computational time. The model was also used to determine the sensitivity of various circuit parameters including typical manufacturing dimensional tolerances and support rod permittivity. By varying the circuit parameters of an accurate model using MAFIA, these sensitivities can be computed for manufacturing concerns, and design optimization previous to fabrication, thus eliminating the need for costly experimental iterations. Several variations were made to a standard helical circuit using MAFIA to investigate the effect that variations on helical tape and support rod width, metallized loading height and support rod permittivity, have on TWT cold-test characteristics.

Kory, Carol L.↗
NASA NTRS1997

Efficient Parallel Kernel Solvers for Computational Fluid Dynamics Applications

Distributed-memory parallel computers dominate today's parallel computing arena. These machines, such as Intel Paragon, IBM SP2, and Cray Origin2OO, have successfully delivered high performance computing power for solving some of the so-called "grand-challenge" problems. Despite initial success, parallel machines have not been widely accepted in production engineering environments due to the complexity of parallel programming. On a parallel computing system, a task has to be partitioned and distributed appropriately among processors to reduce communication cost and to attain load balance. More importantly, even with careful partitioning and mapping, the performance of an algorithm may still be unsatisfactory, since conventional sequential algorithms may be serial in nature and may not be implemented efficiently on parallel machines. In many cases, new algorithms have to be introduced to increase parallel performance. In order to achieve optimal performance, in addition to partitioning and mapping, a careful performance study should be conducted for a given application to find a good algorithm-machine combination. This process, however, is usually painful and elusive. The goal of this project is to design and develop efficient parallel algorithms for highly accurate Computational Fluid Dynamics (CFD) simulations and other engineering applications. The work plan is 1) developing highly accurate parallel numerical algorithms, 2) conduct preliminary testing to verify the effectiveness and potential of these algorithms, 3) incorporate newly developed algorithms into actual simulation packages. The work plan has well achieved. Two highly accurate, efficient Poisson solvers have been developed and tested based on two different approaches: (1) Adopting a mathematical geometry which has a better capacity to describe the fluid, (2) Using compact scheme to gain high order accuracy in numerical discretization. The previously developed Parallel Diagonal Dominant (PDD) algorithm and Reduced Parallel Diagonal Dominant (RPDD) algorithm have been carefully studied on different parallel platforms for different applications, and a NASA simulation code developed by Man M. Rai and his colleagues has been parallelized and implemented based on data dependency analysis. These achievements are addressed in detail in the paper.

Sun, Xian-He↗
NASA NTRS1996

Excitation of Plasma Waves in Aurora by Electron Beams

In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

daSilva, C. E.↗
NASA NTRS1997

Current Collection in a Magnetic Field

It is found that the upper-bound limit for current collection in the case of strong magnetic field from the current is close to that given by the Parker-Murphy formula. This conclusion is consistent with the results obtained in laboratory experiments. This limit weakly depends on the shape of the wire. The adiabatic limit in this case will be easily surpassed due to strong magnetic field gradients near the separatrix. The calculations can be done using the kinetic equation in the drift approximation. Analytical results are obtained for the region where the Earth's magnetic field is dominant. The current collection can be calculated (neglecting scattering) using a particle simulation code. Dr. Singh has agreed to collaborate, allowing the use of his particle code. The code can be adapted for the case when the current magnetic field is strong. The needed dm for these modifications is 3-4 months. The analytical description and essential part of the program is prepared for the calculation of the current in the region where the adiabatic description can be used. This was completed with the collaboration of Drs. Khazanov and Liemohn. A scheme of measuring the end body position is also proposed. The scheme was discussed in the laboratory (with Dr. Stone) and it was concluded that it can be proposed for engineering analysis.

Krivorutsky, E. N.↗
NASA NTRS1998

Analysis of Gamma-Ray Data from Solar Flares in Cycles 21 and 22

One of our primary accomplishments under grant NAGW-35381 was the systematic derivation and compilation, for the first time, of physical parameters for all gamma-ray flares detected by the SMM GRS during its ten year lifetime. The flare parameters derived from the gamma-ray spectra include: bremsstrahlung fluence and best-fit power-law parameters, narrow nuclear line fluence, positron annihilation line fluence, neutron capture line fluence, and an indication of whether or not greater than 10 MeV emissions were present. We combined this compilation of flare parameters with our plots of counting rate time histories and flare spectra to construct an atlas of gamma-ray flare characteristics. The atlas time histories display four energy bands: 56-199 kev, 298526 keV, 4-8 MeV, and 10-25 MeV. These energy bands respectively measure nonrelativistic bremsstrahlung, trans-relativistic bremsstrahlung, nuclear de-excitation, and ultra-relativistic bremsstrahlung. The atlas spectra show the integrated high-energy spectra measured for all GRS flares and dissects them into electron bremsstrahlung, positron annihilation and nuclear emission components. The atlas has been accepted for publication in the Astrophysical Journal Supplements and is currently in press. The atlas materials were also supplied to the Solar Data Analysis Center at Goddard Space Flight Center and were made available through a web site at the University of New Hampshire. Since a uniform methodology was adopted for deriving the flare parameters, this atlas will be very useful for future statistical and correlative studies of solar flares-three independent groups are presently using it to correlate interplanetary energetic particle measurements with our gamma-ray measurements. A better model for the response of the GRS instrument to high energy radiation was also developed. A refined response model was needed because the old model was not adequate for predicting the first and second escape peaks associated with strong nuclear lines nor could it accurately describe the Compton continuum shape. The new response was developed using a GEANT based simulation code and tested against preflight calibration data. The refinement of the response model and the removal of systematic errors now allow more detailed spectral studies of the GRS gamma-ray measurements. This refined response function was supplied to the Solar DAC at Goddard and was also made available via a web site at the University of New Hampshire.

Vestrand, W. Thomas↗
NASA NTRS1997

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling-Wave Tube Structures

Conventional methods used to measure the cold-test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit. It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves-several approximations. With the advent of accurate three- dimensional helical circuit models, these standard approximations can be fully investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the three-dimensional simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using three-dimensional computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1999

Probabilistic Structural Analysis Methods (PSAM) for Select Space Propulsion System Components

Probabilistic Structural Analysis Methods (PSAM) are described for the probabilistic structural analysis of engine components for current and future space propulsion systems. Components for these systems are subjected to stochastic thermomechanical launch loads. Uncertainties or randomness also occurs in material properties, structural geometry, and boundary conditions. Material property stochasticity, such as in modulus of elasticity or yield strength, exists in every structure and is a consequence of variations in material composition and manufacturing processes. Procedures are outlined for computing the probabilistic structural response or reliability of the structural components. The response variables include static or dynamic deflections, strains, and stresses at one or several locations, natural frequencies, fatigue or creep life, etc. Sample cases illustrates how the PSAM methods and codes simulate input uncertainties and compute probabilistic response or reliability using a finite element model with probabilistic methods.

Source record↗
NASA NTRS1998

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling Wave Tube Structures

Conventional methods used to measure the cold- test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. With the advent of accurate three-dimensional (3-D) helical circuit models, these standard approximations can be fully Investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the 3-D simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using 3-D computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1998

Development of a Model for the Night Side Magnetopause using Global Simulations

During the final year of this investigation we have finished several event studies that we considered necessary for the development of a tail magnetopause model and for the calibration of our simulation code. We have not reached the ultimate goal of the project, i.e., the development of an analytical tail magnetopause model. In the course of the investigation we have learned that such a model would be much more complex than we had anticipated. However, the investigations that we conducted towards this goal have led to significant results and discoveries that are of considerable value for understanding the tail magnetopause. These are summarized in the following sections.

Raeder, Joachim↗
NASA NTRS1998

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling-Wave Tube Structures

Conventional methods used to measure the cold-test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit. It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. With the advent of accurate three-dimensional (3-D) helical circuit models, these standard approximations can be fully investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the 3-D simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using 3-D computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1999

Multiphysics Simulation of Active Hypersonic Lip Cooling

This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.

Melis, Matthew E.↗