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,063 records · Page 59

Upwind scheme for solving the Euler equations on unstructured tetrahedral meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a multidimensional linear reconstruction process. The solution gradients required for the higher-order differenes are computed by a novel approach that yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicitly applying a limiter. Solutions are advanced in time by a three-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Transonic solutions are presented for two meshes around the ONERA M6 wing and demonstrate substantial accuracy and insensitivity to mesh size.

Frink, Neal T.↗

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗

Superfluid Helium On-Orbit Transfer (SHOOT) flight demonstration

The Superfluid Helium On-Orbit Transfer (SHOOT) Flight Demonstration was an attached Shuttle payload mounted on a Hitchhiker cross-bay carrier which flew on STS-57 in June of 1993. SHOOT successfully demonstrated the handling and transfer of superfluid helium between two containers, called dewars, in low gravity. SHOOT was a class C payload and for the STS-57 mission was termed a complex secondary payload. The primaries were the retrieval of the EURECA carrier and a collection of modular experiments contained in SPACEHAB. Because the liquid helium was continuously boiling off, SHOOT's activities were scheduled for the first three days of the mission, concurrent with some SPACEHAB experiments, but well before the EURECA retrieval. Control of the SHOOT experiment was highly interactive and originated primarily from the Goddard Payload Operations and Control Center (POCC). Transfer and calibration activities required continuous command windows of up to 50 minutes duration and up to 80 minutes out of each orbit. Occasionally the crew controlled the experiment using the Payload General Support Computer (PGSC) when near-real time control and monitoring was required. SHOOT also placed considerable demands on the orbiter, including a pitch rotation of 3 deg./sec for 15 minutes, and translational burns using both the aft and forward RCS jets to generate accelerations up to 7 milli-g. The basis for these and other requirements are discussed. Interacion with the crew and timing of crew activity during the mission will be detailed. The processing flow of SHOOT at KSC is described with emphasis on the tradeoffs for vertical, as opposed to horizontal, installation in the orbiter. Finally, some lessons learned are presented that are relevant to future cryogenic and Hitchhiker payloads.

Shirron, Peter↗

Superfluid helium sloshing dynamics induced oscillations and fluctuations of angular momentum, force and moment actuated on spacecraft driven by gravity gradient or jitter acceleration associated with slew motion

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by the gravity gradient and jitter accelerations associated with slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with the Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid-vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers have also been derived. Examples are also given for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with slew motion.

Hung, R. J.↗

Direct Replacement of Arbitrary Grid-Overlapping by Non-Structured Grid

A new approach that uses nonstructured mesh to replace the arbitrarily overlapped structured regions of embedded grids is presented. The present methodology uses the Chimera composite overlapping mesh system so that the physical domain of the flowfield is subdivided into regions which can accommodate easily-generated grid for complex configuration. In addition, a Delaunay triangulation technique generates nonstructured triangular mesh which wraps over the interconnecting region of embedded grids. It is designed that the present approach, termed DRAGON grid, has three important advantages: eliminating some difficulties of the Chimera scheme, such as the orphan points and/or bad quality of interpolation stencils; making grid communication in a fully conservative way; and implementation into three dimensions is straightforward. A computer code based on a time accurate, finite volume, high resolution scheme for solving the compressible Navier-Stokes equations has been further developed to include both the Chimera overset grid and the nonstructured mesh schemes. For steady state problems, the local time stepping accelerates convergence based on a Courant - Friedrichs - Leury (CFL) number near the local stability limit. Numerical tests on representative steady and unsteady supersonic inviscid flows with strong shock waves are demonstrated.

Kao, Kai-Hsiung↗

Numerical studies of the surface tension effect of cryogenic liquid helium

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose of performing scientific observation during the normal spacecraft operation is investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics has been based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers, have been derived.

Hung, R. J.↗

NASA GPU Hackathon Yields Significant Code Improvements

