Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differentiable programming”

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 145 records · Page 8

The 3DGRAPE book: Theory, users' manual, examples

A users' manual for a new three-dimensional grid generator called 3DGRAPE is presented. The program, written in FORTRAN, is capable of making zonal (blocked) computational grids in or about almost any shape. Grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. The smoothness for which elliptic methods are known is seen here, including smoothness across zonal boundaries. An introduction giving the history, motivation, capabilities, and philosophy of 3DGRAPE is presented first. Then follows a chapter on the program itself. The input is then described in detail. A chapter on reading the output and debugging follows. Three examples are then described, including sample input data and plots of output. Last is a chapter on the theoretical development of the method.

Sorenson, Reese L.↗

Optimized solution of Kepler's equation

A detailed description is presented of KEPLER, an IBM 360 computer program used for the solution of Kepler's equation for eccentric anomaly. The program KEPLER employs a second-order Newton-Raphson differential correction process, and it is faster than previously developed programs by an order of magnitude.

Kohout, J. M.↗

Recent advances in methods for numerical solution of O.D.E. initial value problems

In the mathematical modeling of physical systems, it is often necessary to solve an initial value problem (IVP), consisting of a system of ordinary differential equations (ODE). A typical program produces approximate solutions at certain mesh points. Almost all existing codes try to control the local truncation error, while the user is really interested in controlling the true or global error. The present investigation provides a review of recent advances regarding the solution of the IVP, giving particular attention to stiff systems. Stiff phenomena are customarily defined in terms of the eigenvalues of the Jacobian. There are, however, some difficulties connected with this approach. It is pointed out that an estimate of the Lipschitz constant proves to be a very practical way to determine the stiffness of a problem.

Bui, T. D.↗

The zinc finger transcription factor Gfi1, implicated in lymphomagenesis, is required for inner ear hair cell differentiation and survival

Gfi1 was first identified as causing interleukin 2-independent growth in T cells and lymphomagenesis in mice. Much work has shown that Gfi1 and Gfi1b, a second mouse homolog, play pivotal roles in blood cell lineage differentiation. However, neither Gfi1 nor Gfi1b has been implicated in nervous system development, even though their invertebrate homologues, senseless in Drosophila and pag-3 in C. elegans are expressed and required in the nervous system. We show that Gfi1 mRNA is expressed in many areas that give rise to neuronal cells during embryonic development in mouse, and that Gfi1 protein has a more restricted expression pattern. By E12.5 Gfi1 mRNA is expressed in both the CNS and PNS as well as in many sensory epithelia including the developing inner ear epithelia. At later developmental stages, Gfi1 expression in the ear is refined to the hair cells and neurons throughout the inner ear. Gfi1 protein is expressed in a more restricted pattern in specialized sensory cells of the PNS, including the eye, presumptive Merkel cells, the lung and hair cells of the inner ear. Gfi1 mutant mice display behavioral defects that are consistent with inner ear anomalies, as they are ataxic, circle, display head tilting behavior and do not respond to noise. They have a unique inner ear phenotype in that the vestibular and cochlear hair cells are differentially affected. Although Gfi1-deficient mice initially specify inner ear hair cells, these hair cells are disorganized in both the vestibule and cochlea. The outer hair cells of the cochlea are improperly innervated and express neuronal markers that are not normally expressed in these cells. Furthermore, Gfi1 mutant mice lose all cochlear hair cells just prior to and soon after birth through apoptosis. Finally, by five months of age there is also a dramatic reduction in the number of cochlear neurons. Hence, Gfi1 is expressed in the developing nervous system, is required for inner ear hair cell differentiation, and its loss causes programmed cell death.

NASA Discipline Developmental Biology↗

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗

FORTRAN 4 computer program for calculating critical speeds of rotating shafts

