Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

An adaptive maneuvering logic computer program for the simulation of one-to-one air-to-air combat. Volume 2: Program description

A detailed description is presented of the computer programs in order to provide an understanding of the mathematical and geometrical relationships as implemented in the programs. The individual sbbroutines and their underlying mathematical relationships are described, and the required input data and the output provided by the program are explained. The relationship of the adaptive maneuvering logic program with the program to drive the differential maneuvering simulator is discussed.

Burgin, G. H.↗

Estimate of Orbiter static aeroelasticity properties via FLEXSTAB

The theory underlying the use of rigid-model experimental data in FLEXSTAB aeroelasticity calculations is summarized. A description is given of the aerodynamic and structural models used in the calculation of Space Shuttle Orbiter aeroelasticity properties during the re-entry phase. Experience gained with the use of pressure, nonlinear force and moment, and area-ratio correction methods in the Orbiter application is described and typical program output data are included. The inability of the program to correctly estimate differential control surface pressures and to accept external corrections for these terms is shown to be a major limitation in the Orbiter application.

Cavin, R. K., III↗

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗

Fourier Series Operating Package

This report presents a computer program for multiplying, adding, differentiating, integrating, "barring" and scalarly multiplying "literal" Fourier series as such, and for extracting the coefficients of specified terms.

Charnow, Milton L.↗

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↗

Mechanical signaling and the cellular response to extracellular matrix in angiogenesis and cardiovascular physiology

Great advances have been made in the identification of the soluble angiogenic factors, insoluble extracellular matrix (ECM) molecules, and receptor signaling pathways that mediate control of angiogenesis--the growth of blood capillaries. This review focuses on work that explores how endothelial cells integrate these chemical signals with mechanical cues from their local tissue microenvironment so as to produce functional capillary networks that exhibit specialized form as well as function. These studies have revealed that ECM governs whether an endothelial cell will switch between growth, differentiation, motility, or apoptosis programs in response to a soluble stimulus based on its ability to mechanically resist cell tractional forces and thereby produce cell and cytoskeletal distortion. Transmembrane integrin receptors play a key role in this mechanochemical transduction process because they both organize a cytoskeletal signaling complex within the focal adhesion and preferentially focus mechanical forces on this site. Molecular filaments within the internal cytoskeleton--microfilaments, microtubules, and intermediate filaments--also contribute to the cell's structural and functional response to mechanical stress through their role as discrete support elements within a tensegrity-stabilized cytoskeletal array. Importantly, a similar form of mechanical control also has been shown to be involved in the regulation of contractility in vascular smooth muscle cells and cardiac myocytes. Thus, the mechanism by which cells perform mechanochemical transduction and the implications of these findings for morphogenetic control are discussed in the wider context of vascular development and cardiovascular physiology.

Non-NASA Center↗

Trajectory Planning and Online Performance Model Estimation for Advanced Air Mobility

We propose a framework for adaptive guidance based on optimal control and estimation using closed-loop performance models. In the proposed approach, parameters of a reduced-order performance model are estimated in real time so that trajectories and guidance commands are replanned with more accurate knowledge of the system's response and performance capabilities. To apply this methodology to flight control systems, we introduce a simple yet expressive performance model modified from a reference linear design model. The proposed framework is applied to guidance of a simulated Advanced Air Mobility class concept aircraft experiencing control effector failures. We demonstrate that, by optimizing performance model parameters in real time, guidance commands can be intelligently adjusted to recover system stability and the performance of a full-order vehicle model, even in the event of effector failures. Furthermore, the reduced-order performance model requires a fraction of the computational cost of the full-order model, facilitating real-time use of adaptive, optimal guidance.

Differential dynamic programming↗

Trajectory Planning and Online Performance Model Estimation for Advanced Air Mobility

