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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

A Framework for Evaluating Distributed Electric Propulsion on the SUSAN Electrofan Aircraft

This work presents a framework for evaluating models and algorithms for Distributed Electric Propulsion (DEP) on the SUSAN Electrofan Aircraft. Throughout the development of the SUSAN aircraft, the performance of various configurations of the aircraft will need to be analyzed. However, the static behavior alone is not sufficient to describe the performance of these configurations. Therefore, simulation with fully integrated subsystem models is required. The proposed framework considers the vehicle aerodynamic, propulsion, and control subsystems. The presented framework automatically generates control laws for any vehicle configuration in response to changes in these subsystems. To compare these different vehicle configurations, various time and frequency domain performance metrics are compared. Three different system modifications are used as cases to evaluate this framework. The first modification integrates the propulsion control system with the flight controller to enable differential thrust without stalling the main engine. This evaluation case is used to validate the framework for aircraft configurations with coupled subsystems. The second modification compares the effect of the vertical tail size on open and closed loop performance. This evaluation case is used to validate the framework for controlling different configurations and tuning towards comparable closed loop performance despite changes to the aircraft's aerodynamic model. The third modification implements two different control allocation schemes. This evaluation case demonstrates the framework's ability to evaluate allocation modifications needed to take advantage of DEP. The first evaluation case is used to show that controller integration enables differential thrust, improving realized wingfan bandwidth by up to 40\% in simulation. The second evaluation case demonstrates that the framework can stabilize the reduced tail size aircraft with closed loop control. The third evaluation case demonstrates that a pseudoinverse control allocation scheme improves lateral velocity settling time by approximately 17~seconds over a symmetric-thrust allocation. These cases show that the framework is useful for evaluating the performance of integrated system designs, enabling analyses of new models and algorithms for the SUSAN distributed electric propulsion vehicle.

Nicholas C Ogden↗
NASA NTRS1994

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
NASA NTRS1992

A new determination of the solar rotation rate

We use 'stackplot' displays to compare observations of the photospheric magnetic field during sunspot cycle 21 with simulations based on the flux-transport model. Adopting nominal rates of diffusion, differential rotation, and meridional flow, we obtain slanted patterns similar to those of the observed field, even when the sources of flux are assigned random longitudes in the model. At low latitudes, the slopes of the nearly vertical patterns of simulated field are sensitive to the rotation rate used in the calculation, and insensitive to the rates of diffusion and flow during much of the sunspot cycle. Good agreement between the observed and simulated patterns requires a synodic equatorial rotation period of 26.75 +/- 0.05 days.

Sheeley, N. R., Jr.↗
NASA NTRS1989

Multiple symbol differential detection of uncoded and trellis coded MPSK

A differential detection for MPSK, which uses a multiple symbol observation interval, is presented and its performance analyzed and simulated. The technique makes use of maximum-likelihood sequence estimation of the transmitted phases rather than symbol-by-symbol detection as in conventional differential detection. As such the performance of this multiple symbol detection scheme fills the gap between conventional (two-symbol observation) differentially coherent detection of MPSK and ideal coherent of MPSK with differential encoding. The amount of improvement gained over conventional differential detection depends on the number of phases, M, and the number of additional symbol intervals added to the observation. What is particularly interesting is that substantial performance improvement can be obtained for only one or two additional symbol intervals of observation. The analysis and simulation results presented are for uncoded and trellis coded MPSK.

Divsalar, Dariush↗
NASA NTRS1982

Evaluation of LANDSAT-D Orbit Determination Using a Filter/Smoother (PREFER)

Simulated range and range rate data for five tracking stations were first generated using batch least squares orbit determination (GTDS). Then GTDS was used (in the differential correction mode) to produce a nominal trajectory which was input to PREFER. The GTDS differential correction (DC) run was made using models which differed from those used to produce the simulated data. These model differences were chosen to be fairly realistic approximations to the errors in the models actually used for operational orbit determination. Several different simulation runs were made with different types of model errors in order to determine the sensitivity to these errors. The nominal trajectory and the simulated measurement data were input to PREFER to produce a smoothed ephemeris file. Numerous runs of PREFER were made in which parameters describing the statistics of the model errors were varied. The likelihood function computed by the Kalman filter determined the ""best'' choice of input parameters. There was strong negative correlation between the likelihood function and the errors in the smoothed ephemeris.

Gibbs, B. P.↗
NASA NTRS1994

Monte Carlo simulation of the terrestrial hydrogen exosphere

