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 379 records · Page 21

Aerodynamic and Aerothermal Simulations of Mars Concept Vehicles Using Overset DPLR

NASA's Descent System Studies (DSS) Program is studying various concept vehicles to enable landing of heavy payloads on the surface of Mars. While it is desirable to run high-fidelity CFD simulations to accurately assess the aerodynamic and aerothermal effects of various design changes during EDL, it is usually difficult to quickly generate high-quality grids suitable for such analyses. One approach to address this bottleneck in mesh generation is through the use oversetting grids. Although the overset approach is efficient and powerful in solving partial differential equations on complex geometries, new users often find it challenging to apply overset concepts for their simulations. For example, generating hyperbolic grids with sufficient overlap; priority in hole-cutting on multiple overlapping grids; and fixes to assemble overlapping viscous grids at the body surface. The objective of this presentation is to introduce a simple process that combines the advantages of near-body, point-matched, structured grids with oversetting background grids suitable for grid alignment. This approach allows for grids that can be sequenced, reclustering of mesh spacing at the wall, and grid alignment with the bow shock. The current methodology is tested on a Mid-L/D configuration using the overset DPLR code.

Overset Grids↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

The use of linked lists in the simulation of controller-structure interaction

An algorithm for the computer simulation of large space structures under active control is considered. Linked lists are used in a matrix data structure to implement the trapezoidal rule on the system differential equations. The use of the trapezoidal rule ensures that the numerical stability is equivalent to the system stability, which is essential for this type of simulation. The sparsity of the system matrices is exploited by the linked lists, and the algorithm efficiently steps through the lists in an orderly fashion. Results of simulations on a NASA large space structure experiment are reported.

Quan, Ralph↗

Heat transfer in the tip region of a rotor blade simulator

In gas turbines, the blades of axial turbine stages rotate in close proximity to a stationary peripheral wall. Differential expansion of the turbine wheel, blades, and the shroud causes variations in the size of the clearance gap between blade tip and stationary shroud. The necessity to tolerate this differential thermal expansion dictates that the clearance gap cannot be eliminated altogether, despite accurate engine machining. Pressure differences between the pressure and suction sides of a blade drives a flow through the clearance gap. This flow, the tip leakage flow, is detrimental to engine performance. The primary detrimental effect of tip leakage flow is the reduction of turbine stage efficiency, and a second is the convective heat transfer associated with the flow. The surface area at the blade tip in contact with the hot working gas represents an additional thermal loading on the blade which, together with heat transfer to the suction and pressure side surface area, must be removed by the blade internal cooling flows. Experimental results concerned with the local heat transfer characteristics on all surfaces of shrouded, rectangular cavities are reported. A brief discussion of the mass transfer system used is given.

Chyu, M. K.↗

Atmospheric Pre-Corrected Differential Absorption Techniques to Retrieve Columnar Water Vapor: Application to AVIRIS 91/95 Data

Water vapor is one of the main forces for weather development as well as for mesoscale air transport processes. The monitoring of water vapor is therefore an important aim in remote sensing of the atmosphere. Current operational systems for water vapor detection use primarily the emission in the thermal infrared (AVHRR, GOES, ATSR, Meteosat) or in the microwave radiation bands (DMSP). The disadvantage of current satellite systems is either a coarse spatial (horizontal) resolution ranging from one to tens of kilometers or a limited insight into the lower atmosphere. Imaging spectrometry on the other hand measures total column water vapor contents at a high spatial horizontal resolution and has therefore the potential of filling these gaps. The sensors of the AVIRIS instrument are capable of acquiring hyperspectral data in 224 bands located in the visible and near infrared at 10 nm resolution. This data includes the information on constituents of the earth's surface as well as of the atmosphere. The optical measurement of water vapor can be performed using sensor channels located in bands or lines of the absorption spectrum. The AVIRIS sensor has been used to retrieve water vapor and with less accuracy carbon dioxide, oxygen and ozone. To retrieve the water vapor amount, the so called differential absorption technique has been applied. The goal of this technique is to eliminate background factors by taking a ratio between channels within the absorption band and others besides the band. Various ratioing methods on the basis of different channels and calculation techniques were developed. The influence of a trace gas of interest on the radiance at the sensor level is usually simulated by using radiative transfer codes. In this study, the spectral transmittance and radiance are calculated by MODTRAN3 simulations with the new DISORT option. The objective of this work is to test the best performing differential absorption techniques for imaging spectrometry of tropospheric water vapor.

