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,243 records · Page 69

Atmospheric Characterization and Further Orbital Modeling of κ Andromeda b

We present κ Andromeda b's photometry and astrometry taken with Subaru/SCExAO+HiCIAO and Keck/NIRC2, combined with recently published SCExAO/CHARIS low-resolution spectroscopy and published thermal infrared photometry to further constrain the companion's atmospheric properties and orbit. The Y/Y−K colors of κ And b are redder than field dwarfs, consistent with its youth and lower gravity. Empirical comparisons of its Y-band photometry and CHARIS spectrum to a large spectral library of isolated field dwarfs reaffirm the conclusion from Currie et al. that it likely has a low gravity but admit a wider range of most plausible spectral types (L0–L2). Our gravitational classification also suggests that the best-fit objects for κ And b may have lower gravity than those previously reported. Atmospheric models lacking dust/clouds fail to reproduce its entire 1–4.7 μm spectral energy distribution (SED), and cloudy atmosphere models with temperatures of ~1700–2000 K better match κ And b data. Most well-fitting model comparisons favor 1700–1900 K, a surface gravity of log(g) ~ 4–4.5, and a radius of 1.3–1.6 R (sub Jup); the best-fit model (Drift-Phoenix) yields the coolest and lowest-gravity values: T (sub eff) = 1700 K and log g = 4.0. An update to κ And b's orbit with ExoSOFT using new astrometry spanning 7 yr reaffirms its high eccentricity (0.77 ± 0.08). We consider a scenario where unseen companions are responsible for scattering κ And b to a wide separation and high eccentricity. If three planets, including κ And b, were born with coplanar orbits, and one of them was ejected by gravitational scattering, a potential inner companion with mass ≲10 M (sub Jup) could be located at ≲25 au.

Exoplanet astronomy↗

Columbia River volcanism - The question of mantle heterogeneity or crustal contamination

It is found that, although crystal fractionation played an important role in producing the chemical characteristics of Columbia River Province basalts displaying a wide range of chemical and isotopic compositions, the isotopic variability calls for the involvement of at least two isotopically distinct components. The major and trace element characteristics of the main volume of the basalts are not consistent with a metasomatized mantle source region, and the presence of a primordial mantel component is not supported by the chemical data. Models of simple binary mixing between a primary magma and Precambrian sialic crustal materials, while satisfying the observed Nd and Sr isotopic variations, fail to account for major trace element abundances. The combination of crustal assimilation and fractional crystalization is found to give a superior fit to the compositional data.

Carlson, R. W.↗

Estimates of M2 Tidal Energy Dissipation from TOPEX/Poseidon Altimeter Data

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

Egbert, Gary D.↗

Thermal Capacitance (Slug) Calorimeter Theory Including Heat Losses and Other Decaying Processes

A mathematical model, termed the Slug Loss Model, has been developed for describing thermal capacitance (slug) calorimeter behavior when heat losses and other decaying processes are not negligible. This model results in the temperature time slope taking the mathematical form of exponential decay. When data is found to fit well to this model, it allows a heat flux value to be calculated that corrects for the losses and may be a better estimate of the cold wall fully catalytic heat flux, as is desired in arc jet testing. The model was applied to the data from a copper slug calorimeter inserted during a particularly severe high heating rate arc jet run to illustrate its use. The Slug Loss Model gave a cold wall heat flux 15% higher than the value of 2,250 W/sq cm obtained from the conventional approach to processing the data (where no correction is made for losses). For comparison, a Finite Element Analysis (FEA) model was created and applied to the same data, where conduction heat losses from the slug were simulated. The heat flux determined by the FEA model was found to be in close agreement with the heat flux determined by the Slug Loss Model.

Hightower, T. Mark↗

Adaptive Pressure Profile Method to Locate the Isolator Shock Train Leading Edge Given Limited Pressure Information

To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.

Robin L Hunt↗

Adaptive Pressure Profile Method to Locate the Isolator Shock Train Leading Edge Given Limited Pressure Information

To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.

Robin Hunt↗

The Atacama Cosmology Telescope: DR4 maps and cosmological parameters