The NASA GPU Hackathon 2020 brought together application developers and computer experts to help get important NASA applications running effectively on graphics processing unit (GPU) nodes. Nine teams of application developers participated in this virtual event, a major impetus for teams to modernize codes of interest for NASA missions to CPU nodes containing GPU accelerators, with a focus on hands-on problem solving. The photo in Figure1 shows 30 of the more than50 participants. The HECC project and NVIDIA jointly organized the event, and HECC provided five Pleiades nodes each with 4 V100 GPUs for teams to use. The virtual event, which took place over four days from September 28–October 7, 2020, used Microsoft Teams and Slack as collaboration tools. Each team consisted of three to six members from NASA Centers and supporting organizations. The teams were paired with one to two mentors from industry, government, and academia. The experience levels of the teams ranged from being GPU novices to advanced CUDA programming experts. OpenACC and the emerging Kokkos API were used in addition to CUDA for GPU programming. During the event, which focused on accelerating AeroSciences and CFD applications, most teams achieved considerable performance improvements on both GPUs and CPUs. For example, a team with no GPU experience completed a first port of a time-critical loop to a GPU. Another team of expert CUDA programmers were able to restructure their algorithm, yielding a factor of five speed-up. And another team sped up some of their CUDA kernels by a factor of 20, which directly translated into their production code. This article highlights some of the many successes resulting from the event.

HECC↗

Coilgun Acceleration Model Containing Interactions Between Multiple Coils

Electromagnetic (EM) accelerators have the potential to fill a performance range not currently being met by conventional chemical and electric propulsion systems by providing a specific impulse of 600-1000 seconds and a thrust-to-power ratio greater than 200 mN/kW. A propulsion system based on EM acceleration of small projectiles has the traditional advantages of using a pulsed system, including precise control over a range of thrust and power levels as well as rapid response and repetition rates. Furthermore, EM accelerators have lower power requirements than conventional electric propulsion systems since no plasma creation is necessary. A coilgun is a specific type of EM device where a high-current pulse through a coil of wire interacts with a conductive projectile via an induced magnetic field to accelerate the projectile. There are no physical or electrical connections to the projectile, which leads to less system degradation and a longer life expectancy. Multi-staging a coilgun by adding multiple turns on a single coil or on the projectile increases the inductance, thus permitting acceleration of the projectile to higher velocities. Previously, a simplified problem of modeling an inductively-coupled, single-coil coilgun using a circuit-based analysis coupled to the one-dimensional momentum equation through Lenz's law was solved; however, the analysis was only conducted on uncoupled coils. The problem is significantly more complicated when multiple, independently-powered coils simultaneously operate and interact with each other and the projectile through induced magnetic fields. This paper presents a multi-coil model developed with the magnetostatic finite element solver QuickField. In the model, mutual inductance values between pairs of conductors were found by first computing the magnetic field energy for different cases where individual coils or multiple coils carry current, then integrating over the entire finite element domain for each case, and finally using the definition of inductive energy storage to solve for the self and mutual inductance. The electric circuit model is coupled to the projectile through Lenz's law, with the coils coupled through mutual inductance but able to be independently triggered at different times to optimize the acceleration profile. This initial model to predict the behavior of a projectile's acceleration through a coupled, multi-coil coilgun increases the potential of building a highly efficient coilgun thruster with key advantages over other EM thruster systems, thus making it a promising candidate for satellite main propulsion or attitude control thrusters.

Liu, Connie↗

Decentralized adaptive control of manipulators - Theory, simulation, and experimentation

The author presents a simple decentralized adaptive-control scheme for multijoint robot manipulators based on the independent joint control concept. The control objective is to achieve accurate tracking of desired joint trajectories. The proposed control scheme does not use the complex manipulator dynamic model, and each joint is controlled simply by a PID (proportional-integral-derivative) feedback controller and a position-velocity-acceleration feedforward controller, both with adjustable gains. Simulation results are given for a two-link direct-drive manipulator under adaptive independent joint control. The results illustrate trajectory tracking under coupled dynamics and varying payload. The proposed scheme is implemented on a MicroVAX II computer for motion control of the three major joints of a PUMA 560 arm. Experimental results are presented to demonstrate that trajectory tracking is achieved despite coupled nonlinear joint dynamics.

Seraji, Homayoun↗

Modelling the solution growth of TGS crystals in low gravity