Schlaepfer, Daniel↗

Optimal Control via Self-Generated Stochasticity

The problem of global maxima of functionals has been examined. Mathematical roots of local maxima are the same as those for a much simpler problem of finding global maximum of a multi-dimensional function. The second problem is instability even if an optimal trajectory is found, there is no guarantee that it is stable. As a result, a fundamentally new approach is introduced to optimal control based upon two new ideas. The first idea is to represent the functional to be maximized as a limit of a probability density governed by the appropriately selected Liouville equation. Then, the corresponding ordinary differential equations (ODEs) become stochastic, and that sample of the solution that has the largest value will have the highest probability to appear in ODE simulation. The main advantages of the stochastic approach are that it is not sensitive to local maxima, the function to be maximized must be only integrable but not necessarily differentiable, and global equality and inequality constraints do not cause any significant obstacles. The second idea is to remove possible instability of the optimal solution by equipping the control system with a self-stabilizing device. The applications of the proposed methodology will optimize the performance of NASA spacecraft, as well as robot performance.

Zak, Michail↗

DYNGEN: A program for calculating steady-state and transient performance of turbojet and turbofan engines

The DYNGEN, a digital computer program for analyzing the steady state and transient performance of turbojet and turbofan engines, is described. The DYNGEN is based on earlier computer codes (SMOTE, GENENG, and GENENG 2) which are capable of calculating the steady state performance of turbojet and turbofan engines at design and off-design operating conditions. The DYNGEN has the combined capabilities of GENENG and GENENG 2 for calculating steady state performance; to these the further capability for calculating transient performance was added. The DYNGEN can be used to analyze one- and two-spool turbojet engines or two- and three-spool turbofan engines without modification to the basic program. A modified Euler method is used by DYNGEN to solve the differential equations which model the dynamics of the engine. This new method frees the programmer from having to minimize the number of equations which require iterative solution. As a result, some of the approximations normally used in transient engine simulations can be eliminated. This tends to produce better agreement when answers are compared with those from purely steady state simulations. The modified Euler method also permits the user to specify large time steps (about 0.10 sec) to be used in the solution of the differential equations. This saves computer execution time when long transients are run. Examples of the use of the program are included, and program results are compared with those from an existing hybrid-computer simulation of a two-spool turbofan.

Sellers, J. F.↗

Mathematical modeling and computer simulation of isoelectric focusing with electrochemically defined ampholytes

A mathematical model of isoelectric focusing at the steady state has been developed for an M-component system of electrochemically defined ampholytes. The model is formulated from fundamental principles describing the components' chemical equilibria, mass transfer resulting from diffusion and electromigration, and electroneutrality. The model consists of ordinary differential equations coupled with a system of algebraic equations. The model is implemented on a digital computer using FORTRAN-based simulation software. Computer simulation data are presented for several two-component systems showing the effects of varying the isoelectric points and dissociation constants of the constituents.

Palusinski, O. A.↗

Simulations and Observations of Heating of Auroral Ion Beams

