Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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 667 records · Page 37

Model Evaluation of Methods for Estimating Surface Emissions and Chemical Lifetimes from Satellite Data

Column densities from satellite retrievals can provide valuable information for estimating emissions and chemical lifetimes objectively across the globe. To better understand the uncertainties associated with these estimates, we test four methods using simulated column densities from a point source: a box model approach, a 2D Gaussian fit, an Inverse Radius fit and an Exponentially-Modified Gaussian fit. The model results were simulated using the WRF and CAMx models for the year 2005, for a single point source outside Atlanta in Georgia, USA with specified emissions and three chemical scenarios: no chemical reactions, 12 h chemical lifetime and 1 h chemical lifetime. No other sources were included in the simulations. We find that the box model provides reliable estimates irrespective of plume speed and plume direction, if the plume speed and the chemical lifetime are known accurately. The 2D Gaussian fit was found to be sensitive to plume speed and direction, and requires omnidirectional dispersion in order to have a decent fit. However, the 2D Gaussian fit is only an approximate fit to the data, and the discrepancies mean that the results are dependent on the geographical domain used for the optimization. An Inverse Radius fit is introduced to correct this issue, which is found to provide improved emissions and lifetime estimates. The Exponentially-Modified Gaussian fit also gave improved estimates. It is however dependent on accurate plume rotation such that reported chemical lifetimes with this method could be significantly underestimated.

Emission inventory↗

Reconstruction of the 1st Space Shuttle (STS-1) entry trajectory

A discussion of the generation of the best estimate trajectory (BET) of the first Space Shuttle Orbiter entry flight is presented. The BET defines a time history of the state, attitude, and atmospheric relative parameters throughout the Shuttle entry from an altitude of approximately 183 km to rollout. The inertial parameters were estimated utilizing a weighted least squares batch filter algorithm. Spacecraft angular rate and acceleration data derived from the Inertial Measurement Unit were utilized to predict the state and attitude which was constrained in a weighted least squares process to fit external tracking data consisting of ground based S-band and C-band data. Refined spacecraft altitude and velocity during and post rollout were obtained by processing artificial altimeter and Doppler data. The BET generation process is discussed. Software and data interface discussions are included. The variables and coordinate systems utilized are defined. STS-1 mission peculiar inputs are summarized. A listing of the contents of the actual BET is provided.

Findlay, J. T.↗

Contouring Trivariate Data

In many applications, the data consist of discrete 3-D points at which one or more parameters are given. To display contours, the data are represented by a continuous function which is evaluated at any point as needed for contouring. Contouring results are presented which are applicable both to artitrarily spaced data and to data which lie on a topologically rectangular three dimensional grid. the contours are assumed to be described by 3-D display lists for viewing on a dynamic color graphics device; that is, they will not simply be projected into 2-D and viewed as a static image on a frame buffer. Dynamic viewing of color contours may be essential to the proper interpretation of results. The gridded data lie on a topologically rectangular grid although two or more nodes may be the same point. Parametric tensor product methods may be used to fit the gridded data and to generate the contours. Rectangular elements are convenient but not necessary. For example, there are other methods which are effective for contouring over tetrahedral elements.

Stead, S. E.↗

Automation of reverse engineering process in aircraft modeling and related optimization problems

During the year of 1994, the engineering problems in aircraft modeling were studied. The initial concern was to obtain a surface model with desirable geometric characteristics. Much of the effort during the first half of the year was to find an efficient way of solving a computationally difficult optimization model. Since the smoothing technique in the proposal 'Surface Modeling and Optimization Studies of Aerodynamic Configurations' requires solutions of a sequence of large-scale quadratic programming problems, it is important to design algorithms that can solve each quadratic program in a few interactions. This research led to three papers by Dr. W. Li, which were submitted to SIAM Journal on Optimization and Mathematical Programming. Two of these papers have been accepted for publication. Even though significant progress has been made during this phase of research and computation times was reduced from 30 min. to 2 min. for a sample problem, it was not good enough for on-line processing of digitized data points. After discussion with Dr. Robert E. Smith Jr., it was decided not to enforce shape constraints in order in order to simplify the model. As a consequence, P. Dierckx's nonparametric spline fitting approach was adopted, where one has only one control parameter for the fitting process - the error tolerance. At the same time the surface modeling software developed by Imageware was tested. Research indicated a substantially improved fitting of digitalized data points can be achieved if a proper parameterization of the spline surface is chosen. A winning strategy is to incorporate Dierckx's surface fitting with a natural parameterization for aircraft parts. The report consists of 4 chapters. Chapter 1 provides an overview of reverse engineering related to aircraft modeling and some preliminary findings of the effort in the second half of the year. Chapters 2-4 are the research results by Dr. W. Li on penalty functions and conjugate gradient methods for quadratic programming problems.

