Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated computing”

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,027 records · Page 57

Comparison of the full potential and Euler formulations for computing transonic airfoil flows

A quantitative comparison between the Euler and full potential formulations with respect to speed and accuracy is presented. The robustness of the codes used is tested by a number of transonic airfoil cases. The computed results are from four transonic airfoil computer codes. The full potential codes use fully implicit iteration algorithms. The first Euler code uses a fully implicit ADI iteration scheme. The second Euler code uses an explicit Runge Kutta time stepping algorithm which is enhanced by a multigrid convergence acceleration scheme. Quantitative comparisons are made using various plots of lift coefficient versus the average mesh spacing along the airfoil. Besides yielding an asymptotic limit to the lift coefficient, these results also demonstrate the truncation error behavior of the various codes. Quantitative conclusions regarding the full potential and Euler formulations with respect to accuracy, speed, and robustness can be presented.

Flores, J.↗

Tests of a Gust-alleviating Flap in the Gust Tunnel

Tests were made to determine the effectiveness of a long-period dynamically overbalanced flap in reducing airplane accelerations due to atmospheric gusts. For two gust shapes, one gust velocity, one forward velocity, and one wing loading, a series of flights was made with the flap locked and was then repeated with the flap free to operate. The records were evaluated by routine methods. The results indicate that the flap reduced the maximum acceleration increment 39 percent for a severe gust but with a representative gust shape (a sharp-edge gust), the reduction was only 3 percent. The results also indicate that the flap tended to reduce the longitudinal stability of the airplane. Computations made of the effectiveness and the action of the flap were in good agreement with the experimental results.

Donely, Philip↗

Electron trapping in evolving coronal structures during a large gradual hard X-ray/radio burst

Gradual hard X-ray/radio bursts are characterized by their long duration, smooth time profile, time delays between peaks at different hard X-ray energies and microwaves, and radiation from extended sources in the low and middle corona. Their characteristic properties have been ascribed to the dynamic evolution of the accelerated electrons in coronal magnetic traps or to the separate acceleration of high-energy electrons in a 'second step' process. The information available so far was drawn from quality considerations of time profiles or even only from the common occurrence of emissions in different spectral ranges. This paper presents model computations of the temporal evolution of hard X-ray and microwave spectra, together with a qualitative discussion of radio lightcurves over a wide spectral range, and metric imaging observations. The basis hypothesis investigated is that the peculiar 'gradual' features can be related to the dynamical evolution of electrons injected over an extended time interval in a coronal trap, with electrons up to relativistic energies being injected simultaneously. The analyzed event (26 April. 1981) is particularly challenging to this hypothesis because of the long time delays between peaks at different X-ray energies and microwave frequencies. The observations are shown to be consistent with the hypothesis, provided that the electrons lose their energy by Coulomb collisions and possibly betatron deceleration. The access of the electrons to different coronal structures varies in the course of the event. The evolution and likely destabilization of part of the coronal plasma-magnetic field configuration is of crucial influence in determining the access to these structures and possibly the dynamical evolution of the trapped electrons through betatron deceleration in the late phase of the event.

Bruggmann, G.↗

Automated Hardware-Identification System

"Compressed symbology" emerging technology involving one- and two-dimensional arrays of surface depressions to form optically readable dots. Patterns more durable and denser than common bar codes. Convey identification data in binary form and read by optoelectric sensors. Computers and compressed-symbology engraving machines they control constitute subsystems of "paperless" hardware-tracking and -identification systems coordinating flows of both identifying information and identified parts themselves, along with ancillary information like work orders. Modifications of software expected to accelerate marking operations, eliminate need for trial or practice marking, and reduce incidence of errors.

Schramm, Harry F., Jr.↗

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗

Numerical computation of transonic flows with shock waves

Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.

Jameson, A.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Simulation of the Solar Energetic Particle Event on 2020 May 29 Observed by Parker Solar Probe

This paper presents a stochastic three-dimensional focused transport simulation of solar energetic particles (SEPs) produced by a data-driven coronal mass ejection (CME) shock propagating through a data-driven model of coronal and heliospheric magnetic fields. The injection of SEPs at the CME shock is treated using diffusive shock acceleration of post-shock suprathermal solar wind ions. A time-backward stochastic simulation is employed to solve the transport equation to obtain the SEP time–intensity profile at any location, energy, and pitch angle. The model is applied to a SEP event on 2020 May 29, observed by STEREO-A close to ∼1 au and by Parker Solar Probe (PSP) when it was about 0.33 au away from the Sun. The SEP event was associated with a very slow CME with a plane-of-sky speed of 337 km s −1 at a height below 6 RS as reported in the SOHO/LASCO CME catalog. We compute the time profiles of particle flux at PSP and STEREO-A locations, and estimate both the spectral index of the proton energy spectrum for energies between ∼2 and 16 MeV and the equivalent path length of the magnetic field lines experienced by the first arriving SEPs. We find that the simulation results are well correlated with observations. The SEP event could be explained by the acceleration of particles by a weak CME shock in the low solar corona that is not magnetically connected to the observers.