Methods for Monte Carlo simulation of planetary exospheres have evolved from early work on the lunar atmosphere, where the regolith surface provides a well defined exobase. A major limitation of the successor simulations of the exospheres of Earth and Venus is the use of an exobase surface as an artifice to separate the collisional processes of the thermosphere from a collisionles exosphere. In this paper a new generalized approach to exosphere simulation is described, wherein the exobase is replaced by a barometric depletion of the major constitents of the thermosphere. Exospheric atoms in the thermosphere-exosphere transition region, and in the outer exosphere as well, travel in ballistic trajectories that are interrupted by collisons with the background gas, and by charge exchange interactions with ionospheric particles. The modified simulator has been applied to the terrestrial hydrogen exosphere problem, using velocity dependent differential cross sections to provide statistically correct collisional scattering in H-O and H-H(+) interactions. Global models are presented for both solstice and equinox over the effective solar cycle range of the F(sub 10.7) index (80 to 230). Simulation results show significant differences with previous terrestrial exosphere models, as well as with the H distributions of the MSIS-86 thermosphere model.

Hodges, R. Richard, Jr.↗
NASA NTRS1990

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

A model is developed that permits the representation of the Galactic disk and bulge, the molecular ring, extragalactic sky, and other Galactic structures in terms of the IR point-source sky. Full characterization is achieved for up to 87 types of Galactic sources in terms of scale heights and absolute magnitudes, and the model is guided by a parallel Monte Carlo simulation of the Galaxy at 12 microns. The model predicts differential and cumulative IR source counts for all bandpasses in the IRAS ranges and can therefore be employed to predict the character of the point-source sky related to IR space observations.

Wainscoat, Richard J.↗
NASA NTRS1995

Investigation of dispersion-relation-preserving scheme and spectral analysis methods for acoustic waves

Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.

Vanel, Florence O.↗
NASA NTRS1999

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur R.↗
NASA NTRS1999

A Viscoelastic Hybrid Shell Finite Element

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur↗
NASA NTRS2007

The Cosmic Battery Revisited

We reinvestigate the generation and accumulation of magnetic flux in optically thin accretion flows around active gravitating objects. The source of the magnetic field is the azimuthal electric current associated with the Poynting-Robertson drag on the electrons of the accreting plasma. This current generates magnetic field loops which open up because of the differential rotation of the flow. We show through simple numerical simulations that what regulates the generation and accumulation of magnetic flux near the center is the value of the plasma conductivity. Although the conductivity is usually considered to be effectively infinite for the fully ionized plasmas expected near the inner edge of accretion disks, the turbulence of those plasmas may actually render them much less conducting due to the presence of anomalous resistivity. We have discovered that if the resistivity is sufficiently high throughout the turbulent disk while it is suppressed interior to its inner edge, an interesting steady-state process is established: accretion carries and accumulates magnetic flux of one polarity inside the inner edge of the disk, whereas magnetic diffusion releases magnetic flux of the opposite polarity to large distances. In this scenario, magnetic flux of one polarity grows and accumulates at a steady rate in the region inside the inner edge and up to the point of equipartition when it becomes dynamically important. We argue that this inward growth and outward expulsion of oppositely-directed magnetic fields that we propose may account for the approx. 30 min cyclic variability observed in the galactic microquasar GRS1915+105.

Contopoulos, Ioannis↗
NASA NTRS2007

Simulated Microgravity Regulates Gene Transcript Profiles of 2T3 Preosteoblasts: Comparison of the Random Positioning Machine and the Rotating Wall Vessel Bioreactor

Microgravity of spaceflight induces bone loss due in part to decreased bone formation by osteoblasts. We have previously examined the microgravity-induced changes in gene expression profiles in 2T3 preosteoblasts using the Random Positioning Machine (RPM) to simulate microgravity conditions. Here, we hypothesized that exposure of preosteoblasts to an independent microgravity simulator, the Rotating Wall Vessel (RWV), induces similar changes in differentiation and gene transcript profiles, resulting in a more confined list of gravi-sensitive genes that may play a role in bone formation. In comparison to static 1g controls, exposure of 2T3 cells to RWV for 3 days inhibited alkaline phosphatase activity, a marker of differentiation, and downregulated 61 genes and upregulated 45 genes by more than two-fold as shown by microarray analysis. The microarray results were confirmed with real time PCR for downregulated genes osteomodulin, bone morphogenic protein 4 (BMP4), runx2, and parathyroid hormone receptor 1. Western blot analysis validated the expression of three downregulated genes, BMP4, peroxiredoxin IV, and osteoglycin, and one upregulated gene peroxiredoxin I. Comparison of the microarrays from the RPM and the RWV studies identified 14 gravi-sensitive genes that changed in the same direction in both systems. Further comparison of our results to a published database showing gene transcript profiles of mechanically loaded mouse tibiae revealed 16 genes upregulated by the loading that were shown to be downregulated by RWV and RPM. These mechanosensitive genes identified by the comparative studies may provide novel insights into understanding the mechanisms regulating bone formation and potential targets of countermeasure against decreased bone formation both in astronauts and in general patients with musculoskeletal disorders.