A FORTRAN 4 computer program, written for the IBM DCS 7094/7044 computer, that calculates the critical speeds of rotating shafts is described. The shaft may include bearings, couplings, extra masses (nonshaft mass), and disks for the gyroscopic effect. Shear deflection is also taken into account, and provision is made in the program for sections of the shaft that are tapered. The boundary conditions at the ends of the shaft can be fixed (deflection and slope equal to zero) or free (shear and moment equal to zero). The fixed end condition enables the program to calculate the natural frequencies of cantilever beams. Instead of using the lumped-parameter method, the program uses continuous integration of the differential equations of beam flexure across different shaft sections. The advantages of this method over the usual lumped-parameter method are less data preparation and better approximation of the distribution of the mass of the shaft. A main feature of the program is the nature of the output. The Calcomp plotter is used to produce a drawing of the shaft with superimposed deflection curves at the critical speeds, together with all pertinent information related to the shaft.

Trivisonno, R. J.↗

A study of the Sabatier-methanation reaction kinetics

The kinetics of the Sabatier methanation reaction, the reduction of carbon dioxide with hydrogen to methane and water, was investigated for 58 percent nickel on kieselguhr catalyst and 20 percent ruthenium on alumina catalyst. Differential rate data from an experimental program were correlated with a power function rate equation both for forward and reverse reactions. The kinetic parameters of activation energy, frequency rate constant and reaction order were determined for the rate equation. The values of these parameters were obtained from an Arrhenius plot of the experimental differential rate data. Also the carbon monoxide side reaction effect was measured and included in the correlation of parameters. The reaction was found to fit the rate equation experimentally within the temperature range 421 K, where the reaction effectively begins, the 800 K where the reaction rate drops and departs from the rate equation form.

Verostko, C. E.↗

Phototransistors Development and their Applications to Lidar

Custom-designed two-micron phototransistors have been developed using Liquid Phase Epitaxy (LPE), Molecular Beam Epitaxy (MBE) and Metal-Organic Chemical Vapor Deposition (MOCVD) techniques under Laser Risk Reduction Program (LRRP). The devices were characterized in the Detector Characterization Laboratory at NASA Langley Research Center. It appears that the performance of LPE- and MBE-grown phototransistors such as responsivity, noise-equivalent-power, and gain, are better than MOCVD-grown devices. Lidar tests have been conducted using LPE and MBE devices under the 2-micrometer CO2 Differential Absorption Lidar (DIAL) Instrument Incubator Program (IIP) at the National Center for Atmospheric Research (NCAR), Boulder, Colorado. The main focus of these tests was to examine the phototransistors performances as compared to commercial InGaAs avalanche photodiode by integrating them into the Raman-shifted Eye-safe Aerosol Lidar (REAL) operating at 1.543 micrometers. A simultaneous measurement of the atmospheric backscatter signals using the LPE phototransistors and the commercial APD demonstrated good agreement between these two devices. On the other hand, simultaneous detection of lidar backscatter signals using MBE-grown phototransistor and InGaAs APD, showed a general agreement between these two devices with a lower performance than LPE devices. These custom-built phototransistors were optimized for detection around 2-micrometer wavelength while the lidar tests were performed at 1.543 micrometers. Phototransistor operation at 2-micron will improve the performance of a lidar system operating at that wavelength. Measurements include detecting hard targets (Rocky Mountains), atmospheric structure consisting of cirrus clouds and boundary layer. These phototransistors may have potential for high sensitivity differential absorption lidar measurements of carbon dioxide and water vapor at 2.05-micrometers and 1.9-micrometers, respectively.

Abedin, M. N.↗

Surface characterization and testing II; Proceedings of the Meeting, San Diego, CA, Aug. 10, 11, 1989

Various papers on surface characterization and testing are presented. Individual topics addressed include: simple Hartmann test data interpretation, optimum configuration of the Offner null corrector, system for phase-shifting interferometry in the presence of vibration, fringe variation and visibility in speckle-shearing interferometry, functional integral representation of rough surfaces, calibration of surface heights in an interferometric optical profiler, image formation in common path differential profilometers, SEM of optical surfaces, measuring surface profiles with scanning tunneling microscopes, surface profile measurements of curved parts, high-resolution optical profiler, scanning heterodyne interferometer with immunity from microphonics, real-time crystal axis measurements of semiconductor materials, radial metrology with a panoramic annular lens, surface analysis for the characterization of defects in thin-film processes, Spacelab Optical Viewport glass assembly optical test program for the Starlab mission, scanning differential intensity and phase system for optical metrology.

