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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Numerical study of vortical flow over a sideslipping delta wing

The three-dimensional vortical viscous flow past a sideslipping delta wing at a high angle of attack is numerically investigated. A computational method for calculating low-speed viscous flowfields is developed. The time-iterative method uses an implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. An algebraic turbulence model is implemented to account for the eddy viscosity. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. Computed results are compared with experimental data.

Hsu, C.-H.↗

NASA NTRS1991

Computational aspects of sensitivity calculations in linear transient structural analysis

The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, W. H.↗

NASA NTRS2005

Studies of Cosmic Ray Modulation and Energetic Particle Propagation in Time-Dependent 3-Dimensional Heliospheric Magnetic Fields

The primary goal of this project was to perform theoretical calculations of propagation of cosmic rays and energetic particles in 3-dimensional heliospheric magnetic fields. We used Markov stochastic process simulation to achieve to this goal. We developed computation software that can be used to study particle propagation in, as two examples of heliospheric magnetic fields that have to be treated in 3 dimensions, a heliospheric magnetic field suggested by Fisk (1996) and a global heliosphere including the region beyond the termination shock. The results from our model calculations were compared with particle measurements from Ulysses, Earth-based spacecraft such as IMP-8, WIND and ACE, Voyagers and Pioneers in outer heliosphere for tests of the magnetic field models. We particularly looked for features of particle variations that can allow us to significantly distinguish the Fisk magnetic field from the conventional Parker spiral field. The computer code will eventually lead to a new generation of integrated software for solving complicated problems of particle acceleration, propagation and modulation in realistic 3-dimensional heliosphere of realistic magnetic fields and the solar wind with a single computation approach.

Zhang, Ming↗

NASA NTRS1982

A vectorized, finite-volume, adaptive-grid algorithm for Navier-Stokes calculations

An adaptive grid, finite-volume method has been used to solve the Navier-Stokes equations for complete (forebody and afterbody) flowfields around blunt bodies. The code, which is applicable for axisymmetric or two-dimensional flows, allows the mesh to adjust during the computation to provide a closer spacing of mesh points in regions of high gradients, thus minimizing the number of required computational points. The solution technique is explicit, utilizing a maximum time-step advancement at each grid point to accelerate convergence to the steady state. The code has been fully vectorized for efficient solution on the CYBER 203 computer. A very flexible rezoning routine is used to concentrate mesh points anywhere in the field, either by a user-defined weighting function or by allowing high gradient regions to adjust the grid. The grid adjustment routine is implicit in nature and represents a very small portion of the total computational cost. Currently, the code runs in approximately 0.000016 seconds per grid point per iteration.

Gnoffo, P. A.↗

NASA NTRS1993

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics, part 2

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗

NASA NTRS1992

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallelize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet -- and perhaps will never -- be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗

NASA NTRS1995

Modeling the vestibulo-ocular reflex of the squirrel monkey during eccentric rotation and roll tilt

Model simulations of the squirrel monkey vestibulo-ocular reflex (VOR) are presented for two motion paradigms: constant velocity eccentric rotation and roll tilt about a naso-occipital axis. The model represents the implementation of three hypotheses: the "internal model" hypothesis, the "gravito-inertial force (GIF) resolution" hypothesis, and the "compensatory VOR" hypothesis. The internal model hypothesis is based on the idea that the nervous system knows the dynamics of the sensory systems and implements this knowledge as an internal dynamic model. The GIF resolution hypothesis is based on the idea that the nervous system knows that gravity minus linear acceleration equals GIF and implements this knowledge by resolving the otolith measurement of GIF into central estimates of gravity and linear acceleration, such that the central estimate of gravity minus the central estimate of acceleration equals the otolith measurement of GIF. The compensatory VOR hypothesis is based on the idea that the VOR compensates for the central estimates of angular velocity and linear velocity, which sum in a near-linear manner. During constant velocity eccentric rotation, the model correctly predicts that: (1) the peak horizontal response is greater while "facing-motion" than with "back-to-motion"; (2) the axis of eye rotation shifts toward alignment with GIF; and (3) a continuous vertical response, slow phase downward, exists prior to deceleration. The model also correctly predicts that a torsional response during the roll rotation is the only velocity response observed during roll rotations about a naso-occipital axis. The success of this model in predicting the observed experimental responses suggests that the model captures the essence of the complex sensory interactions engendered by eccentric rotation and roll tilt.

