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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A study of slow-mode structures in the dayside magnetosheath

Recent observations indicate that a region of enhanced plasma pressure and decreased magnetic field intensity frequently occurs in front of the plasma depletion layer at the dayside magnetopause (Song et al., 1990). This inverse relationship is characteristic of a slow-mode wave. This phenomenon was simulated with a two-dimensional incompressible MHD simulation code. When a normal component of the interplanetary magnetic field is present (Bx not equal to 0), the total magnetic field intensity tends to decrease in front of the depletion layer due to the bending of the magnetic field lines, and the plasma pressure is enhanced in this region. On the other hand, when Bx = 0, this slow-mode structure is not present in the simulation and only the plasma depletion layer is observed.

Lee, L. C.↗
NASA NTRS1998

Results of Two-Stage Light-Gas Gun Development Efforts and Hypervelocity Impact Tests of Advanced Thermal Protection Materials

Gun development efforts to increase the launching capabilities of the NASA Ames 0.5-inch two-stage light-gas gun have been investigated. A gun performance simulation code was used to guide initial parametric variations and hardware modifications, in order to increase the projectile impact velocity capability to 8 km/s, while maintaining acceptable levels of gun barrel erosion and gun component stresses. Concurrent with this facility development effort, a hypervelocity impact testing series in support of the X-33/RLV program was performed in collaboration with Rockwell International. Specifically, advanced thermal protection system materials were impacted with aluminum spheres to simulate impacts with on-orbit space debris. Materials tested included AETB-8, AETB-12, AETB-20, and SIRCA-25 tiles, tailorable advanced blanket insulation (TABI), and high temperature AFRSI (HTA). The ballistic limit for several Thermal Protection System (TPS) configurations was investigated to determine particle sizes which cause threshold TPS/structure penetration. Crater depth in tiles was measured as a function of impact particle size. The relationship between coating type and crater morphology was also explored. Data obtained during this test series was used to perform a preliminary analysis of the risks to a typical orbital vehicle from the meteoroid and space debris environment.

Cornelison, C. J.↗
NASA NTRS1995

A study of workstation computational performance for real-time flight simulation

With recent advances in microprocessor technology, some have suggested that modern workstations provide enough computational power to properly operate a real-time simulation. This paper presents the results of a computational benchmark, based on actual real-time flight simulation code used at Langley Research Center, which was executed on various workstation-class machines. The benchmark was executed on different machines from several companies including: CONVEX Computer Corporation, Cray Research, Digital Equipment Corporation, Hewlett-Packard, Intel, International Business Machines, Silicon Graphics, and Sun Microsystems. The machines are compared by their execution speed, computational accuracy, and porting effort. The results of this study show that the raw computational power needed for real-time simulation is now offered by workstations.

Maddalon, Jeffrey M.↗
NASA NTRS2001

Integration of the QMSFRG Database into the HZETRN Code

Accurate nuclear interaction data bases are needed for describing the transport of space radiation in matter including space craft structures, atmospheres, and tissues. Transport models support the identification and development of new material concepts for human and electronic part protection. Quantum effects are manifested in nuclear reactions in several ways including interference effects between terms in the multiple scattering series, the many-body nuclear wave functions (for e.g. the roles of shell structure and Fermi momentum) and nuclear clustering. The quantum multiple scattering fragmentation model (QMSFRG) is a comprehensive model for generating nuclear interaction databases for galactic cosmic ray (GCR) transport. Other nuclear databases including the NUCFRG model and Monte-Carlo simulation codes such as FLUKA, LAHET, HETC, and GEANT ignore quantum effects. These codes fail to describe many important features of nuclear reactions and are thus inaccurate for the evaluation of materials for radiation protection. Previously we have shown that quantum effects are manifested through constructive interference in forward production spectra, the effects of Fermi momentum on production spectra, cluster nuclei knockout, and the nuclear response function. Quantum effects are especially important for heavy ions with mass numbers less than 20 that dominate radiation transport in human tissues and for the materials that are expected to be superior in space radiation protection. We describe the integration of the QMSFRG model into the HZETRN transport code. Integration milestones include proper treatment of odd-even charge-mass effects in nuclear fragmentation and the momentum distribution of nucleon production from GCR primary heavy ions. We have also modified the two-body amplitudes in the model to include nuclear medium effects. In order to include a comprehensive description of the GCR isotopic composition in materials, we have described the isotopic composition of the GCR by extending the 59-isotope version of HZETRN to an 120-isotope version. The isotopic composition of most primary GCR elements (including H, He, C, N, O, Ne, Mg, Si, Ar, Ca, Cr, and Fe) are included in the extended model. We discuss results for the high-energy neutron composition inside materials, and the charge and mass distribution for benchmark GCR problems.