Greivenkamp, John E.↗

Numerical calculations of high-altitude differential charging: Preliminary results

A two dimensional simulation program was constructed in order to obtain theoretical predictions of floating potential distributions on geostationary spacecraft. The geometry was infinite-cylindrical with angle dependence. Effects of finite spacecraft length on sheath potential profiles can be included in an approximate way. The program can treat either steady-state conditions or slowly time-varying situations, involving external time scales much larger than particle transit times. Approximate, locally dependent expressions were used to provide space charge, density profiles, but numerical orbit-following is used to calculate surface currents. Ambient velocity distributions were assumed to be isotropic, beam-like, or some superposition of these.

Laframboise, J. G.↗

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗

Advanced Weather Awareness and Reporting Enhancements

AWARE (Aviation Weather Awareness and Reporting Enhancements) was a NASA Cooperative Research and Development program conducted jointly by Rockwell Scientific, Rockwell Collins, and NASA. The effort culminated in an enhanced weather briefing and reporting tool prototype designed to integrate graphical and text-based aviation weather data to provide clear situational awareness in the context of a specific pilot, flight and equipment profile. The initial implementation of AWARE was as a web-based preflight planning tool, specifically for general aviation pilots, who do not have access to support such as the dispatchers available for commercial airlines. Initial usability tests showed that for VFR (Visual Flight Rules) pilots, AWARE provided faster and more effective weather evaluation. In a subsequent formal usability test for IFR (Instrument Flight Rules) pilots, all users finished the AWARE tests faster than the parallel DUAT tests, and all subjects graded AWARE higher for effectiveness, efficiency, and usability. The decision analysis basis of AWARE differentiates it from other aviation safety programs, providing analysis of context-sensitive data in a personalized graphical format to aid pilots/dispatchers in their complex flight requirements.

Busquets, Anthony M.↗

MACSYMA's symbolic ordinary differential equation solver

The MACSYMA's symbolic ordinary differential equation solver ODE2 is described. The code for this routine is delineated, which is of interest because it is written in top-level MACSYMA language, and may serve as a good example of programming in that language. Other symbolic ordinary differential equation solvers are mentioned.

Golden, J. P.↗

Modeling of the WSTF frictional heating apparatus in high pressure systems

In order to develop a computer program able to model the frictional heating of metals in high pressure oxygen or nitrogen a number of additions have been made to the frictional heating model originally developed for tests in low pressure helium. These additions include: (1) a physical property package for the gases to account for departures from the ideal gas state; (2) two methods for spatial discretization (finite differences with quadratic interpolation or orthogonal collocation on finite elements) which substantially reduce the computer time required to solve the transient heat balance; (3) more efficient programs for the integration of the ordinary differential equations resulting from the discretization of the partial differential equations; and (4) two methods for determining the best-fit parameters via minimization of the mean square error (either a direct search multivariable simplex method or a modified Levenburg-Marquardt algorithm). The resulting computer program has been shown to be accurate, efficient and robust for determining the heat flux or friction coefficient vs. time at the interface of the stationary and rotating samples.

Skowlund, Christopher T.↗

Studies of the differential absorption rocket experiment

Investigations of the ionosphere, in the rocket program of the Aeronomy Laboratory, include a propagation experiment, the data from which may be analyzed in several modes. This report considers in detail the differential absorption experiment. The sources of error and limitations of sensitivity are discussed. Methods of enhancing the performance of the experiment are described. Some changes have been made in the system and the improvement demonstrated. Suggestions are made for further development of the experiment.

Ginther, J. C.↗

Programming for energy monitoring/display system in multicolor lidar system research

The Z80 microprocessor based computer program that directs and controls the operation of the six channel energy monitoring/display system that is a part of the NASA Multipurpose Airborne Differential Absorption Lidar (DIAL) system is described. The program is written in the Z80 assembly language and is located on EPROM memories. All source and assembled listings of the main program, five subroutines, and two service routines along with flow charts and memory maps are included. A combinational block diagram shows the interfacing (including port addresses) between the six power sensors, displays, front panel controls, the main general purpose minicomputer, and this dedicated microcomputer system.

Alvarado, R. C., Jr.↗