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 379 records · Page 21

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗

Periodic oscillations of the quasi-parallel bow shock as observed with the AMPTE/CCE spacecraft

The AMPTE/CCE spacecraft has, on occasion, encountered the earth's bow shock near apogee (about 9 earth radii). During one encounter, on November 1, 1984, the bow shock was observed to pass back and forth over the spacecraft several times in a periodic fashion. The period of the oscillation was about 20 seconds, and an upstream wave with a similar period was observed in the solar wind just prior to the multiple shock crossings. Minimum variance analysis of the shock crossings shows that the magnetic field carried out a single rotation through 360 deg for each pair of shock crossings. This implies that the shock was responding to the rotation of the field in the upstream wave, and it is shown how the plasma displacement associated with the upstream wave can produce the observed oscillations. These data may provide support for the wavelike structure of quasi-parallel shocks predicted in simulations.

Strangeway, R. J.↗

Time Warp Operating System (TWOS)

Designed to support parallel discrete-event simulation, TWOS is complete implementation of Time Warp mechanism - distributed protocol for virtual time synchronization based on process rollback and message annihilation.

Bellenot, Steven F.↗

A Ka-Band Wide-Bandgap Solid-State Power Amplifier: Architecture Performance Estimates

Motivated by recent advances in wide-bandgap (WBG) gallium nitride (GaN) semiconductor technology, there is considerable interest in developing efficient solidstate power amplifiers (SSPAs) as an alternative to the traveling-wave tube amplifier (TWTA) for space applications. This article documents the results of a study to investigate power-combining technology and SSPA architectures that can enable a 120-W, 40 percent power-added efficiency (PAE) SSPA. Results of the study indicate that architectures based on at least three power combiner designs are likely to enable the target SSPA. The proposed architectures can power combine 16 to 32 individual monolithic microwave integrated circuits (MMICs) with >80 percent combining efficiency. This corresponds to MMIC requirements of 5- to 10-W output power and >48 percent PAE. For the three proposed architectures [1], detailed analysis and design of the power combiner are presented. The first architecture studied is based on a 16-way septum combiner that offers low loss and high isolation over the design band of 31 to 36 GHz. Analysis of a 2-way prototype septum combiner had an input match >25 dB, output match >30 dB, insertion loss <0.02 dB, and isolation >30 dB over the design band. A 16-way design, based on cascading this combiner in a binary fashion, is documented. The second architecture is based on a 24-way waveguide radial combiner. A prototype 24-way radial base was analyzed to have an input match >30 dB (under equal excitation of all input ports). The match of the mode transducer that forms the output of a radial combiner was found to be >27 dB. The functional bandwidth of the radial base and mode transducer, which together will form a radial combiner/divider, exceeded the design band. The third architecture employs a 32-way, parallel-plate radial combiner. Simulation results indicated an input match >24 dB, output match >22 dB, insertion loss <0.23 dB, and adjacent port isolation >20 dB over the design band. All three architectures utilize a low-loss MMIC amplifier module based on commercial MMIC packaging and a custom microstrip-to-rectangular-waveguide transition. The insertion loss of the module is expected to be 0.45 dB over the design band.

Epp, L.↗

Meteorological Controls on Biomass Burning During Santa Ana Events in Southern California

Fires occurring during Santa Ana (SA) events in southern California are driven by extreme fire weather characterized by high temperatures, low humidities, and high wind speeds. We studied the controls on burned area and carbon emissions during two intensive SA burning periods in 2003 and 2007. We therefore used remote sensing data in parallel with fire weather simulations of the Weather and Regional Forecast model. Total carbon emissions were approximately 1800 gigagrams in 2003 and 900 gigagrams in 2007, based on a daily burned area and a fire emission model that accounted for spatial variability in fuel loads and combustion completeness. On a regional scale, relatively strong positive correlations were found between the daily Fosberg fire weather index and burned area/emissions (probability is less than 0.01). Our analysis provides a quantitative assessment of relationships between fire activity and weather during severe SA fires in southern California.