Li, W.↗

Relativistic Iron Emission and Disk Reflection in Galactic Microquasar XTE J1748-288

We report evidence for an Fe K-alpha fluorescence line feature and disk reflection in the very high, high-, and low-state X-ray spectra of the Galactic microquasar XTE J1748 - 288 during its 1998 June outburst. Spectral analyses are made on data gathered throughout the outburst by the Rossi X-Ray Timing Explorer Proportional Counter Array. Gaussian line, relativistic disk emission line, and ionized disk reflection models are fitted to the data. In the very high state the line profile appears strongly redshifted, consistent with disk emission from the innermost stable orbits around a maximally rotating Kerr black hole. In the high state the line profile is less redshifted and increasingly prominent. The low-state line profile is very strong (approx. 0.5 keV equivalent width) and centered at 6.7 +/- 0.10 keV; disk line emission model fits indicate that the inner edge of the disk fluctuates between approx. 20R(sub g) and - approx. 100R(sub g) in this state. The disk reflection fraction is traced through the outburst; reflection from an ionized disk is preferred in the very high and high states, and reflection from a relatively neutral disk is preferred in the low state. We discuss the implications of our findings for the binary system dynamics and accretion flow geometry in XTE J1748 - 288.

Miller, J. M.↗

Relativistic Iron Emission and Disk Reflection in Galactic Microquasar XTE J1748-288

We report evidence for an Fe K(alpha) fluorescence line feature and disk reflection in the very high, high-, and low-state X-ray spectra of the Galactic microquasar XTE J1748-288 during its 1998 June outburst. Spectral analyses are made on data gathered throughout the outburst by the Rossi X-Ray Timing Explorer Proportional Counter Array. Gaussian line, relativistic disk emission line, and ionized disk reflection models are fitted to the data. In the very high state the line profile appears strongly redshifted, consistent with disk emission from the innermost stable orbits around a maximally rotating Kerr black hole. In the high state the line profile is less redshifted and increasingly prominent. The low-state line profile is very strong (approx. 0.5 keV equivalent width) and centered at 6.7 +/- 0.10 keV; disk line emission model fits indicate that the inner edge of the disk fluctuates between approx. 20Rg and approx. 100Rg in this state. The disk reflection fraction is traced through the outburst; reflection from an ionized disk is preferred in the very high and high states, and reflection from a relatively neutral disk is preferred in the low state. We discuss the implications of our findings for the binary system dynamics and accretion flow geometry in XTE J1748-288.

Miller, J. M.↗

Spectra and positions of galactic gamma-ray sources

The UCSD/MIT Hard X-Ray and Low Energy Gamma-Ray Experiment aboard HEAO-1 scanned the galactic center region during three epochs in 1977 and 1978 from 13 to 180 keV. The results are presented from the scanning epoch of 1978 September. Twenty-two known 2 to 10 keV source positions were necessary for an acceptable fit to the data. The spectra of the 16 strongest, least confused sources are all consistent with power laws with photon spectral indices ranging from 2.1 to 7.2. Acceptable fits to thermal bremsstrahlung models are also possible for most sources. No one source in this survey can be extrapolated to higher energy to match the intensity of the gamma-ray continuum as measured by HEAO-1 large field of view detectors, which implies that the continuum is a composite of contributions from a number of sources.

Knight, F. K.↗

