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 1,207 records · Page 67

Analysis of Inter-Moss Loops in the Solar Region with IRIS and SDO AIA: Automatic Event Detection and Characterization

The Interface Region Imaging Spectrograph (IRIS), launched in the summer of 2013, is designed specifically to observe and investigate the transition region and adjacent layers of the solar atmosphere, obtaining images in high spatial, temporal, and spectral resolution. Our particular work is focused on the evolution of inter-moss loops, which have been detected in the lower corona by the Atmospheric Imaging Assembly (AIA) and the High-Resolution Coronal Imager (Hi- C), but are known to have foot points below the transition region. With the high-resolution capabilities of IRIS and its Si IV pass band, which measures activity in the upper chromosphere, we can study these magnetic loops in detail and compare their characteristic length and time scales to those obtained from several AIA image sets, particularly the 171, 193, and 211 pass bands. By comparing the results between these four data sets, one can potentially establish a measure of the ionization equilibrium for the location in question. To explore this idea, we found a large, sit-and-stare observation within the IRIS database that fit our specifications. This data set contained a number of well-defined inter-moss loops (by visual inspection) with a cadence less than or equal to that of AIA (approximately 12 seconds). This particular data set was recorded on October 23, 2013 at 07:09:30, lasting for 3219 seconds with a field of view of 120.6 by 128.1 arcseconds, centered on -53.9 by 59.1 arcseconds from disk center. For ease of comparison, the AIA data has been interpolated to match the IRIS cadence and resolution. In the main portion of the poster, we demonstrate the detection of events, the information collected, and the immediate results to the right, showing the progress of an event with green as the start, blue as the peak, and red as the end. Below here, we demonstrate how pixels are combined to form groups. The 3D results are shown to the right.

event detection↗

Analysis of Inter-Moss Loops in the Solar Region with IRIS and SDO AIA: Automatic Event Detection and Characterization

The Interface Region Imaging Spectrograph (IRIS), launched in the summer of 2013, is designed specifically to observe and investigate the transition region and adjacent layers of the solar atmosphere, obtaining images in high spatial, temporal, and spectral resolution. Our particular work is focused on the evolution of inter-moss loops, which have been detected in the lower corona by the Atmospheric Imaging Assembly (AIA) and the High-Resolution Coronal Imager (Hi- C), but are known to have foot points below the transition region. With the high-resolution capabilities of IRIS and its Si IV pass band, which measures activity in the upper chromosphere, we can study these magnetic loops in detail and compare their characteristic length and time scales to those obtained from several AIA image sets, particularly the 171, 193, and 211 pass bands. By comparing the results between these four data sets, one can potentially establish a measure of the ionization equilibrium for the location in question. To explore this idea, we found a large, sit-and-stare observation within the IRIS database that fit our specifications. This data set contained a number of well-defined inter-moss loops (by visual inspection) with a cadence less than or equal to that of AIA (approximately 12 seconds). This particular data set was recorded on October 23, 2013 at 07:09:30, lasting for 3219 seconds with a field of view of 120.6 by 128.1 arcseconds, centered on -53.9 by 59.1 arcseconds from disk center. For ease of comparison, the AIA data has been interpolated to match the IRIS cadence and resolution. In the main portion of the poster, we demonstrate the detection of events, the information collected, and the immediate results to the right, showing the progress of an event with green as the start, blue as the peak, and red as the end. Below here, we demonstrate how pixels are combined to form groups. The 3D results are shown to the right

IRIS↗

Low Order Equivalent System Identification for the Tu-144LL Supersonic Transport Aircraft

Low order equivalent system models were identified from flight test data for the Tu-144LL supersonic transport aircraft. Flight test maneuvers were executed by Russian and American test Pilots flying the aircraft from Zhukovsky airfield outside Moscow, Russia. Flight tests included longitudinal and lateral / directional maneuvers at supersonic cruise flight conditions. Piloted frequency sweeps and multi-step maneuvers were used to generate data for closed loop low order equivalent system modeling. Model parameters were estimated using a flexible. high accuracy Fourier transform and an equation error / output error (EE/OE) formulation in the frequency domain. Results were compared to parameter estimates obtained using spectral estimation and subsequent least squares fit to frequency response data in Bode plots. Modeling results from the two methods agreed well for both a frequency sweep and multiple concatenated multi-step maneuvers. For a single multi-step maneuvers the EE/OE method gave a better model fit with improved prediction capability. A summary of closed loop low order equivalent system identification results for the Tu-144LL, including estimated parameters, standard errors, and flying qualities level predictions, were computed and tabulated.

Morelli, Eugene A.↗