Solar energetic particles↗

A generalized multiple-input, multiple-output modal parameter estimation algorithm

A new method for experimental determination of the modal parameters of a structure is presented. The method allows for multiple input forces to be applied simultaneously, and for an arbitrary number of acceleration response measurements to be employed. These data are used to form the equations of motion for a damped linear elastic structure. The modal parameters are then obtained through an eigenvalue technique. In conjunction with the development of the equations, an extensive computer simulation study was performed. The results of the study show a marked improvement in the mode shape identification for closely-spaced modes as the number of applied forces is increased. Also demonstrated is the influence of noise on the method's ability to identify accurate modal parameters. Here again, an increase in the number of exciters leads to a significant improvement in the identified parameters.

Craig, R. R., Jr.↗

Current loop coalescence model of solar flares

A computer simulation and theoretical study of the physical characteristics of the explosive coalescence of current-carrying loops is presented. Characteristics of the explosive coalescence include a large impulsive increase of the kinetic energies of electrons and ions, the simultaneous heating and acceleration of electrons and ions in high and low energy ranges, and a break in the energy spectra of electrons and ions. A characteristic double subpeak structure is found in the quasi-periodic oscillations found in the time profiles of the solar flares of June 7, 1980 and November 26, 1982 which can be explained in terms of the coalescence instability of two current loops.

Tajima, T.↗

Particle simulations in magnetospheric plasmas

In view of the recent remarkable advancement of computer technology and simulation software, simulation studies are one of the most powerful academic tools for establishment of quantitative space physics and modelling of our space environment. The complex nature encountered in space plasma physics has motivated considerable development in computer simulations, which have played an essential role in the development of space plasma theory. This report describes research undertaken to understand physical processes involved in plasma waves observed in the magnetospheric plasmas, and associated nonlinear phenomena such as heating, diffusion, and acceleration of particles due to excited waves. The research explains and clarifies the observational data both qualitatively and quantitatively.

Nishikawa, Ken-Ichi↗

Time-optimal Aircraft Pursuit-evasion with a Weapon Envelope Constraint

The optimal pursuit-evasion problem between two aircraft including a realistic weapon envelope is analyzed using differential game theory. Six order nonlinear point mass vehicle models are employed and the inclusion of an arbitrary weapon envelope geometry is allowed. The performance index is a linear combination of flight time and the square of the vehicle acceleration. Closed form solution to this high-order differential game is then obtained using feedback linearization. The solution is in the form of a feedback guidance law together with a quartic polynomial for time-to-go. Due to its modest computational requirements, this nonlinear guidance law is useful for on-board real-time implementation.

Menon, P. K. A.↗

Parallel algorithms for placement and routing in VLSI design

The computational requirements for high quality synthesis, analysis, and verification of very large scale integration (VLSI) designs have rapidly increased with the fast growing complexity of these designs. Research in the past has focused on the development of heuristic algorithms, special purpose hardware accelerators, or parallel algorithms for the numerous design tasks to decrease the time required for solution. Two new parallel algorithms are proposed for two VLSI synthesis tasks, standard cell placement and global routing. The first algorithm, a parallel algorithm for global routing, uses hierarchical techniques to decompose the routing problem into independent routing subproblems that are solved in parallel. Results are then presented which compare the routing quality to the results of other published global routers and which evaluate the speedups attained. The second algorithm, a parallel algorithm for cell placement and global routing, hierarchically integrates a quadrisection placement algorithm, a bisection placement algorithm, and the previous global routing algorithm. Unique partitioning techniques are used to decompose the various stages of the algorithm into independent tasks which can be evaluated in parallel. Finally, results are presented which evaluate the various algorithm alternatives and compare the algorithm performance to other placement programs. Measurements are presented on the parallel speedups available.

Brouwer, Randall Jay↗

Acceleration efficiency in line-driven flows