The experimental growth of triglycine sulfate (TGS) crystals from aqueous solution is modeled here in two dimensions using the PHOENICS finite volume code. Simulations are carried out for steady, impulsive, and periodic accelerations in order to determine tolerable acceleration levels. Scaling arguments are used to estimate the times required for thermal and solutal variations from the initial equilibrium state to be diffusively transported throughout the system, and to obtain order of magnitude information on the relative magnitudes of diffusive and convective transport. The computed concentration fields reflect the features of the concentration distributions found experimentally during experiments conducted aboard Spacelab 3 in 1985.

Nadarajah, Arunan↗

A computational study of advanced exhaust system transition ducts with experimental validation

The current study is an application of CFD to a 'real' design and analysis environment. A subsonic, three-dimensional parabolized Navier-Stokes (PNS) code is used to construct stall margin design charts for optimum-length advanced exhaust systems' circular-to-rectangular transition ducts. Computer code validation has been conducted to examine the capability of wall static pressure predictions. The comparison of measured and computed wall static pressures indicates a reasonable accuracy of the PNS computer code results. Computations have also been conducted on 15 transition ducts, three area ratios, and five aspect ratios. The three area ratios investigated are constant area ratio of unity, moderate contracting area ratio of 0.8, and highly contracting area ratio of 0.5. The degree of mean flow acceleration is identified as a dominant parameter in establishing the minimum duct length requirement. The effect of increasing aspect ratio in the minimum length transition duct is to increase the length requirement, as well as to increase the mass-averaged total pressure losses. The design guidelines constructed from this investigation may aid in the design and manufacture of advanced exhaust systems for modern fighter aircraft.

Wu, C.↗

Physics of solar activity

The aim of the research activity was to increase our understanding of solar activity through data analysis, theoretical analysis, and computer modeling. Because the research subjects were diverse and many researchers were supported by this grant, a select few key areas of research are described in detail. Areas of research include: (1) energy storage and force-free magnetic field; (2) energy release and particle acceleration; (3) radiation by nonthermal electrons; (4) coronal loops; (5) flare classification; (6) longitude distributions of flares; (7) periodicities detected in the solar activity; (8) coronal heating and related problems; and (9) plasma processes.

Sturrock, Peter A.↗

Horizontal tail loads in maneuvering flight

A method is given for determining the horizontal tail loads in maneuvering flight. The method is based upon the assignment of a load-factor variation with time and the determination of a minimum time to reach peak load factor. The tail load is separated into various components. Examination of these components indicated that one of the components was so small that it could be neglected for most conventional airplanes; therefore, the number of aerodynamic parameters needed in this computation of tail loads was reduced to a minimum. In order to illustrate the method, as well as to show the effect of the main variables, a number of examples are given. Some discussion is given regarding the determination of maximum tail loads, maximum pitching accelerations, and maximum pitching velocities obtainable.

Pearson, Henry A↗

Axisymmetric Time-Dependent Computations of Expansion Tube Flow

The goal of this work is to add insight about the flow within expansion tubes by using computational fluid dynamics. This is accomplished by comparing the results of axisymmetric numerical simulations with finite-rate chemistry to data from the HYPULSE expansion tube facility which was previously the NASA Langley expansion tube. The numerical simulations begin at the opening of the primary diaphragm and compute the flow throughout the whole facility and, thus, are able to follow and assess the effect of many of the flow features created during operation of the facility. One particular issue that will be investigated is the effect of boundary layer formation in the acceleration tube on the test gas volume and test gas conditions. Both laminar and turbulent boundary layers will be implemented. The effect of momentary shock reflection off the secondary diaphragm will also be investigated. There is concern that such a reflection will stagnate the test gas and create high levels of dissociated molecules. This is particularly important in propulsion experiments where a freestream composition different from flight conditions may influence ignition and burning data. Several different models of diaphragm rupture will be implemented in order to help understand the importance of this issue.

Wilson, Gregory J.↗

Numerical Simulation of Two-grid Ion Optics Using a 3D Code

A three-dimensional ion optics code has been developed under NASA's Project Prometheus to model two grid ion optics systems. The code computes the flow of positive ions from the discharge chamber through the ion optics and into the beam downstream of the thruster. The rate at which beam ions interact with background neutral gas to form charge exchange ions is also computed. Charge exchange ion trajectories are computed to determine where they strike the ion optics grid surfaces and to determine the extent of sputter erosion they cause. The code has been used to compute predictions of the erosion pattern and wear rate on the NSTAR ion optics system; the code predicts the shape of the eroded pattern but overestimates the initial wear rate by about 50%. An example of use of the code to estimate the NEXIS thruster accelerator grid life is also presented.

