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 451 records · Page 25

Multibody Parachute Flight Simulations for Planetary Entry Trajectories Using "Equilibrium Points"

A method has been developed to reduce numerical stiffness and computer CPU requirements of high fidelity multibody flight simulations involving parachutes for planetary entry trajectories. Typical parachute entry configurations consist of entry bodies suspended from a parachute, connected by flexible lines. To accurately calculate line forces and moments, the simulations need to keep track of the point where the flexible lines meet (confluence point). In previous multibody parachute flight simulations, the confluence point has been modeled as a point mass. Using a point mass for the confluence point tends to make the simulation numerically stiff, because its mass is typically much less that than the main rigid body masses. One solution for stiff differential equations is to use a very small integration time step. However, this results in large computer CPU requirements. In the method described in the paper, the need for using a mass as the confluence point has been eliminated. Instead, the confluence point is modeled using an "equilibrium point". This point is calculated at every integration step as the point at which sum of all line forces is zero (static equilibrium). The use of this "equilibrium point" has the advantage of both reducing the numerical stiffness of the simulations, and eliminating the dynamical equations associated with vibration of a lumped mass on a high-tension string.

Raiszadeh, Ben↗

A real-time subsynoptic, meso, and microscale severe storm forecast system

A system of differential equations is integrated numerically in space and time over several different matrices in an effort to simulate the atmospheric wave structures which organize severe local storms. Preliminary results with case studies indicate that dynamical fields produced by the numerical simulations can be translated into very fine scale space and time zones where severe storm forecast indices can be developed. These fine scale indices are now available in real-time when run on Langley's STAR 100 computer system.

Kaplan, M. L.↗

A mathematical model of post-instability in continuum mechanics

The post-instability can appear in the form of failure of hyperbolicity or in the form of cascade instability (stretching of vorticity in an inviscid fluid or instability of the Navier-Stokes equations). All the cases are characterized by an unlimited decrease of the scale of the motions, in the course of which the derivatives of the corresponding functions tend to infinity, though the functions themselves remain finite. Such an instability shows that the corresponding state of a continuum cannot be properly simulated by smooth functions. Hence, the original model must be corrected by giving up the requirement about differentiability and enlarging the class of functions. Applications to a model of an inviscid fluid and to a model of turbulence are discussed.

Zak, M.↗

Effect of suspension hypokinesia/hypodynamia on glucocorticoid receptor levels in rat hindlimb muscles

Ilyina-Kakueva et. al. (1976) conducted an investigation in which rats were exposed to weightlessness during the Cosmos program. An examination of the rats revealed a marked atrophy of hindlimb muscles. A suspension model has been developed to simulate these weightlessness-induced alterations. In agreement with the Cosmos studies, suspension hypokinesia/hypodynamia (H/H) results in differential atrophy of hindlimb muscles in rats. Recent studies have demonstrated elevated glucocorticoid receptor numbers in the gastrocnemius muscle following immobilization and denervation. One of the objectives of the present investigation was to evaluate the effect of suspension H/H on glucocorticoid receptor levels in rat hindlimb muscles. Another objective was to ascertain whether altered receptor levels reflect the differential nature of hindlimb muscle atrophy during suspension H/H. The obtained findings suggest that differential muscle atrophy resulting from H/H may result from differential alterations of glucocorticoid receptor levels.

Steffen, J. M.↗

DEAN - A program for Dynamic Engine ANalysis

The Dynamic Engine Analysis Program, DEAN, is a FORTRAN code implemented on the IBM/370 mainframe at NASA Lewis Research Center for digital simulation of turbofan engine dynamics. DEAN is an interactive program which allows the user to simulate engine subsystems as well as full engine systems with relative ease. The nonlinear first order ordinary differential equations which define the engine model may be solved by one of four integration schemes, a second order Runge-Kutta, a fourth order Runge-Kutta, an Adams Predictor-Corrector, or Gear's method for still systems. The numerical data generated by the model equations are displayed a specified intervals between which the user may choose to modify various parameters affecting the model equations and transient execution. Following the transient run, versatile graphics capabilities allow close examination of the data. DEAN's modeling procedure and capabilities are demonstrated by generating a model of simple compressor rig.

Sadler, G. G.↗

Statistical analysis of near-wall structures in turbulent channel flow

The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Moser, Robert D.↗

Combined trellis coding and feedforward processing for MSS applications

The idea of using a multiple (more than two) symbol observation interval to improve error probability performance is applied to differential detection of trellis coded MPSK over a mobile satellite (fading) channel. Results are obtained via computer simulation. It is shown that only a slight increase (e.g., one symbol) in the length of the observation interval will provide a significant improvement in bit error probability performance both in AWGN and fading environments.

Divsalar, D.↗

A model of the 8-25 micron point source infrared sky