From EXOSAT to the High Energy Astrophysics Science Archive (HEASARC): X-ray Astronomy Comes of Age

In May 1983 the European Space Agency launched EXOSAT, its first X-ray astronomy observatory. Even though it lasted only 3 short years, this mission brought not only new capabilities that resulted in unexpected discoveries, but also a pioneering approach to operations and archiving that changed X-ray astronomy from observations led by small instrument teams, to an observatory approach open to the entire community through a guest observer program. The community use of the observatory was supported by a small dedicated team of scientists, the precursor to the data center activities created to support e.g. Chandra and XMM-Newton. The new science capabilities of EX OS AT included a 90 hr highly eccentric high earth orbit that allow unprecedented continuous coverage of sources as well as direct communication with the satellite that allowed real time decisions to respond to unexpected events through targets of opportunity. The advantages of this orbit demonstrated by EXOSAT resulted in Chandra and XMM-Newton selecting similar orbits. The three instruments on board the EXOSAT observatory were complementary, designed to give complete coverage over a wide energy band pass of 0.05-50 keY. An onboard processor could be programmed to give multiple data modes that could be optimized in response to science discoveries: These new capabilities resulted in many new discoveries including the first comprehensive study of AGN variability, new orbital periods in X-ray binaries and cataclysmic variables, new black holes, quasi-periodic oscillations from neutron stars and black holes and broad band X-ray spectroscopy. The EXOSAT team generated a well-organized database accessible worldwide over the nascent internet, allowing remote selection of data products, making samples and undertaking surveys from the data. The HEASARC was established by NASA at Goddard Space Flight Center in 1990 as the repository of NASA X-ray and Gamma-ray data. The proven EXOSAT database system became the core of the HEASARC infrastructure. The HEASARC pioneered many concepts now taken for granted including standardized formats using FITS files, restoring data from earlier missions, multi-mission analysis tools and a searchable archive over the world wide web.

White, Nicholas E.↗

Elliptical storm cell modeling of digital radar data

A model for spatial distributions of reflectivity in storm cells was fitted to digital radar data. The data were taken with a modified WSR-57 weather radar with 2.6-km resolution. The data consisted of modified B-scan records on magnetic tape of storm cells tracked at 0 deg elevation for several hours. The MIT L-band radar with 0.8-km resolution produced cross-section data on several cells at 1/2 deg elevation intervals. The model developed uses ellipses for contours of constant effective-reflectivity factor Z with constant orientation and eccentricity within a horizontal cell cross section at a given time and elevation. The centers of the ellipses are assumed to be uniformly spaced on a straight line, with areas linearly related to log Z. All cross sections are similar at different heights (except for cell tops, bottoms, and splitting cells), especially for the highest reflectivities; wind shear causes some translation and rotation between levels. Goodness-of-fit measures and parameters of interest for 204 ellipses are considered.

Altman, F. J.↗

Poisson vs. Gaussian statistics for sparse X-ray data: Application to the soft X-ray spectrometer

Reliable results when fitting X-ray data require proper consideration of the statistics involved. We probe the impact of Gaussian versus Poisson statistics at low count levels using both the standard χ^(2) method and maximum likelihood based on Poisson (C) statistics. The difference is studied and quantified through simulated spectra with known properties. We then test the results through analysis of Mn Kα calibration data taken with the flight spare microcalorimeter for the Hitomi soft X-ray spectrometer. Through comparison with simulations, our results show that the χ^(2) method tends to give overly optimistic estimates of the detector energy resolution, in particular when there are few counts. Given an energy resolution of ∼5 eV and a line with about 100 photons, the line width becomes ∼10% lower in the χ^(2) method than in Poisson statistics. This is a consequence of the uncertainties being dominated by counting statistics, and therefore highlights the need to choose the appropriate fit statistic.

Shinya Yamada↗

X-ray spectra of 12 distant clusters of galaxies

