Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 433 records · Page 24

A numerical study of chemically reacting flow in nozzles

The Space Station uses small rocket motors, called thrusters, for orientation control. Because of the lack of viable design tools for small rockets, the initial thruster design was basically a very small version of a large rocket motor. Thrust measurements of the initial design were lower than predicted. To improve predictions it was decided to develop a version of the RPLUS2D reacting flow code for thruster calculations. RPLUS2D employs an implicit finite volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a coupled and very efficient manner. The combustion processes are modeled by a 9-species, 18 step finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model. The code is extended to handle multiple subsonic inlet conditions where the total mass flow is governed by conditions calculated at the thruster-throat. Results are shown for a thruster design where the overall mixture ratio is hydrogen rich. A calculation of a large area ratio divergent nozzle is also presented.

Vanoverbeke, Thomas J.↗

Nodal Modeling of Tank Pressurization and Draining using a Multi-Node-Ullage Approach

The purpose of the pressurization system in liquid rocket propulsion is to control the pressure in the gas space of the propellant tank (known as the ullage space) and the propellant mass flowrate to the engine. A mathematical model is required to predict the amount of pressurant necessary to ensure that pressure and temperature levels inside the tank remain within acceptable limits and that the propellant pressure leaving the tank satisfies the net positive suction pressure (NPSP) requirement of the pump feeding the engine. Nodal codes typically model tank pressurization and draining using a single node to represent the ullage and a single node to represent the propellant. As the tank drains, the ullage node grows and the propellant node shrinks. The heat transfer between ullage to wall and ullage to propellant is governed by natural convection. Designers of liquid propulsion systems often use empirical correlations to estimate the “Collapse Factor” which represents the ratio of pressurant required with heat transfer and the amount of pressurant required without heat transfer. A single node ullage model of tank pressurization was developed using GFSSP to compute the collapse factor reasonably well and later was used to model tank pressurization during test firing of the FASTRAC rocket engine. The predicted tank pressure compared well with the test data. In the early 1970’s, pressurization and drain tests with liquid methane were performed at NASA Lewis Research Center in a vacuum chamber. A 5 ft diameter spherical aluminum tank was tested to drain from 95% to 5% full using gaseous helium, hydrogen, nitrogen and methane as pressurant. Tests were conducted with different pressurant inlet temperatures and drain times. Measured data include pressurant requirement, amount of pressurant condensed, and ullage and wall temperatures at various heights in the ullage space at the end of draining. A single node GFSSP model was developed to simulate helium pressurization of the methane tank. Predicted helium consumption was 8-23% less than measured. The average error of the six test cases was 16%. A single ullage node with multiple solid node model was developed using Thermal Desktop. Predicted helium consumption compares with the test data within 2%. This paper describes the development of a GFSSP multi-node ullage model of the test configuration and compares the predicted pressurant consumption for both helium and autogenous pressurization using gaseous methane with experimental as well as TD predictions.

Nodal Model↗

Atmospheric Correction of Satellite Imagery Using Modtran 3.5 Code

When performing satellite remote sensing of the earth in the solar spectrum, atmospheric scattering and absorption effects provide the sensors corrupted information about the target's radiance characteristics. We are faced with the problem of reconstructing the signal that was reflected from the target, from the data sensed by the remote sensing instrument. This article presents a method for simulating radiance characteristic curves of satellite images using a MODTRAN 3.5 band model (BM) code to solve the radiative transfer equation (RTE), and proposes a method for the implementation of an adaptive system for automated atmospheric corrections. The simulation procedure is carried out as follows: (1) for each satellite digital image a radiance characteristic curve is obtained by performing a digital number (DN) to radiance conversion, (2) using MODTRAN 3.5 a simulation of the images characteristic curves is generated, (3) the output of the code is processed to generate radiance characteristic curves for the simulated cases. The simulation algorithm was used to simulate Landsat Thematic Mapper (TM) images for two types of locations: the ocean surface, and a forest surface. The simulation procedure was validated by computing the error between the empirical and simulated radiance curves. While results in the visible region of the spectrum where not very accurate, those for the infrared region of the spectrum were encouraging. This information can be used for correction of the atmospheric effects. For the simulation over ocean, the lowest error produced in this region was of the order of 105 and up to 14 times smaller than errors in the visible region. For the same spectral region on the forest case, the lowest error produced was of the order of 10-4, and up to 41 times smaller than errors in the visible region,

Gonzales, Fabian O.↗

The investigation of tethered satellite system dynamics