We present a detailed model for the IR point-source sky that comprises geometrically and physically realistic representations of the Galactic disk, bulge, stellar halo, spiral arms (including the 'local arm'), molecular ring, and the extragalactic sky. We represent each of the distinct Galactic components by up to 87 types of Galactic source, each fully characterized by scale heights, space densities, and absolute magnitudes at BVJHK, 12, and 25 microns. The model is guided by a parallel Monte Carlo simulation of the Galaxy at 12 microns. The content of our Galactic source table constitutes a good match to the 12 micron luminosity function in the simulation, as well as to the luminosity functions at V and K. We are able to produce differential and cumulative IR source counts for any bandpass lying fully within the IRAS Low-Resolution Spectrometer's range (7.7-22.7 microns as well as for the IRAS 12 and 25 micron bands. These source counts match the IRAS observations well. The model can be used to predict the character of the point source sky expected for observations from IR space experiments.

Wainscoat, Richard J.↗

Vacuum-Assisted, Constant-Force Exercise Device

The vacuum-assisted, constant-force exercise device (VAC-FED) has been proposed to fill a need for a safe, reliable exercise machine that would provide constant loads that could range from 20 to 250 lb (0.09 to 1.12 kN) with strokes that could range from 6 to 36 in. (0.15 to 0.91 m). The VAC-FED was originally intended to enable astronauts in microgravity to simulate the lifting of free weights, but it could just as well be used on Earth for simulated weight lifting and other constant-force exercises. Because the VAC-FED would utilize atmospheric/vacuum differential pressure instead of weights to generate force, it could weigh considerably less than either a set of free weights or a typical conventional exercise machine based on weights. Also, the use of atmospheric/ vacuum differential pressure to generate force would render the VAC-FED inherently safer, relative to free weights and to conventional exercise machines that utilize springs to generate forces. The overall function of the VAC-FED would be to generate a constant tensile force in an output cable, which would be attached to a bar, handle, or other exercise interface. The primary force generator in the VAC-FED would be a piston in a cylinder. The piston would separate a volume vented to atmosphere at one end of the cylinder from an evacuated volume at the other end of the cylinder (see figure). Hence, neglecting friction at the piston seals, the force generated would be nearly constant equal to the area of the piston multiplied by the atmospheric/vacuum differential pressure. In the vented volume in the cylinder, a direct-force cable would be looped around a pulley on the piston, doubling the stroke and halving the tension. One end of the direct-force cable would be anchored to a cylinder cap; the other end of the direct-force cable would be wrapped around a variable-ratio pulley that would couple tension to the output cable. As its name suggests, the variable-ratio pulley would contain a mechanism that could be used to vary the ratio between the tension in the direct-force cable and the tension in the output cable. This mechanism could contain gears, pulleys, and/or levers, for example.

Hansen, Christopher P.↗

Observational Effects of Magnetism in O Stars: Surface Nitrogen Abundances

Aims. We investigate the surface nitrogen content of the six magnetic O stars known to date as well as of the early B-type star Tau Sco.. We compare these abundances to predictions of evolutionary models to isolate the effects of magnetic field on the transport of elements in stellar interiors. Methods. We conduct a quantitative spectroscopic analysis of the ample stars with state-of-the-art atmosphere models. We rely on high signal-to-noise ratio, high resolution optical spectra obtained with ESPADONS at CFHT and NARVAL at TBL. Atmosphere models and synthetic spectra are computed with the code CMFGEN. Values of N/H together with their uncertainties are determined and compared to predictions of evolutionary models. Results. We find that the magnetic stars can be divided into two groups: one with stars displaying no N enrichment (one object); and one with stars most likely showing extra N enrichment (5 objects). For one star (Ori C) no robust conclusion can be drawn due to its young age. The star with no N enrichment is the one with the weakest magnetic field, possibly of dynamo origin. It might be a star having experienced strong magnetic braking under the condition of solid body rotation, but its rotational velocity is still relatively large. The five stars with high N content were probably slow rotators on the zero age main sequence, but they have surface N/H typical of normal O stars, indicating that the presence of a (probably fossil) magnetic field leads to extra enrichment. These stars may have a strong differential rotation inducing shear mixing. Our results shOuld be viewed as a basis on which new theoretical simulations can rely to better understand the effect of magnetism on the evolution of massive stars.

Martins, F.↗

Experimental and Numerical Investigation of Adsorption/Desorption in Packed Sorption Beds Under Ideal and Non-Ideal Flows

The importance of the wall effect on packed beds in the adsorption and desorption of carbon dioxide, nitrogen, and water on molecular sieve 5A of 0.127 cm in radius is examined experimentally and with one-dimensional computer simulations. Experimental results are presented for a 22.5-cm long by 4.5-cm diameter cylindrical column with concentration measurements taken at various radial locations. The set of partial differential equations are solved using finite differences and Newman's method. Comparison of test data with the axial-dispersed, non-isothermal, linear driving force model suggests that a two-dimensional model (submitted to Separation Science and Technology) is required for accurate simulation of the average column breakthrough concentration. Additional comparisons of test data with the model provided information on the interactive effects of carrier gas coadsorption with CO2, as well as CO2-H2O interactions.

