Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible 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 199 records · Page 11

Mechanical Backup For Fly-By-Wire Control System

Mechanical device eliminates need for redundant fly-by-wire subsystems. Main components are two linkages. One connected to control column in conventional, reversible control system. Other slides inside first linkage and connected to pilot's control wheel. In addition to aircraft applications, design used in control systems in which computer control desirable but safety backup systems required; for example, in boat rudders, engine controls in boats and automobiles, and controls in construction equipment.

Stewart, Eric C.↗

Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport

Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.

Intriligator, Devrie S.↗

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗

Flat Lens Focusing Demonstrated With Left-Handed Metamaterial

Left-handed metamaterials (LHM's) are a new media engineered to possess an effective negative index of refraction over a selected frequency range. This characteristic enables LHM's to exhibit physical properties never before observed. In particular, a negative index of refraction should cause electromagnetic radiation to refract or bend at a negative angle when entering an LHM, as shown in the figure above on the left. The figure on the right shows that this property could be used to bring radiation to a focus with a flat LHM lens. The advantage of a flat lens in comparison to a conventional curved lens is that the focal length could be varied simply by adjusting the distance between the lens and the electromagnetic wave source. In this in-house work, researchers at the NASA Glenn Research Center developed a computational model for LHM's with the three-dimensional electromagnetic commercial code Microwave Studio, constructed an LHM flat lens, and used it to experimentally demonstrate the reversed refraction and flat lens focusing of microwave radiation.

Wilson, Jeffrey D.↗

Coupled Groups of g-Modes in a Sun with Mixed Core

Groups of linear g-modes can sum to create long-lived nonlinear oscillations in small "hot volumes" very deep in the Sun that help drive the modes. In these volumes (dimensions -10 Mm), the time average rate of He-3 burning doubles as temperature fluctuations exceed 10% and rises by an order of magnitude for fluctuations of 25%. To be consistent with locally large motions, we impose a mixed shell on an otherwise standard solar model before computing g-mode solutions. Mixing in the assumed shell r = (0.10+/-0.03) R(sub sun) is rapid (<<10(exp 6) yr) with slower mixing somewhat beyond. If l is the principal spherical harmonic index, a set of g-modes for any single l less than or equal to 15 with five consecutive radial harmonics can be excited with nearly linear thermal amplitudes, A(sub T) less than or equal to 0.053, throughout the star and a fractional temperature fluctuation in its hot volume of (Delta)T/T less than or equal to 0.18. These thresholds for excitation will become smaller when sets for several values of l are computed simultaneously. There is some evidence for the rotation of g-mode sets in the long solar activity record and g-mode upward wave flux has been suggested to explain the 1.3 yr reversing flows tentatively detected below the Sun's convective envelope (CE). The large local amplitudes needed for excitation implies that g-modes may transport a non-negligible fraction of the solar luminosity, yet their near linear amplitudes outside the hot volume suggests amplitudes over most of the solar surface that would be barely detectable for l > 3. A formalism is presented for summing the g-modes and estimating growth rates under the approximation that modes are strictly linear except in a hot volume which holds only a few percent of mode kinetic energy. Finally over the range 2 less than or equal to l less than or equal to 30, we summed all zonal harmonics, m, for a given l and computed the relative angular orientations that would release the most nuclear energy. This should be close to the physically preferred angular state of such a family and a few examples were displayed.

Wolff, Charles L.↗

Radial rib antenna surface deviation analysis program

A digital computer program was developed which analyzes any radial rib antenna with ribs radiating from a central hub. The program has the capability for calculating the antenna surface contour (reversed pillowing effect), the optimum rib shape for minimizing the rms surface error, and the actual rms surface error. Rib deflection due to mesh tension and catenary cable tension can also be compensated for, and the pattern from which the mesh gores are cut can be determined.

Coyner, J. V., Jr.↗

Computation Of Facilitated Transport of O2 In Hemoglobin

Report describes computations of unsteady facilitated transport of oxygen through liquid membrane of hemoglobin. Used here, "facilitated transport" means diffusion of permeant through membrane in which that diffusion enhanced by reversible chemical reaction between permeant and membrane. In this case, reversible reactions between hemoglobin and oxygen.

Davis, Sanford↗

Efficient GPU-Accelerated MultiSource Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational-wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analyzing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains massive black hole binaries, compact galactic binaries, and a parametrized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble Markov Chain Monte Carlo (MCMC) sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or transdimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data.We recover posterior distributions for all 15 (6) of the injected massive black hole binaries (MBHB) in the LDC2A training (hidden) dataset. We catalog ∼12000 galactic binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA↗

Finite element methods in probabilistic mechanics

Probabilistic methods, synthesizing the power of finite element methods with second-order perturbation techniques, are formulated for linear and nonlinear problems. Random material, geometric properties and loads can be incorporated in these methods, in terms of their fundamental statistics. By construction, these methods are applicable when the scale of randomness is not too large and when the probabilistic density functions have decaying tails. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. Applications showing the effects of combined random fields and cyclic loading/stress reversal are studied and compared with Monte Carlo simulation results.

Liu, Wing Kam↗

General Aviation Task Force report

General aviation is officially defined as all aviation except scheduled airlines and the military. It is the only air transportation to many communities throughout the world. In order to reverse the recent decline in general aviation aircraft produced in the United States, the Task Force recommends that NASA provide the expertise and facilities such as wind tunnels and computer codes for aircraft design. General aviation manufacturers are receptive to NASA's innovations and technological leadership and are expected to be effective users of NASA-generated technologies.

Source record↗

Program Synthesizes UML Sequence Diagrams

A computer program called "Rational Sequence" generates Universal Modeling Language (UML) sequence diagrams of a target Java program running on a Java virtual machine (JVM). Rational Sequence thereby performs a reverse engineering function that aids in the design documentation of the target Java program. Whereas previously, the construction of sequence diagrams was a tedious manual process, Rational Sequence generates UML sequence diagrams automatically from the running Java code.

Barry, Matthew R.↗

Precision in Assembled Discrete Lattice Space Structures for Next-Generation ISAM Applications

Robust autonomous robotic assembly of large-scale space structures has been a long-term and challenging goal to enable higher quality in-space communication and science instrumentation. For optical observatory support structures and antenna structures, the challenge is the strict dimensional precision requirements (generally RMS surface error) driven by the operational radiation wavelength. Early theoretical work in the area linked RMS surface error of a reflector support truss plate to the natural passive dynamic mode frequencies of the structure, as well as to the error distribution in dimensions of constituent structural elements. Prior robotic truss assembly demonstrations focused on designing ultra-high precision structural elements, joints, and robotic actuators. Recently, an alternative space structure assembly strategy based on a programmable matter approach (NASA ARMADAS) was demonstrated. This approach uses lattice building blocks (voxels) that are reversibly, mechanically joined into a bulk lattice structure by robots that locomote in and on the structure itself. Such a system can achieve high-level autonomy with low computation, robustly assemble utilizing inexpensive and imprecise robots, and efficiently build structures several orders of magnitude larger than the assembly robots. However, for instrumentation support structure applications, the resulting precision of these building block-based lattice structures is not well studied. Since they are demonstrated with many more assembly units than prior art trusses, it is unclear whether the same precision design approaches apply. In this work, we study in simulation the effects of voxel geometry, error in distribution, assembly resolution (module size), and assembly geometry on the error of both beam and plate lattice structures. While average RMS error of a plate increases with plate size, increasing plate thickness quickly collapses RMS error towards a limit that is on the order of the error of the constituent parts. We validate our models against previously published precision measurements of built systems. Results from this study will guide manufacturing precision requirements, as well as designs, for future robotically assembled structural applications and establish feasibility for different applications.

Christine E Gregg↗

Unsteady Flow Simulation of High-speed Turbopumps

Computation of high-speed hydrodynamics requires high-fidelity simulation to resolve flow features involving transient flow, cavitation, tip vortex and multiple scales of unsteady fluctuations. One example of this type in aerospace is related to liquid-fueled rocket turbopump. Rocket turbopumps operate under severe conditions at very high rotational speeds typically at thousands of rpm. For example, the Shuttle orbiter low-pressure-fuel-turbopump creates transient flow features associated with reverse flows, tip clearance effects, secondary flows, vortex shedding, junction flows, and cavitation effects. Flow unsteadiness originating from the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the flow liners just upstream of the LPFTP. The reverse flow generated at the tip of the inducer blades travels upstream and interacts with the bellows cavity. Simulation procedure for this type high-speed hydrodynamic problems requires a method for quantifying multi-scale and multi-phase flow as well as an efficient high-end computing strategy. The current paper presents a high-fidelity computational procedure for unsteady hydrodynamic problems using a high-speed liquid-fueled rocket turbopump.

Kiris, Cetin C.↗

High-performance compression of astronomical images

Astronomical images have some rather unusual characteristics that make many existing image compression techniques either ineffective or inapplicable. A typical image consists of a nearly flat background sprinkled with point sources and occasional extended sources. The images are often noisy, so that lossless compression does not work very well; furthermore, the images are usually subjected to stringent quantitative analysis, so any lossy compression method must be proven not to discard useful information, but must instead discard only the noise. Finally, the images can be extremely large. For example, the Space Telescope Science Institute has digitized photographic plates covering the entire sky, generating 1500 images each having 14000 x 14000 16-bit pixels. Several astronomical groups are now constructing cameras with mosaics of large CCD's (each 2048 x 2048 or larger); these instruments will be used in projects that generate data at a rate exceeding 100 MBytes every 5 minutes for many years. An effective technique for image compression may be based on the H-transform (Fritze et al. 1977). The method that we have developed can be used for either lossless or lossy compression. The digitized sky survey images can be compressed by at least a factor of 10 with no noticeable losses in the astrometric and photometric properties of the compressed images. The method has been designed to be computationally efficient: compression or decompression of a 512 x 512 image requires only 4 seconds on a Sun SPARCstation 1. The algorithm uses only integer arithmetic, so it is completely reversible in its lossless mode, and it could easily be implemented in hardware for space applications.

White, Richard L.↗

Experience with k-epsilon turbulence models for heat transfer computations in rotating

This viewgraph presentation discusses geometry and flow configuration, effect of y+ on heat transfer computations, standard and extended k-epsilon turbulence model results with wall function, low-Re model results (the Lam-Bremhorst model without wall function), a criterion for flow reversal in a radially rotating square duct, and a summary.

Tekriwal, Prabbat↗

GRASP/Ada 95: Reverse Engineering Tools for Ada

The GRASP/Ada project (Graphical Representations of Algorithms, Structures, and Processes for Ada) has successfully created and prototyped an algorithmic level graphical representation for Ada software, the Control Structure Diagram (CSD), and a new visualization for a fine-grained complexity metric called the Complexity Profile Graph (CPG). By synchronizing the CSD and the CPG, the CSD view of control structure, nesting, and source code is directly linked to the corresponding visualization of statement level complexity in the CPG. GRASP has been integrated with GNAT, the GNU Ada 95 Translator to provide a comprehensive graphical user interface and development environment for Ada 95. The user may view, edit, print, and compile source code as a CSD with no discernible addition to storage or computational overhead. The primary impetus for creation of the CSD was to improve the comprehension efficiency of Ada software and, as a result, improve reliability and reduce costs. The emphasis has been on the automatic generation of the CSD from Ada 95 source code to support reverse engineering and maintenance. The CSD has the potential to replace traditional prettyprinted Ada source code. The current update has focused on the design and implementation of a new Motif compliant user interface, and a new CSD generator consisting of a tagger and renderer. The Complexity Profile Graph (CPG) is based on a set of functions that describes the context, content, and the scaling for complexity on a statement by statement basis. When combined graphicafly, the result is a composite profile of complexity for the program unit. Ongoing research includes the development and refinement of the associated functions, and the development of the CPG generator prototype. The current Version 5.0 prototype provides the capability for the user to generate CSDs and CPGs from Ada 95 source code in a reverse engineering as well as forward engineering mode with a level of flexibility suitable for practical application. This report provides an overview of the GRASP/Ada project with an emphasis on the current update.

Cross, James H., II↗

Automation of reverse engineering process in aircraft modeling and related optimization problems

During the year of 1994, the engineering problems in aircraft modeling were studied. The initial concern was to obtain a surface model with desirable geometric characteristics. Much of the effort during the first half of the year was to find an efficient way of solving a computationally difficult optimization model. Since the smoothing technique in the proposal 'Surface Modeling and Optimization Studies of Aerodynamic Configurations' requires solutions of a sequence of large-scale quadratic programming problems, it is important to design algorithms that can solve each quadratic program in a few interactions. This research led to three papers by Dr. W. Li, which were submitted to SIAM Journal on Optimization and Mathematical Programming. Two of these papers have been accepted for publication. Even though significant progress has been made during this phase of research and computation times was reduced from 30 min. to 2 min. for a sample problem, it was not good enough for on-line processing of digitized data points. After discussion with Dr. Robert E. Smith Jr., it was decided not to enforce shape constraints in order in order to simplify the model. As a consequence, P. Dierckx's nonparametric spline fitting approach was adopted, where one has only one control parameter for the fitting process - the error tolerance. At the same time the surface modeling software developed by Imageware was tested. Research indicated a substantially improved fitting of digitalized data points can be achieved if a proper parameterization of the spline surface is chosen. A winning strategy is to incorporate Dierckx's surface fitting with a natural parameterization for aircraft parts. The report consists of 4 chapters. Chapter 1 provides an overview of reverse engineering related to aircraft modeling and some preliminary findings of the effort in the second half of the year. Chapters 2-4 are the research results by Dr. W. Li on penalty functions and conjugate gradient methods for quadratic programming problems.

Li, W.↗

Gate-controlled-diodes in silicon-on-sapphire: A computer simulation

The computer simulation of the electrical behavior of a Gate-Controlled Diode (GCD) fabricated in Silicon-On-Sapphire (SOS) was described. A procedure for determining lifetime profiles from capacitance and reverse current measurements on the GCD was established. Chapter 1 discusses the SOS structure and points out the need of lifetime profiles to assist in device design for GCD's and bipolar transistors. Chapter 2 presents the one-dimensional analytical formula for electrostatic analysis of the SOS-GCD which are useful for data interpretation and setting boundary conditions on a simplified two-dimensional analysis. Chapter 3 gives the results of a two-dimensional analysis which treats the field as one-dimensional until the silicon film is depleted and the field penetrates the sapphire substrate. Chapter 4 describes a more complete two-dimensional model and gives results of programs implementing the model.

Gassaway, J. D.↗