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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

NASA NTRS1982

Grid generation by elliptic partial differential equations for a tri-element Augmentor-Wing airfoil

Two efforts to numerically simulate the flow about the Augmentor-Wing airfoil in the cruise configuration using the GRAPE elliptic partial differential equation grid generator algorithm are discussed. The Augmentor-Wing consists of a main airfoil with a slotted trailing edge for blowing and two smaller airfoils shrouding the blowing jet. The airfoil and the algorithm are described, and the application of GRAPE to an unsteady viscous flow simulation and a transonic full-potential approach is considered. The procedure involves dividing a complicated flow region into an arbitrary number of zones and ensuring continuity of grid lines, their slopes, and their point distributions across the zonal boundaries. The method for distributing the body-surface grid points is discussed.

Sorenson, R. L.↗

NASA NTRS2012

High-Performance Modeling and Simulation of Anchoring in Granular Media for NEO Applications

NASA is interested in designing a spacecraft capable of visiting a near-Earth object (NEO), performing experiments, and then returning safely. Certain periods of this mission would require the spacecraft to remain stationary relative to the NEO, in an environment characterized by very low gravity levels; such situations require an anchoring mechanism that is compact, easy to deploy, and upon mission completion, easy to remove. The design philosophy used in this task relies on the simulation capability of a high-performance multibody dynamics physics engine. On Earth, it is difficult to create low-gravity conditions, and testing in low-gravity environments, whether artificial or in space, can be costly and very difficult to achieve. Through simulation, the effect of gravity can be controlled with great accuracy, making it ideally suited to analyze the problem at hand. Using Chrono::Engine, a simulation pack age capable of utilizing massively parallel Graphic Processing Unit (GPU) hardware, several validation experiments were performed. Modeling of the regolith interaction has been carried out, after which the anchor penetration tests were performed and analyzed. The regolith was modeled by a granular medium composed of very large numbers of convex three-dimensional rigid bodies, subject to microgravity levels and interacting with each other with contact, friction, and cohesional forces. The multibody dynamics simulation approach used for simulating anchors penetrating a soil uses a differential variational inequality (DVI) methodology to solve the contact problem posed as a linear complementarity method (LCP). Implemented within a GPU processing environment, collision detection is greatly accelerated compared to traditional CPU (central processing unit)- based collision detection. Hence, systems of millions of particles interacting with complex dynamic systems can be efficiently analyzed, and design recommendations can be made in a much shorter time. The figure shows an example of this capability where the Brazil Nut problem is simulated: as the container full of granular material is vibrated, the large ball slowly moves upwards. This capability was expanded to account for anchors of different shapes and penetration velocities, interacting with granular soils.

Quadrelli, Marco B.↗

NASA NTRS2001

Testing the Preliminary X-33 Navigation System

The X-33 Reusable Launch Vehicle (RLV) must meet the demanding requirements of landing autonomously on a narrow landing strip following a flight that reaches an altitude of up to 200,000 feet and a speed in excess of Mach 9 with significant in-flight energy bleed-off maneuvers. To execute this flight regimen a highly reliable avionics system has been designed that includes three LN-100G Inertial Navigation System/Global Positioning System (INS/GPS) units as the primary navigation system for the X-33. NASA's Marshall Space Flight Center (MSFC) tested an INS/GPS system in real-time simulations to determine the ability of this navigation suite to meet the in flight and autonomous landing requirements of the X-33 RLV. A total of sixty-one open loop tests were performed to characterize the navigation accuracy of the LN-100G. Twenty-seven closed-loop tests were also performed to evaluate the performance of the X-33 Guidance, Navigation and Control (GN&C) algorithms with the real navigation hardware. These closed-loop tests were also designed to expose any integration or operational issues with the real-time X-33 vehicle simulation. Dynamic road tests of the INS/GPS were conducted by Litton to assess the performance of differential and nondifferential INS/GPS hybrid navigation solutions. The results of the simulations and road testing demonstrate that this novel solution is capable of meeting the demanding requirements of take-off, in-flight navigation, and autonomous landing of the X-33 RLV. This paper describes the test environment developed to stimulate the LN-100G and discusses the results of this test effort. This paper also presents recommendations for a navigation system suitable to an operational RLV system.