NASA Discipline Neuroscience↗

NASA NTRS2013

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi-Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post-landing segment of the Flight Software (FSW). An in-depth trade study was completed to determine optimal characteristics of the touchdown detection method resulting in an algorithm monitoring filtered, lever-arm corrected, 200 Hz Inertial Measurement Unit (IMU) vehicle acceleration magnitude data against a tunable threshold using persistence counter logic. Following the design of the algorithm, high fidelity environment and vehicle simulations, coupled with the actual vehicle FSW, were used to tune the acceleration threshold and persistence counter value to result in adequate performance in detecting touchdown and sufficient safety margin against early detection while descending under parachutes. An analytical approach including Kriging and adaptive sampling allowed for a sufficient number of finite element analysis (FEA) impact simulations to be completed using minimal computation time. The combination of a persistence counter of 10 and an acceleration threshold of approximately 57.3 ft/s2 resulted in an impact performance factor of safety (FOS) of 1.0 and a safety FOS of approximately 2.6 for touchdown declaration. An RCS termination acceleration threshold of approximately 53.1 ft/s(exp)2 with a persistence counter of 10 resulted in an increased impact performance FOS of 1.2 at the expense of a lowered under-parachutes safety factor of 2.2. The resulting tuned algorithm was then tested on data from eight Capsule Parachute Assembly System (CPAS) flight tests, showing an experimental minimum safety FOS of 6.1. The formulated touchdown detection algorithm will be flown on the Orion MPCV FSW during the Exploration Flight Test 1 (EFT-1) mission in the second half of 2014.

Daum, Jared↗

NASA NTRS1978

Lunar tidal acceleration obtained from satellite-derived ocean tide parameters

Observation equations for the M2 ocean tide are computed from Geos 3 data for the long periodic variations of the inclination and node of the orbit. M2 ocean tide parameter values C22+ = 3.23 + or - 0.25 cm, epsilon 22+ = 331 + or - 6 deg, and epsilon 42+ = 113 + or - 6 deg are determined. With the assumption of zero solid tide phase lag, the lunar tidal acceleration is mostly (85%) due to the C22+ term in the expansion of the M2 tide with additional small contributions from the O1 and N2 tides. The calculated value for the tidal acceleration in lunar longitude is -27.4 + or - 3 arc sec/sq (100 yr) which is similar to values determined from astronomical data. The mean elements of Geos 3 are presented in tabular form.

Goad, C. C.↗

NASA NTRS2024

Improved Modelling of SEP Event Onset Within the WSA–Enlil–SEPMOD Framework

Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain – usually set at a heliocentric distance of 20–30 R⊙. This means that particle acceleration at coronal mass ejection (CME)-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a “fixed source” of particle injection at the outer boundary of the coronal domain of the coupled WSA–Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is Solar Energetic Particle MODel (SEPMOD), a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD’s fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere – Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft.

Solar eruptions↗

NASA NTRS1994

Peculiar Behavior of He II Disturbances due to Sloshing Dynamics Driven by Jitter Accelerations Associated with Slew Motion in Microgravity

The mathematical formulation of Dewar container sloshing dynamics when partially filled with superfluid liquid He 11 driven by the gravity jitter acceleration associated with slew motion in a microgravity environment is studied. Numerical computation of the sloshing dynamics is based on the non-inertia container bound co-ordinate, and solving time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. This study discloses the capillary effect of the sloshing dynamics governing the liquid-vapour interface fluctuations driven by the gravity jitter acceleration associated with slew motion, which affects the stability of the fluid system in a microgravity environment. The peculiar behaviour of superfluid helium when reacting to the sloshing dynamics driven by the gravity jitter acceleration associated with slew motion in a microgravity environment is also investigated.