We present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013{2016 at 98 and 150 GHz. The maps cover more than 17,000 deg(^2), the deepest 600 deg(^2) with noise levels below 10µK-arcmin. We use the power spectrum derived from almost 6,000 deg(^2) of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, H0. By combining ACT data with large-scale information from WMAP we measure H0 = 67:6±1:1 km/s/Mpc, at 68% confidence, in excellent agreement with the independently measured Planck satellite estimate (from ACT alone we find H0 = 67:9± 1:5 km/s/Mpc). The ΛCDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1σ; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with ΛCDM predictions to within 1.5{2.2σ. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis.

Simone Aiola↗

3 keV to 2 MeV observations of four gamma-ray bursts

Simultaneous X and gamma ray spectra were obtained for four gamma ray burst events using satellite instrumentation. The P78-1 satellite proportional counters recorded X-ray data in the 3-10 keV range, while Pioneer Venus Orbiter and ISEE 3 sensors had a 3 keV lower limit for events. The positively correlated data covered the events GB 790307, 790325, 790504, and 790731. The time histories of the X and gamma ray energies did not correlate well. The post-gamma ray phase of GB 790307 displayed features of simple cooling, and in conjunction with soft X-ray data suggested an origin in an optically thin fireball. However, the data were equally well-fitted by a model of an emitting hot spot on a neutron star. The data negated any concept that gamma ray bursts evolve similarly to classical X-ray bursts.

Laros, J. G.↗

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

System identification or mathematical modelling is utilised in the aerospace community for the 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. The reason for utilising a linear approach has been due to the lack of a proper set of tools for the identification of nonlinear systems. Over the past several decades the controls and biomedical communities have made great advances in developing tools for the identification of nonlinear systems. These approaches are robust and readily applicable to aerospace systems. In this paper, we show the application of one such nonlinear system identification technique, structure detection, for the analysis of F-15B QuietSpike(TradeMark) aeroservoelastic flight test data. Structure detection is concerned with the selection of a subset of candidate terms that best describe the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modelling may be of critical importance for the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion which may save significant development time and costs. The objectives of this study are to demonstrate via analysis of F-15B QuietSpike(TradeMark) aeroservoelastic flight test data for several flight conditions (Mach number) that (i) linear models are inefficient for modelling aeroservoelastic data, (ii) nonlinear identification provides a parsimonious model description whilst providing a high percent fit for cross-validated data and (iii) the model structure and parameters vary as the flight condition is altered.

Kukreja, Sunil L.↗

Structure Computation of Quiet Spike[Trademark] Flight-Test Data During Envelope Expansion

System identification or mathematical modeling is used in the aerospace community for development of simulation models for robust control law design. These models are often described as linear time-invariant processes. Nevertheless, it is well known that the underlying process is often nonlinear. The reason for using a linear approach has been due to the lack of a proper set of tools for the identification of nonlinear systems. Over the past several decades, the controls and biomedical communities have made great advances in developing tools for the identification of nonlinear systems. These approaches are robust and readily applicable to aerospace systems. In this paper, we show the application of one such nonlinear system identification technique, structure detection, for the analysis of F-15B Quiet Spike(TradeMark) aeroservoelastic flight-test data. Structure detection is concerned with the selection of a subset of candidate terms that best describe the observed output. This is a necessary procedure to compute an efficient system description that may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modeling may be of critical importance for the development of robust parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion, which may save significant development time and costs. The objectives of this study are to demonstrate via analysis of F-15B Quiet Spike aeroservoelastic flight-test data for several flight conditions that 1) linear models are inefficient for modeling aeroservoelastic data, 2) nonlinear identification provides a parsimonious model description while providing a high percent fit for cross-validated data, and 3) the model structure and parameters vary as the flight condition is altered.

Kukreja, Sunil L.↗

Molecular hydrogen emission in NGC 7027

The spatial distribution of the emission in the v = 1-0 S(1) line of H2 in the planetary nebula NGC 7027 is presented. The excited H2 molecules do not fill the same volume as the ionized gas and most likely reside near the outer edge of the nebula. The spatial resolution of the data, 5 arcsec, is insufficient to define the location of the molecules precisely. An upper limit to the strength of the v = 2-1 S(1) line is given; it is low enough to preclude simple ultraviolet fluorescence as the source of excitation. A simple shock model can fit the line ratio data, however, and there is enough energy in the expanding nebula to sustain such a shock. A rough estimate of 1-4 solar masses for the mass in the molecular cloud surrounding NGC 7027 is derived.

Beckwith, S.↗

