Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Surface entropy of liquids via a direct Monte Carlo approach - Application to liquid Si

Two methods are presented for a direct Monte Carlo evaluation of the surface entropy S(s) of a liquid interacting by specified, volume-independent potentials. The first method is based on an application of the approach of Ferrenberg and Swendsen (1988, 1989) to Monte Carlo simulations at two different temperatures; it gives much more reliable results for S(s) in liquid Si than previous calculations based on numerical differentiation. The second method expresses the surface entropy directly as a canonical average at fixed temperature.

Wang, Z. Q.↗
NASA NTRS1991

Unsteady surface pressure measurements on a slender delta wing undergoing limit cycle wing rock

An experimental investigation of slender wing limit cycle motion known as wing rock was investigated using two unique experimental systems. Dynamic roll moment measurements and visualization data on the leading edge vortices were obtained using a free to roll apparatus that incorporates an airbearing spindle. In addition, both static and unsteady surface pressure data was measured on the top and bottom surfaces of the model. To obtain the unsteady surface pressure data a new computer controller drive system was developed to accurately reproduce the free to roll time history motions. The data from these experiments include, roll angle time histories, vortex trajectory data on the position of the vortices relative to the model's surface, and surface pressure measurements as a function of roll angle when the model is stationary or undergoing a wing rock motion. The roll time history data was numerically differentiated to determine the dynamic roll moment coefficient. An analysis of these data revealed that the primary mechanism for the limit cycle behavior was a time lag in the position of the vortices normal to the wing surface.

Arena, Andrew S., Jr.↗
NASA NTRS1991

Stability and bifurcation analysis of oscillators with piecewise-linear characteristics - A general approach

A general approach is developed for determining the periodic solutions and their stability of nonlinear oscillators with piecewise-smooth characteristics. A modified harmonic balance/Fourier transform procedure is devised for the analysis. The procedure avoids certain numerical differentiation employed previously in determining the periodic solutions, therefore enhancing the reliability and efficiency of the method. Stability of the solutions is determined via perturbations of their state variables. The method is applied to a forced oscillator interacting with a stop of finite stiffness. Flip and fold bifurcations are found to occur. This led to the identification of parameter ranges in which chaotic response occurred.

Noah, S. T.↗
NASA NTRS1992

Long-term variations in abundance and distribution of sulfuric acid vapor in the Venus atmosphere inferred from Pioneer Venus and Magellan radio occultation studies

Radio occultation experiments have been used to study various properties of planetary atmospheres, including pressure and temperature profiles, and the abundance profiles of absorbing constituents in those planetary atmospheres. However, the reduction of amplitude data from such experiments to determine abundance profiles requires the application of the inverse Abel transform (IAT) and numerical differentiation of experimental data. These two operations preferentially amplify measurement errors above the true signal underlying the data. A new technique for processing radio occultation data has been developed that greatly reduces the errors in the derived absorptivity and abundance profiles. This technique has been applied to datasets acquired from Pioneer Venus Orbiter radio occultation studies and more recently to experiments conducted with the Magellan spacecraft. While primarily designed for radar studies of the Venus surface, the high radiated power (EIRP) from the Magellan spacecraft makes it an ideal transmitter for measuring the refractivity and absorptivity of the Venus atmosphere by such experiments. The longevity of the Pioneer Venus Orbiter has made it possible to study long-term changes in the abundance and distribution of sulfuric acid vapor, H2SO4(g), in the Venus atmosphere between 1979 and 1992. The abundance of H2SO4(g) can be inferred from vertical profiles of 13-cm absorptivity profiles retrieved from radio occultation experiments. Data from 1979 and 1986-87 suggest that the abundance of H2SO4(g) at latitudes northward of 70 deg decreased over this time period. This change may be due to a period of active volcanism in the late 1970s followed by a relative quiescent period, or some other dynamic process in the Venus atmosphere. While the cause is not certain, such changes must be incorporated into dynamic models of the Venus atmosphere. Potentially, the Magellan spacecraft will extend the results of Pioneer Venus Orbiter and allow the continued monitoring of the abundance of distribution of H2SO4(g) in the Venus atmosphere, as well as other interesting atmospheric properties. Without such measurements it will be difficult to address other issues such as the short-term spatial variability of the abundance of H2SO4(g) at similar latitudes in Venus atmosphere, and the identities of particles responsible for large-scale variations observed in NIR images.

Jenkins, J. M.↗
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 NTRS2001

A Leo Satellite Navigation Algorithm Based on GPS and Magnetometer Data

The Global Positioning System (GPS) has become a standard method for low cost onboard satellite orbit determination. The use of a GPS receiver as an attitude and rate sensor has also been developed in the recent past. Additionally, focus has been given to attitude and orbit estimation using the magnetometer, a low cost, reliable sensor. Combining measurements from both GPS and a magnetometer can provide a robust navigation system that takes advantage of the estimation qualities of both measurements. Ultimately, a low cost, accurate navigation system can result, potentially eliminating the need for more costly sensors, including gyroscopes. This work presents the development of a technique to eliminate numerical differentiation of the GPS phase measurements and also compares the use of one versus two GPS satellites.

Deutschmann, Julie↗
NASA NTRS2003

In-Flight MTF Analysis of ADAR 5500 Aircraft Sensor