We propose a framework for adaptive guidance based on optimal control and estimation using closed-loop performance models. In the proposed approach, parameters of a reduced-order performance model are estimated in real time so that trajectories and guidance commands are replanned with more accurate knowledge of the system's response and performance capabilities. To apply this methodology to flight control systems, we introduce a simple yet expressive performance model modified from a reference linear design model. The proposed framework is applied to guidance of a simulated Advanced Air Mobility class concept aircraft experiencing control effector failures. We demonstrate that, by optimizing performance model parameters in real time, guidance commands can be intelligently adjusted to recover system stability and the performance of a full-order vehicle model, even in the event of effector failures. Furthermore, the reduced-order performance model requires a fraction of the computational cost of the full-order model, facilitating real-time use of adaptive, optimal guidance.

Differential dynamic programming↗

Program for the solution of multipoint boundary value problems of quasilinear differential equations

Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.

Source record↗

Cross Section Evaluation for Exclusive Channels of K+Λ and K+Σ0 Electroproduction off Protons Using CLAS Detector Data

In this work, a method for evaluating the cross sections of electroproduction of K+Λ0 and K+Σ0 off protons in the region of invariant masses of final hadrons MK + MY <W < 2.65 GeV (MK and MY being the masses of the kaon and hyperon, respectively) and squares of four-momentum transfers of virtual photons, i.e., photon virtualities 0 < Q2 < 5 GeV2, is developed based on experimental data of these exclusive channels’ cross sections measured by the CLAS detector in Hall B at Jefferson Lab. A set of algorithms has been implemented to evaluate the differential cross sections of these channels, along with their statistical and systematic uncertainties. A program was developed for the evaluation of differential cross sections and structure functions using C++ and Python libraries. An interactive website was created for working with the program, enabling the analysis of one-dimensional and two-dimensional dependences of structure functions and differential cross sections. The evaluation of differential cross sections for the K+Λ and K+Σ0 electroproduction channels is necessary for extracting the structure function σLT from the data on the polarization asymmetry of electroproduction reactions of these final states with longitudinally polarized electrons. The obtained results are also important for the development of realistic Monte Carlo event generators in planning future experiments and for evaluating the efficiency of detecting final particles when extracting reaction cross sections from experimental data.

Golda, A. V.↗

The design and performance of a modem for land mobile satellite communications

The digital implementation and performance of a multipath fade-tolerant differentially coherent 8-PSK-based modem for a land mobile satellite experiment (MSAT-X) program is described. At the receiver, differential detection using a matched filter arrangement and feedforward Doppler compensation is employed to achieve several operational benefits. For this L-band satellite link, the channel utilization is an efficient 1.0 bps/Hz for raw data and 1.5 bps/Hz for coded data. The experimental data analysis shows that the modem functions over the fading channel and performs within 1 dB of theoretical or simulation estimates.

Jedrey, T. C.↗

MGGHAT: Elliptic PDE software with adaptive refinement, multigrid and high order finite elements

MGGHAT (MultiGrid Galerkin Hierarchical Adaptive Triangles) is a program for the solution of linear second order elliptic partial differential equations in two dimensional polygonal domains. This program is now available for public use. It is a finite element method with linear, quadratic or cubic elements over triangles. The adaptive refinement via newest vertex bisection and the multigrid iteration are both based on a hierarchical basis formulation. Visualization is available at run time through an X Window display, and a posteriori through output files that can be used as GNUPLOT input. In this paper, we describe the methods used by MGGHAT, define the problem domain for which it is appropriate, illustrate use of the program, show numerical and graphical examples, and explain how to obtain the software.

Mitchell, William F.↗

Application of digital computer APU modeling techniques to control system design.

Study of the required controls for a H2-O2 auxiliary power unit (APU) technology program for the Space Shuttle. A steady-state system digital computer program was prepared and used to optimize initial system design. Analytical models of each system component were included. The program was used to solve a nineteen-dimensional problem, and then time-dependent differential equations were added to the computer program to simulate transient APU system and control. Some system parameters were considered quasi-steady-state, and others were treated as differential variables. The dynamic control analysis proceeded from initial ideal control modeling (which considered one control function and assumed the others to be ideal), stepwise through the system (adding control functions), until all of the control functions and their interactions were considered. In this way, the adequacy of the final control design over the required wide range of APU operating conditions was established.

Bailey, D. A.↗

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