Mohamadinejad, H.↗

Lidar temperature profiling - Performance simulations of Mason's method

In Mason's method (1975) atmospheric temperatures are inferred from a measure of the Boltzmann distribution of rotational states in one of the vibrational bands of O2. Differential absorption is measured using three tunable, narrowband pulse lasers. The outputs of two are tuned to wavelengths at the centers of absorption lines at either end of a particular branch in the band; the third wavelength is in a region of no absorption. The temperature-altitude profile can be calculated from the ratio of the two line absorption coefficients plus a priori knowledge of the line parameters. In the present paper, computer simulations of various lidar configurations are made, using different line pairs in the atmospheric bands of O2 (approximately 630, 690, and 760 nm). Simulated results are presented for temperature profiles measured from a Space Shuttle lidar.

Schwemmer, G. K.↗

Interfacial Structure and Capacitance of Li-doped Ionic Liquid Electrolytes from Molecular Simulations

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries, supercapacitors, and hybrid energy storage devices. Though Li-salt is often present in these systems, its influence on interfacial properties is largely uncharacterized. We, thereby, present an extensive computational analysis, supported by experimental comparisons, of the properties of a representative set of these electrolytesat an ideal carbon interface as a function of Li-salt doping and voltage. We have performed polarizable molecular (MD) dynamics simulations, using the APPLEP force field, to evaluate electric double layer (EDL) capacitance and distribution of Li+ in the EDL. Differential capacitance exhibits the characteristic camel profile and is insensitive to Li-doping. Li+ localizes in the second molecular layer of the EDL, which is a result of confinement from free energy barriers associated with ion layering. Joint MDelectronic structure computations show the electrochemical window of the electrolytes to be a weak function of Li-doping. Estimates of supercapacitor specific energy are made using the computed window and capacitance. The magnitude and trends in specific energy are in good agreement with experiment.

batteries↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

DIAL simulations for satellite water vapor profile measurements within the BEST project

The project BEST has been proposed by several laboratories of the French Space Agency (CNES) as part of the GEWEX program. Its objective is the study of the energy budget of the tropical system, and aims at measuring different components of the water cycle by coupling several observing systems. As part of this, Service d'Aeronomie (SA) proposed to study a Differential Absorption Lidar (DIAL) system to measure tropical water vapor profiles and to determine humidity convergence fields, when coupled with the Doppler wind lidar. We have performed simulation analyses for water vapor measurement from the selected orbit of 400 km. For these simulations, we have adopted an atmospheric model developed at SA. We have optimized, by a Monte Carlo statistical method, the system parameters and evaluated bias and standard deviation of the retrieved water vapor mixing ratio.

Chazette, P.↗

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗

Ground testing in a simulated auroral environment

Shuttle polar orbit missions are being planned in which astronaut extravehicular activities (EVA) will eventually include external operations or repair and refurbishment of satellites. Recent satellite data and model calculations show that significant spacecraft charging and differential charging occurs in that environment and could endanger equipment and personnel involved in the EVA. These spacecraft charging events take place in response to an intense flux of high energy electrons accompanied by a large drop in the ambient plasma density. The extent of which the auroral environment can be simulated in a large vacuum chamber are determined so that conditions which lead to significant charging and discharging may be observed and controlled. A large chamber would be required in order to simulate the effects of an aurora-like, large cross-section, energetic electron precipitation on a target system the size of an astronaut with a life support system. The precipitating electrons can be simulated with an array of multipactor sources and the background ionization produced by the beam can be limited to the range of the auroral density if the chamber pressure is not significantly above 1 x .000001 Torr. Some small chamber work has shown that the plasma in the beam will diffuse radially at a rate determined by the Bohm diffusion coefficient.

Mcintyre, Bernard↗

Aspects of efficient and reliable multibody system simulation

Multibody system equations can be generated in various forms. All of these may be interpreted as results of two basic approaches, the augmentation- and the elimination-method. The former method yields the descriptor form of the system motion, a set of differential-algebraic equations (DAE), and the latter the state space representation, a minimal set of ordinary differential equations (ODE). Both of these methods are surveyed. Particular emphasis is on the discussion of recursive computational schemes, generating the equations of motion with a number of operations, which is proportional to the number N of system bodies (O(N)-formulations). For simulation purposes one would like to create that set of system equations, which can be generated most efficiently and for which the most efficient and reliable solution techniques are available. Numerical solution techniques for ODE have been studied in great detail and they are well-developed. By contrast, DAE have not been investigated for such a long time. In view of new developments in the latter field the generation of all the equations required for an efficient and reliable solution of DAE describing multibody system motion is discussed. These methods, i.e., an O(N)-formulation and new techniques for solving DAE, are implemented in the SIMPACK code. Its capabilities are illustrated by simulation of multibody robot models.

Schwertassek, R.↗