Patel, Mamta J.↗
NASA NTRS2014

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗
NASA NTRS2014

Solid Rocket Launch Vehicle Explosion Environments

Empirical explosion data from full scale solid rocket launch vehicle accidents and tests were collected from all available literature from the 1950s to the present. In general data included peak blast overpressure, blast impulse, fragment size, fragment speed, and fragment dispersion. Most propellants were 1.1 explosives but a few were 1.3. Oftentimes the data from a single accident was disjointed and/or missing key aspects. Despite this fact, once the data as a whole was digitized, categorized, and plotted clear trends appeared. Particular emphasis was placed on tests or accidents that would be applicable to scenarios from which a crew might need to escape. Therefore, such tests where a large quantity of high explosive was used to initiate the solid rocket explosion were differentiated. Also, high speed ground impacts or tests used to simulate such were also culled. It was found that the explosions from all accidents and applicable tests could be described using only the pressurized gas energy stored in the chamber at the time of failure. Additionally, fragmentation trends were produced. Only one accident mentioned the elusive "small" propellant fragments, but upon further analysis it was found that these were most likely produced as secondary fragments when larger primary fragments impacted the ground. Finally, a brief discussion of how this data is used in a new launch vehicle explosion model for improving crew/payload survival is presented.

Richardson, E. H.↗
NASA NTRSDate not supplied

Reference Command Optimization for the Transition Flight Mode of a Lift plus Cruise Vehicle

Advanced air mobility mainly utilizes vehicles that are capable of vertical takeoff and landing (VTOL) for the simplicity of operation and large-scale deployment. However, VTOL vehicles need specialized trajectory and command design for the transition phase, where the vehicles transition between rotor-borne flight and wing-borne flight. Since VTOL vehicles are commonly designed as over-actuated systems for redundancy, one challenge that arises is actuator ambiguity, where it is unclear how to uniquely command actuators for the VTOL vehicle to track a given trajectory. We propose a method to design optimal reference commands for the transition mode. By formulating an 𝑙 1 -norm cost function on the rotor thrusts of the vehicle, we can achieve economical operation of the rotors such that they only operate when necessary and efficiently utilize the aerodynamics to save energy from reduced rotor actuation. We validate our approach in simulations and show its benefit compared to the commonly used differential flatness-based method.

John L Bullock↗
NASA NTRS2026

Experimental Demonstration of a Supersonic Spike Ejector

Spike diffusers are presented as a means to improve supersonic ejector efficiency via reduction in shock losses. Design, fabrication, and testing of an articulated prototype is described. Experimental results are compared to empirical performance curves representing the loci of optimal one- and two-stage conventional ejectors driven by dry air and discharging to atmosphere. Operated at a compression ratio of ~10 with all inputs equal, the spike ejector is shown to have achieved up to 2.19X and 1.83X the secondary mass throughput of the one- and two-stage references, respectively. This suggests the potential of spike diffusers to reduce ejector operating costs by ~50%.

wind tunnel↗
NASA NTRS2000

A Composite View of Ozone Evolution in the 1995-1996 Northern Winter Polar Vortex Developed from Airborne Lidar and Satellite Observations

The processes which contribute to the ozone evolution in the high latitude northern lower stratosphere are evaluated using a three dimensional model simulation and ozone observations. The model uses winds and temperatures from the Goddard Earth Observing System Data Assimilation System. The simulation results are compared with ozone observations from three platforms: the differential absorption lidar (DIAL) which was flown on the NASA DC-8 as part of the Vortex Ozone Transport Experiment; the Microwave Limb Sounder (MLS); the Polar Ozone and Aerosol Measurement (POAM II) solar occultation instrument. Time series for the different data sets are consistent with each other, and diverge from model time series during December and January. The model ozone in December and January is shown to be much less sensitive to the model photochemistry than to the model vertical transport, which depends on the model vertical motion as well as the model vertical gradient. We evaluate the dependence of model ozone evolution on the model ozone gradient by comparing simulations with different initial conditions for ozone. The modeled ozone throughout December and January most closely resembles observed ozone when the vertical profiles between 12 and 20 km within the polar vortex closely match December DIAL observations. We make a quantitative estimate of the uncertainty in the vertical advection using diabatic trajectory calculations. The net transport uncertainty is significant, and should be accounted for when comparing observations with model ozone. The observed and modeled ozone time series during December and January are consistent when these transport uncertainties are taken into account.

Douglass, A. R.↗