Lomas, James J.↗

NASA NTRS1972

Simulation study of the lift roll coupling problem for hovering VTOL aircraft

The effects of lift-roll coupling on the handling qualities of hovering VTOL aircraft using differential thrust for roll control were assessed in a piloted simulation study in the Ames six-degrees-of-freedom motion simulator. The configuration tested has three vertical thrust sources, one on the roll axis and two laterally displaced from the roll axis, with a thrust distribution of 25%/50%/25%. The outboard thrust sources were modulated to provide roll control whereas all three provided height control. Maximum thrust-to-weight ratio was varied together with a coupling parameter that combined roll-inertia, weight, and engine location. Results showed that handling qualities are affected not only by the occurrence of lift-roll coupling (dependent on both variables) but also by the severity of the coupling (dependent on the coupling parameter alone). However, the advantages of differential thrust for control can be retained with careful design.

Jacobsen, R. A.↗

NASA NTRS1993

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗

NASA NTRS1980

A computer simulation of the midlatitude plasmasphere and ionosphere

A computer model has been developed to simulate species density, temperature, and plasma flow in the ionosphere and plasmasphere. A new approach, the flux preserving method involving a Newton iteration, is used to solve a system of governing equations comprising four second-order partial differentials. Tests are performed to demonstrate that the simulation converges to a stable steady state solution, then steady state ion fluxes are analyzed. Finally, simulations of the collapse of the sunset ionosphere are presented. Comparisons with satellite and radar data show good agreement in a number of cases.

Young, E. R.↗

NASA NTRS1995

Coronal mass ejections and the evolution of the large-scale corona

Coronal mass ejections (CMEs) are dynamic, large-scale events in the solar corona that expel plasma and magnetic fields into the solar wind. The structure of the large-scale corona influences CME occurrence; for example, CMEs are most frequently observed to be disruptions of coronal or helmet streamer configurations on the Sun. We have investigated the evolution of the large-scale corona in the presence of differential rotation, using time-dependent magnetohydrodynamic (MHD) simulations in two and three dimensions. An equilibrium coronal configuration is first computed by superimposing a Parker solar wind flow on an initial potential field, the equations are integrated in time until a steady-state is reached. When differential rotation is introduced, we find that the configuration does not reach a steady state; the closed field regions (helmet streamers) disrupt recurrently, ejecting plasmoids into the solar wind. Our results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated. We will describe the evolution of the polarization brightness (pB) for the resulting configurations, as well as the propagation of the ejected plasmoids in the solar wind, and we will compare the effects of differential rotation with other possible initiation mechanisms, such as the emergence of new magnetic flux.

Linker, J. A.↗

NASA NTRS1982

The condensational instability in the solar transition region and corona

The stability of plasmas at temperatures and densities typical of the solar transition region and corona is investigated using both a linear analysis and nonlinear time-dependent numerical simulations. Growth rates, decay rates, and oscillation frequencies of the perturbations determined from the linear analysis are in excellent agreement with the simulations. The nonlinear regime is characterized by a bifurcation of the plasma into a cool dense condensation surrounded by a hot tenuous corona. The condensation may then be accelerated by forces in the plasma such as those arising from gravity or differential heating. Finally, the results of the detailed simulation show that the transition region is a dynamically stable structure which is the result of the nonlinear evolution of the condensational instability.

Oran, E. S.↗

NASA NTRS1985

Bridgman growth of semiconductors

The purpose of this study was to improve the understanding of the transport phenomena which occurs in the directional solidification of alloy semiconductors. In particular, emphasis was placed on the strong role of convection in the melt. Analytical solutions were not deemed possible for such an involved problem. Accordingly, a numerical model of the process was developed which simulated the transport. This translates into solving the partial differential equations of energy, mass, species, and momentum transfer subject to various boundary and initial conditions. A finite element method with simple elements was initially chosen. This simulation tool will enable the crystal grower to systematically identify and modify the important design factors within her control to produce better crystals.

