Search NASA⌕ Search

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.

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

At least 307 records · Page 17

Communications-imposed pilot workload - A comparison of sixteen estimation techniques

Sixteen potential metrics of mental workload were investigated in regard to their relative sensitivity to communications load and their differential intrusion on primary task performance. A moving-base flight simulator was used to present three cross-country flights to each of 30 subject pilots, each flight varying only in the difficulty of the inherent communications requirements. With the exception of the rating scale measures, which were obtained immediately post-flight, all measures were taken over a seven minute segment of the flight task. The results indicated that both the Modified Cooper-Harper and the workload Multi-descriptor rating scales were reliably sensitive to changes in communications load. Also, the secondary task measure of time estimation and the physiological measure of pupil diameter yielded sensitivity. As expected, those primary task measures which were direct measures of communicative performance were also sensitive to load, while aircraft control primary task measures were not, attesting to the task-specificity of such measures. Finally, the intrusion analysis revealed no differential interference between workload measures.

Casali, J. G.↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Three-Dimensional Finite Element Ablative Thermal Response and Thermostructural Design of Thermal Protection Systems

A finite element ablation and thermal response program is presented for simulation of three-dimensional transient thermostructural analysis. The three-dimensional governing differential equations and finite element formulation are summarized. A novel probabilistic design methodology for thermal protection systems is presented. The design methodology is an eight step process beginning with a parameter sensitivity study and is followed by a deterministic analysis whereby an optimum design can determined. The design process concludes with a Monte Carlo simulation where the probabilities of exceeding design specifications are estimated. The design methodology is demonstrated by applying the methodology to the carbon phenolic compression pads of the Crew Exploration Vehicle. The maximum allowed values of bondline temperature and tensile stress are used as the design specifications in this study.

Dec, John A.↗

Differential Cratering of Synchronously Rotating Satellites by Ecliptic Comets

We use Monte Carlo methods to simulate impacts of ecliptic comets on the synchronously rotating satellites of giant planets. We reconfirm the long-standing prediction that the cratering rate should be much higher on the leading hemispheres than on the trailing hemisphere; indeed we find that previously published analytical formulations modestly underestimate the degree of apex-antapex asymmetry to be expected. We then compare our results to new mapping of impact craters on Ganymede, Callisto, and Triton. Ganymede reveals a pronounced apex-antapex asymmetry that is nonetheless much less than predicted. All of Triton's confirmed impact craters are clustered toward the apex of motion, far exceeding the predicted asymmetry. No asymmetry is observed on Callisto. In each case at least one of our basic assumptions must be wrong. Likely candidates include the following: (i) the surfaces of all but the most sparsely cratered satellites are saturated or nearly saturated with impact craters; (ii) these satellites have rotated nonsynchronously over geological time; (iii) most of the craters are made not by heliocentric (Sun-orbiting) comets and asteroids but rather by planetocentric (planet-orbiting) debris of indeterminate origin; or (iv) pathological endogenic resurfacing has created illusions of structure. Callisto's surface is readily classified as nearly saturated. Ganymede's bright terrains, although less heavily cratered than those of Callisto, can also be explained by crater densities approaching saturation on a world where endogenic processes were active. The leading alternative is nonsynchronous rotation, an explanation supported by the distribution of catenae (crater chains produced by impact of tidally disrupted comets). Triton's craters can be explained by planetocentric debris or by capricious resurfacing, but both hypotheses are inherently improbable.

Zahnle, Kevin↗

Effects of simulated microgravity on arterial nitric oxide synthase and nitrate and nitrite content

The aim of the present work was to investigate the alterations in nitric oxide synthase (NOS) expression and nitrate and nitrite (NOx) content of different arteries from simulated microgravity rats. Male Wistar rats were randomly assigned to either a control group or simulated microgravity group. For simulating microgravity, animals were subjected to hindlimb unweighting (HU) for 20 days. Different arterial tissues were removed for determination of NOS expression and NOx. Western blotting was used to measure endothelial NOS (eNOS) and inducible NOS (iNOS) protein content. Total concentrations of NOx, stable metabolites of nitric oxide, were determined by the chemiluminescence method. Compared with controls, isolated vessels from simulated microgravity rats showed a significant increase in both eNOS and iNOS expression in carotid arteries and thoracic aorta and a significant decrease in eNOS and iNOS expression of mesenteric arteries. The eNOS and iNOS content of cerebral arteries, as well as that of femoral arteries, showed no differences between the two groups. Concerning NOx, vessels from HU rats showed an increase in cerebral arteries, a decrease in mesenteric arteries, and no change in carotid artery, femoral artery and thoracic aorta. These data indicated that there were differential alterations in NOS expression and NOx of different arteries after hindlimb unweighting. We suggest that these changes might represent both localized adaptations to differential body fluid redistribution and other factors independent of hemodynamic shifts during simulated microgravity.