High-degree Gravity Models from GRAIL Primary Mission Data

We have analyzed Ka‒band range rate (KBRR) and Deep Space Network (DSN) data from the Gravity Recovery and Interior Laboratory (GRAIL) primary mission (1 March to 29 May 2012) to derive gravity models of the Moon to degree 420, 540, and 660 in spherical harmonics. For these models, GRGM420A, GRGM540A, and GRGM660PRIM, a Kaula constraint was applied only beyond degree 330. Variance‒component estimation (VCE) was used to adjust the a priori weights and obtain a calibrated error covariance. The global root‒mean‒square error in the gravity anomalies computed from the error covariance to 320×320 is 0.77 mGal, compared to 29.0 mGal with the pre‒GRAIL model derived with the SELENE mission data, SGM150J, only to 140×140. The global correlations with the Lunar Orbiter Laser Altimeter‒derived topography are larger than 0.985 between l = 120 and 330. The free‒air gravity anomalies, especially over the lunar farside, display a dramatic increase in detail compared to the pre‒GRAIL models (SGM150J and LP150Q) and, through degree 320, are free of the orbit‒track‒related artifacts present in the earlier models. For GRAIL, we obtain an a posteriori fit to the S‒band DSN data of 0.13 mm/s. The a posteriori fits to the KBRR data range from 0.08 to 1.5 micrometers/s for GRGM420A and from 0.03 to 0.06 micrometers/s for GRGM660PRIM. Using the GRAIL data, we obtain solutions for the degree 2 Love numbers, k20=0.024615+/-0.0000914, k21=0.023915+/-0.0000132, and k22=0.024852+/-0.0000167, and a preliminary solution for the k30 Love number of k30=0.00734+/-0.0015, where the Love number error sigmas are those obtained with VCE.

Primary 1↗

A Comparative Study of Interferometric Regridding Algorithms

THe paper discusses regridding options: (1) The problem of interpolating data that is not sampled on a uniform grid, that is noisy, and contains gaps is a difficult problem. (2) Several interpolation algorithms have been implemented: (a) Nearest neighbor - Fast and easy but shows some artifacts in shaded relief images. (b) Simplical interpolator - uses plane going through three points containing point where interpolation is required. Reasonably fast and accurate. (c) Convolutional - uses a windowed Gaussian approximating the optimal prolate spheroidal weighting function for a specified bandwidth. (d) First or second order surface fitting - Uses the height data centered in a box about a given point and does a weighted least squares surface fit.

Hensley, Scott↗

JPL pulsar timing observations. I - The Vela pulsar

It is found through an analysis of 12 years of arrival time data for PSR 0833-45, by fitting a model of the pulse period to the data in a least-squares sense, that there exist regions of exceptionally stable behavior in which rms residuals are typically about 40 microsec. In addition, each residual exhibits a significant period parameter second derivative which is found to change little within a given interjump era. Parameter changes in the period P, and the derivatives introduced by each of the four large observed jumps in period, are recomputed to demonstrate that the post-jump decay is nearly exponential in form. It is found that random fluctuations in P are an order of magnitude less than the systematic variations caused by the seond derivative period parameter, whose component associated with the general slowdown of the pulsar can be reliably measured when the observations span 55 years or more.

Downs, G. S.↗

Gamma-ray burst high time-resolution spectral observations made with the Solar Maximum Mission

The gamma-ray bursts of April 19 and 21, 1980, and March 1, 1981, are characterized on the basis of observations obtained at 25-500 keV and time resolution 128 msec using the hard-X-ray-burst spectrometer of the SMM satellite. The data are presented graphically, and the parameters determined by fitting three theoretical models to the data for each phase (rising phase, decay phase, and valleys between pulses) are given in tables. Models used are thin thermal bremsstrahlung with or without absorption lines, thermal synchrotron with or without absorption lines, and power law with exponential cutoff.

Desai, U. D.↗

Low Reynolds number kappa-epsilon modeling with the aid of direct simulation data

The constant C(sub mu) and the near-wall damping function f(sub mu) in the eddy-viscosity relation of the kappa-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f(sub mu) model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.

Rodi, W.↗

Low Reynolds number k-epsilon modelling with the aid of direct simulation data

The constant C sub mu and the near-wall damping function f sub mu in the eddy-viscosity relation of the k-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f sub mu model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.

Rodi, W.↗

ROSAT data analysis with EXSAS

