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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗
NASA NTRS2011

The Angular Momentum of Baryons and Dark Matter Halos Revisited

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive-mesh refinement (AMR), we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r = 0.1R(sub vir). In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/R(sub vir) > 0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons (gas and stars) remains close to that of dark matter halos. Indeed, our simulations indicate that any differential loss of angular momentum amplitude between the two components is minor even though dark matter halos continuously lose between half and two-thirds of their specific angular momentum modulus as they evolve. In light of our results, a substantial revision of the standard theory of disc formation seems to be required. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

Kimm, Taysun↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
NASA NTRS1985

Identification of unsteady aerodynamics and aeroelastic integro-differential systems

The problem of estimating integro-differential models based on test or simulation data is dealt with. The identification techniques proposed for estimating parameters in models described by differential equations need to be considerably extended to deal with the integral terms. Conditions under which the integral terms may be approximated by algebraic values are discussed. The integro-differential models discussed are related to indicial models proposed by aerodynamicists to describe unsteady flow.

Gupta, N. K.↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
NASA NTRS1978

Numerical derivation of linear perturbation model from nonlinear real-time vehicle simulations - A means of validation and increased credibility

A numerical technique has been developed which generates linear perturbation models from nonlinear real-time aircraft vehicle simulations at Langley Research Center. The technique is very general and can be applied to simulations of any system described by nonlinear differential equations. Linear models about various trim conditions have been obtained from simulations of the Boeing 737, Northrop F5-E, Cessna 180, Sikorsky CH-54, and NASA/Army Rotor Systems Research Aircraft (RSRA). One use of the derived linear models is as a validation and verification tool for the nonlinear simulation.

Dieudonne, J. E.↗
NASA NTRS1985

Actual and Simulated Weightlessness Inhibit Osteogenesis in Long Bone Metaphysis by Different Mechanisms

Weightlessness and simulated weightlessness inhibit the rate of periosteal bone formation in long bones. Formation of preosteoblasts is suppressed in periodontal ligament (PDL) of maxillary molars, which suggests a generalized block in osteoblast histogenesis. Growth in length of long bones is decreased by simulated weightlessness, but there are no reliable data on the influence of actual weightlessness on metaphyseal growth. The nuclear size assay for assessing relative numbers of osteoblast precursor cells was utilized in the primary spongiosa of growing long bones subjected to actual and simulated weightlessness. It is found that: (1) Actual weightlessness decreases total number of osteogenic cells and inhibits differentiation of osteoblast precursor cells, (2) Simulated weightlessness suppresses only osteoblast differentation; and (3) The nuclear morphometric assay is an effective means of assessing osteogenic activity in the growing metaphysis or long bones.

Roberts, W. E.↗
NASA NTRS1981

Transient dynamic analysis of high-speed lightly loaded cylindrical roller bearings. 2: Computer program and results

The governing differential equations of motion for a high speed cylindrical roller bearing are programmed for numerical solution and plotted output. This computer program has the capability of performing a two dimensional or three dimensional simulation. Two numerical solutions of the governing differential equations were obtained to simulate the motion of a roller bearing, one for the two dimensional equations of motion and one for the three dimensional equations of motion. Computer generated plots were obtained and present such data as roller/cage interaction forces, roller/race traction forces, roller/race relative slip velocities and cage angular speed over a nondimensional time equivalent to 1.2 revolutions of the inner race. Roller axial displacement, roller skew angle, and skew moment are also plotted for the three dimensional solution. The trajectory of the cage center is plotted for both the two dimensional and three dimensional solutions.

Conry, T. F.↗
NASA NTRS1991

Propellant-remaining modeling

A successful satellite mission is predicted upon the proper maintenance of the spacecraft's orbit and attitude. One requirement for planning and predicting the orbit and attitude is the accurate estimation of the propellant remaining onboard the spacecraft. Focuss is on the three methods that were developed for calculating the propellant budget: the errors associated with each method and the uncertainties in the variables required to determine the propellant remaining that contribute to these errors. Based on these findings, a strategy is developed for improved propellant-remaining estimation. The first method is based on Boyle's law, which related the values of pressure, volume, and temperature (PVT) of an ideal gas. The PVT method is used for the monopropellant and the bipropellant engines. The second method is based on the engine performance tests, which provide data that relate thrust and specific impulse associated with a propellant tank to that tank's pressure. Two curves representing thrust and specific impulse as functions of pressure are then generated using a polynomial fit on the engine performance data. The third method involves a computer simulation of the propellant system. The propellant flow is modeled by creating a conceptual model of the propulsion system configuration, taking into account such factors as the propellant and pressurant tank characteristics, thruster functionality, and piping layout. Finally, a thrust calibration technique is presented that uses differential correction with the computer simulation method of propellant-remaining modeling. Thrust calibration provides a better assessment of thruster performance and therefore enables a more accurate estimation of propellant consumed during a given maneuver.