The tether control law to retrieve the satellite was modified in order to have a smooth retrieval trajectory of the satellite that minimizes the thruster activation. The satellite thrusters were added to the rotational dynamics computer code and a preliminary control logic was implemented to simulate them during the retrieval maneuver. The high resolution computer code for modelling the three dimensional dynamics of untensioned tether, SLACK3, was made fully operative and a set of computer simulations of possible tether breakages was run. The distribution of the electric field around an electrodynamic tether in vacuo severed at some length from the shuttle was computed with a three dimensional electrodynamic computer code.

Lorenzini, E.↗

Software for Collaborative Engineering of Launch Rockets

The Rocket Evaluation and Cost Integration for Propulsion and Engineering software enables collaborative computing with automated exchange of information in the design and analysis of launch rockets and other complex systems. RECIPE can interact with and incorporate a variety of programs, including legacy codes, that model aspects of a system from the perspectives of different technological disciplines (e.g., aerodynamics, structures, propulsion, trajectory, aeroheating, controls, and operations) and that are used by different engineers on different computers running different operating systems. RECIPE consists mainly of (1) ISCRM a file-transfer subprogram that makes it possible for legacy codes executed in their original operating systems on their original computers to exchange data and (2) CONES an easy-to-use filewrapper subprogram that enables the integration of legacy codes. RECIPE provides a tightly integrated conceptual framework that emphasizes connectivity among the programs used by the collaborators, linking these programs in a manner that provides some configuration control while facilitating collaborative engineering tradeoff studies, including design to cost studies. In comparison with prior collaborative-engineering schemes, one based on the use of RECIPE enables fewer engineers to do more in less time.

Stanley, Thomas Troy↗

A coupled radiative-convective-photochemical model of the stratosphere

A one-dimensional atmospheric model is described. The neutral photochemical processes of the troposphere and stratosphere are presented by a reaction set of 87 kinetic processes and 28 photolytic reactions, with vertical transport being taken into account by an eddy diffusion parameter. The vertical temperature profiles are coupled with the chemical calculations by means of a radiative-convective code. The model has the capability of including effects of multiple scattering and diurnal variations in species concentrations. The numerical solution methodology is described in detail, and current input data are given along with results obtained in a sample calculation.

Boughner, R. E.↗

Modelling of capillary optics as a focussing hard X-ray concentrator

The development of glass microcapillary arrays by Kumakhov and co-workers in the Soviet Union suggests that it is possible to construct novel devices for concentrating and focussing X-ray beams. Hypothetical focussing concentrators made up of capillary arrays have been configured to operate over a broad range of X-ray energy and its characteristics have been studied with a Monte Carlo ray tracing code. These modeling results indicate that the concentrators would have significant effective area up to 100 keV while their dimensions are similar to conventional grazing incidence telescopes that operate in the soft X-ray band.

Gorenstein, P.↗

Numerical Predictions of Damage and Failure in Carbon Fiber Reinforced Laminates Using a Thermodynamically-Based Work Potential Theory

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, referred to as enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Consistent characteristic lengths are introduced into the formulation to govern the evolution of the failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs are derived. The theory is implemented into a commercial finite element code. The model is verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared against the experimental results.

composite materials↗

An integrated optimum design approach for high speed prop rotors

The objective is to develop an optimization procedure for high-speed and civil tilt-rotors by coupling all of the necessary disciplines within a closed-loop optimization procedure. Both simplified and comprehensive analysis codes are used for the aerodynamic analyses. The structural properties are calculated using in-house developed algorithms for both isotropic and composite box beam sections. There are four major objectives of this study. (1) Aerodynamic optimization: The effects of blade aerodynamic characteristics on cruise and hover performance of prop-rotor aircraft are investigated using the classical blade element momentum approach with corrections for the high lift capability of rotors/propellers. (2) Coupled aerodynamic/structures optimization: A multilevel hybrid optimization technique is developed for the design of prop-rotor aircraft. The design problem is decomposed into a level for improved aerodynamics with continuous design variables and a level with discrete variables to investigate composite tailoring. The aerodynamic analysis is based on that developed in objective 1 and the structural analysis is performed using an in-house code which models a composite box beam. The results are compared to both a reference rotor and the optimum rotor found in the purely aerodynamic formulation. (3) Multipoint optimization: The multilevel optimization procedure of objective 2 is extended to a multipoint design problem. Hover, cruise, and take-off are the three flight conditions simultaneously maximized. (4) Coupled rotor/wing optimization: Using the comprehensive rotary wing code CAMRAD, an optimization procedure is developed for the coupled rotor/wing performance in high speed tilt-rotor aircraft. The developed procedure contains design variables which define the rotor and wing planforms.