We reexamine the physics of flow driving by line scattering of a continuum radiation source to determine the degree to which such line scattering can heat as well as accelerate the flow. Within the framework of the Sobolev theory for line transfer, we argue that the finite thermal width of the line scattering profile can lead to a significant 'Doppler heating' via photon frequency redistribution within a Sobolev resonance layer. Quantitative computation of this heating shows, however, that it is largely canceled by a corresponding cooling by the diffuse radiation. The resulting reduction in net Doppler heating or cooling means that the overall effect is only of limited importance in the energy balance of line-driven stellar winds. Through simple scaling relations, we compare the effect to other competing heating or cooling terms, including the ion-drag frictional heating recently discussed by Springmann and Pauldrach. We also provide a physical explanation of the unexpected cooling effect, and comment that its near cancellation of the anticipated heating provides another example of the tendency for ideal Sobolev theory to apply to a higher order than expected.

Gayley, Kenneth G.↗

Dominant takeover regimes for genetic algorithms

The genetic algorithm (GA) is a machine-based optimization routine which connects evolutionary learning to natural genetic laws. The present work addresses the problem of obtaining the dominant takeover regimes in the GA dynamics. Estimated GA run times are computed for slow and fast convergence in the limits of high and low fitness ratios. Using Euler's device for obtaining partial sums in closed forms, the result relaxes the previously held requirements for long time limits. Analytical solution reveal that appropriately accelerated regimes can mark the ascendancy of the most fit solution. In virtually all cases, the weak (logarithmic) dependence of convergence time on problem size demonstrates the potential for the GA to solve large N-P complete problems.

Noever, David↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗

One-Dimensional, Two-Phase Flow Modeling Toward Interpreting Motor Slag Expulsion Phenomena

Aluminum oxide slag accumulation and expulsion was previously shown to be a player in various solid rocket motor phenomena, including the Space Shuttle's Reusable Solid Rocket Motor (RSRM) pressure perturbation, or "blip," and phantom moment. In the latter case, such un ]commanded side accelerations near the end of burn have also been identified in several other motor systems. However, efforts to estimate the mass expelled during a given event have come up short. Either bulk calculations are performed without enough physics present, or multiphase, multidimensional Computational Fluid Dynamic analyses are performed that give a snapshot in time and space but do not always aid in grasping the general principle. One ]dimensional, two ]phase compressible flow calculations yield an analytical result for nozzle flow under certain assumptions. This can be carried further to relate the bulk motor parameters of pressure, thrust, and mass flow rate under the different exhaust conditions driven by the addition of condensed phase mass flow. An unknown parameter is correlated to airflow testing with water injection where mass flow rates and pressure are known. Comparison is also made to full ]scale static test motor data where thrust and pressure changes are known and similar behavior is shown. The end goal is to be able to include the accumulation and flow of slag in internal ballistics predictions. This will allow better prediction of the tailoff when much slag is ejected and of mass retained versus time, believed to be a contributor to the widely-observed "flight knockdown" parameter.

Kibbey, Timothy P.↗

Investigations of Pole Erosion Mechanisms in the 12.5 kW HERMeS Hall Thruster Using Numerical Simulations and Ion Velocity Measurements

When the Hall Effect Rocket with Magnetic Shielding (HERMeS) was tested at the NASA GRC’s VF-5 the erosion rate on the pole covers was found to be at least a factor of two higher at 300 V and 20.8 A than at other operating conditions in the 400 V to 600 V range. Simulations using our hydrodynamics code Hall2De do not predict increased erosion rates at 300 V but accurately compute erosion rates similar to the measurements at the other operating conditions. We investigate the source of the discrepancy between measurements and simulations at 300 V using a combination of numerical simulations and experimental measurements of the ion velocity fields in the acceleration region of the thruster. By examining previous simulations that predicted the measured erosion rates for other thrusters and operating conditions, we determine that sputtering of the pole surfaces by high energy ions is the most likely mechanism behind the erosion rates at 300 V. High energy ions sputter the pole surfaces when the acceleration region of the thruster is downstream of the pole surface plane. When the latter occurs, the curvature of the plasma potential contours at the edges of the channel accelerates a small fraction of the high energy ions radially. We find that neither our Hall2De simulations nor the experimental measurements produce ions of sufficiently high energy and flux to the inner pole cover to explain the measured erosion there. Only at the outer edge of the cover is where we find the ions needed to yield simulation results that are comparable to the erosion measurements. We argue that one possible mechanism for the higher erosion is high energy ions that graze the channel corner of the pole and become trapped in a sheath that develops downstream. Within the sheath, we find that the electric field is large enough to turn the ions towards the surface of the pole. A simplified calculation shows that the erosion rates produced by sputtering of trapped ions are similar to the measurements. We also propose that the presence of local oscillations and high-energy cathode ions may be sources of the enhanced erosion We conclude this article by proposing a combination of analysis, simulation, and experimental measurements that can be used to address the validity of these hypotheses.

Mikellides, Ioannis G.↗