Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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 811 records · Page 45

Electrostatic waves due to field-aligned electron beams in the low-latitude boundary layer

Mass-resolved ion, electron, and plasma wave data obtained from several low-latitude boundary layer (LLBL) crossings by the AMPTE CCE satellite are analyzed. The data clearly separate the LLBL from the adjacent magnetosheath and magnetosphere. Attention was focused on wave-particle interactions involving electrons. Electron beams were found to be present in the LLBL during the southward interplanetary magnetic field, along with a simultaneous enhancement of electrostatic waves with parallel polarization. Linear theory analysis shows that for plasma conditions in the LLBL, electron beams are unstable to electrostatic waves that propagate parallel to the local magnetic field, in agreement with observations. A numerical simulation study of the beam-plasma interaction in the LLBL shows that the instability saturates by thermalization of the beam but that a beamlike structure can still remain in the electron distribution for certain initial parameters. It is suggested that peaks in the electron velocity distribution function may be found in the LLBL away from the beam source region.

Peroomian, V.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis

The purpose of this study is the development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates and the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The procedure was also modified to allow coarse parallelization of the solution algorithm. This document is a final report outlining the development and techniques used in the procedure. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Numerical dissipation is used to gain solution stability but is reduced in viscous dominated flow regions. Local time stepping and implicit residual smoothing are used to increase the rate of convergence. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes being generated by the system (TIGG3D) developed earlier under this contract. The grid generation scheme meets the average-passage requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. Pure internal flow solutions were obtained as well as solutions with flow about the cowl/nacelle and various engine core flow conditions. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗

The soliton transform and a possible application to nonlinear Alfven waves in space

The inverse scattering transform (IST) based on the derivative nonlinear Schroedinger (DNLS) equation is applied to a complex time series of nonlinear Alfven wave data generated by numerical simulation. The IST describes the long-time evolution of quasi-parallel Alfven waves more efficiently than the Fourier transform, which is adapted to linear rather than nonlinear problems. When dissipation is added, so the conditions for the validity of the DNLS are not strictly satisfied, the IST continues to provide a compact description of the wavefield in terms of a small number of decaying envelope solitons.

Hada, T.↗

Flight- and ground-test correlation study of BMDO SDS materials: Phase 1 report

The NASA Evaluation of Oxygen Interactions with Materials-3 (EOIM-3) experiment served as a test bed for a variety of materials that are candidates for Ballistic Missile Defense Organization (BMDO) space assets. The materials evaluated on this flight experiment were provided by BMDO contractors and technology laboratories. A parallel ground exposure evaluation was conducted using the FAST atomic-oxygen simulation facility at Physical Sciences, Inc. The EOIM-3 materials were exposed to an atomic oxygen fluence of approximately 2.3 x 10(exp 2) atoms/sq. cm. The ground-exposed materials' fluence of 2.0 - 2.5 x 10(exp 2) atoms/sq. cm permits direct comparison of ground-exposed materials' performance with that of the flight-exposed specimens. The results from the flight test conducted aboard STS-46 and the correlative ground exposure are presented in this publication.

Chung, Shirley Y.↗

Ignition Transient Calculations in the Space Shuttle Solid Rocket Motor

The work presented is part of an effort to develop a multidimensional ignition transient model for large solid propellant rocket motors. On the Space Shuttle, the ignition transient in the slot is induced when the igniter, itself a small rocket motor, is fired into the head-end portion of the main rocket motor. The computational results presented in this paper consider two different igniter configurations. The first configuration is a simulated Space Shuttle RSRM igniter which has one central nozzle that is parallel to the centerline of the motor. The second igniter configuration has a nozzle which is canted at an angle of 45 deg from the centerline of the motor. This paper presents a computational fluid dynamic (CFD) analyses of certain flow field characteristics inside the solid propellant star grain slot of the Space Shuttle during the ignition transient period of operation for each igniter configuration. The majority of studies made to date regarding ignition transient performance in solid rocket motors have concluded that the key parameter to be determined is the heat transfer rate to the propellant surface and hence the heat transfer coefficient between the gas and the propellant. In this paper the heat transfer coefficients, pressure and velocity distributions are calculated in the star slot. In order to validate the computational method and to attempt to establish a correlation between the flow field characteristics and the heat transfer rates a series of cold flow experimental investigations were conducted. The results of these experiments show excellent qualitative and quantitative agreement with the pressure and velocity distributions obtained from the CFD analysis. The CFD analysis utilized a classical pipe flow type correlation for the heat transfer rates. The experimental results provide an excellent qualitative comparison with regard to spatial distribution of the heat transfer rates as a function of nozzle configuration and igniter pressure. The results indicate that from a quantitative point of view that the pipe flow correlation gives reasonably good results. Furthermore, there appears to be a direct correlation between igniter pressure and an average Reynolds number in the star grain slot. This may lead to a simple method for modifying the convection heat transfer correlation. Calculated results of pressure-vs-time for the first 200 msec of motor firing of the Space Shuttle RSRM support the trends shown for the heat transfer rate comparisons between the cold flow CFD and experimental data.

