Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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 541 records · Page 30

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.↗

The Design of a Remote Sensing Data Acquisition Campaign for Precision Agriculture and Some Early Results

In the 1970s NASA and the Department of Agriculture attempted to use the new Landsat MSS system for agricultural purposes. The program had relatively little success. With the advent of differential GPS, yield monitors on harvest equipment and higher spatial resolution remote sensing systems it seemed likely the situation should be reexamined. Therefore, a campaign of data acquisition involving remote sensing and other modalities with dependent research was assembled and funded by the Space Grant Consortia in Alabama and Georgia. The design of the remote sensing data acquisition was driven by the biology and physics of the crop system and limited by the available sensor platforms. Major parameters included crop stage, spatial resolution, seasonal and daily weather conditions, and which portion of the EM spectrum would actually capture the most discriminating information. Joint visible and Near IR with Thermal IR would permit use of existing indices, such as greenness, as well as phenomena driven by the plant' s evapotranspiration. Spatial resolution in the 2-5 meter range was chosen, avoiding many complexities caused by aliasing crop row spacing at, higher resolutions yet finer than the harvester's tightest recording rate. This dictates use of an airborne system. Use of an airborne system also makes scheduling around weather much simpler than use of satellite data. Active video calibration was recognized as essential if quantitative measures were ever to be obtained or reproduced. The system would also have to have onboard geoOF1 Based on these elements 3 data acquisitions have been flown. Seven flight lines were flown twice in 1998 and 16 lines flown in 1999. Total raw data is several GBytes. All of the data has now been geometrically corrected and some preliminary analysis accomplished. The thermal bands have an extremely high correlation with yield. For one@test case with corn, correlation in excess of 0.86 was obtained from a data acquisition two months prior to harvest! Soil images show significant within field variation in clay, soil brightness and emissivity. Light wind has been found to effect the reflectance and temperature of broad leaf crops, including soybeans, cotton and peanuts. Clearly, this work has already demonstrated some very important results. With continued development of the remote sensing technology there is good reason to believe this research will soon be able to help the individual farmer.

Rickmanl, D.↗

Simple Two-Dimensional Corrections for One-Dimensional Pulse Tube Models

One-dimensional oscillating flow models are very useful for designing pulse tubes. They are simple to use, not computationally intensive, and the physical relationship between temperature, pressure and mass flow are easy to understand when used in conjunction with phasor diagrams. They do not possess, however, the ability to directly calculate thermal and momentum diffusion in the direction transverse to the oscillating flow. To account for transverse effects, lumped parameter corrections, which are obtained though experiment, must be used. Or two-dimensional solutions of the differential fluid equations must be obtained. A linear two-dimensional solution to the fluid equations has been obtained. The solution provides lumped parameter corrections for one-dimensional models. The model accounts for heat transfer and shear flow between the gas and the tube. The complex Nusselt number and complex shear wall are useful in describing these corrections, with phase relations and amplitudes scaled with the Prandtl and Valensi numbers. The calculated ratio, a, between a two-dimensional solution of the oscillating temperature and velocity and a one-dimensional solution for the same shows a scales linearly with Va for Va less than 30. In this region alpha less than 0.5, that is, the enthalpy flow calculated with a two-dimensional model is 50% of a calculation using a one-dimensional model. For Va greater than 250, alpha = 0.8, showing that diffusion is still important even when it is confined to a thing layer near the tube wall.

Lee, J. M.↗

Predicting Lightning Initiation using Deep Learning

Lightning occurrence presents safety challenges to people and property. The main challenge with lightning safety is that the majority of guidance is reactive. In other words, lightning has to have already occurred nearby before a person will respond and take shelter. Further, most injuries or fatalities occur as the storm approaches, or as it's moving away, when rainfall may not be present at the time of the flash. Thus, this project develops a physically-based deep learning model to produce lightning probabilities out to 15 minutes. The deep learning model combines a Convolutional Neural Network (CNN) with a Long Short-Term Memory (LSTM) network to capture both the spatial and temporal evolution of storms to predict the probability that lightning initiation will occur in the next 15 minutes. The model combines radar reflectivity, correlation coefficient and differential reflectivity to inferred storm hydrometer type and precipitation phase, which aids in the identification of electrification processes. The model is trained with data from the Geostationary Lightning Mapper (GLM), which is a near infrared sensor onboard the GOES-R series of satellites that measures optical brightness from lightning. This presentation will provide an overview of the project.

Andrew T White↗

Global and Local Helioseismic Studies of Solar Convection Zone Dynamics Using SOI-MDI on SOHO