Carlson, F. M.↗

NASA NTRS1993

Line-of-sight determination in real-time simulations

This paper describes the selection of a method for determining line-of-sight in real-time simulations for the NASA Ames Vertical Motion Simulator (VMS) facility. Five different combinations of terrain representation and line-of-sight determination algorithms were tested. A gridpost terrain format, in conjunction with a Digital Differential Analyzer algorithm, was found to best meet the simulation criteria of high speed, low storage requirements, and accuracy.

Kull, Frederick G., Jr.↗

NASA NTRS1999

Sensitive Procedures for Determining the Permeation Resistance of Chlorinated Polyethylene Against Liquid Propellants

The permeation resistance of chlorinated polyethylene (CPE) used in totally encapsulating chemical protective suits against the aerospace fuels hydrazine, monomethylhydrazine, and uns-dimethylhydrazine was determined by measuring the breakthrough time (BT) and time-averaged vapor transmission rate (VTR) using procedures consistent with ASTM F 739 and ASTM F 1383. Two exposure scenarios were simulated: a 2 hour (h) fuel vapor exposure, and a liquid fuel "splash" followed by a 2 h vapor exposure. To simulate internal suit pressure during operation, a positive differential pressure of 0.3 in. water (75 Pa) on the collection side of the permeation apparatus was used. Using the available data, a model was developed to estimate propellant concentrations inside an air-line fed, totally encapsulating chemical protective suit. Concentrations were calculated under simulated conditions of fixed vapor transmission rate, variable breathing air flow rate, and variable splash exposure area. Calculations showed that the maximum allowable permeation rates of hydrazine fuels through CPE were of the order of 0.05 to 0.08 ng/sq cm min for encapsulating suits with low breathing air flow rates (of the order of 5 scfm or 140 L min-1). Above these permeation rates, the 10 parts-per-billion (ppb) threshold limit value time-weighted average could be exceeded. To evaluate suit performance at 10 ppb threshold-limiting value/time-weighted average level concentrations, use of a sensitive analytical method such as cation exchange high performance liquid chromatography with amperometric detection was found to be essential. The analytical detection limit determines the lowest measurable VTR, which in turn governed the lowest per meant concentration that could be calculated inside the totally encapsulating chemical protective suit.

Waller, Jess M.↗

NASA NTRS2013

UV Lidar Receiver Analysis for Tropospheric Sensing of Ozone

A simulation of a ground based Ultra-Violet Differential Absorption Lidar (UV-DIAL) receiver system was performed under realistic daytime conditions to understand how range and lidar performance can be improved for a given UV pulse laser energy. Calculations were also performed for an aerosol channel transmitting at 3 W. The lidar receiver simulation studies were optimized for the purpose of tropospheric ozone measurements. The transmitted lidar UV measurements were from 285 to 295 nm and the aerosol channel was 527-nm. The calculations are based on atmospheric transmission given by the HITRAN database and the Modern Era Retrospective Analysis for Research and Applications (MERRA) meteorological data. The aerosol attenuation is estimated using both the BACKSCAT 4.0 code as well as data collected during the CALIPSO mission. The lidar performance is estimated for both diffuseirradiance free cases corresponding to nighttime operation as well as the daytime diffuse scattered radiation component based on previously reported experimental data. This analysis presets calculations of the UV-DIAL receiver ozone and aerosol measurement range as a function of sky irradiance, filter bandwidth and laser transmitted UV and 527-nm energy

Pliutau, Denis↗

NASA NTRS2018

Atmospheric Chemistry Modeling Using a Regression Forest Model

Atmospheric chemistry is central to many environmental issues such as air pollution, climate change, and stratospheric ozone loss. Chemistry Transport Models (CTM) are a central tool for understanding these issues, whether for research or for forecasting. 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 CTM. We have investigated a machine learning approach to solving the differential equations instead of solving them numerically. From an annual 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, centennial scale climate simulations etc. We discuss our approches' speed and accuracy, and highlight some potential future directions for improving this approach.

Keller, Christoph A.↗

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