An X-ray spectral analysis of 12 distant clusters of galaxies observed on-axis with the Einstein Observatory imaging proportional counter is presented. Statistically, X-ray spectral data of the highest red shift clusters in the sample are inconsistent with an optically thin plasma emission model plus absorptions only in the Milky Way. Inclusions of excess absorptions in the clusters significantly improve the spectral fits to the data, indicating the possible presence of large amounts of X-ray absorbing cool gas in some of the distant clusters. By using X-ray luminosities determined from observed fluxes only in the 0.8-3.5 keV band and the temperature-to-luminosity correlation for nearby clusters, the estimated temperatures of the hot intracluster medium in these distant clusters constrain the absorption columns in the clusters. The data also show possible positive temperature evolution of the hot ICM with time which, if present, would further increase the absorptions required in fitting the spectra.

Wang, Qingde↗

MODIS Reflective Solar Band Calibration Improvements using Pseudo-Invariant Desert Targets

To provide the best science data quality, an accurate characterization of the response versus scan angle (RVS) function is critical for the MODIS reflective solar bands (RSB) on-orbit calibration. In every MODIS operational scan, the Earth’s surface, referred to here as Earth view (EV), the space view (SV) port, and the onboard calibrators are viewed via a two-sided scan mirror. The RVS is defined as the sensor’s relative response as a function the angle of incidence (AOI) to the scan mirror. Many different approaches have been developed to derive the time-dependent RVS and its look-up table (LUT) applied to MODIS Level 1B (L1B) products since calibration Collection 4. For most MODIS RSB, the on-board calibrators can reasonably track the RVS change with time. In practice, their RVS is derived using data from on-board calibrators and the EV mirror side ratio (for mirror side 2). For Terra bands 1-4, 8-10 and Aqua bands 1-4, 8-9, an enhancement has been employed in Collections 6 and 6.1 (C6/C6.1) by using Earth scene response trending from pseudo-invariant desert sites in addition to the onboard calibrators. The current C6/C6/1 RVS algorithm is focused on fitting the EV data at each AOI over time and then deriving the relative change at different AOI. The EV response trending is currently fitted with multiple segments over time. Alternatively, the EV responses can be fit first as a function of AOI before fitting temporally in order to reduce the dependence on the stability of the desert site. These pre-treatment methods on the EV data provide improvement in the derived calibration coefficients. However, evidence of insufficient calibration is still observed in the MODIS L1B reflectance data, especially in the form of differences between the mirror sides. In this paper, we review the current methodologies that utilize the EV response trends from the pseudo-invariant Libyan desert targets to supplement the gain derived from the onboard calibrators. An improvement is then proposed and investigated such that a sliding window average (SWA) is used to pre-process the raw EV data. The SWA parameters are carefully selected using trade-off studies to accurately track the Earth scene response trending in multiple cases to overcome the reflectance differences between two mirror sides. Calibration results show improvements for both Aqua and Terra MODIS RSB L1B data products. This new adjustment has been included in the recently delivered Collection 7 LUT that will be evident in the L1B products expected to be released in late 2021.

MODIS↗

In Pursuit of the FIP Effect in Stellar Coronae

Observational data have been reduced. Sig Gem data have been fitted using the XSPEC software, though the treatment within this software of the Fe transitions arising from the n=2 levels is thought to be sufficiently inaccurate as to affect the results obtained. Data for the same star obtained by EUVE have been extracted from the archive and are being analyzed in order to try and constrain the coronal DEM. Results of the first sets of model fitting have been presented at an AAS meeting; a more substantial presentation of these and updated results in a main-stream referred journal is in progress.

Drake, Jeremy J.↗

Dynamic model verification of large structural systems

The objective of the present methodology is two-fold: (1) to process test data obtained from either modal survey tests, or slow sine-sweep tests, to extract a set of orthogonal modes best matching the test data while being commensurate with the dynamic model, and (2) to modify submatrices of the dynamic model mass and stiffness matrices to adjust the model to best fit the test data. The method has been implemented using a linear statistical sequential estimator for computation on a CDC computer. Demonstration problems involving Space Shuttle quarter-scale vibration test data and dynamic models have been run. This paper will discuss the general methodology and experience to date.

Lee, L. T.↗

Solar O I 1304-A triplet line profiles

