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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Design of Digital Acquisition for Beam Current Monitor

As a part of the Proton Improvement Plan II (PIP-II) at Fermilab, instrumentation systems are being modernized to take advantage of the higher speeds and ease of use offered by standardized embedded systems like MicroTCA. A rear-transition module (RTM) is being designed to interface with said embedded systems. In each of the four identical channels on the RTM, the differential signal from an alternating-current current transformer (ACCT) transimpedance amplifier will again be amplified by a differential operation-amplifier, then filtered by a low-pass topology. The conditioned signal is then digitized at a maximum of 10MS/s by an analog to digital converter (ADC) integrated circuit. After digitization, the ADC passes the data to an off the shelf AdvancedMC (AMC) Xilinx FPGA module using low voltage differential signals. This paper will describe the simulation of analog circuitry for signal conditioning, simulation of digital signal integrity based on physical design as well as verification of design characteristics critical to signal integrity. This work aims to create a methodology that can be applied to future RTMs requiring application of high-speed digital design principles.

White, Turner [Fermilab] (ORCID:0009000348159256)↗

NASA NTRS1972

Simulation of a steady-state integrated human thermal system.

The mathematical model of an integrated human thermal system is formulated. The system consists of an external thermal regulation device on the human body. The purpose of the device (a network of cooling tubes held in contact with the surface of the skin) is to maintain the human body in a state of thermoneutrality. The device is controlled by varying the inlet coolant temperature and coolant mass flow rate. The differential equations of the model are approximated by a set of algebraic equations which result from the application of the explicit forward finite difference method to the differential equations. The integrated human thermal system is simulated for a variety of combinations of the inlet coolant temperature, coolant mass flow rate, and metabolic rates. Two specific cases are considered: (1) the external thermal regulation device is placed only on the head and (2) the devices are placed on the head and the torso. The results of the simulation indicate that when the human body is exposed to hot environment, thermoneutrality can be attained by localized cooling if the operating variables of the external regulation device(s) are properly controlled.

Hsu, F. T.↗

NASA NTRS1972

Simulation of a steady-state integrated human thermal system.

Hsu, F. T.↗

NASA NTRS1980

A global biogeocenotical biosphere simulation

This model of the D. Forrester type, constructed in differential equations, predicts the food and mineral supply for the factors biosphere population, depending on two socio-economic factors, until about the year 2500. If contemporary rates of natural resources utilization are maintained and there is no management of the environment, food resources will begin to limit human population growth after 2200, and mineral resources will after 2300. A decrease in the biosphere pollution, increase in effective agricultural production, and discovery of new energy sources may forestall or completely avert the onset of a crisis situation. Conservation measures, according to the model, are to a considerable extent realizable only if carried out simultaneously in both areas.

Moiseyev, N. N.↗

NASA NTRS1986

Partitioning and packing mathematical simulation models for calculation on parallel computers

The development of multiprocessor simulations from a serial set of ordinary differential equations describing a physical system is described. Degrees of parallelism (i.e., coupling between the equations) and their impact on parallel processing are discussed. The problem of identifying computational parallelism within sets of closely coupled equations that require the exchange of current values of variables is described. A technique is presented for identifying this parallelism and for partitioning the equations for parallel solution on a multiprocessor. An algorithm which packs the equations into a minimum number of processors is also described. The results of the packing algorithm when applied to a turbojet engine model are presented in terms of processor utilization.

Arpasi, D. J.↗

NASA NTRS1986

Mathematical model partitioning and packing for parallel computer calculation

This paper deals with the development of multiprocessor simulations from a serial set of ordinary differential equations describing a physical system. The identification of computational parallelism within the model equations is discussed. A technique is presented for identifying this parallelism and for partitioning the equations for parallel solution on a multiprocessor. Next, an algorithm which packs the equations into a minimum number of processors is described. The results of applying the packing algorithm to a turboshaft engine model are presented.

Arpasi, Dale J.↗

NASA NTRS1997

Plant and Animal Gravitational Biology