Our joint collaborative analyses of global mode data to characterize the solar differential rotation (e.g. Thompson et al. 1996, Schou et al. 1998), and most recently to detect and analyze temporal variations in angular velocity Omega profiles both within the convection zone and in the deeper radiative interior (e.g. Howe et al 2000a,b; Toomre et al. 2000), have led to a series of fascinating discoveries. These should be pursued further as the solar cycle continues. The physical deductions being made from these studies have been greatly strengthened by utilizing both SOI-MDI and GONG data in order to have two independent observational realizations of Doppler images spanning a five-year interval, using two separate procedures to determine global mode splittings, and then analyzing those splitting data sets using both RLS and SOLA inversion procedures. There are considerable subtleties in the effects of instrumental response functions and calibrations, sensitivity of peak finding algorithms and their mode leakage estimates, and stochastic variations in mode amplitudes that can all contribute to apparent changes in the Omega profiles being inferred from sequences of helioseismic data. We have come to understand the implications of many of these calibration and analysis steps, greatly aided by frequent multi-week collaborative working sessions in our Helioseismic Analysis Facility (HAF) at JILA involving many members of the SO1 dynamics and inversion team, including most of our Co-Is during the summer months when we hold intensive working sessions. Considerable further focused attention is required in a collaborative setting on such global mode issues as we continue studying the changing sun.

Toomre, Juri↗

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

New Boundary Constraints for Elliptic Systems used in Grid Generation Problems

This paper discusses new boundary constraints for elliptic partial differential equations as used in grid generation problems in generalized curvilinear coordinate systems. These constraints, based on the principle of local conservation of thermal energy in the vicinity of the boundaries, are derived using the Green's Theorem. They uniquely determine the so called decay parameters in the source terms of these elliptic systems. These constraints' are designed for boundary clustered grids where large gradients in physical quantities need to be resolved adequately. It is observed that the present formulation also works satisfactorily for mild clustering. Therefore, a closure for the decay parameter specification for elliptic grid generation problems has been provided resulting in a fully automated elliptic grid generation technique. Thus, there is no need for a parametric study of these decay parameters since the new constraints fix them uniquely. It is also shown that for Neumann type boundary conditions, these boundary constraints uniquely determine the solution to the internal elliptic problem thus eliminating the non-uniqueness of the solution of an internal Neumann boundary value grid generation problem.

Kaul, Upender K.↗

Restoration and Archiving of Data from the Plasma Composition Experiment on the International Sun-Earth Explorer One (ISEE 1)

The objective of this project has been to complete the archiving of energetic (10 eV/epsilon - 18 keV/epsilon) ion composition data from the Lockheed Plasma Composition Experiment on the International Sun-Earth Explorer One (ISEE 1) satellite, using a particular data format that had previously been approved by NASA and the NSSDC. That same format, a combination of ion velocity moments and differential flux spectra, had been used in 1991 to archive, at the NSSDC, the first 28 months (the "Prime" period of ISEE investigations) of data from the Lockheed instrument under NASA Contract NAS5-33047. With the completion of this project, the almost 4 1/2-year time span of these unique data is now covered by a very compact set, approximately 1 gigabyte in total, of electronic files with physical quantities, all in ASCII. The files are organized by data type and time of data acquisition, in Universal Time, and named according to year and day of year. Each calendar day has five separate files (five types of data), the lengths of which vary from day to day, depending on the instrument mode of operation. The data format and file structure are described in detail in appendices 1 and 2. The physical medium consists of high-density (6250 cpi) 9-track magnetic tapes, complemented by a set of hardcopy line plots of certain plasma parameters. In this case there are five tapes, to be added to the six previous ones from 1991, and 25 booklets of plots, one per month, to be added to the previous 28. The tapes, including an extra standard-density (1600 cpi) tape with electronic versions of the Data User's Guide and self-guiding VAX/VMS command files, and the hardcopy plots are being boxed for shipment to the NSSDC.

Lennartsson, O. W.↗

Investigation of Vibrational Control of the Bridgman Crystal Growth Technique

The objectives are: Conduct a parametric theoretical and numerical investigation of vibro-convective buoyancy-driven flow in differentially heated cylindrical containers. Investigate buoyant vibro-convective transport regimes in Bridgman-type systems with a focus on the use of vibration to suppress, or control, convection in order to achieve transport control during crystal growth. Assess the feasibility of vibro-convective control as a means of offsetting "g-jitter" effects under microgravity conditions, Exchange information with the experimental group at the General Physics Institute (GPI) of the Russian Academy of Science who are undertaking a complementary experimental program.

Fedoseyev, Alexandre I.↗

Two-Armed, Mobile, Sensate Research Robot

