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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flow fields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane normal to the wing leading edge and compared to 2-D results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that the flow field over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3-D calculation.

Potapczuk, Mark G.↗
NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional ice accretion modeling method is initiated. This first step toward creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flowfields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is intended as a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3D method for a MS-317 swept wing geometry are projected onto a 2D plane normal to the wing leading edge and compared to 2D results for the same geometry. These results indicate that the flowfield over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3D calculation.

Potapczuk, M. G.↗
NASA NTRS1990

Component model reduction via the projection and assembly method

This paper explains the projection and assembly model reduction method which has been used to derive reduced order component models for the Galileo spacecraft. Assembly of reduced order component models produces a reduced order system model which can then be used in multibody simulation codes for efficient run times. The methodology is explained and a proof is given showing the exact reproduction of selected significant system modes. Frequency bounds are obtained for all modes produced when the reduced order components are assembled. The projection and assembly method is demonstrated on two examples. The first example is a simplified model of the Galileo spacecraft, while the second example is the model of the Galileo spacecraft in its early mission configuration. When articulation of components is allowed, the method may not reproduce all system modes precisely. Remedies for this problem are suggested.

Kissel, Glen J.↗
NASA NTRS1990

Collisional electrostatic ion cyclotron waves as a possible source of energetic heavy ions in the magnetosphere

A new mechanism is proposed for the source of energetic heavy ions (NO/+/, O2/+/, and O/+/) found in the magnetosphere. Simulations using a multispecies particle simulation code for resistive current-driven electrostatic ion cyclotron waves show transverse and parallel bulk heating of bottomside ionospheric heavy ion populations. The dominant mechanism for the transverse bulk heating is resonant ion heating by wave-particle ion trapping. Using a linear kinetic dispersion relation for a magnetized, collisional, homogenous, and multiion plasma, it is found that collisional electrostatic ion cyclotron waves near the NO(+), O2(+), and O(+) gyrofrequencies are unstable to field-aligned currents of 50 microA/sq m for a typical bottomside ionosphere.

Providakes, Jason↗
NASA NTRS1990

Magnetic asymmetries of unmagnetized planets

This letter discusses the results produced by three-dimensional hybrid particle code simulations of the solar wind interaction with unmagnetized planets such as Venus and Mars. The solar wind velocity is perpendicular to the IMF in the cases studied. It is found that there are asymmetries in both the magnetic structure and shock location for spherical obstacles ranging in radius from 1000 km to 6000 km. The asymmetries found are due to differences in the electron and ion current paths (diamagnetic behavior). Mass loading of 0(+) was not included in these simulations.

Brecht, Stephen H.↗
NASA NTRS1989

Simulations of light scattering in planetary rings

The current status of efforts to model light scattering in the Uranian and Saturnian rings is briefly reviewed, with an emphasis on the treatment of thin, vertically warped, or clumpy rings. Consideration is given to observations of an opposition surge in the rings (indicating that they are physically thick) and possible reasons for the failure of classical models to predict the forward-scattering characteristics of the B and inner A rings of Saturn, the quadrupole brightness variations in the A ring, the increase in reflectivity with optical depth in the A and B rings, the increase in brightness with solar and observer angle in the B ring, and systematic brightness variations in images of the inner Cassini division. It is pointed out that some of these discrepancies are eliminated by the use of a ray-tracing simulation code which determines the singly scattered light from a layer one or many particles thick.

Dones, L.↗
NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗
NASA NTRS1990

Theoretical and experimental studies of baroclinic processes

The geophysical fluid flow simulation code, GEOSIM, is being used to study the phenomenon of vacillation in the baroclinic annulus. Having verified that the code predicts vacillation for the same points as the experiments, the work is aiming toward explaining the mechanics of vacillation and pointing out some of the sensitivities of the results to the numerical method. Researchers are finding that there is a structural change associated with amplitude vacillation, where the structural changes are in the vertical. The results disagree with the premise of Lindzen et al, that the vacillation is due to constructive and destructive interference of neutral modes with different phase speeds. The researchers are continuing to study the Spacelab 3 Geophysical Fluid Flow Cell (GFFC) results with horizontal temperature gradients and heating from below. GEOSIM has been used to compute a wide range of cases, and these are being compared with the observations. The computations and observations compare well, and the model is being used to extend the results beyond cases studied in the experiments and to study the mechanics and predictability of the flows. The study of fully nonlinear baroclinic instability using the GFFC apparatus is proceeding with the numerical code. While the first instability that occurs is of planetary scale, secondary instabilities consisting of small-scale, penetrative convection occurs where cold fluid flows over a warm surface. The simultaneous modeling of the planetary scale and the convective scale is possible because of the nonhydrostatic formulation of the model. Some of these results have been animated on the Stardent computer, which shows the explosive nature of the small-scale convection.

Miller, Tim↗
NASA NTRS1989

Computer simulation of radiation damage in gallium arsenide