The Modulation Transfer Function (MTF) of an imaging system is one parameter that can be used to describe the spatial resolution or the image quality of an imaging system. MTF is the system's frequency response to an input. The MTF at a specific frequency normally ranges from 0-1, where 0 indicates no frequency response and 1 indicates perfect frequency response. The MTF value at Nyquist frequency is calculated because the Nyquist frequency is the maximum sampling frequency of the system. MTF analysis was performed on data acquired with an ADAR 5500 flown by Positive Systems on September 13-14, 2000 over Brookings, SD. Four edges from ground targets were used as inputs to the system. The system's MTF was calculated using a Matlab-based algorithm, which applied edge detection, numerical differentiation, and Fourier transformation methods to the target images.

Zanoni, Vicki↗
NASA NTRS2001

Estimation of diastolic intraventricular pressure gradients by Doppler M-mode echocardiography

Previous studies have shown that small intraventricular pressure gradients (IVPG) are important for efficient filling of the left ventricle (LV) and as a sensitive marker for ischemia. Unfortunately, there has previously been no way of measuring these noninvasively, severely limiting their research and clinical utility. Color Doppler M-mode (CMM) echocardiography provides a spatiotemporal velocity distribution along the inflow tract throughout diastole, which we hypothesized would allow direct estimation of IVPG by using the Euler equation. Digital CMM images, obtained simultaneously with intracardiac pressure waveforms in six dogs, were processed by numerical differentiation for the Euler equation, then integrated to estimate IVPG and the total (left atrial to left ventricular apex) pressure drop. CMM-derived estimates agreed well with invasive measurements (IVPG: y = 0.87x + 0.22, r = 0.96, P < 0.001, standard error of the estimate = 0.35 mmHg). Quantitative processing of CMM data allows accurate estimation of IVPG and tracking of changes induced by beta-adrenergic stimulation. This novel approach provides unique information on LV filling dynamics in an entirely noninvasive way that has previously not been available for assessment of diastolic filling and function.

Non-NASA Center↗
NASA NTRS2005

Cross Section Sensitivity and Propagated Errors in HZE Exposures

It has long been recognized that galactic cosmic rays are of such high energy that they tend to pass through available shielding materials resulting in exposure of astronauts and equipment within space vehicles and habitats. Any protection provided by shielding materials result not so much from stopping such particles but by changing their physical character in interaction with shielding material nuclei forming, hopefully, less dangerous species. Clearly, the fidelity of the nuclear cross-sections is essential to correct specification of shield design and sensitivity to cross-section error is important in guiding experimental validation of cross-section models and database. We examine the Boltzmann transport equation which is used to calculate dose equivalent during solar minimum, with units (cSv/yr), associated with various depths of shielding materials. The dose equivalent is a weighted sum of contributions from neutrons, protons, light ions, medium ions and heavy ions. We investigate the sensitivity of dose equivalent calculations due to errors in nuclear fragmentation cross-sections. We do this error analysis for all possible projectile-fragment combinations (14,365 such combinations) to estimate the sensitivity of the shielding calculations to errors in the nuclear fragmentation cross-sections. Numerical differentiation with respect to the cross-sections will be evaluated in a broad class of materials including polyethylene, aluminum and copper. We will identify the most important cross-sections for further experimental study and evaluate their impact on propagated errors in shielding estimates.

Heinbockel, John H.↗
NASA NTRS2018

Implementation and Use of a Convergence Test for Iterative Wind Tunnel Balance Load Predictions

A modified and extended version of a convergence test for wind tunnel strain-gage balance load iterations was implemented. The test uses an upper bound of the Lipschitz constant to assess convergence characteristics of balance load predictions if the Iterative Method is applied. Convergence is expected within the use envelope of the balance whenever this upper bound is less than the threshold of one. It is explained in great detail how the convergence test can be applied to the two load iteration equation types that are currently being used in the aerospace testing community. In addition, the application of the test to balances with bi-directional output characteristics is discussed. It is also shown how numerical differentiation can be used in order to obtain partial derivatives that are needed for the calculation of the Lipschitz constant. Finally, machine calibration data of NASA's MC60E six-component force balance is selected to demonstrate both implementation and use of the convergence test.

calibration analysis↗
NASA NTRSDate not supplied

Updates and Correlation of EMU System-Level Model (SINDA EMU)

During United States Extravehicular Activity 80 (US EVA 80), water was observed in the helmet of an Extravehicular Mobility Unit (EMU) during cabin repressurization. One of the primary methods of determining the likely cause of this failure was through a comparison of EVA 80 to other historical EVAs using an analytical approach. The Systems Improved Numerical Differential Analysis EMU model (SINDA EMU) is a system-level model of the EMU that was used in this investigation. SINDA EMU was initially developed and correlated to test data in the 1980s. Since its conception, SINDA EMU has been continually adjusted based on new test data and changes to the EMU design. To support the water in the helmet investigation, SINDA EMU needed to be further updated and recorrelated to ensure accurate results. These changes included changing the primary carbon dioxide (CO2) removal technology, implementing logic to allow for re-evaporation of sweat runoff from the liquid cooling and ventilation garment (LCVG), and improving the transient modeling capabilities. To validate the implementation of these adjustments, SINDA EMU was correlated to test data from the 1990s and human-in-the-loop (HITL) testing from 2014. These updates and correlation efforts proved that SINDA EMU is an effective tool for investigating the EVA 80 water in the helmet failure event.

Noah Lial Andersen↗
NASA NTRS1987

Parallel Algorithm Solves Coupled Differential Equations

Numerical methods adapted to concurrent processing. Algorithm solves set of coupled partial differential equations by numerical integration. Adapted to run on hypercube computer, algorithm separates problem into smaller problems solved concurrently. Increase in computing speed with concurrent processing over that achievable with conventional sequential processing appreciable, especially for large problems.

Hayashi, A.↗