Session TA2 includes short reports covering: (1) The Interaction of Microgravity and Ethylene on Soybean Growth and Metabolism; (2) Structure and G-Sensitivity of Root Statocytes under Different Mass Acceleration; (3) Extracellular Production of Taxanes on Cell Surfaces in Simulated Microgravity and Hypergravity; (4) Current Problems of Space Cell Phytobiology; (5) Biological Consequences of Microgravity-Induced Alterations in Water Metabolism of Plant Cells; (6) Localization of Calcium Ions in Chlorella Cells Under Clinorotation; (7) Changes of Fatty Acids Content of Plant Cell Plasma Membranes under Altered Gravity; (8) Simulation of Gravity by Non-Symmetrical Vibrations and Ultrasound; and (9) Response to Simulated weightlessness of In Vitro Cultures of Differentiated Epithelial Follicular Cells from Thyroid.

Source record↗

NASA NTRS2000

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

NASA NTRS2001

Latitudinal Transport of Angular Momentum by Cellular Flows Observed with MDI

We have analyzed Doppler velocity images from the MDI instrument on SOHO to determine the latitudinal transport of angular momentum by the cellular photospheric flows. Doppler velocity images from 60-days in May to July of 1996 were processed to remove the p-mode oscillations, the convective blue shift, the axisymmetric flows, and any instrumental artifacts. The remaining cellular flows were examined for evidence of latitudinal angular momentum transport. Small cells show no evidence of any such transport. Cells the size of supergranules (30,000 km in diameter) show strong evidence for a poleward transport of angular momentum. This would be expected if supergranules are influenced by the Coriolis force, and if the cells are elongated in an east-west direction. We find good evidence for just such an east-west elongation of the supergranules. This elongation may be the result of differential rotation shearing the cellular structures. Data simulations of this effect support the conclusion that elongated supergranules transport angular momentum from the equator toward the poles, Cells somewhat larger than supergranules do not show evidence for this poleward transport. Further analysis of the data is planned to determine if the direction of angular momentum transport reverses for even larger cellular structures. The Sun's rapidly rotating equator must be maintained by such transport somewhere within the convection zone.

Hathaway, David H.↗

NASA NTRS2002

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

NASA NTRS2013

A Global Model of Meteoric Sodium

A global model of sodium in the mesosphere and lower thermosphere has been developed within the framework of the National Center for Atmospheric Research's Whole Atmosphere Community Climate Model (WACCM). The standard fully interactive WACCM chemistry module has been augmented with a chemistry scheme that includes nine neutral and ionized sodium species. Meteoric ablation provides the source of sodium in the model and is represented as a combination of a meteoroid input function (MIF) and a parameterized ablation model. The MIF provides the seasonally and latitudinally varying meteoric flux which is modeled taking into consideration the astronomical origins of sporadic meteors and considers variations in particle entry angle, velocity, mass, and the differential ablation of the chemical constituents. WACCM simulations show large variations in the sodium constituents over time scales from days to months. Seasonality of sodium constituents is strongly affected by variations in the MIF and transport via the mean meridional wind. In particular, the summer to winter hemisphere flow leads to the highest sodium species concentrations and loss rates occurring over the winter pole. In the Northern Hemisphere, this winter maximum can be dramatically affected by stratospheric sudden warmings. Simulations of the January 2009 major warming event show that it caused a short-term decrease in the sodium column over the polar cap that was followed by a factor of 3 increase in the following weeks. Overall, the modeled distribution of atomic sodium in WACCM agrees well with both ground-based and satellite observations. Given the strong sensitivity of the sodium layer to dynamical motions, reproducing its variability provides a stringent test of global models and should help to constrain key atmospheric variables in this poorly sampled region of the atmosphere.

meteoric flux↗

NASA NTRSDate not supplied

Uncertainty Quantification using Deep Ensembles for Decision Making in Cyber-Physical-Human Systems

