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 289 records · Page 16

Ground resonance analysis using a substructure modeling approach

A convenient and versatile procedure for modeling and analyzing ground resonance phenomena is described and illustrated. A computer program is used which dynamically couples differential equations with nonlinear and time dependent coefficients. Each set of differential equations may represent a component such as a rotor, fuselage, landing gear, or a failed damper. Arbitrary combinations of such components may be formulated into a model of a system. When the coupled equations are formed, a procedure is executed which uses a Floquet analysis to determine the stability of the system. Illustrations of the use of the procedures along with the numerical examples are presented.

Chen, S.-Y.↗

Automated Surface Profilometer

Computer-controlled system saves time and labor. Electromechanical scanning system automatically measures deviation of sheets and plates from flatness. Quickly measures surface profile and detects bumps, bulges, and indentations as small as 0.001 in. in height or depth. X-y table moves specimen to programmed locations under digital linear variable differential transformer (LVDT) while computer records surface-height readings of LVDT. Resulting graphical display shows surface contour.

Chern, E. James↗

Prediction of pressure fluctuation in sounding rockets and manifolded recovery systems

The determination of altitude by means of barometric sensors in sounding rocket applications is discussed. A method for predicting the performance of such sensing systems is needed. A method is developed for predicting the pressure-time response of a volume subjected to subsonic air flow through from one to four passages. The pressure calculation is based on one-dimensional gas flow with friction. A computed program has been developed which solves the differential equations using a self-starting predictor-corrector integration technique. The input data required are the pressure sensing system dimensions, pressure forcing function(s) at the inlet port(s), and a trajectory over the time of analysis (altitude-velocity-time), if the forcing function is trajectory dependent. The program then computes the pressure-temperature history of the gas in the manifold over the time interval specified.

Laudadio, J. F.↗

THE YOUNG SOLAR ANALOGS PROJECT. I. SPECTROSCOPIC AND PHOTOMETRIC METHODS AND MULTI-YEAR TIMESCALE SPECTROSCOPIC RESULTS

This is the first in a series of papers presenting methods and results from the Young Solar Analogs Project, which began in 2007. This project monitors both spectroscopically and photometrically a set of 31 young (300–1500 Myr) solar-type stars with the goal of gaining insight into the space environment of the Earth during the period when life first appeared. From our spectroscopic observations we derive the Mount Wilson S chromospheric activity index (S{sub MW}), and describe the method we use to transform our instrumental indices to S{sub MW} without the need for a color term. We introduce three photospheric indices based on strong absorption features in the blue-violet spectrum—the G-band, the Ca i resonance line, and the Hydrogen-γ line—with the expectation that these indices might prove to be useful in detecting variations in the surface temperatures of active solar-type stars. We also describe our photometric program, and in particular our “Superstar technique” for differential photometry which, instead of relying on a handful of comparison stars, uses the photon flux in the entire star field in the CCD image to derive the program star magnitude. This enables photometric errors on the order of 0.005–0.007 magnitude. We present time series plots of our spectroscopic data for all four indices, and carry out extensive statistical tests on those time series demonstrating the reality of variations on timescales of years in all four indices. We also statistically test for and discover correlations and anti-correlations between the four indices. We discuss the physical basis of those correlations. As it turns out, the “photospheric” indices appear to be most strongly affected by emission in the Paschen continuum. We thus anticipate that these indices may prove to be useful proxies for monitoring emission in the ultraviolet Balmer continuum. Future papers in this series will discuss variability of the program stars on medium (days–months) and short (minutes to hours) timescales.

79 ASTRONOMY AND ASTROPHYSICS↗

The 2nd phase of the LEANDRE program: Water-vapor DIAL measurement