Cucinotta, F. A.↗
NASA NTRS1987

A kinetic study of solar wind mass loading and cometary bow shocks

The kinetic processes associated with solar wind mass loading due to pickup of cometary ions, and the formation of cometary bow shocks were investigated numerically using a hybrid simulation code described by Leroy et al. (1981), in which solar wind protons and heavy cometary ions are treated kinetically, but the electrons are treated as a massless fluid. It was found that the solar wind decelaration and pickup of cometary ions take place through both the macroscopic electromagnetic fields embedded in the solar wind and the microscopic field associated with low-frequency electromagnetic waves that are generated by the unstable velocity distribution function of the cometary ions. The results of the main simulation runs at various cone angles are described and compared with the recent observations of comets Giacobini-Zinner and Halley.

Omidi, N.↗
NASA NTRS1995

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗
NASA NTRS2021

Improvements to the Flight Analysis and Simulation Tool (FAST) and Initial Development of the Genesis Flight Mechanics Simulation for Ascent, Aerocapture, Entry, Descent, and Landing (A2EDL) Trajectory Design

The NASA Engineering and Safety Center (NESC) was requested to provide significant enhancements to the Flight Analysis and Simulation Tool (FAST), a generic, variable-degree-of-freedom, multi-body ascent, aerocapture, entry, descent, and landing (A2EDL) flight simulation code and a key Agency analysis tool. This report details the upgrades to FAST that were completed, describes a change to the assessment’s scope that allowed its goals to be exceeded with a new and innovative approach, and discusses the results of the modified approach.

NASA Engineering and Safety Center; Flight Analysi↗
NASA NTRS2011

Numerical Study Comparing RANS and LES Approaches on a Circulation Control Airfoil

A numerical study over a nominally two-dimensional circulation control airfoil is performed using a large-eddy simulation code and two Reynolds-averaged Navier-Stokes codes. Different Coanda jet blowing conditions are investigated. In addition to investigating the influence of grid density, a comparison is made between incompressible and compressible flow solvers. The incompressible equations are found to yield negligible differences from the compressible equations up to at least a jet exit Mach number of 0.64. The effects of different turbulence models are also studied. Models that do not account for streamline curvature effects tend to predict jet separation from the Coanda surface too late, and can produce non-physical solutions at high blowing rates. Three different turbulence models that account for streamline curvature are compared with each other and with large eddy simulation solutions. All three models are found to predict the Coanda jet separation location reasonably well, but one of the models predicts specific flow field details near the Coanda surface prior to separation much better than the other two. All Reynolds-averaged Navier-Stokes computations produce higher circulation than large eddy simulation computations, with different stagnation point location and greater flow acceleration around the nose onto the upper surface. The precise reasons for the higher circulation are not clear, although it is not solely a function of predicting the jet separation location correctly.

Rumsey, Christopher L.↗
NASA NTRS1993

Hybrid simulations of thin current sheets

