SEARCH · Search NASA

Results for “differentiable simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

NASA NTRS1991

Large liquid rocket engine transient performance simulation system

A simulation system, ROCETS, was designed and developed to allow cost-effective computer predictions of liquid rocket engine transient performance. The system allows a user to generate a simulation of any rocket engine configuration using component modules stored in a library through high-level input commands. The system library currently contains 24 component modules, 57 sub-modules and maps, and 33 system routines and utilities. FORTRAN models from other sources can be operated in the system upon inclusion of interface information on comment cards. Operation of the simulation is simplified for the user by run, execution, and output processors. The simulation system makes available steady-state trim balance, transient operation, and linear partial generation. The system utilizes a modern equation solver for efficient operation of the simulations. Transient integration methods include integral and differential forms for the trapezoidal, first order Gear, and second order Gear corrector equations. A detailed technology test bed engine (TTBE) model was generated to be used as the acceptance test of the simulation system. The general level of model detail was that reflected in the Space Shuttle Main Engine DTM. The model successfully obtained steady-state balance in main stage operation and simulated throttle transients, including engine starts and shutdown. A NASA FORTRAN control model was obtained, ROCETS interface installed in comment cards, and operated with the TTBE model in closed-loop transient mode.

Mason, J. R.↗

NASA NTRS1995

Observational and numerical studies of extreme frontal scale contraction

The general objective of this effort is to increase understanding of how frontal scale contraction processes may create and sustain intense mesoscale precipitation along intensifying cold fronts. The five-part project (an expansion of the originally proposed two-part project) employed conventional meteorological data, special mesoscale data, remote sensing measurements, and various numerical models. First an idealized hydrostatic modeling study of the scale contraction effects of differential cloud cover on low-level frontal structure and dynamics was completed and published in a peer-reviewed journal. The second objective was to complete and publish the results from a three dimensional numerical model simulation of a cold front in which differential sensible heating related to cloud coverage patterns was apparently crucial in the formation of a severe frontal squall line. The third objective was to use a nonhydrostatic model to examine the nonlinear interactions between the transverse circulation arising from inhomogeneous cloud cover, the adiabatic frontal circulation related to semi-geostrophic forcing, and diabatic effects related to precipitation processes, in the development of a density current-like microstructure at the leading edge of cold fronts. Although the development of a frontal model that could be used to initialize such a primitive equation model was begun, we decided to focus our efforts instead on a project that could be successfully completed in this short time, due to the lack of prospects for continued NASA funding beyond this first year (our proposal was not accepted for future funding). Thus, a fourth task was added, which was to use the nonhydrostatic model to test tentative hypotheses developed from the most detailed observations ever obtained on a density current (primarily sodar and wind profiler data). These simulations were successfully completed, the findings were reported at a scientific conference, and the results have recently been submitted to a peer-reviewed journal. The fifth objective was to complete the analysis of data collected during the Cooperative Oklahoma Profiler Studies (COPS-91) field project, which was supported by NASA. The analysis of the mesoscale surface and sounding data, Doppler radar imagery, and other remote sensing data from multi frequency wind profiler, microwave radiometer, and the Radio Acoustic Sounding System has been completed. This study is a unique investigation of processes that caused the contraction of a cold front to a microscale zone exhibiting an undular bore-like structure. Results were reported at a scientific conference and are being prepared for publication. In summary, considerable progress has been achieved under NASA funding in furthering our understanding of frontal scale contraction and density current - gravity wave interaction processes, and in utilizing models and remotely sensed data in such studies.

Koch, Steven E.↗

NASA NTRS2019

Atmospheric Chemistry Modeling and Air Quality Forecasting Using Machine Learning

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O1D to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model.We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry.This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗

NASA NTRS2019

Atmospheric Chemistry Modeling Using Machine Learning

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O1D to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model. We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry. This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗

NASA NTRS2019

Machine Learning Application to Atmospheric Chemistry Modeling

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O (sup 1) D (Deuterium) to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model. We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry. This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗

NASA NTRS1987

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗

NASA NTRS1990

Performance evaluation of land mobile satellite system under vegetative shadowing using differential multiple TCM and QPSK

A comparative analysis by computer simulation of the land mobile satellite system (LMSS) is presented for the uncoded quadrature phase shift keying (QPSK), rate 2/3 8-PSK trellis coded modulation (TCM), and rate 4/6 multiplicity 2 8-PSK TCM schemes. An analytical model is used for simulating the fading channel. The simulation results show that under Rayleigh fades, the TCM designed for the fading channel is superior to uncoded QPSK and to the conventional TCM optimized for the additive white Gaussian noise (AWGN) channel. The performance of differentially detected TCM under AWGN and fading is examined based on the preliminary results. Results of the uncoded QPSK with fading channel modeled with the empirical fade data are given.

Kim, Junghwan↗

NASA NTRS1982

Validation of the NASCAP model using spaceflight data

The NASA Charging Analyzer Program (NASCAP) has been validated in a space environment. Data collected by the SCATHA (Spacecraft Charging at High Altitude) spacecraft has been used with NASCAP to simulate the charging response of the spacecraft ground conductor and dielectric surfaces with considerable success. Charging of the spacecraft ground observed in eclipse, during moderate and severe substorm environments, and in sunlight has been reproduced using the code. Close agreement between both the currents and potentials measured by the SSPM's, and the NASCAP simulated response, has been obtained for differential charging. It is concluded that NASCAP is able to predict spacecraft charging behavior in a space environment.

Stannard, P. R.↗

NASA NTRS1988

Trellis coding techniques for mobile communications

A criterion for designing optimum trellis codes to be used over fading channels is given. A technique is shown for reducing certain multiple trellis codes, optimally designed for the fading channel, to conventional (i.e., multiplicity one) trellis codes. The computational cutoff rate R0 is evaluated for MPSK transmitted over fading channels. Examples of trellis codes optimally designed for the Rayleigh fading channel are given and compared with respect to R0. Two types of modulation/demodulation techniques are considered, namely coherent (using pilot tone-aided carrier recovery) and differentially coherent with Doppler frequency correction. Simulation results are given for end-to-end performance of two trellis-coded systems.

Divsalar, D.↗

NASA NTRS1993

The Communication Link and Error ANalysis (CLEAN) simulator

During the period July 1, 1993 through December 30, 1993, significant developments to the Communication Link and Error ANalysis (CLEAN) simulator were completed and include: (1) Soft decision Viterbi decoding; (2) node synchronization for the Soft decision Viterbi decoder; (3) insertion/deletion error programs; (4) convolutional encoder; (5) programs to investigate new convolutional codes; (6) pseudo-noise sequence generator; (7) soft decision data generator; (8) RICE compression/decompression (integration of RICE code generated by Pen-Shu Yeh at Goddard Space Flight Center); (9) Markov Chain channel modeling; (10) percent complete indicator when a program is executed; (11) header documentation; and (12) help utility. The CLEAN simulation tool is now capable of simulating a very wide variety of satellite communication links including the TDRSS downlink with RFI. The RICE compression/decompression schemes allow studies to be performed on error effects on RICE decompressed data. The Markov Chain modeling programs allow channels with memory to be simulated. Memory results from filtering, forward error correction encoding/decoding, differential encoding/decoding, channel RFI, nonlinear transponders and from many other satellite system processes. Besides the development of the simulation, a study was performed to determine whether the PCI provides a performance improvement for the TDRSS downlink. There exist RFI with several duty cycles for the TDRSS downlink. We conclude that the PCI does not improve performance for any of these interferers except possibly one which occurs for the TDRS East. Therefore, the usefulness of the PCI is a function of the time spent transmitting data to the WSGT through the TDRS East transponder.

Ebel, William J.↗

NASA NTRS1995

The space-time solution element method: A new numerical approach for the Navier-Stokes equations