Non-NASA Center↗

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗

Analysis of Triangular Rollable and Collapsible Composite Booms Under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Analysis of Triangular Rollable and Collapsible (TRAC) Composite Booms under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Thermal Properties of Lunar Regolith Simulants

Various high temperature chemical processes have been developed to extract oxygen and metals from lunar regolith. These processes are tested using terrestrial analogues of the regolith. But all practical terrestrial analogs contain H2O and/or OH-, the presence of which has substantial impact on important system behaviors. We have undertaken studies of lunar regolith simulants to determine the limits of the simulants to validate key components for human survivability during sustained presence on the Moon. Differential Thermal Analysis (DTA) yields information on phase transitions and melting temperatures. Thermo-Gravimetric Analysis (TGA) with Fourier Transform Infrared (FTIR) analysis provides information on evolved gas species and their evolution temperature profiles. The DTA and TGA studies included JSC-1A fine (Johnson Space Center Mare Type 1A simulant), NU-LHT-2M (National Aeronautics and Space Administration (NASA)-- United States Geological Survey (USGS)--Lunar Highlands Type 2M simulant) and its proposed feedstocks: anorthosite; dunite; high quality (HQ) glass and the norite from which HQ glass is produced. As an example, the DTA and TGA profiles for anorthosite follow. The DTA indicates exothermic transitions at 355 and 490 C and endothermic transitions at 970 and 1235 C. Below the 355 C transition, water is lost accounting for approximately 0.1 percent mass loss. Just above 490 C a second type of water is lost, presumably bound in lattices of secondary minerals along with other volatile oxides. Limited TGA-FTIR data is available at the time of this writing. For JSC-1A fine, the TGA-FTIR indicates at least two kinds of water are evolved in the 100 to 500 and the 700 to 900 C ranges. Evolution of carbon dioxide types occurs in the 250 to 545, 545 to 705, and 705 to 985 C ranges. Geologically, the results are consistent with the evolution of "water" in its several forms, CO2 from break down of secondary carbonates and magmatic, dissolved gas and glass recrystallization

Street, Kenneth W., Jr.↗

Aeroelastic stability of coupled flap-lag-torsional motion of helicopter rotor blades in forward flight

A set of coupled periodic nonlinear differential flap-lag-torsional equations of blade motion is described which can simulate a blade having a precone, twist, distributed torsion, root torsion, torque offset, blade-root offsets, and offsets among the elastic axis, the aerodynamic center, and the blade cross-sectional center of mass. It is noted that the aerodynamic loads derived for these equations are applicable to cases of both hover and forward flight. The aeroelastic stability of a hingeless helicopter blade is investigated by linearizing the nonlinear differential equations of motion about a time-dependent equilibrium position of the helicopter in forward flight, using propulsion and moment trim procedures. A comparison of the results with those obtained previously for coupled flap-lag-torsional boundaries in hover indicates that aeroelastic-stability margins are degraded due to forward flight.

Friedmann, P.↗

Boundary value problems in magnetohydrodynamics (and fluid dynamics). I - Radiation boundary condition

In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.

Vanajakshi, T. C.↗

Experimental verification of dynamic simulation

The dynamics model here is a backhoe, which is a four degree of freedom manipulator from the dynamics standpoint. Two types of experiment are chosen that can also be simulated by a multibody dynamics simulation program. In the experiment, recorded were the configuration and force histories; that is, velocity and position, and force output and differential pressure change from the hydraulic cylinder, in the time domain. When the experimental force history is used as driving force in the simulation model, the forward dynamics simulation produces a corresponding configuration history. Then, the experimental configuration history is used in the inverse dynamics analysis to generate a corresponding force history. Therefore, two sets of configuration and force histories--one set from experiment, and the other from the simulation that is driven forward and backward with the experimental data--are compared in the time domain. More comparisons are made in regard to the effects of initial conditions, friction, and viscous damping.

Yae, K. Harold↗

MicroRNA Expression Profiles in Cultured Human Fibroblasts in Space