Chattopadhyay, Aditi↗

Electron impact emission cross sections needed in modeling UV observations of the upper atmospheres of the planets

In recognition of the need for accurate collision cross-sections suitable for modeling spectroscopic observations of the planetary systems and calibrate flight instruments, NASA-JPL has determined the laboratory emission cross-sections of H, H2, N2, SO2, etc. Resonance excitation, predissociation, and non-Franck-Condon band intensity systems are determined by configuration interaction and must be taken into account in analyses of spacecraft observations. Attention is given to the application of the cross-sectional data obtained to both electron energy loss codes and models of Voyager, IUE, Galileo, and HST observations of the solar system.

Ajello, J. M.↗

Refactoring and Componentizing Legacy Codes: GMI Case Study

As Earth science applications grow increasingly complex over time, maintenance of such software becomes a significant challenge. In general, the software subsystems were designed, implemented and integrated (over an extended period of time) by different groups that did not follow the same software design standard. The objective was understandably to obtain scientific results quickly, but this goal was often achieved at the expense of good software practices. With that pattern of development, the cost of incremental changes to the applications grow ominously large. This increase in effort, tends to squeeze development resources and largely prevent any systematic attempt to reintroduce better software practices and thereby enhance the maintainability of the software. In this paper, we present the process we used to incrementally refactor and componentize the Global Modeling Initiative code to make it more understandable, maintainable, flexible and extensible. The resulting product not only preserves the scientific integrity of the original code but also extends the capabilities of the code for several applications.

Kouatchou, Jules↗

Modeling of the Convection and Interaction of Ring Current, Plasmaspheric and Plasma Sheet Plasmas in the Inner Magnetosphere

Distinctive sources of ions reside in the plasmasphere, plasmasheet, and ring current regions at discrete energies constitute the major plasma populations in the inner/middle magnetosphere. They contribute to the electrodynamics of the ionosphere-magnetosphere system as important carriers of the global current system, in triggering; geomagnetic storm and substorms, as well as critical components of plasma instabilities such as reconnection and Kelvin-Helmholtz instability at the magnetospheric boundaries. Our preliminary analysis of in-situ measurements shoves the complexity of the plasmas pitch angle distributions at particularly the cold and warm plasmas, vary dramatically at different local times and radial distances from the Earth in response to changes in solar wind condition and Dst index. Using an MHD-ring current coupled code, we model the convection and interaction of cold, warm and energetic ions of plasmaspheric, plasmasheet, and ring current origins in the inner magnetosphere. We compare our simulation results with in-situ and remotely sensed measurements from recent instrumentation on Geotail, Cluster, THEMIS, and TWINS spacecraft.

Fok, Mei-Ching↗

Wind-US Results for the AIAA 2nd Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 2nd Propulsion Aerodynamics Workshop. The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a 25deg conical nozzle for a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which the 25deg conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume shock (NPR=4.0). A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of these cases across various codes and models. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations were examined.

Wind-US↗

A Numerical Investigation of the Extinction of Low Strain Rate Diffusion Flames by an Agent in Microgravity

Our goal has been to investigate the influence of both dilution and radiation on the extinction process of nonpremixed flames at low strain rates. Simulations have been performed by using a counterflow code and three radiation models have been included in it, namely, the optically thin, the narrowband, and discrete ordinate models. The counterflow flame code OPPDIFF was modified to account for heat transfer losses by radiation from the hot gases. The discrete ordinate method (DOM) approximation was first suggested by Chandrasekhar for solving problems in interstellar atmospheres. Carlson and Lathrop developed the method for solving multi-dimensional problem in neutron transport. Only recently has the method received attention in the field of heat transfer. Due to the applicability of the discrete ordinate method for thermal radiation problems involving flames, the narrowband code RADCAL was modified to calculate the radiative properties of the gases. A non-premixed counterflow flame was simulated with the discrete ordinate method for radiative emissions. In comparison with two other models, it was found that the heat losses were comparable with the optically thin and simple narrowband model. The optically thin model had the highest heat losses followed by the DOM model and the narrow-band model.

Puri, Ishwar K.↗

Pre-engineering Spaceflight Validation of Environmental Models and the 2005 HZETRN Simulation Code

The HZETRN code has been identified by NASA for engineering design in the next phase of space exploration highlighting a return to the Moon in preparation for a Mars mission. In response, a new series of algorithms beginning with 2005 HZETRN, will be issued by correcting some prior limitations and improving control of propagated errors along with established code verification processes. Code validation processes will use new/improved low Earth orbit (LEO) environmental models with a recently improved International Space Station (ISS) shield model to validate computational models and procedures using measured data aboard ISS. These validated models will provide a basis for flight-testing the designs of future space vehicles and systems of the Constellation program in the LEO environment.