This paper is one of a series of papers describing the development of a new numerical method for the Navier-Stokes equations. Unlike conventional numerical methods, the current method concentrates on the discrete simulation of both the integral and differential forms of the Navier-Stokes equations. Conservation of mass, momentum, and energy in space-time is explicitly provided for through a rigorous enforcement of both the integral and differential forms of the governing conservation laws. Using local polynomial expansions to represent the discrete primitive variables on each cell, fluxes at cell interfaces are evaluated and balanced using exact functional expressions. No interpolation or flux limiters are required. Because of the generality of the current method, it applies equally to the steady and unsteady Navier-Stokes equations. In this paper, we generalize and extend the authors' 2-D, steady state implicit scheme. A general closure methodology is presented so that all terms up through a given order in the local expansions may be retained. The scheme is also extended to nonorthogonal Cartesian grids. Numerous flow fields are computed and results are compared with known solutions. The high accuracy of the scheme is demonstrated through its ability to accurately resolve developing boundary layers on coarse grids. Finally, we discuss applications of the current method to the unsteady Navier-Stokes equations.

Scott, James R.↗

NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗

NASA NTRS2001

Using EEG To Detect and Monitor Mental Fatigue

This project aims to develop EEG-based methods for detecting and monitoring mental fatigue. Mental fatigue poses a serious risk, even when performance is not apparently degraded. When such fatigue is associated with sustained performance of a single type of cognitive task it may be related to the metabolic energy required for sustained activation of cortical areas specialized for that task. The objective of this study was to adapt EEG to monitor cortical energy over a long period of performance of a cognitive task. Multielectrode event related potentials (ERPs) were collected every 15 minutes in nine subjects who performed a mental arithmetic task (algebraic sum of four randomly generated negative or positive digits). A new problem was presented on a computer screen 0.5 seconds after each response; some subjects endured for as long as three hours. ERPs were transformed to a quantitative measure of scalp electrical field energy. The average energy level at electrode P3 (near the left angular gyrus), 100-300 msec latency, was compared over the series of ERPs. For most subjects, scalp energy density at P3 gradually fell over the period of task performance and dramatically increased just before the subject was unable to continue the task. This neural response can be simulated for individual subjects using, a differential equation model in which it is assumed that the mental arithmetic task requires a commitment of metabolic energy that would otherwise be used for brain activities that are temporarily neglected. Their cumulative neglect eventually requires a reallocation of energy away from the mental arithmetic task.

Montgomery, Leslie↗

NASA NTRS2006

Practical Aspects of the Equation-Error Method for Aircraft Parameter Estimation

Various practical aspects of the equation-error approach to aircraft parameter estimation were examined. The analysis was based on simulated flight data from an F-16 nonlinear simulation, with realistic noise sequences added to the computed aircraft responses. This approach exposes issues related to the parameter estimation techniques and results, because the true parameter values are known for simulation data. The issues studied include differentiating noisy time series, maximum likelihood parameter estimation, biases in equation-error parameter estimates, accurate computation of estimated parameter error bounds, comparisons of equation-error parameter estimates with output-error parameter estimates, analyzing data from multiple maneuvers, data collinearity, and frequency-domain methods.

Morelli, Eugene a.↗

NASA NTRS1987

Time-partitioning simulation models for calculation on parallel computers

A technique allowing time-staggered solution of partial differential equations is presented in this report. Using this technique, called time-partitioning, simulation execution speedup is proportional to the number of processors used because all processors operate simultaneously, with each updating of the solution grid at a different time point. The technique is limited by neither the number of processors available nor by the dimension of the solution grid. Time-partitioning was used to obtain the flow pattern through a cascade of airfoils, modeled by the Euler partial differential equations. An execution speedup factor of 1.77 was achieved using a two processor Cray X-MP/24 computer.

Milner, Edward J.↗

NASA NTRS1987

Time-partitioning simulation models for calculation of parallel computers

Milner, Edward J.↗

NASA NTRS1973

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.↗