In this paper and its companion, Differential Equation Approximation Using Gradient-Boosted Quantile Regression, Robison et al., we examine an approach to quantifying model uncertainty with the aim of increasing the trustworthiness of computational models in human-machine interactions. In Differential Equation Approximation Using Gradient-Boosted Quantile Regression, we focus on gradient-boosted decision trees, while in this one, we give more details about deep ensembles. Uncertainty quantification is crucial for building trustworthy autonomous decision-making agents in human-machine teams. There are two types of uncertainties: aleatoric and epistemic. The former is related to the inherent stochasticity (noise) of the process, whereas the latter is associated with the lack of knowledge or representation capability of models, such as neural networks. By lack of knowledge, we mean the model’s inability to accurately predict outputs for all possible inputs. The aleatory uncertainty can be estimated fairly easily with, for example, filters, whereas epistemic uncertainty is challenging to compute. This paper uses deep ensembles to quantify both aleatory and epistemic uncertainty. It can act as an uncertainty-aware surrogate transition model for decision-making frameworks. "Uncertainty-aware" means that the surrogate transition model should make predictions along with confidence in those predictions. In the context of decision-making, the transition models are ordinary differential equations (ODEs). Since ODEs can be simulated to make one-step or multi-step predictions, a good surrogate model for them should perform reasonably well in both modes. In a multi-step approach, the trajectory sampling method TS∞ was used to propagate uncertainty over multiple steps. The cartpole dynamical system was selected to demonstrate the ability of deep ensembles as good surrogate transition models for decision-making frameworks. The deep ensembles modeled the dynamics of cartpole ODEs and made uncertainty-aware predictions in single-step and multi-step transition modes.

CPH systems↗

DOE OSTI2025

Distinguishing fissile uranium isotopes using an active well neutron coincidence counter

Proposed thorium-based nuclear fuel cycles are likely to require quantification and verification of 233 U within nuclear material. Because of their similar fission cross sections, active neutron nondestructive assay (NDA) systems may respond similarly to 233 U and 235 U. Traditional safeguards equipment has been optimized for 235 U and 238 U quantification associated with conventional uranium/plutonium fuel cycles and may not be directly applicable to 233 U quantification when mixed with other actinides. This work used models of the large volume active well coincidence counter (LV-AWCC) at Oak Ridge National Laboratory to evaluate the performance of this neutron NDA system to differentiate fissile uranium isotopes. The models were developed to simulate NDA system performance in response to a number of triangular radiation signature training device sources within the central cavity or well. This work predicted that the LV-AWCC can effectively differentiate 233 U from 235 U in certain modes of operation. In active mode, the LV-AWCC with the cadmium liner results in different doubles count rates between the fissile isotopes for a given fissile uranium mass. Without the cadmium liner, the uranium isotopes provide a statistically indistinguishable doubles count rate response for the fissile masses considered in this work (up to approximately 150 g). The cadmium liner serves to harden the neutron interrogation spectrum, which better exploits the notable difference in the 233 U and 235 U fission cross sections at approximately 1 eV. In passive mode, the two fissile isotopes exhibit different doubles and singles count rates regardless of liner presence because the passive source strength of 233 U is approximately 2 orders of magnitude stronger than that of 235 U due to the shorter half-life and correspondingly higher (α, n) yield. We conclude that using neutron interrogation in the LV-AWCC, two measurements are needed to quantify 233 U content in mixed uranium items. The first measurement is used to determine the total fissile uranium mass using a mode that cannot distinguish fissile isotopes (i.e., where a similar response is observed for both fissile uranium isotopes such as active doubles without cadmium or using a thermal neutron interrogation source). In conclusion, the second measurement is used to determine the 233 U content by using a differentiating technique (e.g., passive doubles, passive doubles to singles ratio, active doubles with cadmium).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Kinetics of particles with short-range interactions

Self-assembly is one of the grand challenges of the 21st century – as the devices and materials we would like to build become too complex or small-scale for top-down manufacturing to be efficient, it is increasingly important to find ways to create these through bottom-up, dynamical approaches. Many particles used in self-assembly have very short-ranged attractive interactions, making simulations expensive or impossible. This proposal develops a set of conceptual and computational tools to study the dynamics of self-assembly for particles with short-ranged interactions, harnessing ideas in differential and computational geometry, and stochastic analysis, to accelerate simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2005