A version of the binary-collision simulation code MARLOWE was used to study the spatial characteristics of radiation damage in proton and electron irradiated gallium arsenide. Comparisons made with the experimental results proved to be encouraging.

Stith, John J.↗
NASA NTRS1990

Modelling the solution growth of TGS crystals in low gravity

The experimental growth of triglycine sulfate (TGS) crystals from aqueous solution is modeled here in two dimensions using the PHOENICS finite volume code. Simulations are carried out for steady, impulsive, and periodic accelerations in order to determine tolerable acceleration levels. Scaling arguments are used to estimate the times required for thermal and solutal variations from the initial equilibrium state to be diffusively transported throughout the system, and to obtain order of magnitude information on the relative magnitudes of diffusive and convective transport. The computed concentration fields reflect the features of the concentration distributions found experimentally during experiments conducted aboard Spacelab 3 in 1985.

Nadarajah, Arunan↗
NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict three-dimensional ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗
NASA NTRS1991

Reconnection voltage as a function of IMF clock angle

Magnetic reconnection between the IMF and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. Both analytic modeling and analysis of geophysical data have shown that this coupling process should be a sensitive function of the clock angle of the IMF. Results are presented from a three-dimensional, MHD, global numerical simulation code for the reconnection voltage between the closed geomagnetic field and the IMF as a function of the IMF clock angle. These results are consistent with a sin(theta/2) functional behavior.

Fedder, J. A.↗
NASA NTRS1991

Structure of medium Mach number quasi-parallel shocks - Upstream and downstream waves

The transition from steady low-Mach-number to unsteady high-Mach-number quasi-parallel shocks was investigated by performing large-scale 1D hybrid code simulations at increasing Mach numbers. It was found that only at very low Mach number shocks the steepening is limited by upstream phase-standing whistlers, as predicted by the classical theory (Tidman and Northrop, 1968). In the intermediate region of Mach numbers between 1.5 and 3.5, a very diverse behavior is observed. Backstreaming ions generate fast magnetosonic waves which dominate the upstream, with wavelengths longer than phase-standing whistlers. At increasing Mach numbers, the phase and group velocities of the dominant waves are reduced until they point back toward the shock; when there is sufficient energy flux in these waves, they lead to unsteady shock behavior and eventually to shock reformation.

Krauss-Varban, D.↗
NASA NTRS1992

Electrostatic Particle-In-Cell Code For Hypercube Computer

Code simulates two-dimensional motions of plasma particles in self-consistent electrostatic field and externally applied magnetic field developed for execution on Mark IIIfp hypercube computer. Based on generalization of one-dimensional particle-in-cell algorithm applicable to many different parallel computing architectures. Intermediate product of continuing effort to speed particle-in-cell computations by taking advantage of distributed-memory parallel computers like those of hypercube class.

Ferraro, Robert D.↗
NASA NTRS1991

Studies of the charging of a thin dust layer in a plasma

The present study is designed to extend Wilson's model to include refinements in the actual operational design of the simulation code itself and to include additional physics in order to assess the resulting electric field structure and grain charges. Most of the initial phases of the work have focussed on the fact that the code, though in principle reliable, suffers from severe time constraints, requiring as must as 10-50 hours of CPU time for typical runs on a VAX computer. It is likely that by staging the code into steps in which early steps are carried out crudely and later steps successively fine tuned the total CPU time can be cut significantly. Preliminary results indicate the possibility to cut CPU time to as much as 1/3 or 1/5 the former time. The original Wilson code has been carefully examined, tested and documented. Shortcuts and staging systems have been introduced and tested which should make the code produce more accurate results with a finer grid system than would have been feasible heretofor. Continued testing and implementation of these refinements will be carried out at the author's home institution in the near future. In addition to the effort at reducing CPU time to a manageable amount, work is in progress to include a dust-size distribution instead of assuming all dust grains have the same size.

Peterson, Lennart R.↗
NASA NTRS1992

Statistical discrete particle simulation studies

A particle simulation code was developed that is suited to Multiple Instruction Multiple Data multiprocessor computers. The resulting code is designed to be portable between a large class of multiprocessor computer architectures, including hypercubes (e.g., Intel iPSC/860), shared memory machines (e.g., Cray, SGI, DASH), and simple uniprocessors (e.g., SUN, VAX, IBM PC). This code provides a powerful general framework that can be adapted to specific applications, such as thermomechanical modeling, general 3-D geometry support, or visualization support via a codeveloped CPlot data reduction and visualization system. Chemistry models were implemented and validated with comparison to continuum solutions of thermochemically relaxing gas mixtures and with experimental results from the high speed flow about a circular cylinder. Furthermore, the exchange of energy between translational and internal modes was modeled. New models were developed, combining greater efficiency than earlier phenomenological models. They offer greater physical detail by addressing vibrational energy as a discretely distributed quantity.

Mcdonald, Jeffrey D.↗
NASA NTRS1991

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble averaged results for a widely used 1-D turbulence approximation.

J. Riaz↗