The Anthropomorphic Robotic Testbed (ART) is an experimental prototype of a partly anthropomorphic, humanoid-size, mobile robot. The basic ART design concept provides for a combination of two-armed coordination, tactility, stereoscopic vision, mobility with navigation and avoidance of obstacles, and natural-language communication, so that the ART could emulate humans in many activities. The ART could be developed into a variety of highly capable robotic assistants for general or specific applications. There is especially great potential for the development of ART-based robots as substitutes for live-in health-care aides for home-bound persons who are aged, infirm, or physically handicapped; these robots could greatly reduce the cost of home health care and extend the term of independent living. The ART is a fully autonomous and untethered system. It includes a mobile base on which is mounted an extensible torso topped by a head, shoulders, and two arms. All subsystems of the ART are powered by a rechargeable, removable battery pack. The mobile base is a differentially- driven, nonholonomic vehicle capable of a speed >1 m/s and can handle a payload >100 kg. The base can be controlled manually, in forward/backward and/or simultaneous rotational motion, by use of a joystick. Alternatively, the motion of the base can be controlled autonomously by an onboard navigational computer. By retraction or extension of the torso, the head height of the ART can be adjusted from 5 ft (1.5 m) to 6 1/2 ft (2 m), so that the arms can reach either the floor or high shelves, or some ceilings. The arms are symmetrical. Each arm (including the wrist) has a total of six rotary axes like those of the human shoulder, elbow, and wrist joints. The arms are actuated by electric motors in combination with brakes and gas-spring assists on the shoulder and elbow joints. The arms are operated under closed-loop digital control. A receptacle for an end effector is mounted on the tip of the wrist and contains a force-and-torque sensor that provides feedback for force (compliance) control of the arm. The end effector could be a tool or a robot hand, depending on the application.

Engelberger, J. F.↗

The physical properties of coronal holes

Coronal hole data obtained by Skylab during the period May 1973 to February 1974 are reviewed for the region from the photosphere out through the low corona. Coronal holes have a significantly lower density and temperature than the typical background corona, possess an open divergent magnetic field structure, and are the source of the high-speed solar wind stream which give rise to recurrent geomagnetic storms. The photospheric, chromospheric, and coronal signatures of coronal holes are described, and the evolution of coronal holes is examined, with attention to the lifetimes, the differential rotation, and the areas covered by the holes. The association of coronal holes to large-scale magnetic fields and to disk activity is considered.

Bohlin, J. D.↗

A physical basis for remote rock mapping of igneous rocks using spectral variations in thermal infrared emittance

Results of a theoretical investigation of the relation between spectral features in the 8-12 micrometer region and rock type are presented. Data on compositions of a suite of rocks and measurements of their spectral intensities in 8.2-10.9 and 9.4-12.1 micrometer bands published by Vincent (1973) were subjected to various quantitative procedures. There was no consistent direct relationship between rock group names and the relative spectral intensities. However, there is such a relationship between the Thornton-Tuttle (1960) Differentiation Index and the relative spectral intensities. This relationship is explicable on the basis of the change in average Si-O bond length which is a function of the degree of polymerization of the SiO4 tetrahedra of the silicate minerals in the igneous rocks.

Walter, L. S.↗

A differential CDM model for fatigue of unidirectional metal matrix composites

A multiaxial, isothermal, continuum damage mechanics (CDM) model for fatigue of a unidirectional metal matrix composite volume element is presented. The model is phenomenological, stress based, and assumes a single scalar internal damage variable, the evolution of which is anisotropic. The development of the fatigue damage model, (i.e., evolutionary law) is based on the definition of an initially transversely isotropic fatigue limit surface, a static fracture surface, and a normalized stress amplitude function. The anisotropy of these surfaces and function, and therefore the model, is defined through physically meaningful invariants reflecting the local stress and material orientation. This transversely isotropic model is shown, when taken to it's isotropic limit, to directly simplify to a previously developed and validated isotropic fatigue continuum damage model. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation in attempting to characterize a class of composite materials, and (2) the capability of the formulation in predicting anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Also, specific material parameters representing an initial characterization of the composite system SiC/Ti 15-3 and the matrix material (Ti 15-3) are reported.

Arnold, S. M.↗

Merging IceSAT GLAS and Terra MODIS Data in Order to Derive Forest Type Specific Tree Heights in the Central Siberian Boreal Forest