In the auroral zone, quasi-static parallel electric fields produce beams of ionospheric ions (e.g., H(+), He(+) and 0(+)), which flow outward into the magnetosphere, providing a significant source of ions for the ring current and plasma sheet. Because the velocities to which these beams are accelerated is dependent on the mass of the ions, differential flows between the various ion species can develop which are unstable to an ion-ion streaming instability. Particle simulations and observations from DE 1 are used to investigate the heating of the ion beams produced by this instability. It is shown that there is net transfer of energy from the light ions to the heavy ions, with the heavy ions reaching maximum velocities near the beam velocity of the light ions. Bulk heating of the heavy ions occurs when their relative density is low while high-energy tails are produced when their relative density is high. The heating is primarily parallel to the magnetic field if the difference in the heavy and light ion beam velocities is subsonic while both perpendicular and parallel heating can occur if it is supersonic. In the latter case, very strong heating of an intermediate ion's species such as He(+) can also occur. Comparison with observations shows features consistent with heating via the ion-ion instability including perpendicular heating in the supersonic regime and parallel heating in the subsonic regime and a change in the heating between these regimes as the ratio of the H(+) beam speed to the local sound speed is observed to decrease. This heating is, however, not always observed in association with enhanced wave emissions. This lack of waves is attributed to reabsorption of the waves as the ions become heated.

Winglee, R. M.↗

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges II: Numerical Simulations and Comparisons with Experiments

We present a model that quantifies the magnitude of the ion-acoustic turbulence (IAT) in the plume of hollow cathodes and its effect on the resistivity and ion heating. The model takes the form of a partial differential equation (PDE) that can be solved concurrently with the equations of motion for a partially ionized plasma already included in our numerical code for the simulation of the plasma discharge in hollow cathodes, OrCa2D. We also determine that self-induced magnetic fields are not negligible in hollow cathodes operating at large discharge currents and implement in our code Ampere’s law and modifications to Ohm’s law that account for this effect. Numerical simulations that employed these models show large improvements in our agreement with experimental measurements with respect to a previous model, which assumed complete saturation of the IAT and did not account for the growth stage of the waves. In particular, the model is able to accurately predict the location and magnitude of the maximum resitivity to the electron current along the cathode centerline.

Lopez Ortega, Alejandro↗

Design study of shaft face seal with self-acting lift augmentation. 5: Performance in simulated gas turbine engine operation

The feasibility and the noncontact operation of the self-acting seal was demonstrated over a range of simulated gas turbine engine conditions from 200 to 500 ft/sec sliding speed. Sealed pressure differentials were 50 to 300 psi and sealed temperatures were 150 to 1200 F. Low leakage (about 1/10 that of conventional labyrinth seals) was exhibited in two endurance runs (200 and 338 hr) at 400 ft/sec, 200 psi and 1000 F (gas temperature). For these endurance runs, the self-acting pad wear was less than 3.8 micrometers (0.00015 in.); this low wear was attributed to the noncontact operation of the primary seal. Operating problems identified were fretting wear of the secondary seal and erosion of the primary seal by hard particles.

Ludwig, L. P.↗

Active control of potential of the geosynchronous satellites ATS-5 and ATS-6

The ATS-6 geosynchronous satellite carries two cesium ion thrusters as part of the technology experiment package. Their effectiveness in controlling the spacecraft potential was shown during normal operation of the thrusters. Auroral particle detectors were used to estimate spacecraft potential, amount of cesium propellant backflow to the spacecraft, and electron signature of the thruster. Operation of the thruster clamps the spacecraft to about -10 v in the presence of a wide range of ambient particle flux. Some positive ions from the thruster are returned to the spacecraft during this period. The plasma bridge neutralizer of the thruster is also effective in controlling the spacecraft potential. Laboratory simulation tests of the ATS-5 ion thruster thermionic emitter neutralizer showed that negative differential charging on the spacecraft surface probably prevents the emitter from completely discharging the spacecraft in flight.

Goldstein, R.↗

Simulation study of an automatic trim system for reducing the control forces on a light twin after an engine failure

An automatic trim system for reducing the control forces after an engine failure on a light twin has been investigated on the Langley General Aviation Simulator. The system schedules open-loop trim tab deflections as a function of differential propeller slipstream dynamic pressure and freestream dynamic pressure. The system is described and the airplane-system static and dynamic characteristics are documented. Three NASA research pilots evaluated the effectiveness of the system for takeoff and landing maneuvers. A variety of off-nominal system characteristics were studied. The system was judged to be generally beneficial, providing a 2 to 3 point improvement in pilot rating for the tasks used in the evaluations.

