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 757 records · Page 42

The new challenge of computational aeroscience

This paper discusses NASA's Computational Aerosciences (CAS) Project of the High Performance Computing and Communications Program (HPCCP). The project is aimed at developing advanced, multidisciplinary simulation capabilities for aerospace vehicle and propulsion system design. It is also aimed at overcoming computational performance barriers by accelerating the development of parallel computer technology. The goals and approach of the CAS Project are described and the challenges to its implementation are addressed. Specific vehicle class simulations to be demonstrated and the principal mutidisciplinary modeling approaches to be emphasized are described. The computational speed and memory requirements for representative multidisciplinary applications are estimated. Finally, the state of parallel computer technology including programming issues and the results of performance measurements are explored.

Bailey, F. R.↗

Pitch-angle scattering of cometary ions - Computer simulations

The relationships between the rate of pitch-angle scattering induced by wave-particle interactions and two parameters, alpha, and the cometary ion interjection rate, Lambda-j, are studied with the aid of one-dimensional electromagnetic hybrid computer simulations of homogeneous plasmas. The simulations are done with relatively weak ion injection rates commensurate with those in the distant environment of the Comet Halley. In response to the linear temporal growth of the fluctuating magnetic field energy, the injected ions pitch-angle scatter toward isotropy in both the quasi-parallel and quasi-perpendicular regimes. For the injection of cometary oxygen ions, the simulations show pitch-angle scattering rates that increase as the square root of the injection rate. The oxygen ion pitch-angle scattering rates are roughly twice as large at perpendicular injection, again in qualitative agreement with observations.

Gary, S. P.↗

A parallel row-based algorithm for standard cell placement with integrated error control

A new row-based parallel algorithm for standard-cell placement targeted for execution on a hypercube multiprocessor is presented. Key features of this implementation include a dynamic simulated-annealing schedule, row-partitioning of the VLSI chip image, and two novel approaches to control error in parallel cell-placement algorithms: (1) Heuristic Cell-Coloring; (2) Adaptive Sequence Length Control.

Sargent, Jeff S.↗

Investigation of radiative interaction in laminar flows using Monte Carlo simulation

The Monte Carlo method (MCM) is employed to study the radiative interactions in fully developed laminar flow between two parallel plates. Taking advantage of the characteristics of easy mathematical treatment of the MCM, a general numerical procedure is developed for nongray radiative interaction. The nongray model is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. To validate the Monte Carlo simulation for nongray radiation problems, the results of radiative dissipation from the MCM are compared with two available solutions for a given temperature profile between two plates. After this validation, the MCM is employed to solve the present physical problem and results for the bulk temperature are compared with available solutions. In general, good agreement is noted and reasons for some discrepancies in certain ranges of parameters are explained.

Liu, Jiwen↗

Particle Acceleration and Emission in Relativistic Jets

Shock wave acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

Particle Acceleration and Radiation associated with Magnetic Field Generation from Relativistic Collisionless Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

Particle Acceleration in Relativistic Jets Due to Weibel Instability

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

Flare model chromospheres and photospheres

Homogeneous plane-parallel model atmospheres for solar flares have been constructed to approximately simulate observations of flares. The wings of the Ca II lines have been used to derive flare upper photosphere models, which indicate temperature increases of about 100 K over the temperature distribution in the pre-existing facula at a height of 300 km. In the case of flares covering sunspots the temperature rise seems to occur much higher in the atmosphere. We find that with increasing flare importance the heights of the upper chromosphere and transition region occur lower in the solar atmosphere, accounting for the factor of 60-600 increase in pressure in these regions relative to the quiet sun. The Ca II line profiles agree with observations only by assuming a macrovelocity distribution that increases with height. Also, the chromospheric parts of flares appear to be highly inhomogeneous. We show that shock and particle heated flare models do not agree with the observations and propose a thermal response model for flares. In particular, it appears that heating in the photosphere is an essential aspect of flares.

Machado, M. E.↗

Anomalous transport by magnetohydrodynamic Kelvin-Helmholtz instabilities in the solar wind-magnetosphere interaction

The high latitude, or downstream flank, and dayside low latitude magnetospheric boundaries are modeled in an MHD simulation of Kelvin-Helmholtz instablities in a compressible plasma for parallel and transverse configurations. Detailed attention is given to the nonlinear consequences of the instabilities for several different Alfven and sound Mach numbers in both configurations. Emphasis is given to the anomalous transport of momentum and energy by the Kelvin-Helmholtz instabilities across the magnetospheric boundary, which is crucial in evaluating the instabilities' contribution to magnetospheric convection. It is concluded that the anomalous tangential stress at the magnetospheric boundary caused by the instability may be of the order of 1 percent of the magnetosheath flow momentum adjacent to the boundary, and gives a contribution to the convection potential drop over the polar cap of 10-30 kV.

Miura, A.↗

Solar reflection from interacting and shadowing cloud elements

