Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable simulations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Propellant-remaining modeling

A successful satellite mission is predicted upon the proper maintenance of the spacecraft's orbit and attitude. One requirement for planning and predicting the orbit and attitude is the accurate estimation of the propellant remaining onboard the spacecraft. Focuss is on the three methods that were developed for calculating the propellant budget: the errors associated with each method and the uncertainties in the variables required to determine the propellant remaining that contribute to these errors. Based on these findings, a strategy is developed for improved propellant-remaining estimation. The first method is based on Boyle's law, which related the values of pressure, volume, and temperature (PVT) of an ideal gas. The PVT method is used for the monopropellant and the bipropellant engines. The second method is based on the engine performance tests, which provide data that relate thrust and specific impulse associated with a propellant tank to that tank's pressure. Two curves representing thrust and specific impulse as functions of pressure are then generated using a polynomial fit on the engine performance data. The third method involves a computer simulation of the propellant system. The propellant flow is modeled by creating a conceptual model of the propulsion system configuration, taking into account such factors as the propellant and pressurant tank characteristics, thruster functionality, and piping layout. Finally, a thrust calibration technique is presented that uses differential correction with the computer simulation method of propellant-remaining modeling. Thrust calibration provides a better assessment of thruster performance and therefore enables a more accurate estimation of propellant consumed during a given maneuver.

Torgovitsky, S.↗
NASA NTRS1990

Simulation of the collapse and dissipation of Langmuir wave packets

Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.

Newman, D. L.↗
NASA NTRS2002

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗
NASA NTRS2016

Reserves and Trade Jointly Determine Exposure to Food Supply Shocks

While a growing proportion of global food consumption is obtained through international trade, there is an ongoing debate on whether this increased reliance on trade benefits or hinders food security, and specifically, the ability of global food systems to absorb shocks due to local or regional losses of production. This paper introduces a model that simulates the short-term response to a food supply shock originating in a single country, which is partly absorbed through decreases in domestic reserves and consumption, and partly transmitted through the adjustment of trade flows. By applying the model to publicly-available data for the cereals commodity group over a 17 year period, we find that differential outcomes of supply shocks simulated through this time period are driven not only by the intensification of trade, but as importantly by changes in the distribution of reserves. Our analysis also identifies countries where trade dependency may accentuate the risk of food shortages from foreign production shocks; such risk could be reduced by increasing domestic reserves or importing food from a diversity of suppliers that possess their own reserves. This simulation-based model provides a framework to study the short-term, nonlinear and out-of-equilibrium response of trade networks to supply shocks, and could be applied to specific scenarios of environmental or economic perturbations.

supply↗
DOE OSTI2024

Towards reverse mode automatic differentiation of Kokkos-based codes

Derivative computation is a key component of optimization, sensitivity analysis, uncertainty quantification, and the solving of nonlinear problems. Automatic differentiation (AD) is a powerful technique for evaluating such derivatives, and in recent years, has been integrated into programming environments such as Jax, PyTorch, and TensorFlow to support derivative computations needed for training of machine learning models, facilitating wide-spread use of these technologies. The C++ language has become the de facto standard for scientific computing due to numerous factors, yet language complexity has made the wide-spread adoption of AD technologies for C++ difficult, hampering the incorporation of powerful differentiable programming approaches into C++ scientific simulations. This is exacerbated by the increasing emergence of architectures, such as GPUs, with limited memory capabilities and requiring massive thread-level concurrency. C++ AD tools must effectively use these environments to bring novel scientific simulations to next-generation DOE experimental and observational facilities. In this project, we investigated source transformation-based automatic differentiation using LLVM compiler infrastructure to automatically generate portable and efficient gradient computations of Kokkos-based code. We have demonstrated that our proposed strategy is feasible by investigating the usage of a prototype LLVM-based source transformation tool to generate gradients of simple functions made of sequences of simple Kokkos parallel regions. Speedups of up to 500x compared to Sacado were observed on NVIDIA V100 GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Processing CCD Images to Detect Transits of Earth-Sized Planets: Maximizing Sensitivity While Achieving Reasonable Downlink Requirements