Nealy, John E.↗

Mixed-field GCR Simulations for Radiobiological Research Using Ground Based Accelerators

Space radiation is comprised of a large number of particle types and energies, which have differential ionization power from high energy protons to high charge and energy (HZE) particles and secondary neutrons produced by galactic cosmic rays (GCR). Ground based accelerators such as the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL) are used to simulate space radiation for radiobiology research and dosimetry, electronics parts, and shielding testing using mono-energetic beams for single ion species. As a tool to support research on new risk assessment models, we have developed a stochastic model of heavy ion beams and space radiation effects, the GCR Event-based Risk Model computer code (GERMcode). For radiobiological research on mixed-field space radiation, a new GCR simulator at NSRL is proposed. The NSRL-GCR simulator, which implements the rapid switching mode and the higher energy beam extraction to 1.5 GeV/u, can integrate multiple ions into a single simulation to create GCR Z-spectrum in major energy bins. After considering the GCR environment and energy limitations of NSRL, a GCR reference field is proposed after extensive simulation studies using the GERMcode. The GCR reference field is shown to reproduce the Z and LET spectra of GCR behind shielding within 20% accuracy compared to simulated full GCR environments behind shielding. A major challenge for space radiobiology research is to consider chronic GCR exposure of up to 3-years in relation to simulations with cell and animal models of human risks. We discuss possible approaches to map important biological time scales in experimental models using ground-based simulation with extended exposure of up to a few weeks and fractionation approaches at a GCR simulator.

Kim, Myung-Hee Y.↗

Dynamical experiments on models of colliding disk galaxies

Collisions between galaxies can induce large morphological changes in the participants and, in the case of colliding disk galaxies, bridges and tails are often formed. Observations of such systems indicate a wide variation in color (see Larson and Tinsley, 1978) and that some of the particpants are experiencing enhanced rates of star formation, especially in their central regions (Bushouse 1986, 1987; Kennicutt et al., 1987, Bushouse, Lamb, and Werner, 1988). Here the authors describe progress made in understanding some of the dynamics of interacting galaxies using N-body stellar dynamical computer experiments, with the goal of extending these models to include a hydrodynamical treatment of the gas so that a better understanding of globally enhanced star formation will eventually be forthcoming. It was concluded that close interactions between galaxies can produce large perturbations in both density and velocity fields. The authors measured, via computational experiments that represent a galaxy's stars, average radial velocity flows as large as 100 km/sec and 400 percent density increases. These can occur in rings that move outwards through the disk of a galaxy, in roughly homologous inflows toward the nucleus, and in off center, non-axisymmetric regions. Here the authors illustrate where the gas is likely to flow during the early stages of interaction and in future work they plan to investigate the fate of the gas more realistically by using an N-body/Smoothed Particle Hydrodynamics code to model both the stellar and gaseous components of a disk galaxy during a collision. Specifically, they will determine the locations of enhanced gas density and the strength and location of shock fronts that form during the interaction.

Gerber, Richard A.↗

Overview of the Graphical User Interface for the GERMcode (GCR Event-Based Risk Model)

The descriptions of biophysical events from heavy ions are of interest in radiobiology, cancer therapy, and space exploration. The biophysical description of the passage of heavy ions in tissue and shielding materials is best described by a stochastic approach that includes both ion track structure and nuclear interactions. A new computer model called the GCR Event-based Risk Model (GERM) code was developed for the description of biophysical events from heavy ion beams at the NASA Space Radiation Laboratory (NSRL). The GERMcode calculates basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at NSRL for the purpose of simulating space radiobiological effects. For mono-energetic beams, the code evaluates the linear-energy transfer (LET), range (R), and absorption in tissue equivalent material for a given Charge (Z), Mass Number (A) and kinetic energy (E) of an ion. In addition, a set of biophysical properties are evaluated such as the Poisson distribution of ion or delta-ray hits for a specified cellular area, cell survival curves, and mutation and tumor probabilities. The GERMcode also calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle. The contributions from primary ion and nuclear secondaries are evaluated. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections, and has been used by the GERMcode for application to thick target experiments. The GERMcode provides scientists participating in NSRL experiments with the data needed for the interpretation of their experiments, including the ability to model the beam line, the shielding of samples and sample holders, and the estimates of basic physical and biological outputs of the designed experiments. We present an overview of the GERMcode GUI, as well as providing training applications.

Kim, Myung-Hee Y.↗