Santa Ana (SA)↗

Experimental and Finite Element Modeling of Near-Threshold Fatigue Crack Growth for the K-Decreasing Test Method

The experimental methods to determine near-threshold fatigue crack growth rate data are prescribed in ASTM standard E647. To produce near-threshold data at a constant stress ratio (R), the applied stress-intensity factor (K) is decreased as the crack grows based on a specified K-gradient. Consequently, as the fatigue crack growth rate threshold is approached and the crack tip opening displacement decreases, remote crack wake contact may occur due to the plastically deformed crack wake surfaces and shield the growing crack tip resulting in a reduced crack tip driving force and non-representative crack growth rate data. If such data are used to life a component, the evaluation could yield highly non-conservative predictions. Although this anomalous behavior has been shown to be affected by K-gradient, starting K level, residual stresses, environmental assisted cracking, specimen geometry, and material type, the specifications within the standard to avoid this effect are limited to a maximum fatigue crack growth rate and a suggestion for the K-gradient value. This paper provides parallel experimental and computational simulations for the K-decreasing method for two materials (an aluminum alloy, AA 2024-T3 and a titanium alloy, Ti 6-2-2-2-2) to aid in establishing clear understanding of appropriate testing requirements. These simulations investigate the effect of K-gradient, the maximum value of stress-intensity factor applied, and material type. A material independent term is developed to guide in the selection of appropriate test conditions for most engineering alloys. With the use of such a term, near-threshold fatigue crack growth rate tests can be performed at accelerated rates, near-threshold data can be acquired in days instead of weeks without having to establish testing criteria through trial and error, and these data can be acquired for most engineering materials, even those that are produced in relatively small product forms.

Smith, Stephen W.↗

Lightforce Photon-Pressure Collision Avoidance: Efficiency Analysis in the Current Debris Environment and Long-Term Simulation Perspective

This work provides an efficiency analysis of the LightForce space debris collision avoidance scheme in the current debris environment and describes a simulation approach to assess its impact on the long-term evolution of the space debris environment. LightForce aims to provide just-in-time collision avoidance by utilizing photon pressure from ground-based industrial lasers. These ground stations impart minimal accelerations to increase the miss distance for a predicted conjunction between two objects. In the first part of this paper we will present research that investigates the short-term effect of a few systems consisting of 20-kilowatt-class lasers directed by 1.5-meter-diameter telescopes using adaptive optics. The results found such a network of ground stations to mitigate more than 85 percent of conjunctions and could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. While these are impressive numbers that indicate LightForce's utility in the short-term, the remaining 15 percent of possible collisions contain (among others) conjunctions between two massive objects that would add large amount of debris if they collide. Still, conjunctions between massive objects and smaller objects can be mitigated. Hence, we choose to expand the capabilities of the simulation software to investigate the overall effect of a network of LightForce stations on the long-term debris evolution. In the second part of this paper, we will present the planned simulation approach for that effort. For the efficiency analysis of collision avoidance in the current debris environment, we utilize a simulation approach that uses the entire Two Line Element (TLE) catalog in LEO for a given day as initial input. These objects are propagated for one year and an all-on-all conjunction analysis is performed. For conjunctions that fall below a range threshold, we calculate the probability of collision and record those values. To assess efficiency, we compare a baseline (without collision avoidance) conjunction analysis with an analysis where LightForce is active. Using that approach, we take into account that collision avoidance maneuvers could have effects on third objects. Performing all-on-all conjunction analyses for extended period of time requires significant computer resources; hence we implemented this simulation utilizing a highly parallel approach on the NASA Pleiades supercomputer.

space debris mitigation↗

Simulation Schiaparelli's Entry and Comparison to Aerothermal Flight Data