Inductive Learning Approaches for Improving Pilot Awareness of Aircraft Faults

Neural network flight controllers are able to accommodate a variety of aircraft control surface faults without detectable degradation of aircraft handling qualities. Under some faults, however, the effective flight envelope is reduced; this can lead to unexpected behavior if a pilot performs an action that exceeds the remaining control authority of the damaged aircraft. The goal of our work is to increase the pilot s situational awareness by informing him of the type of damage and resulting reduction in flight envelope. Our methodology integrates two inductive learning systems with novel visualization techniques. One learning system, the Inductive Monitoring System (IMS), learns to detect when a simulation includes faulty controls, while two others, Inductive Classification System (INCLASS) and multiple binary decision tree system (utilizing C4.5), determine the type of fault. In off-line training using only non-failure data, IMS constructs a characterization of nominal flight control performance based on control signals issued by the neural net flight controller. This characterization can be used to determine the degree of control augmentation required in the pitch, roll, and yaw command channels to counteract control surface failures. This derived information is typically sufficient to distinguish between the various control surface failures and is used to train both INCLASS and C4.5. Using data from failed control surface flight simulations, INCLASS and C4.5 independently discover and amplify features in IMS results that can be used to differentiate each distinct control surface failure situation. In real-time flight simulations, distinguishing features learned during training are used to classify control surface failures. Knowledge about the type of failure can be used by an additional automated system to alter its approach for planning tactical and strategic maneuvers. The knowledge can also be used directly to increase the pilot s situational awareness and inform manual maneuver decisions. Our multi-modal display of this information provides speech output to issue control surface failure warnings to a lesser-used communication channel and provides graphical displays with pilot-selectable !eve!s of details to issues additional information about the failure. We also describe a potential presentation for flight envelope reduction that can be viewed separately or integrated with an existing attitude indicator instrument. Preliminary results suggest that the inductive approach is capable of detecting that a control surface has failed and determining the type of fault. Furthermore, preliminary evaluations suggest that the interface discloses a concise summary of this information to the pilot.

Spikovska, Lilly↗

NASA NTRS1984

A study of the application of differential techniques to the global positioning system for a helicopter precision approach

The results of a simulation study to define the functional characteristics of a airborne and ground reference GPS receiver for use in a Differential GPS system are doumented. The operations of a variety of receiver types (sequential-single channel, continuous multi-channel, etc.) are evaluated for a typical civil helicopter mission scenario. The math model of each receiver type incorporated representative system errors including intentional degradation. The results include the discussion of the receiver relative performance, the spatial correlative properties of individual range error sources, and the navigation algorithm used to smooth the position data.

Mccall, D. L.↗

DOE OSTI2021

A Continuous-Discontinuous Galerkin Method for Electromagnetic Simulations Based on an All-Frequency Stable Formulation

In this paper, a potential-based partial-differential formulation, called the all-frequency stable formulation, is presented for the accurate and robust simulation of electromagnetic problems at all frequencies. Due to its stability from (near) dc to microwave frequencies, this formulation can be applied to simulate wide-band and multiscale problems without encountering the infamous low-frequency breakdown issue or applying basis function decompositions such as the tree-cotree splitting technique. To provide both efficient and flexible numerical solutions to the electromagnetic formulation, a mixed continuous-discontinuous Galerkin (CDG) method is proposed and implemented. In regions with homogeneous media, the continuous Galerkin method is employed to avoid the introduction of duplicated degrees of freedom (DoFs) on the elemental interfaces, while on the interfaces of two different media, the discontinuous Galerkin method is applied to permit the jump of the normal components of the electromagnetic fields. Numerical examples are provided to validate and demonstrate the proposed numerical solver for problems in a wide electromagnetic spectrum.

Yan, Su↗

DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