Mapping of boreal forest's type, biomass, and other structural parameters are critical for understanding of the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. We believe the nature of the forest structure information available from MISR and GLAS can be used to help identify forest type, age class, and estimate above ground biomass levels beyond that now possible with MODIS alone. The ground measurements will be used to develop relationships between remote sensing observables and forest characteristics and provide new information for understanding forest changes with respect to environmental change. Lidar is a laser altimeter that determines the distance from the instrument to the physical surface by measuring the time elapsed between the pulse emission and the reflected return. Other studies have shown that the returned signal may identify multiple returns originating from trees, building and other objects and permits the calculation of their height. Studies using field data have shown that lidar data can provide estimates of structural parameters such as biomass, stand volume and leaf area index and allows remarkable differentiation between primary and secondary forest. NASA's IceSAT Geoscience Laser Altimeter System (GLAS) was launched in January 2003 and collected data during February and September of that year. This study used data acquired over our study sites in central Siberia to examine the GLAS signal as a source of forest height and other structural characteristics. The purpose of our Siberia project is to improve forest cover maps and produce above-ground biomass maps of the boreal forest in Northern Eurasia from MODIS by incorporating structural information inherent in the Terra MISR and ICESAT Geoscience Laser Altimeter System (GLAS) instruments. A number of forest cover classifications exist for the boreal forest. We believe the limiting factor in these products is the lack of structural information, particularly in the vertical dimension. The emphasis of this project is to improve upon satellite maps of boreal forest structure parameters (i.e. height and biomass) using temporal, multi-angle, and vertical profile information of GLAS data. The existing and near future lidar data is useful for demonstrating these techniques and pursuing current estimates. Future lidar missions may be several years in the future, so we will work other new data sets that may aide in biomass estimates such as ALOS PALSAR We will continue this work to produce an accurate map of current above ground forest phytomass/carbon storage possible for the study area. We plan to develop, test, and integrate remote sensing methods for extracting forest canopy structure measures. We are compiling our field measurements and will compare them with the remote sensing methods where possible. We also be able to produce a realistic error bound on the remotely sensed carbon estimates.

Ranson, K. Jon↗

Imaging Asteroid 4 Vesta Using the Framing Camera

The Framing Camera (FC) onboard the Dawn spacecraft serves a dual purpose. Next to its central role as a prime science instrument it is also used for the complex navigation of the ion drive spacecraft. The CCD detector with 1024 by 1024 pixels provides the stability for a multiyear mission and its high requirements of photometric accuracy over the wavelength band from 400 to 1000 nm covered by 7 band-pass filters. Vesta will be observed from 3 orbit stages with image scales of 227, 63, and 17 m/px, respectively. The mapping of Vesta s surface with medium resolution will be only completed during the exit phase when the north pole will be illuminated. A detailed pointing strategy will cover the surface at least twice at similar phase angles to provide stereo views for reconstruction of the topography. During approach the phase function of Vesta was determined over a range of angles not accessible from earth. This is the first step in deriving the photometric function of the surface. Combining the topography based on stereo tie points with the photometry in an iterative procedure will disclose details of the surface morphology at considerably smaller scales than the pixel scale. The 7 color filters are well positioned to provide information on the spectral slope in the visible, the depth of the strong pyroxene absorption band, and their variability over the surface. Cross calibration with the VIR spectrometer that extends into the near IR will provide detailed maps of Vesta s surface mineralogy and physical properties. Georeferencing all these observation will result in a coherent and unique data set. During Dawn s approach and capture FC has already demonstrated its performance. The strong variation observed by the Hubble Space Telescope can now be correlated with surface units and features. We will report on results obtained from images taken during survey mode covering the whole illuminated surface. Vesta is a planet-like differentiated body, but its surface gravity and escape velocity are comparable to those of other asteroids and hence much smaller than those of the inner planets or

Keller, H. Uwe↗

Modeling Atmosphere-Ocean Radiative Transfer: A PACE Mission Perspective

The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth's sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates.

Chowdhary, Jacek↗

Mercury - Internal structure and thermal evolution

Mercury's thermal evolution and internal structure are modeled based on the planet's gross physical properties (which imply a high metallic iron content) and predictions for its chemistry made from the Lewis-Cameron model of condensation of the primitive solar nebula (which implies that Mercury may be composed only of those materials that condensed at temperatures near that of metallic iron condensation in the cooling nebula). Various heat sources, initial temperatures, and thermal conductivities are considered for a homogeneous model and a differentiated two-layer model. Density distributions are calculated from the mean density and estimates of the present-day temperature. The moment of inertia and the hydrostatic value of the second degree harmonic coefficient of Mercury's gravity field are found for the differentiated and undifferentiated models. These results should be useful for preliminary interpretation of the Mariner 10 measurements of Mercury's gravitational field.

Siegfried, R. W., II↗

A comparative study of theoretical methods on Goertler instability

Goertler vortices arise in boundary layers along concave surfaces due to centrifugal effects and these vortices in combination with Tollmien-Schlichting (TS) waves and crossflow vortices may play an important role in triggering early transition. There have been two distinct theoretical approaches to solve the problem. In the classical, normal-mode approach (NMA) the linear development of Goertler vortices is reduced to an eigenvalue problem. In the other method, initiated by Hall, the governing partial differential equations for the perturbations which are parabolic in the streamwise direction, are solved as an initial value problem. This method predicts multiple neutral curves depending on how and where the mean flow is perturbed, unlike the NMA where a unique but different neutral curve is predicted by each investigator. The present paper attempts to compare the two techniques and shows that mutually compatible results can be obtained when physically consistent assumptions are made.

Kalburgi, Vijay↗