Torgovitsky, S.↗
NASA NTRS1990

Simulation of the collapse and dissipation of Langmuir wave packets

Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.

Newman, D. L.↗
NASA NTRS2002

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗
NASA NTRS2016

Reserves and Trade Jointly Determine Exposure to Food Supply Shocks

While a growing proportion of global food consumption is obtained through international trade, there is an ongoing debate on whether this increased reliance on trade benefits or hinders food security, and specifically, the ability of global food systems to absorb shocks due to local or regional losses of production. This paper introduces a model that simulates the short-term response to a food supply shock originating in a single country, which is partly absorbed through decreases in domestic reserves and consumption, and partly transmitted through the adjustment of trade flows. By applying the model to publicly-available data for the cereals commodity group over a 17 year period, we find that differential outcomes of supply shocks simulated through this time period are driven not only by the intensification of trade, but as importantly by changes in the distribution of reserves. Our analysis also identifies countries where trade dependency may accentuate the risk of food shortages from foreign production shocks; such risk could be reduced by increasing domestic reserves or importing food from a diversity of suppliers that possess their own reserves. This simulation-based model provides a framework to study the short-term, nonlinear and out-of-equilibrium response of trade networks to supply shocks, and could be applied to specific scenarios of environmental or economic perturbations.

supply↗
NASA NTRS2000

Processing CCD Images to Detect Transits of Earth-Sized Planets: Maximizing Sensitivity While Achieving Reasonable Downlink Requirements