Jenkins, Rhonald M.↗

FFT Computation with Systolic Arrays, A New Architecture

The use of the Cooley-Tukey algorithm for computing the l-d FFT lends itself to a particular matrix factorization which suggests direct implementation by linearly-connected systolic arrays. Here we present a new systolic architecture that embodies this algorithm. This implementation requires a smaller number of processors and a smaller number of memory cells than other recent implementations, as well as having all the advantages of systolic arrays. For the implementation of the decimation-in-frequency case, word-serial data input allows continuous real-time operation without the need of a serial-to-parallel conversion device. No control or data stream switching is necessary. Computer simulation of this architecture was done in the context of a 1024 point DFT with a fixed point processor, and CMOS processor implementation has started.

Boriakoff, Valentin↗

An Experimental and Computational Investigation of Oscillating Airfoil Unsteady Aerodynamics at Large Mean Incidence

A major challenge in the design and development of turbomachine airfoils for gas turbine engines is high cycle fatigue failures due to flutter and aerodynamically induced forced vibrations. In order to predict the aeroelastic response of gas turbine airfoils early in the design phase, accurate unsteady aerodynamic models are required. However, accurate predictions of flutter and forced vibration stress at all operating conditions have remained elusive. The overall objectives of this research program are to develop a transition model suitable for unsteady separated flow and quantify the effects of transition on airfoil steady and unsteady aerodynamics for attached and separated flow using this model. Furthermore, the capability of current state-of-the-art unsteady aerodynamic models to predict the oscillating airfoil response of compressor airfoils over a range of realistic reduced frequencies, Mach numbers, and loading levels will be evaluated through correlation with benchmark data. This comprehensive evaluation will assess the assumptions used in unsteady aerodynamic models. The results of this evaluation can be used to direct improvement of current models and the development of future models. The transition modeling effort will also make strides in improving predictions of steady flow performance of fan and compressor blades at off-design conditions. This report summarizes the progress and results obtained in the first year of this program. These include: installation and verification of the operation of the parallel version of TURBO; the grid generation and initiation of steady flow simulations of the NASA/Pratt&Whitney airfoil at a Mach number of 0.5 and chordal incidence angles of 0 and 10 deg.; and the investigation of the prediction of laminar separation bubbles on a NACA 0012 airfoil.

Capece, Vincent R.↗

SALT: The Simulator for the Analysis of LWP Timing

With the emergence of new processor architectures that are highly multithreaded, and support features such as full/empty memory semantics and split-phase memory transactions, the need for a processor simulator to handle these features becomes apparent. This paper describes such a simulator, called SALT.

simulation↗

PPM Receiver Implemented in Software

A computer program has been written as a tool for developing optical pulse-position- modulation (PPM) receivers in which photodetector outputs are fed to analog-to-digital converters (ADCs) and all subsequent signal processing is performed digitally. The program can be used, for example, to simulate an all-digital version of the PPM receiver described in Parallel Processing of Broad-Band PPM Signals (NPO-40711), which appears elsewhere in this issue of NASA Tech Briefs. The program can also be translated into a design for digital PPM receiver hardware. The most notable innovation embodied in the software and the underlying PPM-reception concept is a digital processing subsystem that performs synchronization of PPM time slots, even though the digital processing is, itself, asynchronous in the sense that no attempt is made to synchronize it with the incoming optical signal a priori and there is no feedback to analog signal processing subsystems or ADCs. Functions performed by the software receiver include time-slot synchronization, symbol synchronization, coding preprocessing, and diagnostic functions. The program is written in the MATLAB and Simulink software system. The software receiver is highly parameterized and, hence, programmable: for example, slot- and symbol-synchronization filters have programmable bandwidths.

Gray, Andrew↗

Characterizing a Wake-Free Safe Zone for the Simplified Aircraft-Based Paired Approach Concept

The Federal Aviation Administration (FAA) has proposed a concept of operations geared towards achieving increased arrival throughput at U.S. Airports, known as the Simplified Aircraft-based Paired Approach (SAPA) concept. In this study, a preliminary characterization of a wake-free safe zone (WFSZ) for the SAPA concept has been performed. The experiment employed Monte-Carlo simulations of varying approach profiles by aircraft pairs to closely-spaced parallel runways. Three different runway lateral spacings were investigated (750 ft, 1000 ft and 1400 ft), along with no stagger and 1500 ft stagger between runway thresholds. The paired aircraft were flown in a leader/trailer configuration with potential wake encounters detected using a wake detection surface translating with the trailing aircraft. The WFSZ is characterized in terms of the smallest observed initial in-trail distance leading to a wake encounter anywhere along the approach path of the aircraft. The results suggest that the WFSZ can be characterized in terms of two primary altitude regions, in ground-effect (IGE) and out of ground-effect (OGE), with the IGE region being the limiting case with a significantly smaller WFSZ. Runway stagger was observed to only modestly reduce the WFSZ size, predominantly in the OGE region.