ion optics↗

Improved Pyrotechnic Shock Scaling Method Based on Shock Response Spectrum and Statistical Energy Analysis

It has been more than forty years since the pyrotechnic shock scaling method was introduced. The scaling method estimates the attenuation of the Shock Response Spectrum (SRS) based on the distance from the source, structural configurations, types of structural joints and interfaces, and intervening structure [1]. The method has been successfully used in the spacecraft community countless times and still is frequently used to develop pyrotechnic shock requirements at various levels of assembly of a spacecraft. However, since the method was derived empirically from a limited set of shock test data [2], the aerospace community has been looking for an alternative approach. In this paper, a computational scaling method based on Statistical Energy Analysis (SEA) is re-visited. SEA is traditionally used as a method for investigating the diffusion of acoustic and vibratory energies of a system at steady-state [3]. Previously, an approach using the SRS as an acceleration constraint condition to the SEA model was introduced to estimate the attenuation from a shock source [4,5]. In the current investigation, additional examples are provided to further validate the approach. In addition, responses in the time domain are produced by the Local Modal Phase Reconstruction (LMPR) approach which was recently introduced to the community [6].

Fernandez, Juan P.↗

Analysis of Launch Vehicle Liftoff Debris: Historical Perspective from Space Shuttle and Application to Artemis I