Hung, R. J.↗

NASA NTRS1976

Condensation of water vapor in rarefaction waves. I - Homogeneous nucleation

A detailed theoretical investigation has been made of the condensation of water vapor/carrier gas mixtures in the nonstationary rarefaction wave generated in a shock tube. It is assumed that condensation takes place by homogeneous nucleation. The equations of motion together with the nucleation rate and the droplet growth equations were solved numerically by the method of characteristics and Lax's method of implicit artificial viscosity. It is found that, for the case considered, the condensation wave formed by the collapse of the metastable nonequilibrium state is followed by a shock wave generated by the intersection of characteristics of the same family. The expansion is practically isentropic up to the onset of condensation. The condensation front accelerates in the x,t plane. The results of the computations for a chosen case of water vapor/nitrogen mixture are presented by plotting variations of pressure, nucleation rate, number density of critical clusters, and condensate mass-fraction along three particle paths. Some consideration is given to homogeneous condensation experiments conducted in a shock tube. Although a direct comparison of the present theoretical work and these experiments is not possible, several worthwhile interpretative features have resulted nevertheless.

Sislian, J. P.↗

NASA NTRS1978

Platform motion for fighter simulations

The reported investigation is concerned with aspects of 'realistic' cue generation and evaluation. 'Realistic' cues are defined to mean cues that are representative of an aircraft, that are of sufficient magnitude to be sensed, that are in the proper direction arriving at the proper time, and yet that are within the presentation capabilities of the washout/synergistic platform system. Attention is given to the aircraft mathematical model, computer implementation, cockpit hardware, visual display, the motion system, vertical acceleration, longitudinal acceleration and pitch rate, and lateral acceleration and roll rate. The motion cues for fighter aircraft produced by Langley nonlinear adaptive washout/synergistic motion base system are found to be unacceptable in the roll axis representation. However, a method for producing potentially realistic roll cues is identified.

Parrish, R. V.↗

NASA NTRS1984

Survey of NASA research on crash dynamics

Ten years of structural crash dynamics research activities conducted on general aviation aircraft by the National Aeronautics and Space Administration (NASA) are described. Thirty-two full-scale crash tests were performed at Langley Research Center, and pertinent data on airframe and seat behavior were obtained. Concurrent with the experimental program, analytical methods were developed to help predict structural behavior during impact. The effects of flight parameters at impact on cabin deceleration pulses at the seat/occupant interface, experimental and analytical correlation of data on load-limiting subfloor and seat configurations, airplane section test results for computer modeling validation, and data from emergency-locator-transmitter (ELT) investigations to determine probable cause of false alarms and nonactivations are assessed. Computer programs which provide designers with analytical methods for predicting accelerations, velocities, and displacements of collapsing structures are also discussed.

Thomson, R. G.↗

NASA NTRS2002

Non-Rocket Missile Rope Launcher

The method, installation, and estimation for delivering payload and missiles into outer space are presented. This method uses, in general, the engines and straight or closed-loop cables disposed on a planet surface. The installation consists of a space apparatus, power drive stations located along trajectory of the apparatus, the cables connected to the apparatus and to the power stations, a system for suspending the cable, and disconnected device. The drive stations accelerate the apparatus up to hypersonic speed. The estimations and computations show the possibility of making these projects a reality in a short period of time (see attached project: launcher for missiles and loads). The launch will be very cheap $1-$2 per LB. We need only light strong cable, which can be made from artificial fibers, whiskers, nanotubes, which exist in industry and scientific laboratories.

Bolonkin, Alexander↗

NASA NTRS2004

Test of the Equivalence Principle in an Einstein Elevator