The European Space Agency recently flew an entry, descent, and landing demonstrator module called Schiaparelli that entered the atmosphere of Mars on the 19th of October, 2016. The instrumentation suite included heatshield and backshell pressure transducers and thermocouples (known as AMELIA - Atmospheric Mars Entry and Landing Investigations and Analysis) and backshell radiation and direct heat flux-sensing sensors (known as COMARS (Combined Aerothermal and Radiometer Sensors Instrument Package) and ICOTOM (narrow band radiometers)). Due to the failed landing of Schiaparelli, only a subset of the flight data was transmitted before and after plasma black-out. The goal of this paper is to present comparisons of the flight data with calculations from NASA simulation tools, DPLR (Data Parallel Line Relaxation) / NEQAIR (NonEQuilibrium AIr Radiation) and LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) / HARA (High-temperature Aerothermodynamic RAdiation ). DPLR and LAURA are used to calculate the flowfield around the vehicle and surface properties, such as pressure and convective heating. The flowfield data are passed to NEQAIR and HARA to calculate the radiative heat flux. Comparisons will be made to the COMARS total heat flux, radiative heat flux and pressure measurements. Results will also be shown against the reconstructed heat flux which was calculated from an inverse analysis of the AMELIA thermocouple data performed by Astrium. Preliminary calculations are presented in this abstract.

Brandis, A. M.↗

Circumbinary Disk Accretion into Spinning Black Hole Binaries

Supermassive black hole binaries are likely to accrete interstellar gas through a circumbinary disk.Shortly before merger, the inner portions of this circumbinary disk are subject to general relativisticeffects. To study this regime, we approximate the spacetime metric of close orbiting black holes bysuperimposing two boosted Kerr-Schild terms. After demonstrating the quality of this approximation,we carry out very long-term general relativistic magnetohydrodynamic simulations of the circumbinarydisk. We consider black holes with spin dimensionless parameters of magnitude 0.9, in one simula-tion parallel to the orbital angular momentum of the binary, but in another anti-parallel. These arecontrasted with spinless simulations. We find that, for a fixed surface mass density in the inner cir-cumbinary disk, aligned spins of this magnitude approximately reduce the mass accretion rate by 14%and counter-aligned spins increase it by 45%, leaving many other disk properties unchanged.

Federico G Lopez Armengol↗

On the Presence and Thermalization of Cold Ions in the Exhaust of Antiparallel Symmetric Reconnection

Using fully kinetic 2.5 dimensional particle-in-cell simulations of anti-parallel symmetric magnetic reconnection, we investigate how initially cold ions are captured by the reconnection process, and how they evolve and behave in the exhaust. We find that initially cold ions can remain cold deep inside the exhaust. Cold ions that enter the exhaust downstream of active separatrices, closer to the dipolarization front, appear as cold counter-streaming beams behind the front. In the off-equatorial region, these cold ions generate ion-acoustic waves that aid in the thermalization both of the incoming and outgoing populations. Closest to the front, due to the stronger magnetization, the ions can remain relatively cold during the neutral plane crossing. In the intermediate exhaust, the weaker magnetization leads to enhanced pitch angle scattering and reflection. Cold ions that enter the exhaust closer to the X line, at active separatrices, evolve into a thermalized exhaust. Here, the cold populations are heated through a combination of thermalization at the separatrices and pitch angle scattering in the curved magnetic field around the neutral plane. Depending on where the ions enter the exhaust, and how long time they have spent there, they are accelerated to different energies. The superposition of separately thermalized ion populations that have been accelerated to different energies form the hot exhaust population.

Cecilia Norgren↗

Parallelized Carbon Oxidation and Surface Recession Model in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that often consist of a carbon-based reinforcement and a polymeric matrix. During Earth re-entry, they endure high-temperature oxidation and surface recession. Oxidation is an important mechanism for ablation, sometimes leading to the weakening, spallation, or failure of the oxidized fibers at the surface and in the char layer. However, more details are required including accurate material properties of the fiber microstructure, whether the fibers recede homogenously or localized at pits, and the role of pyrolysis outgassing in order to evaluate the role oxidation plays in the degradation and failure mechanisms of these materials. In this work, we demonstrate a parallelizable oxidation-driven ablation model developed for detailed, large-scale simulations in the DSMC code SPARTA. We also develop robust algorithms which handle the conservation of the surface state after an ablation step. Finally, we verify our model for both simple and more complex chemistry as well as microstructures with reaction statistics, oxidation depth calculations, and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