We have performed end-to-end laboratory and numerical simulations to demonstrate the capability of differential photometry under realistic operating conditions to detect transits of Earth-sized planets orbiting solar-like stars. Data acquisition and processing were conducted using the same methods planned for the proposed Kepler Mission. These included performing aperture photometry on large-format CCD images of an artificial star fields obtained without a shutter at a readout rate of 1 megapixel/sec, detecting and removing cosmic rays from individual exposures and making the necessary corrections for nonlinearity and shutterless operation in the absence of darks. We will discuss the image processing tasks performed `on-board' the simulated spacecraft, which yielded raw photometry and ancillary data used to monitor and correct for systematic effects, and the data processing and analysis tasks conducted to obtain lightcurves from the raw data and characterize the detectability of transits. The laboratory results are discussed along with the results of a numerical simulation carried out in parallel with the laboratory simulation. These two simulations demonstrate that a system-level differential photometric precision of 10-5 on five- hour intervals can be achieved under realistic conditions.

Earth-size planets↗
NASA NTRSDate not supplied

Differential Gene Expression in A Cross-Feeding Two-Species Model Microbial Community Under Simulated Microgravity and Deep-Space Radiation

A long-term goal of space biology is to understand interspecies microbial interactions in space. Presently, little is known about the combined effect of microgravity and ionizing radiation on bacterial community response when species are interdependent through exchange of metabolites in fluid medium (cross-feeding). Microgravity is expected to slow interspecies mass transfer and growth in cross-feeding communities in the low-shear, diffusion-limited environment, while ionizing radiation may influence stress response to direct (DNA damage) and indirect damage (ROS). Using a well-understood, two-species (Escherichia coli and Salmonella enterica) microbial community engineered to be a model for studying cross-feeding, we simulated galactic cosmic rays (GCRsim) and microgravity to test the hypothesis: exposure to ionizing radiation causes cell damage or stress, altering transcriptomic community responses in metabolically interdependent cells, which is exacerbated by microgravity. We expect to see differential gene expression between cross-feeding and non-cross-feeding communities. We measured GCRsim effects on growth and gene expression in well-mixed versus simulated-microgravity conditions and in cross-feeding and non-cross-feeding medium. Microbial cultures were inoculated into liquid medium in rotating wall vessels (RWV) with different rotation rates: 5 RPM (simulated microgravity) and 50 RPM (well-mixed). The E. coli-S. enterica consortium, under simulated microgravity, were exposed to 500 mGy of Simplified 5-ion Galactic Cosmic Ray Simulation for 2 hours at Brookhaven National Lab. We harvested samples 40 minutes after irradiation for extraction and sequencing (NASA GeneLab). Here we present the differential gene expression analysis results, which reveal altered transcriptomic community responses, even where growth rate differences are not observed. Gene expression of these actively metabolizing microbial communities in GCRsim may illuminate molecular mechanisms of microbial interactions in space. Understanding how microbial community gene expression, metabolism, and other cellular processes are influenced by spaceflight stressors can inform the use of microbes in human life support for low Earth orbit missions and beyond.

microgravity↗
NASA NTRS1984

Effects of suspension on tissue levels of glucocorticoid receptors

Differential muscle responses can be simulated by hypokinetic/hypodynamic (H/H) suspension of rats with complete unloading of the hindlimb muscles. Since mechanism(s) underlying these atrophic effects were not clearly elucidated, experiments were initiated to investigate a possible role for glucocorticoids in the physiological and biochemical responses to H/H. The principal objective was to assess the potential for alterations in peripheral responsiveness to glucocorticoids in response to H/H. Studies have initially focused on the determination of tissue levels of glucocorticoid receptors as one index of hormonal sensitivity at the cellular level. Four hindlimb muscles (soleus, gastrocnemius, plantaris and EDL), previously demonstrated to exhibit differential responses to H/H, were investigated. Receptor levels in other glucocorticoid sensitive tissues (heart, liver, and kidney) were determined. Male rats (180-200g) were suspended for 7 or 14 days, sacrificed by cervical dislocation, and the tissues excised.

Steffen, J. M.↗
NASA NTRS1987

Simulation of AFC for a DMSK Receiver

Low bit rates suitable for low-cost satellite communications. LOOP computer program written to simulate automatic frequency control (AFC) subsystem of differential minimum-shift keying (DMSK) receiver with bit rate of 2,400 baud. AFC simulated by LOOP is first-order loop configuration with first-order resistance-and-capacitance filter. Written in MS SuperSoft FORTRAN 77.

Davarian, Faramaz↗
NASA NTRS1998

X-33 Integrated Test facility, Extended Range Simulation

In support of the X-33 Single Stage To Orbit program, NASA Dryden Flight Research Center was selected to provide continuous communications coverage of the X-33 vehicle from launch, through landing at Malmstrom Air Force Base, Montana and Michaels Army Air Field, Utah. An extensive real-time range simulation capability is being developed to ensure successful communications with the autonomous X-33 vehicle. This paper will provide an overview of the various levels of simulation being developed to support the X-33 extended range subsystems. These subsystems include the Flight Termination System, L-Band command uplink subsystem and the S-Band telemetry downlink subsystem. In addition, the radar model developed provides continuous azimuth, elevation and range information based on the flight trajectory. The Dynamic Ground Station Analysis model developed by NASA Goddard Space Flight Center, calculate the received signal strength at each ground station. This model takes into consideration Radio Frequency (RF) link parameters such as frequency, antenna gain, space loss, plasma effects and the vehicle's position and attitude at any point in time during the flight path. All three RF links are then attenuated based on this calculated level and the RF signals are sent into telemetry receivers to emulate remote sites, or the power incident on the vehicle from uplinked signals. The best source received telemetry data is then passed back to the Launch and Mission Control Monitoring System (LMCMS) resident in the Operations Control Center. The LMCMS also provides the range simulation system the uplink command combined with differential GPS corrections. Later stages will require the progressive integration of actual range hardware with this simulation effort, leading to communication between telemetry, uplink and FTS antennas at NASA Dryden Flight Research Center, with vehicle antennas mounted on the Walter C. Williams Research Aircraft Integration Facility (RAIF). Decommutated Pulse Code Modulated (PCM) data is displayed on one of the four monitors that comprise the Range Safety Officer's (RSO) station. Also displayed are instantaneous impact prediction models, and Federal Aviation Administration (FAA) data for notification of other traffic in the area. Aside from initiating the flight termination command and validating communication links, the RSO station with the range simulation will be used to provide both range control and range safety officers training. The training is necessary to perform their respective functions with greater levels of confidence prior to first flight.

Sharma, Ashley↗