For the x-ray observatory ROSAT, data from survey and pointed mission phases taken with different focal plane instruments and according to a complex mission timeline have to be handled. Data analysis therefore puts high demands on appropriate software tools. With EXSAS - the Extended Scientific Analysis System developed with an effort of 20 man years by the German ROSAT Scientific Data Center - a comfortable system for the reduction of data from the ROSAT x-ray and XUV instruments has been made available. EXSAS comprises a large collection of application modules as typically required in analyzing data of this wavelength regime and runs as a specific context in the wide-spread ESO-MIDAS environment. EXSAS, completely written in FORTRAN 77, takes full advantage of all the standards used in MIDAS and therefore, reflects the same portability (different UNIX installations and VMS). If required, the FORTRAN code also enables users to adapt the software in an easy way to their specific needs. To maintain independence from the specifics of different operating systems also on the data input side, all ROSAT data redistributed in the widely accepted FITS format. Although EXSAS has been developed specifically for data analysis of the ROSAT instruments, its structural design is sufficiently general to serve equally well also data from other X-ray and XUV instruments. EXSAS analysis modules are grouped into 4 application packages dealing with Data Preparation and Instrument Correction, Spatial Analysis, Spectral Analysis and Timing Analysis. A special EXSAS header, read and updated by each application, maintains the general information transfer on the origin, the history and the parameter space of the data stored in tables and images. About 100 genuine commands (most of which offer several additional options) allow to interactively explore the functionality of the system. Up to now 40 institutes all over the world have requested the EXSAS software. Maintenance and regular updates of the software and the comprehensive documentation are provided by the ROSAT Scientific Data Center at Garching.

Zimmermann, H. U.↗

Opus: A Coordination Language for Multidisciplinary Applications

Data parallel languages, such as High Performance fortran, can be successfully applied to a wide range of numerical applications. However, many advanced scientific and engineering applications are multidisciplinary and heterogeneous in nature, and thus do not fit well into the data parallel paradigm. In this paper we present Opus, a language designed to fill this gap. The central concept of Opus is a mechanism called ShareD Abstractions (SDA). An SDA can be used as a computation server, i.e., a locus of computational activity, or as a data repository for sharing data between asynchronous tasks. SDAs can be internally data parallel, providing support for the integration of data and task parallelism as well as nested task parallelism. They can thus be used to express multidisciplinary applications in a natural and efficient way. In this paper we describe the features of the language through a series of examples and give an overview of the runtime support required to implement these concepts in parallel and distributed environments.

Chapman, Barbara↗

Cyclotron Resonant Scattering Feature Simulations II. Description of the CRSF Simulation Process

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons o_ quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems.Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time.Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua.Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗

Development of an Optical Model for the MEDLI2 Radiometer

We present an optical model of the MEDLI2 radiometer, which was flown onboard the Mars 2020 vehicle during atmospheric entry on February 18, 2021. The proposed model takes into account the radiometer geometry and optical phenomena such as reflection, refraction, and absorption to estimate the signal loss inherent in the design of the radiometer. The model serves as a filter, which can be used to convert predicted total radiative heat fluxes on the radiometer window into predicted measured heat fluxes, enabling more accurate comparisons of the radiometer flight data with predicted radiometer measurements. Optical properties of the radiometer’s sapphire window are inferred from spectral transmissivity measurements. A modified Lorentz Oscillator Model is used to represent the dielectric function of the window and is fit to the transmission data. The resulting optical model of the radiometer is validated against experimental data used to measure the dependence of the radiometer output voltage on the incident angle of the light impinging on it. Comparison of predicted radiometer measurements with the MEDLI2 flight data show good agreement between the predicted and measured peak heat fluxes, though offset at different times. A linear increase in the signal loss due to radiometer optics alone is also demonstrated, suggesting that the changing shock angle over the course of the trajectory impacts the signal loss.

MEDLI2↗

Low Order Equivalent System for the Tu-144LL Supersonic Transport Aircraft

Low order equivalent system models were identified from flight test data for the Tu- 144LL supersonic transport aircraft. Flight test maneuvers were executed by Russian and American test pilots flying the aircraft from Zhukovsky airfield outside Moscow, Russia. Flight tests included longitudinal and lateral/directional maneuvers at supersonic cruise flight conditions. Piloted frequency sweeps and multi-step maneuvers were used to Generate data for p closed loop low order equivalent system modeling Model parameters were estimated using a flexible, high accuracy Fourier transform and an equation error output error (EE/OE) formulation in the frequency domain. Results were compared to parameter estimates obtained using spectral estimation and subsequent least squares fit to frequency response data in Bode plots. Modeling results from the two methods a-reed well for both a frequency sweep and multiple concatenated multi-step maneuvers. For a single multi-step maneuver, the EE/OE method gave a better model fit with improved prediction capability. A summary of closed loop low order equivalent system identification results for the Tu-144LL. including estimated parameters, standard errors, and flying qualities level predictions, were computed and tabulated.