Two-temperature interpretation of dissociation rate data for N2 and O2

The existing experimental data on dissociation of nitrogen and oxygen obtained using shock-tubes during the 1960's and 1970's are reinterpreted using the two-temperature thermo-chemical model developed recently in order to determine the rate coefficients consistent with the model. In this model, the vibrational-electronic temperature is calculated by integrating a separate conservation equation accounting for the suppression of vibrational energy during dissociation due to preferential removal of high vibrational states. The rate coefficient is assumed to be a function of the geometrically-averaged temperature between the translational-rotational temperature and the vibrational-electronic temperature. By comparing the computed overall and species densities with the experimental data, the rate coefficient values most consistent with the model, and their ranges of uncertainty, are deduced for dissociation of N2 through collisions with N2 or N, and for O2 through collisions with O2, O or N2. It is seen that a single set of such rate coefficients fit all existing experimental data closely. According to the two-temperature model, density and species density are insensitive to the rate coefficients, and so the rate coefficients so determined are uncertain to within a factor of at least 1.5.

Park, Chul↗

Determination of the Charon/Pluto Mass Ratio from Center-of-Light Astrometry

The Charon/Pluto mass ratio is a fundamental but poorly known parameter of the two-body system. Previous values for the mass ratio have ranged from 0.0837 plus or minus 0.0147 (Null et al., 1993, Astron. J. 105, 2319-2335) to 0.1566 plus or minus 0.0035 (Young et al., 1994, Icarus 108,186-199). We report here a new determination of the Charon/Pluto mass ratio, using five sets of groundbased images taken at four sites in support of Pluto occultation predictions. Unlike the Null et al. and Young et A determinations, where the centers of light for Pluto and Charon could be determined separately, this technique examines the motion of the center of light of the blended Pluto-Charon image. We compute the offsets of the observed center-of-light position of Pluto-Charon from the ephemeris position of the system and fit these offsets to a model of the Pluto-Charon system. The least-squares fits to the five data sets agree within their errors, and the weighted mean mass ratio is 0.117 plus or minus 0.006. The effects of errors in the Charon light fraction, semimajor axis, and ephemeris have been examined and are equal to only a small fraction of the formal error from the fit. This result is intermediate between those of Null et al., and Young et al. and matches a new value of 0.124 plus or minus 0.008 by Null and Owen (1996, Astron. J. 111, 1368-1381). The mass ratio and resulting individual masses and densities of Pluto and Charon are consistent with a collisional origin for the Pluto-Charon system.

Foust, Jeffrey A.↗

The Hopkins Ultraviolet Telescope: The Final Archive

The Hopkins Ultraviolet Telescope (HUT) was a 0.9 m telescope and moderate-resolution (Delta)lambda equals 3 A) far-ultraviolet (820-1850 Å) spectrograph that flew twice on the space shuttle, in 1990 December (Astro-1, STS-35) and 1995 March (Astro-2, STS-67). The resulting spectra were originally archived in a nonstandard format that lacked important descriptive metadata. To increase their utility, we have modified the original datareduction software to produce a new and more user-friendly data product, a time-tagged photon list similar in format to the Intermediate Data Files (IDFs) produced by the Far Ultraviolet Spectroscopic Explorer calibration pipeline. We have transferred all relevant pointing and instrument-status information from locally-archived science and engineering databases into new FITS header keywords for each data set. Using this new pipeline, we have reprocessed the entire HUT archive from both missions, producing a new set of calibrated spectral products in a modern FITS format that is fully compliant with Virtual Observatory requirements. For each exposure, we have generated quicklook plots of the fully-calibrated spectrum and associated pointing history information. Finally, we have retrieved from our archives HUT TV guider images, which provide information on aperture positioning relative to guide stars, and converted them into FITS-format image files. All of these new data products are available in the new HUT section of the Mikulski Archive for Space Telescopes (MAST), along with historical and reference documents from both missions. In this article, we document the improved data-processing steps applied to the data and show examples of the new data products.

Hopkins↗

Atmospheric particulate analysis using angular light scattering