Guerreiro, Nelson M.↗

Pair Cascades and Deathlines in Offset Magnetic Dipole Fields

We investigate electron-positron pair cascades in a dipole magnetic field whose axis is offset from the neutron star center. In such a field geometry, the polar cap is displaced from the neutron star symmetry axis and the field line radius of curvature is modified. Using the modified parallel electric field near the polar cap of an offset dipole, we simulate pair cascades to determine the pair deathlines and pair multiplicities as a function of the offset parameter. We find that the pair multiplicity can change dramatically with a modest offset, with a significant increase on one side of the polar cap. Lower pair deathlines allow a larger fraction of the pulsar population, that include old and millisecond pulsars, to produce cascades with high multiplicity.

Harding, Alice↗

Analysis and Simulation of the Simplified Aircraft-Based Paired Approach Concept With the ALAS Alerting Algorithm in Conjunction With Echelon and Offset Strategies

This report presents analytical and simulation results of an investigation into proposed operational concepts for closely spaced parallel runways, including the Simplified Aircraft-based Paired Approach (SAPA) with alerting and an escape maneuver, MITRE?s echelon spacing and no escape maneuver, and a hybrid concept aimed at lowering the visibility minima. We found that the SAPA procedure can be used at 950 ft separations or higher with next-generation avionics and that 1150 ft separations or higher is feasible with current-rule compliant ADS-B OUT. An additional 50 ft reduction in runway separation for the SAPA procedure is possible if different glideslopes are used. For the echelon concept we determined that current generation aircraft cannot conduct paired approaches on parallel paths using echelon spacing on runways less than 1400 ft apart and next-generation aircraft will not be able to conduct paired approach on runways less than 1050 ft apart. The hybrid concept added alerting and an escape maneuver starting 1 NM from the threshold when flying the echelon concept. This combination was found to be effective, but the probability of a collision can be seriously impacted if the turn component of the escape maneuver has to be disengaged near the ground (e.g. 300 ft or below) due to airport buildings and surrounding terrain. We also found that stabilizing the approach path in the straight-in segment was only possible if the merge point was at least 1.5 to 2 NM from the threshold unless the total system error can be sufficiently constrained on the offset path and final turn.

Torres-Pomales, Wilfredo↗

Near-Body Grid Adaption for Overset Grids

A solution adaption capability for curvilinear near-body grids has been implemented in the OVERFLOW overset grid computational fluid dynamics code. The approach follows closely that used for the Cartesian off-body grids, but inserts refined grids in the computational space of original near-body grids. Refined curvilinear grids are generated using parametric cubic interpolation, with one-sided biasing based on curvature and stretching ratio of the original grid. Sensor functions, grid marking, and solution interpolation tasks are implemented in the same fashion as for off-body grids. A goal-oriented procedure, based on largest error first, is included for controlling growth rate and maximum size of the adapted grid system. The adaption process is almost entirely parallelized using MPI, resulting in a capability suitable for viscous, moving body simulations. Two- and three-dimensional examples are presented.

Buning, Pieter G.↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

Simulation and control of a 20 kHz spacecraft power system

A detailed computer representation of four Mapham inverters connected in a series, parallel arrangement has been implemented. System performance is illustrated by computer traces for the four Mapham inverters connected to a Litz cable with parallel resistance and dc receiver loads at the receiving end of the transmission cable. Methods of voltage control and load sharing between the inverters are demonstrated. Also, the detailed computer representation is used to design and to demonstrate the advantages of a feed-forward voltage control strategy. It is illustrated that with a computer simulation of this type, the performance and control of spacecraft power systems may be investigated with relative ease and facility.

Wasynczuk, O.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

Temporal Precedence Checking for Switched Models and its Application to a Parallel Landing Protocol

This paper presents an algorithm for checking temporal precedence properties of nonlinear switched systems. This class of properties subsume bounded safety and capture requirements about visiting a sequence of predicates within given time intervals. The algorithm handles nonlinear predicates that arise from dynamics-based predictions used in alerting protocols for state-of-the-art transportation systems. It is sound and complete for nonlinear switch systems that robustly satisfy the given property. The algorithm is implemented in the Compare Execute Check Engine (C2E2) using validated simulations. As a case study, a simplified model of an alerting system for closely spaced parallel runways is considered. The proposed approach is applied to this model to check safety properties of the alerting logic for different operating conditions such as initial velocities, bank angles, aircraft longitudinal separation, and runway separation.

Duggirala, Parasara Sridhar↗