Morelli, Eugene A.↗

Parametrization of electron impact ionization cross sections for CO, CO2, NH3 and SO2

The electron impact ionization and dissociative ionization cross section data of CO, CO2, CH4, NH3, and SO2, measured in the laboratory, were parameterized utilizing an empirical formula based on the Born approximation. For this purpose an chi squared minimization technique was employed which provided an excellent fit to the experimental data.

Srivastava, Santosh K.↗

Stability Analysis of the Planetary System Orbiting Upsilon Andromedae

We present results of long-term numerical orbital integrations designed to test the stability of the three-planet system orbiting Upsilon Andromedae and short-term integrations to test whether mutual perturbations among the planets can be used to determine planetary masses. Our initial conditions are based on the latest fits to the radial velocity data obtained by the planet-search group at Lick Observatory. The new fits result in significantly more stable systems than did the initially announced planetary parameters. An analytic analysis of the star and the two outer planets shows that this subsystem is Hill stable up to five. Our integrations involving all three planets show that the system is stable for at least 100 Myr for up to four. In our simulations, we still see a secular resonance between the outer two planets and in some cases large oscillations in the eccentricity of the inner planet.

Lissauer, Jack J.↗

F-15B Quiet Spike(TradeMark) Aeroservoelastic Flight-Test Data Analysis

System identification is utilized in the aerospace community for development of simulation models for robust control law design. These models are often described as linear, time-invariant processes and assumed to be uniform throughout the flight envelope. Nevertheless, it is well known that the underlying process is inherently nonlinear. Over the past several decades the controls and biomedical communities have made great advances in developing tools for the identification of nonlin ear systems. In this report, we show the application of one such nonlinear system identification technique, structure detection, for the an alysis of Quiet Spike(TradeMark)(Gulfstream Aerospace Corporation, Savannah, Georgia) aeroservoelastic flight-test data. Structure detectio n is concerned with the selection of a subset of candidate terms that best describe the observed output. Structure computation as a tool fo r black-box modeling may be of critical importance for the development of robust, parsimonious models for the flight-test community. The ob jectives of this study are to demonstrate via analysis of Quiet Spike(TradeMark) aeroservoelastic flight-test data for several flight conditions that: linear models are inefficient for modelling aeroservoelast ic data, nonlinear identification provides a parsimonious model description whilst providing a high percent fit for cross-validated data an d the model structure and parameters vary as the flight condition is altered.

Kukreja, Sunil L.↗

Decadal-Scale Responses in Middle and Upper Stratospheric Ozone From SAGE II Version 7 Data

Stratospheric Aerosol and Gas Experiment (SAGE II) version 7 (v7) ozone profiles are analyzed for their decadal-scale responses in the middle and upper stratosphere for 1991 and 1992-2005 and compared with those from its previous version 6.2 (v6.2). Multiple linear regression (MLR) analysis is applied to time series of its ozone number density vs. altitude data for a range of latitudes and altitudes. The MLR models that are fit to the time series data include a periodic 11 yr term, and it is in-phase with that of the 11 yr, solar UV (Ultraviolet)-flux throughout most of the latitude/ altitude domain of the middle and upper stratosphere. Several regions that have a response that is not quite in-phase are interpreted as being affected by decadal-scale, dynamical forcings. The maximum minus minimum, solar cycle (SClike) responses for the ozone at the low latitudes are similar from the two SAGE II data versions and vary from about 5 to 2.5% from 35 to 50 km, although they are resolved better with v7. SAGE II v7 ozone is also analyzed for 1984-1998, in order to mitigate effects of end-point anomalies that bias its ozone in 1991 and the analyzed results for 1991-2005 or following the Pinatubo eruption. Its SC-like ozone response in the upper stratosphere is of the order of 4%for 1984-1998 vs. 2.5 to 3%for 1991-2005. The SAGE II v7 results are also recompared with the responses in ozone from the Halogen Occultation Experiment (HALOE) that are in terms of mixing ratio vs. pressure for 1991-2005 and then for late 1992- 2005 to avoid any effects following Pinatubo. Shapes of their respective response profiles agree very well for 1992-2005. The associated linear trends of the ozone are not as negative in 1992-2005 as in 1984-1998, in accord with a leveling off of the effects of reactive chlorine on ozone. It is concluded that the SAGE II v7 ozone yields SC-like ozone responses and trends that are of better quality than those from v6.2.

Remsberg, E. E.↗