As a follow-on of the backscattered lidar, a differential absorption lidar (LEANDRE 2) is now being developed as part of the LEANDRE program for airborne meteorological studies. The primary measurement objective of LEANDRE 2 is water vapor. Pressure and temperature measurements are aimed at a second stage. The goals are to obtain a horizontal resolution of a few hundred meters for a vertical resolution of less than a hundred meters, with an absolute accuracy of 10 percent for humidity measurement. As compatibility is an important feature between the 2 first phases of LEANDRE, most of the LEANDRE 1 sub-system will be used and adapted for LEANDRE 2. For example, detection electronics, central computer, detectors and telescope will be the same. However, important modifications have to be done on the laser source, and spectral control has to be added. Most of the work is thus devoted to those developments, and the status is presented here.

Quaglia, P.↗

Random harmonic analysis program, L221 (TEV156). Volume 1: Engineering and usage

A digital computer program capable of calculating steady state solutions for linear second order differential equations due to sinusoidal forcing functions is described. The field of application of the program, the analysis of airplane response and loads due to continuous random air turbulence, is discussed. Optional capabilities including frequency dependent input matrices, feedback damping, gradual gust penetration, multiple excitation forcing functions, and a static elastic solution are described. Program usage and a description of the analysis used are presented.

Miller, R. D.↗

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗

Application of Aerosol Hygroscopicity Measured at the Atmospheric Radiation Measurement Program's Southern Great Plains Site to Examine Composition and Evolution

A Differential Mobility Analyzer/Tandem Differential Mobility Analyzer (DMA/TDMA) was used to measure submicron aerosol size distributions, hygroscopicity, and occasionally volatility during the May 2003 Aerosol Intensive Operational Period (IOP) at the Central Facility of the Atmospheric Radiation Measurement Program's Southern Great Plains (ARM SGP) site. Hygroscopic growth factor distributions for particles at eight dry diameters ranging from 0.012 micrometers to 0.600 micrometers were measured throughout the study. For a subset of particle sizes, more detailed measurements were occasionally made in which the relative humidity or temperature to which the aerosol was exposed was varied over a wide range. These measurements, in conjunction with backtrajectory clustering, were used to infer aerosol composition and to gain insight into the processes responsible for evolution. The hygroscopic growth of both the smallest and largest particles analyzed was typically less than that of particles with dry diameters of about 0.100 micrometers. It is speculated that condensation of secondary organic aerosol on nucleation mode particles is largely responsible for the minimal hygroscopic growth observed at the smallest sizes considered. Growth factor distributions of the largest particles characterized typically contained a nonhygroscopic mode believed to be composed primarily of dust. A model was developed to characterize the hygroscopic properties of particles within a size distribution mode through analysis of the fixed size hygroscopic growth measurements. The performance of this model was quantified through comparison of the measured fixed size hygroscopic growth factor distributions with those simulated through convolution of the size-resolved concentration contributed by each of the size modes and the mode-resolved hygroscopicity. This transformation from sizeresolved hygroscopicity to mode-resolved hygroscopicity facilitated examination of changes in the hygroscopic properties of particles within a size distribution mode that accompanied changes in the sizes of those particles. This model was used to examine three specific cases in which the sampled aerosol evolved slowly over a period of hours or days.

Gasparini, Roberto↗

Assembly and Delivery of Rabbit Capsules for Irradiation of Silicon Carbide Cladding Tube Specimens in the High Flux Isotope Reactor