Microgravity, or an altered gravity environment from the static 1g, has been shown to influence global gene expression patterns and protein levels in living organisms. However, it is unclear how these changes in gene and protein expressions are related to each other or are related to other factors regulating such changes. A different class of RNA, the small non-coding microRNA (miRNA), can have a broad effect on gene expression networks by mainly inhibiting the translation process. Previously, we investigated changes in the expression of miRNA and related genes under simulated microgravity conditions on the ground using the NASA invented bioreactor. In comparison to static 1 g, simulated microgravity altered a number of miRNAs in human lymphoblastoid cells. Pathway analysis with the altered miRNAs and RNA expressions revealed differential involvement of cell communication and catalytic activity, as well as immune response signaling and NGF activation of NF-kB pathways under simulated microgravity condition. The network analysis also identified several projected networks with c- Rel, ETS1 and Ubiquitin C as key factors. In a flight experiment on the International Space Station (ISS), we will investigate the effects of actual spaceflight on miRNA expressions in nondividing human fibroblast cells in mostly G1 phase of the cell cycle. A fibroblast is a type of cell that synthesizes the extracellular matrix and collagen, the structural framework for tissues, and plays a critical role in wound healing and other functions. In addition to miRNA expressions, we will investigate the effects of spaceflight on the cellular response to DNA damages from bleomycin treatment.

Wu, Honglu↗

Evaluation of coated columbium alloy heat shields for space shuttle thermal protection system application. Volume 2, phase 2: Subsize heat shield and small size TPS evaluation

Initially a trade study was conducted of seven heat shield configurations. These were evaluated for structural reliability, fabricability, weight, inspectability, and refurbishability. Two concepts, a tee-stiffened and an open corrugation, were selected as offering the most potential for system success. Fourteen subsize heat shields of a full scale section were fabricated from C-129Y and Cb-752 and silicide coated with R-512E. These subsize panels were subjected to a simulated flight profile representing temperature, local surface pressures, and applied pressure differential loads. All corrugated panels of both alloys sustained 100 cycles without structural or coating failure. All Cb-752/R-512E panels performed well with one panel being successfully repaired after 66 cycles and completing 100 cycles. As a result of this evaluating the Cb-752/R-512E system was selected for hardware application during the subsequent phases. In addition, the tee-stiffened configuration was selected for further development and application in Phase III. This selection was based on an overall assessment of relative weight, cost, and structural performance of the tee-stiffened and open corrugation TPS.

Black, W. E.↗

Numerical solutions of atmospheric flow over semielliptical simulated hills

Atmospheric motion over obstacles on plane surfaces to compute simulated wind fields over terrain features was studied. Semielliptical, two dimensional geometry and numerical simulation of flow over rectangular geometries is also discussed. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale were solved by a finite difference technique. The mechanism of flow separation induced by a semiellipse is the same as flow over a gradually sloping surface for which the flow separation is caused by the interaction between the viscous force, the pressure force, and the turbulence level. For flow over bluff bodies, a downstream recirculation bubble is created which increases the aspect ratio and/or the turbulence level results in flow reattachment close behind the obstacle.

Shieh, C. F.↗

Vibrational relaxation of diatomic molecules in solids at low temperatures

The application of a hemiquantal method to the specific problem of the vibrational relaxation of a diatomic molecule embedded in a one dimensional lattice is presented. The vectorization of a CYBER 205 algorithm which integrates the 1,000 to 10,000 simultaneous hemiquantal differential equations is examined with comments on optimization. Results of the simulations are briefly discussed.

Halcomb, L. L.↗

The identification of a distributed parameter model for a flexible structure

A computational method is developed for the estimation of parameters in a distributed model for a flexible structure. The structure we consider (part of the RPL experiment) consists of a cantilevered beam with a thruster and linear accelerometer at the free end. The thruster is fed by a pressurized hose whose horizontal motion effects the transverse vibration of the beam. The Euler-Bernoulli theory is used to model the vibration of the beam and treat the hose thruster assembly as a lumped or point mass dashpot spring system at the tip. Measurements of linear acceleration at the tip are used to estimate the hose parameters (mass, stiffness, damping) and a Voigt-Kelvin viscoelastic structural damping parameter for the beam using a least squares fit to the data. Spline based approximations are considered to the hybrid (coupled ordinary and partial differential equations) systems; theoretical convergence results and numerical studies with both simulation and actual experimental data obtained from the structure are presented and discussed.

Banks, H. T.↗

Mathematical modeling of a class of multibody flexible spacecraft structures

A mathematical model for a general multibody flexible spacecraft is obtained. The generic spacecraft considered consists of a flexible central body to which a number of flexible multibody structures are attached. The coordinate systems used in the derivation allow effective decoupling of the translational motion of the entire spacecraft from its rotational motion about its center of mass. The derivation assumes that the deformations in the bodies are only due to elastic motions. The dynamic model derived is a closed-form vector-matrix differential equation. The model developed can be used for analysis and simulation of many realistic spacecraft configurations.

Kelkar, Atul, G.↗