High-performance computing for flight vehicles; Proceedings of the Symposium, Washington, Dec. 7-9, 1992

The present conference discusses high-performance computing systems for flight vehicles, large-scale simulations on high-performance flight computers and software, multidisciplinary and design/optimization applications of computers, computational electromagnetics and acoustics, the simulation of aircraft powerplant turbomachinery and reacting flows, and flow calculations on parallel machines. Also discussed are direct flow simulation Monte Carlo methods, structural mechanics sensitivity and fracture calculations on parallel machines, grid-generation and advanced algorithms for CFD, advanced solid-mechanics and structures applications, and advancements in flow visualization technology and neural networks.

Noor, Ahmed K.↗

Ion-beam-driven electrostatic ion cyclotron instabilities

Results are presented of a particle simulation study of the electrostatic ion-cyclotron (EIC) instability driven by a parallel ion beam. The results of this simulation study demonstrate the nonlinear consequences of nonresonant EIC waves destabilized by an ion beam parallel to the magnetic field for the case of a large beam velocity. As a consequence of the instability, it is shown that the beam ions are heated strongly in the perpendicular direction and suffer a strong anomalous friction via EIC waves which leads to the beam slowing down. Simulation results indicate that the anomalous slowing down of beam ions by EIC waves is much larger than that from the classical electron-drag, and the perpendicular collision frequency measured from perpendicular beam heating is as large as that from Bohm diffusion. It is concluded that the ion beam driven EIC wave is a viable mechanism for the transfer of ion parallel beam energy to the ion perpendicular energy.

Okuda, H.↗

Cyclic behavior at quasi-parallel collisionless shocks

Large scale one-dimensional hybrid simulations with resistive electrons have been carried out of a quasi-parallel high-Mach-number collisionless shock. The shock initially appears stable, but then exhibits cyclic behavior. For the magnetic field, the cycle consists of a period when the transition from upstream to downstream is steep and well defined, followed by a period when the shock transition is extended and perturbed. This cyclic shock solution results from upstream perturbations caused by backstreaming gyrating ions convecting into the shock. The cyclic reformation of a sharp shock transition can allow ions, at one time upstream because of reflection or leakage, to contribute to the shock thermalization.

Burgess, D.↗

A multiarchitecture parallel-processing development environment

A description is given of the hardware and software of a multiprocessor test bed - the second generation Hypercluster system. The Hypercluster architecture consists of a standard hypercube distributed-memory topology, with multiprocessor shared-memory nodes. By using standard, off-the-shelf hardware, the system can be upgraded to use rapidly improving computer technology. The Hypercluster's multiarchitecture nature makes it suitable for researching parallel algorithms in computational field simulation applications (e.g., computational fluid dynamics). The dedicated test-bed environment of the Hypercluster and its custom-built software allows experiments with various parallel-processing concepts such as message passing algorithms, debugging tools, and computational 'steering'. Such research would be difficult, if not impossible, to achieve on shared, commercial systems.

Townsend, Scott↗

Near-specular reflection of ions at quasi-parallel shocks

One-dimensional hybrid simulations and a semianalytical model of the shock front are employed to investigate the source regions in the incident ion phase space of reflected and transmitted ions, the relative importance of the electric and magnetic forces in the reflection of incident ions, and how these characteristics change in time. The phase space origin of reflected particles and the fraction of incident ions that reflect are found to depend on the electromagnetic field structure of the shock front at the time the ions encounter it. The reflection fraction is maximized when the electric field along the shock normal (Ex) and the noncoplanar magnetic field (By) are at their maximum values. When Ex and By are large, the reflection process produces a beam that is cooler, more dense, and closer to specular than when they are small.

Mckean, M. E.↗