An attempt is made to relate the reflection from a broken cloud field to that from a plane parallel cloud of the same optical properties by using Monte Carlo simulations for a very simple model of an array of regularly spaced, infinitely long bar clouds. Reflection from the array is represented by an effective cloud fraction, and it is shown that for many cases this may be obtained by treating the cloud elements to be noninteracting azimuthally symmetric reflectors. The effective cloud fraction is a function of solar zenith angle, the aspect ratio of the individual cloud elements, and the normal or zenith cloud fraction of the array. Model computations for reflection from an array of randomly distributed identical cylinders are also presented as an example of a possible application of this technique.

HARSHVARDHAN↗

Zero-Lift Wave Drag of Complex Aircraft Configurations

WAVDRAG calculates supersonic zero-lift wave drag of complex aircraft configurations. Numerical model of aircraft used throughout design process from concept to manufacturing. Incorporates extended geometric input capabilities to permit use of more accurate mathematical model. Engineer defines aircraft components as fusiform or nonfusiform in terms of traditional parallel contours or nonintersecting contours in any direction. Laterally asymmetric configurations simulated. WAVDRAG calculates total drag and wave-drag coefficient of specified aircraft configuration.

Craidon, C., B.↗

Anomalous transport by Kelvin-Helmholtz instabilities

A simulation of magnetohydrodynamic Kelvin-Helmholtz instabilities has been performed for parallel and transverse configurations, modeling high latitude (or downstream flanks) and dayside low latitude magnetospheric boundaries. In the parallel configuration, a super-Alfvenic and transsonic shear flow develops into small eddies, which strongly compresses, twists, and hence amplifies the magnetic field by the dynamo action with an amplification factor of M(A)/2. In the nonlinear stage, however large the initial Alfven mach number M(A) may be, the magnetic field amplified and twisted by the hydromagnetic flow vortices reacts back upon the flow evolution, and the flow vortices cascade into smaller structures. In the transverse configuration, the instability leads to the formation of a fast shock discontinuity from an initially subfast shear flow. Anomalous tangential stress by the instability in the transverse configuration reaches 1 percent of the magnetosheath momentum flux, but for the parallel configuration, the anomalous transport is 2-3 times larger than the anomalous transport in the transverse configuration. The anomalous transport for both configurations satisfies the requirement of the viscous-like interaction at the magnetospheric boundary.

Miura, A.↗

Electron acceleration in stochastic double layers

Transversely localized double layers evolve randomly in turbulent regions of strongly magnetized plasma carrying current along the magnetic field. Results from numerical simulations and spacecraft observations in the auroral plasma indicate that the parallel electric field in such regions is microscopically intermittent or stochastic. The implications of stochastic double layer fields on electron acceleration will be discussed in terms of a statistical process involving ensemble averages over test particle motion. A Fokker-Planck equation can be derived for the electron phase space density, which depends on the mean and rms amplitudes of the double layers, the mean double layer density, and the initial electron velocity distribution. It is shown that the resulting electron acceleration is very sensitive to the ratio of the initial electron energy to the rms double layer amplitude. When this ratio is large, the acceleration process differs little from that expected in a dc electric field. When it is small, stochastic heating competes with directed acceleration. Evidence for both cases can be found in the auroral ionosphere in association with so-called inverted-V precipitation and collimated edge precipitation.

Lotko, William↗

Synchronization trigger control system for flow visualization

The use of cinematography or holographic interferometry for dynamic flow visualization in an internal combustion engine requires a control device that globally synchronizes camera and light source timing at a predefined shaft encoder angle. The device is capable of 0.35 deg resolution for rotational speeds of up to 73 240 rpm. This was achieved by implementing the shaft encoder signal addressed look-up table (LUT) and appropriate latches. The developed digital signal processing technique achieves 25 nsec of high speed triggering angle detection by using direct parallel bit comparison of the shaft encoder digital code with a simulated angle reference code, instead of using angle value comparison which involves more complicated computation steps. In order to establish synchronization to an AC reference signal whose magnitude is variant with the rotating speed, a dynamic peak followup synchronization technique has been devised. This method scrutinizes the reference signal and provides the right timing within 40 nsec. Two application examples are described.

Chun, K. S.↗

Neural nets on the MPP

The Massively Parallel Processor (MPP) is an ideal machine for computer experiments with simulated neural nets as well as more general cellular automata. Experiments using the MPP with a formal model neural network are described. The results on problem mapping and computational efficiency apply equally well to the neural nets of Hopfield, Hinton et al., and Geman and Geman.

Hastings, Harold M.↗

Real time animation of space plasma phenomena

In pursuit of real time animation of computer simulated space plasma phenomena, the code was rewritten for the Massively Parallel Processor (MPP). The program creates a dynamic representation of the global bowshock which is based on actual spacecraft data and designed for three dimensional graphic output. This output consists of time slice sequences which make up the frames of the animation. With the MPP, 16384, 512 or 4 frames can be calculated simultaneously depending upon which characteristic is being computed. The run time was greatly reduced which promotes the rapid sequence of images and makes real time animation a foreseeable goal. The addition of more complex phenomenology in the constructed computer images is now possible and work proceeds to generate these images.

Jordan, K. F.↗

Oxygen acoustic solitary waves in a magnetized plasma

Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.

Qian, S.↗