Estimates were made of the full-disk line profiles of the solar O I resonant triplet at 1304 A, using data from the SMM and the Orbiting Solar Observatory 8 for the time period 1975-1985. The observed line profiles are parameterized using a nonlinear least squares fit of the data by a simple empirical representation of the general line shape. The derived full-disk line shape parameters may be used, in combination with an appropriate value for the line-integrated full-disk solar flux, in the analysis of planetary or cometary O I 1304-A observations under any solar conditions.

Gladstone, G. R.↗

Using Decision Trees to Detect and Isolate Simulated Leaks in the J-2X Rocket Engine

The goal of this work was to use data-driven methods to automatically detect and isolate faults in the J-2X rocket engine. It was decided to use decision trees, since they tend to be easier to interpret than other data-driven methods. The decision tree algorithm automatically "learns" a decision tree by performing a search through the space of possible decision trees to find one that fits the training data. The particular decision tree algorithm used is known as C4.5. Simulated J-2X data from a high-fidelity simulator developed at Pratt & Whitney Rocketdyne and known as the Detailed Real-Time Model (DRTM) was used to "train" and test the decision tree. Fifty-six DRTM simulations were performed for this purpose, with different leak sizes, different leak locations, and different times of leak onset. To make the simulations as realistic as possible, they included simulated sensor noise, and included a gradual degradation in both fuel and oxidizer turbine efficiency. A decision tree was trained using 11 of these simulations, and tested using the remaining 45 simulations. In the training phase, the C4.5 algorithm was provided with labeled examples of data from nominal operation and data including leaks in each leak location. From the data, it "learned" a decision tree that can classify unseen data as having no leak or having a leak in one of the five leak locations. In the test phase, the decision tree produced very low false alarm rates and low missed detection rates on the unseen data. It had very good fault isolation rates for three of the five simulated leak locations, but it tended to confuse the remaining two locations, perhaps because a large leak at one of these two locations can look very similar to a small leak at the other location.

Schwabacher, Mark A.↗

Non-Contact Thrust Stand Calibration Method for Repetitively-Pulsed Electric Thrusters

A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.

Wong, Andrea R.↗

GOES-17 Magnetometer On-Orbit Calibration Design and Results

The Geostationary Operational Environmental Satellite-R (GOES-R) is a series of four NOAA weather satellites. Each spacecraft flies a Magnetometer instrument which includes two boom mounted fluxgate sensors located 6.3 and 8.5 meters from the spacecraft. Sensor readings are corrupted by both internal (bias and scale factor) and external (spacecraft emitted fields and misalignments) error sources. In order to estimate the bias, an on-orbit calibration is performed, slewing the spacecraft about multiple axes and then performing a fit of the data to determine zero offset and misalignment terms. The first satellite in the series, GOES-16 was launched in November 2016. Its magnetometer calibration maneuver was found to be inadequate providing insufficient rotation about orthogonal axes. The uncertainty in the estimated biases was found to be in the 3-5 nanotesla range, exceeding the magnetometer accuracy requirement of 1.0 nanotesla. Based on lessons learned, the maneuver was redesigned to provide multiple 360 degree rotations about two orthogonal axes, providing optimal visibility into bias and misalignment. Varying ambient fields over the duration of the calibration maneuvers, corrupts the ability to accurately determine parameters of interest. To address this issue, rather than assuming a constant field for the duration of the maneuver, a non-linear least squares filter was employed that uses a series of spline functions to better fit the ambient field. The number of splines can be varied to minimize the least squares filter residuals. Simulations were performed using 360 days of GOES-16 data to bound the uncertainty in the biases for given residuals of the fit to the ambient field. Following the launch of the second satellite in the series, GOES-17 in March 2018, the new calibration maneuver and enhanced bias estimation algorithm were exercised twice. Based on the residuals of the fit, the uncertainty of the estimated zero offset was reduced by an order of magnitude from 3-5 nanotesla to approximately 0.4 nanotesla per axis. This paper provides details of the maneuver design and the innovative algorithm used to estimate the GOES-17 biases and misalignments.

GOES-17↗