Using the light scattering matrix elements measured by a polar nephelometer, a procedure for estimating the characteristics of atmospheric particulates was developed. A theoretical library data set of scattering matrices derived from Mie theory was tabulated for a range of values of the size parameter and refractive index typical of atmospheric particles. Integration over the size parameter yielded the scattering matrix elements for a variety of hypothesized particulate size distributions. A least squares curve fitting technique was used to find a best fit from the library data for the experimental measurements. This was used as a first guess for a nonlinear iterative inversion of the size distributions. A real index of 1.50 and an imaginary index of -0.005 are representative of the smoothed inversion results for the near ground level atmospheric aerosol in Tucson.

Hansen, M. Z.↗

Determination of Chemical Kinetic Rate Constants of a Model for Carbothermal Processing of Lunar Regolith Simulant Using Methane

We have previously developed a chemical conversion model of the carbothermal processing of lunar regolith using methane to predict the rate of production of carbon monoxide. In this carbothermal process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Hydrogen is carried away by the exiting stream of gases and carbon is deposited on the melt surface. The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. In our model, we assume that the flux of carbon deposited is equal to the product of the surface reaction rate constant gamma and the concentration of methane adjacent to the melt surface. Similarly, the rate of consumption of carbon per unit volume in the melt is equal to the product of the melt reaction rate constant k and the concentrations of carbon and metal oxide in the melt. In this paper, we describe our effort to determine gamma and k by comparison of the predictions from our model with test data obtained by ORBITEC (Orbital Technologies Corporation). The concentration of methane adjacent to the melt surface is a necessary input to the model. It is inferred from the test data by a mass balance of methane, adopting the usual assumptions of the continuously-stirred-tank-reactor model, whereby the average concentration of a given gaseous species equals its exit concentration. The reaction rates gamma and k have been determined by a non-linear least-squares fit to the test data for the production of carbon monoxide and the fraction of the incoming methane that is converted. The comparison of test data with our model predictions using the determined chemical kinetic rate constants provides a consistent interpretation of the process over the full range of temperatures, pressures, and methane flow rates used in the tests, thereby increasing our confidence to use the model for scale-up purposes.

Balasubramaniam, R↗

In Hot Pursuit of the Hidden Companion of Carinae: An X-Ray Determination of the Wind Parameters

We present X-ray spectral fits to a recently obtained Chandra grating spectrum of eta Carinae, one of the most massive and powerful stars in the Galaxy and which is strongly suspected to be a colliding wind binary system. Hydrodynamic models of colliding winds are used to generate synthetic X-ray spectra for a range of mass-loss rates and wind velocities. They are then fitted against newly acquired Chandra grating data. We find that due to the low velocity of the primary wind (approximately 500 km per sec), most of the observed X-ray emission appears to arise from the shocked wind of the companion star. We use the duration of the lightcurve minimum to fix the wind momentum ratio at eta = 0.2. We are then able to obtain a good fit to the data by varying the mass-loss rate of the companion and the terminal velocity of its wind. We find that M(sub 2) is approximately 10(exp -5) solar mass per yr and upsilon(sub infinity(sub 2)) is approximately 3000 km per sec. With observationally determined values of approximately 500- 700 km per sec for the velocity of the primary wind, our fit implies a primary mass-loss rate of M(sub 1) approximately 2.5 x 10(exp -4) solar mass per yr. This value is smaller than commonly inferred, although we note that a lower mass-loss rate can reduce some of the problems noted when a value as high as 10(exp -3) solar mass per yr is used. The wind parameters of the companion are indicative of a massive star which may or may not be evolved. The line strengths appear to show slightly sub-solar abundances, although this needs further confirmation. It also appears that the primary star is responsible for the nebula.

Pittard, J. M.↗

Solar neutron emissivity during the large flare on 1982 June 3

For the solar neutron event on June 3, 1982, it is shown here that the combined SMM Gamma Ray Spectrometer and Jungfraujoch neutron monitor data require a time-extended emission of neutrons at the sun with energies of 100 MeV to about 2 GeV. The solar neutron emissivity spectrum is shown to have a strong downward curvature or truncation between 2 and 4 GeV. A Bessel function and truncated power law give acceptable fits to the observational data, but only the power law can explain the rapid rise of the neutron monitor count rate. The integrated emissivity of neutrons above E(n) of 100 MeV is strongly constrained at 8 x 10 to the 28th neutrons/sr and is essentially independent of neutron spectral shape. At neutron energies of about 100 MeV, good agreement is found for both spectral forms with observations of neutron decay protons.

Chupp, E. L.↗