Neutron irradiation of silicon carbide (SiC)-based fuel cladding under a high radial heat flux presents a critical challenge for SiC cladding concepts in light water reactors (LWRs). Fission heating in the fuel provides a high heat flux through the cladding, which, combined with the degraded thermal conductivity of SiC under irradiation, results in a large temperature gradient through the thickness of the cladding. The strong temperature dependence of swelling in SiC creates a complex stress profile in SiCbased cladding tubes as a result of differential swelling. The Nuclear Science User Facilities (NSUF) Program within the US Department of Energy Office of Nuclear Energy is supporting research efforts to improve the scientific understanding of the effects of irradiation on SiC cladding tubes. Ultimately, the results of this project will provide experimental validation of multi-physics models for SiC-based fuel cladding during LWR operation. The first objective of this project is to irradiate tube specimens using a previously developed design that allows for irradiation testing of miniature SiC tube specimens subjected to a high radial heat flux. The previous “rabbit” capsule design uses the gamma heating in the core of the High Flux Isotope Reactor (HFIR) to drive a high heat flux through the cladding tube specimens. A compressible aluminum foil allows for a constant thermal contact conductance between the cladding tubes and the rabbit housing despite swelling of the SiC tubes. To allow separation of the effects of irradiation from those due to differential swelling under a high heat flux, a new design was developed under the NSUF program. This design allows for irradiation of similar SiC cladding tube specimens without a high radial heat flux. This report briefly describes the irradiation experiment design concepts, summarizes the irradiation test matrix, and reports on the successful delivery of six rabbit capsules to the HFIR. Rabbits of both low and high heat flux configurations have been assembled, welded, evaluated, and delivered to the HFIR along with a complete quality assurance fabrication package. These rabbits contain a wide variety of specimens including monolith tubes, SiC fiber SiC matrix (SiC/SiC) composites, duplex specimens (inner composite, outer monolith), and specimens with a variety of metallic or ceramic coatings on the outer surface. The rabbits are targeted for insertion during HFIR cycle 475, which is scheduled for September 2017.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inelastic, nonlinear analysis of stiffened shells of revolution by numerical integration

This paper describes the latest addition to the STARS system of computer programs, STARS-2P, for the plastic, large deflection analysis of axisymmetrically loaded shells of revolution. The STARS system uses a numerical integration scheme to solve the governing differential equations. Several unique features for shell of revolution programs that are included in the STARS-2P program are described. These include orthotropic nonlinear kinematic hardening theory, a variety of shell wall cross sections and discrete ring stiffeners, cyclic and nonproportional mechanical and thermal loading capability, the coupled axisymmetric large deflection elasto-plastic torsion problem, an extensive restart option, arbitrary branching capability, and the provision for the inelastic treatment of smeared stiffeners, isogrid, and waffle wall constructions. To affirm the validity of the results, comparisons with available theoretical and experimental data are presented.

Levine, H. S.↗

Hes6 is required for actin cytoskeletal organization in differentiating C2C12 myoblasts

Hes6 is a member of the hairy-enhancer-of-split family of transcription factors that regulate proliferating cell fate in development and is known to be expressed in developing muscle. Here we investigate its function in myogenesis in vitro. We show that Hes6 is a direct transcriptional target of the myogenic transcription factors MyoD and Myf5, indicating that it is integral to the myogenic transcriptional program. The localization of Hes6 protein changes during differentiation, becoming predominantly nuclear. Knockdown of Hes6 mRNA levels by siRNA has no effect on cell cycle exit or induction of myosin heavy chain expression in differentiating C2C12 myoblasts, but F-actin filament formation is disrupted and both cell motility and myoblast fusion are reduced. The knockdown phenotype is rescued by expression of Hes6 cDNA resistant to siRNA. These results define a novel role for Hes6 in actin cytoskeletal dynamics in post mitotic myoblasts.

60 APPLIED LIFE SCIENCES↗

Field-Programmable Gate Array Computer in Structural Analysis: An Initial Exploration

This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.

Singleterry, Robert C., Jr.↗

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗

A proposed university/industry/government cooperative aero educational program

A program is proposed to provide mutual aid between the beleaguered Aero Engineering degree-granting departments in the United States and United States industry/government. It is proposed that any accrued cost differentials be supported by the latter groups, although all programs proposed will provide true mutual assistance. The program suggested is designed to provide assured industrial experience for Aero faculty members, combat engineering obsolescence for industry engineers, provide assured summer work for Aero undergraduates, encourage industry/government input in curriculum and course planning, and funnel significant laboratory devices into the university system.

Brodsky, R. F.↗

pyomo.dae: a modeling and automatic discretization framework for optimization with differential and algebraic equations

