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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Hot Chips and Hot Interconnects for High End Computing Systems

I will discuss several processors: 1. The Cray proprietary processor used in the Cray X1; 2. The IBM Power 3 and Power 4 used in an IBM SP 3 and IBM SP 4 systems; 3. The Intel Itanium and Xeon, used in the SGI Altix systems and clusters respectively; 4. IBM System-on-a-Chip used in IBM BlueGene/L; 5. HP Alpha EV68 processor used in DOE ASCI Q cluster; 6. SPARC64 V processor, which is used in the Fujitsu PRIMEPOWER HPC2500; 7. An NEC proprietary processor, which is used in NEC SX-6/7; 8. Power 4+ processor, which is used in Hitachi SR11000; 9. NEC proprietary processor, which is used in Earth Simulator. The IBM POWER5 and Red Storm Computing Systems will also be discussed. The architectures of these processors will first be presented, followed by interconnection networks and a description of high-end computer systems based on these processors and networks. The performance of various hardware/programming model combinations will then be compared, based on latest NAS Parallel Benchmark results (MPI, OpenMP/HPF and hybrid (MPI + OpenMP). The tutorial will conclude with a discussion of general trends in the field of high performance computing, (quantum computing, DNA computing, cellular engineering, and neural networks).

Saini, Subhash↗
NASA NTRS2016

Electron Scale Structures and Magnetic Reconnection Signatures in the Turbulent Magnetosheath

Collisionless space plasma turbulence can generate reconnecting thin current sheets as suggested by recent results of numerical magnetohydrodynamic simulations. The Magnetospheric Multiscale (MMS) mission provides the first serious opportunity to verify whether small ion-electron-scale reconnection, generated by turbulence, resembles the reconnection events frequently observed in the magnetotail or at the magnetopause. Here we investigate field and particle observations obtained by the MMS fleet in the turbulent terrestrial magnetosheath behind quasi-parallel bow shock geometry. We observe multiple small-scale current sheets during the event and present a detailed look of one of the detected structures. The emergence of thin current sheets can lead to electron scale structures. Within these structures, we see signatures of ion demagnetization, electron jets, electron heating, and agyrotropy suggesting that MMS spacecraft observe reconnection at these scales.

Yordanova, E.↗
NASA NTRS2017

Structure and Properties of the Foreshock at Venus

The interaction of the solar wind with Venus is dominated by the planet's ionosphere that acts as an obstacle to the flow resulting in an induced magnetosphere and bow shock much smaller than their terrestrial counterparts. This study presents a 3-D electromagnetic hybrid (kinetic ions and fluid electrons) simulation of the solar wind interaction with an unmagnetized obstacle to examine the structure and properties of the Cytherean foreshock during periods of near radial IMF, that is, when it lies upstream of the ionosphere. The interaction between the backstreaming ions and the solar wind results in the generation of two classes of ULF waves: (1) parallel propagating sinusoidal waves with periods approx. 20-30 s and (2) highly oblique fast magnetosonic waves. The joint nonlinear evolution of these waves results in the formation of structures called foreshock cavitons with dimensions comparable to the size of the planet. Foreshock cavitons are also present in the terrestrial foreshock. The excavation of plasma and magnetic field from their cores leads to lower average densities and magnetic field strengths in the foreshock. As in the case of Earth, this excavation results in the formation of a fast magnetosonic pulse/shock at the edge of the foreshock named the foreshock compressional boundary. Also similar to Earth is the formation of spontaneous hot flow anomalies (SHFAs) as foreshock cavitons approach the bow shock. The size and properties of SHFAs at Venus are comparable to those at Earth, and their existence has recently been established at Mars and Venus in a companion paper.

Omidi, N.↗
NASA NTRS2021

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of thePatch-workcode [2], for which several algorithms presented here were co-developed. Its purpose is to create a“multipatch” scheme compatible with numerical simulations of arbitrary equations of motion at any dis-cretization order in space and time. In thePatchworkframework, the global simulation is comprised of anarbitrary number of moving, local meshes, or “patches”, which are free to employ their own resolution, co-ordinate system/topology, physics equations, reference frame, and in our new approach, numerical method.Each local patch exchanges boundary data with a single global patch on which all other patches residethrough a client-router-server parallelization model. In generalizingPatchworkto be compatible witharbitrary order time integration,PatchworkMHDandPatchworkWavehave significantly improved theinterpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in theoriginalPatchworkcode. Furthermore, we extendPatchworkto bemultimethodby allowing multiplestate vectors to be updated simultaneously, with each state vector providing its own interpatch interpolationand transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partialdifferential equations. We demonstrate our changes through the implementation of a scalar wave toy-modelthat is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

Dennis B Bowen↗
NASA NTRS2022

Whistler Waves Generated by Nongyrotropic and Gyrotropic Electron Beams During Asymmetric Guide Field Reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

magnetic reconnection↗
NASA NTRSDate not supplied

Transients and Dayside Magnetospheric Interactions

Observations at and deep within the Earth’s dayside magnetosphere provide considerable evidence indicating that the solar wind-magnetosphere interaction is often unsteady. Bursty reconnection resulting in the formation of flux transfer events and transient auroral brightenings is common. Fluctuations in solar wind parameters may trigger some bursts of reconnection, while other bursts may result from intrinsic magnetopause instabilities. Large amplitude solitons often travel along the dayside magnetopause. Deeper within the dayside magnetosphere, quasi-periodic transient compressions energize radiation belt electrons and ring current ions. Most of these perturbations occur behind the quasi-parallel bow shock, indicating that they are driven by density structures generated by kinetic processes within the foreshock. This presentation compares and contrasts the amplitudes and occurrence patterns of density and dynamic pressure perturbations intrinsic to the solar wind with those generated within the foreshock, specifies how they move across the magnetosphere, examines recent global MHD simulation results for their interaction with the bow shock, magnetosheath, and magnetopause. It then outlines new paths toward quantifying the significance of bursty reconnection to the overall interaction, via new imaging and analysis techniques.

D. G. Sibeck↗
NASA NTRS1993

The Generation of Oblique Magnetosonic Waves

One of the outstanding issues regarding the excitation of magnetosonic waves has been the observational evidence that obliquely propagating waves are dominant at comets and planetary foreshocks despite the predictions of linear theory that maximum growth occurs at parallel propagation. To address this issue, we have conducted a detailed linear theory using a beam-ring distribution function. The results have shown that such distribution functions are associated with four separate instabilities. Two of these instabilities are similar to the right hand resonant ion/ion and the non-resonant instabilities which are also present in the case of a field aligned beam. The other two instabilities are associated with the presence of ring part of the distribution function. One of these instabilities excites magnetosonic waves with maximum growth in the oblique directions. The other excites Alfven waves with maximum growth in the oblique directions. The importance of the former instability is that it may explain the oblique nature of the magnetosonic waves observed at planetary foreshocks and comets. In order to understand the nonlinear properties of these various instabilities, we have also conducted 2-D hybrid (particle ions, fluid electrons) simulations. We have found that when the beam density is sufficiently large so that the non-resonant instability has the largest growth rate, these waves dominate the initial wave power in the system. At later times, however, the obliquely propagating magnetosonic waves become dominant. Another important finding was that when the two instabilities which excite magnetosonic waves have the largest growth rates, the system is dominated by the obliquely propagating magnetosonic waves. This is despite the fact that the largest growth rate for one of the instabilities occurs in the parallel direction. The exact cause of this is not currently understood and is under investigation.

Source record↗
NASA NTRS1984

NASA welding assessment program

A long duration test has been conducted for comparing various methods of attaching electrical interconnects to solar cells for near Earth orbit spacecraft. Representative solar array modules have been thermally cycled for 36,000 cycles between -80 and +80 C on this JPL and NASA Lewis Research Center sponsored work. This test simulates the environmental stress of more than 6 years on a near Earth spacecraft as it cycles in and out of the Earth's shadow. Evaluations of the integrity of these modules were made by visual and by electrical examinations before starting the cycling and then at periodic intervals during the cycling tests. Modules included examples of parallel gap and of ultrasonic welding, as well as soldering. The materials and fabrication processes are state of the art, suitable for forming large solar arrays of spacecraft quality. The modules survived his extensive cycling without detectable degradation in their ability to generate power under sunlight illumination.

Stofel, E. J.↗
NASA NTRS1990

Particle injection and acceleration at earth's bow shock - Comparison of upstream and downstream events

The injection and acceleration of thermal solar wind ions at the quasi-parallel earth's bow shock during radial interplanetary magnetic field conditions is investigated. Active Magnetospheric Particle Tracer Explorers/Ion Release Module satellite observations of complete proton spectra, and of heavy ion spectra above 10 keV/Q, made on September 12, 1984 near the nose of the shock, are presented and compared to the predictions of a Monte Carlo shock simulation which includes diffusive shock acceleration. It is found that the spectral observations are in good agreement with the predictions of the simulation when it is assumed that all accelerated ions originate in the solar wind and are injected into the acceleration mechanism by thermal leakage from the downstream plasma. The efficiency, which is determined directly from the downstream observations, is high, with at least 15 percent of the solar wind energy flux going into accelerated particles. The comparisons allow constraints to be placed on the rigidity dependence of the scattering mean free path and suggest that the upstream solar wind must be slowed substantially by backstreaming accelerated ions prior to undergoing a sharp transition in the viscous subshock.

Ellison, Donald C.↗
NASA NTRS1994

High Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

In order to predict the dynamic response of a flexible structure in a fluid flow, the equations of motion of the structure and the fluid must be solved simultaneously. In this paper, we present several partitioned procedures for time-integrating this focus coupled problem and discuss their merits in terms of accuracy, stability, heterogeneous computing, I/O transfers, subcycling, and parallel processing. All theoretical results are derived for a one-dimensional piston model problem with a compressible flow, because the complete three-dimensional aeroelastic problem is difficult to analyze mathematically. However, the insight gained from the analysis of the coupled piston problem and the conclusions drawn from its numerical investigation are confirmed with the numerical simulation of the two-dimensional transient aeroelastic response of a flexible panel in a transonic nonlinear Euler flow regime.

Felippa, C. A.↗
NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗
NASA NTRS1994

Shuttle Upper Atmosphere Mass Spectrometer Experimental Flight Results

Calibrated pressure measurements for species with mass-to-charge ratios up to 50 amu/e(-) were obtained trom the shuttle upper atmosphere mass spectrometer experiment during re-entry on the STS-35 mission. The principal experimental objective is to obtain measurements of freestream density in the hypersonic rarefied flow flight regime. Data were collected from 180 to about 87 km. However, data above 115 km were contaminated from a source of gas emanating from pressure transdueers connected in parallel to the mass spectrometer. At lower altitudes, the pressure transducer data are compared to the mass spectrometer total pressure with excellent agreement. Near the orifice entrance, a significant amount of CO2 was generated from chemical reactions. The freestream density in the rarefied flow flight regime is calculated using an orifice pressure coefficient model based upon direct simulation Monte Carlo results. This density, when compared with the 1976 U.S. Standard Atmosphere model, exhibits the wavelike nature seen on previous flights using accelerometry. Selected spectra are presented at higher altitudes (320 km) showing the effects of the ingestion of gases from a forward fuselage fuel dump.

Blanchard, R. C.↗
NASA NTRS2002

Load Balancing Strategies for Multi-Block Overset Grid Applications

The multi-block overset grid method is a powerful technique for high-fidelity computational fluid dynamics (CFD) simulations about complex aerospace configurations. The solution process uses a grid system that discretizes the problem domain by using separately generated but overlapping structured grids that periodically update and exchange boundary information through interpolation. For efficient high performance computations of large-scale realistic applications using this methodology, the individual grids must be properly partitioned among the parallel processors. Overall performance, therefore, largely depends on the quality of load balancing. In this paper, we present three different load balancing strategies far overset grids and analyze their effects on the parallel efficiency of a Navier-Stokes CFD application running on an SGI Origin2000 machine.

Djomehri, M. Jahed↗
NASA NTRS2009

A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners

A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state-of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they should behave as resonant liners while the construction of others suggest that they should behave as broadband attenuators. In each case the finite element theory is observed to predict the proper attenuation trend.

Watson, W. R.↗
NASA NTRS2011

EOS MLS Level 2 Data Processing Software Version 3

This software accepts the EOS MLS calibrated measurements of microwave radiances products and operational meteorological data, and produces a set of estimates of atmospheric temperature and composition. This version has been designed to be as flexible as possible. The software is controlled by a Level 2 Configuration File that controls all aspects of the software: defining the contents of state and measurement vectors, defining the configurations of the various forward models available, reading appropriate a priori spectroscopic and calibration data, performing retrievals, post-processing results, computing diagnostics, and outputting results in appropriate files. In production mode, the software operates in a parallel form, with one instance of the program acting as a master, coordinating the work of multiple slave instances on a cluster of computers, each computing the results for individual chunks of data. In addition, to do conventional retrieval calculations and producing geophysical products, the Level 2 Configuration File can instruct the software to produce files of simulated radiances based on a state vector formed from a set of geophysical product files taken as input. Combining both the retrieval and simulation tasks in a single piece of software makes it far easier to ensure that identical forward model algorithms and parameters are used in both tasks. This also dramatically reduces the complexity of the code maintenance effort.

Livesey, Nathaniel J.↗
NASA NTRS2014

Exploration of the Theoretical Physical Capacity of the John F. Kennedy International Airport Runway System

A design study was completed to explore the theoretical physical capacity (TPC) of the John F. Kennedy International Airport (KJFK) runway system for a northflow configuration assuming impedance-free (to throughput) air traffic control functionality. Individual runways were modeled using an agent-based, airspace simulation tool, the Airspace Concept Evaluation System (ACES), with all runways conducting both departures and arrivals on a first-come first-served (FCFS) scheduling basis. A realistic future flight schedule was expanded to 3.5 times the traffic level of a selected baseline day, September 26, 2006, to provide a steady overdemand state for KJFK runways. Rules constraining departure and arrival operations were defined to reflect physical limits beyond which safe operations could no longer be assumed. Safety buffers to account for all sources of operational variability were not included in the TPC estimate. Visual approaches were assumed for all arrivals to minimize inter-arrival spacing. Parallel runway operations were assumed to be independent based on lateral spacing distances. Resulting time intervals between successive airport operations were primarily constrained by same-runway and then by intersecting-runway spacing requirements. The resulting physical runway capacity approximates a theoretical limit that cannot be exceeded without modifying runway interaction assumptions. Comparison with current KJFK operational limits for a north-flow runway configuration indicates a substantial throughput gap of approximately 48%. This gap may be further analyzed to determine which part may be feasibly bridged through the deployment of advanced systems and procedures, and which part cannot, because it is either impossible or not cost-effective to control. Advanced systems for bridging the throughput gap may be conceptualized and simulated using this same experimental setup to estimate the level of gap closure achieved.

Neitzke, Kurt W.↗
NASA NTRS1984

Analytical description of an imperfect Fabry-Perot etalon

A novel method is presented for the investigation of Fabry-Perot etalon defects, incorporating the effects of microsurface imperfection, bowing, and departure from parallelism, as well as the effect of finite aperture on the etalon, into the perturbation function. The tolerance imperfections of reflective layers can be incorporated into the transfer coefficients of a noncorrugated thin-film system, and any defect can be investigated either separately or in combination with others through the use of the appropriate analytical or numerical simulation of the corrugated function. In the context of numerical computation, the present method permits the simulation and evaluation of different specific situations. The expressions derived for the transfer coefficients of a corrugated, equidistant thin film system may be applicable to other branches of applied optics.

Prikryl, I.↗
NASA NTRS1995

CDMA with interference cancellation for multiprobe missions

Code division multiple-access spread spectrum has been proposed for use in future multiprobe/multispacecraft missions. This article considers a general parallel interference-cancellation scheme that significantly reduces the degradation effect of probe (user) interference but with a lesser implementation complexity than the maximum-likelihood technique. The scheme operates on the fact that parallel processing simultaneously removes from each probe (user) the total interference produced by the remaining most reliably received probes (users) accessing the channel. The parallel processing can be done in multiple stages. The proposed scheme uses tentative decision devices with different optimum thresholds at the multiple stages to produce the most reliably received data for generation and cancellation of probe/spacecraft interference. The one-stage interference cancellation was analyzed for two types of tentative decision devices, namely, hard and null zone decisions. Simulation results are given for one- and two-stage interference cancellation for equal as well as unequal received power probes.

Divsalar, D.↗