A one-dimensional, hybrid simulation code is used to study current sheets with a nonzero normal magnetic field B(sub z) and a dawn-to-dusk electric field E(sub y). Such configurations are dependent upon only two parameters: we use the normalized normal magnetic field B-normalized (sub z) = B(sub z)/(4(pi)(n(sub b)) (v(exp 2 sub T))(exp 1/2) and normalized electric field V-normalized (sub D) = (1/V(sub T)(cE(sub y)/B(sub z)), where V(sub T) is the thermal velocity of ions prior to their interaction with the current sheet and n(sub b) is the number density outside the current sheet (at the simulation boundary). A third parameter that is relevant to the motion of particles in current sheets is kappa(sub A), the value of kappa = (R(sub min)/rho(sub max))(exp 1/2) for particles of average energy. We find that if either B-normalized (sub z) is close to or greater than 1, or if kappa(sub A) is close to 1, a rotational mode develops in which the z = 0 current rotates with the ion sense about the normal magnetic field, while for small values of both B-normalized (sub z) or kappa(sub A), the configuration is quasi-steady. To achieve values of kappa(sub A) of the order of or larger than 1, we decrease the value of V-normalized (sub D) uniformly. We find that the magnetic field fluctuations and particle distribution functions are similar in many respects to what was observed in the day 240, 1986, Active Magnetospheric Particle Tracer Explorer (AMPTE)/CCE current disruption event, an event that appears to be located at the site of initiation of current disruption and related particle energization.

Burkhart, G. R.↗
NASA NTRS1997

Numerical Assessment of Four-Port Through-Flow Wave Rotor Cycles with Passage Height Variation

The potential for improved performance of wave rotor cycles through the use of passage height variation is examined. A Quasi-one-dimensional CFD code with experimentally validated loss models is used to determine the flowfield in the wave rotor passages. Results indicate that a carefully chosen passage height profile can produce substantial performance gains. Numerical performance data are presented for a specific profile, in a four-port, through-flow cycle design which yielded a computed 4.6% increase in design point pressure ratio over a comparably sized rotor with constant passage height. In a small gas turbine topping cycle application, this increased pressure ratio would reduce specific fuel consumption to 22% below the un-topped engine; a significant improvement over the already impressive 18% reductions predicted for the constant passage height rotor. The simulation code is briefly described. The method used to obtain rotor passage height profiles with enhanced performance is presented. Design and off-design results are shown using two different computational techniques. The paper concludes with some recommendations for further work.

Paxson, D. E.↗
NASA NTRS1993

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS1994

Novel high-gain, improved-bandwidth, finned-ladder V-band Traveling-Wave Tube slow-wave circuit design

The V-band frequency range of 59-64 GHz is a region of the millimeter-wave spectrum that has been designated for inter-satellite communications. As a first effort to develop a high-efficiency V-band Traveling-Wave Tube (TWT), variations on a ring-plane slow-wave circuit were computationally investigated to develop an alternative to the more conventional ferruled coupled-cavity circuit. The ring-plane circuit was chosen because of its high interaction impedance, large beam aperture, and excellent thermal dissipation properties. Despite these advantages, however, low bandwidth and high voltage requirements have, until now, prevented its acceptance outside the laboratory. In this paper, the three-dimensional electrodynamic simulation code MAFIA (solution of MAxwell's Equation by the Finite-Integration-Algorithm) is used to investigate methods of increasing the bandwidth and lowering the operating voltage of the ring-plane circuit. Calculations of frequency-phase dispersion, beam on-axis interaction impedance, attenuation and small-signal gain per wavelength were performed for various geometric variations and loading distributions of the ring-plane TWT slow-wave circuit. Based on the results of the variations, a circuit termed the finned-ladder TWT slow-wave circuit was designed and is compared here to the scaled prototype ring-plane and a conventional ferruled coupled-cavity TWT circuit over the V-band frequency range. The simulation results indicate that this circuit has a much higher gain, significantly wider bandwidth, and a much lower voltage requirement than the scaled ring-plane prototype circuit, while retaining its excellent thermal dissipation properties. The finned-ladder circuit has a much larger small-signal gain per wavelength than the ferruled coupled-cavity circuit, but with a moderate sacrifice in bandwidth.

Kory, Carol L.↗
NASA NTRS1993

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS2010

Calculation of Dose Deposition in 3D Voxels by Heavy Ions

The biological response to high-LET radiation is very different from low-LET radiation, and can be partly attributed to the energy deposition by the radiation. Several experiments, notably detection of gamma-H2AX foci by immunofluorescence, has revealed important differences in the nature and in the spatial distribution of double-strand breaks (DSB) induced by low- and high-LET radiations. Many calculations, most of which are based on amorphous track models with radial dose, have been combined with chromosome models to calculate the number and distribution of DSB within nuclei and chromosome aberrations. In this work, the Monte-Carlo track structure simulation code RITRACKS have been used to calculate directly the energy deposition in voxels (3D pixels). A cubic volume of 5 micrometers of side was irradiated by 1) 450 (1)H+ ions of 300 MeV (LET is approximately 0.3 keV/micrometer) and 2) by 1 (56)Fe26+ ion of 1 GeV/amu (LET is approximately 150 keV/micrometer). In both cases, the dose deposited in the volume is approximately 1 Gy. All energy deposition events are recorded and dose is calculated in voxels of 20 micrometers of side. The voxels are then visualized in 3D by using a color scale to represent the intensity of the dose in a voxel. This simple approach has revealed several important points which may help understand experimental observations. In both simulations, voxels which receive low dose are the most numerous, and those corresponding to electron track ends received a dose which is in the higher range. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. The distribution of the voxels shows major differences for the (56)Fe26+ ion. The track structure can still be seen, and voxels with much higher dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and may be responsible for DSB that are more difficult to repair. By applying a threshold on the dose visualization, voxels corresponding to electron track ends are evidenced and the spatial distribution of voxels is very similar to the distribution of DSB observed in gamma H2AX experiments, even if no chromosomes have been included in the simulation. Furthermore, this work has shown that a significant dose is deposited in voxels corresponding to electron track ends. Since some delta-rays from iron ion can travel several millimeters, they may also be of radiobiological importance.

Plante, Ianik↗
NASA NTRS2003

A Nonlinear Model for Fuel Atomization in Spray Combustion

Most gas turbine combustion codes rely on ad-hoc statistical assumptions regarding the outcome of fuel atomization processes. The modeling effort proposed in this project is aimed at developing a realistic model to produce accurate predictions of fuel atomization parameters. The model involves application of the nonlinear stability theory to analyze the instability and subsequent disintegration of the liquid fuel sheet that is produced by fuel injection nozzles in gas turbine combustors. The fuel sheet is atomized into a multiplicity of small drops of large surface area to volume ratio to enhance the evaporation rate and combustion performance. The proposed model will effect predictions of fuel sheet atomization parameters such as drop size, velocity, and orientation as well as sheet penetration depth, breakup time and thickness. These parameters are essential for combustion simulation codes to perform a controlled and optimized design of gas turbine fuel injectors. Optimizing fuel injection processes is crucial to improving combustion efficiency and hence reducing fuel consumption and pollutants emissions.

Liu, Nan-Suey↗
NASA NTRS1991

Evolution of magnetic flux ropes associated with flux transfer events and interplanetary magnetic clouds

Spacecraft observations suggest that flux transfer events and interplanetary magnetic clouds may be associated with magnetic flux ropes which are magnetic flux tubes containing helical magnetic field lines. In the magnetic flux ropes, the azimuthal magnetic field is superposed on the axial field. The time evolution of a localized magnetic flux rope is studied. A two-dimensional compressible MHD simulation code with a cylindrical symmetry is developed to study the wave modes associated with the evolution of flux ropes. It is found that in the initial phase both the fast magnetosonic wave and the Alfven wave are developed in the flux rope. After this initial phase, the Alfven wave becomes the dominant wave mode for the evolution of the magnetic flux rope and the radial expansion velocity of the flux rope is found to be negligible. Numerical results further show that even for a large initial azimuthal component of the magnetic field, the propagation velocity along the axial direction of the flux rope remains the Alfven velocity. It is also found that the localized magnetic flux rope tends to evolve into two separate magnetic ropes propagating in opposite directions. The simulation results are used to study the evolution of magnetic flux ropes associated with flux transfer events observed at the earth's dayside magnetopause and magnetic clouds in the interplanetary space.

Wei, C. Q.↗
NASA NTRS2000

Unsteady Aerodynamic Effects on the Flight Characteristics of an F-16XL Configuration

Unsteady aerodynamic models based on windtunnel forced oscillation test data and analyzed with a fuzzy logic algorithm arc incorporated into an F-16XL flight simulation code. The reduced frequency needed in the unsteady models is numerically calculated by using a limited prior time history of state variables in a least-square sense. Numerical examples arc presented to show the accuracy of the calculated reduced frequency. Oscillatory control inputs are employed to demonstrate the differences in the flight characteristics based on unsteady and quasi-steady aerodynamic models. Application of the unsteady aerodynamic models is also presented and the results are compared with one set of F16XIL longitudinal maneuver flight data. It is shown that the main differences in dynamic response are in the lateral-directional characteristics, with the quasi-steady model being more stable than the flight vehicle, while the unsteady model being more unstable. Similar conclusions can also be made in a simulated rapid sideslipping roll. To improve unsteady aerodynamic modeling, it is recommended to acquire test data with coupled motions in pitch, roll and yaw.

Wang, Zhongjun↗
NASA NTRS2016

A New Code SORD for Simulation of Polarized Light Scattering in the Earth Atmosphere

We report a new publicly available radiative transfer (RT) code for numerical simulation of polarized light scattering in plane-parallel atmosphere of the Earth. Using 44 benchmark tests, we prove high accuracy of the new RT code, SORD (Successive ORDers of scattering). We describe capabilities of SORD and show run time for each test on two different machines. At present, SORD is supposed to work as part of the Aerosol Robotic NETwork (AERONET) inversion algorithm. For natural integration with the AERONET software, SORD is coded in Fortran 90/95. The code is available by email request from the corresponding (first) author or from ftp://climate1.gsfc.nasa.gov/skorkin/SORD/.

AERONET↗