We describe pyomo.dae, an open source Python-based modeling framework that enables high-level abstract specification of optimization problems with differential and algebraic equations. The pyomo.dae framework is integrated with the Pyomo open source algebraic modeling language, and is available at http://www.pyomo.org. One key feature of pyomo.dae is that it does not restrict users to standard, predefined forms of differential equations, providing a high degree of modeling flexibility and the ability to express constraints that cannot be easily specified in other modeling frameworks. Other key features of pyomo.dae are the ability to specify optimization problems with high-order differential equations and partial differential equations, defined on restricted domain types, and the ability to automatically transform high-level abstract models into finite-dimensional algebraic problems that can be solved with off-the-shelf solvers. Moreover, pyomo.dae users can leverage existing capabilities of Pyomo to embed differential equation models within stochastic and integer programming models and mathematical programs with equilibrium constraint formulations. Collectively, these features enable the exploration of new modeling concepts, discretization schemes, and the benchmarking of state-of-the-art optimization solvers.

97 MATHEMATICS AND COMPUTING↗

Single-Crystal Sapphire Optical Fiber Sensor Instrumentation

This report summarizes technical progress on the program “Single-Crystal Sapphire Optical Fiber Sensor Instrumentation,” funded by the National Energy Technology Laboratory of the U.S. Department of Energy, and performed by the Center for Photonics Technology of the Bradley Department of Electrical and Computer Engineering at Virginia Tech. This project was completed in three phases, each with a separate focus. Phase I of the program, from October 1999 to April 2002, was devoted to development of sensing schema for use in high temperature, harsh environments. Different sensing designs were proposed and tested in the laboratory. Phase II of the program, from April 2002 to April 2009, focused on bringing the sensor technologies, which had already been successfully demonstrated in the laboratory, to a level where the sensors could be deployed in harsh industrial environments and eventually become commercially viable through a series of field tests. Also, a new sensing scheme was developed and tested with numerous advantages over all previous ones in Phase II. Phase III of the program, September 2009 to December 2013, focused on development of the new sensing scheme for field testing in conjunction with materials engineering of the improved sensor packaging lifetimes. In Phase I, three different sensing principles were studied: sapphire air-gap extrinsic Fabry-Perot sensors; intensity-based polarimetric sensors; and broadband polarimetric sensors. Black body radiation tests and corrosion tests were also performed in this phase. The outcome of the first phase of this program was the selection of broadband polarimetric differential interferometry (BPDI) for further prototype instrumentation development. This approach is based on the measurement of the optical path difference (OPD) between two orthogonally polarized light beams in a single-crystal sapphire disk. At the beginning of Phase II, in June 2004, the BPDI sensor was tested at the Wabash River coal gasifier facility in Terre Haute, Indiana. Due to business conditions at industrial partner and several logistical problems, this field test was not successful. An alternative high-temperature sensing system using sapphire wafer-based extrinsic Fabry-Perot interferometry was then developed as a significant improvement over the BPDI solution. From June 2006 to June 2008, three consecutive field tests were performed with the new sapphire wafer sensors at the TECO coal gasifier in Tampa, Florida. One of the sensors survived in the industrial coal gasifier for 7 months, over which time the existing thermocouples were replaced twice. The outcome of these TECO field tests suggests that the sapphire wafer sensor has very good potential to be commercialized. However packaging and sensor protection issues need additional development. During Phase III, several major improvements in the design and fabrication process of the sensor have been achieved through experiments and theoretical analysis. Studies on the property of the key components in the sensor head, including the sapphire fiber and sapphire wafer, were also conducted, for a better understanding of the sensor behavior. A final design based on all knowledge and experience has been developed, free of any issues encountered during the entire research. Sensors with this design performed well as expected in lab long-term tests, and were deployed in the sensing probe of the final coal-gasifier field test. Sensor packaging and protection was improved through materials engineering through testing of packaging designs in two blank probe packaging tests at Eastman Chemical in Kingsport, TN. Performance analysis of the blank probe packaging resulted in improve package designs culminating in a 3rd generation probe packaging utilized for the full field test of the sapphire optical sensor and materials designed sensor packaging.

47 OTHER INSTRUMENTATION↗