Human exploration-class launch vehicles are inherently prone to debris due to the extreme environments generated during pre-launch operations, liftoff, and flight. The use of cryogenic propellants often requires thermal protection system (TPS) coatings, typically foam, to maintain the propellant conditions in the tank and prevent an accumulation ice on the external surface of the vehicle. Some ice growth is to be expected at umbilical interfaces, vents, flanges, or brackets where it is difficult to apply TPS. This ice may come loose at any time due to wind on the launch pad, structural vibration and acoustics after rocket ignition, or aerodynamic forces during flight. This phenomena is especially apparent on vehicles with no TPS, such as the Saturn V rockets used in the Apollo Program, see Figure 1. During propellant tanking, the thermal contraction of the underlying substrate may generate cracks in the TPS (Figure 1). Chunks of TPS can release due to the expansion of ingested gas from cryopumping or from aerodynamic forces if the crack creates an offset surface. Most foams will also have a certain amount of “popcorning” where small pieces of foam will pop off during flight because of the differential between the static surface pressure and the pressure of the gas trapped in the foam cell structure. There are a number of other coating or closeout materials that may be shed from the vehicle and become debris. During pre-launch operations and liftoff, the vehicle may also be exposed to debris originating from the launch pad or ground support equipment. This debris is separate from foreign object debris, or FOD, which is not intended to be present and is strictly controlled through operations and maintenance procedures. In this case, debris is generated from hardware and materials that are necessary for launch and are subject to the intense vibration, acoustics, and direct plume impingement of the launch environment. Examples include ice from umbilicals, tape and tie wraps that protect cables, and rust or corrosion from the launch platform. While NASA has historically been aware of debris as a potential issue that could cause a failure resulting in loss of mission, loss of vehicle, or loss of crew, the likelihood and severity of that risk was not always well understood or given sufficient weight in program and flight decisions. After the Space Shuttle Columbia accident (STS-107), the investigation found that foam TPS debris shed from the external tank was the proximate cause of the damage to the orbiter wing. Six previous observations of debris released from the foam ramp that covered the bipod connecting the forward end of the orbiter to the external tank resulted in minor changes or were determined to be accepted flight risks. Two occurrences of bipod ramp foam loss were not identified until the STS-107 investigation. Despite the damage inflicted by these debris strikes, the Shuttle Program Requirements Control Board deemed the vehicle safe to fly. During the Return to Flight effort following the Columbia disaster, NASA Engineering developed a process for the assessment of debris transport, impact, and damage tolerance to support independent assessments of risk by NASA Safety and Mission Assurance (S&MA). Under this system, each element (vehicle or ground system) defines a catalog of all expected debris based on launch history, component testing, or analysis. Debris transport analysis (DTA) is conducted using the debris catalog characteristics and potential flow transport mechanisms (e.g., vehicle aerodynamics, gravity, wind, plume-driven). The predicted debris impact locations and velocities are provided to the hardware owners, who use available test data and analysis to determine whether each component can withstand the impacts. In cases where the element hardware may be severely damaged or fail, the options are to mitigate the debris source through some change in design or operation, or to work with S&MA to try to characterize the probability of the impact and damage for program risk acceptance. Because of the differences in debris characteristics and transport, the DTA has been divided between the Liftoff and Ascent regimes. The development and application of Liftoff DTA methodology from the Shuttle Program to the current Artemis Program is the subject of this paper. Liftoff DTA covers the time from the start of pre-launch operations at the launch pad, up until the vehicle clears the launch tower and there is no longer any interaction with ground systems. Debris transport during this period is broadly classified as either gravity, wind, and plume-entrained (GWPE) or plume driven (PD). GWPE debris is generally lower speed, travelling in a forward-to-aft direction. PD transport includes flow features from the rocket ignition transient, as well as plume impingement and recirculation that occur as the vehicle lifts off the launch platform. In these cases, the debris typically moves in an aft-to-forward direction at higher speeds. The applicable transport mechanisms must be considered for each piece of debris depending on the material, and release location and time. For example, rust or metallic debris from the tower could fall (GWPE) and impact the vehicle before landing on the launch platform deck where it could be also be transported by plume impingement (PD). However, falling ice (GWPE) from an umbilical is unlikely to survive impact with the vehicle or launch platform and be available for PD transport. Modeling of debris transport is accomplished using a set of DTA tools which simulate debris trajectories subject to a reference frame acceleration (i.e., gravity) and aerodynamic drag. Where the trajectory encounters a solid surface, the debris is allowed to rebound with a specified coefficient of restitution. The drag is calculated by interpolating the fluid state at each point in the debris trajectory from high-fidelity computational fluid dynamics (CFD) simulations of the launch vehicle and pad. The CFD data may either be static (steady state or time averaged), typically for GWPE transport, or dynamic (time-accurate) for PD flow features like the ignition transient. Examples of the CFD flow field solutions for the Space Launch System (SLS) rocket and launch pad are shown in Figure 2. Typical SLS debris trajectory predictions from DTA are illustrated in Figure 3. The final version of this paper will include a more detailed examination of the Liftoff DTA process developed during the Shuttle Program, and how it has been augmented and applied to the SLS rocket under the Artemis Program. Comparisons with debris observations from the Artemis I launch will demonstrate validation of the tools and methodology.

Debris↗

ICME: The Design of Fit-for-Purpose Materials

With the increased emphasis on reducing the cost and time to market of new materials, ICME (Integrated Computational Materials Engineering) has become a fast growing discipline within materials science and engineering. The vision of ICME is compelling in many respects, not only for the value added in reducing time to market for new products with advanced, tailored materials, but also for enhanced efficiency and performance of these materials. Although the challenges and barriers (both technical and cultural) are formidable, substantial cost, schedule, and technical benefits can result from broad development, implementation, and validation of ICME principles. ICME is an integrated approach to the design of products, and the materials that comprise them, by linking material and structural models at multiple time and length scales. NASA’s Transformational Tools and Technology (TTT) Project sponsored a study (performed by a team led by Pratt & Whitney) in 2016 to define the potential 25-year future state required for integrated multiscale modeling of materials and systems (e.g., load-bearing structures) to accelerate the pace and reduce the expense of innovation in future aerospace and aeronautical systems. This talk will briefly review ICME, the findings of the 2040 Vision study, and discuss NASA’s TTT 2040 implementation activities: with special emphasis on recent accomplishments. The 2040 study, NASA CR 2018- 219771, envisions the development of a cyber-physical-social ecosystem comprised of experimentally verified and validated computational models, tools, and techniques, along with the associated digital tapestry, that impacts the entire supply chain to enable cost-effective, rapid, and revolutionary design of fit-for-purpose materials, components, and systems by enabling the engineer to not only “design-with-the” material but also concurrently “design-the” material.

multiscale analysis↗