Stewart, E. C.↗

Chemical contamination remote sensing

A ground mobile laser test bed system was assembled to assess the feasibility of detection of various types of chemical contamination using Differential Scattering (DISC) and Differential Absorption (DIAL) Lidar techniques. Field experiments with the test bed system using chemical simulants were performed. Topographic reflection and range resolved DIAL detection of vapors as well as DISC detection of aerosols and surface contamination were achieved. Review of detection principles, design of the test bed system, and results of the experiments are discussed.

Carrico, J. P.↗

Quasi-geostrophic free mode models of long-lived Jovian eddies: Forcing mechanisms and crucial observational tests

Observations of Jupiter and Saturn long-lived eddies, such as Jupiter's Great Red Spot and White Ovals, are presently compared with laboratory experiments and corresponding numerical simulations for free thermal convection in a rotating fluid that is subject to horizontal differential heating and cooling. Difficulties in determining the essential processes maintaining and dissipating stable eddies, on the basis of global energy budget studies, are discussed; such difficulties do not arise in considerations of the flow's potential vorticity budget. On Jupiter, diabatically forced and transient eddy-driven flows primarily differ in the implied role of transient eddies in transporting potential vorticity across closed geostrophic streamlines in the time mean.

Read, P. L.↗

Dynamics of emerging active region flux loops

The buoyant rise of a magnetic flux loop arising from a single perturbed segment of a toroidal flux ring lying slightly beneath the base of the convection zone is studied by way of numerical simulations. We have considered flux loop evolution assuming both solid-body rotation, and differential rotation consistent with recent results from helioseismology. Our major results are presented, and we offer some speculations on the decay of active regions, based on the results of our studies. We speculate that as plasma in the tube attempts to establish hydrostatic equilibrium along the field lines after the flux emergence has taken place, the tube field strength at some intermediate depths below the surface becomes sufficiently small at the surface portions of the tube (which have cooled and undergone convective collapse) become dynamically disconnected from those portions near the base of the convection zone. The surface proportions of the emerged flux tubes are then transported by motions near the photosphere, such as supergranular convection and meridional flow.

Fan, Y.↗

Shear-induced opening of the coronal magnetic field

This work describes the evolution of a model solar corona in response to motions of the footpoints of its magnetic field. The mathematics involved is semianalytic, with the only numerical solution being that of an ordinary differential equation. This approach, while lacking the flexibility and physical details of full MHD simulations, allows for very rapid computation along with complete and rigorous exploration of the model's implications. We find that the model coronal field bulges upward, at first slowly and then more dramatically, in response to footpoint displacements. The energy in the field rises monotonically from that of the initial potential state, and the field configuration and energy appraoch asymptotically that of a fully open field. Concurrently, electric currents develop and concentrate into a current sheet as the limiting case of the open field is approached. Examination of the equations shows rigorously that in the asymptotic limit of the fully open field, the current layer becomes a true ideal MHD singularity.

Wolfson, Richard↗

Satellite Angular Rate Estimation From Vector Measurements

This paper presents an algorithm for estimating the angular rate vector of a satellite which is based on the time derivatives of vector measurements expressed in a reference and body coordinate. The computed derivatives are fed into a spacial Kalman filter which yields an estimate of the spacecraft angular velocity. The filter, named Extended Interlaced Kalman Filter (EIKF), is an extension of the Kalman filter which, although being linear, estimates the state of a nonlinear dynamic system. It consists of two or three parallel Kalman filters whose individual estimates are fed to one another and are considered as known inputs by the other parallel filter(s). The nonlinear dynamics stem from the nonlinear differential equation that describes the rotation of a three dimensional body. Initial results, using simulated data, and real Rossi X ray Timing Explorer (RXTE) data indicate that the algorithm is efficient and robust.

Azor, Ruth↗