The scientific goal of the experiment is to test the equality of gravitational and inertial mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate of fall of bodies from their compositions. The measurement is accomplished by measuring the relative displacement (or equivalently acceleration) of two falling bodies of different materials which are the proof masses of a differential accelerometer spinning about an horizontal axis to modulate a possible violation signal. A non-zero differential acceleration appearing at the signal frequency will indicate a violation of the Equivalence Principle. The goal of the experiment is to measure the Eotvos ratio og/g (differential acceleration/common acceleration) with a targeted accuracy that is about two orders of magnitude better than the state of the art (presently at several parts in 10(exp 13). The analyses carried out during this first grant year have focused on: (1) evaluation of possible shapes for the proof masses to meet the requirements on the higher-order mass moment disturbances generated by the falling capsule; (2) dynamics of the instrument package and differential acceleration measurement in the presence of errors and imperfections; (3) computation of the inertia characteristic of the instrument package that enable a separation of the signal from the dynamics-related noise; (4) a revised thermal analysis of the instrument package in light of the new conceptual design of the cryostat; (5) the development of a dynamic and control model of the capsule attached to the gondola and balloon to define the requirements for the leveling mechanism (6) a conceptual design of the leveling mechanism that keeps the capsule aligned before release from the balloon; and (7) a new conceptual design of the customized cryostat and a preliminary valuation of its cost. The project also involves an international cooperation with the Institute of Space Physics (IFSI) in Rome, Italy. The group at IFSI is in charge of prototyping the differential accelerometer and carrying out precursor laboratory measurements. During this grant year, our partners analyzed and then designed a new prototype of differential accelerometer that has several characteristics in common with the flight accelerometer at this point of the instrument development. The highlights of these activities are documented in a section of this report.

Shapiro, Irwin I.↗

NASA NTRS2009

Validation of Finite Element Crash Test Dummy Models for Predicting Orion Crew Member Injuries During a Simulated Vehicle Landing

A series of crash tests were conducted with dummies during simulated Orion crew module landings at the Wright-Patterson Air Force Base. These tests consisted of several crew configurations with and without astronaut suits. Some test results were collected and are presented. In addition, finite element models of the tests were developed and are presented. The finite element models were validated using the experimental data, and the test responses were compared with the computed results. Occupant crash data, such as forces, moments, and accelerations, were collected from the simulations and compared with injury criteria to assess occupant survivability and injury. Some of the injury criteria published in the literature is summarized for completeness. These criteria were used to determine potential injury during crew impact events.

Tabiei, Al↗

NASA NTRS1992

Optimal Estimation of Rain Rate Profiles from Single-Frequency Radar Echoes

The significant ambiguities inherent in the determination of a particular vertical rain intensity profile from a given time profile of radar echo powers measured by a downward-looking (spaceborne or airborne) radar at a single attenuating frequency are well-documented. Indeed, one already knows that by appropriately varying the parameters of the reflectivity-rain-rate (Z - R) and/or attenuation-rain-rate (k - R) relationships, one can produce several substantially different hypothetical rain rate profiles which would have the same radar power profile. Imposing the additional constraint that the path-averaged rain-rate be a given fixed number does reduce the ambiguities but falls far short of eliminating them. While we now know how to generate as many mutually ambiguous rain-rate profiles from a given profile of received radar reflectivities as we like, there remains to produce a quantitative measure to assess how likely each of these profiles is, what the appropriate 'average' profile should be, and what the 'variance' of these multiple solutions is. Of course, in order to do this, one needs to spell out the stochastic constraints that can allow us to make sense of the words 'average' and 'variance' in a mathematically rigorous way. Such a quantitative approach would be particularly well-suited for such systems as the proposed Precipitation Radar of the Tropical Rainfall Measuring Mission (TRMM). Indeed, one would then be able to use the radar reflectivities measured by the TRMM radar from one particular look in order to estimate the most likely rain-rate profile that would have produced the measurements, as well as the uncertainty in the estimated rain-rates as a function of range. Such an optimal approach is described in this paper.

deep space flight computer real-time advanced spac↗