We have performed end-to-end laboratory and numerical simulations to demonstrate the capability of differential photometry under realistic operating conditions to detect transits of Earth-sized planets orbiting solar-like stars. Data acquisition and processing were conducted using the same methods planned for the proposed Kepler Mission. These included performing aperture photometry on large-format CCD images of an artificial star fields obtained without a shutter at a readout rate of 1 megapixel/sec, detecting and removing cosmic rays from individual exposures and making the necessary corrections for nonlinearity and shutterless operation in the absence of darks. We will discuss the image processing tasks performed `on-board' the simulated spacecraft, which yielded raw photometry and ancillary data used to monitor and correct for systematic effects, and the data processing and analysis tasks conducted to obtain lightcurves from the raw data and characterize the detectability of transits. The laboratory results are discussed along with the results of a numerical simulation carried out in parallel with the laboratory simulation. These two simulations demonstrate that a system-level differential photometric precision of 10-5 on five- hour intervals can be achieved under realistic conditions.

Earth-size planets↗
NASA NTRSDate not supplied

Differential Gene Expression in A Cross-Feeding Two-Species Model Microbial Community Under Simulated Microgravity and Deep-Space Radiation

A long-term goal of space biology is to understand interspecies microbial interactions in space. Presently, little is known about the combined effect of microgravity and ionizing radiation on bacterial community response when species are interdependent through exchange of metabolites in fluid medium (cross-feeding). Microgravity is expected to slow interspecies mass transfer and growth in cross-feeding communities in the low-shear, diffusion-limited environment, while ionizing radiation may influence stress response to direct (DNA damage) and indirect damage (ROS). Using a well-understood, two-species (Escherichia coli and Salmonella enterica) microbial community engineered to be a model for studying cross-feeding, we simulated galactic cosmic rays (GCRsim) and microgravity to test the hypothesis: exposure to ionizing radiation causes cell damage or stress, altering transcriptomic community responses in metabolically interdependent cells, which is exacerbated by microgravity. We expect to see differential gene expression between cross-feeding and non-cross-feeding communities. We measured GCRsim effects on growth and gene expression in well-mixed versus simulated-microgravity conditions and in cross-feeding and non-cross-feeding medium. Microbial cultures were inoculated into liquid medium in rotating wall vessels (RWV) with different rotation rates: 5 RPM (simulated microgravity) and 50 RPM (well-mixed). The E. coli-S. enterica consortium, under simulated microgravity, were exposed to 500 mGy of Simplified 5-ion Galactic Cosmic Ray Simulation for 2 hours at Brookhaven National Lab. We harvested samples 40 minutes after irradiation for extraction and sequencing (NASA GeneLab). Here we present the differential gene expression analysis results, which reveal altered transcriptomic community responses, even where growth rate differences are not observed. Gene expression of these actively metabolizing microbial communities in GCRsim may illuminate molecular mechanisms of microbial interactions in space. Understanding how microbial community gene expression, metabolism, and other cellular processes are influenced by spaceflight stressors can inform the use of microbes in human life support for low Earth orbit missions and beyond.

microgravity↗
NASA NTRS1984

Effects of suspension on tissue levels of glucocorticoid receptors

Differential muscle responses can be simulated by hypokinetic/hypodynamic (H/H) suspension of rats with complete unloading of the hindlimb muscles. Since mechanism(s) underlying these atrophic effects were not clearly elucidated, experiments were initiated to investigate a possible role for glucocorticoids in the physiological and biochemical responses to H/H. The principal objective was to assess the potential for alterations in peripheral responsiveness to glucocorticoids in response to H/H. Studies have initially focused on the determination of tissue levels of glucocorticoid receptors as one index of hormonal sensitivity at the cellular level. Four hindlimb muscles (soleus, gastrocnemius, plantaris and EDL), previously demonstrated to exhibit differential responses to H/H, were investigated. Receptor levels in other glucocorticoid sensitive tissues (heart, liver, and kidney) were determined. Male rats (180-200g) were suspended for 7 or 14 days, sacrificed by cervical dislocation, and the tissues excised.

Steffen, J. M.↗
NASA NTRS1987

Simulation of AFC for a DMSK Receiver

Low bit rates suitable for low-cost satellite communications. LOOP computer program written to simulate automatic frequency control (AFC) subsystem of differential minimum-shift keying (DMSK) receiver with bit rate of 2,400 baud. AFC simulated by LOOP is first-order loop configuration with first-order resistance-and-capacitance filter. Written in MS SuperSoft FORTRAN 77.

Davarian, Faramaz↗
DOE OSTI2025

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum- inspired techniques with classical methods. It focuses on fixed- point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed- point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum-inspired techniques with classical methods. It focuses on fixed-point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed-point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

X-33 Integrated Test facility, Extended Range Simulation

In support of the X-33 Single Stage To Orbit program, NASA Dryden Flight Research Center was selected to provide continuous communications coverage of the X-33 vehicle from launch, through landing at Malmstrom Air Force Base, Montana and Michaels Army Air Field, Utah. An extensive real-time range simulation capability is being developed to ensure successful communications with the autonomous X-33 vehicle. This paper will provide an overview of the various levels of simulation being developed to support the X-33 extended range subsystems. These subsystems include the Flight Termination System, L-Band command uplink subsystem and the S-Band telemetry downlink subsystem. In addition, the radar model developed provides continuous azimuth, elevation and range information based on the flight trajectory. The Dynamic Ground Station Analysis model developed by NASA Goddard Space Flight Center, calculate the received signal strength at each ground station. This model takes into consideration Radio Frequency (RF) link parameters such as frequency, antenna gain, space loss, plasma effects and the vehicle's position and attitude at any point in time during the flight path. All three RF links are then attenuated based on this calculated level and the RF signals are sent into telemetry receivers to emulate remote sites, or the power incident on the vehicle from uplinked signals. The best source received telemetry data is then passed back to the Launch and Mission Control Monitoring System (LMCMS) resident in the Operations Control Center. The LMCMS also provides the range simulation system the uplink command combined with differential GPS corrections. Later stages will require the progressive integration of actual range hardware with this simulation effort, leading to communication between telemetry, uplink and FTS antennas at NASA Dryden Flight Research Center, with vehicle antennas mounted on the Walter C. Williams Research Aircraft Integration Facility (RAIF). Decommutated Pulse Code Modulated (PCM) data is displayed on one of the four monitors that comprise the Range Safety Officer's (RSO) station. Also displayed are instantaneous impact prediction models, and Federal Aviation Administration (FAA) data for notification of other traffic in the area. Aside from initiating the flight termination command and validating communication links, the RSO station with the range simulation will be used to provide both range control and range safety officers training. The training is necessary to perform their respective functions with greater levels of confidence prior to first flight.

Sharma, Ashley↗
DOE OSTI2024

Implementing Ordinary Differential Equation Solvers in Rust Programming Language for Modeling Vehicle Powertrain Systems: Preprint

Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.

ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1979

Numerical calculations of high-altitude differential charging: Preliminary results

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

Laframboise, J. G.↗
DOE OSTI2025

Understanding the Role of Hydroxyl Functionalization in Linear Poly(Ethylenimine) for Oxidation‐Resistant Direct Air Capture of CO 2

Aminopolymer-based adsorbents are a prominent class of materials being used for direct air capture of CO 2 at the industrial scale. However, improving their working lifetime, specifically by increasing their resilience to oxidative degradation, remains an ongoing challenge. Toward this end, functionalization of aminopolymers with non-amine functionalities such as hydroxyls has emerged in recent years as a promising strategy toward improving adsorbent lifetime. Although there is a growing body of work demonstrating the effectiveness of this approach and investigating the origin of this improved stability, studies to date have primarily focused on branched aminopolymer systems such as branched poly(ethylenimine). In this work, hydroxyl-functionalized linear poly(ethylenimine) is used to continue to probe the underlying protective mechanism of this strategy. A combination of thermogravimetric analysis, NMR relaxometry, differential scanning calorimetry, and computational simulations is used to better understand the relationship between the extent of chemical functionalization, physical properties, and adsorbent performance.

amine-based adsorbent↗